Ik ben bezig met een tool in python om bepaalde infrastructuur beter te herkennen. hiervoor is de 25cm resolutie infrarood dataset ideaal. Ik wil deze graag benaderen via python (waarin de rest van de tool geprogrammeerd is) maar het lukt me niet om de WMTS op een handige manier te benaderen. In essentie heb ik een shapefile (polygon) waarvan ik de IR luchtfoto van wil dowloaden. Kan dit en wat is hier de handigste manier voor?
Ik zou dit zelf met deze Python bibliotheken oplossen:
- Fiona: https://fiona.readthedocs.io/
- Shapely: https://shapely.readthedocs.io/ en
- OWSLib: https://owslib.readthedocs.io/
Met Fiona kun je een Shapefile inlezen.
Met Shapely kun je een “bounding box” van de ingelezen Polygon bepalen, met de envelope() functie.
Met OWSLib kun je web-services bevragen als WMTS. Maar dan moet je de tile-nummers en zoomlevel bepalen. Mogelijk kun je dan beter die resulterende ‘envelope’ (bbox coördinaten) gebruiken om uit de IR Luchtfoto WMS een image te halen. Maar dan moet wel de lengte/breedte verhouding kloppen, plus daar zijn beperkingen aan qua opvraagbare size.
Dit in grote lijnen.