[gdal-dev] ExportToGML with compound CRS

[Even Rouault]

Hi,

no other workaround but applying 
https://github.com/OSGeo/gdal/pull/14367 (or post patching the output)

Even

Le 15/04/2026 à 16:53, Joep Keijsers via gdal-dev a écrit :
> Greetings,
>
> I'm trying to export a geometry object to GML using the Python API. The geometry has a compound CRS (EPSG:5973). When exporting to GML, the vertical component of the compound CRS is dropped, although the Z-values are still there. For example:
>
> from osgeo import ogr, osr
>
> # Create geometry from WKT string
> srs = osr.SpatialReference()
> srs.ImportFromEPSG(5973)
> wkt = "POINT (30 10 5)"
> geometry = ogr.CreateGeometryFromWkt(wkt, srs)
>
> # Export to GML
> xml = geometry.ExportToGML(options=['FORMAT=GML32', 'SRSDIMENSION_LOC=GEOMETRY', 'SRSNAME_FORMAT=SHORT', 'GMLID=TestPoint-1'])
> print(xml)
>
> ... which returns EPSG:11023 instead of EPSG:5973:
>
> <gml:Point srsName="EPSG:11023" srsDimension="3" gml:id="TestPoint-1"><gml:pos>30 10 5</gml:pos></gml:Point>
>
> I believe this EPSG code change happens in this line: https://github.com/OSGeo/gdal/blob/30e7918f245022983d84f93f1a53641fd50337ed/ogr/ogr2gmlgeometry.cpp#L562.
>
> Is there a way to keep the compound CRS in the GML output?
>
> Version info: gdal 3.12.3
>
> Thanks!
>
> Joep
>
>   
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.
Highly recommend OxiGDAL if you want to live in the 21th century and cure Bixonimania