地平坐標系(英语:Horizontal coordinate system),是天球坐標系統中的一種,以觀測者所在地為中心點,所在地的地平線作為基礎平面,將天球適當的分成能看見的上半球和看不見(被地球本身遮蔽)的下半球。上半球的頂點(最高點)稱為天頂,下半球的頂點(最低點)稱為地底。
地平坐標系統由兩個夾角來定義一個天體位置的極座標:
- 高度角(Altitude, Alt)或仰角又稱地平緯度,是天體和觀測者所在地的地平線的夾角。
- 方位角(Azimuth, Az)又稱地平經度,是沿着地平線測量的角度. 一般文獻指稱的方位角是以正北方為0度起點, 順時鐘向東方測量. 但對於某些觀星者或航海家而言, 定義以南方為0度起點的方向角, 有其方便性. 因此以下將以 A_N 及 A_S 分別代表(以正北為0度的)北方位角及(以正南為0度的)南方位角。
因此地平坐標系有時也被稱為高度/方位(Alt/Az)坐標系統。
簡略的觀測
地平坐標系統是固定在地球上而不是恆星,所以天體出現在天球上的高度和方位會隨著時間,在天球上不停的改變。另一方面,因為基礎平面是觀測者所在地的地平面,所以相同的天體在相同的時間從不同的位置觀察,也會有不同的高度和方位。
地平坐標系在測量天體的出沒上非常的好用,當一個天體的高度為0°,就表示他位於地平線上。此時若其高度增加,就代表上升;若高度減少,便是下降。然而天球上所有天體的運動都受到由西向東的周日運動支配,所以與其笨拙的去觀察高度是增加或減少,不如改為觀察天體的方位更容易來判斷是上升或是下降:
*當天體的方位在0°~180°之間(北方—東方—南方,亦即子午線之東)是上升。
*當天體的方位在180°~360°之間(南方—西方—北方,亦即子午線之西)是下降。
但在下面的特殊位置則例外:
*在北極點,因為天頂就是北天極,所有的方向都是南方,所以無法定出方位,但這並不造成問題,因為所有天體的高度無論任何時間都不會改變,即既不升高也不降低,只繞北極星以逆時針轉動。(头朝下感觉天星是顺时针转,抬头望天,才看见天星逆时针转)
*在南極,地面上所有方向都是北方,也會有與北極相同情況,只是所有星星皆繞天頂的南天極順時針轉動。
*在赤道,位於極點的天體會固定不動的永遠停留在地平線上的那一個點。(但實際上由於天極很接近地平線,在該處天體未必能直接看到)
需要注意的是:前面所考慮的衹是理論上的幾何地平,即不考慮地球大氣層對天體位置的影響,讓觀測者的地平線完全以理想的海平面構成。因為地球有弧度,實際上看見的視地平面會隨著觀測者的高度增加而降低(出現負值)。另一方面大氣層也會將地平線下半度的天體折射到地平線上。
與赤道坐標系的互換
只要知道觀測者的地理坐標與時間,就可以將地平坐標轉換成赤道坐標,或是反過來將赤道坐標轉換成地平坐標。
在以下公式中,以A代表方位,a代表高度。
以\alpha表示赤經,\delta表示赤緯,H表示時角。
φ為觀測者所在地的纬度。
不管赤緯或地理緯度, 都是以北極點為+90°,在赤道是0°,南極點是-90°。
天體時角與本地恆星時
在地平座標轉換前一般會先計算天體的本地時角 (Local Hour Angle, LHA) (或稱地方時角)。天體的本地時角H為觀測時通過本地子午圈的天球經線的赤經值\alpha_L 與天體赤經 \alpha的差值 (H \equiv \alpha_L - \alpha), 也代表星體所在的赤經線與南方子午線在赤道面的夾角. 由於方位角是以南方(或北方)為基準, 所以用時角來轉換到方位角頗為直覺. 上述的 \alpha_L 正式名稱為本地恆星時 (Local Sidereal Time, LST, LST \equiv \alpha_L). 想像當地球以穩定的自轉速度旋轉時, 在每個恆星日, 南方子午線上會陸續通過赤經為 0^h, 1^h, ..., \alpha^h, ..., \alpha_L^h, ..., 24^h (\equiv 0^h) 東昇西落的星星, 就可想像 \alpha_L 可以當成觀測本地的一個時鐘, 上面顯示的時鐘刻度就是本地恆星時 (LST), 換算成一小時 15 度, 也就是觀測地經線相對於天球赤道起點 (春分點, \alpha_0 = 0^h) 的旋轉角度. 而天體的時角就代表從天體中天時刻到觀測時刻所經歷的時間或轉動的角度. 顯然, 本地恆星時由觀測時間t及觀測地經度\lambda決定 (\alpha_L=\alpha_L(t, \lambda)). 所以, 天體的時角也由天體的赤經\alpha及t, \lambda共同決定,故H有時也會寫成H (\alpha, t, \lambda) 或 H (\alpha), 代表赤經在特定觀測時地的替代表示方式. 這也是為什麼在空間座標轉換時, (星體座標)會用去除時地標誌的天體時角(及赤緯)來代替其赤經(及赤緯)的原因. 總結上述說明, 星體的時角與本地恆星時的關係及計算公式為:
:\begin{align}
H(\alpha, t, \lambda) &= LST(t, \lambda) - \alpha &&\equiv \alpha_L - \alpha \\
LST(jd,\lambda) &= GST(jd) + \lambda &&\equiv \alpha_L (jd, \lambda)\\
GST(jd) &= 241.3872 + 360.9856091 \times (jd-2440000.5) &&\equiv LST(jd, 0^\circ) \\
GHA(\alpha, t) &\equiv H(\alpha, t, 0^\circ) &&= GST(t) - \alpha \\
H(\alpha, t, \lambda) &= GHA(\alpha, t) + \lambda &&= GST(t) + \lambda - \alpha
\end{align}
如上所示, LST 可由 GST 加計本地地理經度求得. 其中, GST 為格林威治恆星時, 亦即 0 度經線上之觀測站的 LST. 上述公式中, \theta_0 =241.3872 (度)代表在參考曆元 T_0 = 2440000.5 JD (儒略日, Julian Date) (相當於1968/5/24.0) 時, 經過格林威治本初子午線的遙遠恆星的赤經.
R_0 =360.9856091 (度/太陽日)代表一天 (一個平太陽日) 之內地球轉動的度數. 乘以 t (用儒略日 jd 表示) 與 T_0 的差值, 代表至觀測時間 t 總共新增的轉動度數.
當然, 這些角度都要調整到 [0, 360] 或 [-180, +180] 的範圍. 由 LST 就可以知道觀測時通過本地子午線的星體的赤經了.
一般導航用的天文年鑑或曆書 (almanac), 並無法把主要天文導航天體(如太陽, 月亮, 行星, 及約 57 顆導航用亮星), 在所有城市的本地時角, 都刊印出來, 僅能列印他們在格林威治所觀測到的天體時角, 即格林威治時角 (Greenwich Hour Angle, GHA), 再由領航員從 LHA = GHA + \lambda 的關係中加計經度推算出 LHA. 因此, 上列公式也把 GHA 的相關式子列出來做為參考.
- 單位與慣例: 上列公式中的經度 \lambda, 以東經為正, 西經為負, 故亦稱為東經度, 可表示為 \lambda_E. 某些地區習慣取西經為正, 東經為負, 稱為西經度. 為示區別, 可表示為 \lambda_W (\lambda_W = - \lambda_E). 本頁所稱經度概以東經度為主, 即 \lambda \equiv \lambda_E. 而公式中的夾角皆以角度(degree)為單位, 若改成弳度 (radian) 或時間, 則可用 360^\circ \equiv 2\pi \text{ rad} \equiv 24^h 換算.
此外, 對同一觀測目標 (\alpha), 在同一觀測地 (\lambda)而言:
:\begin{align}
GST(jd) &= \theta_0 + R_0 \times (jd-T_0) &&= LST(jd,\lambda) - \lambda &&= H(\alpha, t,\lambda) - \alpha -\lambda \\
\Delta ST &\triangleq GST(jd1) - GST(jd0) &&= LST(jd1,\lambda) - LST(jd0,\lambda) \\
&= H(\alpha, jd1, \lambda) - H(\alpha, jd0, \lambda) && \triangleq \Delta H \\
\Delta ST &= R_0 \times (jd1-jd0) &&= \Delta t \times R_0 &&= \Delta H \\
\Delta t &\triangleq jd1-jd0 &&= \Delta ST / R_0 &&= \Delta H / R_0 \\
jd1 &= jd0 + \Delta ST / R_0 &&= jd0 + \Delta H /R_0 \\
\end{align}
也就時說, 在同一觀測地, 恆星時差(\Delta ST)與天體時角差(\Delta H)是相同的, 且都跟觀測時間差(\Delta t)成正比. 只不過恆星時鐘與太陽時鐘的時間長度及速度不一樣, 地球公轉一周看到遠處恆星的次數比看到近處太陽的次數正好多1次. 所以, 恆星時鐘比太陽時鐘走得快一點. 若要把恆星時差換算成手錶上的時差(平太陽時), 就必須多除以 R_0 這個係數 (R_0 \approx (365.25+1)/365.25 \times 360^o). 在許多有關天文事件時間 (jd1) 或時差(duration) (Δ t) 的計算問題上 (如日出時刻、日落時刻、星體中天時刻、曙暮光始末時刻、日照時間或白天長度)
, 要記得用這個比例常數來調整兩種不同時間的刻度. 例如, 1 恆星日 (\Delta ST = 360^0) 的時間長度大約相當於 \Delta t = 23^h 56^m 4.09^s (平太陽時), 與一天(太陽日)的長度差了約 4 分鐘.
赤道坐標轉為地平坐標
赤道坐標轉為地平坐標時, 可以透過以下的關係, 由天體的赤經 (\alpha) 及赤緯 (\delta), 求得天體的方位角(A) 及高度角 (a)。
:\begin{align}
Z_h &= \sin a &&= \sin \phi \cdot \sin \delta + \cos \phi \cdot \cos \delta \cdot \cos H \\
X_h &= \cos A \cdot \cos a &&= -\cos \phi \cdot \sin \delta + \sin \phi \cdot \cos \delta \cdot \cos H \\
Y_h &= \sin A \cdot \cos a &&= \cos \delta \cdot \sin H \\
\end{align}
根據以上關係式, a = arcsin(Z_h). A 則可由 X_h, Y_h 求得.
有種方式是把 Y_h, X_h 相除後消去\cos a項,而化簡為\tan A = Y_h/X_h, 再用 \arctan (Y_h/X_h) 來求 A。但是, \arctan(Y_h/X_h) 使用的反正切函數的值域只在[-90, 90] 度之間, 無法完整涵蓋 [0, 360] (或 [-180,+180]) 度的方位角. 而在 0 到 360 (或 [-180,+180]) 度之間, \tan 值相同的角度有兩個 (\tan A = \tan (A \mp 180) = Y_h/X_h). 例如45°和225°是完全不同的方位, 但正切值相同。因此, 必須根據X_h 及 Y_h 的正負符號, 決定方位角落在哪個象限. 如果這些同值的角度落在非值域的第二及第三象限, 即 X 值為負時, \arctan(Y/X) 必須 +/-180 度, 才會得到正確的 A. 若為 X=0 (Y/X 為無限大) 的特殊狀況, 則依 Y_h 的正負符號, 定義其方位角為 +90 或 -90 度。若 X, Y 皆為 0 (即天體在天頂), 則可依習慣定義方位。
:\begin{align}
\tan A &
= \frac{\cos \delta \cdot \sin H}{-\cos \phi \cdot \sin \delta + \sin \phi \cdot \cos \delta \cdot \cos H}
&& = \frac{Y_h}{X_h} \\
&= \frac{\sin H}{-\cos \phi \cdot \tan \delta + \sin \phi \cdot \cos H} \\
A_0 &\triangleq \arctan \frac{Y_h}{X_h} && \in [-90,90] \\
A &= \begin{cases}
A_0 & \qquad X_h > 0 & \in [-90,90] \\
A_0 + 180^\circ & \qquad Y_h \ge 0 , X_h 0 , X_h = 0 \\
- 90^\circ & \qquad Y_h
其實不少程式語言(如 C, C++, Java, Python) 都有提供一個叫做 ATAN2(Y,X) (或 ATAN2(X,Y)) 的反三角函數 (atan2是已將象限納入考量的反正切函數), 可算出 \arctan(Y/X) 的值, 並根據 (X,Y) 的正負號判斷所屬象限, 從而決定 (X, Y) 向量與 X 軸的夾角, 讓他的值域涵蓋 360 度角. 這對決定方位角非常方便, 省掉自己編寫程式碼來判斷象限的麻煩. 至於高度角 a 的求解, 可令第一個公式等號右邊的值為 Z, 用 a = \arcsin(Z) 求 a 值即可, 不必再做調整. 因為, \arcsin 的值域為正負 90 度, 正好對應地平線上下夾角 (這狀況同樣適用於之後在計算赤經赤緯時對應北南半球緯度).
兩種方位角
需要特別注意的是, 上面計算出來的方位角 A 其實指的是以南方為0度向西遞增的方位角, 而不是一般文獻指稱的, 以北方為0度, 向東遞增的方位角. 這種一般文獻上所稱的 (北)方位角 若表示成A_N, 則與上列計算出來的 A, 或特意表示成 A_S 的 南方位角, 兩者相差正好 180 度, 可以用 A_N = A_S + 180 計算出來, 並調整到 0~360 度即可. 由於很多人不明白其間的差異, 因此由其他文獻上抄錄來的公式, 常因公式中某些項目的正負符號與其他來源(如維基網頁)不同, 而誤以為錯誤, 甚至錯誤更改維基百科的公式而不自知 (可察看本頁歷史編輯紀錄). 其算出的結果也可能與預期有 180 度的差異. 所以, 參照不同來源公式時, 必須小心. 而之所以會有人定義這種南方為零的南方向角, 主要是一些北半球的觀星者平時觀測的星體以南方星體為主. 因此, 以南方為零度方位, 有其方便性.
赤平座標轉換之矩陣轉換式
上列公式並不容易理解其來由, 若移項重新整理, 並刻意以 A_s 提醒此方位角為南方位角, 則可得:
:\begin{align}
X_h &= \cos a \cdot \cos A_S &&= \sin \phi \cdot \cos \delta \cdot \cos H -\cos \phi \cdot \sin \delta \\
Y_h &= \cos a \cdot \sin A_S &&= \cos \delta \cdot \sin H \\
Z_h &= \sin a &&= \cos \phi \cdot \cos \delta \cdot \cos H + \sin \phi \cdot \sin \delta \\
\end{align}
其矩陣形式則為:
:
\begin{bmatrix}
X_h \\
Y_h \\
Z_h
\end{bmatrix}
=
\begin{bmatrix}
\cos a \cdot \cos A_S \\
\cos a \cdot \sin A_S \\
\sin a \\
\end{bmatrix}
=
\begin{bmatrix}
\sin \phi & 0 & - \cos \phi \\
0 & 1 & 0 \\
\cos \phi & 0 & \sin \phi
\end{bmatrix}
\times
\begin{bmatrix}
\cos \delta \cdot \cos H \\
\cos \delta \cdot \sin H \\
\sin \delta
\end{bmatrix}
A_S = atan2(Y_h, X_h), a = \arcsin(Z_h)
其中, 最右邊要被轉換的行向量表示赤道極座標 (\alpha, \delta) 或 (H, \delta) (H = LST - \alpha) 投影在赤道面某選定直角座標的三個分量 (X_e, Y_e, Z_e), 等號左邊的轉換後所得行向量表示地平極座標 (A_S, a) 投影在地平面某選定直角座標的三個分量 (X_h, Y_h, Z_h). 中間的轉換矩陣代表將赤道座標沿著子午線由天球北極 (Z 軸) 轉向赤道面 (X軸) 轉動 90-\phi 度角的座標旋轉矩陣. 這樣的矩陣式說明了原公式的直覺意義, 對於需要時常計算的觀星者,航海家或天文計算程式員而言比較不必硬記, 也較不容易弄錯.
地平坐標轉為赤道坐標
上列矩陣轉換公式, 也讓地平座標轉赤道座標變得容易. 事實上, 只要把轉換對象調換, 並進行逆轉換即可. 換句話說, 前式的兩個行向量只要互相調換, 並把原來的轉換矩陣變成他的逆矩陣 (inverse matrix) 即可得到反向轉換公式. 有趣的是, 座標轉換的逆矩陣也是他的轉置矩陣 (transpose matrix), 也就是行列互換的矩陣, 因此並不需要費力去求原轉換矩陣的逆矩陣. 因此, 我們可以輕易得到:
:
\begin{bmatrix}
X_e \\
Y_e \\
Z_e
\end{bmatrix}
=
\begin{bmatrix}
\cos \delta \cdot \cos H \\
\cos \delta \cdot \sin H \\
\sin \delta
\end{bmatrix}
=
\begin{bmatrix}
\sin \phi & 0 & \cos \phi \\
0 & 1 & 0 \\
- \cos \phi & 0 & \sin \phi
\end{bmatrix}
\times
\begin{bmatrix}
\cos a \cdot \cos A_S \\
\cos a \cdot \sin A_S \\
\sin a \\
\end{bmatrix}
亦即
:\begin{align}
X_e &= \cos \delta \cdot \cos H &&= \sin \phi \cdot \cos a \cdot \cos A_S + \cos \phi \cdot \sin a \\
Y_e &= \cos \delta \cdot \sin H &&= \cos a \cdot \sin A_S \\
Z_e &= \sin \delta &&= - \cos \phi \cdot \cos a \cdot \cos A_S + \sin \phi \cdot \sin a \\
\end{align}
H = atan2(Y_e, X_e), \alpha = LST(t,\lambda) - H(\alpha), \delta = \arcsin(Z_e)
.
其中, LST(t, \lambda) 為觀測者所在經度 \lambda 於觀測時間 t 的本地恆星時.
方位角以北為零時的赤平座標相互轉換
前面已提到, 實用上有兩種方位角, 前面計算的其實是南方位角 A_S. 一般官方文獻所提的方位角為北方位角 A_N. 為了避免混淆, 以下將使用 A_N時的座標轉換公式也一併列出, 以便相互對照.
赤道轉北地平
赤道轉地平, 求 A_N 的方法除了用先前方法算出 A_S 再加 180 度之外, 也可以將原來的轉換公式中的 X_h 跟 Y_h等號右側方程式都加負號, 並把等號左側的 A_S 改成 A_N 即可. 其結果是:
:\begin{align}
X_{hN} &= \cos a \cdot \cos A_N &&= - \sin \phi \cdot \cos \delta \cdot \cos H +\cos \phi \cdot \sin \delta &&= - X_h\\
Y_{hN} &= \cos a \cdot \sin A_N &&= - \cos \delta \cdot \sin H &&= - Y_h\\
Z_{hN} &= \sin a &&= \cos \phi \cdot \cos \delta \cdot \cos H + \sin \phi \cdot \sin \delta &&= + Z_h\\
\end{align}
A_N = atan2(Y_{hN}, X_{hN}), a = \arcsin(Z_{hN}), 或者
A_N = atan2(-Y_{h}, -X_{h}), a = \arcsin(Z_{h}).
同時:
:\begin{align}
\tan A_N &= \frac{Y_{hN}}{X_{hN}} = \frac{-Y_{h}}{-X_{h}} = \frac{Y_{h}}{X_{h}} = \tan A_S = \tan A \\
\end{align}
兩者相除後,除正負號的區別外,形式完全一樣,已無法區分這裡的方位角是南方位角或北方位角。且已失去判斷象限的訊息,必須由分子分母的正負來輔助判斷。這跟之前討論如何由 \tan A (即 \tan A_S) 求 A 的情況一樣。
有興趣者可以把他轉成矩陣轉換式, 會發現這樣的轉換是經過兩道轉換手續, 即先轉成原先的南地平, 再把 X 軸轉 180 度, 也就是 X 值跟 Y 值都取負號.
:
\begin{bmatrix}
X_{hN} \\
Y_{hN} \\
Z_{hN}
\end{bmatrix}
=
\begin{bmatrix}
\cos a \cdot \cos A_N \\
\cos a \cdot \sin A_N \\
\sin a \\
\end{bmatrix}
=
\begin{bmatrix}
-1 & 0 & 0 \\
0 & -1 & 0 \\
0 & 0 & 1
\end{bmatrix}
\times
\begin{bmatrix}
\sin \phi & 0 & - \cos \phi \\
0 & 1 & 0 \\
\cos \phi & 0 & \sin \phi
\end{bmatrix}
\times
\begin{bmatrix}
\cos \delta \cdot \cos H \\
\cos \delta \cdot \sin H \\
\sin \delta
\end{bmatrix}
北地平轉赤道
要得到使用 A_N 時的地平轉赤道座標轉換公式, 只要將 A_S = A_N - 180代入原來的南地平轉赤道的轉換公式即可. 此代換會得出, \cos(A_S) = \cos(A_N - 180) = -\cos(A_N), \sin(A_S) = \sin(A_N - 180) = -\sin(A_N), 因此, 有以下轉換公式:
:\begin{align}
X_{eN} &= \cos \delta \cdot \cos H &&= - \sin \phi \cdot \cos a \cdot \cos A_N + \cos \phi \cdot \sin a &&=X_e\\
Y_{eN} &= \cos \delta \cdot \sin H &&= - \cos a \cdot \sin A_N &&=Y_e\\
Z_{eN} &= \sin \delta &&= \cos \phi \cdot \cos a \cdot \cos A_N + \sin \phi \cdot \sin a &&=Z_e\\
\end{align}
H = atan2(Y_{eN}, X_{eN}), \alpha = LST(t,\lambda) - H(\alpha), \delta = \arcsin(Z_{eN}).
其中, LST(t, \lambda) 為觀測者所在經度 \lambda 於觀測時間 t 的本地恆星時.
注意,(X_{eN}, Y_{eN}, Z_{eN}) 其實與 (X_{e}, Y_{e}, Z_{e}) 是相同的兩組赤道座標,只是以 A_S 及 A_N 表達時,形式不同而已。
比較 (X_{eN}, Y_{eN}, Z_{eN}) 跟 (X_{hN}, Y_{hN}, Z_{hN}) 會發現, 兩者的轉換公式長得完全一樣, 不同的只是符號的代換. 把 H 跟 \delta 分別與 A_N 跟 a 互相代換就會得到另一組轉換公式. 這是因為赤道轉北地平的轉換矩陣(即矩陣式中的兩個矩陣相乘)是對稱矩陣, 所以它的逆矩陣 (已知等於轉置矩陣) 跟原轉換矩陣是一樣的. 所以, 除了符號互相替換之外, 公式的形式完全相同. 這個有趣的結果可以有兩個應用. 第一,是可以由兩個方向的轉換矩陣或轉換公式的形式是否一樣來判斷公式裡的方位角到底是不是以北方為零度的方位角. 第二,如果採用 A_N 為方位角, 則撰寫轉換程式碼時其實只需要寫一個函數.
赤道座標與地平座標轉換之應用
赤道座標與地平座標之轉換, 牽涉到觀測時間 t 或 GST(t) 或 LST(t,\lambda), 觀測位置(\lambda, \phi), 觀測天體座標 (\alpha, \delta) 或 (H(\alpha, t, \lambda), \delta) 和觀測者地面量測的視角及數據(a, A). 透過這些相依關係, 只要固定某些變數或進行相關測量, 就可以求得其他感興趣的變數, 進行預測或量測. 這在天體追蹤, 觀測活動規劃, 個人位置定位, 天文導航 (celestial navigation) 等方面, 應用極為廣泛. 舉例而言:
- 天體位置追蹤: 由觀測時地 t, (\lambda, \phi) 及天體座標 (\alpha, \delta) 求其高度及方位 (a, A).
** 標定天體的方位及高度, 以規劃適當的觀測場地, 設置觀測儀器, 進行觀測活動.
** 太陽位置追蹤及太陽能應用:
*** 將地面或太空船上的太陽光電模組、太陽熱能模組對準陽光方向, 以獲得最大日照強度及能量.
*** 日蝕時, 由太陽及月球相對方位及高度繪製模擬過程. 提供太陽能設備特殊處理所需的位置資訊.
- 事件時程預測: 由觀測地位置 (\lambda, \phi) 及天體座標 (\alpha, \delta) 求天體特殊事件發生的時刻 t, 持續時間 \Delta t 及方位高度 (a, A).
** 可由Z_h (\equiv Z_{hN}) 推算時間 (隱藏在時角變數H及恆星時GST內), 由 Z_{e} 或 (Z_{eN}) 推算方位.
** 星起(rise)、星落(set)、中天(transit)時刻、可觀測時間(duration).
** 曙光開始及暮光終止(twilight)時間: 決定最佳觀測或活動時間窗口.
** 太陽活動(日出日落)時間及太陽能應用:
*** 計算日出(sunrise)、日落(sunset)、中天時間、方位, 及時啟動及關閉太陽能裝置, 適時啟動備用儲能系統或其他發電系統.
*** 計算日照時間, 預測太陽光電全年發電量, 規劃與其他發電裝置及儲能系統的最佳搭配.
*** 由太陽及月亮位置, 預測日蝕發生時間. 提前預警太陽能設備進行特殊處理所需的時間資訊.
*** 預測人造衛星進入地影時間.
- 觀測者地理位置定位(positioning)及天文導航 (celestial navigation) 應用:
** 定位: 由一個或多個已知星體 (\alpha, \delta) (或 (H(\alpha, t, \lambda), \delta)) 在不同時間 t 的量測高度與方位 (a, A) 決定觀測者(或船舶及飛行器)所在地理位置 (\lambda, \phi). 理論上, 同一時間三個星體或同一星體三個不同時間的測量, 可以決定觀測者所在的經緯度.
*** 原理: 所有位置中, 會觀測到天體的高度角為 a 的所有經緯度的集合, 是一個以天體星下點地理位置 GP 為中心(半徑為餘高, co-altitude, 90^\circ - a) 的圓圈, 稱為等高度角圈 (簡稱等高圈) (circle of equal altitude) 或稱該觀測高度角對應的位置圈. 觀測者的經緯度必在其中. 若觀測兩個天體的兩個高度角, 可以繪出兩個位置圈, 其交點僅剩兩個, 觀測者的位置必在其中之一. 通常這時就可以輕易排除其一, 而猜到正確經緯度. 若再觀察另一天體並繪出其位置圈, 則三個天體的三個位置圈的交點為唯一, 即為觀測者所在的經緯度.
*** 公式: 利用地平座標中的參數 Z_h 的角度與位置關係, 給定高度角 a 及時間 (t, 即可算出天體等高度角圈 (circle of equal altitude), 即位置圈 (Circle of Positions, COP) 中的所有經緯度. 位置圈中的一小段弧線, 約可視為一直線, 稱為位置線 (Line of Positions, LOP). LOP 是天文導航中常用來決定船舶或飛機位置的重要資訊. 用越多的 LOP 的交點所決定的經緯度越精確.
** 導航: 由假設的地理位置 (AP, Assumed Position) (\lambda_{AP}, \phi_{AP}) 計算天體在 AP 應有的高度角及方位 (a_{AP}, A_{AP}), 並將計算值與實際量測到的天體高度角 (a) 比較, 推定航行器的適當航向及航程. 這是用 LOP 概念的一種直覺導航方法, 稱為截距法 (intercept method, IM). 可用來修正航位推測 (Dead reckoning) 的誤差.
- 天體地面位置投影及軌跡推測:
** 計算天體於特定時間 t 向地心投影的地面投影點(稱為星下點 Substellar point)的地理位置 (GP, Geographic Position) (以經緯度表示).
*** GP(\alpha, \delta, t) \triangleq (\lambda_{GP}, \phi_{GP}) = ( - GHA(\alpha, t), \delta) = (\alpha - GST(t), \delta).
*** GHA: 天體的格林威治時角 (向西為正向東為負).
*** 此處經度(\lambda_{GP})為東經度(\lambda\equiv\lambda_E). 若以西經度(\lambda_W)表示, 則 -GHA 前的負號應去除, 改為 +GHA.
** 計算天體及人造衛星的地面軌跡 (ground track).
** 繪製太陽所定義出來的晨昏線 (circle of illumination, Terminator (solar)).
** 繪製日全蝕的地面的全蝕路徑 (path of totality).
- 天體軌道測定 (orbit determination): 由移動天體或人造衛星在不同時間 t 的地面觀測數據 (a, A) (或 (\alpha, \delta)) 及雷達測距等數據, 推算天體的位置及速度向量 (狀態向量), 並透過軌道力學公式, 推測天體的 軌道要素 (), 以便預測其後續在任意時間點 t 的位置 ((\alpha, \delta) 或 (a,A)).
日出方程式
天體的出沒及持續時間,可以由前述Z_h中的高度角a及時角H(\alpha,t,\lambda)的關係推導出來。其中, 時角隱藏時間(t (jd))、觀測地經度(\lambda)、觀測天體赤經(\alpha), 是計算與時間相關的問題時,會被檢視的變數。
天體出沒高度的修正因素
天體出沒於地平線時,按理說其高度角a_0應該為0^\circ。但由於某些天體視直徑(d)很可觀,或觀測地的氣候條件(如大氣折射效應, Refraction)影響天體出沒時間,或周遭地理狀況(如海拔及障礙物)影響出沒高度,或觀察對象為較近地球的天體,地表及地心視差不容忽略,因此,天體出沒時的高度角a_0未必為零。
- 比如說,太陽的視直徑 d 約有半度角。因此,當太陽中心還在地平面下半個視直徑的時候,太陽圓盤的上沿 (superior limb) 就已經碰到地平面,日出就算發生了。要計算這種狀況,就必須令 a_0 = -{d}/{2}。其中的負號,表示還在地平面以下,出沒就發生的狀況。
- 又比如,天體在地平面之下時,由於大氣折射的關係,視位置會被拉高某個角度 R,而導致還沒浮出地平面,就已經被看到。那麼,就必須令 a_0 = -R。
- 若兩效應都考慮,就應把這兩個視角度加總。一般計算日出日落時間時,a_0 = -{d}/{2} - R = -50' 是很常用的修正條件。
- 觀測站如果位於高海拔,其效應等同降低地平面高度,會提早看到真正地平面下的天體。若該角度為 \eta_1,可令 a_0 = -\eta_1 予以補償。
- 如果有天體浮出地平面後還看不到的情況,比如被地平面上的高山或建築擋到了,而該障礙物的視角度為 \eta_2,那就該在高度角上加計這樣的角度 a_0 = +\eta_2。
- 當觀測天體為月球時,因地月平均距離D_{EM}僅約為地球半徑R_E的60倍(約直徑的30倍)。這比其它天體都要近得多。所以,由地心觀看的月球高度角與由地表觀測的高度角,會出現顯著的地心視差 (geocentric parallax) P。這個視差在月球出沒於地平線時達最大P_{max} = \arcsin(R_E/D_{EM}) \approx 1/60 \text{rad} = 0.955^\circ = 57'。亦即月亮中心的地心高度角必需昇到 57'時,地表的觀察者才會由地平線看到其中心昇上地平線。故計算月出及月沒時間時,令a_0 = +P_{max},可讓因地心視差之因素而造成的月出月沒時間延遲,得以修正。又因 57' 的視差約為月亮直徑 (32') 的兩倍。故地心人看到月昇大約2個月亮高時,地表人才看得到月亮升起。
- 有時,我們關心的並非天體實體而是它引發的效應所發生的時間。比如,曙光(太陽的光線而不是太陽實體)開始出現或暮光結束的時間。這時,可以設定一個慣用的高度角 TW (twilight angle),來計算此類事件的起始及終止時間,而令 a_0 = -TW。由於曙暮光計算所用的角度遠大於其他修正角度,因此其他修正項基本上可以忽略。
所以,計算天體(含太陽及月亮)的實體或虛像出現或隱沒時,高度角一般由以下幾個可選的參數決定其原理類似前述的截距法。但所有修正都由飛行電腦及導航機構自動完成。
在實際運用上,天體導航有時會與其他導航系統結合,截長補短,以提高導航的精確度。例如,恆星追蹤儀可能跟慣性導航系統 (INS, Inertial Navigation System) 結合,形成所謂的星光慣性導航系統 (Stellar Inertial Navigation Systems)。
太阳的位置
在地平坐標系統中,有好幾種方法可以計算太陽的視位置。
完整和精確的計算方法可以參考比利時天文學家簡米斯的天文計算(Astronomical Algorithms)
下面是一種簡單的近似計算法的例子:
已知:
- 在一年中的日期和當天的時間
- 觀測者的地理坐標(經度、緯度)和時區
以下的公式可以算出太陽的赤緯:
*\delta = -23.45^\circ \cdot \cos \left(\frac{360^\circ}{365} \cdot \left( N + 10 \right)\right)
** 此處的N是自1月1日開始的天數。
- 用以下步驟可以計算太陽時角 H。此處的真時角是觀測者因為地球的自轉與太陽之間相對應的角度。
* 令hh:mm*為觀測者由計時器所得到的時間。
* 將時與分結合成一個變數T = hh + mm*/60,單位為時。
* hh:mm*是官方(公眾)在時區中所使用的時間,與觀測者所需要的地方時(真正以太陽的位置定出的時間)是不同的,T必須依照經度來修正+(經度/15-時區),這是時區內的標準時間和觀測者在地的真太阳時之間的差異。
** 如果在夏季有使用日光節約時間(或稱夏令時),還要從官方時間減一小時,才是當地的標準時。
** 當天的均時差也要加入,由於T的單位是時,均時差需要除以60,由分轉為時之後才能併入。
** 現在已經可以算出太陽的時角了。事實上這個角度是由下面的算式直接得到:
*H =(12 - T)\cdot 15
** 由於T是以時來計算,而地球每小時轉動15度,所以H的單位是度。如果要轉成徑度量,只要成上2π/360就可以了。
- 使用赤道座標轉地平座標的公式,由 (H, \delta) 計算太陽的高度與方位:(A_S, a) 或 (A_N, a)。
** 這樣的訊息,可以用來自動追蹤太陽,讓以太陽驅動的能源裝置,如地面或太空中的太陽能光電模組,能對準太陽光來向,獲取最大的日照強度及太陽能。
天體的位置
其他移動的天體,如行星、彗星、小行星、月亮、行星的衛星、人造衛星、太空船等,也可以透過解克卜勒方程式,求得它們在其軌道面的即時位置,再透過適當的座標轉換,將軌道面座標轉換到赤道座標,再轉換到水平座標。這樣就可以預測它們的方位角及高度角,再以人工或自動的方式,加以追蹤。(太陽系行星位置計算及行星軌道要素
可參考 NASA 外部鏈結.)
參考文獻
外部鏈結
*[https://docs.kde.org/trunk5/en/kstars/kstars/ai-skycoords.html Kstars, Jason Harris, Celestial Coordinate Systems]
*[https://edu.kde.org/kstars/index.php 適合Linux/KDE和Windows/OsX的天象儀程式KStars]
*[https://promenade.imcce.fr/en/pages3/367.html Rise and Set of a Celestial Body]
*[https://ssd.jpl.nasa.gov/planets/approx_pos.html JPL, NASA, Approximate Positions of the Planets]
*[https://www.youtube.com/watch?v=fn9xMkNUMmY Celestial Navigation (+Intercept Methods)]
*[https://www.youtube.com/watch?v=G4DRBi66cOA WWII U.S. Army Air Force Training Film - Celestial Navigation]
评论 (0)