在R的terra包中修改多边形内指定值的栅格数据

可复现示例 #

先创建测试数据，方便验证：

library(terra)
# 生成随机值栅格
r <- rast(nrow=10, ncol=10, vals=sample(1:25, 100, replace=TRUE))
# 创建覆盖中心区域的多边形
p <- vect("POLYGON ((-50 50, 50 50, 50 -50, -50 -50, -50 50))", crs=crs(r))

正确解决方案 #

方法1：利用单元格索引精准赋值

先提取多边形覆盖的单元格索引，再针对特定值修改：

# 获取多边形对应的栅格单元格索引
cell_idx <- cells(r, p)[, 1]
# 仅修改索引内值为21的单元格
r[cell_idx][r[cell_idx] == 21] <- 7

方法2：用ifel函数实现条件赋值（更简洁）

terra::ifel可以直接根据条件批量修改栅格值，结合mask限定区域：

# 仅在多边形区域内，将值为21的单元格替换为7，其余保持不变
r <- ifel(mask(r, p) == 21, 7, r)

你之前代码报错的原因 #

你尝试的myraster[mypolygon][myraster[mypolygon] == 21] <- 7逻辑有问题：

• 第一次myraster[mypolygon]返回的是多边形内所有单元格值的向量
• 第二次子集[myraster[mypolygon] ==21]只筛选出其中等于21的元素，此时你要赋值的元素数量远小于myraster[mypolygon]的总长度，R无法匹配赋值长度，因此抛出Error: [[] value is too long的错误。

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