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])
