使用terra::predict()直接预测SpatRaster失败，求技术支持

在处理大规模栅格数据时，无法将SpatRaster转换为data.frame后进行预测，尝试使用terra包的predict()直接调用随机森林模型（ranger/randomForest）时失败，且原本正常运行的旧脚本也出现相同问题，怀疑是terra近期版本更新导致兼容性问题。

可复现代码 #

library(terra)
library(ranger)

r <- rast(ncols=95, nrows=90, nlyrs=5, crs="+proj=longlat",
          xmin=5.74167, xmax=6.5333, ymin=49.44167, ymax=50.19167,
          names=c("band1","band2","band3","band4","band5"))

# 填充随机值
set.seed(123)
r$band1 <- runif(n = ncell(r))
r$band2 <- init(r, "cell")
r$band3 <- sample(1:10, ncell(r), replace = T)
r$band4 <- sample(1:10, ncell(r), replace = T)
r$band5 <- sample(1:10, ncell(r), replace = T)

# 提取训练样本
extr_values <- spatSample(r, 10)
extr_values$response <- as.factor(sample(c(1,2), nrow(extr_values), replace = T))

# 构建ranger模型
rf_model <- ranger(response ~., data = extr_values, num.trees = 50)

# 尝试预测（失败）
prediction <- predict(r, rf_model,  na.rm = TRUE)$predictions

解决方案 #

1. 使用自定义预测函数显式处理模型调用 #

terra新版本对模型的预测接口兼容性要求更严格，需要显式定义如何将栅格数据传入模型的预测函数。针对不同随机森林包，自定义适配函数：

适配ranger模型

# 自定义预测函数，处理ranger的输入输出逻辑
rf_predict_fun <- function(model, data) {
  predict(model, data = data)$predictions
}

# 执行预测
prediction <- predict(r, rf_model, fun = rf_predict_fun, na.rm = TRUE)

适配randomForest模型

library(randomForest)
rf_model_rf <- randomForest(response ~., data = extr_values, ntree = 50)

# 适配randomForest的预测函数
rf_predict_fun_rf <- function(model, data) {
  predict(model, newdata = data)
}

prediction <- predict(r, rf_model_rf, fun = rf_predict_fun_rf, na.rm = TRUE)

2. 检查变量名匹配性 #

确认栅格的图层名称（names(r)）与训练数据中的特征列名完全一致，包括大小写、拼写，避免因名称不匹配导致模型无法定位对应特征。

3. 降级terra到兼容版本 #

如果上述方法无效，可尝试降级到之前能正常运行的terra版本：

# 使用pak包安装指定版本
install.packages("pak")
pak::pak("rspatial/terra@1.7-39") # 替换为你之前正常使用的版本号

4. 显式指定分类模型类型 #

针对分类任务，在模型训练和预测时明确指定类型，避免terra自动推断出错：

# 训练ranger模型时明确指定分类类型
rf_model <- ranger(response ~., data = extr_values, num.trees = 50, classification = TRUE)

# 自定义函数中明确返回类别
rf_predict_fun <- function(model, data) {
  predict(model, data = data, type = "response")$predictions
}

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