QGIS中用GRASS v.voronoi.skeleton生成复杂几何中心线的问题

针对你遇到的问题——简单几何用v.voronoi.skeleton生成中心线效果良好，但复杂分叉几何会出现Voronoi剖分偏向单侧、仅生成单路径中心线的情况，可通过以下方案解决：

一、预处理：智能分割复杂几何 #

将多路径的复杂几何拆分为独立单路径要素，是最直接的解决思路：

• 基于初步中轴线的自动分割：先用QGIS「提取中心线」工具（Processing工具箱 → 矢量几何工具 → 提取中心线）生成初步中轴线，从中提取分叉节点；再用「分割线」工具，以这些节点为分割点，将原几何拆分为多个独立的单路径要素，之后对每个要素单独运行v.voronoi.skeleton。
• 拓扑检测分割：若原数据是面要素，先转线（「多边形转线」），再用「拓扑检查器」找出分叉/自相交节点，手动或自动分割为单分支要素。批量处理可借助Python脚本结合shapely库自动检测分叉点：

  from shapely.geometry import LineString, Point, MultiPoint
  from shapely.ops import split
  
  def split_line_at_junctions(line):
      junctions = []
      coords = list(line.coords)
      for i in range(1, len(coords)-1):
          # 通过向量叉积判断夹角，识别分叉点
          prev_vec = (coords[i][0]-coords[i-1][0], coords[i][1]-coords[i-1][1])
          next_vec = (coords[i+1][0]-coords[i][0], coords[i+1][1]-coords[i][1])
          # 叉积绝对值大于阈值则判定为分叉，可根据数据精度调整阈值
          if abs(prev_vec[0]*next_vec[1] - prev_vec[1]*next_vec[0]) > 0.1:
              junctions.append(coords[i])
      if junctions:
          split_points = MultiPoint([Point(p) for p in junctions])
          return split(line, split_points)
      return [line]
  

• 自动化流程：用QGIS「图形模型器」把「分割→中心线生成」的流程串起来，实现批量处理。

二、调整GRASS工具参数 #

无需分割几何，通过调整v.voronoi.skeleton或换用GRASS其他工具优化结果：

• 优化v.voronoi.skeleton阈值：减小threshold参数（距离阈值），默认值可能过大导致分叉路径被合并。尝试将阈值从0.5调整为0.1左右，让工具保留更多细节，识别出两条独立路径的骨架。
• 换用v.skeleton工具：GRASS的v.skeleton专门针对多边形生成骨架线，支持prune参数（开启后修剪冗余分支），搭配prune_threshold调整修剪强度，能更好地保留分叉路径的中心线。

三、替代工具方案 #

如果上述方法效果不佳，可尝试更适配复杂几何的工具：

• QGIS原生「提取中心线」工具：基于CGAL库开发，对分叉路径的支持更智能，无需手动分割即可生成多分支中心线。调整「偏移」和「简化容差」参数适配你的数据精度。
• PostGIS函数：将数据导入PostGIS后，使用ST_ApproximateMedialAxis函数生成中轴线，该函数对复杂多边形的分叉处理逻辑更完善：

  SELECT ST_ApproximateMedialAxis(geom) AS centerline FROM your_dataset;
  

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