使用空间多边形更改栅格值

Question

要更改 SpatialPoints 下的栅格值，您只需使用 [ 即可。

r <- raster(system.file("external/test.grd", package="raster"))
rTp <- rasterToPoints(r, spatial = T)

set.seed(666)
rTpS <- rTp[sample(1:length(rTp), 500),]
plot(r)
plot(rTpS, add = TRUE, pch = ".")

现在，空间点下方的栅格值更改如下：

r1 <- r
r1[rTpS] <- 99999
plot(r1)

具有 SpatialPoints 的栅格（左）和具有更改值的栅格（右）：

因此，我假设多边形的工作方式相同。因此从这个栅格数据集生成一些多边形：

rst <- as.integer(stretch(r, 0, 10))
rTpol <- rasterToPolygons(rst, dissolve = T)
rTpol <- rTpol[rTpol$layer > 3,]

plot(r)
plot(rTpol, add = T)

rst[rTpol[,]] <- 100
plot(rst)

# also tried 
# rst[rTpol] <- 100

但是，当我尝试将栅格值更改为低于某些 SpatialPolygons 时，多边形的 ID（或其他内容）不知何故。然而，它应该是 100。

具有 SpatialPolygons 的栅格（左）和具有奇怪变化的值的栅格（右）：

因此，我的问题是如何更改多边形内的栅格值。

Answer 1

编辑（两年后）：

此问题已按照 Hijmans 先生在评论中的承诺进行修复。
感谢您使用这种简单、很棒且方便的地理数据处理方式！

原答案：

终于，我自己找到了答案。

函数rasterize()可以帮助解决这个问题。

r <- raster(system.file("external/test.grd", package="raster"))
rst <- as.integer(stretch(r, 0, 10))
plot(rst)
rTpol <- rasterToPolygons(rst, dissolve = TRUE)
rTpol <- rTpol[rTpol$layer > 3,]
plot(rTpol)

rpol <- rasterize(rTpol, rst, field = 100, update = TRUE)
plot(rpol)

rasterize() 将 field = 值放入从多边形生成的栅格单元中，而 update = TRUE 获取栅格并使用字段值更新与多边形对应的所有单元。

结果如下所示：

使用空间多边形更改栅格值

Change raster values using spatial polygons

plot

r

raster

spatial

r-raster

编辑（两年后）：

原答案：