行内のポリゴンを有効なポリゴンに変換する方法は？

Question

多角形をより小さな部分に分割するために、shapelyとpythonを使用しています。台形、x軸に平行な辺を持つ平行四辺形、x軸に平行な辺を持つ三角形です。初期データは.gdsファイルから得られ、タプルx、y座標のリストのように表示されます。私の目的のために、私はhereのような違いの方法を使いました。ポリゴンが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 

そこで問題は、ポリゴンの正しい作成のための船体と穴にポイントのリストを変換する方法を最も単純な方法とは何かありますか？

ありがとうございました！

Answer 1

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