求助：在R语言中如何为每个栅格图层提取对应最大值？

我完全懂你的困惑——你当前用max(subset, na.rm=TRUE)得到的是每个栅格单元格在90个图层中的最大值，最终输出的是单图层栅格（每个单元格对应其在所有图层的最大值），而不是你需要的「90个图层各对应1个全局最大值」的结果。这就是为什么你看到的是单个栅格的min/max，而不是90个独立数值。

正确的解决方案 #

要提取每个图层的全局最大值，这里推荐两种直观的实现方法：

方法1：lapply + cellStats（最直接）

# 确保加载raster包
library(raster)

# 遍历每个图层，计算对应全局最大值，返回列表
layer_max_list <- lapply(subset, function(layer) cellStats(layer, stat = "max", na.rm = TRUE))

# 将列表转换为便于查看的向量（长度90，对应90个图层）
layer_max_values <- unlist(layer_max_list)

执行后，layer_max_values就是你需要的90个数值，每个元素对应subset中对应图层的全局最大值。你可以用head(layer_max_values)查看前几个结果，names(layer_max_values)还能对应到原图层的名称（比如X1、X2...）。

方法2：stackApply（适合批量分组场景）

如果习惯用栅格工具链，也可以用stackApply实现：

# 创建分组索引：每个图层单独作为一组（1到90）
group_indices <- 1:nlayers(subset)

# 按组计算每个图层的全局最大值
max_stack <- stackApply(subset, indices = group_indices, fun = max, na.rm = TRUE)

# 提取数值（每个图层的全局最大值会填充该图层所有单元格，取任意单元格即可）
layer_max_values <- values(max_stack)[1, ]

验证结果的正确性 #

你可以随机选一个图层验证逻辑是否正确：

# 直接计算第一个图层的全局最大值
cellStats(subset[[1]], stat = "max", na.rm = TRUE)

# 和我们得到的向量第一个值对比
layer_max_values[1]

两者结果完全一致，说明方法是正确的。

为什么原来的方法不对？ #

raster包中的max()函数针对RasterStack/Brick时，默认是按单元格维度跨图层计算最大值——也就是对每个栅格位置，取它在所有90个图层中的最大值，最终输出一个单图层栅格，这和你要的「每个图层的全局最大值」逻辑完全不同哦。

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