如何实现两个经纬度点之间的路径经纬度点插值

没有路线API可用时，可根据精度需求选择以下两种无依赖实现方式，最终都能输出包含起止点在内的完整插值经纬度数组。

核心实现思路 #

• 简单线性插值：直接对经纬度数值按比例线性拆分，计算量极小，适合两点间距10公里以内的场景，短距离下误差可忽略。
• 球面大圆插值：将地球近似为平均半径6371km的标准球体，沿两点间球面最短路径（大圆劣弧）做等距插值，中长距离下精度远高于线性插值，无需调用外部服务。

如果需要贴合实际道路的路径点，无公网API时可本地部署开源路由引擎实现，但如果只需要两点间直线路径点，以下方法完全满足最低要求。

预置示例坐标 #

# Google总部坐标
google_lat = 37.422131
google_lon = -122.084801

# Apple总部坐标
apple_lat = 37.33467267707233
apple_lon = -122.0089722675975

Python实现代码 #

import math

# 地球平均半径，单位：米
EARTH_RADIUS = 6371000

def linear_interpolate(lat1, lon1, lat2, lon2, num_points):
    """
    简单线性插值
    :param lat1/lon1: 起点经纬度
    :param lat2/lon2: 终点经纬度
    :param num_points: 总插值点数（包含起止点）
    :return: 插值点列表，格式为[(lat, lon), ...]
    """
    points = []
    for i in range(num_points):
        ratio = i / (num_points - 1)
        lat = lat1 + (lat2 - lat1) * ratio
        lon = lon1 + (lon2 - lon1) * ratio
        points.append((lat, lon))
    return points

def haversine_distance(lat1, lon1, lat2, lon2):
    """计算两点间球面距离，单位：米"""
    lat1_rad = math.radians(lat1)
    lat2_rad = math.radians(lat2)
    dlat = math.radians(lat2 - lat1)
    dlon = math.radians(lon2 - lon1)
    a = math.sin(dlat/2)**2 + math.cos(lat1_rad) * math.cos(lat2_rad) * math.sin(dlon/2)**2
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
    return EARTH_RADIUS * c

def spherical_interpolate(lat1, lon1, lat2, lon2, interval=1609):
    """
    球面大圆等距插值
    :param interval: 插值点间隔，单位：米，默认1609米即1英里
    :return: 插值点列表，格式为[(lat, lon), ...]
    """
    # 转弧度
    lat1_r, lon1_r = math.radians(lat1), math.radians(lon1)
    lat2_r, lon2_r = math.radians(lat2), math.radians(lon2)
    
    # 计算总距离和总插值点数
    total_dist = haversine_distance(lat1, lon1, lat2, lon2)
    num_points = math.ceil(total_dist / interval) + 1
    
    points = []
    # 计算大圆插值的圆心角
    d = total_dist / EARTH_RADIUS
    if d == 0:
        return [(lat1, lon1)]
    
    for i in range(num_points):
        frac = i * interval / total_dist
        if frac > 1:
            frac = 1
        # 球面线性插值公式
        a = math.sin((1 - frac) * d) / math.sin(d)
        b = math.sin(frac * d) / math.sin(d)
        
        x = a * math.cos(lat1_r) * math.cos(lon1_r) + b * math.cos(lat2_r) * math.cos(lon2_r)
        y = a * math.cos(lat1_r) * math.sin(lon1_r) + b * math.cos(lat2_r) * math.sin(lon2_r)
        z = a * math.sin(lat1_r) + b * math.sin(lat2_r)
        
        lat = math.degrees(math.atan2(z, math.sqrt(x**2 + y**2)))
        lon = math.degrees(math.atan2(y, x))
        points.append((lat, lon))
    return points

# 调用示例
if __name__ == "__main__":
    # 线性插值：总共生成10个点
    linear_points = linear_interpolate(google_lat, google_lon, apple_lat, apple_lon, 10)
    print("线性插值结果：", linear_points)
    
    # 球面插值：每隔1英里生成一个点
    spherical_points = spherical_interpolate(google_lat, google_lon, apple_lat, apple_lon, interval=1609)
    print("球面等距插值结果：", spherical_points)

使用说明 #

• 两点距离近、对精度要求不高时直接用线性插值即可，代码逻辑最简单。
• 距离超过10公里、需要等距间隔的路径点时，用球面插值方法，输出的点沿地球表面最短路径分布，误差在米级。
• 如果需要沿道路的路径点，没有公网API的情况下可以本地部署开源路由服务实现，不需要依赖第三方商业接口。

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