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