# NOT RUN {
library(sf)
communes <- gisco_get_communes(country = "Belgique")
plot(
st_geometry(communes),
col = c("black", "gold", "brown2"),
border = "grey90"
)
title(
main = "Communes on Belgium (2016)",
sub = gisco_attributions(copyright = FALSE),
line = 1,
cex.sub = 0.8
)
# }
Run the code above in your browser using DataLab