JED 2026年1月Dave Adamy空间几何专栏深读
站得更高,通常能看到更远的地平线。但当观察者换成卫星,“更远”至少有两种意思:卫星到地平线点的直线距离,以及从卫星正下方沿地表走到地平线的距离。它们来自同一个三角形,却不能在计算和比较时相互替换。
Dave Adamy在JED 2026年1月的《Space EW (Part 8) – Horizon Calculation》中,用球形地球模型讨论这一问题。[1] 本文围绕基本几何展开,重点解释距离、角度、单位和舍入,说明原文的计算能够回答什么,以及为什么几何可见不能直接等同于传感器已经获得有效观测。
原刊延续空间几何系列,把问题从“两个位置之间有多远”推进到“它们之间的直线是否被地球遮住”。这是一个可见性问题。在本文采用的理想模型中,地球是半径固定的球,传播按直线处理,地面点位于球面上;地形起伏和实际传播环境暂不加入。
这些条件不是为了把真实地球说成绝对光滑,而是让基本关系能够清楚地推导。如果一条视线穿过球体,球面上的相应点就不属于此时的几何可见范围。边界情况则是视线刚好与球面相切,接触点便是这条剖面上的地平线点。
这种“看见”尚未涉及设备是否指向那里、是否具有足够灵敏度,或者是否得到可用图像。NASA关于宇航员摄影的说明也区分正下方的天底位置与倾斜拍摄的实际图像中心。[2] 位置几何提供一个观察范围的背景,并不代替具体仪器的视场说明。

图1 卫星视界的球面剖面。编辑依据原刊几何关系重绘;红线为卫星至地平线点的直线距离,金色弧线为星下点至地平线的地表弧长,二者起点不同。
图1中的O是地心,S是卫星位置。连接O和S,这条线与近侧球面的交点P就是本模型中的星下点。高度h从P量到S,地心距r从O量到S,地球半径R则从O量到球面。这三个长度的关系非常直接:地心距等于地球半径加上高度。
再从卫星向球面作一条切线,切点记为T。S到T的线段长度是L,表示卫星与地平线点之间的直线距离;P到T沿球面的短弧长度是s,表示沿地表的距离。两者不只是“直线与曲线”的区别,还具有不同的起点,不能放在同一栏里而省略定义。
在三维球面上,不同方向的地平线点构成可见球冠的边界,图中只截取其中一个包含地心的平面。于是,原本看起来复杂的空间关系,可以先化为平面三角形和圆弧。原刊也提醒示意图不按实际卫星高度比例绘制,读者不应从印刷图上的尺长反推数据。
本文改用R、r、L、s等符号,是为了避免原刊中同一个小写字母在不同位置承担不同含义。原文三角形里的a表示地球半径,表格里的a却表示轨道半长轴。两处各有定义,但若把它们直接拼成一段计算,就很容易代入错误的长度。
圆的半径与切线在切点处垂直。因此,O、T、S构成的三角形在T处为直角,斜边是卫星到地心的距离r,两个直角边分别为地球半径R和直线距离L。原刊用勾股定理得到到地平线的直线距离,关键不在公式复杂,而在于哪一条边是斜边。
从这个关系立即可以检查一个结果是否基本合理:L应小于r,因为它是直角边;L又大于高度h,因为从卫星斜向地平线的距离,比直接向下到星下点更长。这些判断只依赖本模型的几何关系,可以在输入数值以前使用。
公式还可展开成高度形式。展开后含有高度的平方项,不能在轨道尺度下未经检查就省略。只有在高度远小于地球半径时,忽略这一项才可能构成有用近似。原刊的示例高度为4,189千米,已经不是可以随意按很低观察点处理的尺度。
所有长度必须使用同一单位。若半径以千米表示,高度也应以千米代入,结果才以千米给出。把米与千米混在平方项中,程序仍可能返回一个数字,但那个数字失去了正确的物理含义;公式正确并不能防止单位输入错误。

