常用緯度差異極值符號表達式

李忠美，李厚樸，邊少鋒，2

1.海軍工程大學 導航工程系，湖北 武漢430033；2.中國科學院 測量與地球物理研究所，湖北 武漢430077

1 引 言

大地緯度是測量和地球科學計算中最常用的一種緯度，但是在測量和地圖投影理論推導中，為滿足某種投影性質，也常會用到其他5種輔助緯度（地心緯度、歸化緯度、等距離緯度、等角緯度和等面積緯度）的概念［1－4］。隨著空間技術和計算機技術在大地測量及地圖學中的應用和發展，研究大地緯度及其他5種輔助緯度間的關系以及它們之間的差異問題具有更加重要的實用價值。對于這一問題，國內外許多學者進行了深入研究，取得了顯著成果。文獻［4］推算出了大地緯度、地心緯度、歸化緯度之間的關系以及它們的拉格朗日共軛級數展開式；文獻［5—11］人工推導出了等角緯度、等面積緯度及等距離緯度展至sin 8B展開式；文獻［12—16］借助具有強大數學分析功能的計算機代數系統［17－19］推算出等面積緯度、等角緯度及等距離緯度的偏心率e的冪級數展開式，發現和糾正了人工推導的正解公式中某些項的偏差，推導出的系數具有更高的精確度。

從目前來看，前人對這一領域做了很多卓有成效的工作，但是卻鮮有文獻將這幾種常用緯度進行系統比較。為豐富對這一問題的研究，使人們對這幾種常用緯度形成較直觀的認識，本文著重研究了常用緯度間的差異極值問題，推導出了常用緯度間差異極值點及對應差異極值的符號表達式，最后以CGCS2000參考橢球［20］為例，對常用緯度間的差異進行了數值分析和對比。

2 常用緯度定義

根據文獻［3］可知，歸化緯度u、地心緯度φ是大地測量學中常用緯度，記橢球偏心率為e，則它們關于大地緯度B的表達式為

又由文獻［13］可知，等距離緯度ψ、等角緯度φ和等面積緯度?是地圖學中常用的3種緯度，它們與大地緯度B的關系式為

式中的系數如下

3 常用緯度間的差異極值符號表達式

由于6種常用緯度相差微小，實際應用中經常會用到它們的差值表達式。為了解常用緯度間的差異情況，特別是差異極值點及其對應差異極值，可基于它們之間的差值表達式推導出常用緯度間的差異極值符號表達式。

3.1 輔助緯度與大地緯度差異極值符號表達式

通過對輔助緯度與大地緯度的差值表達式進行求導，即可推算出輔助緯度與大地緯度差異極值點的解析表達式。以歸化緯度為例

略去推導過程，可得歸化緯度與大地緯度差異一階導數為零的點

將式（5）代入歸化緯度與大地緯度差值表達式的二階導數進行檢驗，可知二階導數不為零，即式（5）所示的大地緯度即為所求的差異極值點，將它代入歸化緯度與大地緯度的差值表達式，即可推算出差異極值

為方便不同的輔助緯度與大地緯度差異極值間的比較，可進一步將式（5）、式（6）展開成偏心率e的冪級數形式

類似的，對地心緯度與大地緯度的差值表達式進行求導，如式（8）所示

略去推導過程，可得地心緯度與大地緯度差異一階導數為零時對應的大地緯度

類似的，將式（9）、式（10）展開為關于偏心率e的冪級數形式

經檢驗，式（9）所示的大地緯度為極值點，將其代入地心緯度與大地緯度差值表達式，即可推算出地心緯度與大地緯度差異的極值

等距離緯度與大地緯度的關系式為關于偏心率e冪級數形式，可通過對它們的差值表達式求導，并利用迭代法推導出等距離緯度與大地緯度的差異極值點及對應的差異極值

即

將cos 2B項作為等式左邊項，其他項移到等式右邊，并對等式兩端同時取反三角函數

由式（13）可以看出，由于地球橢球偏心率很小，在B∈［0，π／2］范圍內，若要滿足緯度差值表達式的導數為零，即等式左邊首項近似為零，也就是差異極值點在B＝π／4附近。

將初值B0＝π／4代入式（14）進行迭代，經過1次迭代后的結果為

將B1代入式（14）繼續進行第2次迭代，依此類推，可以發現第4次迭代結果與第3次迭代結果相比，它們擴展至e10的展開式中各項系數已不再發生變化，故迭代終止。略去推導過程，得差異極值點的符號表達式，如式（16）所示

