country <- get_map("country")
tm_shape(country) + tm_polygons()
division <- get_map("division")
tm_shape(division) + tm_polygons()
bd_plot("district",type = "interactive")
## ℹ tmap mode set to "view".
bd_plot("upazila")
## ℹ tmap mode set to "plot".
Kurukpata_union_data <- get_bandarban_map(upazila = "Alikadam", union = "Kurukpata")
datatable(sf::st_drop_geometry(Kurukpata_union_data))
sylhet <- get_divisions(divisions = "Sylhet",level = "upazila")
# single division
ggplot(data = sylhet) +
geom_sf() +
xlab("")+ ylab("")+
theme_minimal()
sylhet_chittagong_dhaka <- get_divisions(divisions = c("Sylhet", "Chittagong", "Dhaka"),level = "upazila")
ggplot(data = sylhet_chittagong_dhaka) +
geom_sf() +
xlab("")+ ylab("")+
theme_minimal()
population <- bangladeshPLUS::pop_district_2011[, c("district", "population")]
district <- get_map("district")
map_data <- dplyr::left_join(district, population, by = c("District" = "district"))
ggplot(data = map_data) +
geom_sf(aes(fill = population))+
theme_void()+
viridis::scale_fill_viridis(trans = "log", name="Population", labels = scales::unit_format(unit = "M", scale = 1e-6)) +
labs(
title = "Bangladesh Population Map",
subtitle = "Population & Housing Census 2011",
caption = "Data Source: BBS"
)
amtali <- bd_search("amtali", level = "union", as.is = T, coordinates = T)
knitr::kable(amtali, format = "html")
|
Division
|
District
|
Upazila
|
Union
|
lat
|
lon
|
|
Barisal
|
Barguna
|
Amtali
|
Amtali
|
22.07556
|
90.24699
|
|
Chittagong
|
Rangamati
|
Baghai Chhari
|
Amtali
|
23.10559
|
92.18832
|
|
Dhaka
|
Gopalganj
|
Kotali Para
|
Amtali
|
22.98264
|
90.03070
|
|
Chittagong
|
Khagrachhari
|
Matiranga
|
Amtali
|
23.11701
|
91.88545
|