Berechnen von Linienschnittpunkten auf einem Parallelogrammgitter ⇐ Python

[Anonymous]

• Ich muss ~2 Millionen Koordinaten berechnen (64-Bit-Float-Tupel)
• Sie sind die Schnittpunkte von Linien, die ein parallelogrammförmiges Gitter bilden (siehe Abbildung unten)
• Die Form des Parallelogramms ist bei jeder Iteration anders.
• Ich habe alle Grenzen Koordinaten (in Grün) stehen mir zur Verfügung.
• Die von Ost nach West verlaufenden Linien haben konstante Breitengrade.

Ich bin sicher, dass dies schon oft gelöst wurde, aber leider weiß ich nicht, nach welchen Schlüsselwörtern ich suchen sollte, um die Lösung zu finden. Logische Ansätze, die ich bisher in Betracht gezogen habe:

• Berechnen Sie die Steigungen der Nord-Süd-Linie (d. h.: (y2 – y1) / (x2 – x1)). Lösen Sie nach x auf (d. h.: Längengrad der Ost-West-Linie am Schnittpunkt).
• Berechnen Sie die Längen der Ost-West-Linie. Subtrahieren Sie die gesamten konstanten Breiten (d. h. 6x schmale Streifen). Teilen Sie den Rest durch 6. Wenden Sie das Ergebnis als Längsversatz für jeden Punkt auf der Ost-West-Linie an.
• Ich habe versucht, die Funktion „split()“ von Shapely zu verwenden, um das parallelogrammförmige Polygon mit Linienzeichenfolgen zu teilen, aber es wurde schnell klar, dass es ein Albtraum wäre, die zurückgegebene GeometryCollection zu verarbeiten und sie sortieren zu müssen.

Ich bin sicher, ihr wisst davon Allerdings viel bessere Möglichkeiten.
Fragen:

• Was ist der beste mathematische Berechnungsansatz, um dies zu erreichen? und
• Gibt es Python-Bibliotheken, mit denen ich die Dinge beschleunigen könnte?
• Bonus: Ich muss dies wahrscheinlich als Multiprozess ausführen, wenn die Zeit gekommen ist, den Datensatz zu generieren. Können Sie bestätigen, dass diese Art von Berechnungen am besten auf den CPU-Kernen und nicht auf den GPU-Kernen durchgeführt werden?