將式（16）代入ψ－B表達式中，并對結果進行冪級數展開，可推導出等距離緯度與大地緯度差異的極值，如式（17）所示

由于等角緯度、等面積緯度與大地緯度的關系式形式復雜，求導后不便于解算，又由文獻［13］知，它們擴展至sin 10B的級數展開式也已具有相當高的精度，等角緯度和等面積緯度的正解展開式的精度均優于10－8″，故可采用與等距離緯度類似的方法，推導出它們與大地緯度的差異極值點及對應的差異極值符號表達式，最后結果列于表1。

表1 輔助緯度與大地緯度差異極值符號表達式Tab.1 The symbolic expressions of difference extrema between auxiliary and geodetic latitudes

由表1可以看出，各輔助緯度與大地緯度的差異極值點均在B＝π／4右側，且略有不同，各輔助緯度與大地緯度差異極值點按從小到大的順序依次為：歸化緯度、等面積緯度、等距離緯度、等角緯度和地心緯度。輔助緯度中與大地緯度差異極值的絕對值較大的是地心緯度以及等角緯度，且它們與大地緯度的差異極值十分接近；等距離緯度次之，而歸化緯度與大地緯度的差異極值的絕對值最小。

3.2 輔助緯度之間差異極值符號表達式

為系統地比較各常用緯度之間的差異，除對輔助緯度與大地緯度間的差異進行分析以外，需對輔助緯度之間的差異進行分析。以下分別對大地測量中常用的地心緯度與歸化緯度間的差異，地圖學中常用的等距離緯度、等角緯度以及等面積緯度之間的差異進行分析。與推導輔助緯度和大地緯度間差異極值點類似，通過對輔助緯度之間的差值表達式進行求導（以地心緯度與歸化緯度間為例）

可得

經檢驗，式（19）中所示大地緯度為地心緯度與歸化緯度差異極值點，將它代入地心緯度與歸化緯度的差值表達式，推算出兩者的差異極值

為方便不同極值間的比較，將式（19）、式（20）展開為關于偏心率e的冪級數形式，可得差異極值點及對應極值的級數展開式

與地心緯度和歸化緯度不同，地圖投影理論中常用的3種緯度關于大地緯度的函數形式復雜，導致它們的差值表達式不便于求導，因此可借助計算機代數系統，將它們間的差值表達式進行級數展開，并按正弦函數的倍角形式將展開結果進行整理，略去推導過程，可得（以等面積緯度與等角緯度的差異極值點為例）

式中系數為

對式（22）進行求導并采用迭代法，即可推算出等面積緯度與等角緯度之間的差異極值點。類似的，可推導出等角緯度、等面積緯度分別與等距離緯度間的差異極值點及對應極值，最后結果列于表2。

表2 等角緯度、等面積緯度及等距離緯度間差異極值符號表達式Tab.2 The symbolic expressions of difference extrema between conformal，authalic and equidistant latitudes

由表2可以看出，地圖學中常用的3種緯度間的差異極值點均在B＝π／4右側，且存在微小差異。結合表1可知，它們之間差異極值的絕對值均小于它們與大地緯度差異極值的絕對值，其中，等距離緯度與等面積緯度間差異極值的絕對值最小，而等面積緯度與等角緯度差異極值的絕對值最大。

4 算例分析

由于地球橢球偏心率很小，且不同參考橢球的偏心率非常接近。為使人們對各緯度的差異在數值上有一個直觀的認識，下面以CGCS2000參考橢球（e＝0.081 819 191 042 8）為例，對6種常用緯度間的差異進行數值比較與分析。

4.1 輔助緯度與大地緯度的差異極值

為明確了解各輔助緯度與大地緯度的差異情況，可繪制出大地緯度B∈［0°，90°］范圍內，各輔助緯度與大地緯度的差異曲線圖，如圖1所示。表3為大地緯度B∈［0°，90°］時，每隔15°對應的各輔助緯度與大地緯度間的差異。

圖1 輔助緯度與大地緯度的差異曲線圖Fig.1 Chart of differences between auxiliary and geodetic latitudes

表3 輔助緯度與大地緯度的差異Tab.3 Values of differences between auxiliary and geodetic latitudes

由圖1及表3可以看出，各輔助緯度與大地緯度差異的絕對值先變大后變小，且在大地緯度45°附近處，各輔助緯度與大地緯度的差異出現極值。

根據已推導出的各輔助緯度與大地緯度的差異極值點及對應的差異極值的符號表達式，可以計算出B∈［0°，90°］時，各輔助緯度與大地緯度間的差異極值，如表4所示。

