Vietnam_Map
hoangphuluu
9/14/2020
R Markdown
This is an R practice for illustrating map document:
#Get Vietnam geography by provinces:
vietnam <- getData("GADM", country = "Vietnam", level = 1)
head(vietnam)## GID_0 NAME_0 GID_1 NAME_1 VARNAME_1 NL_NAME_1
## 1 VNM Vietnam VNM.1_1 An Giang An Giang <NA>
## 12 VNM Vietnam VNM.2_1 B<U+1EA1>c Liêu Bac Lieu <NA>
## 23 VNM Vietnam VNM.3_1 B<U+1EAF>c Giang Bac Giang <NA>
## 34 VNM Vietnam VNM.4_1 B<U+1EAF>c K<U+1EA1>n Bac Kan <NA>
## 45 VNM Vietnam VNM.5_1 B<U+1EAF>c Ninh Bac Ninh <NA>
## 56 VNM Vietnam VNM.6_1 B<U+1EBF>n Tre Ben Tre <NA>
## TYPE_1 ENGTYPE_1 CC_1 HASC_1
## 1 T<U+1EC9>nh Province <NA> VN.AG
## 12 T<U+1EC9>nh Province <NA> VN.BL
## 23 T<U+1EC9>nh Province <NA> VN.BG
## 34 T<U+1EC9>nh Province <NA> VN.BK
## 45 T<U+1EC9>nh Province <NA> VN.BN
## 56 T<U+1EC9>nh Province <NA> VN.BRnames(vietnam)## [1] "GID_0" "NAME_0" "GID_1" "NAME_1" "VARNAME_1"
## [6] "NL_NAME_1" "TYPE_1" "ENGTYPE_1" "CC_1" "HASC_1"vietnam_df <- vietnam %>% fortify()## Regions defined for each Polygonshead(vietnam_df)## long lat order hole piece id group
## 1 105.3745 10.24604 1 FALSE 1 1 1.1
## 2 105.3362 10.23442 2 FALSE 1 1 1.1
## 3 105.3154 10.26842 3 FALSE 1 1 1.1
## 4 105.3120 10.27393 4 FALSE 1 1 1.1
## 5 105.3065 10.26894 5 FALSE 1 1 1.1
## 6 105.3013 10.27698 6 FALSE 1 1 1.1unique(vietnam_df$piece)## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
## [18] 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34
## [35] 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
## [52] 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68
## [69] 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85
## [86] 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102
## [103] 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119
## [120] 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136
## [137] 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153
## [154] 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170
## [171] 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187
## [188] 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204
## [205] 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221
## [222] 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238
## [239] 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255
## [256] 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272
## [273] 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289
## [290] 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306
## [307] 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323
## [324] 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340
## [341] 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357
## [358] 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374
## [375] 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391
## [392] 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408
## [409] 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425
## [426] 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442
## [443] 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459
## [460] 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476
## [477] 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493
## [494] 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510
## [511] 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527
## [528] 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544
## [545] 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561
## [562] 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578
## [579] 579 580 581
## 581 Levels: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 ... 581theme_set(theme_minimal())
m1 <- vietnam_df %>%
ggplot(aes(x = long, y = lat, group = group))+
geom_polygon(aes(fill = id), color = "grey30", show.legend = FALSE) +
labs(x = NULL, y = NULL)
m1
k1 <- ggplot() +
geom_polygon(data = vietnam_df,
aes(long, lat, group = group, fill = id),
show.legend = FALSE, color = "grey")
k1 + coord_map("albers", lat0 = 30, lat1 = 40)
#Province_names
vn_province <- vietnam$NAME_1
vn_province %>% head()## [1] "An Giang" "B<U+1EA1>c Liêu" "B<U+1EAF>c Giang" "B<U+1EAF>c K<U+1EA1>n" "B<U+1EAF>c Ninh" "B<U+1EBF>n Tre"#Simulating the income per households by provinces:
set.seed(1)
mydf <- data.frame(Tinh = vn_province, id = 1:length(vn_province),
share = runif(length(vn_province), 20, 60))
#Join data:
vietnam_df_income <- merge(vietnam_df, mydf, by = "id", all.x = TRUE)
m2 <- vietnam_df_income %>%
ggplot(aes(x = long, y = lat, group = group))+
geom_polygon(aes(fill = share), color = "grey30") +
labs(x = NULL, y = NULL)
m2
library(viridis)## Loading required package: viridisLitem2 + scale_fill_viridis(direction = -1, option = "D", "income")
#Get data map by districts, then level = 2
#Get data map by communes, then level = 3
vietnam_by_districts <- getData("GADM", country = "Vietnam", level = 2)
vietnam_by_districts %>% head()## GID_0 NAME_0 GID_1 NAME_1 NL_NAME_1 GID_2 NAME_2
## 1 VNM Vietnam VNM.1_1 An Giang <NA> VNM.1.1_1 An Phú
## 4 VNM Vietnam VNM.1_1 An Giang <NA> VNM.1.2_1 Ch<U+1EE3> M<U+1EDB>i
## 5 VNM Vietnam VNM.1_1 An Giang <NA> VNM.1.3_1 Châu Ð<U+1ED1>c
## 6 VNM Vietnam VNM.1_1 An Giang <NA> VNM.1.4_1 Châu Phú
## 7 VNM Vietnam VNM.1_1 An Giang <NA> VNM.1.5_1 Châu Thành
## 8 VNM Vietnam VNM.1_1 An Giang <NA> VNM.1.6_1 Long Xuyên
## VARNAME_2 NL_NAME_2 TYPE_2 ENGTYPE_2 CC_2 HASC_2
## 1 An Phu <NA> Huy<U+1EC7>n District <NA> VN.TT.AL
## 4 Cho Moi <NA> Huy<U+1EC7>n District <NA> VN.BD.AL
## 5 Chau Doc <NA> Thành ph<U+1ED1> City <NA> VN.KG.AB
## 6 Chau Phu <NA> Huy<U+1EC7>n District <NA> VN.HP.AD
## 7 Chau Thanh <NA> Huy<U+1EC7>n District <NA> VN.GL.AK
## 8 Long Xuyen <NA> Thành ph<U+1ED1> City <NA> VN.HP.ALunique(vietnam_by_districts$NAME_1)## [1] "An Giang" "B<U+1EA1>c Liêu" "B<U+1EAF>c Giang"
## [4] "B<U+1EAF>c K<U+1EA1>n" "B<U+1EAF>c Ninh" "B<U+1EBF>n Tre"
## [7] "Bà R<U+1ECB>a - Vung Tàu" "Bình Ð<U+1ECB>nh" "Bình Duong"
## [10] "Bình Phu<U+1EDB>c" "Bình Thu<U+1EAD>n" "C<U+1EA7>n Tho"
## [13] "Cà Mau" "Cao B<U+1EB1>ng" "Ð<U+1EAF>k L<U+1EAF>k"
## [16] "Ð<U+1EAF>k Nông" "Ð<U+1ED3>ng Nai" "Ð<U+1ED3>ng Tháp"
## [19] "Ðà N<U+1EB5>ng" "Ði<U+1EC7>n Biên" "Gia Lai"
## [22] "H<U+1EA3>i Duong" "H<U+1EA3>i Phòng" "H<U+1EAD>u Giang"
## [25] "H<U+1ED3> Chí Minh" "Hà Giang" "Hà N<U+1ED9>i"
## [28] "Hà Nam" "Hà Tinh" "Hoà Bình"
## [31] "Hung Yên" "Khánh Hòa" "Kiên Giang"
## [34] "Kon Tum" "L<U+1EA1>ng Son" "Lai Châu"
## [37] "Lâm Ð<U+1ED3>ng" "Lào Cai" "Long An"
## [40] "Nam Ð<U+1ECB>nh" "Ngh<U+1EC7> An" "Ninh Bình"
## [43] "Ninh Thu<U+1EAD>n" "Phú Th<U+1ECD>" "Phú Yên"
## [46] "Qu<U+1EA3>ng Bình" "Qu<U+1EA3>ng Nam" "Qu<U+1EA3>ng Ngãi"
## [49] "Qu<U+1EA3>ng Ninh" "Qu<U+1EA3>ng Tr<U+1ECB>" "Sóc Trang"
## [52] "Son La" "Tây Ninh" "Th<U+1EEB>a Thiên Hu<U+1EBF>"
## [55] "Thái Bình" "Thái Nguyên" "Thanh Hóa"
## [58] "Ti<U+1EC1>n Giang" "Trà Vinh" "Tuyên Quang"
## [61] "Vinh Long" "Vinh Phúc" "Yên Bái"# filter Hanoi :
hanoi <- vietnam_by_districts[vietnam_by_districts$NAME_1 == "Hà Nội", ]
# Transform Hanoi to a dataframe:
hanoi_df <- hanoi %>% fortify()## Regions defined for each Polygons# Hanoi map:
#m5 <- hanoi_df %>%
# ggplot(aes(x = long, y = lat, group = group))+
# geom_polygon(aes(fill = id), color = "grey30", show.legend = FALSE) +
# labs(x = NULL, y = NULL, title = "Districts of Hanoi",
# caption = "map")
#m5
#m5 + scale_fill_manual(values = c('#a6cee3','#fb9a99','#b2df8a','#33a02c','#fb9a99','#e31a1c',
#'#1f78b4','#ff7f00','#cab2d6','#6a3d9a','#ffff99','#b15928',
# '#33a02c','#fb9a99','#e31a1c', '#fdbf6f','#ff7f00','#cab2d6','#6a3d9a',
#'#6a3d9a','#1f78b4','#b2df8a','#33a02c','#fb9a99','#e31a1c','#fdbf6f',
#'#ff7f00','#cab2d6','#b15928','#ffff99','#b15928','#b2df8a'))