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

问题描述

我是空间分析的新手,我对R中的栅格图层有疑问。 我正在使用以下库:

library(sp)
library(raster)

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

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

任何帮助将不胜感激!

编辑:这些点以坐标给出(即,每个点都由具有该点坐标的向量组成)。所有坐标都存储在CSV文件中,我已使用SpatialPoint中的函数SpatialPoint将其转换为library(sp)类。

解决方法

在询问R问题时,请提供一个最小的,独立的,可复制的示例。像这样:

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

要获取xy的{​​{1}}的单元格编号:

r