Re: Dank 🙂, loop nu een beetje vast met Shapefiles

Bericht van: Rutger (Meppel) , 30-12-2017 09:26 

Dank voor de bruikbare tips Rutger en Bram!

Ik ben met Shapefiles nu een heel eind verder geraakt & heb nu onderstaand resultaat - ook een stuk sneller t.o.v. Basemap.

Maar ik loop nu vast; delen van land worden 'overspoeld' door het water (zie afbeelding).

Ik gebruik pyshp & gebruik de code van dit artikel: http://www.geophysique.be/2013/02/12/matplotlib-basemap-tutorial-10-shapefiles-unleached-continued/

Ik gebruik shapefiles van België, Frankrijk, Duitsland & de Noordzee om eea in te tekenen. Die van de Noordzee lijkt echter niet goed 'dicht' en als gevolg overspoelen de Randstad en de Wadden 🙂 

Hoe hebben jullie dit opgelost bij het gebruiken van Shapefiles (bv. die kaartjes van Rutger)? Andere shapefile gebruikt (ik gebruik de IHO Sea Areas en dan een met QGIS gemaakte uitsnede van het Noordzee-vlak (zodat alle andere zeeën wegvallen ivm snelheid: http://www.marineregions.org/downloads.php#iho)

Dank alvast 🙂 

Als je die Shapefile in QGIS bekijkt zie je dat die wel goed is, daar ligt het niet aan. Ik zou liever de meren (zie bv Grevelingenmeer) er ook nog bij hebben etc, maar dat is niet het probleem. Het lijkt mij, op basis van je afbeelding, dat je de 'inner rings' van de polygonen of negeert of als 'outer rings' behandeld. Bijvoorbeeld dat stuk bij de Alblasserwaard, dat is duidelijk een 'inner ring'.

Zo ziet die "World_Seas_IHO_v2.shp" eruit, daar zit niet een dergelijk gat in.

Ik gebruik altijd andere data en methoden, afhankelijk van het doel. Voor visualiseren met Bokeh is GeoJSON bijvoorbeeld een handiger formaat dan een ESRI Shapefile. Als je iets als OGR gebruikt maakt het formaat niet zo heel veel uit verder.

Voor Nederland kan je bij het CBS of Kadaster wel terecht voor grenzen. Internationaal zijn er ook wel allerlei bronnen, of de lokale geografische dienst of globale datasets zoals bijvoorbeeld Natural Earth of GADM.

Als snelheid een probleem is kan je naast een uitsnede, en versimpelen zoals Bram al aangeeft, ook de Matplotlib collectie opslaan met de pickle module. Het geeft ook wat nadelen in de zin dat je die bestanden moet beheren, maar het voordeel is dat je het Python object opslaat, en dus niet elke keer het process "Shapefile -> Matplotlib collectie" door hoeft. Ik vermijd zoiets liever, maar als je bijvoorbeeld elke dag duizenden afbeelding moet maken helpt zoiets zeker.

 

 


Python interpolatie: wat doe ik verkeerd?   ( 1261)
Christian (Oosterhout gld) -- 27-12-2017 21:58
@Rutger (Meppel) en eigenlijk ook @mods   ( 823)
Christian (Oosterhout gld) -- 27-12-2017 22:48
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 734)
Saskia (Diepenveen) ( 5m) -- 27-12-2017 23:09
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 787)
Rutger (Meppel) -- 27-12-2017 23:38
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 699)
Christian (Oosterhout gld) -- 28-12-2017 02:44
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 657)
Bram (Oldebroek/Bilthoven) -- 28-12-2017 10:32
Dank 🙂, loop nu een beetje vast met Shapefiles   ( 745)
Christian (Oosterhout gld) -- 28-12-2017 19:09
Re: Dank 🙂, loop nu een beetje vast met Shapefiles   ( 618)
Christian (Oosterhout gld) -- 28-12-2017 22:43
Re: Dank 🙂, loop nu een beetje vast met Shapefiles   ( 537)
Bram (Oldebroek/Bilthoven) -- 30-12-2017 08:52
Re: Dank 🙂, loop nu een beetje vast met Shapefiles   ( 809)
Rutger (Meppel) -- 30-12-2017 09:26
Bram, Rutger, bedankt voor jullie hulp!   ( 780)
Christian (Oosterhout gld) -- 31-12-2017 14:48
Re: Bram, Rutger, bedankt voor jullie hulp!   ( 636)
Rutger (Meppel) -- 31-12-2017 16:07
Re: Bram, Rutger, bedankt voor jullie hulp!   ( 530)
Bram (Oldebroek/Bilthoven) -- 31-12-2017 16:53
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 726)
Christian (Oosterhout gld) -- 29-12-2017 00:46
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 671)
Christian (Oosterhout gld) -- 28-12-2017 02:42
Re: @Rutger (Meppel) en eigenlijk ook @mods   ( 654)
Saskia (Diepenveen) ( 5m) -- 28-12-2017 22:21