matlab 多边形 自相交,高效提取MultiPolygon中自相交特征生成的所有子多边形

从包含相当大数量(约20000)可能部分重叠的多边形的shapefile开始，我需要提取通过交叉它们的不同“边界”而产生的所有子多边形。

在实践中，从一些模拟数据开始：

library(tibble)

library(dplyr)

library(sf)

ncircles

rmax

x_limits

y_limits

set.seed(100)

xy

id = paste0("id_", 1:ncircles),

x = runif(ncircles, min(x_limits), max(x_limits)),

y = runif(ncircles, min(y_limits), max(y_limits))) %>%

as_tibble()

polys %

st_buffer(runif(ncircles, min = 1, max = 20))

plot(polys[1])

我需要派生一个包含ALL的sf或sp多边形，并且只需要交叉点生成的多边形，例如：

(请注意，颜色仅用于举例说明预期结果，其中每个“不同颜色”区域是一个单独的多边形，不覆盖任何其他多边形)

我知道我可以通过一次分析一个多边形，识别并保存所有交叉点然后“擦除”那些区域形成完整的多边形并继续循环来解决问题，但这很慢。

我觉得应该有一个更有效的解决方案，但我无法弄清楚，所以任何帮助将不胜感激！

(欢迎基于sf和sp的解决方案)

更新：

最后，我发现即使是“一次一个多边形”，任务也远非简单！我真的在努力解决这个显而易见的“简单”问题！任何提示？即使是一个缓慢的解决方案或提示开始正确的道路将不胜感激！

更新2：

也许这会澄清一些事情：所需的功能类似于这里描述的功能：

更新3：

我将赏金奖给了@ shuiping-chen(谢谢！)，他的回答正确地解决了所提供的示例数据集上的问题。然而，“方法”一般被推广到“四倍”或“n-uple”交叉点可能的情况。如果我管理的话，我会在未来几天尝试解决这个问题并发布一个更通用的解决方案！