表4 輔助緯度與大地緯度差異極值Tab.4 Difference extrema between auxiliary and geodetic latitudes

由圖1及表4可以看出，輔助緯度與大地緯度的差異極值點均在B＝45°右側，在B＝45°05′46″處，地心緯度與大地緯度的差異極值最大，差異絕對值最大達到11′32.73″。等距離緯度與大地緯度差異極值的絕對值次之，為8′39.55″，而歸化緯度與大地緯度的差異最小，它們的差異絕對最大值為5′46.36″。等角緯度和地心緯度這兩種緯度與大地緯度的差異十分接近，它們與大地緯度的差異曲線近乎重合。

4.2 輔助緯度間的差異極值

為全面分析在CGCS2000參考橢球中常用緯度間的差異極值問題，可繪制出大地緯度B∈［0°，90°］時大地測量中常用的地心緯度與歸化緯度間差異圖，及地圖學中常用等角緯度、等面積緯度及等距離緯度這3種緯度間的差異圖，如圖2、圖3所示。

圖2 地心緯度與歸化緯度差異圖Fig.2 Chart of difference between geocentric and reduced latitudes

圖3 等角緯度、等面積緯度、等距離緯度間的差異Fig.3 Chart of differences between conformal，equivalent and equidistant latitudes

表5為大地緯度B∈［0°，90°］時，每隔15°對應的輔助緯度間的差異。

表5 輔助緯度間的差異Tab.5 Values of differences between auxiliary latitudes

由圖2、圖3及表5可以看出，在B∈［0°，90°］范圍內，地心緯度與歸化緯度間差異的絕對值和等角緯度、等面積緯度與等距離緯度間差異的絕對值都是先變大后變小，且這4種差異極值均出現在大地緯度45°附近。根據3.2中已推導出的極值符號表達式，可計算出輔助緯度間差異極值點及對應差異極值，如表6所示。

表6 輔助緯度間差異極值Tab.6 Difference extrema between auxiliary latitudes

由圖2、圖3及表6可以看出，地心緯度與歸化緯度的差異極值點為B＝45°08′39″，對應差異極值為－5′46.36″；地圖學中常用的3種緯度間，等距離緯度與等角緯度間差異極值最大，為3′50.47″，而等距離緯度與等面積緯度間差異最小，它們差異的極值為－57.68″。

5 結 論

本文對6種常用緯度間的差異進行了研究，推導出常用緯度間的差異極值點及對應極值的符號表達式，并以CGCS2000橢球為例，對它們進行了數值分析和對比，得出如下結論。

（1）輔助緯度與大地緯度差異的絕對值隨大地緯度的增大先變大后變小，且在B＝π／4右側，差異出現極值。其中，等角緯度與大地緯度差異極值的絕對值最大，極值符號表達式的首項為－e2／2；歸化緯度與大地緯度差異極值絕對值最小，極值符號表達式的首項為－e2／4。

（2）地圖學中常用的3種緯度間，等面積緯度與等角緯度差異極值最大，極值符號表達式首項為e2／6；等距離緯度與等面積緯度差異極值最小，極值符號表達式首項為－e2／24。

（3）將常用緯度間差異極值點及對應極值表示為符號形式，并統一展開為偏心率e的冪級數形式，該表達式易于比較分析，一定程度上豐富了測量及地圖學數學分析理論。

［1］ LüZhiping，ZHANG Jianjun，QIAO Shubo.Basis of Geodesy［M］.Beijing：Publishing House of Liberation Army，2005.（呂志平，張建軍，喬書波.大地測量學基礎［M］.北京：解放軍出版社，2005.）

［2］ HUA Tang.The Mathematical Basis of Chart［M］.Beijing：Haichao Press，1985.（華棠.海圖數學基礎［M］.北京：海潮出版社，1985.）

［3］ WANG Rui，LI Houpu.Symbolic Expressions of the Track Calculation Methods Based on the Ellipsoidal Model［J］.Acta Geodaetica et Cartographica Sinica，2010，39（2）：151－155.（王瑞，李厚樸.基于地球橢球模型的符號形式的航跡計算法［J］.測繪學報，2010，39（2）：151－155.）

［4］ XIONG Jie.Ellipsoidal Geodesy［M］.Beijing：Publishing House of Liberation Army，1988：18－23.（熊介.橢球大地測量學［M］.北京：解放軍出版社，1988：18－23.）

［5］ ADAMS O S.Latitude Developments Connected with Geodesy and Cartography with Tables，Including a Table for Lambert Equal－Area Meridional Projection［M］.Washington D C：Government Printing Office，1921.

