我在R区用土处理环境层.我想修改落入多边形内的栅格的值.我知道classify函数可以更新栅格的所有值，但我不知道如何仅将其应用于多边形的区域.我还发现a nice and neat way可以改变多边形内的所有值，如下所示:
myraster[mypolygon] <- 7
但我只需要更新一个特定的值，而不是所有的，例如，将所有包含值21的单元格更改为值7.
我try 过，但没有成功:
myraster[mypoligon][myraster[mypoligon] == 21] <- 7
得到消息:
错误:[`[`]值太长
作为替代方案，我正在考虑获得与多边形相交的栅格，然后对值进行分类，并将这个小栅格与原始栅格合并，但我想知道是否存在更快的过程.
PS.我在下面贴了一个可复制的例子作为答案.
这里有一种使用cover到达那里的方法.
示例数据
library(terra)
p <- vect(system.file("ex/lux.shp", package="terra"))[3,]
r <- rast(system.file("ex/elev.tif", package="terra"))
让我们将p覆盖的区域r更新为400，但仅当该区域的值为400时.>
x <- mask(r, p)
y <- clamp(x, -Inf, 400)
z <- cover(y, r)
如果你想把一个特定的值，比如482，改成1000，你可以
yy <- classify(x, cbind(482, 1000))
zz <- cover(yy, r)
以及一系列的价值观
yyy <- classify(x, cbind(480, 500, 1000))
zzz <- cover(yyy, r)