Python 如何将直线中与自身接触的多边形转换为有效多边形?
我使用shapely和python将多边形划分为更小的部分:梯形、边平行于x轴的平行四边形和边平行于x轴的三角形。初始数据来自.gds文件,如元组x、y坐标列表所示。出于我的目的,我使用了上述的差分方法。但当多边形有以下情况时,我会得到堆栈,例如:Python 如何将直线中与自身接触的多边形转换为有效多边形?,python,polygon,shapely,Python,Polygon,Shapely,我使用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,
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)
(请参阅)
谢谢你的快速回复!在我写这篇文章之前,我已经读过了,但是我对我通常检查多边形
x,y=shape.external.coords.xy的方法感到困惑,当然没有内部点。所以你清除了我的思想。谢谢
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