Raster tuğla nasıl döngü yapılır?

Question

Fraksiyonel kar örtüsünü içeren bir SCA (nrow = 108, ncol = 132, nlayers = 365) raster tuğlasına sahibim. Bundan 8 katmandan 46'şer yığın yapmak ve bu 46 yığından maksimum fraksiyonel kar örtüsünü hesaplamak istiyorum. Bunu nasıl yapabilirim?

Answer 1

Senden isteyebilirsiniz düşünüyorum bu şekilde o: Bir döngü istedi (ama bu R, döngüler önlemek için deneyin) Eğer yapabilirsem

library(raster) 
# example data 
sca <- brick(nrow=108,ncol=132,nl=365) 
values(sca) <- runif(ncell(sca)*nlayers(sca)) 

# indices grouping sets of 8 
i <- rep(1:ceiling(365/8), each=8) 
# the last period is not a complete set of 8 days 
i <- i[1:nlayers(sca)] 

x <- stackApply(sca, i, max) 

for (i in 1:nlayers(sca)) { 
    x <- sca[[i]] 
    # etc. 
}