趣味大地测量学/主题 01 · 地球椭球理论/第 1 集

第一次算出一个数

地图上的距离,到底是怎么算出来的

武汉 30.5928°N 114.3055°E北京 39.9042°N 116.4074°E

EP 013:53直观入门
第一次算出一个数 · 封面
在 B 站看在 YouTube 看点击封面开始播放

本集速览

  • 真实贴地曲线最符合物理直觉,却没有解析公式可算,且实地逐点丈量无法推广到全球任意两点。
  • 问题发生重构:不再执着于「哪条曲线最真实」,而是转向「哪条曲线最可用」
  • 模型降级:主动将起伏地球简化为半径 6371 km 的光滑圆球,舍弃局部地形以换取可算性。
  • 球面两点间最短路径是大圆弧;弧长由半径乘以圆心角决定,圆心角通过空间向量点积求得。
  • 代入武汉与北京实算得出 1052.7 公里——与现代地图软件测距完全吻合,揭示其背后的圆球假设。

哪条线最可用▶ 0:00

当我们在地图软件中搜索两地之间的距离——例如从武汉到北京——屏幕上会给出一个清晰确定的数字。但在一个真实起伏的地球表面上,两点之间的「距离」究竟该如何定义,又是怎么算出来的?

在三维空间中连接地球表面的两点,几何上有三条最自然的候选路径:

  1. 穿地心的直线(弦):在三维欧几里得空间里,连接两点的直线段在数学上绝对最短。但这条线必须直接穿透地球内部的地壳与地幔;除非人类能打通一条贯穿地底的真空隧道,否则它无法作为地面通行与测量的基准。
  2. 平滑的几何圆弧:在几何面上拉出一条规则的平滑弧线。它虽然不穿过地心,但悬在空中或埋在地下,与起伏的真实地表存在偏差。
  3. 贴合地表的起伏曲线:严格贴着地表山川起伏走出的实际轨迹。从直觉上看,这无疑是最真实的一条路线——任何人在地面上移动,实际走过的都是这条起伏的路径。

但「真实」与「可用」完全是两回事。

贴地曲线虽然最符合现实,却面临一个致命的数学困境:真实地球表面千沟万壑,山峰、峡谷与平原错落交织。这种复杂的地形特征是局部的、不规则的,物理上没有任何公式能够解析表达这整条起伏曲线的数学方程,更无法对其直接积分求长。

贴地曲线沿着起伏的地表,右上角放大出剖面的锯齿
图版 1 · 视频 0:32贴着地面的那条线(橙)确实最真实,但它每一段都在起伏——放大之后是一条无法写成公式的锯齿。

如果无法用公式计算,能不能退而求其次:跳过数学推导,直接派人带着测绳或仪器在地面上实地量一遍,把数值记录下来备查?

对于一条特定的铁路线或公路,工程勘测确实会实地丈量。但在全球尺度上,地球表面存在无穷多个点对组合。如果任意两地之间的距离都必须依赖人力先在崇山峻岭间实地量一次,不仅成本极其高昂,在逻辑上也根本无法推广。

因此,大地测量学在出发的一刻,就必须把问题做一次根本性的重构:面对复杂的物理现实,我们不该执着于寻找「哪条曲线最真实」,而应该寻找「哪条曲线最可用」。

换一个有公式的表面▶ 0:57

科学与工程中的「可用」,通常必须同时满足两个硬性条件:

  1. 能算(可计算性):对象的几何形状必须拥有简洁的数学解析式,使得任意两点间的距离可以通过几何推导得出,而不是每次都依赖实地勘测;
  2. 足够接近(有效近似):简化后的模型不能脱离物理现实太远,计算出来的数值必须在误差允许范围内具备实际的工程与参考价值。

如果用这两个标准去重新审视前面的三条候选线,平滑的几何圆弧之所以脱颖而出,并不是因为它完美还原了地表的起伏,而是因为它背靠一个拥有完美数学公式的几何表面——球体

这一步在方法论上完成了一个关键动作:模型降级

面对复杂到无法建立方程的真实地表,大地测量学没有选择死磕每一个局部的山头与峡谷,而是主动退后一步——用一个高度简化的、数学上完全光滑的规则几何体来替代真实的物理地球。我们丢弃了千变万化的局部地形,换取了全局的、确定性的可计算性。

