如何在corLocal函数中忽略栅格NA值并解决结果为空问题

问题原因 #

你遇到的问题主要来自两个点：

• corLocal默认的na.rm参数无法直接处理窗口内有效样本数不足的情况（比如窗口内r1的非NA值少于2个，或r1/r2的窗口值全部相同），这种情况下Spearman秩相关系数无法计算，直接返回NA。
• corLocal没有自动同步r1的NA区域，导致结果全NA。

解决方案 #

步骤1：自定义局部Spearman计算函数 #

用focalPair（corLocal的底层依赖函数）自定义逻辑，先过滤NA并检查样本有效性，再计算相关系数：

library(raster)

# 自定义局部Spearman相关计算函数
local_spearman <- function(x, y) {
  # 筛选r1非NA的有效配对
  valid_idx <- !is.na(x)
  x_valid <- x[valid_idx]
  y_valid <- y[valid_idx]
  
  # 仅当有效样本数≥2，且两组值都存在变异时计算相关
  if (length(x_valid) >= 2 && length(unique(x_valid)) > 1 && length(unique(y_valid)) > 1) {
    return(cor(x_valid, y_valid, method = "spearman", use = "complete.obs"))
  } else {
    return(NA)
  }
}

步骤2：计算局部相关并同步NA区域 #

用focalPair执行窗口计算，再用mask将结果的NA区域与r1对齐：

# 3x3窗口计算（可根据需求调整为5x5等更大窗口）
result <- focalPair(r1, r2, fun = local_spearman, window = matrix(1, 3, 3))

# 同步r1的NA区域到结果
result <- mask(result, r1)

可选调整：增大窗口尺寸 #

如果r1的NA分布较分散，3x3窗口有效样本不足，可以尝试更大的窗口（比如5x5）：

result <- focalPair(r1, r2, fun = local_spearman, window = matrix(1, 5, 5))
result <- mask(result, r1)

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