De boel is uitgeplozen en herleidt tot het volgende:
Sinds versie 3.7.0 gaat GDAL ervan uit dat geometrie van CityGML in 3D is conform de CityGML standaard:
(OGC City Geography Markup Language (CityGML) 3.0 Conceptual Model Users Guide):
Spatial properties of all CityGML feature types are represented using the geometry classes defined in ISO 19107. Spatial representations can have 0-, 1-, 2-, or 3-dimensional extents depending on the respective feature type and Levels of Detail (LOD). The LOD concept is discussed in [ug-levels-of-detail-section] and [ug-geometry-lod-section]. With only a few exceptions, all geometries must use 3D coordinate values.
Voor de een aantal lijngeometrie elementen wordt het attribuut srsDimension niet gebruikt. GDAL gaat er dan van uit dat de poslist steeds x y z waarden bevat ipv xy waarden.
Het toevoegen van srsDimension="2" aan het posList element lost het probleem op.
Voorbeeld: in bgt_ondersteunendWegdeel komt een element voor met ifentificatie:
<imgeo:identificatie>
<imgeo:NEN3610ID>
<imgeo:namespace>NL.IMGeo</imgeo:namespace>
<imgeo:lokaalID>G0274.6d2f97620c5346af9f092ee268d2c6a2</imgeo:lokaalID>
</imgeo:NEN3610ID>
</imgeo:identificatie>
Deze bevat een element:
<imgeo:geometrie2dOndersteunendWegdeel><gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="urn:ogc:def:crs:EPSG::28992"><gml:exterior><gml:LinearRing><gml:posList srsDimension="2">187888.305 444056.518 …
Hier wordt de srsDimension expliciet genoemd zoals vereist door de standaard.
Daarnaast is er ook ook:
<imgeo:kruinlijnOndersteunendWegdeel><gml:LineString xmlns:gml="http://www.opengis.net/gml"><gml:posList>187888.305 444056.518 …
Hier wordt de srsDimension niet expliciet genoemd hetgeen leidt tot foutief inlezen (GDAL verwacht een posList met x y z ipv x y).
Wijzigen tiot
<imgeo:kruinlijnOndersteunendWegdeel><gml:LineString xmlns:gml="http://www.opengis.net/gml"><gml:posList srsDimension="2">187888.305 444056.518 …
leidt wel tot correct inlezen.
Vraag aan de beheerders van BGT en/ of PDOK: Kan de gml conform de standaarden uitgeleverd worden met het srsDimension attribuut op alle geometrie?