library(bangladesh)
district <- get_map("district")
upazila <- get_map("upazila")
plot(upazila[1])

plot(district[1])

library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
# Create subset for Kushtia
sd <- district %>%
  filter(District %in% c("Kushtia"))

plot(sd[1])

# Create subset for Jatrabari
subset_upazila <- upazila %>%
  filter(Upazila %in% c("Jatrabari"))

plot(subset_upazila[1])