如何用Python计算MULTILINESTRING中点与点之间的距离？

你当前得到的小数值是经纬度坐标的度数差——因为EPSG:4326是地理坐标系（用经纬度表示位置），直接调用Point.distance()计算的是球面的度数距离，并非实际物理距离，所以这些数值没有直观意义。以下是两种可行的解决方案：

方案一：转换到投影坐标系（推荐，精度更高） #

地理坐标系不适合直接计算距离，应先转换到投影坐标系（比如UTM坐标系，以米为单位）。UTM会根据数据所在位置自动匹配对应带号，geopandas可自动完成判断。

修正后的完整代码 #

import geopandas as gpd
from shapely.geometry import MultiLineString, Point

# 读取gdb文件
rivers_main = gpd.read_file("你的河流文件.gdb")

# 自动转换到UTM投影（单位为米）
rivers_main_utm = rivers_main.to_crs(rivers_main.estimate_utm_crs())

def calculate_segment_distances(line_geom):
    distance_values = []
    # 处理MultiLineString的所有子线段
    if isinstance(line_geom, MultiLineString):
        lines = list(line_geom.geoms)
    else:
        lines = [line_geom]
    
    for line in lines:
        points = [Point(xy) for xy in line.coords]
        for idx in range(len(points) - 1):
            point1 = points[idx]
            point2 = points[idx + 1]
            # 此时计算的距离单位为米
            distance = point1.distance(point2)
            distance_values.append(distance)
    return distance_values

# 计算第一条河流的各段距离
river_dists = calculate_segment_distances(rivers_main_utm.loc[0, 'geometry'])
for dist in river_dists:
    print(f"线段距离（米）: {dist:.2f}")

方案二：直接计算球面距离（无需转投影） #

如果不想转换投影，可使用geopy库的geodesic函数计算两点间的球面距离（基于WGS84椭球，精度可靠）。

代码示例 #

import geopandas as gpd
from shapely.geometry import MultiLineString
from geopy.distance import geodesic

# 读取gdb文件（保持EPSG:4326）
rivers_main = gpd.read_file("你的河流文件.gdb")

def calculate_geodesic_distances(line_geom):
    distance_values = []
    if isinstance(line_geom, MultiLineString):
        lines = list(line_geom.geoms)
    else:
        lines = [line_geom]
    
    for line in lines:
        coords = list(line.coords)
        for idx in range(len(coords) - 1):
            # geodesic接受(lat, lon)格式的坐标
            point1 = (coords[idx][1], coords[idx][0])
            point2 = (coords[idx+1][1], coords[idx+1][0])
            distance = geodesic(point1, point2).meters
            distance_values.append(distance)
    return distance_values

# 计算第一条河流的各段距离
river_dists = calculate_geodesic_distances(rivers_main.loc[0, 'geometry'])
for dist in river_dists:
    print(f"线段距离（米）: {dist:.2f}")

关键说明 #

• 原代码中point1和point2未定义，属于语法错误，已在上述代码中修正。
• 针对MultiLineString，原代码仅处理了第一个子线段，上述代码会遍历所有子线段并计算每一段的距离。
• UTM投影适合小范围区域计算，精度更高；geodesic适合跨区域或全球范围的场景。

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