相交线以使用带有QGIS的Python进行交叉？

我有一组代表公交线路的线路。一些线重叠并且走相同的路。

我能够提取节点。

但是我对仅提取这样的交叉点感兴趣：

我怎样才能做到这一点？我正在寻找QGIS或Python的方法。

我尝试了GDAL Python 的交集方法，但这基本上只返回了顶点。

如果有两条线交叉，则QGIS的“ 线交叉点”方法会返回交叉点。但是，如果两条公交线路在同一条道路上走得很远，这并不能给我他们指向合并的位置。

节点：

您需要两件事，即折线的端点（无中间节点）和交点。还有一个问题，某些折线端点也是交点：

一个解决方案是使用Python以及Shapely和Fiona模块

1）读取shapefile：

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

2）找到线的终点（如何获得折线的终点？）：

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）计算相交（使用itertools模块通过层中的几何对进行迭代）。一些交点的结果是多点，我们需要一个点列表：

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）消除端点和交点之间的重复项（如您在图中看到的）

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

5）保存生成的shapefile

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

最后结果：

线段

如果你也想在节点之间的段，你需要“平面化”（平面图形，没有棱角相互交叉）的shape文件。这可以通过Shapely 的unary_union函数完成

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