Is natuurlijk voor een grensbewoner nuttig dat er ook wat van het buitenland op staat. Heb je helemaal gelijk. Sorry ik was even gefocussed op de programmeeruitdaging om precies op de grens te knippen.

En natuurlijk is er verschil bij gebruik van een bbox of polygoon met weinig of veel punten. Ik hoop dat osmosis zo slim is om eenmalig voor een polygoon een bbox te berekenen en vervolgens bij isinarea() aanroepen eerst de lat,lon tegen de bbox te houden. Scheelt echt veel. En dan pas de polygoon zelf langs te gaan als dat dan nog moet.

Om 100.000 's te localiceren gebruikte mijn programma 1:12 m:s met benelux.gpx (24424 punten) en 0:13 m:s voor jouw polygoon(60 punten).

Maar met een truc kan ik toch mijn benelux.gpx gebruiken en kom dan op 0:27. Dus maar twee keer zoveel tijd nodig en haartje precies.

Hierop zie je drie polygonen. In blauw de preciese benelux.gpx en dan nog een outherpolygoon (rood en 60 ptn) en een innerpolygoon (groen en 100 ptn).

De truc is dan dat ik eerst controleer of die buiten het outherpolygoon valt en zo ja dan ben ik klaar. Daarna of die binnen het innerpolygoon valt en zoja dan ben ik ook klaar. Slecht voor lat,lon’s in de grensstreek moet ik dan de polygoon waar het om gaat langs.

Die geofabrik osm files gebruik ik ook wel. Alleen maak ik er geen kaarten van dus valt me ook niet op dat het een knoeiboel is aan de grenzen.