使用terra包保存加载SpatRasters时遇外部指针无效错误求助

SpatRaster对象依赖外部指针关联底层栅格数据，这类指针仅在创建它的当前R会话中有效。直接用save()保存、load()加载时，指针无法跨会话复用，就会触发"external pointer is not valid"错误。

提供两种实用方法，按需选择：

方法1：用terra::wrap()/unwrap()序列化栅格对象 #

该方法可以在R环境内完整保存栅格状态，无需导出到外部文件：

library(terra)
x <- terra::rast(xmin=-110, xmax=-80, ymin=40, ymax=70, ncols=30, nrows=30)
values(x) <- 1:ncell(x)
r <- c(x, x, x, x)
r <- list(r, r, r, r)

# 序列化栅格列表（嵌套列表需逐层处理）
r_wrapped <- lapply(r, function(sub_list) lapply(sub_list, wrap))
save(r_wrapped, file = "test_wrapped.Rdata")

# 加载并反序列化
rm(list=ls(all=TRUE))
load("test_wrapped.Rdata")
r_unwrapped <- lapply(r_wrapped, function(sub_list) lapply(sub_list, unwrap))

# 验证正常输出
r_unwrapped

方法2：导出为标准栅格文件后重新读取 #

如果需要长期存储或跨工具共享，推荐将栅格导出为GeoTIFF等标准格式：

library(terra)
x <- terra::rast(xmin=-110, xmax=-80, ymin=40, ymax=70, ncols=30, nrows=30)
values(x) <- 1:ncell(x)
r <- c(x, x, x, x)
r <- list(r, r, r, r)

# 创建存储目录
dir.create("raster_files", showWarnings = FALSE)

# 批量导出所有栅格
for (i in seq_along(r)) {
  for (j in seq_along(r[[i]])) {
    writeRaster(
      r[[i]][[j]],
      filename = paste0("raster_files/raster_", i, "_", j, ".tif"),
      overwrite = TRUE
    )
  }
}

# 清空环境后重新读取构建列表
rm(list=ls(all=TRUE))
library(terra)

r_loaded <- list()
for (i in 1:4) {
  sub_list <- list()
  for (j in 1:4) {
    sub_list[[j]] <- rast(paste0("raster_files/raster_", i, "_", j, ".tif"))
  }
  r_loaded[[i]] <- sub_list
}

# 验证正常
r_loaded

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