球面三角

天文學，特別是球面天文學需要球面三角學的知識。球面三角中，常要用到角度和圓弧的度量關係。從平面三角學我們知道，一圓周的1/360，叫做1度的弧。

角和弧的量度單位，常用的有兩種：

弧度：長度和半徑相等的圓弧所對的圓心角，叫做1弧度（rad）。

由於一圓周的長度等於2π個圓半徑的弧長，根據以上弧度的定義，得到弧度和度的關係如下：

2πrad=360°；

1rad= 360/2π =57.3°= 3438′= 206265′′；

或者 1°=1/57.3 rad；

1′=（1/60）°=1/3438 rad；

1′′=（1/60）′=1/206265 rad。

如果一個角的值以弧度表示時為θ，那麼以度表示時其值為57.3°×θ；以角分表示時為3438′×θ；以角秒表示時為206265′′×θ。為了方便起見，我們用符號θ°，θ′，θ′′表示一個角的度數、角分數、角秒數。

θ°=57.3°θ，θ′=3438′θ，θ′′=206265θ′′。

當角度很小時，角度的正弦或正切常可以近似地用它所對的弧來表示。

例如：sin1′′≈tan1′′≈1′′=1/206265 rad

由此得：1rad=206265′′=206265 sin1′′

根據相同的理由，得：sinθ′′≈tanθ′′≈θ′′= θ/206265=θsin1′′

上式常寫為：θ=θ′′sin1′′

球面上的圓：從立體幾何學得知，通過球心的平面截球面所得的截口是一個圓，叫做大圓；不通過球心的平面截球面所得的截口也是一個圓，叫做小圓。通過球面上不在同一直徑兩端的兩個點，能做並且只能做一個大圓。

例如通過圖1中的任意兩點A和B，也僅可以做一個大圓ABC。A、B兩點間的大圓弧（小於180°的那段弧）可以用線長、也可以用角度計量，在天文上常用角度來計量，叫做A、B間的角距，記為⌒AB（⌒應該畫在AB的上方，下同） ，它等於大圓弧⌒AB所對的中心角∠AOB。

球面上圓的極：設⌒ABC為球面上的一個任意圓（圖2），它所在的平面為MABC，又設PP’為垂直於平面MABC的球直徑，則它的兩個端點P和P’叫做圓⌒ABC的極。如果用一句話來表達，可以這樣説：垂直於球面上一已知圓（不論大圓或小圓）所在平面的球直徑的端點，叫做這個圓的極。

球面上某一圓的極和這個圓上任一點的角距，叫做極距。可以證明，極到圓上各點的角距都是相等的；如果所討論的圓是一個大圓的話，則極距為90°。

球面角：兩個大圓弧相交所成的角，叫做球面角。它們的交點叫做球面角的頂點。大圓弧本身叫做球面角的邊。圖3繪出了兩個相交的大圓弧⌒PA和⌒PB，O為球心，⌒PA所在的平面為POA，⌒PB所在的平面為POB，兩者的交線為OP。球面角∠APB用POA和POB所構成的兩面角來量度。在圖3中做以P為極的大圓⌒QQ’，設⌒PA（或其延線）和⌒QQ’相交於A’，⌒PB（或其延線）和⌒QQ’相交於B’，則由於P為⌒QQ’的極，所以OP垂直於平面QQ’，因而也垂直於OA’和OB’，所以∠A’OB’就是平面POA和POB所構成的兩面角。即：球面角∠APB可以用∠A'OB'量度，又因為∠A'OB'可以用 A'B'量度，所以最後得到的球面角∠APB是以⌒A'B'弧量度的。

從上面的討論可以概括出下述結果：如果以球面角的頂點為極作大圓，則球面角的邊或其延長線在這個大圓上所截取的那個弧段便是球面角的數值。

球面三角形：把球面上的三個點用三個大圓弧聯結起來，所圍成的圖形叫做球面三角形。這三個大圓弧叫做球面三角形的邊，通常用小寫拉丁字母a、b、c表示；這三個大圓弧所構成的角叫做球面三角形的角，通常用大寫拉丁字母A、B、C表示，並且規定：A角和a邊相對，B角和b邊相對，C角和c邊相對（如圖4所示）。三個邊和三個角合稱球面三角形的六個元素。 ^[1] 

極三角形：設球面三角形ABC各邊a、b、c的極分別為A'、B'、C'（圖5），並設弧⌒AA'、⌒BB'、⌒CC'都小於90°，則由通過A'、B'、C'的大圓弧構成的球面三角形A'B'C'叫做原球面三角形的極三角形。

極三角形和原三角形有着非常密切的關係，這種關係存在着兩條定理。

