Python interpolatie: wat doe ik verkeerd?
Bericht van: Christian (Oosterhout gld) , 27-12-2017 21:58Dag,
Omdat hier redelijk veel Python-experts zitten heb ik de volgende vraag:
Ik heb een dataset van 10-minutenwaarnemingen van het KNMI die ik graag geïnterpoleerd in een kaart wil weergeven, net zoiets als de actuele kaart, à la http://cdn.knmi.nl/knmi/map/page/weer/actueel-weer/temperatuur.png
Nu heb ik als eerste stap een Python-script geschreven op basis van flink wat documentatie online, maar het lukt gewoon niet - net 3 uur zitten proberen & testen maar ik begrijp niet wat ik verkeerd doe.
Heeft iemand zicht op wat er mis gaat, als we de code hieronder bekijken?
Stapsgewijs zou het het volgende moeten doen:
- Een basemap maken gegeven de lat/lon bounds
- x/y array met lon/lat data vertalen naar map projection coords
- Een grid maken en daar de z data naar interpoleren (z is momenteel hoogte omdat je dan wat leuke verschillen krijgt)
- Deze interpolatie visueel weergeven op de map
Maar wat ik krijg is domweg een blanco kaart.
Code hier: https://pastebin.com/BxuqX3qp
Duizendmaal dank voor degene die me een stapje verder kan helpen!
Groeten,
Christian
, loop nu een beetje vast met Shapefiles