光滑圆球上标出两点
图版 2 · 视频 1:27起伏没有了。代价是地形被整个丢掉,收获是任意两点之间的距离从此可以计算。

这是一笔明码标价的科学交易。在整门大地测量学的演进史中,我们还会反复看到这个动作:真实地表 → 圆球(本集)→ 旋转椭球(后续章节)→ 投影平面。每一次模型的切换,都是在「计算成本」与「描述精度」之间重新寻找平衡。

早在两千多年前的古希腊时期,学者们就已经把脚下的世界抽象为一个光滑的圆球。由此,我们建立起本课程的第一个基准地球模型:一个地球平均半径取为 6371 公里的光滑正球体。

在这个理想的球面上,两点之间的最短路径和距离该如何计算?问题就此进入了纯粹而优雅的球面几何世界。

球上的最短 · 大圆弧▶ 1:29

在一只光滑的圆球表面上,连接任意两点的路径有无穷多条。哪一条才是真正的最短路径?

我们可以借助一个物理直觉:想象球面上没有任何摩擦阻力,在两点之间紧绷一根弹性橡皮筋。在张力的作用下,橡皮筋会自然贴紧球面滑移,最终停留在一条长度最短的几何轨迹上。

在几何学中,这条路径被称为大圆弧

大圆与小圆

什么是大圆?

在三维空间中,球面上不重合且非对跖(不关于球心对称)的任意两点,与球心三者可以唯一确定一个平面。当这个平面切过球体时,截面与球面的交线就是一个圆。

  • 大圆(Great Circle):由经过球心的平面截球面所得的圆。由于截面恰好平分球体,大圆的半径等于球体本身的半径 ,是球面上所能画出的最大圆周;大圆上的任意一段圆弧,就称为大圆弧
  • 小圆(Small Circle):由不经过球心的平面截球面所得的圆(例如赤道以外的所有纬线圈),其半径

在球面上,半径越大的圆弯曲程度越平缓。因此,连接两点的大圆弧在所有贴附球面的曲线中弯曲度最小,它正是球面上两点之间的测地线(最短路径)。

这一步只用中学几何就能证到底,橡皮筋的直觉背后是一条干净的推理:深入大圆弧为什么是最短的

两点与球心确定的平面截出大圆,弧与两条半径构成扇形
图版 3 · 视频 1:52两点 + 球心 → 一个平面 → 截出大圆。要算的量在这张图上已经全了:两条半径和它们之间的夹角。

一个反直觉的现象:沿纬线走并非最短

大圆弧有一件打破直觉的事:如果两座城市恰好位于同一条纬线上(例如同处北纬 40° 附近的北京与马德里),沿着纬线正东或正西方向走,并不是最短路线

在常见的平面长方形地图(如墨卡托投影)上,同一条纬线表现为一条水平横线,给人一种「直着走最近」的视觉错觉。但在真实球面上,除赤道以外的所有纬线圈都是小圆。沿着小圆行走,相当于一直在绕着一个较小的半径侧向转弯,路线被无形拉长了。

真正最短的大圆弧,实际上会向高纬度地区弯曲拱起,甚至切入极地边缘。在现实中的民航远程航线上(例如从北京飞往纽约或欧洲),飞机在平面航线图上往往呈现为一条深拱向北极的高纬度弧线——那条看似绕了大远路的弯曲航迹,正是为了追求航程最短、油耗最省的大圆航线。

到底绕远了多少,把两条路都算一遍就知道:硬核沿纬线飞到底绕远了多少

地图上纬线路径被否掉,弯向高纬的大圆弧才是最短
图版 4 · 视频 2:15同纬度两地,沿纬线走(直线)反而更长;弯向高纬的那条才是大圆弧。

化归为平面几何

最短的路径既然已经锁定为大圆弧,它的长度该怎么计算?

由于两点和大圆弧完全包含在过球心的那个截面内,三维球面上的弧长计算,在此刻直接化归为最基础的平面扇形几何

在截面圆盘中,两点与球心构成了半径为 的圆弧,其长度 直接遵循中学几何公式:

其中 是已知的球半径(), 是从球心看两点所张开的圆心角(以弧度为单位)。

