从定义的间隔中选择栅格堆栈图层

时间:2016-06-28 14:12:59

标签: r raster

我有一个包含744个图层的光栅堆栈s。其中,744,我需要制作具有第18至24层(即18层中的6层)的新堆栈S1,然后在距离18(即42至48)24的间隙处制作下6层。所以我需要在24的间隙继续挑选6层,直到我到达堆栈的末尾。所以我所做的是:

x<-data.frame(seq(18, 744, by = 24))
x$v2<-seq(24, 744, by = 24)
S1<-stack(s[[x[1,1]:x[1,2]]])
for(i in 2:nrow(x)){
  S1<-stack(S1,s[[x[i,1]:x[i,2]]])
}

但是,我想知道是否有更清洁,更好的方法来做到这一点。

1 个答案:

答案 0 :(得分:2)

你可以做到

i <- rep(18:24, 31) + rep(0:30, each=7) * 24
S2 <- S1[[i]]