¿Intersección de líneas para obtener cruces usando Python con QGIS?

Tengo un conjunto de líneas que representan líneas de autobús. Algunas de las líneas se superponen y toman los mismos caminos.

Soy capaz de extraer los nodos.

Sin embargo, estoy interesado en extraer solo cruces como este:

¿Cómo puedo hacer esto? Estoy buscando formas con QGIS o Python.

Probé el método de intersección de GDAL Python, pero esto básicamente me devuelve solo los vértices.

El método de intersecciones de línea de QGIS me devuelve los cruces si se cruzan dos líneas. Sin embargo, en el caso de que dos líneas de autobús recorran gran parte de su ruta en el mismo camino, no me da que señalen dónde se fusionan.

Los nodos:

Desea dos cosas, los puntos finales de las polilíneas (sin nodos intermedios) y los puntos de intersección. Hay un problema adicional, algunos puntos finales de polilíneas también son puntos de intersección:

Una solución es utilizar Python y los módulos Shapely y Fiona.

1) Lea el archivo de forma:

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

2) Encuentre los puntos finales de las líneas ( ¿cómo se obtendrían los puntos finales de una polilínea? ):

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) Calcule las intersecciones (iterando a través de pares de geometrías en la capa con el módulo itertools ). El resultado de algunas intersecciones son MultiPoints y queremos una lista de puntos:

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) Elimine los duplicados entre los puntos finales y los puntos de intersección (como puede ver en las figuras)

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

5) Guarde el archivo de forma resultante

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

Resultado final:

Los segmentos de línea

Si desea también los segmentos entre los nodos, debe "planarizar" ( Gráfico plano , sin bordes cruzados) su archivo de forma. Esto se puede hacer por el unary_union función del bien proporcionado

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