如何将在一条线上接触自身的多边形转换为有效的多边形?

How to convert polygon that touches itself in a line into valid polygon?

我正在使用 shapely 和 python 将多边形分成更小的部分:梯形、边平行于 x 轴的平行四边形和边平行于 x 轴的三角形。初始数据来自 .gds 文件,并以 x,y 坐标元组列表的形式呈现。为了我的目的,我使用了 . But I get stack when polygon has holes 中描述的差异方法,例如:

from shapely.geometry import Polygon
points = [(0.0, -1.0), (0.0, 2.0), (3.0, 2.0), (3.0, 1.0), (1.0, 1.0), (1.0, 0.0), (2.0, 0.0), (2.0, 1.0), (3.0, 1.0), (3.0,-1.0)]
poly = Polygon(points)
print(poly.is_valid)
#Self-intersection at or near point 2 1
#False

所以问题是如何将此点列表转换为外壳和孔以正确创建多边形的最简单方法是什么?

谢谢!

标准技巧是使用 poly.buffer(0)(参见 the shapely manual)。

polyb = poly.buffer(0)
print(shapely.geometry.mapping(polyb))
# {'type': 'Polygon', 'coordinates': (((0.0, -1.0), (0.0, 2.0), (3.0, 2.0), (3.0, 1.0), (3.0, -1.0), (0.0, -1.0)), ((2.0, 1.0), (1.0, 1.0), (1.0, 0.0), (2.0, 0.0), (2.0, 1.0)))}
print(polyb.is_valid)
# True