library(geofabrik) # Downloads OSM data from geofabrik
# customize this script
# Check all geofabrick links
dir.create("OSM")
head(urlgeo)
names(urlgeo)
# convert to data.table
# search africa
af <- grep("Africa", urlgeo$Subregion)
urlgeo[af, ]
if (download_osm) {
file_down <- paste0("OSM/", OSM_region)
get_osm(region = OSM_region, file = file_down)
unzip(
zipfile = paste0(gsub(" ", "", file_down),
".shp.zip"),
exdir = "network"
)
# This approach may be faster with QGIS
# let us read
roads <- st_read("network/gis_osm_roads_free_1.shp")
# let us select the streets that we need
st <- c(
"motorway", "motorway_link", "trunk", "trunk_link",
"primary", "primary_link", "secondary", "secondary_link",
"tertiary", "tertiary_link"
)
roads <- roads[roads$fclass %in% st, ]
# save roads on network
saveRDS(roads, "network/roads.rds")
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.