QGIS ile Python kullanarak geçiş yapmak için kesişen çizgiler mi var?

Otobüs hatlarını temsil eden bir dizi hattım var. Bazı çizgiler üst üste biniyor ve aynı yolları kullanıyor.

Düğümleri çıkarabiliyorum.

Ancak ben sadece böyle geçişleri ayıklamak ilgileniyorum:

Bunu nasıl yapabilirim? QGIS veya Python ile yollar arıyorum.

GDAL Python'dan kesişim yöntemini denedim ama bu temelde bana sadece köşeleri döndürüyor.

Çizgi Intersections yöntemi iki satır çapraz eğer QGIS dan bana geçişleri döndürür. Ancak iki otobüs hattının aynı yol üzerindeki güzergahlarının bir kısmına gitmesi durumunda, bana birleştikleri yere işaret etmiyorlar.

Düğümler:

İki şey istiyorsunuz, çoklu çizgilerin bitiş noktaları (ara düğümler olmadan) ve kesişim noktaları. Ek bir sorun var, bazı çoklu hatların uç noktaları da kesişim noktalarıdır:

Çözüm Python ve Shapely ve Fiona modüllerini kullanmaktır

1) şekil dosyasını okuyun:

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

2) Çizgilerin uç noktalarını bulun (bir çoklu çizginin uç noktalarını nasıl elde edersiniz? ):

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) Kavşakları hesaplayın ( itertools modülü ile katmandaki geometri çiftleri üzerinden yineleme ). Bazı kavşakların sonucu MultiPoints'tir ve bir nokta listesi istiyoruz:

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) Bitiş noktaları ile kesişme noktaları arasındaki kopyaları ortadan kaldırın (şekillerde gördüğünüz gibi)

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

5) Ortaya çıkan şekil dosyasını kaydedin

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}})

Son sonuç:

Çizgi parçaları

Düğümler arasındaki segmentleri de istiyorsanız, şekil dosyanızı "planarize" ( Düzlemsel Grafik , hiçbir kenar birbirini geçmez ). Bu Shapely'nin unary_union işlevi tarafından yapılabilir

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