R计算重叠部分（多边形交叉点）的快速方式

Question

我正在寻找一种方法来计算基础多边形被另一个（选择）多边形覆盖的百分比。

以此为例。我想计算红色区域的红色和蓝色区域。

group_a <- data.frame(x = c(1, 4, 1),
                      y = c(2, 4, 4),
                      group = "base")
group_b <- data.frame(x = c(2, 5, 2),
                      y = c(3, 3, 5),
                      group = "selection")

dat <- rbind(group_a, group_b)
dat
x y     group
1 1 2      base
2 4 4      base
3 1 4      base
4 2 3 selection
5 5 3 selection
6 2 5 selection

library(ggplot2)

ggplot(dat, aes(x = x, y = y, fill = group)) + geom_polygon(alpha = 0.5)

一种方法是使用rgdal - 和rgeos - 包将data.frames更改为Polygon，然后使用gArea和gIntersection查找必要的价值观。但是，我认为这可能不是最快的方法。你知道更快的方法来计算区域/交叉点吗？如果有RCpp方法，我会特别高兴。

我希望我没有问过多;和往常一样：感谢您的任何想法/提示！

编辑：这就是我使用其他软件包的方式：

library(rgdal)
library(rgeos)
base <- group_a[, c("x", "y")]
sel <- group_b[, c("x", "y")]

base_pol <- Polygons(list(Polygon(base)), "base")
sel_pol <- Polygons(list(Polygon(sel)), "sel")

shape <- SpatialPolygons(list(base_pol, sel_pol))
plot(shape)

base_area <- gArea(shape["base"])
intersections <- gIntersection(shape["base"], shape["sel"])

gArea(intersections)/base_area
# [1] 0.4027778

加速基准测试：正在运行microbenchmark，我获得了以下结果：

OvPerc <- function(base, sel) {
  base <- base[, c("x", "y")]
  sel <- sel[, c("x", "y")]

  base_pol <- Polygons(list(Polygon(base)), "base")
  sel_pol <- Polygons(list(Polygon(sel)), "sel")

  shape <- SpatialPolygons(list(base_pol, sel_pol))
  # plot(shape)

  base_area <- gArea(shape["base"])
  intersections <- gIntersection(shape["base"], shape["sel"])

  return(gArea(intersections)/base_area)
}

library(microbenchmark)
microbenchmark(OvPerc(group_a, group_b), times = 1000)
# Unit: milliseconds
#                    expr      min       lq     mean   median       uq      max neval
# OvPerc(group_a, group_b) 3.680386 3.970394 4.932217 4.157952 4.745398 67.13696  1000