<div dir="ltr">Beste Osgeo.nl lijst,<div><br></div><div>Allereerst dank dat jullie naar mijn verzoek zouden willen kijken: mijn naam is Rein van 't Veer en ik werk bij Geodan, gedetacheerd bij het Kadaster. <span style="line-height:1.5">In teamverband werk ik aan het project GRID: een vernieuwingsproject om geodata beter beschikbaar te krijgen.</span></div><div><span style="line-height:1.5"><br></span></div><div><span style="line-height:1.5">Eén van de doelen van het GRID-project is om de onlangs vrijgegeven Digitale Kadastrale Kaart (DKK, handig per provincie te downloaden vanaf </span><a href="https://www.pdok.nl/nl/producten/pdok-downloads/basis-registratie-kadaster/kadastrale-kaart">https://www.pdok.nl/nl/producten/pdok-downloads/basis-registratie-kadaster/kadastrale-kaart</a>)<span style="line-height:1.5"> op andere manieren te ontsluiten - als REST API en als Linked Data. Maar we lopen er even tegenaan dat de dataset, gedownload als GML de nodige uitdagingen bevat. Dit probleem leg ik graag voor aan deze open groep, aangezien het antwoord erop, naar ik hoop, voor anderen waardevol kan zijn die de DKK willen gebruiken.</span></div><div><span style="line-height:1.5"><br></span></div><div><span style="line-height:1.5">Zo bevat de 'laag' Percelen per perceelsobject twee typen geometrie: een punt-object en (naar ik in de ruwe GML kan zien), een polygoongeometrie. We zijn natuurlijk vooral geïnteresseerd in de laatste, maar deze wordt door QGIS (via GDAL) gevoeglijk genegeerd. GDAL lijkt maar moeizaam om te kunnen gaan met meer dan één geometrie per object.</span></div><div><span style="line-height:1.5"><br></span></div><div><span style="line-height:1.5">Vraag is nu: hoe krijg ik de polygonen (de </span>brk_krt:begrenzingPerceel) <span style="line-height:1.5">geconverteerd naar GeoJSON? En hoe kan ik niet de punten, maar de polygonen vanuit GML openen in QGIS? GDAL lijkt helaas een voorkeur te hebben voor de puntgeometrieën. De pagina </span><a href="http://www.gdal.org/drv_geojson.html" style="line-height:1.5">http://www.gdal.org/drv_geojson.html</a><span style="line-height:1.5"> geeft wel enige aanwijzing voor conversie met behulp van de omgevingsvariabele </span><span style="line-height:1.5">GEOMETRY_AS_COLLECTION=YES, maar dit lijkt niet het gewenste effect te hebben op bijvoorbeeld een opdracht als: </span></div><div><br></div><div><div>$ ogrinfo Perceel.gml</div><div>Had to open data source read-only.</div><div>INFO: Open of `Perceel.gml'</div><div>      using driver `GML' successful.</div><div>1: Perceel (Point)</div></div><div><br></div><div>$ ogr2ogr -f GeoJSON Perceel.NO.geojson Perceel.gml</div><div><br></div><div>$ ogrinfo Perceel.NO.geojson</div><div>ERROR 4: GeoJSON Driver doesn't support update.</div><div>Had to open data source read-only.</div><div>INFO: Open of `Perceel.NO.geojson'</div><div>      using driver `GeoJSON' successful.</div><div>1: OGRGeoJSON (Point)</div><div><br></div><div>$ export GEOMETRY_AS_COLLECTION=YES</div><div><div><br></div><div><span style="line-height:1.5">$ ogrinfo Perceel.YES.geojson</span><br></div></div><div><div>ERROR 4: GeoJSON Driver doesn't support update.</div><div>Had to open data source read-only.</div><div>INFO: Open of `Perceel.geojson'</div><div>      using driver `GeoJSON' successful.</div><div>1: OGRGeoJSON (Geometry Collection)</div></div><div><br></div><div><div>De Geometry Collection bevat vervolgens enkel punten. <span style="line-height:1.5">GDAL onderkent geen polygonen in de dataset. Hoe kom ik er op een elegante (automatische) manier achter welke typen geometrie er in de dataset zitten? Iemand enig idee hoe dit benaderd kan worden? Ik heb wel een vergelijkbare melding gevonden op </span><a href="http://gis.stackexchange.com/questions/37217/why-are-the-polygons-in-my-gml-not-rendered" style="line-height:1.5">http://gis.stackexchange.com/questions/37217/why-are-the-polygons-in-my-gml-not-rendered</a><span style="line-height:1.5"> maar het handmatig aanpassen van het GDAL .gfs is niet echt een elegante oplossing. Moet ik met Sax of een ander olifantenmiddel de puntgeometrieën eruit filteren? Of kan er een andere library benut worden zoals de Java(Script) Topology Suite? </span></div></div><div><br></div><div>Alles suggesties zijn welkom, alvast veel dank,</div><div>Rein van 't Veer</div></div>