Weld segmen garis individual menjadi satu LineString menggunakan Shapely

Saya menggunakan Shapely in python dan saya diberi MultiLineStringdengan banyak Linestringobjek. Saya dapat menjamin bahwa semua LineStringobjek adalah garis-garis sederhana dengan hanya 2 simpul dan mereka semua adalah bagian dari satu garis tunggal (tanpa cabang).

Saya ingin "menghubungkan-titik" dan membuat satu LineString. Apakah saya perlu menulis metode pengelasan rekursif untuk ini atau apakah ada cara yang lebih cepat?

Anda dapat menggunakan shapely's ops.linemergeuntuk mencapai hal ini:

from shapely import geometry, ops

# create three lines
line_a = geometry.LineString([[0,0], [1,1]])
line_b = geometry.LineString([[1,1], [1,0]])
line_c = geometry.LineString([[1,0], [2,0]])

# combine them into a multi-linestring
multi_line = geometry.MultiLineString([line_a, line_b, line_c])
print(multi_line)  # prints MULTILINESTRING ((0 0, 1 1), (1 1, 2 2), (2 2, 3 3))

# you can now merge the lines
merged_line = ops.linemerge(multi_line)
print(merged_line)  # prints LINESTRING (0 0, 1 1, 2 2, 3 3)

# if your lines aren't contiguous
line_a = geometry.LineString([[0,0], [1,1]])
line_b = geometry.LineString([[1,1], [1,0]])
line_c = geometry.LineString([[2,0], [3,0]])

# combine them into a multi-linestring
multi_line = geometry.MultiLineString([line_a, line_b, line_c])
print(multi_line)  # prints MULTILINESTRING ((0 0, 1 1), (1 1, 1 0), (2 0, 3 0))

# note that it will now merge only the contiguous portions into a component of a new multi-linestring
merged_line = ops.linemerge(multi_line)
print(merged_line)  # prints MULTILINESTRING ((0 0, 1 1, 1 0), (2 0, 3 0))

Saya pikir Anda bisa melakukannya dengan Shapely menggunakan metode shapely.ops.linemerge ().

Sepertinya bisa mengambil daftar baris sebagai input dan menggabungkannya. Saya menggunakan metode 'polygonize' sebelumnya dan itu mengambil daftar garis.

Lihatlah dokumen di sini: http://toblerity.org/shapely/manual.html#shapely.ops.linemerge

shapely.ops.linemerge()gagal untuk beberapa baris saya jadi saya harus melakukannya secara manual. Tampaknya gagal untuk garis yang "kembali" ke dirinya sendiri, yaitu, melewati titik yang sama lebih dari sekali. Untuk kasus saya, saya tahu bahwa garis berada dalam urutan yang benar sehingga mudah untuk menulis fungsi kecil untuk menggabungkannya.

from shapely.geometry import LineString
from typing import List


def merge_lines(lines: List[LineString]) -> LineString:
    last = None
    points = []
    for line in merged_line:
        current = line.coords[0]

        if last is None:
            points.extend(line.coords)
        else:
            if last == current:
                points.extend(line.coords[1:])
            else:
                print('Skipping to merge {} {}'.format(last, current))
                return None
        last = line.coords[-1]
    return LineString(points)

Semoga ini bisa membantu seseorang

shapely.ops.linemergeberfungsi jika garis bersebelahan ("tips" bertepatan dengan "ekor" garis konstituen), tetapi jika mereka tidak bersebelahan (jika ada celah antara ujung dan ekor) ia mengembalikan MultiLineString lainnya. Jika garis konstituen Anda tertata dengan baik (dengan satu garis berakhir di dekat awal baris berikutnya) tetapi memiliki celah tip-to-tail, Anda dapat mengekstrak koordinat dan menggunakannya membuat garis sederhana baru. Pendekatan ini juga berfungsi untuk multi-garis yang dibuat dari subline yang lebih kompleks (yaitu subline dengan lebih dari dua poin).

import shapely

# Make a MultiLineString to use for the example
inlines = shapely.geometry.MultiLineString(
    [shapely.geometry.LineString([(0,0),(0,0.9)]), 
     shapely.geometry.LineString([(0,1),(1,1)])]
)

# Put the sub-line coordinates into a list of sublists
outcoords = [list(i.coords) for i in inlines]

# Flatten the list of sublists and use it to make a new line
outline = shapely.geometry.LineString([i for sublist in outcoords for i in sublist])