Skip to content

Instantly share code, notes, and snippets.

@dieghernan
Created March 16, 2022 09:48
Show Gist options
  • Save dieghernan/d257a7690795077e1182db38b0f49536 to your computer and use it in GitHub Desktop.
Save dieghernan/d257a7690795077e1182db38b0f49536 to your computer and use it in GitHub Desktop.
library(ggplot2)
library(sf)
download.file(paste0("https://inspire.navarra.es/services/CP/wfs?",
"&request=Getfeature",
"&typenames=CP:CadastralParcel",
"&SRSname=EPSG:25830",
"&bbox=609758.580940066,4740059.01117674,611758.580940066,4742059.01117674"),
"cp.gml")
cp <- sf::st_read("cp.gml")
ggplot(cp) +
geom_sf()
# Buildings
download.file(paste0("https://inspire.navarra.es/services/BU/wfs?",
"&request=Getfeature",
"&typenames=BU:Building",
"&SRSname=EPSG:25830",
"&bbox=609758.580940066,4740059.01117674,611758.580940066,4742059.01117674"),
"bu.gml")
bu <- sf::st_read("bu.gml")
ggplot(bu) +
geom_sf()
# Addresses
download.file(paste0("https://inspire.navarra.es/services/AD/wfs?",
"&request=Getfeature",
"&typenames=AD:Address",
"&SRSname=EPSG:25830",
"&bbox=609758.580940066,4740059.01117674,611758.580940066,4742059.01117674"),
"ad.gml")
ad <- sf::st_read("ad.gml")
ggplot(ad) +
geom_sf()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment