如何通过调整 Math.acos() 输入值精度来优化地理大圆距离计算？

本文共计834个文字，预计阅读时间需要4分钟。

直接使用 `Math.acos()` 计算大圆距离时，若两点几乎重合或几乎对齐（经度差小于180°且符号相反），输入值可能超出范围（如 1.0000000000000002 或 -1.0000000000000002），导致结果为 `NaN`。这不是公式错误，而是由于浮点数误差引起的。关键不是避免使用 `acos`，而是确保安全截断输入值。

用 clamp 技术预处理 acos 输入值

Math.acos(x) 只接受 x ∈ [-1, 1]。当球面余弦公式结果因浮点误差超出该区间时，应主动“拉回”到合法范围，而非抛异常或跳过：

• 定义安全函数：const safeAcos = x => Math.acos(Math.max(-1, Math.min(1, x)));
• 它等价于“硬限幅”：小于 -1 → 设为 -1；大于 1 → 设为 1；中间值不变
• 这比 isNaN(x) ? 0 : Math.acos(x) 更合理——因为 x=1.0000000000000002 的物理含义就是“夹角为 0”，对应距离为 0，而 acos(1)=0

使用更稳定的球面余弦公式变体

标准 Haversine 或球面余弦公式都可能在极近或极远时放大浮点误差。