讀取 raster 檔案
在 R 中讀取 raster 檔案最簡單的方法,是使用 raster 套件的 raster() 函式。你只要把 raster 的檔名(包含副檔名)當作第一個引數 x 傳入即可。
raster() 會依據檔名的副檔名,對部分檔案型別使用 raster 套件內建的讀取函式;否則就交由 rgdal 套件的 readGDAL() 來讀取。之所以不直接用 readGDAL(),是因為 raster() 會回傳 RasterLayer 物件,較為方便。
常見的 raster 檔案格式是 GeoTIFF,副檔名為 .tif 或 .tiff。我們已從 US census 下載了中位收入的 raster,並放在你的工作目錄中。
現在就來看看並把它讀進來吧。
本練習屬於課程
在 R 視覺化地理空間資料
練習說明
- 使用
dir()檢視你的工作目錄內容。 - 再次使用
dir()檢視nyc_grid_data資料夾中的內容。 - 使用
raster(),將中位收入的 raster 讀入變數income_grid,並傳入該.tif檔案的完整路徑。 - 使用
summary()確認該 raster 以RasterLayer儲存。 - 使用
plot()檢視 raster 的內容。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
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