# NOT RUN {
library(sf)
Base <- esp_get_munic(region = c("Castilla y Leon"))
SAN <-
esp_get_munic(
region = c("Castilla y Leon"),
munic = c("^San ", "^Santa ")
)
plot(st_geometry(Base), col = "cornsilk", border = "grey80")
plot(st_geometry(SAN),
col = "firebrick3",
border = NA,
add = TRUE)
# }
Run the code above in your browser using DataLab