VDV/IDF Driver doesn't create Layer "Node" with Z-Values (3D Data)

OSGeo / gdal

GDAL is an open source MIT licensed translator library for raster and vector geospatial data formats.

https://gdal.org

Other

4.9k stars 2.55k forks source link

VDV/IDF Driver doesn't create Layer "Node" with Z-Values (3D Data) #964

Closed benehuber closed 6 years ago

benehuber commented 6 years ago

Expected behavior and actual behavior.

We try to convert Data from the Austrian official open government street graph in "INTREST Data format" with Z Values (3D) in a Geopackage. Conversition works fine, but Layer "Node" doesnt have Z-Values. Layers "LinkCoordinate" and "Node" are (as apspected) Point Layers with Z-Values.

Steps to reproduce the problem.

Download IDF File from http://open.gip.gv.at/ogd/3_routingexport_ogd.zip (https://www.data.gv.at/katalog/dataset/intermodales-verkehrsreferenzsystem-osterreich-gip-at-beta)
convert it to Geopackage with ogr2ogr -f "GPKG" 3_routingexport_ogd.gpkg 3_routingexport_ogd.txt -a_srs EPSG:4326 -zfield Z -dim XYZ
Geopackage will be created
Layers LinkCoordinate and Node are Point Layers with Z-Values
Layer Node will be created without Z-Values

Operating system

Windows 7/ Server 2012 R2

GDAL version and provenance

"release-1911-x64-gdal-2-3-1-mapserver-7-2-0" from gisinternals

benehuber commented 6 years ago

Is there a possibility to check the result?

rouault commented 6 years ago

Is there a possibility to check the result?

What do you mean ? Test the updated master version ? Yes, you can build from source. Or on Windows, use builds at http://gisinternals.com/development.php (the 1900 ones, since the others seem to be not up to date)

benehuber commented 6 years ago

Thank you very much for your quick response. I tried the conversion with http://download.gisinternals.com/sdk/downloads/release-1911-x64-gdal-mapserver.zip (File Date 2018-10-04 21:11:07) and got Geopackage where layer "line" was defined with z-values. But when I select one feature (line) not all vertices have a z-values! Two points at the beginning and one point at the end have a z-values. All vertices between have 0 as z-value. In Screenshot you can see the list of vertices of one line in the editing mode of Arcgis Pro. Therefore, it was necessary to contert the link layer from GPKG to ESRI FGDB.

rouault commented 6 years ago

@benehuber You're right. I've just pushed a fix for those z=0 components

benehuber commented 6 years ago

Thank you for the fix. Now the points have Z-values. But we have another Problem: there must something wrong with the Geopackage. We are not able to copy the LINK layer in an ESRI filegeodatabase because of following error:

Before it was no Problem to copy the geopackage layer to a ESRI filegeodatabase (FGDB). If we convert the IDF File to ESRI Shape Files we are also able to import them to FGDB.