ĐỒ THỊ THỊ TRƯỜNG

GDP CÁC TỈNH MIỀN TÂY NAM BỘ

VIỆT NAM

1. Kết nối dữ liệu

dulieu <-read.csv("http://solieu.vip/csv/mtay.csv")
head(dulieu)
##     ï..Tinh    GDP       lo        la id   PCI   FDI   HANG
## 1  An Giang  79083 105.1740 10.518005  1 58.44   1.8 Hang 2
## 2  Bac Lieu  32567 105.4841  9.226981 10 56.58 368.1 Hang 1
## 3   Ben Tre 123098 106.4613 10.045397 12 63.46 403.5 Hang 3
## 4    Ca Mau  34521 104.9583  8.856421 27 59.32  43.9 Hang 2
## 5   Can Tho 345892 105.4079 10.177783 17 56.51  37.7 Hang 1
## 6 Dong Thap  45672 105.6247 10.624494 33 66.42   7.3 Hang 3

2. Lấy bản đồ về

2.1 Bản đồ Việt Nam

library(raster)
## Loading required package: sp
vietnam <-getData("GADM", country="Vietnam", level=1)
library(ggplot2)
ggplot() + geom_polygon(data=vietnam,aes(x=long,y=lat,group=group,fill=id))
## Regions defined for each Polygons

2.2 Bản đồ miền tây

mientay <-vietnam[c(1,2,6,12,13,18,24,33,39,51,58,59,61),]
ggplot() + geom_polygon(data=mientay,aes(x=long,y=lat,group=group,fill=id))
## Regions defined for each Polygons

3. Chuyển mã

library(tidyverse)
## -- Attaching packages ------------------------------------------------------- tidyverse 1.3.0 --
## v tibble  2.1.3     v dplyr   0.8.3
## v tidyr   1.0.0     v stringr 1.4.0
## v readr   1.3.1     v forcats 0.4.0
## v purrr   0.3.3
## -- Conflicts ---------------------------------------------------------- tidyverse_conflicts() --
## x tidyr::extract() masks raster::extract()
## x dplyr::filter()  masks stats::filter()
## x dplyr::lag()     masks stats::lag()
## x dplyr::select()  masks raster::select()
mtay<-fortify(mientay)
## Regions defined for each Polygons

4. Kết nối dữ liệu

id <- c(1, 10, 12, 17, 27, 33, 4, 47, 5, 54, 55, 56,58)
dulieu2 <-data.frame(dulieu,id)
noi <- merge(mtay, dulieu2, by="id", all.x = TRUE)

5. Vẽ đồ thị

5.1 Vẽ bản đồ hành chính miền tây

m1 <-ggplot(data=noi,aes(x = long, y = lat,group = group))+ geom_polygon(aes(fill = id), color = "white") + labs(x = NULL, y =  NULL) + theme(legend.position = "none")
m1

5.2 Thêm địa điểm thủ phủ tỉnh

m2 <- m1 + geom_point(aes(x=lo, y=la), col="red")
m2

5.3 Thêm tên tỉnh

###{r} ###m3 <- m2 + geom_text(aes(x=lo, y=la, label= Tinh), vjust=- 1.2, col="green")## ###m3 ###

5.4 Thêm giá trị GDP cho tỉnh

###{r} ###m4 <- m3 + geom_text(aes(x=lo, y=la, label= noi$FDI), vjust= 1.2,col="#35c908") ###m4 ###

** Như vậy là chúng ta đã xong bản đồ GDP của các tỉnh miền tây rồi **