如何在R中用grp函数分组计算不同年份同月栅格均值并导出图像

步骤1：加载必要包 #

首先确保安装并加载raster包（若处理大尺寸数据，推荐用新版terra包，后文附替代方案）：

library(raster)

步骤2：读取并堆叠栅格文件 #

先修正你代码里list.files的语法错误，再批量读取栅格并堆叠：

# 列出指定路径下所有.tiff文件（带完整路径）
rast_files <- list.files(path = "F:/Test/", pattern = "\\.tiff$", full.names = TRUE)
# 堆叠所有栅格
stk <- stack(rast_files)

步骤3：从文件名提取月份分组 #

你原代码用rep的方式不适用，因为不同月份的年份数量不一致（1-5月有13年数据，6-12月只有12年）。直接从文件名提取月份作为分组依据更可靠：

# 获取文件名（去掉路径）
file_names <- basename(rast_files)
# 提取月份：文件名中第7-8位是月份（如index_200801.tiff中的"01"）
months <- as.integer(substr(file_names, 7, 8))

步骤4：计算各月份均值 #

用stackApply按月份分组计算均值：

# 按月份分组计算栅格均值，得到包含12个图层的栈（对应1-12月）
monthly_means <- stackApply(stk, indices = months, fun = mean)

步骤5：导出每个月的均值栅格 #

循环导出每个图层为独立的TIFF文件，文件名统一格式：

# 创建输出目录（如果不存在）
output_dir <- "F:/Test/Monthly_Means/"
if (!dir.exists(output_dir)) dir.create(output_dir)

# 逐个导出12个月的均值栅格
for (i in 1:nlayers(monthly_means)) {
  # 生成带两位月份的文件名（如mean_month_01.tiff）
  output_file <- paste0(output_dir, "mean_month_", sprintf("%02d", i), ".tiff")
  # 导出栅格，允许覆盖已存在文件
  writeRaster(monthly_means[[i]], filename = output_file, format = "GTiff", overwrite = TRUE)
}

新版terra包替代方案（更高效） #

如果你处理的栅格数据量较大，推荐使用terra包替代raster：

library(terra)
# 读取并堆叠栅格
rast_files <- list.files(path = "F:/Test/", pattern = "\\.tiff$", full.names = TRUE)
stk <- rast(rast_files)
# 提取月份分组
months <- as.integer(substr(basename(rast_files), 7, 8))
# 计算月度均值
monthly_means <- tapp(stk, months, mean)
# 创建输出目录并导出
output_dir <- "F:/Test/Monthly_Means/"
dir.create(output_dir, showWarnings = FALSE)
writeRaster(monthly_means, filename = paste0(output_dir, "mean_month_", sprintf("%02d", 1:12), ".tiff"), overwrite = TRUE)

为什么你原来的分组方法行不通？ #

rep(1:12, rep(??,12))需要手动指定每个月份对应的年份数量，但你的数据到2020年5月就结束了，1-5月有13年数据（2008-2020），6-12月只有12年（2008-2019），手动填数容易出错且不灵活。直接从文件名提取月份的方法能自动适配数据范围，更鲁棒。

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