问题描述
我是空间分析的新手,我对R中的栅格图层有疑问。 我正在使用以下库:
library(sp)
library(raster)
我有一个栅格图层,表示某个地理区域的气候变量,我有一些矢量,表示该区域的点。
如何检查两个点是否在栅格图层的同一像元中? 我知道如何提取每个点的气候值,但是我不知道如何获取包含某个点的单元格。
任何帮助将不胜感激!
编辑:这些点以坐标给出(即,每个点都由具有该点坐标的向量组成)。所有坐标都存储在CSV文件中,我已使用SpatialPoint
中的函数SpatialPoint
将其转换为library(sp)
类。
解决方法
在询问R问题时,请提供一个最小的,独立的,可复制的示例。像这样:
library(raster)
r <- raster()
xy <- cbind(1:5,1:5)
要获取xy
的{{1}}的单元格编号:
r