基于输入距离计算两点间经纬度坐标的函数问题排查及求解

一、如何找到距起点恰好300km处的坐标 #

要计算起点沿两点连线方向移动指定距离后的坐标，需基于**大圆航线（Great Circle）**的数学逻辑：

1. 先计算起点到终点的方位角（即两点连线的方向角）；
2. 以起点为原点，沿该方位角移动指定距离，通过球面三角公式推导新的经纬度。

二、你的Python函数问题分析 #

函数输出结果始终一致，核心问题有两个：

1. 未使用传入的distance_km参数：函数计算中间点时用的是两点间的总距离distance，而非用户传入的目标距离distance_km，导致无论传什么值，计算逻辑都固定；
2. 中间点坐标计算公式错误：原公式完全不符合球面移动的推导逻辑，正确逻辑应基于方位角和目标移动距离来计算。

修正后的函数代码 #

import math

def find_point(lon1, lat1, lon2, lat2, distance_km):
    # 地球平均半径（单位：km）
    EARTH_RADIUS = 6371.0
    
    # 将经纬度转换为弧度
    lat1_rad = math.radians(lat1)
    lon1_rad = math.radians(lon1)
    lat2_rad = math.radians(lat2)
    lon2_rad = math.radians(lon2)
    
    # 计算起点到终点的方位角（弧度）
    delta_lon = lon2_rad - lon1_rad
    y = math.sin(delta_lon) * math.cos(lat2_rad)
    x = math.cos(lat1_rad)*math.sin(lat2_rad) - math.sin(lat1_rad)*math.cos(lat2_rad)*math.cos(delta_lon)
    bearing = math.atan2(y, x)
    
    # 计算目标距离对应的弧度值
    angular_distance = distance_km / EARTH_RADIUS
    
    # 计算目标点的纬度（弧度）
    target_lat_rad = math.asin(
        math.sin(lat1_rad)*math.cos(angular_distance) + 
        math.cos(lat1_rad)*math.sin(angular_distance)*math.cos(bearing)
    )
    
    # 计算目标点的经度（弧度）
    target_lon_rad = lon1_rad + math.atan2(
        math.sin(bearing)*math.sin(angular_distance)*math.cos(lat1_rad),
        math.cos(angular_distance) - math.sin(lat1_rad)*math.sin(target_lat_rad)
    )
    
    # 转换回十进制度数
    target_lat = math.degrees(target_lat_rad)
    target_lon = math.degrees(target_lon_rad)
    
    return target_lon, target_lat

代码说明 #

• 先正确计算起点到终点的方位角，确保移动方向与两点连线一致；
• 使用传入的distance_km计算对应的弧度距离，让目标点位置随参数变化；
• 采用标准球面三角公式推导目标点经纬度，完全符合大圆航线的移动逻辑；
• 针对你提供的示例坐标，调用find_point(16.36920394, 48.20342942, 35.93935074, 31.95158219, 300)即可得到距起点300km处的正确坐标。

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