我正在使用shapely和python将多边形分成较小的部分:梯形,边与x轴平行的平行四边形和边与x轴平行的三角形。初始数据来自.gds文件,并以元组x,y坐标列表的形式显示。为了我的目的,我使用了这里描述的差异方法。但是当多边形有孔时,我会得到堆栈,例如:
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)
(请参阅标准手册)。
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
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句