<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thanks for the quick patch! After building from the PR branch, I can confirm that it correctly returns the compound CRS for both approaches.</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Joep</div>
<div id="Signature">
<table cellspacing="0" cellpadding="0" border="0">
<tbody>
<tr>
<td>
<div id="x_LSI_marker">
<div style="font-size: 0px;"> </div>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Rahkonen Jukka <jukka.rahkonen@maanmittauslaitos.fi><br>
<b>Sent:</b> 15 April 2026 17:28<br>
<b>To:</b> Joep Keijsers <joep.keijsers@sogelink.com>; gdal-dev@lists.osgeo.org <gdal-dev@lists.osgeo.org>; Even Rouault <even.rouault@spatialys.com><br>
<b>Subject:</b> Re: [gdal-dev] ExportToGML with compound CRS</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">[You don't often get email from jukka.rahkonen@maanmittauslaitos.fi. Learn why this is important at
<a href="https://aka.ms/LearnAboutSenderIdentification">https://aka.ms/LearnAboutSenderIdentification</a> ]<br>
<br>
Hi<br>
<br>
Even made a fix already but here is a command for testing with the old school CLI.<br>
<br>
ogr2ogr -f gml -a_srs EPSG:5973 compound.gml :memory: -sql "select 1 as id, st_geomfromtext('POINTZ (1 1 1)') as geom" -nln test<br>
<br>
"boundedBy" seems to receive<br>
srsName="urn:ogc:def:crs:EPSG::5973 but the geometry itself srsName="urn:ogc:def:crs:EPSG::11023"<br>
<br>
-Jukka Rahkonen-<br>
<br>
________________________________________<br>
Lähettäjä: gdal-dev <gdal-dev-bounces@lists.osgeo.org> käyttäjän Even Rouault via gdal-dev <gdal-dev@lists.osgeo.org> puolesta<br>
Lähetetty: Keskiviikko 15. huhtikuuta 2026 18.25<br>
Vastaanottaja: Joep Keijsers <joep.keijsers@sogelink.com>; gdal-dev@lists.osgeo.org <gdal-dev@lists.osgeo.org><br>
Aihe: Re: [gdal-dev] ExportToGML with compound CRS<br>
<br>
<br>
HUOM! Ulkoinen lähettäjä. Älä avaa linkkejä tai liitteitä, ellet tunnista lähettäjää.<br>
OBS! Extern avsändare. Öppna inte länkar eller bilagor om du inte känner igen avsändaren.<br>
NOTE! External sender. Do not open links or attachments unless you recognize the sender.<br>
<br>
<br>
<br>
<br>
Hi,<br>
<br>
no other workaround but applying<br>
<a href="https://github.com/OSGeo/gdal/pull/14367">https://github.com/OSGeo/gdal/pull/14367</a> (or post patching the output)<br>
<br>
Even<br>
<br>
Le 15/04/2026 à 16:53, Joep Keijsers via gdal-dev a écrit :<br>
> Greetings,<br>
><br>
> 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:<br>
><br>
> from osgeo import ogr, osr<br>
><br>
> # Create geometry from WKT string<br>
> srs = osr.SpatialReference()<br>
> srs.ImportFromEPSG(5973)<br>
> wkt = "POINT (30 10 5)"<br>
> geometry = ogr.CreateGeometryFromWkt(wkt, srs)<br>
><br>
> # Export to GML<br>
> xml = geometry.ExportToGML(options=['FORMAT=GML32', 'SRSDIMENSION_LOC=GEOMETRY', 'SRSNAME_FORMAT=SHORT', 'GMLID=TestPoint-1'])<br>
> print(xml)<br>
><br>
> ... which returns EPSG:11023 instead of EPSG:5973:<br>
><br>
> <gml:Point srsName="EPSG:11023" srsDimension="3" gml:id="TestPoint-1"><gml:pos>30 10 5</gml:pos></gml:Point><br>
><br>
> I believe this EPSG code change happens in this line: <a href="https://github.com/OSGeo/gdal/blob/30e7918f245022983d84f93f1a53641fd50337ed/ogr/ogr2gmlgeometry.cpp#L562">
https://github.com/OSGeo/gdal/blob/30e7918f245022983d84f93f1a53641fd50337ed/ogr/ogr2gmlgeometry.cpp#L562</a>.<br>
><br>
> Is there a way to keep the compound CRS in the GML output?<br>
><br>
> Version info: gdal 3.12.3<br>
><br>
> Thanks!<br>
><br>
> Joep<br>
><br>
><br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
> gdal-dev@lists.osgeo.org<br>
> <a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
--<br>
<a href="http://www.spatialys.com/">http://www.spatialys.com/</a><br>
My software is free, but my time generally not.<br>
Highly recommend OxiGDAL if you want to live in the 21th century and cure Bixonimania<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
gdal-dev@lists.osgeo.org<br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
<br>
</div>
</span></font></div>
</body>
</html>