Free Great Circle Calculator

地球表面两点之间的最短路径并非地图上常见的那样呈直线，而是沿着一条特殊的弧线——大圆（Great Circle）。Great Circle Distance Calculator（大圆距离计算器）正是基于球面或椭球模型，快速计算这条最短路径的长度。无论是航空航线规划、航海导航，还是卫星通信覆盖分析，掌握大圆距离都具有重要价值。本文详细阐述大圆的几何定义、球面距离计算公式、Vincenty 椭球算法及其迭代过程，并通过具体示例展示两种模型的差异，帮助读者全面理解最短路径的计算原理。

大圆的定义与几何意义

从几何上看，大圆是球面与一个通过球心的平面相交所形成的封闭曲线。该平面将球体均分为两个半球，且大圆是球面上能画出的最大圆。对于球面上任意两个不同的点，通常可以确定唯一的大圆（唯一的例外是当两点恰好位于球心的对径位置，即互为对跖点时，存在无穷多个大圆通过它们）。大圆上较短的弧段即为连接两点的最短球面路径，称之为大圆距离。这条路径是球面上的测地线——即曲面上的“直线”。当我们使用 Shortest Path on Earth Calculator（地球最短路径计算器）时，实际就是在寻找这条测地线。

由于三维空间中的直线在球面上无法存在，大圆便成为两点间最自然的“最短路线”。这一特性在航空与航海中被广泛应用：飞机航线在地图上看似弯曲，实则遵循的是大圆路径。Sphere Distance Calculator（球面距离计算器）即采用球面三角公式计算该距离，而更精确的 Geodesic Distance Calculator（测地线距离计算器）则考虑地球的真实形状，使用椭球模型计算。

球面上的大圆距离计算

假设地球是一个半径为 $R$ 的标准球体。已知两点 $A$ 和 $B$ 的经纬度坐标分别为 $(\phi_A,\lambda_A)$ 与 $(\phi_B,\lambda_B)$（北纬为正、南纬为负；东经为正、西经为负），则大圆距离 $d$ 可由球面余弦公式直接给出：

$$d = R \cdot \arccos\!\big( \sin\phi_A \sin\phi_B + \cos\phi_A \cos\phi_B \cos(\Delta\lambda) \big)$$

其中 $\Delta\lambda = \lambda_A - \lambda_B$ 为经度差。该公式基于球面三角形的余弦定理，适用于任意球面上的两点。当两点接近对跖点时，反余弦函数的自变量接近 $-1$，数值稳定性可能下降，此时可改用半正矢公式替代。但在绝大多数情况下，上述公式已足够可靠且易于计算。使用 Sphere Distance Calculator 时，只需输入经纬度，工具即自动代入此公式得出结果。

椭球模型与 Vincenty 距离算法

真实地球并非完美球体，而是两极稍扁、赤道略鼓的旋转椭球体。为了更精确地描述地球形状，全球定位系统（GPS）采用 WGS84 参考椭球模型，其主要参数为：

• 长半轴 $a = 6\,378\,137.8\ \text{m}$
• 短半轴 $b \approx 6\,356\,752.314\ \text{m}$
• 扁率 $f = 1 / 298.257223563$

在此模型上计算大圆距离（更准确地说是测地线距离）需要采用 Vincenty 公式。Vincenty Distance Calculator（文森特距离计算器）以此为内核，通过迭代求解两点间的精确距离。其主要步骤如下：

1. 计算归化纬度（reduced latitude）：

$$U_A = \arctan\!\big((1-f)\tan\phi_A\big), \quad U_B = \arctan\!\big((1-f)\tan\phi_B\big)$$

2. 设初始值 $\lambda = L = \lambda_A - \lambda_B$（经度差）。

3. 在每次迭代中，计算以下中间变量：

$$\sin\sigma = \sqrt{(\cos U_B \sin\lambda)^2 + (\cos U_A \sin U_B - \sin U_A \cos U_B \cos\lambda)^2}$$ $$\cos\sigma = \sin U_A \sin U_B + \cos U_A \cos U_B \cos\lambda$$ $$\sigma = \operatorname{atan2}(\sin\sigma, \cos\sigma)$$ $$\sin\alpha = \frac{\cos U_A \cos U_B \sin\lambda}{\sin\sigma}$$ $$\cos^2\alpha = 1 - \sin^2\alpha$$

当 $\cos^2\alpha = 0$（即两点位于赤道）时，公式不适用，需另作处理。

$$\cos(2\sigma_m) = \cos\sigma - \frac{2 \sin U_A \sin U_B}{\cos^2\alpha}$$ $$C = \frac{f}{16} \cos^2\alpha \big(4 + f(4 - 3\cos^2\alpha)\big)$$ $$\lambda_{\text{new}} = L + (1 - C) f \sin\alpha \big( \sigma + C \sin\sigma [\cos(2\sigma_m) + C \cos\sigma (2\cos^2(2\sigma_m) - 1)] \big)$$

重复迭代直到 $|\lambda_{\text{new}} - \lambda|$ 小于预设阈值（通常 $10^{-12}$ 弧度），每次迭代后将 $\lambda$ 更新为 $\lambda_{\text{new}}$。

4. 收敛后，计算最终距离：

$$D = b \big( A\sigma - B \sin\sigma \cos(2\sigma_m) + C \sin\sigma \cos(2\sigma_m) \big)$$

其中系数 $A$、$B$ 等由中间变量进一步导出。完整的表达式在实际计算中由程序自动完成，用户无需手动迭代。

Vincenty 公式在绝大多数情况下能达到亚毫米级精度，但在处理对跖点或近乎对跖的点时收敛可能不畅。因此，此计算器内部设定最多 100 次迭代，若未收敛则输出当前精度值供参考。Shortest Distance Between Two Points on Earth Calculator（地球两点最短距离计算器）通常会同时提供球面与 Vincenty 两种结果，便于对比。

球面与椭球结果差异示例

为了直观展现两种模型的差异，考虑一个跨越半球的典型航线：纽约（40.71° N, 74.01° W）至悉尼（33.87° S, 151.21° E）。使用球面公式（取 $R = 6371\ \text{km}$）计算得到大圆距离为 $15988.75\ \text{km}$。若采用 WGS84 椭球模型并运行 Vincenty 算法，所得距离会略短约 $0.1\%$（具体缩短约 $15\ \text{km}$）。这表明，对于长距离航线，球面近似已经相当不错；但在需要极高精度的科研或精密导航中，椭球修正是必不可少的。

再以洛杉矶（33.94° N, 118.41° W）至伦敦（51.47° N, 0.45° W）的航线为例，球面与椭球计算结果的差异小于较长距离的 $0.25\%$。正是这种微小偏差累积，使得航线看似“绕道”格陵兰岛，实际上却是空间中最短的路径。

地图投影与视觉扭曲

大多数世界地图采用墨卡托投影，其特点是保持角度正确，但面积和形状在高纬度地区被严重放大。因此，大圆路径在平面地图上呈现显著的弯曲，而地图上的直线（等角航线）反而更长。Great Circle Distance Calculator 可以帮助我们直观地理解这一现象：当你输入两点的坐标并查看结果时，背后的几何原理正是测地线。使用该工具，只需选择适合的模型（球面或椭球），即可立即获得两点间的真正最短距离。

无论是航空公司的飞行计划、跨洋航线的设计，还是业余无线电的通联距离估算，大圆距离计算器都提供了准确、便捷的解决方案。掌握其计算公式与适用条件，能让用户在最短路径规划中做出更明智的决策。