李厚樸，邊少鋒

1.海軍工程大學(xué)導(dǎo)航工程系，湖北武漢430033；2.海島（礁）測繪技術(shù)國家測繪地理信息局重點實驗室，山東 青島266510；3.中國科學(xué)院測量與地球物理研究所，湖北武漢430077

不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換

李厚樸1，2，邊少鋒1，3

1.海軍工程大學(xué)導(dǎo)航工程系，湖北武漢430033；2.海島（礁）測繪技術(shù)國家測繪地理信息局重點實驗室，山東 青島266510；3.中國科學(xué)院測量與地球物理研究所，湖北武漢430077

為避免不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間傳統(tǒng)間接變換繁瑣的計算過程，利用子午線弧長、等量緯度和等面積緯度函數(shù)間變換的直接展開式，建立了相應(yīng)投影坐標(biāo)間的直接變換模型，無需計算大地緯度即可完成變換。本文導(dǎo)出公式均為含參考橢球第一偏心率的符號形式，可解決兩類投影在不同參考橢球下的變換問題。算例分析表明與傳統(tǒng)間接變換模型相比，本文建立的直接變換模型提高了計算效率和計算精度，可供實際使用。

地圖投影；正軸圓柱投影；正軸圓錐投影；不同變形性質(zhì)；直接變換

1 引 言

正軸圓柱投影和正軸圓錐投影是在測量和地圖制圖中應(yīng)用非常廣泛的兩類重要投影，航空圖、海圖、世界全圖和區(qū)域地圖等的繪制主要采用這兩類投影方式［1－2］，在實際生產(chǎn)中經(jīng)常會遇到這兩類投影間的變換問題。文獻［3—6］對這一問題進行了深入研究，取得了卓有成效的成果，但對于橢球情形下不同變形性質(zhì)投影間的變換，則是通過間接變換法實現(xiàn)的，計算思路可簡單概括為“圓柱→橢球→圓錐”，由于需要反解出大地緯度，計算過程相當(dāng)復(fù)雜，并且有的公式表現(xiàn)為具體的數(shù)值形式，僅能解決某一特定參考橢球下的變換問題，不便于推廣使用。

子午線弧長、等量緯度和等面積緯度函數(shù)是不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影中的3類重要變量。這3類變量都是大地緯度的函數(shù)，它們和大地緯度間的正反算問題已得到國內(nèi)外學(xué)者的深入研究并已有效解決［7－14］，但對于它們之間的變換，傳統(tǒng)上則是通過解算大地緯度間接實現(xiàn)的，不僅計算非常繁瑣，而且不便于理論分析。筆者借助具有強大符號運算功能的計算機代數(shù)系統(tǒng)Mathematica［15］，對這一問題進行了新的研究，推導(dǎo)出它們之間變換的直接展開式［16－18］，為實現(xiàn)不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換提供了新的思路。本文利用文獻［16］導(dǎo)出的子午線弧長、等量緯度和等面積緯度函數(shù)間變換的直接展開式，深入研究了不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換問題，建立了相應(yīng)投影坐標(biāo)間的直接變換模型。

2 子午線弧長、等量緯度和等面積緯度

函數(shù)間變換的直接展開式［16］

2.1 子午線弧長和等量緯度間變換的直接展開式

子午線弧長X變換至等量緯度q的直接展開式為

式中，R＝a（1－e2）K0，a為參考橢球長半軸，e為參考橢球第一偏心率，系數(shù)為

等量緯度q變換至子午線弧長X的直接展開式為

式中，系數(shù)為

2.2 子午線弧長和等面積緯度函數(shù)間變換的直接展開式

子午線弧長X變換至等面積緯度函數(shù)F的直接展開式為

式中，系數(shù)為

等面積緯度函數(shù)F變換至子午線弧長X的直接展開式為

式中，R′2＝a2（1－e2）A，系數(shù)為

2.3 等量緯度和等面積緯度函數(shù)間變換的直接展開式

等量緯度q變換至等面積緯度函數(shù)F的直接展開式為

式中，系數(shù)為

等面積緯度函數(shù)F變換至等量緯度q的直接展開式為

式中，系數(shù)為

3 不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換

3.1 不同變形性質(zhì)的正軸圓柱投影和正軸圓錐投影

3.1.1 不同變形性質(zhì)的正軸圓柱投影

等距離正軸圓柱投影的坐標(biāo)公式為

等角正軸圓柱投影的坐標(biāo)公式為

等面積正軸圓柱投影的坐標(biāo)公式為

式（13）～式（15）中，l為經(jīng)差；r0i（i＝1，2，3）為基準(zhǔn)緯度處的圓柱半徑。

