如何在PyQGIS中连接线图层边缘至最近线图层并删除交点后多余线段

场景1：线边缘连接至最近线图层 #

GUI操作步骤（新手优先） #

• 提取待连接线路的端点：打开提取节点工具，选择需要连接的线图层，设置提取类型为仅端点，得到端点矢量图层。
• 生成最短连线：打开最近邻分析工具，输入刚才生成的端点图层，参考图层选择目标线图层，勾选生成最近点和创建连线，运行后直接得到端点到目标线的最短连线。
• 合并线路（可选）：如果需要把原线和连线合并成一个图层，用合并矢量图层工具，选择原线图层和连线图层即可。

补充：你之前尝试失败可能是没筛选仅端点，或者参考图层选错了，确保目标线图层是你要连接的那一个。

PyQGIS批量处理脚本 #

适合需要重复操作的场景，直接复制到QGIS Python控制台运行：

from qgis.core import QgsProject
import processing

# 替换成你的图层名称
source_layer_name = "待连接线图层"
target_layer_name = "目标线图层"

# 获取图层对象
source_lines = QgsProject.instance().mapLayersByName(source_layer_name)[0]
target_lines = QgsProject.instance().mapLayersByName(target_layer_name)[0]

# 提取仅端点
extract_result = processing.run(
    "native:extractvertices",
    {"INPUT": source_lines, "TYPE": 1, "OUTPUT": "memory:端点图层"}
)
endpoints_layer = extract_result["OUTPUT"]

# 生成端点到目标线的连线
nn_result = processing.run(
    "native:nearestneighbouranalysis",
    {
        "INPUT": endpoints_layer,
        "REFERENCE_LAYER": target_lines,
        "GENERATE_NEAREST_POINTS": True,
        "CREATE_CONNECTIONS": True,
        "OUTPUT": "memory:最短连线图层"
    }
)
connections_layer = nn_result["OUTPUT"]

# 将结果添加到画布
QgsProject.instance().addMapLayer(endpoints_layer)
QgsProject.instance().addMapLayer(connections_layer)

场景2：线边缘匹配对齐 #

GUI操作步骤（新手优先） #

• 提取两个线图层的端点：分别对两个线图层用提取节点工具，设置提取类型为仅端点，得到端点图层A（待匹配线的端点）和端点图层B（目标匹配线的端点）。
• 对齐端点：打开对齐要素工具，选择端点图层A作为输入，对齐目标选择端点图层B，对齐方式选最近点，运行后即可将A的端点移动到B的最近端点位置。
• 更新原线路：如果需要让原线的端点同步更新，用修改线工具手动拖动原线端点到对齐后的位置，或者用替换几何图形工具批量替换（需要先关联原线和端点的ID）。

简化批量方案 #

如果不需要保留原端点，直接用最近邻分析获取目标坐标后移动原线：

1. 用最近邻分析得到端点A到端点B的坐标映射。
2. 用字段计算器给原线图层添加目标起点X/Y和目标终点X/Y字段，关联端点的目标坐标。
3. 用PyQGIS脚本批量更新线的起止点：

from qgis.core import QgsProject, QgsGeometry, QgsPointXY

# 替换成你的图层名称
line_layer_name = "待匹配线图层"
line_layer = QgsProject.instance().mapLayersByName(line_layer_name)[0]

# 开启编辑
line_layer.startEditing()
for feature in line_layer.getFeatures():
    # 假设字段已添加：target_start_x, target_start_y, target_end_x, target_end_y
    start_point = QgsPointXY(feature["target_start_x"], feature["target_start_y"])
    end_point = QgsPointXY(feature["target_end_x"], feature["target_end_y"])
    new_geom = QgsGeometry.fromPolylineXY([start_point, end_point])
    line_layer.changeGeometry(feature.id(), new_geom)
# 保存编辑
line_layer.commitChanges()

内容的提问来源于stack exchange，提问作者Ayan