Remove ads
来自维基百科,自由的百科全书
地平坐标系(英语:Horizontal coordinate system),是天球坐标系统中的一种,以观测者所在地为中心点,所在地的地平线作为基础平面,将天球适当的分成能看见的上半球和看不见(被地球本身遮蔽)的下半球。上半球的顶点(最高点)称为天顶,下半球的顶点(最低点)称为地底。
地平坐标系统由两个夹角来定义一个天体位置的极座标:
因此地平坐标系有时也被称为高度/方位(Alt/Az)坐标系统。
地平坐标系统是固定在地球上而不是恒星,所以天体出现在天球上的高度和方位会随着时间,在天球上不停的改变。另一方面,因为基础平面是观测者所在地的地平面,所以相同的天体在相同的时间从不同的位置观察,也会有不同的高度和方位。
地平坐标系在测量天体的出没上非常的好用,当一个天体的高度为0°,就表示他位于地平线上。此时若其高度增加,就代表上升;若高度减少,便是下降。然而天球上所有天体的运动都受到由西向东的周日运动支配,所以与其笨拙的去观察高度是增加或减少,不如改为观察天体的方位更容易来判断是上升或是下降:
但在下面的特殊位置则例外:
需要注意的是:前面所考虑的只是理论上的几何地平,即不考虑地球大气层对天体位置的影响,让观测者的地平线完全以理想的海平面构成。因为地球有弧度,实际上看见的视地平面会随着观测者的高度增加而降低(出现负值)。另一方面大气层也会将地平线下半度的天体折射到地平线上。
只要知道观测者的地理坐标与时间,就可以将地平坐标转换成赤道坐标,或是反过来将赤道坐标转换成地平坐标。
在以下公式中,以代表方位,代表高度。
不管赤纬或地理纬度, 都是以北极点为+90°,在赤道是0°,南极点是-90°。
在地平座标转换前一般会先计算天体的本地时角 (Local Hour Angle, LHA) (或称地方时角)。天体的本地时角为观测时通过本地子午圈的天球经线的赤经值 与天体赤经 的差值 (), 也代表星体所在的赤经线与南方子午线在赤道面的夹角. 由于方位角是以南方(或北方)为基准, 所以用时角来转换到方位角颇为直觉. 上述的 正式名称为本地恒星时 (Local Sidereal Time, LST, ). 想像当地球以稳定的自转速度旋转时, 在每个恒星日, 南方子午线上会陆续通过赤经为 , ..., , ..., 东昇西落的星星, 就可想像 可以当成观测本地的一个时钟, 上面显示的时钟刻度就是本地恒星时 (), 换算成一小时 15 度, 也就是观测地经线相对于天球赤道起点 (春分点, ) 的旋转角度. 而天体的时角就代表从天体中天时刻到观测时刻所经历的时间或转动的角度. 显然, 本地恒星时由观测时间及观测地经度决定 (). 所以, 天体的时角也由天体的赤经及共同决定,故有时也会写成 或 , 代表赤经在特定观测时地的替代表示方式. 这也是为什么在空间座标转换时, (星体座标)会用去除时地标志的天体时角(及赤纬)来代替其赤经(及赤纬)的原因. 总结上述说明, 星体的时角与本地恒星时的关系及计算公式为:
如上所示, LST 可由 GST 加计本地地理经度求得. 其中, GST 为格林威治恒星时, 亦即 0 度经线上之观测站的 LST. 上述公式中, 241.3872 (度)代表在参考历元 2440000.5 JD (儒略日, Julian Date) (相当于1968/5/24.0) 时, 经过格林威治本初子午线的遥远恒星的赤经. 360.9856091 (度/太阳日)代表一天 (一个平太阳日) 之内地球转动的度数. 乘以 (用儒略日 jd 表示) 与 的差值, 代表至观测时间 总共新增的转动度数. 当然, 这些角度都要调整到 [0, 360] 或 [-180, +180] 的范围. 由 LST 就可以知道观测时通过本地子午线的星体的赤经了.
一般导航用的天文年鉴或历书 (almanac), 并无法把主要天文导航天体(如太阳, 月亮, 行星, 及约 57 颗导航用亮星), 在所有城市的本地时角, 都刊印出来, 仅能打印他们在格林威治所观测到的天体时角, 即格林威治时角 (Greenwich Hour Angle, GHA), 再由领航员从 的关系中加计经度推算出 . 因此, 上列公式也把 GHA 的相关式子列出来做为参考.
此外, 对同一观测目标 (), 在同一观测地 ()而言:
也就时说, 在同一观测地, 恒星时差()与天体时角差()是相同的, 且都跟观测时间差()成正比. 只不过恒星时钟与太阳时钟的时间长度及速度不一样, 地球公转一周看到远处恒星的次数比看到近处太阳的次数正好多1次. 所以, 恒星时钟比太阳时钟走得快一点. 若要把恒星时差换算成手表上的时差(平太阳时), 就必须多除以 这个系数 (). 在许多有关天文事件时间 (jd1) 或时差(duration) (Δ t) 的计算问题上 (如日出时刻、日落时刻、星体中天时刻、曙暮光始末时刻、日照时间或白天长度) [1], 要记得用这个比例常数来调整两种不同时间的刻度. 例如, 1 恒星日 () 的时间长度大约相当于 (平太阳时), 与一天(太阳日)的长度差了约 4 分钟.
赤道坐标转为地平坐标时, 可以透过以下的关系, 由天体的赤经 () 及赤纬 (), 求得天体的方位角() 及高度角 ()。
根据以上关系式, . 则可由 求得.
有种方式是把 相除后消去项,而化简为, 再用 来求 。但是, 使用的反正切函数的值域只在[-90, 90] 度之间, 无法完整涵盖 [0, 360] (或 [-180,+180]) 度的方位角. 而在 0 到 360 (或 [-180,+180]) 度之间, 值相同的角度有两个 (). 例如45°和225°是完全不同的方位, 但正切值相同。因此, 必须根据 及 的正负符号, 决定方位角落在哪个象限. 如果这些同值的角度落在非值域的第二及第三象限, 即 X 值为负时, 必须 +/-180 度, 才会得到正确的 . 若为 X=0 (Y/X 为无限大) 的特殊状况, 则依 的正负符号, 定义其方位角为 +90 或 -90 度。若 X, Y 皆为 0 (即天体在天顶), 则可依习惯定义方位。
其实不少编程语言(如 C, C++, Java, Python) 都有提供一个叫做 ATAN2(Y,X) (或 ATAN2(X,Y)) 的反三角函数 (atan2是已将象限纳入考量的反正切函数), 可算出 的值, 并根据 (X,Y) 的正负号判断所属象限, 从而决定 (X, Y) 向量与 X 轴的夹角, 让他的值域涵盖 360 度角. 这对决定方位角非常方便, 省掉自己编写程式码来判断象限的麻烦. 至于高度角 的求解, 可令第一个公式等号右边的值为 , 用 求 值即可, 不必再做调整. 因为, 的值域为正负 90 度, 正好对应地平线上下夹角 (这状况同样适用于之后在计算赤经赤纬时对应北南半球纬度).
需要特别注意的是, 上面计算出来的方位角 其实指的是以南方为0度向西递增的方位角, 而不是一般文献指称的, 以北方为0度, 向东递增的方位角. 这种一般文献上所称的 (北)方位角 若表示成, 则与上列计算出来的 , 或特意表示成 的 南方位角, 两者相差正好 180 度, 可以用 计算出来, 并调整到 0~360 度即可. 由于很多人不明白其间的差异, 因此由其他文献上抄录来的公式, 常因公式中某些项目的正负符号与其他来源(如维基网页)不同, 而误以为错误, 甚至错误更改维基百科的公式而不自知 (可察看本页历史编辑纪录). 其算出的结果也可能与预期有 180 度的差异. 所以, 参照不同来源公式时, 必须小心. 而之所以会有人定义这种南方为零的南方向角, 主要是一些北半球的观星者平时观测的星体以南方星体为主. 因此, 以南方为零度方位, 有其方便性.
上列公式并不容易理解其来由, 若移项重新整理, 并刻意以 提醒此方位角为南方位角, 则可得:
其矩阵形式则为:
,
其中, 最右边要被转换的行向量表示赤道极座标 或 () 投影在赤道面某选定直角座标的三个分量 , 等号左边的转换后所得行向量表示地平极座标 投影在地平面某选定直角座标的三个分量 . 中间的转换矩阵代表将赤道座标沿着子午线由天球北极 (Z 轴) 转向赤道面 (X轴) 转动 90- 度角的座标旋转矩阵. 这样的矩阵式说明了原公式的直觉意义, 对于需要时常计算的观星者,航海家或天文计算程式员而言比较不必硬记, 也较不容易弄错.
上列矩阵转换公式, 也让地平座标转赤道座标变得容易. 事实上, 只要把转换对象调换, 并进行逆转换即可. 换句话说, 前式的两个行向量只要互相调换, 并把原来的转换矩阵变成他的逆矩阵 (inverse matrix) 即可得到反向转换公式. 有趣的是, 座标转换的逆矩阵也是他的转置矩阵 (transpose matrix), 也就是行列互换的矩阵, 因此并不需要费力去求原转换矩阵的逆矩阵. 因此, 我们可以轻易得到:
亦即
, , .
其中, 为观测者所在经度 于观测时间 的本地恒星时.
前面已提到, 实用上有两种方位角, 前面计算的其实是南方位角 . 一般官方文献所提的方位角为北方位角 . 为了避免混淆, 以下将使用 时的座标转换公式也一并列出, 以便相互对照.
赤道转地平, 求 的方法除了用先前方法算出 再加 180 度之外, 也可以将原来的转换公式中的 跟 等号右侧方程式都加负号, 并把等号左侧的 改成 即可. 其结果是:
, , 或者 , . 同时:
两者相除后,除正负号的区别外,形式完全一样,已无法区分这里的方位角是南方位角或北方位角。且已失去判断象限的讯息,必须由分子分母的正负来辅助判断。这跟之前讨论如何由 (即 ) 求 的情况一样。
有兴趣者可以把他转成矩阵转换式, 会发现这样的转换是经过两道转换手续, 即先转成原先的南地平, 再把 X 轴转 180 度, 也就是 X 值跟 Y 值都取负号.
要得到使用 时的地平转赤道座标转换公式, 只要将 代入原来的南地平转赤道的转换公式即可. 此代换会得出, , , 因此, 有以下转换公式:
, , .
其中, 为观测者所在经度 于观测时间 的本地恒星时. 注意, 其实与 是相同的两组赤道座标,只是以 及 表达时,形式不同而已。
比较 跟 会发现, 两者的转换公式长得完全一样, 不同的只是符号的代换. 把 跟 分别与 跟 互相代换就会得到另一组转换公式. 这是因为赤道转北地平的转换矩阵(即矩阵式中的两个矩阵相乘)是对称矩阵, 所以它的逆矩阵 (已知等于转置矩阵) 跟原转换矩阵是一样的. 所以, 除了符号互相替换之外, 公式的形式完全相同. 这个有趣的结果可以有两个应用. 第一,是可以由两个方向的转换矩阵或转换公式的形式是否一样来判断公式里的方位角到底是不是以北方为零度的方位角. 第二,如果采用 为方位角, 则撰写转换程式码时其实只需要写一个函数.
赤道座标与地平座标之转换, 牵涉到观测时间 或 或 , 观测位置(), 观测天体座标 () 或 ( 和观测者地面量测的视角及数据. 透过这些相依关系, 只要固定某些变数或进行相关测量, 就可以求得其他感兴趣的变数, 进行预测或量测. 这在天体追踪, 观测活动规划, 个人位置定位, 天文导航 (celestial navigation) 等方面, 应用极为广泛. 举例而言:
天体的出没及持续时间,可以由前述中的高度角及时角的关系推导出来。其中, 时角隐藏时间()、观测地经度()、观测天体赤经(), 是计算与时间相关的问题时,会被检视的变数。
天体出没时,按理说其高度角应该为。但由于不同天体视直径()及观测地的气候条件(如大气折射效应, Refraction)或周遭地理状况之不同(如海拔及障碍物),未必为零。
所以,计算天体(含太阳)的实体或虚像出现或隐没时,高度角一般由以下几个可选的参数决定[1]:
其中,
曙暮光的标准并非依据客观的太阳物理数据,而是依据不同人群受日光出现之影响程度所订定的经验标准。因此有民用曙暮光 (civil twilight)、航海曙暮光 (nautical twilight) 及天文曙暮光 (astronomical twilight) 的订定。天文观测一般不愿受日光影响,所以希望太阳于地平线下 18 度以下的期间才进行观测。船只航行只要清楚看得到海天之交的弧线就可以安心航行,所以用较宽松的海平面下 12 度,作为航海曙暮光的标准。至于一般人民,只要天将亮夜未深之际进行作息即可,所以民用曙暮光的标准是更宽松的海平面下 6 度。
给定出没时的天体高度角后,假设 为对应的时角,则计算日(星)出、日(星)落、曙暮光等事件的时间可摘要如下:
此处, 天体出现时的天体时角为,隐没时的时角为。由于 且 ,故满足 且在 [-180,180] 或 [0,360] 的角度有两个, 即 。因此,其中一个为天体出现时的时角,另一个为隐没时的时角。
若无永昼永夜之类的极端情况,则可由分别先求出星体出没的本地恒星时,及格林威治恒星时。并以当日子夜零点的GST () 为基准,求恒星时差异。当然也可以由两个 LST 求恒星时差异。随后将恒星时差异调为太阳时差异。必要时将以日为单位的太阳时差异,以每日24小时换算为时:分:秒,就可求出星体出没的时间。以下公式概括所有步骤:
天体出现在地平线上的时间,即从上升到下沉所经历的时间,其实就是相当于两倍 的时间。如果天体是指太阳的话,这就是白天的时间或日照时间。对其他天体而言,就是最长可观测的时间 (不考虑曙暮光的影响的话)。
天体出没的方位角,显然只跟天体所在的赤纬平行圈及观测地纬度有关,与出没时间无关。所有赤纬相同的天体都在同一天球赤纬平行圈上,该平行圈与地平线交于相同的两点,其方位即天体出或没时的方位。给定出没时的天体高度角 及观测地纬度 ,假设 及 分别为对应的南方位角及北方位角,则该方位角可计算如下 (如果有出没状况的话):
另外, 也可以直接由 的公式, 求解日出日落时的北方位角,
以上, 的公式,不只可以求解日升日落(星升星落)(即 ) 的状况,也可以求一般状况下的天体高度及方位角。
这时常应用在天文航海或天体导航 (Celestial Navigation) 中的截距法 (Intercept Method)。在此法中,领航员为了确定自身的经纬度,必须计算在某个假设的地理位置 (AP, Assumed Position) 上,导航天体会被观测到的计算高度 (通常记为 , 相当于公式中的 ) 及方位角 (北方位角通常记为 , 相当于公式中的 ),再与实际观测到的高度角 (通常记为 ) 比较,来决定应该将 AP 往天体的地理位置 (GP, Geographic Position, 即天体在地球表面的星下点位置) 方向前进或后退若干海里,以得到航行器本身实际的经纬度。
如果观测到的高度角大于在 AP 的计算高度 (),代表航行器比假设位置更靠近天体。这就好像面对远方导航的灯塔(天体)时,越靠近灯塔,所测量到的塔顶高度角越大。在此情况下,就应把 AP 的座标点,沿计算方位角 Zn 的方向,往天体(GP)方向推进,以推得航行器本身的位置。反之,若 ,就应把 AP 的座标往 Zn 的相反方向后退,远离天体,来推得航行器本身的位置。这种直观的位置修正法则,一般简记为 HoMoTo 法则 (if HO is MOre than Hc, move AP TOward GP for a fix.)
以上所选的 AP 通常是利用其他航位推测法 (如航速及水流/风的速度) 所获得的大略位置,故与实际位置所推算的导航天体方位角及高度角不会相差太多。但透过星体导航,可以得到更精确的位置,因为星体的位置可以精密推算得到,变成所有航行器都可以参考的 '灯塔'。至于 AP 往前推进(或后退)的距离,称为截距距离 (intercept distance),可以由 推得。这是因为观测地与 GP 的角距离为 , 等于天体的天顶距 (ZD, Zenith Distance)。所以,观测地与 AP 的角距离为 。而 1 海里 (nautical mile, nm) 的定义为地表航行 1 分角的平均距离 (约 1.852 km)。故可换算截距距离为高度角差距乘以每一分角度一海里(或每一度60海里)。
火箭、人造卫星、弹道导弹、航天飞机、太空船也可以用明亮的恒星执行导航任务。执行此类导航的装置通常称为恒星追踪仪 (Star Tracker),或恒星追踪器、跟星仪、星象仪、星光探测器等等。多数的恒星追踪仪内部存有已知位置的亮星赤经赤纬数据库。航行途中,追踪仪透过一个或一个以上的相机镜头或望远镜获取视野内的星空影像,并将所拍摄的影像与数据库里的亮星比对,辨识出视野内的恒星或星座。并从镜头感光器上的恒星亮点与镜头的相对角度,求出飞行器的姿态(altitude)及定位。就像人类从天空辨识到北斗七星及相对高度方位,从而了解所在位置一样。
对于固定轨道的弹道导弹而言,则可预先计算每个预订时刻在预定位置可以观察到的亮星,及其相对于飞行路径的假想地平面的高度角及方位角。并从实际的观测值与计算值的差距,产生错误修正讯号,以修正其飞行轨迹,最终将导弹导引至目的地。[3] 其原理类似前述的截距法。但所有修正都由飞行电脑及导航机构自动完成。
在实际运用上,天体导航有时会与其他导航系统结合,截长补短,以提高导航的精确度。例如,恒星追踪仪可能跟惯性导航系统 (INS, Inertial Navigation System) 结合,形成所谓的星光惯性导航系统 (Stellar Inertial Navigation Systems)。
在地平坐标系统中,有好几种方法可以计算太阳的视位置。
完整和精确的计算方法可以参考比利时天文学家简米斯的天文计算(Astronomical Algorithms)
下面是一种简单的近似计算法的例子:
已知:
以下的公式可以算出太阳的赤纬:
其他移动的天体,如行星、彗星、小行星、月亮、行星的卫星、人造卫星、太空船等,也可以透过解开普勒方程式,求得它们在其轨道面的即时位置,再透过适当的座标转换,将轨道面座标转换到赤道座标,再转换到水平座标。这样就可以预测它们的方位角及高度角,再以人工或自动的方式,加以追踪。(太阳系行星位置计算及行星轨道要素 [4] 可参考 NASA 外部链接.)
Seamless Wikipedia browsing. On steroids.
Every time you click a link to Wikipedia, Wiktionary or Wikiquote in your browser's search results, it will show the modern Wikiwand interface.
Wikiwand extension is a five stars, simple, with minimum permission required to keep your browsing private, safe and transparent.