r - 如何在每一层应用光栅堆栈对象？

Question

我想在堆栈的每一层上应用函数并返回一个堆栈。

library(raster)
r1 <- raster(vals=1:20,nrow=4,ncols=5)
r2 <- raster(vals=2:21,nrow=4,ncols=5)
r3 <- raster(vals=3:22,nrow=4,ncols=5)
stk <- stack(r1,r2,r3)
func <- function(x){
  calc(x,function(y){
    ifelse(y>10,0,y)
  })
}

我使用lapply在每一层上应用func

lapply(stk,func)
Error in (function (classes, fdef, mtable)  : 
unable to find an inherited method for function ‘calc’ for signature ‘"integer", "function"’ 

看起来lapply切片stk为stk[1]stk[2]...，但是实际层是 stk[[1]],stk[[2]]....

并使用stackApply:

stackApply(stk,indices=1:3,fun=func)
Error in FUN(newX[, i], ...) : unused argument (na.rm = na.rm)

Answer 1

看起来您可以单独使用calc:

stack(calc(stk, function(y) ifelse(y > 10, 0, y)))
# class       : RasterStack 
# dimensions  : 4, 5, 20, 3  (nrow, ncol, ncell, nlayers)
# resolution  : 72, 45  (x, y)
# extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
# names       : layer.1, layer.2, layer.3 
# min values  :       0,       0,       0 
# max values  :      10,      10,      10