如何在R包函数输出列表中集成terra对象并保证可移植性？

问题背景 #

我正在更新R包virtualspecies，将其依赖从raster切换到terra。此前该包函数的输出是附加了virtualspecies类的列表，其中包含1-2个RasterLayer或RasterStack对象，可通过saveRDS()导出并保留功能。现在需在输出列表中集成terra的SpatRaster对象，要求对象自包含且能跨R会话/电脑传输。我考虑了三种方案：

1. 使用terra的wrap()函数存储，需频繁封装/解封装，且proxy=FALSE可能导致大栅格或低内存环境下的问题；
2. 使用writeRaster()保存栅格并记录路径，但对象不再自包含，且跨系统移植复杂度高；
3. 优先用wrap()，内存不足时切换为writeRaster()并记录路径。

请问哪种方案最优？

方案分析与推荐 #

方案1：纯wrap()实现 #

• 优点：完全自包含，saveRDS()导出后跨会话/电脑直接可用，无需额外文件依赖；贴合原有包的使用习惯，用户无需处理外部文件。
• 缺点：频繁的封装/解封装会增加代码复杂度；当栅格数据量极大时，proxy=FALSE强制加载数据到内存，可能触发内存不足错误。

方案2：纯writeRaster()实现 #

• 直接排除。该方案彻底破坏对象自包含性，用户导出后必须同时拷贝栅格文件和RDS文件，跨系统时路径格式差异（Windows vs Unix）还会导致加载失败，大幅提升用户使用成本，不符合原有包的设计逻辑。

方案3：动态切换（优先wrap()，内存不足时降级为文件存储） #

• 这是最优方案，兼顾了自包含性和内存友好性：

  • 常规场景（栅格数据量小、内存充足）下，用wrap()封装SpatRaster，保持对象自包含，用户体验和原有包一致；
  • 检测到大栅格或内存不足时（可通过terra::ncell()估算数据量，或捕获wrap()的内存错误），自动调用writeRaster()将栅格保存到指定目录（或允许用户自定义路径），并在列表中存储路径信息+类型标记；
  • 加载时（比如为virtualspecies类实现自定义反序列化逻辑，或封装专用加载函数），自动识别标记：若为封装的SpatRaster则直接unwrap()，若为文件路径则用terra::rast()加载。

• 注意事项：

  • 需在包中实现自定义序列化/反序列化逻辑，比如为virtualspecies类添加serialize和unserialize方法，或封装save_virtualspecies()/load_virtualspecies()函数替代原生saveRDS()/readRDS()，确保跨环境兼容性；
  • 若使用自定义路径存储栅格，需提示用户导出时同步拷贝RDS文件和栅格文件，避免文件丢失；
  • 内存检测逻辑要稳健，可结合memory.limit()（Windows）或memuse::memuse()估算所需内存，提前判断是否需要降级到文件存储。

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