使用套件取得資料
從檔案讀取空間資料是把空間資料載入 R 的一種方式,但也有一些套件提供常用的空間資料。例如,rnaturalearth 套件提供來自 Natural Earth 的資料,這是一個高解析度世界地圖的來源,包含海岸線、州界與聚落地點。事實上,第 2 章所用的資料就是來自這裡。
你將要在 New York County(也就是曼哈頓區)的普查分區層級檢視家戶所得中位數。不過在此之前,你需要知道各普查分區的邊界。R 的 tigris 套件提供一種簡單的方法,能依照美國普查的地理分區下載並匯入 shapefile。你會使用 tracts() 函式下載分區邊界;此外,tigris 也提供 states()、counties()、places() 等多個函式,對應普查所定義的各種地理層級。
現在就來取得這些分區的空間資料吧。
本練習屬於課程
在 R 視覺化地理空間資料
練習說明
- 以
state = "NY"、county = "New York"、cb = TRUE呼叫tracts(),並將結果存成nyc_tracts。 - 對
nyc_tracts使用summary(),確認回傳物件是SpatialPolygonsDataFrame。 - 使用
plot()繪製nyc_tracts,檢查其內容。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
library(sp)
library(tigris)
# Call tracts(): nyc_tracts
# Call summary() on nyc_tracts
# Plot nyc_tracts