在 R 中的空间对象中创建孔

Question

我需要在 spatialPolygonsDataFrame 中基于一个完全包含的较小多边形创建一个洞。

这说明了使用 WKT 创建的多边形的期望结果：

require(rgeos)
p = readWKT("POLYGON((1 1,5 1,5 5,1 5,1 1),(2 2,2 3,3 3,3 2,2 2))")
p = SpatialPolygonsDataFrame(p, data.frame(x=1))
plot(p, col='grey')

现在假设我们有以下没有任何孔几何形状的多边形，我们如何插入孔（即没有 readWKT）？我猜是 rgeos 函数，但我无法确定解决方案。

poly = readWKT("POLYGON((1 1,5 1,5 5,1 5,1 1))")
hole = readWKT("POLYGON((2 2,2 3,3 3,3 2,2 2))")

提前致谢。

Answer 1

library(rgeos)
gDifference(poly, hole)

应该这样做。