栅格堆栈像元中的底值到像元值的特定分位数r

Question

我正在努力处理栅格堆栈的地板像元值。这里是一些示例数据：

library(raster)

r <- raster(ncol=10, nrow=10)
s <- stack(lapply(1:365, function(i) setValues(r, runif(100, 0, 1))))
# adding NAs
s[[2]][sample(100, 25, TRUE)] <- NA
s[sample(100, 25, TRUE)] <- NA

我有一整年的时间序列栅格堆栈。在时间序列中，某些像元已完全设置为NA以及某些单个值。 我正在尝试将每个单元格中的值下限为该单元格的较低.025分位数。这是我的功能达到的程度：

f <- function(i) {
  v_q <- quantile(i,probs = .025 ,na.rm=TRUE)
  if (length(v_q) > 0) {
    if(i[v_q] < v_q) {
      i[v_q] <- v_q
    }
    i
  }
  i  
}

x <- calc(s, fun = f)

将其应用于calc时，会不断出现以下错误：

Error in .calcTest(x[1:5], fun, na.rm, forcefun, forceapply) : 
  cannot use this function

有什么想法吗？

Answer 1

不确定该错误，但是您的函数看起来有些奇怪。

首先，让我们添加一个随机种子以提高可重复性：

library(raster)

r <- raster(ncol=10, nrow=10)
set.seed(42)
s <- stack(lapply(1:365, function(i) setValues(r, runif(100, 0, 1))))
# adding NAs
s[[2]][sample(100, 25, TRUE)] <- NA
s[sample(100, 25, TRUE)] <- NA

calc在堆栈或砖上基本上适用于跨层的值向量。

所以这应该起作用：

f <- function(i) {
  q <- quantile(i,probs = .025 ,na.rm=TRUE)

  i[i<q] <- q
  i

}

x <- calc(s,fun=f)

并进行测试：

q <- calc(s,fun = function(x) quantile(x,0.025,na.rm=T))

par(mfrow=c(1,2))

plot(s[[1]] < q,main="s[[1]] < q")
plot(x[[1]] == q,main="x[[1]] == q")