如何使用grp函数计算不同年份同月栅格文件的平均值

Question

我有从2008年到2020年的每月栅格数据（.tiff）。文件命名如下：

index_200801
index_200802
..
..
..
index_202005

我想计算不同年份相同月份的平均值。

我在 R 中使用以下代码：

rast <- list.files((path="F:/Test/", pattern=".tiff$", full.names=TRUE) 
stk = stack(rast)
grp = rep(1:12, rep(??,12)

但是我不明白如何对文件进行分组

谁能帮我解决这个问题并导出情节。

在输出中，应该有 12 张图像，每个月一张，即 Jan（具有 2008 年 1 月、2009 年 1 月...2020 年 1 月的平均值）..同样适用于 2 月...12 月。

Answer 1

你可以做到

ff <- list.files((path="F:/Test/", pattern=".tiff$", full.names=TRUE) 
grp <- substr(basename(ff), 11, 12)

library(terra)
x <- rast(ff)
y <- tapp(x, grp, mean)