在 R 中访问栅格图层的像元答案

【问题标题】：Access cells of a raster layer in R在 R 中访问栅格图层的像元
【发布时间】：2020-12-13 04:29:09
【问题描述】：

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

library(sp)
library(raster)

我有一个栅格图层表示某个地理区域的气候变量，我有一些向量表示该区域中的点。

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

任何帮助将不胜感激！

编辑：点是在坐标中给出的（即每个点都由一个带有该点坐标的向量组成）。所有坐标都存储在一个 CSV 文件中，我已使用 library(sp) 中的函数 SpatialPoint 将其转换为 SpatialPoint 类。

【问题讨论】：

【解决方案1】：

这是一个最小的、独立的、可重现的示例：

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

要为r 获取xy 的单元格编号：

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

【讨论】：