读取栅格文件
在 R 中读取栅格文件，最简单的方法是使用 raster 包中的 raster() 函数。您只需将栅格文件的文件名（包含扩展名）作为第一个参数 x 传入即可。
raster() 函数会根据文件名中的扩展名，优先使用 raster 包的原生函数读取特定类型的文件；否则会将读取工作交给 rgdal 包中的 readGDAL()。不直接使用 readGDAL() 的好处是，raster() 会返回一个 RasterLayer 对象。
一种常见的栅格文件类型是 GeoTIFF，其扩展名为 .tif 或 .tiff。我们已经从 US census 下载了一份家庭收入中位数的栅格，并放在您的工作目录中。
让我们来查看并将其读入。
本练习是课程的一部分
在 R 中可视化地理空间数据
练习说明
- 使用
dir()查看您的工作目录。 - 再次使用
dir()查看nyc_grid_data目录中的内容。 - 使用
raster()，传入该.tif文件的完整路径，将家庭收入中位数的栅格读入变量income_grid。 - 使用
summary()验证该对象为RasterLayer。 - 使用
plot()查看并验证栅格的内容。
交互式实操练习
通过完成这段示例代码来试试这个练习。
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