load the required packages
library(dplyr) # general data wrangling
library(raster) # raster data
library(sf) # spatial data handling
library(tmap) # thematic mapping
download and prepare the data
## extract boundaries and data from OpenStreetMap
library(osmdata)
## define a bounding box for area of interest
osm_bbox <- c(76.