每6个栅格图层的平均值和此文件的名称

Question

我在一个文件夹中组合了28个栅格文件。我想做3个第一个栅格的平均值，然后是另外3个栅格的平均值等等。然后，我应该获得9个栅格。我想根据月份来命名。我希望我的第一个栅格图层命名为Tmin-2010-01-15，我的第二个栅格图层命名为“Tmin-2010-02-15”，...，以及我的最后一个和9个栅格名为'Tmin-2010-09-15 ”。我编写了以下代码。我没有任何错误消息，但我看不到任何结果，所以我不知道代码的逻辑是否良好以及如何处理。

感谢您的帮助！

if

我堆积了我的栅格

$status

我定义参数n =栅格图层数，m =我想要收集的栅格数量

setwd("/Users/lisemontefiore/Documents/Data/Climate Data/test2") 

我定义了函数创建的文件名，并编写了栅格

bands <-list.files(path="xxx",pattern="*.tif",full.names=TRUE)
data_stack <- stack(bands) 

Answer 1

您可以这样做，而不是计算：

 b <- stackApply(b, group3, mean)

但没有输出，因为你没有调用writeRaster 你需要这样的东西：

months <- seq(from=as.Date("2010-01-15"), to= as.Date("2010-09-15"), by='month')
ff <- paste0('Tmin', months, '.tif')

for (i in seq_along(months)) {
    writeRaster(b[[i]]], ff[[i]])
}