访问 R 中栅格图层的单元格

Access cells of a raster layer in R

我对空间分析完全陌生,我对 R 中的栅格图层有疑问。 我正在使用以下库:

library(sp)
library(raster)

我有一个栅格层代表某个地理区域的气候变量,我有一些矢量代表这个区域的点。

如何检查两个点是否在栅格图层的同一个像元中? 我知道如何提取每个点的气候值,但我不知道如何获取包含某个点的单元格。

非常感谢任何帮助!

编辑:这些点以坐标形式给出(即每个点都由一个带有该点坐标的向量组成)。所有坐标都存储在 CSV 文件中,我使用 library(sp).

中的函数 SpatialPoint 将其转换为 SpatialPoint class

这是一个最小的、self-contained、可重现的例子:

library(raster)
r <- raster()
xy <- cbind(1:5, 1:5)

r获取xy的手机号:[=​​14=]

cellFromXY(r, xy)
[1] 32222 31863 31504 31145 30786