# NOT RUN {
library(sf)
GreatCities <-
gisco_get_urban_audit(level = "GREATER_CITIES")
unique(GreatCities$CNTR_CODE)
countries <- gisco_get_countries()
plot(
st_geometry(countries) ,
col = "grey10",
xlim = c(-15, 25),
ylim = c(35, 60)
)
plot(st_geometry(GreatCities),
add = TRUE,
col = "darkgoldenrod3")
box()
title(
main = "Greater Cities on Europe (2020)",
sub = gisco_attributions(copyright = FALSE),
cex.main = 0.8,
cex.sub = 0.7,
line = 1
)
# }
Run the code above in your browser using DataLab