在R terra中基于多边形字段用值向量更新栅格值的问题

错误原因 #

你碰到的Error: [[] value is too long，核心问题是r[lux]返回的是栅格对象而非像元值的向量，直接给栅格对象赋值长度不匹配的向量就会触发报错。另外nrow(r[lux])得到的是栅格的行数，并非该区域实际包含的像元数量，这也是导致长度不匹配的关键。

正确实现方法 #

方法1：通过像元索引赋值 #

先提取多边形覆盖的像元索引，再用索引精准赋值：

library(terra)
p <- vect(system.file("ex/lux.shp", package="terra"))
r <- rast(p, res=0.01)
values(r) <- 1

lux <- p[p$NAME_1 == "Luxembourg",]
# 获取多边形覆盖的像元索引
cell_idx <- cells(r, lux)[,1]
# 给对应像元赋值向量
r[cell_idx] <- 1:length(cell_idx)

方法2：使用replace函数 #

利用terra::replace直接替换指定区域的像元值：

# 生成和目标区域像元数匹配的向量
val_vec <- 1:length(cells(r, lux)[,1])
r <- replace(r, lux, val_vec)

方法3：结合extract提取像元后赋值 #

先提取区域内的像元位置，再构造匹配向量赋值：

# 提取多边形覆盖的像元（返回矩阵，第一列为像元索引）
lux_cells <- extract(r, lux, cells=TRUE)[,1]
r[lux_cells] <- seq_along(lux_cells)

验证结果 #

可以通过以下代码查看更新后的栅格值：

# 查看指定区域的前几个像元值
head(values(r[lux]))

内容的提问来源于stack exchange，提问作者we need a Mat. Stat.