基于两点、距离与角度计算第三点经纬度的技术求助

你的核心问题是生成Point C的逻辑错误：当前代码直接将A→B的方位角与指定角度相加，混淆了「地理方位角」和「两点连线夹角」的概念，导致最终计算出的∠ABC与预期不符。

比如你传入105度时，代码计算的是「从正北顺时针出发，先走到A→B的方向，再转105度的方位」，但你实际需要的是「在Point B处，BC与BA连线形成105度夹角」的方向，二者的角度基准完全不同。

正确思路是：

1. 计算B→A的方位角（即从Point B指向Point A的地理方位，等于A→B的方位角加180度）
2. 以B→A的方位角为基准，加上/减去指定角度，得到Point C的目标方位角
3. 从Point B出发，沿目标方位角移动指定距离，得到Point C

以下是修正后的代码：

from geopy import Point
from geopy.distance import geodesic
import math

def extend_angle_line(pointA, pointB, distance, direction=1, angle=90):
    # 计算A→B的方位角
    ab_bearing = calculate_bearing(pointA, pointB)
    # 计算B→A的方位角（作为夹角的基准方向）
    ba_bearing = (ab_bearing + 180) % 360
    # 计算BC的目标方位角：direction=1为逆时针转angle，direction=-1为顺时针转angle
    bc_bearing = (ba_bearing + angle * direction) % 360
    # 从B出发沿目标方位移动指定距离
    extended_point = geodesic().destination(pointB, bc_bearing, distance)
    # 验证角度（可选）
    calculated_angle = calculate_angle(pointA, pointB, extended_point)
    print(f"预期角度: {angle}°, 实际计算角度: {calculated_angle:.1f}°")
    return extended_point

def calculate_bearing(pointA, pointB):
    lat1, lon1 = math.radians(pointA.latitude), math.radians(pointA.longitude)
    lat2, lon2 = math.radians(pointB.latitude), math.radians(pointB.longitude)
    diff_lon = lon2 - lon1
    
    x = math.sin(diff_lon) * math.cos(lat2)
    y = math.cos(lat1) * math.sin(lat2) - math.sin(lat1) * math.cos(lat2) * math.cos(diff_lon)
    
    initial_bearing = math.atan2(x, y)
    bearing = (math.degrees(initial_bearing) + 360) % 360
    return bearing

def calculate_angle(pointA, pointB, pointC):
    a = calculate_distance(pointB, pointC)
    b = calculate_distance(pointA, pointC)
    c = calculate_distance(pointA, pointB)
    
    # 余弦定理计算∠ABC
    cos_B = (a**2 + c**2 - b**2) / (2 * a * c)
    # 处理浮点精度导致的cos_B超出[-1,1]的情况
    cos_B = max(min(cos_B, 1.0), -1.0)
    angle = math.degrees(math.acos(cos_B))
    return angle

def calculate_distance(pointA, pointB):
    return geodesic((pointA.latitude, pointA.longitude), 
                   (pointB.latitude, pointB.longitude)).km

# 测试用例
pointA = Point(46.9540700,7.4474400)
pointB = Point(46.9560700,7.4494400)

# 测试90度
print("测试90度：")
c1 = extend_angle_line(pointA, pointB, 0.25, 1, 90)
# 测试105度
print("\n测试105度：")
c2 = extend_angle_line(pointA, pointB, 0.25, 1, 105)

1. 基准方向调整：将夹角的基准从A→B的方位角改为B→A的方位角，确保计算的是∠ABC的角度
2. 浮点精度处理：在calculate_angle中添加了对cos_B的范围限制，避免因地理计算的浮点误差导致math.acos报错
3. 函数命名优化：将extend_perpendicular_line改为extend_angle_line，更贴合功能实际（支持任意角度，不仅限于垂直）

测试后你会发现，实际计算出的角度会非常接近预期值（误差在0.1度以内，属于地理计算的正常精度范围）。

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