Đọc một shapefile
Shapefile là một trong những cách phổ biến nhất để chia sẻ dữ liệu không gian và có thể đọc vào R dễ dàng bằng readOGR() từ gói rgdal. readOGR() có hai đối số quan trọng: dsn và layer. Chính xác bạn truyền gì cho các đối số này phụ thuộc vào loại dữ liệu bạn đang đọc. Bạn đã học trong video rằng với shapefile, dsn nên là đường dẫn đến thư mục chứa các tệp cấu thành shapefile và layer là tên tệp của shapefile cụ thể (không kèm phần mở rộng).
Với bản đồ của bạn, bạn muốn ranh giới khu phố. Chúng tôi đã tải về Neighborhood Tabulation Areas, do Thành phố New York định nghĩa, từ Open Data Platform of the Department of City Planning. Tệp tải về ở dạng kho lưu trữ zip và chúng tôi đã giải nén và đặt kết quả vào thư mục làm việc của bạn.
Bạn sẽ dùng hàm dir() từ base R để xem nội dung thư mục làm việc, rồi đọc shapefile vào R.
Bài tập này là một phần của khóa học
Trực quan hóa dữ liệu địa lý trong R
Hướng dẫn bài tập
- Dùng
dir()không truyền đối số để tìm tên thư mục chứa shapefile. - Dùng
dir(), truyền vào đường dẫn đến thư mục shapefile, để xem các tệp bên trong. - Giờ bạn đã biết thư mục và tên tệp. Dùng
readOGR()để đọc shapefile khu phố vào một đối tượng tên làneighborhoods. - Kiểm tra nội dung bằng cách gọi
summary()trênneighborhoods. - Kiểm tra nội dung bằng cách vẽ
neighborhoods.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
library(sp)
library(rgdal)
# Use dir() to find directory name
# Call dir() with directory name
# Read in shapefile with readOGR(): neighborhoods
# summary() of neighborhoods
# Plot neighborhoods