rgeos::gDifference
を使用して、大きいポリゴンから小さなポリゴンを差し引くことで、これを行うことができます。
例:
まず、我々は我々が今、私たち
a1 = gArea(p1)
## [1] 320000
a2 = gArea(p2)
## [1] 150000
小さくなる見つけるために
gArea
を使用することができます
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)
に証明するためにいくつかのダミーポリゴンを作ります交差点をより大きなものから取り除くことができます
if (a1>a2){
p3 = gDifference(p1,p2)
p4 = p2
} else {
p3 = gDifference(p2,p1)
p4 = p1
}
plot(p3)
2つのポリゴンが一緒にこの
plot(p4, add=T)
次のようになります。このようにしてポリゴン
出典
2016-12-05 04:21:19
dww