在R中提取区域逐年最大值时掩码后仅得NA值的问题

你遇到的这个问题核心原因其实很明确——max()函数默认不会自动忽略NA值。因为掩码后的栅格里，海洋区域全是NA，陆地只有部分有效值，直接调用max(values(r))时，函数只要检测到存在NA就会返回NA，哪怕有大量非NA的陆地数据。

下面是具体的解决办法和优化建议：

1. 快速修复：给max()添加na.rm=TRUE参数 #

修改你现有的sapply代码，在max里开启忽略NA的设置：

Gmax <- sapply(unstack(land), function(r){max(values(r), na.rm=TRUE)})

这样函数就会跳过所有海洋区域的NA值，只计算陆地区域的最大值。

2. 更高效的替代方案：用calc()直接处理多层栅格 #

在raster包中，处理多层栅格的层统计有专门的calc()函数，比unstack+sapply的组合更简洁高效，也更符合栅格数据的处理逻辑：

# 直接对land栅格砖计算每一层的最大值（自动忽略NA）
Gmax <- calc(land, fun = function(x) max(x, na.rm = TRUE))
# 如果需要转成向量格式，再加一步
Gmax <- as.vector(Gmax)

3. 额外验证项（避免隐性问题） #

虽然你说绘图显示掩码已经生效，但还是可以快速检查一下原始数据和掩码图层的投影是否一致——如果投影不匹配，有时候会出现看似掩码生效但实际数据被错误遮蔽的情况：

# 检查两个图层的投影信息
projection(Model66)
projection(b)

# 如果投影不一致，先转换掩码图层的投影
if (projection(Model66) != projection(b)) {
  b <- spTransform(b, projection(Model66))
  land <- mask(Model66, b)
}

按上面的方法修改后，应该就能得到你想要的145年陆地全局最大值了。

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