Je voudrais importer tous les rasters de mosaïques multiples (chaque mosaïque est composée de 7 couches) du dossier indiqué dans mon R. Ensuite, accédez-y en tant que rasters multicouches uniques.
J'ai donc fait:
# read all mosaics named "mos....img" in R
raster_data<-list.files(path=getwd(), pattern="mos.*.img$")
# read files as rasters
s <- stack(raster_data)
# check my imported rasters p.ex. raster n°8 from "s" raster stack
s[[8]]
et mon raster s [[8]] ne contient qu'une seule couche, donc toute la mosaïque n'a pas été importée!
nlayers(s[[8]])
[[1]]
Si je lis chaque mosaïque séparément, cela fonctionne:
# read 1 mosaic (composed by 7 bands)
mosaic1<-brick("mosaic1.img")
# extract one band
band4<-subset(mosaic1, 4)
Pourquoi l'outil "empiler" n’importe pas des mosaïques entières mais seulement une bande de la mosaïque et comment est-il possible de l’arranger?