将多边形复制区域(重叠)分配给仅一个多边形

时间:2016-12-05 01:50:32

标签: r overlap polygons

我有一个shapefile(可以在这里下载http://www.car.gov.br/publico/municipios/downloads?sigla=MA(任何人都可以)),其中每个多边形代表一个农村私有财产,因为所有者提交了它。但是,多边形重叠。对于我需要运行的分析,一个区域不能被计算两次(即,在两个属性中)。因此,在重叠区域,我想将这个重复区域分配给面积最小的多边形,并将其从另一个区域中删除。

我查看了其他问题,例如this one。但是没有人能真正解决这个问题。

1 个答案:

答案 0 :(得分:3)

您可以使用rgeos::gDifference从较大的多边形中减去较小的多边形来完成此操作。

一个例子:

首先,我们制作一些虚拟多边形来演示

library(sp)
library(rgeos)
p1 <- Polygon(matrix(c(
  100, 100,
  100, 500,
  900, 500,
  900, 100,
  100, 100), ncol=2, byrow = T))
p2 <- Polygon(matrix(c(
  50, 400,
  50, 600,
  800, 600,
  800, 400,
  50, 400), ncol=2, byrow = T))
p1 <- SpatialPolygons(list(Polygons(list(p1), "p1")))
p2 <- SpatialPolygons(list(Polygons(list(p2), "p2")))
plot(p1)
plot(p2, add=T)

enter image description here

我们可以使用gArea来查找哪个更小

a1 = gArea(p1)
## [1] 320000
a2 = gArea(p2)
## [1] 150000

现在我们可以删除较大多边形的交点,如下所示:

if (a1>a2){
  p3 = gDifference(p1,p2)
  p4 = p2
} else {
  p3 = gDifference(p2,p1)
  p4 = p1
}
plot(p3)

enter image description here

两个多边形一起看起来像这样

plot(p4, add=T)

enter image description here