Itération sur des fichiers NetCDF
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:
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
