Bonjour,

Je travaille avec des gros jeux de données de type raster (Fichiers NetCDF4 sur les USA avec données climatiques quotidiennes par pixel de 1km²). Je souhaiterais faire la même chose sur chacun, mais pour cela, il faudrait que je puisse changer le nom du raster à chaque itération, avec le code ce sera sûrement plus clair :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
library("ncdf")
library("ncdf4")
library("raster")
 
for (i in (1990:2013)) {   ## Pour chaque année j'ai un fichier en 3 dimensions (latitude, longitude et jours) avec 365 bandes correspondant aux 365 jours
 
file<-paste("tmin_",i,".nc4", sep="")
file2<-brick(file)
early<-stack(file2[[1:120]])   ## Sélections de la fin de l'hiver (début d'année i) (jours 1 à 120)
late<-stack(file2[[305:365]])  ## Sélection du début de l'hiver suivant (fin d'année i) (jours 305 à 365)
 
}
Je voudrais pouvoir créer des stacks pour chaque hiver de la façon suivante : stack(late1990, early1991), mais pour cela, il faudrait que je puisse nommer différemment mes stacks de sortie de boucle à chaque itération...

Est-ce que quelqu'un aurait une solution ?
Merci d'avance !
Eva