R栅格时间序列分解错误答案

【问题标题】：R raster timeseries decomposition errorR栅格时间序列分解错误
【发布时间】：2018-01-10 13:04:31
【问题描述】：

我正在尝试分解每月时间序列堆栈（维度：336、221、74256、420）以获得随机分量。我正在使用以下内容：

tsfun <- function(x) {  
  my.ts = na.remove(ts(x, start=c(1982,1), frequency=12)) 
  my.decomp = decompose(my.ts, type="multi") 
  my.exp = my.decomp$random 
}
random <- calc(ras_in, fun=tsfun)

但我收到以下错误

.calcTest(x[1:5], fun, na.rm, forcefun, forceapply) 中的错误：无法使用此功能

有什么建议吗？谢谢一百万

【问题讨论】：

使用包含至少 2-3 个 NA 值的 sn-p 数据会更容易。所以我们可以在发布之前测试我们的理论。它不需要全部，只需dput(x) 整个内容并提取一个块，但在发布之前确保该块也失败。

标签： r time-series r-raster decomposition

【解决方案1】：

我认为您的功能无效。当我测试它时，我得到：

tsfun(1:24)
#Error in na.remove(ts(x, start = c(1982, 1), frequency = 12)) : 
#  could not find function "na.remove"

这是一个适合我的版本：

tsfun2 <- function(x) {  
  my.ts = na.omit(ts(x, start=c(1982,1), frequency=12)) 
  decompose(my.ts, type="multi")$random 
}

tsfun2(1:24)
#     Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
#1982  NA  NA  NA  NA  NA  NA   1   1   1   1   1   1
#1983   1   1   1   1   1   1  NA  NA  NA  NA  NA  NA

【讨论】：