3.1.2 不同變形性質(zhì)的正軸圓錐投影

等距離正軸圓錐投影的坐標(biāo)公式為

式中

等角正軸圓錐投影的坐標(biāo)公式為

式中

式中，exp（－a2q）＝e－a2q，e＝2.718 281 828為自然對數(shù)的底。

等面積正軸圓錐投影的坐標(biāo)公式為

式中

式（16）～式（21）中，ρsi、αi、Ci（i＝1，2，3）為相應(yīng)投影的投影常數(shù)，計算公式見文獻［3］。

3.2 不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換模型

3.2.1 等距離正軸圓柱投影和等角正軸圓錐投影間的直接變換模型

由式（13）可反解得X、l

將式（22）代入式（18），顧及式（19），可得等距離正軸圓柱投影直接變換至等角正軸圓錐投影的坐標(biāo)關(guān)系式為

式中，q可根據(jù)式（1）表示為

由式（18）可反解得ρ2、δ2，將其代入式（19）可反解得q、l

將式（25）代入式（13），并顧及式（3），可得等角正軸圓錐投影直接變換至等距離正軸圓柱投影的坐標(biāo)關(guān)系式為

3.2.2 等距離正軸圓柱投影和等面積正軸圓錐投影間的直接變換模型

將式（22）代入式（20），并顧及式（21），可得等距離正軸圓柱投影直接變換至等面積正軸圓錐投影的坐標(biāo)關(guān)系式為

式中，F(xiàn)可根據(jù)式（5）表示為

由式（20）可反解得ρ3、δ3，將其代入式（21）可反解得F、l

將式（29）代入式（13），并顧及式（7），可得等面積正軸圓錐投影直接變換至等距離正軸圓柱投影的坐標(biāo)關(guān)系式為

3.2.3 等角正軸圓柱投影和等距離正軸圓錐投影間的直接變換模型

由式（14）可反解得q、l

將式（31）代入式（16），并顧及式（17），可得等角正軸圓柱投影直接變換至等距離正軸圓錐投影的坐標(biāo)關(guān)系式為

式中，X可根據(jù)式（3）表示為

由式（16）可反解得ρ1、δ1，將其代入式（17）可反解得X、l

將式（34）代入式（14），并顧及式（1），可得等距離正軸圓錐投影直接變換至等角正軸圓柱投影的坐標(biāo)關(guān)系式為

3.2.4 等角正軸圓柱投影和等面積正軸圓錐投影間的直接變換模型

將式（31）代入式（20），并顧及式（21），可得等角正軸圓柱投影直接變換至等面積正軸圓錐投影的坐標(biāo)關(guān)系式為

式中，F(xiàn)可根據(jù)式（9）表示為

將式（29）代入式（14），并顧及式（11），可得等面積正軸圓錐投影直接變換至等角正軸圓柱投影的坐標(biāo)關(guān)系式為

3.2.5 等面積正軸圓柱投影和等距離正軸圓錐投影間的直接變換模型

由式（15）可反解得

將式（39）代入式（16），并顧及式（17），可得等面積正軸圓柱投影直接變換至等距離正軸圓錐投影的坐標(biāo)關(guān)系式為

式中，X可根據(jù)式（7）表示為

將式（34）代入式（15），顧及式（5），可得等距離正軸圓錐投影直接變換至等面積正軸圓柱投影的坐標(biāo)關(guān)系式為

3.2.6 等面積正軸圓柱投影和等角正軸圓錐投影間的直接變換模型

將式（39）代入式（18），并顧及式（19），可得等面積正軸圓柱投影直接變換至等角正軸圓錐投影的坐標(biāo)關(guān)系式為

式中，q可根據(jù)式（11）表示為

將式（25）代入式（18），并顧及式（9），可得等角正軸圓錐投影直接變換至等面積正軸圓柱投影的坐標(biāo)關(guān)系式為

4 算例分析

為驗證本文建立的不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換模型的有效性與快捷性，同時為了與文獻［3］給出的傳統(tǒng)間接變換模型進行比較，以等角正軸圓錐投影變換至等面積正軸圓柱投影為例，選用CGCS2000橢球常數(shù)［19－20］a＝6 378 137m、1／f＝298.257 222 101進行了計算分析，具體分析思路如下。

