The World Database on Protected Areas (WDPA) là một cơ sở dữ liệu mở về các khu bảo tồn trên toàn cầu. WDPA là một dự án hợp tác giữa UN Environment Programme và International Union for Conservation of Nature (IUCN) và được quản lý bởi UN Environment Programme World Conservation Monitoring Centre (UNEP-WCMC). Dữ liệu WDPA được cập nhật hàng tháng và có thể được tải trực tiếp từ trang web hoặc qua thư viện wdpar trên R.
Nếu tải dữ liệu trực tiếp từ trang web thì dữ liệu sẽ có phạm vi toàn cầu với dung lượng rất lớn, còn nếu tải dữ liệu qua wdpar thì có thể tải dữ liệu riêng lẻ của từng quốc gia. Trong bài viết này, mình sẽ chia sẻ cách dùng wdpar để tải dữ liệu khu bảo tồn ở Việt Nam.
Tải thư viện wdpar và các thư viện khác
library(wdpar) # tải dữ liệu khu bảo tồn
library(tidyverse) # xử lý và trực quan hóa dữ liệu
library(sf) # xử lý dữ liệu không gian dạng vector
library(stars) # xử lý dữ liệu không gian dạng raster
library(tmap) # hiển thị dữ liệu không gian
# webdriver::install_phantomjs() # có thể phải cài PhantomJS để chạy các lệnh của wdpar
# Dữ liệu khu bảo tồn Việt Nam
vietnam <- wdpa_fetch(
x = "Vietnam", # Tên nước hoặc mã nước theo ISO3 (ví dụ VNM)
wait = TRUE,
download_dir = "." # Thư mục lưu file, "." nghĩa là sẽ lưu ở thư mục làm việc hiện tại (working directory)
)
vietnam
Simple feature collection with 191 features and 28 fields
Geometry type: GEOMETRY
Dimension: XY
Bounding box: xmin: 102.1522 ymin: 8.429265 xmax: 109.4629 ymax: 23.17159
Geodetic CRS: WGS 84
# A tibble: 191 × 29
WDPAID WDPA_PID PA_DEF NAME ORIG_NAME DESIG DESIG_ENG DESIG_TYPE IUCN_CAT INT_CRIT MARINE REP_M_AREA REP_AREA
* <dbl> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <dbl> <dbl>
1 902735 902735 1 Bau Sau Bau Sau Rams… Ramsar S… Internati… Not Rep… (i)(ii)… 0 0 138.
2 7877 7877 1 Cat Ba Cát Bà Vườn… National… National II Not App… 1 0 174.
3 7878 7878 1 Cuc Phuong Cúc Phươ… Vườn… National… National II Not App… 0 0 224.
4 10110 10110 1 Bach Ma Bạch Mã Vườn… National… National II Not App… 0 0 375.
5 10111 10111 1 Con Dao Côn Đảo Vườn… National… National II Not App… 1 42.7 58.3
6 10113 10113 1 Ba Vi Vườn quố… Vườn… National… National II Not App… 0 0 97.0
7 10114 10114 1 Ngoc Trao Ngoc Trao Cult… Cultural… National V Not App… 0 0 3
8 10117 10117 1 Den Hung Đền Hùng Khu … Cultural… National V Not App… 0 0 5.38
9 10127 10127 1 Bac Son Bắc Sơn Khu … Nature R… National V Not App… 0 0 10.1
10 10187 10187 1 Ban dao Son… Bán đảo … Khu … Nature R… National Not Rep… Not App… 0 0 25.9
# ℹ 181 more rows
# ℹ 16 more variables: NO_TAKE <chr>, NO_TK_AREA <dbl>, STATUS <chr>, STATUS_YR <int>, GOV_TYPE <chr>, OWN_TYPE <chr>,
# MANG_AUTH <chr>, MANG_PLAN <chr>, VERIF <chr>, METADATAID <int>, SUB_LOC <chr>, PARENT_ISO <chr>, ISO3 <chr>,
# SUPP_INFO <chr>, CONS_OBJ <chr>, geometry <GEOMETRY [°]>
Hiển thị toàn bộ dữ liệu
tmap_mode(mode = "view") # thiết lập chế độ tương tác
tm_shape(vietnam) +
tm_polygons(id = "ORIG_NAME")