如何计算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))