如何用R计算月度NDVI栅格

Question

NDVI RasterStack

我有一个15天的1981-2015 NDVI RasterStack。

我需要使用15天的数据来计算每月NDVI。

我想知道如何计算1981-2015年每月的同名MM栅格的均值

感谢您的帮助！非常感谢。

names XYYYY.MM.DD

Answer 1

我最近一直在为您解决相同的问题，这也应该对您有用。

您要创建一个单独的变量，其中包含图层名称中的日期。

#this removes the "X" character from the name leaving only the dates

layer_name <- sub('.', '', names(NDVI_stack))

#install.packages("lubridate")
library(lubridate)

layer_name <- ymd(layer_name)

#Create an indices to prepare it for stackApply, which takes the means for all the days of the month within each year.
indices <- format(as.Date(layer_name, format = "%Y.%m.%d"), format = "%Y.m") 

NDVI_mean <- stackApply(NDVI_stack, indices, mean)