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:
- 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()
- 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.