使用R语言terra包修改SpatRaster范围时原对象同步变更的问题

核心差异：terra与raster的对象拷贝机制不同 #

1. terra包的浅拷贝+写时复制策略

当执行r2 = r1创建SpatRaster对象副本时，terra采用的是浅拷贝：

• r2和r1会共享同一个底层元数据结构（包含范围、投影、行列数等信息），这是terra为优化内存占用和执行效率设计的默认行为。
• 因此修改r2的范围ext(r2) = ext(r2)/2时，实际修改的是两者共享的元数据，会同步影响到r1。
• 但修改数值时，terra会触发**写时复制（Copy-on-Write）**机制：只有当你对r2的数值进行修改操作时，才会真正创建独立的数据副本，所以不会影响r1的数值。

2. raster包的深拷贝机制

raster包中执行r2 = r1时，会直接创建一个完全独立的深拷贝：

• r2的元数据和数值数据都会和r1彻底分离，修改r2的范围或数值都只会作用于自身，不会牵连原对象r1。

解决方法：创建独立的SpatRaster对象 #

如果需要在terra中避免元数据共享，使用deepcopy()函数创建完全独立的对象：

library(terra)

r1 = rast(nrow = 2, ncol = 2)
values(r1) = c(0, 1, 2, 3)
r2 = deepcopy(r1)  # 执行深拷贝
ext(r2) = ext(r2)/2
ext(r2)  # 输出减半后的范围
ext(r1)  # 原范围保持不变

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