Denmark Map

library(osmextract)

## Data (c) OpenStreetMap contributors, ODbL 1.0. https://www.openstreetmap.org/copyright.
## Any product made from OpenStreetMap must cite OSM as the data source.
## Geofabrik data are taken from https://download.geofabrik.de/
## For usage details of bbbike data see https://download.bbbike.org/osm/
## OpenStreetMap_fr data are taken from http://download.openstreetmap.fr/

mypath <- "C://users//np83zg//OneDrive - Aalborg Universitet//Skrivebord"
# mydata_line <- oe_get(place="Denmark",download_directory=mypath)
# mydata_poly <- oe_get(place="Denmark", layer = "multipolygons",download_directory=mypath,stringsAsFactors = FALSE)

mypath <- "C://users//np83zg//OneDrive - Aalborg Universitet//Skrivebord//denmark"
setwd(mypath)
mydata_line <- oe_read(file="geofabrik_denmark-latest.gpkg")
mydata_poly <- oe_read(file="geofabrik_denmark-latest.gpkg",layer = "multipolygons",stringsAsFactors = FALSE)

library(tmap)
ht <- c("primary","secondary","primary_link","secondary_link")
osm_major_roads = mydata_line[mydata_line$highway %in% ht, ]
tm_shape(osm_major_roads) + tm_lines(col="grey29") 

## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1

index1 <- which(mydata_poly$admin_level==7)
index2 <- which(mydata_poly$landuse=="residential")
tmap_mode("plot")

## tmap mode set to plotting

tm_shape(mydata_poly[index1,"admin_level"]) + tm_fill("seashell3")+tm_borders("white", lwd = 1) + 
tm_shape(osm_major_roads[,"highway"])+ tm_lines(col="grey29") + 
tm_shape(mydata_poly[index2,"landuse"]) + tm_fill("black")