RGeo - 修复自相交

Question

我有一堆具有自相交的多边形，这在进一步后处理它们时会导致一些错误（特别是 - 我无法计算这些多边形与其他多边形的相交面积）。这是一个破碎多边形的例子：

{
  "type": "MultiPolygon",
  "coordinates": [
    [
      [
        [
          6.881057785381658,
          46.82373306675715
        ],
        [
          6.857171686909481,
          46.81861230543794
        ],
        [
          6.857354659059071,
          46.81856788926046
        ],
        [
          6.856993473052509,
          46.82693029065604
        ],
        [
          6.8612894138116785,
          46.83422796373707
        ],
        [
          6.86720955648855,
          46.835636765630476
        ],
        [
          6.871281147359957,
          46.83078486366309
        ],
        [
          6.871573291317274,
          46.8306215963777
        ],
        [
          6.877608228639841,
          46.82771553607934
        ],
        [
          6.877758462659651,
          46.82772313420989
        ],
        [
          6.877852632482749,
          46.827735617670285
        ],
        [
          6.880928107931434,
          46.82630213148064
        ],
        [
          6.8810399979122305,
          46.82622029042867
        ],
        [
          6.881117606743071,
          46.826115612819855
        ],
        [
          6.881057785381658,
          46.82373306675715
        ]
      ]
    ]
  ]
}

这就是地图上的样子 - 如您所见，两条多边形边相交。RGeo引发错误，指向交点坐标（I猜）： => "Geos::GEOSException: TopologyException: Input geom 0 is invalid: Self-intersection at or near point 6.8573510795579145 46.818650764080992 at 6.8573510795579145 46.818650764080992"。所以，我至少拥有它。

我的问题是：有没有办法自动修复那个交叉点？我读到，一个可能的解决方案是插入 2 个具有自交坐标的相似点。但问题是 - 多边形有特定的顺序，我不知道在哪里插入这些点。

另外，也许有一些现有的工具可以帮助解决这个问题......

Answer 1

我将使用的解决方案是postgis 的 postgres ST_MakeValid 选项，如果这是您的一个选项，您可以按照以下方式执行某些操作ST_AsText(ST_MakeValid(geom_column))，或者如果您宁愿传递文本，这里是一个使用prepair中显示的领结示例的示例：

select ST_AsText(ST_MakeValid(ST_GeomFromText('POLYGON((0 0, 0 10, 10 0, 10 10, 0 0))')));
st_astext                         
-----------------------------------------------------------
MULTIPOLYGON(((0 0,0 10,5 5,0 0)),((5 5,10 10,10 0,5 5)))
(1 row)

如果您对此不感兴趣，您可以导出这些几何图形并使用prepair等工具来转换它们。总结一下它在幕后的工作原理，它将把这些“领结”分割成多个多边形，然后将它们制成一个多边形。相同类型的修复将应用于多边形。