library(bangladesh)

country <- get_map("country")
division <- get_map("division")
district <- get_map("district")
upazila <- get_map("upazila")
union <- get_map("union")
plot(division[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
Khulna_division <- division %>%
  filter(Division %in% c("Khulna"))

plot(Khulna_division[1])

plot(district[1])

library(dplyr)

Kushtia_district <- district %>%
  filter(District %in% c("Kushtia"))

plot(Kushtia_district[1])

plot(upazila[1])

library(dplyr)
subset_upazila <- upazila %>%
  filter(Upazila %in% c("Daulatpur", "Bheramara", "Khoksa","Kumarkhali", "Kushtia", "Mirpur"))

plot(subset_upazila[1])