如何从日尺度Raster Brick计算年度总和与月度均值（含闰年处理）

你的核心问题是直接取日均值再乘固定天数会忽略闰年/平年的天数差异，正确的逻辑应该是先计算每年各月的实际日数据总和，再对这些年度月度总和求均值，具体步骤如下：

1. 生成年份+月份的分组索引
   首先提取图层对应的日期，生成包含年份和月份的索引，用来区分不同年份的同一月份：

   # 提取图层名称并转换为标准日期格式（去掉开头的X，替换点为横杠）
   layer_dates <- as.Date(sub("^X", "", names(brk)), format = "%Y.%m.%d")
   # 生成"年-月"格式的索引，如"1971-01"
   year_month_idx <- format(layer_dates, "%Y-%m")
   

2. 计算每年各月的日数据总和
   使用stackApply按年-月分组，对每个月的所有日数据求和，得到每年每个月的实际总量：

   # 计算每年各月的总和，每个图层对应某一年的某个月
   annual_monthly_sums <- stackApply(brk, indices = year_month_idx, fun = sum)
   

3. 计算同月份的年度总和均值
   对第一步得到的年度月度总和图层，按月份分组求均值，得到多年平均的月度总量：

   # 从年度月度总和的图层名称中提取月份索引
   sum_layer_dates <- as.Date(paste0(names(annual_monthly_sums), "-01"), "%Y-%m-%d")
   month_idx <- format(sum_layer_dates, "%m")
   # 对同一月份的年度总和求均值，最终得到12个图层（对应1-12月）
   mean_monthly_total <- stackApply(annual_monthly_sums, indices = month_idx, fun = mean)
   

逻辑说明 #

这种方法完全避开了闰年天数问题：

• 第一步的年-月索引确保每个分组都是某一年的完整月份，不管当月是28/29/30/31天，求和结果都是该月的实际总量；
• 第二步对同月份的年度总量取均值，直接得到多年平均的月度总数值，不需要再额外乘天数。

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