定理1：如果一球面三角形為另一球面三角形的極三角形，則另一球面三角形也為這一球面三角形的極三角形。這條定理很容易證明，請讀者自證。

定理2：極三角形的邊和原三角形的對應角互補；極三角形的角和原三角形的對應邊互補。

證明：B'是b的極（圖5），C'是c的極，所以有：

⌒B'E=⌒C'D=90°；⌒B'E+⌒C'D=180°。

即⌒B'C'+⌒DE=180°

但由定理1，A是⌒B'C'的極，故有⌒DE=A，將此式以及⌒B'C'=a'代人上式，便得到

a'+A=180° （1.1）

（1.1）式即定理2的前半的證明。定理2的後半不需證明；因為實際上，它只是定理1和定理2的前半的一個推論。

1．球面三角形兩邊之和大於第三邊。

證明：將球面三角形ABC的頂點和球心O連結起來（圖6），由立體幾何得知：三面角的兩個面角之和大於第三個面角，即∠AOB+∠BOC>；∠AOC。故c+a>b。同理a+b>c，a+c>a。

推理：球面三角形兩邊之差小於第三邊。

2．球面三角形三邊之和大於0°而小於360°。

證明：因為a，b，c均為正，故a+b+c>0°，又由立體幾何得知凸多面角各面角之和小於360°，因此∠AOB+∠BOC+∠COA<360°；0

3．球面三角形三角之和大於180 °而小於540 °。

證明：由極三角形和原三角形的關係得：a'+A = 180°， b'+B = 180°，c'+C = 180°，即 A+B+C = 540°―(a'+b'+c')。

但根據定理2有：0°

所以上式化為180°

除了上述三個基本性質以外，還有兩個重要的基本性質；對於這兩個性質，我們只寫出結果，而不給出證明。

4．若球面三角形的兩邊相等，則這兩邊的對角也相等。反之，若兩角相等，則這兩角的對邊也相等。

5．在球面三角形中，大角對大邊，大邊對大角。

下面我們要推導出六個基本公式，它們全是針對三個邊都小於90°的球面三角形導出的，但是能夠證明所得公式適用於任何球面三角形。 ^[2] 

取球面三角形ABC，將各頂點與球心°連結，可得一球心三面角O-ABC（圖7）。

過頂點A做b、c邊的切線，分別交OC，OB的延長線於N、M，由此得到兩個平面直角三角形OAM、OAN和兩個平面普通三角形△OMN、△AMN。

在平面三角形OMN中，應用平面三角的餘弦定理，得MN²=OM²+ON²－2OM·ONcosa。

同理，在平面△AMN中，得MN²=AM²+AN²－2AM·ANcosA．