等角正軸圓錐投影的投影區(qū)域為18°N≤B≤54°N，雙標(biāo)準(zhǔn)緯線分別取為B1＝27°N，B2＝45°N，經(jīng)差變化范圍為0°≤l≤30°，等面積正軸圓柱投影的基準(zhǔn)緯度取為0°。取定大地緯度B和經(jīng)差l，分別代入式（15）和式（18）可得等面積正軸圓柱投影坐標(biāo)（x3，y3）和等角正軸圓錐投影坐標(biāo)（X2，Y2），將（X2，Y2）分別代入文獻［3］給出的間接變換模型和本文建立的直接變換模型式（45）可得變換后的等面積正軸圓柱投影坐標(biāo)，依次記為（x′3，y′3）和（x″3，y″3），分別與（x3，y3）相減，可得間接變換模型的計算誤差（Δx′3，Δy′3）和直接變換模型的計算誤差（Δx″3，Δy″3）。

投影區(qū)域包含的點數(shù)可由大地緯度和經(jīng)差方向的分辨率確定，取分辨率為1°×1°、0.5°× 0.5°、0.1°×0.1°3種情況，間接變換模型和直接變換模型的計算用時分別記為t1、t2，如表1所示，分辨率為0.1°×0.1°時兩種模型的計算誤差統(tǒng)計情況如表2所示。

表1 間接變換模型和直接變換模型的計算用時Tab.1 The calculation time of the indirect and direct transformation models s

表2 分辨率為0.1°×0.1°時間接變換模型和直接變換模型的計算誤差Tab.2 Statistics of the calculation errors of the indirect and direct transformation models at the resolution of 0.1°×0.1° m

由表1可以看出，當(dāng)分辨率為1°×1°（共37× 31個點）時，間接變換模型用時為4.391s，直接變換模型用時僅為0.469s，約為間接變換模型用時的10.7%；當(dāng)分辨率為0.1°×0.1°（共361× 301個點）時，間接變換模型用時為275.391s，直接變換模型用時僅為45.359s，約為間接變換模型用時的16.5%。由表2可以看出，傳統(tǒng)間接變換模型的橫坐標(biāo)計算誤差和直接變換模型的橫坐標(biāo)計算誤差基本一致，均在10－9m量級；本文建立的直接變換模型的縱坐標(biāo)計算精度相比于傳統(tǒng)間接變換模型提高了4個數(shù)量級。

5 結(jié) 論

本文研究了不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的變換問題，主要結(jié)論如下：

（1）為避免傳統(tǒng)間接變換“圓柱→橢球→圓錐”繁瑣的計算過程，利用子午線弧長、等量緯度和等面積緯度函數(shù)間變換的直接展開式，建立了不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換模型，無需反解大地緯度即可實現(xiàn)投影變換。

（2）本文建立的直接變換模型為含參考橢球第一偏心率的符形式，可解決兩類投影在不同參考橢球下的變換問題，適合于計算機編程計算。

（3）以等角正軸圓錐投影變換至等面積正軸圓柱投影為例，對本文建立的直接變換模型的有效性和快捷性進行了分析。結(jié)果表明，與傳統(tǒng)間接變換模型相比，本文建立的直接變換模型提高了計算效率和計算精度。

［1］ LI Guozao，YANG Qihe，HU Dingquan.Map Projections［M］.Beijing：PLA Press，1993.（李國藻，楊啟和，胡定荃.地圖投影［M］.北京：解放軍出版社，1993.）

［2］ SUN Da，PU Yingxia.Map Projections［M］.Nanjing：Nangjing University Press，2005.（孫達(dá)，蒲英霞.地圖投影［M］.南京：南京大學(xué)出版社，2005.）

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

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

［5］ SNYDER J P.Map Projections－a Working Manual［M］.Washington：U.S.Government Printing Office，1987.

［6］ YANG Xiaomei，YANG Qihe.A Reseach on the Transformation Models between Coordinates of Cylindrical and Conic Projections［J］.Bulletin of Surveying and Mapping，1999（2）：5－8.（楊曉梅，楊啟和.圓錐投影和圓柱投影坐標(biāo)變換模型研究［J］.測繪通報，1999（2）：5－8.）

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

［8］ 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 Institute of Surveying and Mapping，1985（2）：64－75.（孫群，楊啟和.底點緯度解算以及等量緯度和面積函數(shù)反解問題的探討［J］.解放軍測繪學(xué)院學(xué)報，1985（2）：64－75.）

［9］ YANG Qihe，YANG Xiaomei.Three Kinds of Latitude Functions in Surveying and Mapping and Their Inverse Solution Using the Linear Interpolation Method［J］.Acta Geodaetica et Cartographica Sinica，1997，26（1）：92－93.（楊啟和，楊曉梅.測量和地圖學(xué)中應(yīng)用的三種緯度函數(shù)及其反解變換的線性插值方法［J］.測繪學(xué)報，1997，26（1）：92－93.）

