应使用安全截断函数处理acos输入：const safeAcos = x => Math.acos(Math.max(-1, Math.min(1, x)))，避免浮点误差导致NaN；推荐鲁棒球面余弦公式计算cosσ后调用safeAcos；对cosσ>0.999999的情况启用平面近似优化。
直接用 Math.acos() 计算大圆距离时，若两点几乎重合或几乎对跖（经度差≈180°且纬度符号相反），输入值可能略超 [-1, 1] 范围（如 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 或球面余弦公式都可能在极近或极远时放大浮点误差。推荐用以下形式（基于余弦定理，但结构更鲁棒）：
- 设
φ1, φ2为两点纬度（弧度），λ1, λ2为经度（弧度） - 计算：
const cosσ = Math.sin(φ1)*Math.sin(φ2) + Math.cos(φ1)*Math.cos(φ2)*Math.cos(λ2−λ1); - 立即应用：
const σ = safeAcos(cosσ); - 再乘地球平均半径（如 6371e3 米）得距离
该形式比先算 Δλ 再调用 cos(Δλ) 略少一次三角运算，且 cos(λ2−λ1) 在 Δλ 接近 π 时本身精度尚可（cos 在 ±π 处导数为 0，误差传播弱）。
对极接近点（cosσ ≈ 1）做短路优化（可选）
当两点非常接近（如 cosσ > 0.999999），acos 在 1 附近导数趋于无穷（d(acos x)/dx = −1/√(1−x²)），微小输入误差会导致角度误差被剧烈放大。此时可改用平面近似或 Haversine 的小角度展开：
- 判断：
if (cosσ > 0.999999) { const Δφ = φ2 − φ1; const Δλ = λ2 − λ1; const a = Δφ*Δφ + Math.cos(φ1)*Math.cos(φ2)*Δλ*Δλ; return 6371e3 * Math.sqrt(a); } - 这是球面正交投影的一阶近似（单位：弧度），对百米级距离误差 acos
- 阈值
0.999999对应夹角约 0.001 rad ≈ 63 km，可根据精度需求调整
验证边界行为的简单测试用例
写几个关键测试确保鲁棒性：
-
safeAcos(1.0 + 1e-15)→ 返回0（不是NaN） -
safeAcos(-1.0 - 1e-15)→ 返回Math.PI - 同一点：
φ1=φ2=0.5, λ1=λ2=0.3→cosσ ≈ 1.0→ 距离 ≈ 0 - 对跖点：
φ1=0.1, φ2=-0.1, λ1=0, λ2=Math.PI→cosσ ≈ -1.0→ 距离 ≈ 20015 km
不复杂但容易忽略