因此OM²+ON²－2OM·ONcosa = AM²+AN²－2AM·ANcosA，即2OM·ONcosa = (ON²－AN²）+(OM²－AM²）+2AM·ANcosA=OA²+OA²+2AM·ANcosA。或cosa=(OA/ON)(OA/OM) +(AN/ON)(AM/OM）·cosA。

將OA/ON=cosb，OA/OM=cosc，AN/ON=sinb，AM/OM = sinc代入上式，便得到cosa=cosbcosc+sinb sinccosA （1.2）

（1.2）式是a邊的餘弦公式，其他兩個邊的餘弦公式在形式上和（1.2）式完全一樣，可以用依次輪換邊和角的字母的方法而得出，所以從（1.2）式得到b邊的餘弦公式為cosb=cosccosa+sincsinacosB （1.3）

由（1.3）式得到c邊的餘弦公式為cosc=cosacosb+sina sinbcosC （1.4）

（1.2）式、（1.3）和（1.4）式合稱邊的餘弦公式，可以用文字表達為：球面三角形任意邊的餘弦等於其他兩邊餘弦的乘積加上這兩邊的正弦及其夾角餘弦的連乘積。

設球面三角形ABC的極三角形為A'B'C'，則按照（1.2）式有cosa'=cosb' cosc'+sinb' sinc' cosA'。

因為 a'=180°－A，b'=180°－B，c'=180°－C，A' =180°－a，所以上式化為－cosA=cosBcosC－sinBsinC cosa。即 cosA=－cosBcosC+sinBsinCcosa （1.5）

利用輪換變更字母法，可以得出B角和C角的餘弦公式。（1.5）式就是角的餘弦公式，可用文字表達為：球面三角形任一角的餘弦等於其它兩角餘弦的乘積冠以負號加上這兩角的正弦及其夾邊餘弦的連乘積。

取球面三角形ABC，做球心三面角O-ABC。過C點做OAB平面的垂線交此平面於D（圖8），再從D向OA、OB引垂線DE，DF。連接CE和CF；由此得四個平面三角形OEC、OFC、CDE、CDF。因CD垂直於平面OAB，DE⊥OA，所以OA⊥CE；同理OB⊥CF，因此，四個平面三角形OEC、OFC、CDE、CDF都是直角三角形，並且有∠CED=A，∠CFD=B。

從圖8可得sina/sinA =(CF/OC)/(CD/CE)=(CF·CE)/(OC·CD)，sinb/sinB =(CE/OC)/(CD/CF)=(CF·CE)/(OC·CD)。因得 sina/sinA= sinb/sinB。利用輪換變更字母法，可以得出其它兩個類似的式子，最後得sina/sinA= sinb/sinB= sinc/sinC（1.6）

（1.6）式就是正弦公式，文字表達為：球面三角形各邊的正弦和對角的正弦成正比。

由邊的餘弦公式有：cosa =cosb cosc + sinb sinc cosA；cosb=cosc cosa + sinc sina cosB。

第二個式子可以改寫為：sinc sina cosB=cosb－cosc cosa。

將第一個式子代入上式的右邊，得sinc sina cosB=cosb－cosc（cosb cosc+sinb sinc cosA）=cosb－cosb (cosc)²－sinb sinc cosc cosA=cosb (sinc)²－sinb sinc cosC cosA。

將上式兩端各除以sinc，便得到sinacosB=cosbsinc－sinbcosccosA （1.7）

同理，得sinacosC=cosc sinb－sinc cosb cosA （1.8）

其他類似的式子可以從（1.7）式或（1.8）式，利用輪換變更字母法得出。（1.7）式或（1.8）式都是第一五元素公式，它具有一定的規律，但是它的文字表達式很繁瑣，因此這裏不寫出來。

利用極三角形和原三角形的關係（定理2），可以導出下列兩個公式：

sinA cosb=cosB sinC+sinB cosC cosa （1.9）

sinA cosc=cosC sinB+sinC cosB cosa （1.10）

其它類似的式子可以從（1.9）式或（1.10）式，利用輪換變更字母法得出，（1.9）或（1.10）式都是第二五元素公式。它的文字表達式也沒有必要寫出來。

把第一五元素公式和正弦公式聯合起來，可以導出球面三角形中相鄰的四個元素的關係式，即：

cotA sinC=－cosCcosb+sinb cota （1.11）

cotA sinB=－cosB cosc+sinc cota （1.12）

其他類似的式子，可以從（1.11）式或（1.12）式利用輪換變更字母法得出。

有一個角等於90°的球面三角形叫做直角球面三角形。設球面三角形ABC中，C=90°，且cosC=0，sinC=1，將它們代入以上各公式，經過適當的變換，可得下列常用的直角三角形公式：

cosc = cosa cosb， sinb = sinc sinB

sina = sinc sinA， sinb = tana cotA

sina = tanbcotB， cosc = cotAcotB

cosB = tana cotc， cosA=tanb cotc

cosB = cosb sinA， cosA=cosa sinB ——（1.13）

為了便於記憶這十個直角三角形公式，聶比爾提出了一條很有用的定則。除掉直角C，用（90°－a）和（90°－b）分別代替夾直角的兩個邊a和b，然後把所得的五個元素依序排成一個圓（如圖9所示）；這樣，每個元素有兩個相鄰元素和兩個相對元素。聶比爾定則為：每個元素的餘弦等於兩相鄰元素的餘切的乘積或者等於兩相對元素的正弦的乘積。例如，當所選元素為C時根據定則的前半得cosc=cotAcotB，這就是（1.13） 式裏的第六式。根據定則的後半得cosc=sin（90°－b)sin（90 °－a)=cosa cosb，這就是（1.13）式裏的第一式。 ^[3] 

全球網格模型在全球空間數據庫、全球性問題研究和基於空間信息集成等方面具有潛在的應用價值，國內外學者對之開展了大量工作。全球網格構建主要採用下面3種方式(Whiteetal.，1992，1998，1999；Sahretal.，2003)：第1種方式是在球面柏拉圖立體上直接剖分球面，如Sadourny等(1968)的剖分模型以及LianSong等(2000)所提出的小圓弧剖分模型等；第2種是循環剖分柏拉圖立體的弦，然後映射到球面，如Dutton(2000，1996，1998)所提出的QTM模型等；第3種均勻剖分展開在柏拉圖立體各面的地圖投影，如Snyder(1989)多面體投影族、Fekete(1990a，b)提出的SQT模型(該模型實質上是一種Gnomonic投影剖分(Snyder&Voxland，1989))、Fuller-Gray投影(Gray，1995，1994)等，然後反投影到球面。3種方法各有優缺點。其中Sadourny及Dutton的QTM模型頂點在球面上的分佈具有確切的映射位置，但邊界模糊，缺乏明確的數學定義，直接影響模型的實際應用。Dutton利用ZOT投影將正八面體投影到平面上，轉換為由平面直角三角組成的正方形網格，通過遞歸求解找到經緯度所對應的三角格子，座標轉換相對容易；而GoodChild等(1992)以及Otoo等(1993)基於一定的數學假設推導出QTM擴展投影方程，將經緯度座標快速轉換到平面上，座標計算簡單快速。LianSong所設計的模型通過調整小圓弧空間關係構建等面積三角剖分，計算方式複雜，座標轉換困難，計算量大。第3種剖分方式直接剖分投影三角面，單點座標轉換方便，但剖分幾何意義不明確，有的能保證等面積特性，但等角變異大，如基於Snyder等面積投影的ISEA剖分模型邊界扭曲嚴重，投影計算複雜，投影和逆投影轉換沒有顯式求解方法；有的能保證等角形變，如Gnomonic投影(即SQT剖分模型)，但面積變異太大；Fuller-Gray投影能保證等角和等面積特性相對均衡，但座標計算較複雜(投影座標反算到經緯度座標無顯式解，需要通過數字逼近方法計算得到)。Goodchild，White和Clarke等(1995)分別對全球網格剖分模型的設計和選擇標準進行了研究，指出全球網格模型評價和設計原則是一個多要素相互獨立的多指標體系，已知剖分模型中不存在一種所有評價指標最好的模型，每個模型或者一些指標性能很好，但另外一些指標相對較差。 ^[4] 

根據空間信息表達、組織和管理的應用特點，本文選擇剖分等形狀、等面積、座標轉換計算簡單快速作為全球格網主要評價指標。按該三指標體系，已有的剖分模型的指標性能不均衡，如Fuller-Gray和ISEA剖分相對等形狀和等面積，但座標轉換繁瑣，而且ISEA邊界扭曲嚴重；QTM座標轉換相對容易，但形狀和麪積的分佈存在較大差異。本文設計了等角比例投影，該投影可支持包括正八面體，正二十面體等多種柏拉圖立體，球面特徵曲線(緯度圈，大圓弧等)上的均勻剖分點均勻地投影到投影面中。該投影的數學幾何模型簡單，逆計算方便。按照該投

影性質，然後基於正二十面體構建了EARPIH剖分模型。該模型幾何意義明確，較好地滿足了近似規則、等邊、等面積的要求；剖分邊界相對平滑，扭曲不明顯；同時球面點與格網的關係轉換為投影與投影三角格網間的線性關係。與其他剖分模型相比，三指標相對均衡，是一種比較理想的剖分模型。

全球格網中只有原始球面三角能夠保證嚴格的等形狀特性，比如正八面體的球面剖分或者正二十面體球面剖分，其他任意深度的剖分將破壞等邊性，等面積性。因此只能儘量保證近似等邊性以及近似等面積性。由於全球網格模型構建方法種類繁多，包括球面直接剖分、球面弦剖分以及投影平面剖分等方式，難以直接進行投影方程的面積形變和角度形變分析。因多數剖分的邊界定義模糊，計算複雜，本文將以剖分單元的3個頂點為計算依據，採用類似Dutton的方法直接利用每個球面三角單元對應的平面三角作為標準進行QTM、ISEA、SQT和EARPIH模型的比較，當剖分越深，平面三角的弦邊與球面三角的弧邊越接近，球面三角也越接近平面三角。

在正二十面體等角比例投影基礎上，本文構建了EARPIH正二十面體球面剖分模型。該剖分模型有兩種產生方法：一種是投影平面均分的三角格網反投影到球面產生；第二種方法直接均勻剖分球面特徵弧段，結點的幾何意義明確，但邊界定義模糊，必須藉助第一種方法解決。

全球網格模型的空間幾何均勻性直接與正多面體的階數相關，階數越大，剖分越相對均勻，如基於正八面體的QTM與基於正二十面體的EARPIH，ISEA，SQT等模型相比，等面積性和等邊性存在較大差異，性能相對較差。而後3種剖分基於正二十面體，階數最高，有關指標相對接近。與ISEA相比，EARPIH的主要優勢在於計算相對簡單，投影方程和逆投影方程具有顯式解，而ISEA計算複雜，投影和逆投影轉換隻能通過數字計算方式，且ISEA的剖分邊界存在嚴重的扭曲。EARPIH的等面積性好於SQT。總之，以等邊性比較，SQT最好，EARPIH好於ISEA，QTM最差；在等面積性方面，ISEA最好，EARPIH次之，SQT第三，QTM最差；計算方法方面，QTM最簡單，SQT與EARPIH涉及三角函數，計算複雜度接近，而ISEA最為複雜。因此EARPIH基本滿足了設計標準，三個指標相對均衡。 ^[5]