Hoi Lennart,
Je zou eens naar onderstaande post kunnen kijken. Dit gaat niet over een straal, maar over een bounding box via WFS.
Hieronder staat de R code die ik daar als voorbeeld geef.
Tip: kijk ook eens naar de library httr
(van Hadley Wickham) met de functies parse_url
en build_url
. Hiermee kun je al je parameters in een overzichtelijke list opgeven.
HTH,
Egge-Jan
library(httr)
library(sf)
library(tmap)
url <- parse_url("https://geodata.nationaalgeoregister.nl/nwbwegen/wfs")
url$query <- list(service = "WFS",
version = "2.0.0",
request = "GetFeature",
typename = "nwbwegen:wegvakken",
BBOX = "154500,462500,155500,463500",
outputFormat = "application/json")
request <- build_url(url)
nwb_wegen_amersfoort <- st_read(request)
tmap_mode("view")
tm_shape(nwb_wegen_amersfoort)+tm_lines(col="red", lwd = 4)+tm_format("NLD")