Terra自定义FGT函数计算栅格全局统计结果异常原因咨询

针对你用terra包实现Alpha=0的FGT指数时返回结果为1（预期0.5）的问题，核心是函数计算逻辑与理论定义不符，以下是常见错误点：

1. 贫困像素的判定逻辑颠倒 #

Alpha=0的FGT指数本质是贫困像素数与有效总像素数的比值，可达性贫困通常定义为「旅行时间超过阈值的像素为贫困」。如果你的函数误将「小于等于阈值」的像素判定为贫困，会导致统计的贫困像素数等于总像素数，最终返回1。

• 错误示例：sum(r <= threshold)
• 正确逻辑（根据贫困定义调整符号）：sum(r > threshold, na.rm=TRUE)

2. 未正确处理NA值 #

若测试栅格存在缺失值，函数中总像素数与贫困像素数的统计口径不一致会导致结果偏差：

• 如果总像素数用ncell(r)（包含NA），但贫困像素数统计时未排除NA（na.rm=FALSE），会把NA误判为贫困像素；
• 反之，若总像素数用length(na.omit(values(r)))，但贫困像素数统计时包含NA，也会导致比值异常。

3. 栅格值提取或计算的语法错误 #

使用terra包时，若直接对栅格对象进行统计而非提取其数值，可能触发非预期的计算逻辑：

• 错误示例：直接对栅格对象执行sum(r > threshold)（部分场景下会误统计栅格元数据而非像素值）；
• 正确方式：先提取像素值r_vals <- values(r)，再基于数值向量计算。

4. 阈值计算错误 #

若阈值（全局中位数）的计算未处理NA或提取值错误，会导致判定基准偏差：

• 错误示例：threshold <- median(r)（未指定na.rm=TRUE，若存在NA会得到无效阈值）；
• 正确方式：threshold <- median(values(r), na.rm=TRUE)

快速验证步骤 #

1. 手动提取栅格数值：r_vals <- na.omit(values(your_raster))
2. 计算正确阈值：threshold <- median(r_vals)
3. 统计贫困像素数：poor_count <- sum(r_vals > threshold)
4. 计算预期比值：poor_count / length(r_vals)
   将上述结果与自定义函数的每一步输出对比，即可定位错误环节。

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