如何在R中保存含多图层的SpatRaster对象并重新加载？

针对你的Climate SpatRaster对象，以下几种方法可以确保保存后重新加载时完全保留原格式、变量名及所有元数据：

方法1：使用多波段GeoTIFF格式（推荐） #

GeoTIFF是GIS领域通用格式，能完美保留SpatRaster的变量名、分辨率、坐标参考系等所有元数据，同时支持跨软件读取。

保存代码 #

library(terra)
# 写入多波段GeoTIFF，显式指定变量名确保一致性
writeRaster(Climate, "Climate.tif", overwrite = TRUE, names = names(Climate))

加载代码 #

# 直接读取即可恢复完整对象
Climate_reloaded <- rast("Climate.tif")

验证变量名：

names(Climate_reloaded)  # 应与原对象names(Climate)完全一致

方法2：使用R原生RDS格式 #

RDS是R专属的二进制序列化格式，会完整保存整个SpatRaster对象的所有属性，适合仅在R环境内使用的场景。

保存代码 #

saveRDS(Climate, "Climate.rds")

加载代码 #

Climate_reloaded <- readRDS("Climate.rds")

加载后对象的所有属性（包括变量名、时间信息、元数据）都会和原对象完全一致。

方法3：修复NetCDF格式的变量名问题 #

如果必须使用NetCDF格式，需要在保存时显式指定变量名参数，避免默认的自动命名（如band1、band2）。

保存代码 #

writeRaster(Climate, "Climate.nc", overwrite = TRUE, 
            varname = names(Climate),  # 指定NetCDF变量名
            longname = names(Climate)) # 可选，设置变量长描述

加载代码 #

Climate_reloaded <- rast("Climate.nc")

此时加载后的变量名将与原对象完全匹配。

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