合併不同 CRS/投影的資料
每個空間物件都會附帶一個座標參考系統(CRS)。一般來說,這會在匯入資料時設定,並會直接從空間檔案讀取。neighborhoods 和 nyc_tracts 的座標系統資訊就是這樣取得的。
sp 與 raster 這兩個套件都有 proj4string() 函式,可回傳被呼叫物件的 CRS。
嘗試用不同 CRS 的空間資料一起作業,有點像同時使用一個以英里為單位、另一個以公里為單位的資料集。它們量的是同一件事,但數值無法直接比較。
我們來看看這兩個多邊形物件。
本練習屬於課程
在 R 視覺化地理空間資料
練習說明
- 對
neighborhoods呼叫proj4string(),再對nyc_tracts呼叫一次。*確認兩個字串不同。* - 查看
neighborhoods的coordinates()的head(),並對nyc_tracts重複一次。*你看出問題了嗎?nyc_tracts的 x 座標大約在 -70,而neighborhoods則在 1,000,000 附近!* - 先繪製
neighborhoods,再以col = "red"與add = TRUE繪製nyc_tracts,把它們疊加在上面。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts