无法在 R 中的栅格堆栈上执行 median()

Question

我试图在R中找到一组光栅的每个单元格中位数，但是当我使用median函数时，我得到以下错误:

Error in if (any(is.na(x))) return(x[FALSE][NA]) : 
  argument is not interpretable as logical

然而，我可以在我的堆栈上执行max, mean和calc(r, median)，所以我怀疑这不是数据格式(tif)的问题。我使用的代码是:

数据加载:

f = list.files("./")
r = stack(f)

做总结:

median(r) # does not work
mean(r) # works
max(r) # works
calc(r, median) # works

作为calc正在工作，我不需要一个答案继续，但我想了解为什么median不工作。这是标准行为吗?

Answer 1

根据参考手册p. 201:

以下总结方法可用于栅格*对象:mean, max, min, range, prod, sum, any, all

median不是Raster*对象可用的摘要方法。这就是为什么你需要使用calc函数