如何在Python中从LineString序列获取6个等间距点？

核心思路 #

直接对坐标用np.linspace或np.arange做线性插值会出错——因为LineString的"等间距"是沿路径的长度分布，不是坐标的线性分布。正确做法是先计算路径的累计长度，再按总长度的比例（0%、20%、40%、60%、80%、100%，对应6个点）找到对应位置。

具体实现（基于Shapely库） #

假设你用Shapely处理几何数据（GIS场景常用库），以下是可直接运行的算法：

1. 计算单条LineString的总长度

   from shapely.geometry import LineString
   import numpy as np
   
   # 替换成你的LineString数据
   line = LineString([(0,0), (1,1), (2,0), (3,1), (4,0)])
   total_length = line.length
   

2. 生成等间距的距离节点
   要6个点，距离节点取总长度的0、0.2、0.4、0.6、0.8、1倍：

   distance_steps = np.linspace(0, total_length, num=6)
   

3. 沿路径插值获取对应点
   用interpolate方法根据距离获取路径上的坐标：

   equidistant_points = []
   for dist in distance_steps:
       point = line.interpolate(dist)
       equidistant_points.append((point.x, point.y))
   
   print("等间距点坐标：", equidistant_points)
   

错误原因说明 #

你之前的方法直接对坐标数组插值，相当于把x、y轴分别做线性分割，完全忽略了路径的弯曲长度——比如路径中某段折线的实际长度远大于坐标差，导致25%、75%这类中间分位点的位置严重偏离真实路径。

批量处理LineString序列 #

如果是多条LineString组成的序列，循环遍历处理即可：

line_strings = [LineString([(0,0), (1,2)]), LineString([(3,1), (5,4), (6,2)])]
all_equidistant_points = []

for line in line_strings:
    total_len = line.length
    dist_steps = np.linspace(0, total_len, 6)
    points = [(line.interpolate(d).x, line.interpolate(d).y) for d in dist_steps]
    all_equidistant_points.append(points)

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