開始使用免費開始

合併不同 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

編輯並執行程式碼