Opvragen enkelbestemming ruimtelijke plannen bedrijf WMS

Hi!
Ik probeer voor een project een dataset/API te maken waarin alle gebieden staan waarbij enkelbestemming bedrijf is. Ik wil het liefst kunnen zien of er in de straal van een aantal kilometer van een coördinaat gebied is met deze enkelbestemming. Ik ben niet echt thuis in deze wereld en loop dus tegen een paar dingen aan:

  1. Ik gebruik nu de WMS-interface om data op te vragen. Ik gebruik hiervoor de link die hier staat: Geo services - PDOK. Als ik heb geprobeerd om als eerste test een jpeg op te vragen, maar ik krijg alleen wit terug:
from owslib.wms import WebMapService
wms = WebMapService("https://geodata.nationaalgeoregister.nl/plu/wms")
print(list(wms.contents)) # Lijst met onder andere de LAYER: Enkelbestemming
print(wms["Enkelbestemming"].queryable) # 1
print(wms["Enkelbestemming"].boundingBox) # BoundingBox (-14247.901, -10814.1739999994, 278026.09, 897413.893000001, 'EPSG:28992')
print(wms["Enkelbestemming"].boundingBoxWGS84) # BoundingBox (2.6760616721056922, 47.874168536166074, 7.360312726930821, 56.05651865088286)
print(wms["Enkelbestemming"].crsOptions) # Opties: ['EPSG:28992', 'EPSG:4326', 'EPSG:25832', 'EPSG:3035', 'EPSG:3034', 'EPSG:25831', 'EPSG:3857', 'EPSG:4258']
print(wms["Enkelbestemming"].styles) # Opties: Enkelbestemming met een legenda
print(wms.getOperationByName('GetMap').formatOptions) # een hele hoop format options, ik heb gekozen voor image/jpeg


img = wms.getmap(   layers=["Enkelbestemming"],
                    styles=["Enkelbestemming"],
                    srs='EPSG:28992',
                    bbox=[-14247.901, -10814.1739999994, 278026.09, 897413.893000001], # het hele gebied volgends de bounding box
                    size=(800, 800),
                    format='image/jpeg',
                    )

# Opslaan van plaatje
out = open('image.jpeg', 'wb')
out.write(img.read())
out.close()

  1. Ik ben niet erg thuis in deze wereld, maar wil er graag meer over leren. Is dit de handigste manier om dit probleem aan te pakken, of bestaat er al een API/dataset waarin je dit soort data kan opvragen en bijvoorbeeld data-analyse op kan doen?

Sorry als ik een hele gekke vraag stel. Zoals ik al zei, ben ik niet erg thuis in deze wereld.

Allereerst: je probeert nu een wms plaatje van 800 x 800 pixels te maken van dit gebied. Dat gebied is “te ver uitgezoomd” voor deze wms. Dat is waarom je een wit plaatje krijgt.

Dus zonder de code uit te proberen, dit is mogelijk een manier om de wms te bevragen, maar wat je daarvoor nog nodig hebt is een verder ingezoomde bounding box. Dus een van een paar kilometer rond een punt waarin je geïnteresseerd bent.

En de tweede vraag: https://www.pdok.nl/downloads/-/article/ruimtelijke-plannen je kan ook de atom downoaden. Bijvoorbeeld van de enkelbestemming. Daar kan je analyses op doen, als het je niet om figuren gaat.

Bedankt! Ik ga de data downloaden en daar mee aan de slag!

Graag wijs ik jou ook op de DSO-API’s, in dit geval Ruimtelijke plannen opvragen - Ontwikkelaarsportaal.
Middels deze REST-API worden alle ruimtelijke plannen (IMRO) ontsloten welke ook na in werking treden van de Omgevingswet (1 januari 2024) juridische werking hebben. Deze services zijn gebaseerd op de actuele planvoorraad van Ruimtelijkeplannen.nl, en zijn dus altijd actueel. De services en atomfeed van PDOK voor de ruimtelijke plannen worden maandelijks geactualiseerd.

en mocht je meer willen weten over de DSO-API’s dan attendeer ik jou graag op dit bericht;

2 likes