基于C#的GPS两点方位角计算实现咨询

我帮你把这个方位角计算的代码补全并修正输出范围的问题，直接看完整实现：

public static double Bearing(IPointGps pt1, IPointGps pt2)
{
    // 重要提示：如果你的IPointGps存储的是角度值（如十进制度数），必须先转弧度
    // 取消下面注释启用转换：
    // double lat1 = Math.PI * pt1.Latitude / 180;
    // double lat2 = Math.PI * pt2.Latitude / 180;
    // double lon1 = Math.PI * pt1.Longitude / 180;
    // double lon2 = Math.PI * pt2.Longitude / 180;

    double lat1 = pt1.Latitude;
    double lat2 = pt2.Latitude;
    double deltaLon = pt2.Longitude - pt1.Longitude;

    double x = Math.Cos(lat1) * Math.Sin(lat2) - Math.Sin(lat1) * Math.Cos(lat2) * Math.Cos(deltaLon);
    double y = Math.Sin(deltaLon) * Math.Cos(lat2);

    // 初始计算的方位角范围是 -π 到 π
    double bearingRadians = Math.Atan2(y, x);

    // 转换为 0 ≤ 结果 < 2π 的目标范围
    bearingRadians = (bearingRadians + 2 * Math.PI) % (2 * Math.PI);

    return bearingRadians;
}

关键细节说明 #

• 弧度制必须注意：.NET的Math类三角函数参数都是基于弧度的！如果你的IPointGps接口存储的是角度（比如常用的十进制GPS坐标），一定要先执行注释里的弧度转换，否则计算结果会完全错误。
• 解决Atan2负值问题：Math.Atan2(y, x)的返回范围是-π到π，通过(值 + 2π) % 2π的操作，既能把负数转为正区间的角度，也能保证正数不会超出2π（比如刚好等于2π时，取模后会回到0）。
• 公式有效性：你选用的是球面几何中计算两点间初始方位角的标准公式，逻辑是正确的，只要单位匹配就能得到准确结果。

如果需要返回0°到360°的角度值，只需要在最后把弧度转成角度：

return bearingRadians * 180 / Math.PI;

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