如何计算R中多边形重叠的面积？

Question

有谁知道如何计算R中2个或更多个多边形之间的共同区域？我希望这种计算的输出是该重叠区域的新多边形的坐标。 干杯

Answer 1

编辑：这些天我会在'raster'包中使用'intersect'，'cover'，'erase'，'union'和相关函数。他们努力工作以保持顶级对象和属性。

ORIG： 您可以将rgeos包与其gIntersection函数一起使用。成对和交叉路口之间的连续呼叫将使您到达那里。见

library(rgeos)
?gIntersection

您需要进入sp包中的“SpatialPolygons”结构以获取最终坐标。看小插图（“sp”）。

Answer 2

我想我会添加我最终使用的解决方案 - joinPolys包中的PBSmapping函数。

实施例

library(PBSmapping)
p1 <- data.frame(PID=rep(1, 4), POS=1:4, X=c(1,1,6,6), Y=c(1,3,3,1))
p2 <- data.frame(PID=rep(2, 5), POS=1:5, X=c(4,4,8,8,6), Y=c(2,4,4,2,1))
p3 <- joinPolys(p1,p2)
x11()
par(mar=c(3,3,1,1))
plot(1,1,ylim=c(0,5),xlim=c(0,9), t="n", xlab="", ylab="")
polygon(p1$X, p1$Y, border=2)
polygon(p2$X, p2$Y)
polygon(p3$X, p3$Y, col=rgb(0,0,1,0.2))