如何在terra中覆盖栅格文件？解决writeRaster同名报错问题

当尝试直接覆盖正在被terra关联的栅格文件时，即使设置overwrite=TRUE也会报错，这是因为terra为避免数据损坏，禁止源文件与目标文件为同一文件的操作。以下是几种可行的解决方法：

方法1：将栅格数据加载到内存后再覆盖 #

读取栅格时强制加载到内存，断开与原文件的关联，之后即可正常覆盖写入：

library(terra)

x <- rast(ncol=10, nrow=10, nlyr=1)
x <- init(x, "cell")
x <- spatSample(x, ncell(x), "random", as.raster=TRUE)

f <- file.path(tempdir(), "test.tif")
writeRaster(x, f)

# 读取时直接加载到内存
y <- rast(f, mem=TRUE)
# 也可读取后调用readAll将数据全部载入内存
# y <- readAll(y)

# 添加新变量（示例操作）
y$new_var <- y[[1]] * 2

writeRaster(y, f, overwrite=TRUE)

方法2：通过临时文件中转 #

先将修改后的栅格写入临时文件，再替换原文件：

library(terra)

x <- rast(ncol=10, nrow=10, nlyr=1)
x <- init(x, "cell")
x <- spatSample(x, ncell(x), "random", as.raster=TRUE)

f <- file.path(tempdir(), "test.tif")
writeRaster(x, f)

y <- rast(f)
# 添加新变量
y$new_var <- y[[1]] * 2

# 创建临时文件路径
temp_f <- file.path(tempdir(), "temp_test.tif")
writeRaster(y, temp_f)

# 删除原文件并替换
file.remove(f)
file.rename(temp_f, f)

方法3：直接修改原文件（适用于添加图层场景） #

如果需求仅为添加新变量（图层），可直接用addLayer()函数在原文件上追加，无需覆盖整个文件：

library(terra)

x <- rast(ncol=10, nrow=10, nlyr=1)
x <- init(x, "cell")
x <- spatSample(x, ncell(x), "random", as.raster=TRUE)

f <- file.path(tempdir(), "test.tif")
writeRaster(x, f)

# 直接向原文件添加新图层
new_layer <- x[[1]] * 2
addLayer(f, new_layer, overwrite=TRUE)

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