Đọc một tệp raster
Các tệp raster được đọc vào R dễ nhất bằng hàm raster() từ gói raster. Bạn chỉ cần truyền tên tệp (kèm phần mở rộng) của raster làm đối số đầu tiên x.
Hàm raster() sử dụng một số hàm gốc của gói raster để đọc các loại tệp nhất định (dựa trên phần mở rộng trong tên tệp), và nếu không thì sẽ chuyển việc đọc tệp cho readGDAL() từ gói rgdal. Lợi ích của việc không gọi trực tiếp readGDAL() là raster() trả về một đối tượng RasterLayer.
Một loại tệp raster phổ biến là GeoTIFF, với phần mở rộng .tif hoặc .tiff. Chúng tôi đã tải về một raster thu nhập trung vị từ US census và đặt nó trong thư mục làm việc của bạn.
Hãy xem và đọc nó 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()để xem bên trong thư mục làm việc của bạn. - Dùng
dir()lần nữa để xem bên trong thư mụcnyc_grid_data. - Dùng
raster()để đọc raster thu nhập trung vị vào biếnincome_gridbằng cách truyền đường dẫn đầy đủ tới tệp.tif. - Dùng
summary()để kiểm tra raster được lưu dưới dạngRasterLayer. - Dùng
plot()để kiểm tra nội dung của raster.
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(raster)
# Call dir()
# Call dir() on the directory
# Use raster() with file path: income_grid
# Call summary() on income_grid
# Call plot() on income_grid