Extract national roads from Microsoft's dataset (Bing Maps) in R

microsoft_roads.R

library("sf")
library("readr")
library("geojsonsf")

## download all roads for Europe
data = "https://usaminedroads.blob.core.windows.net/road-detections/Europe-Full.zip"
options(timeout = 6000)
download.file(data, "Europe-Full.zip")
unzip("Europe-Full.zip")

## read data in chunks (and select only in Czech Republic)
f = function(x, pos) subset(x, country == "CZE")
df = read_tsv_chunked("Europe-Full.tsv", DataFrameCallback$new(f),
                      col_names = c("country", "geom"),
                      chunk_size = 100000, col_types = "fc")
geom = geojson_sfc(df$geom)
write_sf(geom, "CZE.gpkg") # save as geopackage