Ich arbeite mit der täglichen Beobachtung von Klimadaten, die in .nc-Dateien organisiert sind. Ich habe sie mit dem Stack-Befehl des Raster-Pakets gelesen. Jede Datei (die einem Jahr entspricht) ist ein RasterStack-Element mit den folgenden Merkmalen:

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) 

Jede Ebene ist das Raster der Werte eines Tages.
Ich möchte die Schichten summieren, um die monatlichen Werte zu berechnen. Ich glaube, die Lösung sollte calc oder stackApply {raster} verwenden, aber ich konnte keinen Weg finden, von Layer x zu Layer y zu summieren oder den RasterStack vor der Summierung zu subset- zen.

Ich habe ein @ vorbereitexample file mit nur 12 Schichten (um die Größe zu reduzieren).

Ich weiß nicht genau, wie ich einen Code vorschlagen soll, sorry, aber es sollte ungefähr so aussehen:

library(raster)
setwd("myfolder")
data<-stack(mydata.nc)

datasum<- stackApply(data, ??? ,fun=sum)

Vielen Dan