［10］ CRAIG R.Auxiliary Latitude Formulas：Finding the Coefficients Numerically and Symbolically［C］∥Proceedings of Wolfram Technology Conference.America：Champaign，2006.

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

［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ì)兵.等距離緯度等量緯度和等面積緯度展開式［J］.測繪學(xué)報，2007，36（2）：218－223.）

［13］ 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］.武漢大學(xué)學(xué)報：信息科學(xué)版，2008，33（6）：623－626.）

［14］ 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級數(shù)法推演［J］.海洋測繪，2008，28（3）：18－23.）

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

［16］ LI Houpu.The Research of the Precise Computation Theory and Its Application Based on Computer Algebra System for Geodetic Coordinate System［D］.Wuhan：Naval University of Engineering，2010.（李厚樸.基于計算機代數(shù)系統(tǒng)的大地坐標(biāo)系精密計算理論及其應(yīng)用研究［D］.武漢：海軍工程大學(xué)，2010.）

［17］ LI Houpu，LIU Min，KONG Haiying，et al.The Direct Expansions for Transformations between Meridian Arc and Authalic Latitude Function［J］.Hydrographic Surveying and Charting，2011，31（1）：17－19.（李厚樸，劉敏，孔海英，等.子午線弧長和等面積緯度函數(shù)變換的直接展開式［J］.海洋測繪，2011，31（1）：17－19.）

［18］ LI Houpu，BIAN Shaofeng，CHEN Liangyou.The Direct Calculating Formulae for Transformations between Authalic Latitude Function and Isometic Latitude［J］.Geomatics and Information Science of Wuhan University，2011，36（7）：843－846.（李厚樸，邊少鋒，陳良友.等面積緯度函數(shù)和等量緯度變換的直接解算公式［J］.武漢大學(xué)學(xué)報：信息科學(xué)版，2011，36（7）：843－846.）

［19］ 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.（陳俊勇.中國現(xiàn)代大地基準(zhǔn)：中國大地坐標(biāo)系統(tǒng)2000（CGCS2000）及其框架［J］.測繪學(xué)報，2008，37（3）：269－271.）

［20］ CHENG Pengfei，WEN Hanjiang，CHENG Yingyan，et al.Parameters of the CGCS2000Ellipsoid and Comparisons with GRS80and WGS84［J］.Acta Geodaetica et Cartographica Sinica，2009，38（6）：189－194.（程鵬飛，文漢江，成英燕，等.2000國家大地坐標(biāo)系與GRS80和WGS84的比較［J］.測繪學(xué)報，2009，38（6）：189－194.）

The Direct Transformations between Normal Cylindrical and Conic Projections with Different Distortion Properties

LI Houpu1，2，BIAN Shaofeng1，3
1.Department of Navigation，Naval University of Engineering，Wuhan430033，China；2.Key Laboratory of Surveying and Mapping Technology on Island and Reef，National Administration of Surveying，Mapping and Geoinformation，Qingdao 266510，China；3.Institute of Geodesy and Geophysics，Chinese Academy of Sciences，Wuhan 430077，China

In order to avoid the fussy processes of traditional indirect transformations between normal cylindrical and normal conic projections with different distortion properties，the direct transformation models between the corresponding projection coordinates were established using the direct expansions of transformations between meridian arc，isometric latitude and authalic latitude function.The formulas are symbolical expressions which include the first eccentricity of the referenced ellipsoid and could solve the transformation problems when different reference ellipsoids are used.Numerical examples show that the direct transformation modes have improved the computation efficiency and accuracy compared to the traditional indirect transformation ones，and they could satisfy practical use.

map projection；normal cylindrical projection；normal conic projection；different distortion properties；direct transformation

LI Houpu（1985—），male，PhD，lecturer，majors in geodesy and satellite navigation.

LI Houpu，BIAN Shaofeng.The Direct Transformations between Normal Cylindrical and Conic Projections with Different Distortion Properties［J］.Acta Geodaetica et Cartographica Sinica，2012，41（4）：536－542.（李厚樸，邊少鋒.不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換［J］.測繪學(xué)報，2012，41（4）：536－542.）

P282.1

A

1001－1595（2012）04－0536－07

國家973計劃（2012CB719902）；國家自然科學(xué)基金（41071295；40904018）；海島（礁）測繪技術(shù)國家測繪地理信息局重點實驗室資助項目（2010B04）；水下測控技術(shù)重點實驗室延伸性發(fā)展基金（YS0C261102）

宋啟凡）

2011－08－29

2011－12－21

李厚樸（1985—），男，博士，講師，主要研究方向為大地測量和衛(wèi)星導(dǎo)航。

E－mail：lihoupu1985＠126.com