<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) 对浮点误差敏感，方位角抖动明显。
C++ Code Review Master
下载
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
这不是代码写错了，而是球面模型在极点退化所致。应对方式只有两种：
- 业务上明确限制输入范围（例如不支持纬度 > ±85°）
- 改用更鲁棒的库，比如
GeographicLib的Direct或Inverse函数，它内部做了极点处理
自己硬补特判容易引入新 bug，不如早做输入校验。
别忽略地球椭球模型带来的偏差
上面用的是球面近似（WGS84 简化版），对几百公里内的点误差通常
如果项目要求亚度级精度：
- 不要手写公式，用
GeographicLib（C++ 原生支持） - 调用其
Geodesic::Inverse，返回的azi1就是标准正北顺时针方位角 - 注意链接时加
-lGeographicLib，头文件是#include <geographiclib></geographiclib>
球面公式够用就别折腾，但一旦发现实测方位和计算值系统性偏移 0.3° 以上，大概率是模型问题，不是算错。
