如何高效提取SpatRaster中值大于200的像元位置与属性信息？

针对大体积多层SpatRaster，避免使用as.array()或as.data.frame()带来的性能损耗，可以用以下高效方法提取符合条件的像元数据：

修正后的可复现代码 #

首先修正示例代码中的笔误（原代码中* *应为*）：

library(terra)
f <- system.file("ex/logo.tif", package="terra")
r <- rast(f, f) * seq(1,2, by = 0.2)

高效提取步骤 #

1. 创建逻辑掩码，标记各图层中数值大于200的像元

mask <- r > 200

2. 获取满足条件的像元索引（cell）和对应图层编号

matches <-```r
matches <- which(mask, cells = TRUE)

返回的matches是一个两列矩阵：第一列为像元的全局索引，第二列为对应的图层编号。

3. 提取对应像元的x、y坐标

coords <- xyFromCell(r, matches[, 1])

4. 提取对应像元的数值

values <- extract(r, matches[, 1], layer = matches[, 2])

5. 合并为最终数据框

result <- cbind(coords, layer = matches[, 2], value = values)

方法优势 #

该方法仅处理满足条件的像元，无需将整个大栅格数据加载到内存，大幅提升处理大体积SpatRaster的性能。

内容的提问来源于stack exchange，提问作者ahmathelte