Lange smalle percelen selecteren (vermoedelijke wegen)

Gepubliceerd door Marloes Decraemer op

Soms vind je een methode die werkt als snelle oplossing maar waarvan je vermoed (of hoopt?) dat er nog een betere manier in de wereld bestaat… In dit geval hadden we een snelle oplossing nodig om wegen die niet als openbaar domein in het GRB opgenomen zijn te identificeren. We waren dus op zoek naar lange smalle percelen.

Na wat google- en zoekwerk kregen we enkele ideetjes (waaronder deze) en werkten we de volgende methode uit:

Het idee is om de oppervlakte van het perceel zoals het is te vergelijken met wat de oppervlakte van het perceel zou zijn als het een lange smalle vorm had, berekend vanaf de omtrek.

Er van uitgaande dat een weg tussen de 1 en 10 meter breed is berekende we de fictieve oppervlakte als volgt:

De omtrek van het perceel (LENGTE) delen we door twee en verminderen we met de voorop genomen breedte. Zo bepalen we de enkele lengte van de kandidaat weg. Wanneer we dit vermenigvuldigen met de breedte van de weg krijgen we de oppervlakte van ons fictief wegperceel.
CALC1 : ((“LENGTE” /2)-1)*1
CALC2 : ((“LENGTE” /2)-10)*10

Hierna vergelijken we de reële oppervlakte van de percelen met de calculaties als volgt:
“OPPERVL” > “CALC1” and “OPPERVL” < “CALC2”

Het resultaat is niet perfect, maar geeft een goede indicatie en een vlotte start om verder op te werken. speel vooral eens met smaller percelen dan 10

Heb je een beter idee of kennis van een handige functie (in ArcGIS pro, Qgis of FME) die dit voor je kan doen, laat het zeker weten in de comments!

Met dank aan Bart Boute!


4 reacties

Mus Luc · 9 maart 2023 op 07:40

In de kadasterdatabank zoeken op aard = weg

    Marloes Decraemer · 9 maart 2023 op 08:37

    Dat is een heel fijne suggestie en zeer nuttig voor wie toegang heeft tot de kadasterdatabank, dankje! Als intercommunale hebben we spijtig genoeg geen toegang tot de legger.

Kris Taillaert · 14 maart 2023 op 14:26

Lukt idd ook via FME door als volgt te werk te gaan:
1 Neem als reader de ADP’s
2 doe een BoundingBoxReplacer (met vb “lengte” als lange zijde en “breedte” als korte zijde)
3 gebruik een expressionevaluator om ’t “verschil” weer te geven (@Value(lengte)-@Value(breedte))
4 gebruik vervolgens een AttributeRangeFilter (behouden van ’t verschil tussen vb 400 en 10000)
5 doe een FeatureMerger om ’t resultaat aan de originele percelen te koppelen (join vb on OIDN)
6 ’t resultaat geeft (bij Brugge) zo’n 85 percelen
(best wat spelen met de laagste waarde in de filter)

    Marloes Decraemer · 28 april 2023 op 16:19

    Dankje Kris,!

Geef een reactie

Avatar plaatshouder

Het e-mailadres wordt niet gepubliceerd. Vereiste velden zijn gemarkeerd met *

Deze site gebruikt Akismet om spam te verminderen. Bekijk hoe je reactie-gegevens worden verwerkt.