R Notebook

https://drive.google.com/file/d/1icDRT8jQ86A7Ya_5-n5vIuHapiqsjHYu/view?usp=sharing

# Install if not already installed
if (!require(leaflet)) install.packages("leaflet")

## Loading required package: leaflet

library(leaflet)
if (!require(sf)) install.packages("sf") 

## Loading required package: sf

## Linking to GEOS 3.13.1, GDAL 3.11.4, PROJ 9.7.0; sf_use_s2() is TRUE

library(sf)

# Google Drive file ID
id <- "1icDRT8jQ86A7Ya_5-n5vIuHapiqsjHYu"
# Download & unzip
zip <- tempfile(fileext = ".zip")
dir <- tempdir()
download.file(sprintf("https://drive.google.com/uc?export=download&id=%s", id), zip, mode = "wb")
unzip(zip, exdir = dir)
# Read shapefile
shp <- list.files(dir, "hotosm_bgd_health_facilities_points_shp\\.shp$", 
            full.names = TRUE, recursive = TRUE)
health_facilities <- st_read(shp[1])

## Reading layer `hotosm_bgd_health_facilities_points_shp' from data source 
##   `C:\Users\hp\AppData\Local\Temp\RtmpuglG08\hotosm_bgd_health_facilities_points_shp.shp' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 6929 features and 14 fields
## Geometry type: POINT
## Dimension:     XY
## Bounding box:  xmin: 88.10011 ymin: 20.63178 xmax: 92.58573 ymax: 26.41615
## Geodetic CRS:  WGS 84

names(health_facilities)

##  [1] "name"       "name_en"    "amenity"    "building"   "healthcare"
##  [6] "healthca_1" "operator_t" "capacity_p" "addr_full"  "addr_city" 
## [11] "source"     "name_bn"    "osm_id"     "osm_type"   "geometry"

leaflet(health_facilities) %>%  addTiles() %>%  
  addCircleMarkers(
    popup = ~name)