r - 从 Raster 对象获取 XYZ 列表

Question

我知道要从 XYZ 列表(其中 x 和 y 是坐标的数字向量，z 是数字矩阵)获取栅格，我可以直接使用函数 raster。但是如果我有一个光栅并且我想获得一个 XYZ 列表怎么办？ raster 包上是否已经开发了一个简单的功能？

当然，我知道我可以在栅格对象上建立一个索引不同层次的列表，但我发现我必须做一些逆向来实现它，所以我想知道是否有一种简单有效的方法.

谢谢。

Answer 1

您可以考虑使用 rasterToPoints:

library(raster)
r <- raster()
values(r) <- 1:ncell(r)

xyz <- rasterToPoints(r)

head(xyz)
           x    y layer
#[1,] -179.5 89.5     1
#[2,] -178.5 89.5     2
#[3,] -177.5 89.5     3
#[4,] -176.5 89.5     4
#[5,] -175.5 89.5     5
#[6,] -174.5 89.5     6

但请注意，该函数只返回不是 NA 的单元格的值并且这会返回一个矩阵；不是 list 。但这可能更实用。否则请参阅@Val 的解决方案...