Bij het samenvoegen van polygonen uit de top10nl gpkg kwam ik erachter dat de hoekpunten van aangrenzende polygonen niet altijd exact hetzelfde zijn. Afgerond op 3 cijfers zijn de punten exact hetzelfde, dus met de data uit de GML-leveringen zal het wel ok zijn. Maar in de GPKG is iets raars met de afronding gebeurd.
Het bestand heb ik hier gedownload:
https://service.pdok.nl/brt/topnl/atom/downloads/top10nl_Compleet.gpkg
Hier de coordinaten van 4 punten van polygonen a en b uit de laag top10nl_wegdeel_vlak die identiek zouden moeten zijn:
a0 123937.67299999999522697 488709.71799999999348074,
b2 123937.67299999999522697 488709.71799999900395051,
a4 123934.82399999999324791 488705.20899999997345731,
b5 123934.82400000099733006 488705.20899999898392707,
a5 123935.63499999999476131 488706.58399999997345731,
b4 123935.63500000200292561 488706.58399999898392707,
a6 123936.36400000000139698 488707.70500000001629815,
b3 123936.36400000000139698 488707.70499999797903001,
Hier linkjes naar de betreffende polygonen uit mijn voorbeeld, maar hetzelfde geldt voor vele andere punten in de gpkg:
polygon a
https://api.pdok.nl/brt/top10nl/ogc/v1/collections/wegdeel_vlak/items?limit=1000&lokaal_id=128028701
polygon b
https://api.pdok.nl/brt/top10nl/ogc/v1/collections/wegdeel_vlak/items?limit=1000&lokaal_id=128028700
Nu weet ik dat binaire dataformaten vaak dit soort verschijnselen vertonen ver achter de komma, maar ik zou verwachten dat die afwijking dan steeds exact hetzelfde is voor elke identieke input.
Is dit probleem bekend? En kan het worden opgelost bij de bron?
Ik heb overigens QGIS gebruikt om de coordinaten te bepalen die ik hierboven post. Maar voor het samenvoegen van de polygonen gebruik ik direct sqlite+spatialite.
Groet,
Raymond
). En zou je eigenlijk ook de gelijkheid van geometrie moeten testen met een tolerantie, ipv ordinaten als getal met elkaar… maar dan spreekt de landmeter/geodeet in mij…