# NOT RUN {
library(sf)
NL <- gisco_get_countries(country = "NL")
AirP_NL <- gisco_get_airports(country = "NL")
Ports <- gisco_get_ports()
# Intersect with NL
PortsNL <- st_intersection(Ports, NL)
plot(st_geometry(NL), col = "wheat")
plot(
st_geometry(PortsNL),
pch = 22,
col = "forestgreen",
add = TRUE,
cex = 0.8
)
plot(
st_geometry(AirP_NL),
pch = 20,
col = "steelblue",
add = TRUE,
cex = 1.2
)
legend(
"topright",
legend = c("Port", "Airport"),
col = c("forestgreen", "steelblue"),
cex = 0.9,
bty = "n",
pch = c(22, 20),
pt.cex = c(1, 1.5),
y.intersp = 2
)
title(
main = "Transport Network on the Nethelands",
sub = gisco_attributions(),
line = 1,
cex.sub = 0.7,
font.sub = 3
)
# }
Run the code above in your browser using DataLab