ラスターレンガをループする方法は？

Question

私は小雪を含むSCA（nrow = 108、ncol = 132、nlayers = 365）のラスタブリックを持っています。私はこれから8枚のレイヤーをそれぞれ46枚ずつ積み重ねて、これらのスタックから最大の積雪量を計算したいと思います。私はこれをどうやって行うことができますか？

Answer 1

は、私はあなたがやりたいことだと思い、このようこと：あなたがループを望んでいた（これはR、ループを回避しようで）あなたができる場合は

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. 
}