我有一堆具有自交集的多边形,这会导致在进一步处理它们时出现一些错误(特别是-我无法计算这些多边形与其他多边形的相交面积)。下面是一个破碎多边形的示例:
{
  "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抛出一个错误,指向交叉点坐标(我猜):=> "Geos::GEOSException: TopologyException: Input geom 0 is invalid: Self-intersection at or near point 6.8573510795579145 46.818650764080992 at 6.8573510795579145 46.818650764080992"。所以,至少我有它。
我的问题是:有没有办法自动修复那个交叉点?我读到,一种可能的解决方案是插入两个坐标为自交的类似点。但问题是-多边形有一个特定的顺序,我不知道在哪里插入这些点。
另外,也许有一些现有的工具可以帮助解决这个问题。

发布于 2018-02-21 23:22:33
我使用的解决方案是postgis's ST_MakeValid option for postgres,如果这是您的选择,您可以按照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等工具进行转换。总结一下这在幕后是如何工作的,它将把这些“领结”分成多个多边形,然后这些多边形将被制作成一个多多边形。相同类型的修复将应用于多个多边形。
https://stackoverflow.com/questions/48785494
复制相似问题