半径已经确定,弧长公式也极其简单。求解两地距离的全部问题,最终收敛为一个纯粹的几何任务:如何根据已知的经纬度坐标,求出两点与球心之间的圆心角

实算与回身 · 1052.7▶ 2:31

有了「弧长 = 半径 × 圆心角」的基础,计算两点距离的核心就在于求解夹角

虽然球面三角学中有现成的球面余弦定理,但在现代测量学与计算机计算中,更直观、通用且便于推广的方法是空间向量法。两条路其实通向同一个式子:深入点积法与球面余弦定理是同一个式子

第一步:建立空间直角坐标系

我们将光滑的圆球地球放入三维空间直角坐标系 中:

  • 坐标原点 定在地球球心;
  • 轴沿着地轴指向北极;
  • 轴穿过赤道与本初子午线(0° 经线)的交点;
  • 轴在赤道平面内与 轴垂直指向东经 90°,构成标准的右手坐标系。

球面上任意一点的位置由大地纬度 与大地经度 确定。根据三维空间中的三角投影关系,该点的空间直角坐标 为:

其中 为纬度, 为经度。从几何上看, 是点在赤道平面上的投影半径,再沿经度分解为 则是沿 轴的垂向高度。这套转换关系在主题 04 讲坐标系统时还会再次登场。

这三行不必硬记,拆成两个平面直角三角形就能一眼看出来:深入从经纬度到直角坐标的那三行

如果将球半径归一化为单位长度(),由球心指向该点的单位向量,即为该点的方向向量

第二步:向量点积求圆心角

当武汉与北京两地的方向向量 确定之后,两向量从球心出发所张开的夹角 ,可以通过高中数学的向量点积公式直接求解:

由于 均为单位向量(模长为 1),分母为 1,点积直接等于两向量夹角的余弦值:

球心出发的两个方向向量与坐标轴,右侧三张公式卡
图版 5 · 视频 2:58经纬度 → 方向向量 → 点积 → 夹角。三步推导全程只用到高中的向量与三角。

通过反余弦函数 ,即可直接解出以弧度为单位的圆心角

第三步:代入实算

现在将武汉与北京的经纬度坐标代入上述算法:

  • 武汉
  • 北京

计算步骤如下:

  1. 分别算出两地的方向向量
  2. 计算点积
  3. 求反余弦得到圆心角 弧度(换算为角度约为 9.47°);
  4. 乘以地球平均半径 ,得出两地大圆弧距离:

想自己跑一遍的,这条链十行代码就能复现:硬核十行代码复现 1052.7

左侧三维坐标系中的两个向量,右侧圆心角与 1052.7 公里的结果
图版 6 · 视频 3:14第一个数落地。整门课的基准数字,这是第一个。

回身:地图软件背后的球假设

打开地图软件搜索武汉到北京的直线距离,屏幕上给出的正是 1052.7 公里

这个每天被无数人查询、从未被怀疑过的数字,本质上就是我们刚才算出的那段大圆弧。换句话说,地图软件在计算两点距离时,底层默认使用的正是半径 6371 公里的正球体模型

至此,第一笔交易达成了:我们主动舍弃局部地形,换来了全局可算的确定性,在复杂的物理现实中拿到了第一个可用、可信的基准数字。

说明
武汉30.5928°N, 114.3055°E全课统一取点
北京39.9042°N, 116.4074°E全课统一取点
球半径 R6371 km地球平均半径
圆心角 θ0.16524 rad ≈ 9.47°由点积求出
大圆弧长 s1052.7 km本集算出的那个数
穿地心的弦1051.5 km比大圆弧只短 1.2 km

拓展

深入
硬核

动手试一试

交互
  • 大圆弧探索器

下一个问题

然而,这一整套干净的计算,自始至终都押在一个未经检验的前提上:地球是一只半径 6371 公里的正球体。

如果追问下去,两个更根本的问号立刻浮出水面:人类站在平坦的地面上,凭什么断定大地是球形的?而 6371 公里这个决定了全局尺度的半径,历史上又是谁第一次量出来的?

那时候没有卫星,希腊人的活动范围也只有地中海附近一小块——他们怎么知道地球是圆的,甚至还知道地球有多大?

趣味大地测量学GEODESY本页可独立阅读,与视频内容一致出品方标志

鄂ICP备2025098393号-2