要更改 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 的栅格(左)和具有奇怪变化值的栅格(右):
因此,我的问题是如何更改多边形内的栅格值。
最佳答案
编辑(两年后):
这个问题已按照 Hijmans 先生在评论中的 promise 得到修复。
感谢您提供这种简单、出色且方便的地理数据处理方式!
原始答案:
最后,我自己找到了答案。
函数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
则获取栅格并更新所有与具有字段值的多边形相对应的单元格。
关于r - 使用空间多边形更改栅格值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/38609065/