Bản đồ Việt Nam bằng phần mềm R (phần 2)

2. Kỹ thuật vẽ:

a) Vẽ bản đồ tỉnh Cà Mau (mstinh 18).

# Muốn vẽ bản đồ cần package sp

>Library(sp)

# Gọi tập tin VNM_adm1.rds trong ổ đĩa C, tên du1 trong R (không phân chia theo huyện)

> du1=readRDS("c:/VNM_adm1.rds")

> plot(du1[c(18),],col="orange",main=" Tỉnh Cà Mau")

Kết quả

Muốn phân chia theo ranh giới các huyện ta gọi tập tin VNM_adm2.rds. Nhưng phải xem mã số các huyện trong tập tin  VNM_adm2.csv (mã số các huyện từ 154 đến 162)

# Gọi tập tin VNM_adm2.rds trong ổ đĩa C vào R tên là du2

>du2=readRDS("c:/VNM_adm2.rds")

>plot(du2[c(154,155,156,157,158,159,160,161,162),],col="red",main=" Tỉnh Cà Mau")

Kết quả

b) Bản đồ các tỉnh đồng bằng sông Cửu Long

>plot(du1[c(3,7,13,17,18,20,31,33,34,40,53,60,61,63),],col="blue",main="Các tỉnh ĐBSLL")

Kết quả

c) Nước Việt Nam

Chỉ đơn giản lệnh sau:

> plot(du1[c(1:64),])

Đôi khi trong Hội nghị cần nhấn mạnh đến một tỉnh/huyện nào đó về một chủ đề cụ thể, ta có thể cho màu tỉnh/huyện đó để trực quan hơn. Việc làm này chỉ cần copy bản đồ sang phần mềm Untitled Paint (có sẳn trong máy) và tô màu theo ý.

Ví dụ: tỉnh Cà Mau màu đỏ

Kết quả

d) Bản đồ huyện xã (huyện Cái Nước và TT Cái Nước tỉnh Cà Mau).

Cần tới tập tin VNM_adm3.rds.

Mở tập tin VNM_adm3.csv ta thấy huyện Cái Nước, tỉnh Cà Mau có các mã số xã, thị trấn từ 2083,2084,2085,2086,2087,2088,2089,2090,2091,2092, 2093 ta làm như sau:

>plot(du3[c(2083,2084,2085,2086,2087,2088,2089,2090,2091,2092,2093),],col="red",main="Huyện Cái Nước - Cà Mau")

Cố nhiên có thể vẽ được bản đồ cho từng xã, thị trấn.

Vẽ TT Cái Nước

>plot(du3[c(2085),],col="3",main="TT Cái Nước, huyện Cái Nước - Cà Mau")

Kết quả

Tóm lại: để vẽ được bản đồ Việt Nam theo ý muốn ta chỉ cần thực hiện các nội dung sau:

1. Gọi package sp vận hành trong môi trường R.

2. Tải dữ liệu VNM_adm0.rds; VNM_adm1.rds; VNM_adm2.rds; VNM_adm3.rds từ website http://gadm.org/ về máy để vẽ bản đồ.

3. Tải dữ liệu VNM_adm0.csv; VNM_adm1.csv; VNM_adm2.csv; VNM_adm3.csv từ website http://gadm.org/ về máy để tìm mã số tỉnh, huyện, xã/thị trấn để vẽ bản đồ theo vị trí mình mong muốn.

BÙI VĂN DỦ