将多边形重叠部分的面积分配给一个多边形

Question

我有一个shapefile文件（可在此处下载：http://www.car.gov.br/publico/municipios/downloads?sigla=MA），其中每个多边形代表一个农村私有财产，由业主提交。然而，这些多边形存在重叠。为了进行分析，一个区域不能被计算两次（即在两个财产中都存在）。因此，在存在重叠的区域，我想将该重复区域分配给面积最小的多边形，并从其他多边形中删除它。

我查看了其他类似问题，例如 这个 ，但没有真正解决这个问题。

- Liz

Answer 1

你可以使用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)

我们可以使用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)

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

plot(p4, add=T)

- dww