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 :
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...
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) }
Est-ce que quelqu'un aurait une solution ?
Merci d'avance !
Eva
Partager