Opened 8 years ago
Closed 8 years ago
#3161 closed defect (worksforme)
different results by v.in.ogr --ui and vector import wizzard (data shift)
| Reported by: | hellik | Owned by: | |
|---|---|---|---|
| Priority: | blocker | Milestone: | 7.2.0 |
| Component: | Vector | Version: | svn-trunk |
| Keywords: | vector import | Cc: | |
| CPU: | x86-64 | Platform: | MSWindows 7 |
Description
follwoing location
g.proj -p -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1
and vector data (shapefile to import):
INFO: Open of `Gesamtgewaessernetz_v11_Tirol_epsg31287.shp'
using driver `ESRI Shapefile' successful.
Layer name: Gesamtgewaessernetz_v11_Tirol_epsg31287
Metadata:
DBF_DATE_LAST_UPDATE=2016-09-16
Geometry: Measured Line String
Feature Count: 5095
Extent: (118287.824753, 197521.746410) - (833235.119008, 544697.312403)
Layer SRS WKT:
PROJCS["MGI_Austria_Lambert",
GEOGCS["GCS_MGI",
DATUM["Militar_Geographische_Institute",
SPHEROID["Bessel_1841",6377397.155,299.1528128]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",49],
PARAMETER["standard_parallel_2",46],
PARAMETER["latitude_of_origin",47.5],
PARAMETER["central_meridian",13.33333333333333],
PARAMETER["false_easting",400000],
PARAMETER["false_northing",400000],
UNIT["Meter",1]]
KURZRID: Integer64 (10.0)
RL: String (2.0)
VERSION: String (5.0)
GEW_ID: String (15.0)
GEW_NAME: String (120.0)
OBJEKT: String (6.0)
STAND: Date (10.0)
EMASST: Integer64 (10.0)
SHAPE_LEN: Real (19.11)
GEW_TYP: String (30.0)
GEW_GRKL: String (30.0)
GEW_LAENGE: Real (19.10)
GEW_NAME2: String (120.0)
srs of location and data for import are matching.
using the vector import wizzard
v.import input=D:\temp\Gesamtgewaessernetz_v11_Tirol layer=Gesamtgewaessernetz_v11_Tirol_epsg31287 output=Gesamtgewaessernetz_v11_Tirol_epsg31287_vimportwizzard WARNING: All available OGR layers will be imported into vector map <Gesamtgewaessernetz_v11_Tirol> Importing <D:\temp\Gesamtgewaessernetz_v11_Tirol> ... Check if OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287> contains polygons... Standard Treiber / Datenbank ist: Treiber: sqlite Datenbank: $GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db Importing 5095 features (OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287>)... ----------------------------------------------------- Erstelle Topologie für die Vektorkarte <Gesamtgewaessernetz_v11_Tirol_epsg31287_vimportwizzard@PERMANENT>... Registriere Primitive... 5095 primitives registered 1505379 Vertices registriert Erzeuge Flächen... 0 areas built 0 isles built Füge Inseln hinzu... Füge Zentroide hinzu... Anzahl der Knoten: 10108 Anzahl der Primitive: 5095 Anzahl der Punkte: 0 Anzahl der Linien: 5095 Anzahl der Grenzen: 0 Anzahl der Zentroide: 0 Anzahl der Flächen: 0 Anzahl der Inseln: 0 Reprojecting <Gesamtgewaessernetz_v11_Tirol_epsg31287_vimportwizzard>... Reprojecting primitives ... Erstelle Topologie für die Vektorkarte <Gesamtgewaessernetz_v11_Tirol_epsg31287_vimportwizzard@data>... Registriere Primitive... 5095 primitives registered 1505379 Vertices registriert Erzeuge Flächen... 0 areas built 0 isles built Füge Inseln hinzu... Füge Zentroide hinzu... Anzahl der Knoten: 10108 Anzahl der Primitive: 5095 Anzahl der Punkte: 0 Anzahl der Linien: 5095 Anzahl der Grenzen: 0 Anzahl der Zentroide: 0 Anzahl der Flächen: 0 Anzahl der Inseln: 0
you can see a reprojecting
Reprojecting primitives ...
although srs of location and import vector data matches
then starting v.in.ogr
v.in.ogr input=D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp layer=Gesamtgewaessernetz_v11_Tirol_epsg31287 output=Gesamtgewaessernetz_v11_Tirol_epsg31287_vinogrui Check if OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287> contains polygons... Importing 5095 features (OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287>)... ----------------------------------------------------- Erstelle Topologie für die Vektorkarte <Gesamtgewaessernetz_v11_Tirol_epsg31287_vinogrui@data>... Registriere Primitive... 5095 primitives registered 1505379 Vertices registriert Erzeuge Flächen... 0 areas built 0 isles built Füge Inseln hinzu... Füge Zentroide hinzu... Anzahl der Knoten: 10108 Anzahl der Primitive: 5095 Anzahl der Punkte: 0 Anzahl der Linien: 5095 Anzahl der Grenzen: 0 Anzahl der Zentroide: 0 Anzahl der Flächen: 0 Anzahl der Inseln: 0
v.info map=Gesamtgewaessernetz_v11_Tirol_epsg31287_vimportwizzard@data
+----------------------------------------------------------------------------+
| Name: Gesamtgewaessernetz_v11_Tirol_epsg31287_vimportwizzard |
| Mapset: data |
| Location: locat_epsg31287 |
| Datenbank: D:\grassdata |
| Titel: |
| Maßstab: 1:1 |
| Name des Erzeugers:hkmyr |
| Organisation: |
| Quelldatum: Fri Sep 16 21:51:42 2016 |
| Timestamp (first layer): none |
|----------------------------------------------------------------------------|
| Kartenformat: native |
|----------------------------------------------------------------------------|
| Kartenart:: Vektor (level: 2) |
| |
| Anzahl der Punkte: 0 Anzahl der Zentroide: 0 |
| Anzahl der Linien: 5095 Anzahl der Grenzen: 0 |
| Anzahl der Flächen: 0 Anzahl der Inseln: 0 |
| |
| Map is 3D: No |
| Anzahl der dblinks: 1 |
| |
| Projektion: MGI / Austria Lambert |
| |
| N: 544693.22594229 S: 197500.09732215 |
| E: 832927.81730718 W: 117982.48429679 |
| |
| Digitalisierungs-Schwellwert:: 0 |
| Kommentar: |
| |
+----------------------------------------------------------------------------+
v.info map=Gesamtgewaessernetz_v11_Tirol_epsg31287_vinogrui@data
+----------------------------------------------------------------------------+
| Name: Gesamtgewaessernetz_v11_Tirol_epsg31287_vinogrui |
| Mapset: data |
| Location: locat_epsg31287 |
| Datenbank: D:\grassdata |
| Titel: |
| Maßstab: 1:1 |
| Name des Erzeugers:hkmyr |
| Organisation: |
| Quelldatum: Fri Sep 16 21:53:48 2016 |
| Timestamp (first layer): none |
|----------------------------------------------------------------------------|
| Kartenformat: native |
|----------------------------------------------------------------------------|
| Kartenart:: Vektor (level: 2) |
| |
| Anzahl der Punkte: 0 Anzahl der Zentroide: 0 |
| Anzahl der Linien: 5095 Anzahl der Grenzen: 0 |
| Anzahl der Flächen: 0 Anzahl der Inseln: 0 |
| |
| Map is 3D: No |
| Anzahl der dblinks: 1 |
| |
| Projektion: MGI / Austria Lambert |
| |
| N: 544697.31240282 S: 197521.74641048 |
| E: 833235.11900777 W: 118287.82475345 |
| |
| Digitalisierungs-Schwellwert:: 0 |
| Kommentar: |
| |
+----------------------------------------------------------------------------+
there is a data shift/reprojection by vector import wizzard!
import via v.in.ogr --ui is correct.
tested with
GRASS Version: 7.3.svn GRASS SVN revision: r69490 Build date: 2016-09-15 Build platform: x86_64-w64-mingw32 GDAL: 2.1.1 PROJ.4: 4.9.3 GEOS: 3.5.0 SQLite: 3.14.1 Python: 2.7.5 wxPython: 2.8.12.1 Platform: Windows-8-6.2.9200 (OSGeo4W)
attaching a screenshot of the data shift.
Attachments (1)
Change History (9)
by , 8 years ago
| Attachment: | check_vimport_vinogrui.png added |
|---|
follow-ups: 2 3 4 comment:1 by , 8 years ago
Replying to hellik:
follwoing location
g.proj -p -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1
The towgs84 parameters could be a problem. GDAL uses similar towgs84 parameters as default while GRASS uses as default
towgs84 : 653,-212,449,0,0,0,0
Since
and vector data (shapefile to import):
INFO: Open of `Gesamtgewaessernetz_v11_Tirol_epsg31287.shp' using driver `ESRI Shapefile' successful. Layer name: Gesamtgewaessernetz_v11_Tirol_epsg31287 Metadata: DBF_DATE_LAST_UPDATE=2016-09-16 Geometry: Measured Line String Feature Count: 5095 Extent: (118287.824753, 197521.746410) - (833235.119008, 544697.312403) Layer SRS WKT: PROJCS["MGI_Austria_Lambert", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",49], PARAMETER["standard_parallel_2",46], PARAMETER["latitude_of_origin",47.5], PARAMETER["central_meridian",13.33333333333333], PARAMETER["false_easting",400000], PARAMETER["false_northing",400000], UNIT["Meter",1]]
does not report towgs84 parameters, GRASS might decide that
srs of location and data for import are matching.
srs of location and data for import are not matching.
Is the output of
g.proj -p
and
g.proj -p georef="D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp"
identical?
If not, as v.import assumes, you might need to provide the EPSG code and the datum transformation parameters as input to v.import (which makes v.import more complicated as it was intended to be).
If the towgs84 parameters as reported by g.proj for the current location and for the input shapefile are not identical, this could be regarded as a bug in GRASS since the vast majority of GIS data is imported via r.in.gdal/v.in.ogr and GDAL/OGR automatically picks a preferred set of datum transformation parameters which GRASS should detect.
In this case the datum transformation parameters selected by GDAL are
towgs : 577.326,90.129,463.919,5.137,1.474,5.297,2.4232
which is similar but not equal to those used in your location
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
In the end, the user has to decide if reprojection is needed or if the source srs and the location srs are similar enough to allow for using v.in.ogr -o (override projection check). After investigating the input data with g.proj georef=<data to be imported> and ogrinfo/gdalinfo.
comment:2 by , 8 years ago
Replying to mmetz:
Replying to hellik:
follwoing location
g.proj -p -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1The towgs84 parameters could be a problem. GDAL uses similar towgs84 parameters as default while GRASS uses as default
towgs84 : 653,-212,449,0,0,0,0
the location was created by the location wizzard by EPSG:31287 and then the transformation for the region of Austria was used.
trying
g.proj -p epsg=31287
in a EPSG:4326 and a EPSG:31287 location (both locations are created by the EPSG code, I get
g.proj -p epsg=31287 -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 towgs84 : 577.326,90.129,463.919,5.137,1.474,5.297,2.4232 no_defs : defined -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1
so GRASS seems to choose the same +towgs84 parameter as the underlying gdal/ogr.
Since
and vector data (shapefile to import):
INFO: Open of `Gesamtgewaessernetz_v11_Tirol_epsg31287.shp' using driver `ESRI Shapefile' successful. Layer name: Gesamtgewaessernetz_v11_Tirol_epsg31287 Metadata: DBF_DATE_LAST_UPDATE=2016-09-16 Geometry: Measured Line String Feature Count: 5095 Extent: (118287.824753, 197521.746410) - (833235.119008, 544697.312403) Layer SRS WKT: PROJCS["MGI_Austria_Lambert", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",49], PARAMETER["standard_parallel_2",46], PARAMETER["latitude_of_origin",47.5], PARAMETER["central_meridian",13.33333333333333], PARAMETER["false_easting",400000], PARAMETER["false_northing",400000], UNIT["Meter",1]]does not report towgs84 parameters, GRASS might decide that
srs of location and data for import are matching.
srs of location and data for import are not matching.
Is the output of
g.proj -p
g.proj -p -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1
the location was created by the location wizzard by EPSG:31287 and then the transformation for the region of Austria was used.
and
g.proj -p georef="D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp"identical?
in an EPSG:4326 location I get:
g.proj -p georef=D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp Versuche mit OGR zu öffnen... ...erfolgreich beendet. -PROJ_INFO------------------------------------------------- name : MGI_Austria_Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined -PROJ_EPSG------------------------------------------------- epsg : 4326 -PROJ_UNITS------------------------------------------------ unit : Meter units : Meters meters : 1
and in an EPSG:31287 location I get
g.proj -p georef=D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp Versuche mit OGR zu öffnen... ...erfolgreich beendet. -PROJ_INFO------------------------------------------------- name : MGI_Austria_Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : Meter units : Meters meters : 1
If not, as v.import assumes, you might need to provide the EPSG code and the datum transformation parameters as input to v.import (which makes v.import more complicated as it was intended to be).
now declaring EPSG code and datum_trans in v.import:
v.import --verbose epsg=31287 datum_trans=2 input=D:\temp\Gesamtgewaessernetz_v11_Tirol layer=Gesamtgewaessernetz_v11_Tirol_epsg31287 output=checkimport Creating temporary location for <D:\temp\Gesamtgewaessernetz_v11_Tirol>... Projektionsinformationen aktualisiert WARNING: All available OGR layers will be imported into vector map <Gesamtgewaessernetz_v11_Tirol> -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1 Importing <D:\temp\Gesamtgewaessernetz_v11_Tirol> ... Die Projektionsinformationen des Eingabedatensatzes und der aktuellen Location scheinen übereinzustimmen. Check if OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287> contains polygons... Using native format Standard Treiber / Datenbank ist: Treiber: sqlite Datenbank: $GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db Importing 5095 features (OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287>)... ----------------------------------------------------- Erstelle Topologie für die Vektorkarte <checkimport@PERMANENT>... Registriere Primitive... 5095 primitives registered 1505379 Vertices registriert Erzeuge Flächen... 0 areas built 0 isles built Füge Inseln hinzu... Füge Zentroide hinzu... Die Topologie wurde erstellt. Anzahl der Knoten: 10108 Anzahl der Primitive: 5095 Anzahl der Punkte: 0 Anzahl der Linien: 5095 Anzahl der Grenzen: 0 Anzahl der Zentroide: 0 Anzahl der Flächen: 0 Anzahl der Inseln: 0 Reprojecting <checkimport>... Eingabe der Projektionsparameter: +proj=lcc +lat_1=49 +lat_2=46 +lat_0=47.5 +lon_0=13.33333333333333 +x_0=400000 +y_0=400000 +no_defs +a=6377397.155 +rf=299.1528128 +towgs84=577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 Eingabe des Einheitenfaktors: 1 Ausgabe der Projektionsparameter: +proj=lcc +lat_1=49 +lat_2=46 +lat_0=47.5 +lon_0=13.33333333333333 +x_0=400000 +y_0=400000 +no_defs +a=6377397.155 +rf=299.1528128 +towgs84=577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 Ausgabe des Einheitenfaktors: 1 Using native format Reprojecting primitives ... Erstelle Topologie für die Vektorkarte <checkimport@data>... Registriere Primitive... 5095 primitives registered 1505379 Vertices registriert Erzeuge Flächen... 0 areas built 0 isles built Füge Inseln hinzu... Füge Zentroide hinzu... Die Topologie wurde erstellt. Anzahl der Knoten: 10108 Anzahl der Primitive: 5095 Anzahl der Punkte: 0 Anzahl der Linien: 5095 Anzahl der Grenzen: 0 Anzahl der Zentroide: 0 Anzahl der Flächen: 0 Anzahl der Inseln: 0
If the towgs84 parameters as reported by g.proj for the current location and for the input shapefile are not identical, this could be regarded as a bug in GRASS since the vast majority of GIS data is imported via r.in.gdal/v.in.ogr and GDAL/OGR automatically picks a preferred set of datum transformation parameters which GRASS should detect.
In this case the datum transformation parameters selected by GDAL are
towgs : 577.326,90.129,463.919,5.137,1.474,5.297,2.4232which is similar but not equal to those used in your location
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232In the end, the user has to decide if reprojection is needed or if the source srs and the location srs are similar enough to allow for using v.in.ogr -o (override projection check). After investigating the input data with g.proj georef=<data to be imported> and ogrinfo/gdalinfo.
for me it's not clear, as I've choosen GUI -> import common formats and not GUI -> import of common formats with reprojection. in this case I'm sure/assume my data is in the same srs and import should just works.
it's quite confusing as v.in.ogr --ui imports the data correctly without the need of override projection check.
for me it's a regression, as I have to type v.in.ogr --ui to get correct results for a simple task of import data in the same srs as the location.
comment:3 by , 8 years ago
Replying to mmetz:
In this case the datum transformation parameters selected by GDAL are
towgs : 577.326,90.129,463.919,5.137,1.474,5.297,2.4232which is similar but not equal to those used in your location
towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
g.proj epsg=31287 datum_trans=-1 --- 1 Used in whole hermannskogel region towgs84=653.000,-212.000,449.000 Default 3-Parameter Transformation (May not be optimum for older datums; use this only if no more appropriate options are available.) --- 2 Used in Austria towgs84=577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 Accuracy approx. 1.5m --- 3 Used in Czech Republic towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 --- 4 Used in Slovakia towgs84=485.021,169.465,483.839,7.786342,4.397554,4.102655,0 --- 5 Used in Slovenia towgs84=426.9,142.6,460.1,4.91,4.49,-12.42,17.1
mmhh, g.proj gives slightly different +towgs84 parameters as you mention it for GDAL/OGR.
System Info GRASS Version: 7.3.svn GRASS SVN revision: r69510 Build date: 2016-09-17 Build platform: x86_64-w64-mingw32 GDAL: 2.1.1 PROJ.4: 4.9.3 GEOS: 3.5.0 SQLite: 3.14.1 Python: 2.7.5 wxPython: 2.8.12.1 Platform: Windows-8-6.2.9200 (OSGeo4W)
follow-up: 5 comment:4 by , 8 years ago
Replying to mmetz:
In the end, the user has to decide if reprojection is needed or if the source srs and the location srs are similar enough to allow for using v.in.ogr -o (override projection check). After investigating the input data with g.proj georef=<data to be imported> and ogrinfo/gdalinfo.
as mentioned before, there is no override projection check needed by import via v.in.ogr --ui.
if you want to try
- data can be downloaded river net;
- then ogr2ogr -s_srs EPSG:31254 -t_srs EPSG:31287
- create location by EPSG:31287 and datum_trans 2
- then import GUI -> import common formats
follow-up: 6 comment:5 by , 8 years ago
Replying to hellik:
Replying to mmetz:
In the end, the user has to decide if reprojection is needed or if the source srs and the location srs are similar enough to allow for using v.in.ogr -o (override projection check). After investigating the input data with g.proj georef=<data to be imported> and ogrinfo/gdalinfo.
as mentioned before, there is no override projection check needed by import via v.in.ogr --ui.
if you want to try
- data can be downloaded river net;
- then ogr2ogr -s_srs EPSG:31254 -t_srs EPSG:31287
- create location by EPSG:31287 and datum_trans 2
- then import GUI -> import common formats
I tested. The location has the same parameters like yours:
> g.proj -p -PROJ_INFO------------------------------------------------- name : MGI / Austria Lambert datum : hermannskogel ellps : bessel proj : lcc lat_1 : 49 lat_2 : 46 lat_0 : 47.5 lon_0 : 13.33333333333333 x_0 : 400000 y_0 : 400000 no_defs : defined towgs84 : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232 -PROJ_EPSG------------------------------------------------- epsg : 31287 -PROJ_UNITS------------------------------------------------ unit : meter units : meters meters : 1
Now test if the SRS of the shapefile is identical
> v.in.ogr in=Gesamtgewaessernetz_v11_Tirol_epsg31287 layer=Gesamtgewaessernetz_v11_Tirol -j Projection of input dataset and current location appear to match
fine, now v.import
> v.import in=Gesamtgewaessernetz_v11_Tirol_epsg31287 out=Gesamtgewaessernetz_v11_Tirol_epsg31287 [...] Input <Gesamtgewaessernetz_v11_Tirol_epsg31287> successfully imported without reprojection
All fine, no reprojection. What is your output of v.in.ogr -j?
follow-up: 7 comment:6 by , 8 years ago
Replying to mmetz:
All fine, no reprojection. What is your output of
v.in.ogr -j?
v.in.ogr -j input=D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp Die Projektionsinformationen des Eingabedatensatzes und der aktuellen Location scheinen übereinzustimmen.
they are also matching here.
will test it with today's winGRASS daily.
comment:7 by , 8 years ago
Replying to hellik:
Replying to mmetz:
All fine, no reprojection. What is your output of
v.in.ogr -j?v.in.ogr -j input=D:\temp\Gesamtgewaessernetz_v11_Tirol\Gesamtgewaessernetz_v11_Tirol_epsg31287.shp Die Projektionsinformationen des Eingabedatensatzes und der aktuellen Location scheinen übereinzustimmen.they are also matching here.
will test it with today's winGRASS daily.
now tested with
System Info GRASS Version: 7.3.svn GRASS SVN revision: r69517 Build date: 2016-09-19 Build platform: x86_64-w64-mingw32 GDAL: 2.1.1 PROJ.4: 4.9.3 GEOS: 3.5.0 SQLite: 3.14.1 Python: 2.7.5 wxPython: 2.8.12.1 Platform: Windows-8-6.2.9200 (OSGeo4W)
v.import input=D:\temp\Gesamtgewaessernetz_v11_Tirol layer=Gesamtgewaessernetz_v11_Tirol_epsg31287 output=Gesamtgewaessernetz_v11_Tirol_epsg31287 WARNING: All available OGR layers will be imported into vector map <data_epsg31287> Check if OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287> contains polygons... Importing 5095 features (OGR layer <Gesamtgewaessernetz_v11_Tirol_epsg31287>)... ----------------------------------------------------- Erstelle Topologie für die Vektorkarte <Gesamtgewaessernetz_v11_Tirol_epsg31287@data>... Registriere Primitive... 5095 primitives registered 1505379 Vertices registriert Erzeuge Flächen... 0 areas built 0 isles built Füge Inseln hinzu... Füge Zentroide hinzu... Anzahl der Knoten: 10108 Anzahl der Primitive: 5095 Anzahl der Punkte: 0 Anzahl der Linien: 5095 Anzahl der Grenzen: 0 Anzahl der Zentroide: 0 Anzahl der Flächen: 0 Anzahl der Inseln: 0 Input <D:\temp\Gesamtgewaessernetz_v11_Tirol> successfully imported without reprojection
now it works also without reprojection.
it seems some temporary instability. closing ticket as worksforme.
thanks for testing.
comment:8 by , 8 years ago
| Resolution: | → worksforme |
|---|---|
| Status: | new → closed |
screenhot of data shift: green correctly imported, red false