Vẽ bản đồ địa hình 3D với rasterVis

Gọi thư viện

library(rasterVis)
library(rgl)
library(raster)

Vẽ bản đồ 3D

Mở tập tin bản đồ địa hình

Các bạn có thể tải tập tin DEM từ địa chỉ sau: http://marcchoisy.free.fr/srtm90.tif. Hoặc từ gói srtmVN của tác giả Marc Choisy. Nguồn bản đồ CGIAR.

Đây là bản đồ về cao độ của Việt Nam phần đất liền.

DEM_VN <- raster("map/srtm90.tif")

Vẽ bản đồ 3D

Để vẽ bản đồ, ta sử dụng lệnh plot3D() của gói rasterVis.

plot3D(DEM_VN, zfac = 0.2, col = rainbow)

Vẽ bản đồ 3D cho 1 tỉnh

Ví dụ cho tỉnh Lâm Đồng. Đầu tiên ta cần tải ranh giới tỉnh Lâm Đồng với gói gadmVN.

lamdong_shp <- subset(gadmVN::gadm(resolution = "high"), province == "Lam Dong")

Cắt raster DEM_VN theo ranh giới tỉnh Lâm Đồng

lamdong_dem <- crop(DEM_VN, lamdong_shp)
lamdong_dem <- mask(lamdong_dem, lamdong_shp)

Vẽ bản đồ

plot3D(lamdong_dem, zfac = 0.5, col = rainbow)