Fix ungültige polygon in Formschöne

Formschöne definiert ein Polygon als ungültig, wenn jedes seiner Segmente überschneiden, einschließlich der Segmente, sind colinear. Viele software-Pakete erstellen Sie eine region oder einen Bereich mit einer "Aussparung", wie hier gezeigt, die colinear Segmente:

Fix ungültige polygon in Formschöne

>>> pp = Polygon([(0,0), (0,3), (3,3), (3,0), (2,0), 
                  (2,2), (1,2), (1,1), (2,1), (2,0), (0,0)])
>>> pp.is_valid
WARNING:shapely.geos:Self-intersection at or near point 2 0
False

Natürlich, den Ausschnitt, die implementiert werden können, die nativ in Formschönen, oder dieselbe geometrie umgesetzt werden können, wie zwei gültige Polygone, aber wenn ich nur die Liste der Punkte, die oben gezeigt, gibt es eine leicht zu "beheben" dieses (erstellen von geometrie aus dieser Liste der Punkte)?

InformationsquelleAutor jpcgt | 2013-12-30
  1. 37

    Ich habe eine Lösung gefunden, die funktioniert für den speziellen Fall gegeben:

    >>> pp2 = pp.buffer(0)
>>> pp2.is_valid
True
>>> pp2.exterior.coords[:]
[(0.0, 0.0), (0.0, 3.0), (3.0, 3.0), (3.0, 0.0), (2.0, 0.0), (0.0, 0.0)]
>>> pp2.interiors[0].coords[:]
[(2.0, 1.0), (2.0, 2.0), (1.0, 2.0), (1.0, 1.0), (2.0, 1.0)]
    • Ich finde, dass dieser buffer( 0 ) - Methode funktioniert im Allgemeinen für die Befestigung Polygone mit deckungsgleiche Linien unter sub-Polygone. Dies ist sehr nützlich. Kennen Sie einen Ort, wo dieser trick ist offiziell empfohlen/sanktioniert?
    • Es ist in der Anleitung! Aber nicht angekündigte überhaupt.
    • Durch diese Methode nicht behoben mein Problem: ich hatte ein paar selbst schneidende Multipolygone.
    • Sie müssen ein wenig vorsichtig mit der buffer(0) Technik. Wir hatten bowtie Fällen, in denen es zerstört den größten Teil des Polygons und linke nur eine kleine bowtied Ecke. YMMV.
    InformationsquelleAutor jpcgt
Schreibe einen Kommentar