如何在谷歌地图经纬度线段上获取指定数量的点？

嘿，这个问题我之前做地图相关开发时也踩过坑！平面几何的直线公式在经纬度坐标系里完全行不通，毕竟地球是个近似球体，得用球面插值的思路来解决，我给你梳理下具体方法和思路：

核心思路 #

经纬度本质是球面坐标系的坐标，所以我们不能直接在二维平面上插值，得把坐标转换到三维空间（地球的球面）上做线性插值，再转换回经纬度；或者直接用球面几何的大圆插值公式来计算中间点。

具体实现步骤（球面插值法） #

1. 把经纬度转换成弧度 #

所有球面计算的公式基本都基于弧度，转换公式很简单：

弧度 = 角度 × π / 180

比如把起点(lat1, lon1)和终点(lat2, lon2)都转换成弧度值lat1_rad, lon1_rad, lat2_rad, lon2_rad。

2. 转换为三维笛卡尔坐标 #

把球面坐标转成三维空间中的点，公式如下（R是地球半径，一般取6371km，要是只需要相对位置，R取1也不影响最终结果）：

• x = R × cos(lat_rad) × cos(lon_rad)
• y = R × cos(lat_rad) × sin(lon_rad)
• z = R × sin(lat_rad)

3. 三维空间线性插值+归一化 #

假设要取X个点（包含起点和终点），那么每个点的插值参数t从0到1均匀分布：t = i/(X-1)，其中i从0到X-1（t=0对应起点，t=1对应终点）。

先对三维坐标做线性插值：

• x_interp = x1 + t*(x2 - x1)
• y_interp = y1 + t*(y2 - y1)
• z_interp = z1 + t*(z2 - z1)

插值后的点可能不在球面上，所以要归一化：把(x_interp, y_interp, z_interp)除以它的模长，得到单位球面上的点，再乘以R（如果之前用了R的话）。

4. 转换回经纬度 #

把归一化后的三维坐标转回经纬度，公式：

• lat_rad = arcsin(z / R)
• lon_rad = arctan2(y, x)
  再把弧度转成角度：角度 = 弧度 × 180 / π，就得到了插值点的经纬度。

伪代码示例（Python） #

import math

def get_interpolated_points(lat1, lon1, lat2, lon2, num_points):
    # 1. 转弧度
    lat1_rad = math.radians(lat1)
    lon1_rad = math.radians(lon1)
    lat2_rad = math.radians(lat2)
    lon2_rad = math.radians(lon2)
    
    # 2. 转三维笛卡尔坐标
    earth_radius = 6371000  # 地球半径，单位米
    x1 = earth_radius * math.cos(lat1_rad) * math.cos(lon1_rad)
    y1 = earth_radius * math.cos(lat1_rad) * math.sin(lon1_rad)
    z1 = earth_radius * math.sin(lat1_rad)
    
    x2 = earth_radius * math.cos(lat2_rad) * math.cos(lon2_rad)
    y2 = earth_radius * math.cos(lat2_rad) * math.sin(lon2_rad)
    z2 = earth_radius * math.sin(lat2_rad)
    
    interpolated_points = []
    for i in range(num_points):
        t = i / (num_points - 1)
        # 3. 线性插值
        x = x1 + t * (x2 - x1)
        y = y1 + t * (y2 - y1)
        z = z1 + t * (z2 - z1)
        
        # 归一化到球面上
        norm = math.sqrt(x**2 + y**2 + z**2)
        x /= norm
        y /= norm
        z /= norm
        
        # 4. 转回经纬度
        lat_rad = math.asin(z)
        lon_rad = math.atan2(y, x)
        
        lat = math.degrees(lat_rad)
        lon = math.degrees(lon_rad)
        interpolated_points.append((lat, lon))
    
    return interpolated_points

额外说明 #

• 如果需要更高精度（比如考虑地球是椭球体而非完美球面），可以使用文森特椭圆体插值公式，不过对于大多数地图应用场景，球面插值的精度已经足够。
• 要是两点是对跖点（比如北极和南极附近），插值路径会有歧义，这时候需要指定具体的行进方向（比如沿东经还是西经）。

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