栅格到缺少值的点

Question

我正在使用软件包rasterToPoints()中的函数raster。默认情况下，此函数忽略NA值。

raster <- raster(matrix(c(1,2,NA, 1,2,3, 1,2,3), nrow = 3, ncol = 3 ))
raster_points <- rasterToPoints(raster, na.rm = FALSE)

因此，我放置了一个具有200个像元的栅格，并且得到了只有150行的data.frame。如何在最终data.frame中获得带有空白单元格的点？

Answer 1

您可以结合使用xyFromCell，values和cbind。这将以矩阵形式返回每个像元的中心坐标（初始化栅格时，除非另有定义，否则它假定为1 x 1范围）及其值（包括NAs），以矩阵形式：

cbind(xyFromCell(raster, 1:ncell(raster)), values(raster))