Gộp dữ liệu từ các CRS/phép chiếu khác nhau
Mỗi đối tượng không gian đều có một hệ quy chiếu tọa độ (CRS) đi kèm. Thông thường, điều này được thiết lập khi nhập dữ liệu và sẽ được đọc trực tiếp từ các tệp không gian. Đây chính là cách neighborhoods và nyc_tracts nhận thông tin hệ tọa độ của chúng.
Cả hai gói sp và raster đều có hàm proj4string() trả về CRS của đối tượng được gọi.
Cố gắng làm việc với dữ liệu không gian dùng các CRS khác nhau cũng giống như làm việc với một tập dữ liệu dùng dặm và tập khác dùng kilômét. Chúng đo cùng một thứ, nhưng các con số không thể so sánh trực tiếp.
Hãy cùng xem hai đối tượng đa giác của chúng ta.
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
- Gọi
proj4string()trênneighborhoods, rồi gọi lại trênnyc_tracts. Xác nhận rằng hai chuỗi là khác nhau. - Xem
head()củacoordinates()choneighborhoodsvà lặp lại chonyc_tracts. Bạn có thấy vấn đề không?nyc_tractscó tọa độ x khoảng -70, cònneighborhoodskhoảng 1.000.000! - Vẽ
neighborhoods, sau đó vẽnyc_tractsvớicol = "red"vàadd = TRUEđể chồng lên trên.
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)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts