如何用另一SpatRaster的指定行替换目标SpatRaster对应行的值？

错误原因 #

直接用values(x[1,,,,drop=FALSE]) <- values(y[1,,,,drop=FALSE])触发错误，核心问题在于x[1,,,,drop=FALSE]返回的是临时生成的SpatRaster子集副本，它和原对象x没有内存层面的绑定关系。你试图通过values()给这个临时子集赋值时，本质上只是修改了这个短命的临时对象，原对象x完全没被改动。terra包为了避免这种无效操作，就直接抛出了错误来提醒你。而先把子集存到中间变量再操作，其实是先创建了一个可修改的独立副本，之后再把修改后的副本赋值回原对象的对应位置，这时候是实实在在修改了x，所以能成功。

直接修改原对象的解决方案 #

这里提供几种无需中间变量就能直接完成替换的方法：

方法1：通过行索引定位单元格赋值 #

利用rows()获取目标行的单元格索引，直接操作原对象的values()矩阵：

library(terra)

# 示例数据
x <- rast(nrow=5, ncol=5, nlyr=2)
values(x) <- 1:ncell(x)
y <- rast(nrow=5, ncol=5, nlyr=2)
values(y) <- ncell(x):1

# 获取第1行的单元格索引
row_cells <- rows(x, 1)
# 直接给x的对应行赋值
values(x)[row_cells, ] <- values(y)[row_cells, ]

方法2：直接将y的行子集赋值给x的行子集 #

直接提取y的目标行子集，赋值给x的对应行子集，terra会自动处理内存关联：

# 直接替换x的第1行为y的第1行
x[1,,,,drop=FALSE] <- y[1,,,,drop=FALSE]

方法3：用setValues()精准赋值 #

如果需要更精细的单元格控制，可以用cells()获取目标行的单元格ID，再用setValues()赋值：

cell_ids <- cells(x, row=1)
x <- setValues(x, values(y)[cell_ids, ], cells=cell_ids)

以上三种方法都能直接修改原对象x，不需要额外的中间变量，也不会触发错误。

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