[GRASS-dev] [GRASS GIS] #3161: different results by v.in.ogr --ui and vector import wizzard (data shift)
GRASS GIS
trac at osgeo.org
Fri Sep 16 13:24:30 PDT 2016
#3161: different results by v.in.ogr --ui and vector import wizzard (data shift)
---------------------------+-------------------------
Reporter: hellik | Owner: grass-dev@…
Type: defect | Status: new
Priority: blocker | Milestone: 7.2.0
Component: Vector | Version: svn-trunk
Keywords: vector import | CPU: x86-64
Platform: MSWindows 7 |
---------------------------+-------------------------
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 at 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 at 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 at 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 at 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 at 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.
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3161>
GRASS GIS <https://grass.osgeo.org>
More information about the grass-dev
mailing list