与 R 中的值匹配的像素的栅格网格位置/坐标

Question

有没有办法提取与特定值匹配的像素的网格位置或(最好是具有明确范围的栅格)点/质心坐标？我几乎有一个非常低效的工作流转换为矩阵并使用 which(mtrx == max(mtrx), arr.ind = TRUE) 来获取矩阵位置，但是这个 (a) 丢失了地理空间信息和 (b) 导致数据在矩阵转换过程中旋转 90 度，这两者都需要额外的代码才能使其工作并显着减慢计算速度。是否有任何人都知道的等效栅格工作流程？

Answer 1

示例数据:

library(raster)
set.seed(0)
r <- raster(ncols=10, nrows=10)
r[] <- sample(50, 100, replace=T)

现在做:

p <- rasterToPoints(r, function(x) x == 11)

得到

       x   y layer
[1,]   18  81    11
[2,] -126  63    11
[3,]  -90  45    11
[4,]   54 -63    11

如果你想要具有最大值的单元格

vmax = maxValue(r)
p <- rasterToPoints(r, function(x) all.equal(x, vmax)

(不使用@data@max)