Сумма слоев растрового стека в R
Я работаю с ежедневным наблюдением за климатическими данными, организованными в файлах .nc. Я прочитал их, используя команду стека растрового пакета. Каждый файл (соответствующий году) является элементом RasterStack со следующими характеристиками:
class : RasterStack
dimensions : 360, 720, 259200, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.5, 0.5 (x, y)
extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
каждый слой - это растр значений дня.
Я хотел бы суммировать слои, чтобы рассчитать месячные значения. Я считаю, что решение должно использовать calc или stackApply {растр}, но я не смог найти способ суммирования от слоя x до слоя y или способ поднабора RasterStack перед суммированием.
Я подготовилфайл примера всего 12 слоев (чтобы уменьшить размер).
Я не знаю точно, как предложить код, извините, но он должен выглядеть примерно так:
library(raster)
setwd("myfolder")
data<-stack(mydata.nc)
datasum<- stackApply(data, ??? ,fun=sum)
Спасибо