R为球形地球半径,h为卫星距球面的高度,r为地心距,L为卫星至地平线点的直线距离。四者均为长度,本文统一用千米;要求R大于零、h不小于零,且采用球面相切与直线传播模型。
直线距离算出以后,还需要确定球面上的弧长。图1把地心O处的夹角记为θ,它位于指向卫星的方向与指向地平线点的半径之间。在同一个直角三角形中,邻边是R,斜边是r,因此可以由余弦关系得到θ。
原刊采用正弦定理和反正弦求地心角。对于此处的直角三角形,两种路线在保留足够精度时相容。本文选择余弦写法,是因为已知R和r时可直接计算,不必先把L舍入成整数再带入下一步;这属于等价整理,不是对原文基本几何的否定。
圆弧长度等于半径乘以对应的弧度数。这里最常见的错误,是把以度为单位的角度直接乘上半径。度和弧度描述同一个角,却使用不同的数值标度;如果计算器输出的是度,就必须先做转换,或者按整圈周长的比例计算。
原文使用大圆周长40,030千米,并乘以地心角占360度的比例。以R等于6,371千米计算的周长约为40,030.17千米,两者非常接近。它们在本例只造成很小的结果差别,但对复算而言,仍应说明到底沿用哪一种常量口径。

θ为地心角,采用弧度;s为星下点至地平线点沿大圆方向的球面短弧长,单位与R相同。若使用度数,应先乘以π/180转换为弧度,不能将度数直接代入s等于R乘θ。
同一幅图还包含卫星位置S处的角。它是从指向地心的天底方向转到地平线视线所夹的角,本文记为α,称为离天底角。因为T处为直角,α与地心角θ互余。在本例中,一个约为37.11度,另一个约为52.89度,二者不是同一个角。
原刊将相关角度称作elevation angle,并在后文补上“from nadir”。如果只翻译成“仰角”,就容易让读者误以为它是在地面相对水平面测得的角。按图2的顶点和边线来看,该角位于卫星,而不是星下点;因此,更明确的说法是说明顶点在何处、从什么方向量起。
地面观测者通常谈论的仰角,以当地水平面为参考。ESA Navipedia的局部东、北、天坐标说明,将这一角度与观察位置的局部坐标联系起来。[3] 在本文的理想地平线点T处,卫星视线就在切平面内,所以地面仰角为零;这与卫星处的37.11度没有冲突。
检查角度时,可以依次问三个问题:角的顶点在哪,零方向是哪条线,数值用度还是弧度。只要其中一个条件不同,相同的“仰角”字样就可能代表不同量。本文保留原刊数值关系,同时把命名中的参考方向明确写出,避免让术语代替几何定义。

