Geography assignment 02

#Data

library(bangladesh)

country <- get_map("country")
division <- get_map("division")
district <- get_map("district")
upazila <- get_map("upazila")
union <- get_map("union")

plot(country[1])

plot(division[1])

plot(district[1])

plot(upazila[1])

plot(union[1])

Division Map

plot(division[1], main = "Division Map of Bangladesh")

library(sf)

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

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

chittagong <- division %>%
  filter(Division =="Chittagong")

plot(chittagong[1])

rangamati_district <- district %>%
  filter(Division =="Chittagong")

plot(rangamati_district[1])

chittagong_district <- district %>%
  filter(District =="Rangamati")

plot(chittagong_district[1])

rangamati_upazila <- upazila %>% 
  filter(District == "Rangamati")

plot(rangamati_upazila[1])

rangamati_upazila <- upazila %>%
  filter(Upazila %in% c("Barkal", "Langadu"))

plot(rangamati_upazila[1])