如何为Terra包的Raster对象添加计算列并保存？

不用非得把栅格转成data.frame才能加新列，terra包的rast对象支持直接操作，不管是单个新增还是批量加20列都很方便，下面分两种场景说明：

一、直接对rast对象操作 （推荐） #

1. 单个添加新图层 #

直接通过赋值就能给原栅格加新图层，和你在data.frame里加列逻辑一致：

library(terra)

# 加载示例数据
f <- system.file("extdata/asia.tif", package = "tidyterra")
xx <- rast(f)

# 直接计算并添加
xx$col1 <- cells(xx) * 5
xx$col2 <- cells(xx) * 20
xx$col3 <- cells(xx) * 50

# 查看结果，能看到新增的图层
xx

这里cells(xx)直接获取像元编号，和你转data.frame后的cell列完全对应。

2. 批量添加20个图层 #

如果要加20列，用循环或者lapply批量生成后合并就行：

# 定义20个乘数，比如从5到100每次加5
multipliers <- seq(5, 100, by = 5)

# 批量生成新图层
new_layers <- lapply(multipliers, function(m) cells(xx) * m)

# 给新图层命名，比如col1到col20
names(new_layers) <- paste0("col", 1:length(multipliers))

# 合并到原栅格
xx <- c(xx, rast(new_layers))

# 检查图层数量，原图层+20个新图层
nlyr(xx)

二、从已有的data.frame导回栅格 #

如果你已经在data.frame里完成了计算，也能把新增列转成栅格合并回去：

# 你的原有代码
rast_df <- as.data.frame(xx, na.rm = TRUE, cells = TRUE)
rast_df$col1 <- with(rast_df, cell*5)
rast_df$col2 <- with(rast_df, cell*20)
rast_df$col3 <- with(rast_df, cell*50)

# 提取所有新增的col开头的列
new_cols <- rast_df[, grepl("^col", names(rast_df))]

# 把data.frame转成和原栅格匹配的rast对象
new_rast <- rast(xx, nlyr = ncol(new_cols))
values(new_rast) <- new_cols
names(new_rast) <- names(new_cols)

# 合并到原栅格
xx <- c(xx, new_rast)

三、保存最终栅格 #

完成所有添加后，用writeRaster保存成你需要的格式：

# 保存为GeoTIFF，允许覆盖已存在的文件
writeRaster(xx, "asia_with_new_columns.tif", overwrite = TRUE)

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