如何计算点到线段中点的连线与该线段的夹角？

我来帮你排查代码里的问题啦～首先你的需求是计算p1到p2-p3线段中点的连线，和p2-p3线段本身的夹角，但代码里的向量选得不对，这是结果出错的核心原因！

先梳理正确的向量逻辑 #

要计算两条线的夹角，我们需要取两条线上的同向/反向向量：

• 第一条线：p1 ↔ 中点M，我们可以取向量 M → p1（坐标差为 p1 - M）
• 第二条线：p2 ↔ p3，我们可以取向量 p2 → p3（坐标差为 p3 - p2）

你的代码里错误地把第二条向量写成了 p3 → p1，这完全不是p2-p3线段的方向，结果自然不对啦！

修正后的代码 #

先补上math模块的导入，然后调整向量的定义：

import math

def calculateAngle(point1, point2, point3):
    lon1, lat1 = point1
    lon2, lat2 = point2
    lon3, lat3 = point3
    
    # 计算p2-p3的中点M
    latCenter = (lat2 + lat3) / 2
    lonCenter = (lon2 + lon3) / 2
    
    # 构造正确的两个向量
    # 向量1：中点M指向p1
    xV1 = lat1 - latCenter
    yV1 = lon1 - lonCenter
    # 向量2：p2指向p3（代表p2-p3线段的方向）
    xV2 = lat3 - lat2
    yV2 = lon3 - lon2
    
    # 计算点积和模长
    dotProduct = xV1 * xV2 + yV1 * yV2
    magnitudeV1 = math.sqrt(xV1 ** 2 + yV1 ** 2)
    magnitudeV2 = math.sqrt(xV2 ** 2 + yV2 ** 2)

    if magnitudeV1 == 0 or magnitudeV2 == 0:
        return 0.0, 0.0
    
    # 计算夹角（处理数值溢出问题）
    cosTheta = dotProduct / (magnitudeV1 * magnitudeV2)
    cosTheta = max(min(cosTheta, 1), -1)
    theta = math.acos(cosTheta)
    angleInDegrees = math.degrees(theta)
    
    return angleInDegrees, theta

额外注意点 #

• 经纬度的平面近似问题：你直接把经纬度当作平面坐标计算，如果是小范围的地理区域（比如城市内部），这个近似可以接受；但如果是跨越大范围（比如跨国家），平面坐标的误差会很大，这时候需要用球面坐标系的夹角计算方法。
• 向量方向不影响夹角：不管你取M→p1还是p1→M，p2→p3还是p3→p2，最终计算的夹角结果都是一样的，因为余弦值是偶函数。

附上你提供的示意图：

备注：内容来源于stack exchange，提问作者mikezang