α为卫星处从天底方向量起的离天底角,θ为地心角;此式二者均以弧度表示,π/2为直角。若二者改用度表示,关系是两角之和为90度。该式不是任意地面点的仰角公式。
原刊给出一个3小时周期的圆轨道例子,采用高度4,189千米、地心距10,560千米与地球半径6,371千米。本文将这些公开的整数作为复算输入,不另行根据轨道周期反推参数。高度加半径恰好等于地心距,说明这三个输入之间是一致的。
按前面的相切三角形,直线距离约为8,421.636千米,地心角约为52.892356度。以半径乘弧度计算,地表弧长约为5,881.362千米;若沿用原刊40,030千米周长,则约为5,881.336千米。整数级表达均约为5,881千米,与原刊正文相符。
需要留意的是,原文将中间比值写成0.798,紧接着给出反正弦为52.89度。若真的把已经舍入的0.798重新输入计算器,结果约为52.94度。52.89度与保留较多中间精度的计算相符,不能把印刷上缩短的小数当成后续计算实际使用的完整输入。
下表把原刊展示值与编辑复算并列,目的在于解释计算口径,而不是增加观测精度。输入本来就是整数千米,本文列出更多小数只是为了显示运算与舍入如何衔接,并不表示卫星位置或地球模型具有相应的小数级准确度。
表1 原刊示例的数值与同一输入下的复算
| 项目 | 原刊展示 | 编辑复算 |
|---|---|---|
| 卫星至地平线直线距离 | 8,422 km | 约8,421.636 km |
| 地心角 | 52.89° | 约52.892356° |
| 星下点至地平线弧长 | 正文5,881 km | 约5,881.362 km |
| 同一例子的表格弧长 | 5,880 km | 与上述统一口径略有差别 |
原刊表1把轨道周期、高度、半长轴、直线距离和地表距离放在一起。周期单位是分钟,其余距离采用千米。它分左右两组排列,右半部分是较长周期的后续行,不是另一套物理模型。若从抽取文字复制数据,应先与实际表格版面核对列的对应关系。
在180分钟这一行,原表的高度和地心距与正文例子相同,但地表距离写为5,880千米,正文则为5,881千米。按本文统一使用的整数输入复算,结果约为5,881.36千米。这个差异应当说明,而不是悄悄改掉后再声称原文始终一致。
表中其他行也存在少量千米级差别。仅凭印刷整数,不能确定作者生成整张表时是否使用了隐藏精度、不同常量或其他取整过程。因此,本文不把差异全部归为同一种原因,也不把原表重新包装成高精度计算结果。使用表格前,仍应确认它提供的是哪一种近似。
原刊说明,表格中的半长轴列可以与不同轨道形状联系,但高度列和视界距离按圆轨道讨论。圆轨道的地心距保持不变,因此在球形模型中,相应的地平线距离也保持不变。注意,这不意味着地表上始终是同一片区域可见,卫星位置仍在随时间变化。
椭圆轨道则不能把半长轴直接当作任意时刻的地心距。本文的相切三角形仍然是一种瞬时几何关系,但输入应是当时的地心距或高度,而不是未经区分的平均尺度。原文关于“适用于任意轨道形状”的半长轴说明,不能扩大为整张高度与视界表都通用。
另一个容易混淆的量是轨道周期。它描述绕行一周的时间,不直接等于某个地面点处于视野内的时间。原刊把过顶飞越的可见时长留到后续讨论,也表明本篇回答的是空间边界,而不是完整的时间过程。
因此,用这篇文章学习几何时,可以先把卫星位置固定在一个瞬间,再判断可见边界。若进一步讨论连续覆盖,需要另行描述位置如何变化以及观察条件如何限制。这些问题不能由单个最大地平线距离自动回答,本文也不据此生成轨道安排或具体任务窗口。
原刊沿用上期例子,分别比较地面点与地平线的直线距离,以及对应的地表距离。它的思路是用同一类量作判断:从卫星出发的直线距离,对照卫星到地平线的直线距离;从星下点出发的弧长,对照星下点到地平线的弧长。
在本文限定的同一球面和同一卫星位置下,这两类比较可以相互核对。但如果把某个地面点的海拔、其他观察位置或不同传播模型混入,原有对应关系就需要重新检查。公式的适用对象是关系成立的组成部分,不能因为两个结果都以千米表示就直接比较。
尤其不能把直线距离较大的数字当成较大的“地表覆盖半径”。前者从卫星开始,后者从星下点开始。图1已经显示二者的起点差异;表1中8,422与5,881两个相近但明显不同的数值,恰好说明同一问题可以有不同而同时正确的距离答案。
还可以用极端情况检查公式。当高度趋近于零,地心距趋近于地球半径,直线距离和地表弧长都趋近于零。这符合光滑球面上一个无高度观察点的几何极限,并不意味着站在真实地面上的人永远没有视野;现实中的眼睛高度、地形与环境尚未放入这个理想化条件。
当高度不断增大,地心角趋近于直角,星下点至地平线的弧长趋近大圆周长的四分之一。直线距离则仍可继续增大。这个差别说明“离地更远”不等于“沿地表能看得无限远”,瞬时球面可见性的几何上限与传播路径长度不是同一件事。
真正需要带走的不是一个脱离上下文的最远距离,而是一套清楚的提问方式:从哪里量到哪里,沿什么路径,用什么角度,以及结果能够证明到哪一步。几何计算由此成为理解观察条件的基础,而不是对实际探测、识别或持续覆盖能力的笼统保证。
[1] Dave Adamy. Space EW (Part 8) – Horizon Calculation. JED,2026年1月,第40—41页。正文与表1数值分别保留,并按其给定的R为6,371 km、r为10,560 km复算。
[2] NASA Johnson Space Center,Gateway to Astronaut Photography of Earth,FAQ中的nadir、oblique与图像中心说明。
[3] J. Sanz Subirana等,Transformations between ECEF and ENU coordinates,ESA Navipedia,2011。用于局部坐标与仰角定义。
[4] 示意图、等价公式与数值复算由编辑完成,采用球形地球与直线传播模型。