Memotong garis untuk mendapatkan persimpangan menggunakan Python dengan QGIS?

Saya memiliki satu set garis yang mewakili jalur bus. Beberapa jalur tumpang tindih dan mengambil jalan yang sama.

Saya dapat mengekstrak node.

Namun saya tertarik mengekstraksi penyeberangan seperti ini:

Bagaimana saya bisa melakukan ini? Saya mencari cara dengan QGIS atau Python.

Saya mencoba metode persimpangan dari GDAL Python tetapi ini pada dasarnya mengembalikan saya hanya simpul.

The Jalur Persimpangan Metode dari QGIS kembali saya penyeberangan jika dua garis silang. Namun dalam kasus bahwa dua jalur bus pergi jauh dari rute mereka di jalan yang sama, itu tidak memberi saya titik di mana mereka bergabung.

Simpul:

Anda menginginkan dua hal, titik akhir dari polyline (tanpa titik tengah) dan titik persimpangan. Ada masalah tambahan, beberapa titik akhir polyline juga merupakan titik persimpangan:

Solusinya adalah dengan menggunakan Python dan modul Shapely dan Fiona

1) Baca shapefile:

from shapely.geometry import Point, shape
import fiona
lines = [shape(line['geometry']) for line in fiona.open("your_shapefile.shp")]

2) Temukan Poin akhir dari garis ( bagaimana cara mendapatkan titik akhir dari polyline? ):

endpts = [(Point(list(line.coords)[0]), Point(list(line.coords)[-1])) for line  in lines]
# flatten the resulting list to a simple list of points
endpts= [pt for sublist in endpts  for pt in sublist] 

3) Hitung persimpangan (iterasi melalui pasangan geometri di lapisan dengan modul itertools ). Hasil dari beberapa persimpangan adalah MultiPoints dan kami ingin daftar poin:

import itertools
inters = []
for line1,line2 in  itertools.combinations(lines, 2):
  if  line1.intersects(line2):
    inter = line1.intersection(line2)
    if "Point" == inter.type:
        inters.append(inter)
    elif "MultiPoint" == inter.type:
        inters.extend([pt for pt in inter])
    elif "MultiLineString" == inter.type:
        multiLine = [line for line in inter]
        first_coords = multiLine[0].coords[0]
        last_coords = multiLine[len(multiLine)-1].coords[1]
        inters.append(Point(first_coords[0], first_coords[1]))
        inters.append(Point(last_coords[0], last_coords[1]))
    elif "GeometryCollection" == inter.type:
        for geom in inter:
            if "Point" == geom.type:
                inters.append(geom)
            elif "MultiPoint" == geom.type:
                inters.extend([pt for pt in geom])
            elif "MultiLineString" == geom.type:
                multiLine = [line for line in geom]
                first_coords = multiLine[0].coords[0]
                last_coords = multiLine[len(multiLine)-1].coords[1]
                inters.append(Point(first_coords[0], first_coords[1]))
                inters.append(Point(last_coords[0], last_coords[1]))

4) Hilangkan duplikat antara titik akhir dan titik persimpangan (seperti yang Anda lihat pada gambar)

result = endpts.extend([pt for pt in inters  if pt not in endpts])

5) Simpan shapefile yang dihasilkan

from shapely.geometry import mapping
# schema of the shapefile
schema = {'geometry': 'Point','properties': {'test': 'int'}}
# creation of the shapefile
with fiona.open('result.shp','w','ESRI Shapefile', schema) as output:
    for i, pt in enumerate(result):
        output.write({'geometry':mapping(pt), 'properties':{'test':i}})

Hasil akhir:

Segmen garis

Jika Anda ingin juga segmen antara node, Anda perlu "planarize" ( Grafik Planar , tidak ada ujung yang saling bersilangan) shapefile Anda. Ini dapat dilakukan dengan fungsi unary_union dari Shapely

from shapely.ops import unary_union
graph = unary_union(lines)