Intersection de lignes pour obtenir des croisements à l'aide de Python avec QGIS?

J'ai un ensemble de lignes représentant des lignes de bus. Certaines lignes se chevauchent et empruntent les mêmes routes.

Je peux extraire les nœuds.

Cependant, je suis intéressé à extraire uniquement les croisements comme celui-ci:

Comment puis-je faire ceci? Je cherche des moyens avec QGIS ou Python.

J'ai essayé la méthode d'intersection de GDAL Python mais cela ne me renvoie essentiellement que les sommets.

La méthode Line Intersections de QGIS me renvoie les croisements si deux lignes se croisent. Cependant, dans le cas où deux lignes de bus parcourent une grande partie de leur itinéraire sur la même route, cela ne me donne pas le point de fusionner.

Les nœuds:

Vous voulez deux choses, les points d'extrémité des polylignes (sans nœuds intermédiaires) et les points d'intersection. Il y a un problème supplémentaire, certains points d'extrémité de polylignes sont également des points d'intersection:

Une solution consiste à utiliser Python et les modules Shapely et Fiona

1) Lisez le fichier de formes:

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

2) Trouvez les points d'extrémité des lignes ( comment obtenir les points d'extrémité d'une polyligne? ):

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) Calculez les intersections (itération à travers des paires de géométries dans la couche avec le module itertools ). Le résultat de certaines intersections est MultiPoints et nous voulons une liste de points:

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) Éliminez les doublons entre les points d'extrémité et les points d'intersection (comme vous pouvez le voir sur les figures)

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

5) Enregistrez le fichier de formes résultant

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

Résultat final:

Les segments de ligne

Si vous souhaitez également les segments entre les nœuds, vous devez "planariser" ( graphique planaire , aucune arête ne se croise) votre fichier de formes. Cela peut être fait par la unary_union fonction de Shapely

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