［6］ SNYDER J P.Map Projections：A Working Manual［M］.Washington D C：Government Printing Office，1987.

［7］ YANG Qihe.The Theory and Method of Map Projection［M］.Beijing：PLA Press，1989.（楊啟和.地圖投影變換原理與方法［M］.北京：解放軍出版社，1989.）

［8］ YANG Q H，SNYDER J P，TOBLER W R.Map Projection Transformation： Principles and Application ［M ］.London：Taylor ＆Francis，2000.

［9］ YANG Qihe，Yang Xiaomei.Three Kinds of Latitudes in Surveying and Mapping and the Liner Interpolation Method for Inverse Solution［J］.Acta Geodaetica et Cartographica Sinica，1997，26（1）：92－95.（楊啟和，楊曉梅.測量和地圖學中應用的三種緯度函數及其反解變換的線性插值方法［J］.測繪學報，1997，26（1）：92－95.）

［10］ WU Zhongxing，YANG Qihe.Principle of Mathematical Cartography［M］.Beijing：Surveying and Mapping Press，1989.（吳忠性，楊啟和.數學制圖學原理［M］.北京：測繪出版社，1989.）

［11］ SUN Qun，YANG Qihe.The Research on the Computation of the Foot－point Latitude and the Inverse Solution of Isometric Latitude and Area Function［J］.Journal of PLA of Surveying and Mapping，1985（2）：64－75.（孫群，楊啟和.底點緯度解算以及等量緯度和面積函數反解問題的探討［J］.解放軍測繪學院學報，1985（2）：64－75.）

［12］ BIAN Shaofeng，JI Bing.The Expansions of Rectifying Latitude，Conformal Latitude and Authalic Latitude［J］.Acta Geodaetica et Cartographica Sinica，2007，36（2）：218－223.（邊少鋒，紀兵.等距離緯度等量緯度和等面積緯度展開式［J］.測繪學報，2007，36（2）：218－223.）

［13］ BIAN S F，LI H P.Mathematical Analysis in Cartography by Means of Computer Algebra System［C］∥Carlos Bateira.Cartography：A Tool for Spatial Analysis.Croatia：InTech，2012.

［14］ LI Houpu.The Research on the Precise Calculation Theory of Geodetic Coordinate System and Its Application Based on Computer Algebra System［D］.Wuhan：Naval University of Engineering，2010.（李厚樸.基于計算機代數系統的大地坐標系精密計算理論及其應用研究［D］.武漢：海軍工程大學，2010.）

［15］ WANG Rui，LI Houpu.The Derivation of the Inverse Expansions for Auxiliary Latitudes by Lagrange Series Method［J］.Hydrographic Surveying and Charting，2008，28（3）：18－23.（王 瑞，李 厚 樸.輔 助 緯 度 反 解 公 式的Lagrange級數法推演［J］.海洋測繪，2008，28（3）：18－23.）

［16］ LI Houpu，BIAN Shaofeng.Derivation of Inverse Expansions for Auxiliary Latitudes by Hermite Interpolation Method［J］. Geomatics and Information Science of Wuhan University，2008，33（6）：623－626.（李厚樸，邊少鋒.輔助緯度反解公式的Hermite插值法新解［J］.武漢大學學報：信息科學版，2008，33（6）：623－626.）

［17］ BIAN Shaofeng，XU Jiangning.Computer Algebra System and Mathematical Analysis in Geodesy ［M］. Beijing：National Defense Industry Press，2004.（邊少鋒，許江寧.計算機代數系統與大地測量數學分析［M］.北京：國防工業出版社，2004.）

［18］ BIAN S F，CHEN Y B.Solving and Inverse Problem of a Meridian Arc in Terms of Computer Algebra System［J］.Journal of Surveying Engineering，2006，132（2）：1－4.

［19］ BIAN Shaofeng，CHEN Liangyou，JI Bing.Application of Mathematica in Geodesy［J］.Geomatics Science and Engineering，2004，24（2）：25－28.（邊少鋒，陳良友，紀兵.Mathematica軟件在大地測量中的應用［J］.測繪科學與工程，2004，24（2）：25－28.）

［20］ CHEN Junyong.Chinese Modern Geodetic Datum：Chinese Geodetic Coordinate System 2000（CGCS2000）and Its Frame［J］.Acta Geodaetica et Cartographica Sinica，2008，37（3）：269－271.（陳俊勇.中國現代大地基準：中國大地坐標系統2000（CGCS2000）及其框架［J］.測繪學報，2008，37（3）：269－271.）