在R中计算时序栅格趋势报错问题求助

问题场景 #

使用terra包的app()函数计算时间序列栅格的线性趋势时，触发如下错误：

Error: Not compatible with requested type: [type=list; target=double].

原始代码与数据信息：

library(terra)
trend <- function(x){
  if(is.na(x[1])){rep(NA,length(x))}
  else{
 summary(lm(as.numeric(x)~c(2015:2100)))$coefficients[2]
   }
}
# 数据信息
data
#> class       : SpatRaster
#> dimensions  : 40, 81, 86  (nrow, ncol, nlyr)
#> resolution  : 0.5, 0.5  (x, y)
#> extent      : 46.5, 87, 35.5, 55.5  (xmin, xmax, ymin, ymax)
#> coord. ref. : lon/lat WGS 84
#> source(s)   : memory
#> names       :      y_2015,       y_2016,       y_2017,       y_2018,       y_2019,       y_2020, ...
#> min values  : 0.00000e+00, 0.000000e+00, 0.000000e+00, 0.000000e+00, 0.000000e+00, 0.000000e+00, ...
#> max values  : 6.49618e-10, 6.122992e-10, 5.359198e-10, 6.447496e-10, 6.671828e-10, 6.464858e-10, ...
#> time (years): 2015 to 2100

# 执行计算触发错误
data |> app(trend)

错误原因 #

1. 返回值类型/长度不匹配：app()函数要求每个栅格单元的处理函数返回单个数值（用于生成单波段栅格），但原始函数中：
   • 第一个分支返回长度为86的NA向量（与栅格列数一致）
   • 第二个分支返回带名称属性的lm系数，虽为单个值，但与向量分支的返回类型不统一
2. NA判断逻辑不完善：仅检查第一个元素是否为NA，无法覆盖栅格单元部分值为NA的情况，可能导致后续lm拟合失败。

修正后的代码 #

library(terra)

trend <- function(x) {
  # 检查是否全为NA，或非NA值不足2个（无法拟合线性模型）
  if (all(is.na(x)) || sum(!is.na(x)) < 2) {
    return(NA_real_)  # 返回双精度类型的NA，匹配数值要求
  } else {
    # 用数据框传递变量更规范；直接提取lm系数并转为纯数值
    mod <- lm(as.numeric(x) ~ year, data = data.frame(x = x, year = 2015:2100))
    return(as.numeric(coef(mod)[2]))
  }
}

# 计算趋势栅格
trend_raster <- app(data, trend)

关键修正点 #

• 统一返回值为单个双精度数值（NA_real_是标准数值型NA），确保与app()的类型要求匹配
• 完善NA判断逻辑，避免因部分NA或全NA导致的拟合错误
• 直接使用lm()的coef()提取系数，比summary()更高效；用as.numeric()去除系数的名称属性，确保返回纯数值
• 用数据框传递自变量，避免潜在的长度不匹配问题，代码更易读

内容的提问来源于stack exchange，提问作者ChuDingjin