C++如何计算两个地理坐标点（经纬度）的方位

<p>atan2能根据dx、dy符号准确判断象限，返回[-π, π]弧度值，而atan仅返回[-π/2, π/2]且无法区分象限；地理方位角需转换为[0, 360)度：fmod(atan2(dx, dy) * 180.0 / M_PI + 360.0, 360.0)，输入经纬度须先转弧度，极区需校验纬度范围或改用GeographicLib。</p>

用 atan2 实现方位角计算，别用 atan

直接用 atan 会丢象限信息，导致正北和正南都返回同一个值。必须用 atan2(dy, dx)，它能根据坐标差的符号自动判断 quadrant（象限），结果范围是 [-π, π]。

实际地理中，方位角定义为：从正北方向顺时针旋转到目标点的角度，范围是 [0, 360) 度。所以得把 atan2 的输出先转成角度，再做偏移和模运算：

• dx = sin(Δλ) * cos(φ₂)（东西向分量，已考虑纬度缩放）
• dy = cos(φ₁) * sin(φ₂) - sin(φ₁) * cos(φ₂) * cos(Δλ)（南北向分量）
• 最终方位角：fmod(atan2(dx, dy) * 180.0 / M_PI + 360.0, 360.0)

注意经纬度单位和顺序：输入必须是弧度，不是度

C++ 标准数学函数（sin、cos、atan2）只认弧度。如果原始数据是度（比如 lat = 39.9042, lon = 116.4074），必须先转：

• φ₁ = lat1 * M_PI / 180.0
• φ₂ = lat2 * M_PI / 180.0
• Δλ = (lon2 - lon1) * M_PI / 180.0

漏掉这步，结果完全错误，且不容易排查——因为数值看起来“像那么回事”，只是方向全反或偏 90°。

立即学习“C++免费学习笔记（深入）”；

当两点纬度接近极点时，dy 可能趋近于 0

公式里 dy 在极区（比如两个点都在北纬 85° 附近）可能非常小，导致 atan2(dx, dy) 对浮点误差敏感，方位角抖动明显。

这不是代码写错了，而是球面模型在极点退化所致。应对方式只有两种：

• 业务上明确限制输入范围（例如不支持纬度 > ±85°）
• 改用更鲁棒的库，比如 GeographicLib 的 Direct 或 Inverse 函数，它内部做了极点处理

自己硬补特判容易引入新 bug，不如早做输入校验。

别忽略地球椭球模型带来的偏差

上面用的是球面近似（WGS84 简化版），对几百公里内的点误差通常

如果项目要求亚度级精度：

• 不要手写公式，用 GeographicLib（C++ 原生支持）
• 调用其 Geodesic::Inverse，返回的 azi1 就是标准正北顺时针方位角
• 注意链接时加 -lGeographicLib，头文件是 #include <geographiclib></geographiclib>

球面公式够用就别折腾，但一旦发现实测方位和计算值系统性偏移 0.3° 以上，大概率是模型问题，不是算错。