地圖投影計算機代數(shù)分析研究進展

邊少鋒，李厚樸，李忠美

海軍工程大學(xué)導(dǎo)航工程系，湖北 武漢 430033

地圖投影計算機代數(shù)分析研究進展

邊少鋒，李厚樸，李忠美

海軍工程大學(xué)導(dǎo)航工程系，湖北 武漢 430033

地圖投影是現(xiàn)代地圖學(xué)的重要組成部分，涉及大量的橢圓函數(shù)冪級數(shù)展開、隱函數(shù)復(fù)合函數(shù)微分、橢圓積分、復(fù)變函數(shù)運算等一系列煩瑣的數(shù)學(xué)分析過程，人工推導(dǎo)不但費時費力，而且容易出錯，有時由于難以忍受的復(fù)雜性等各種原因，甚至根本無法實現(xiàn)。本文主要從橢球各緯度間正反解符號表達式、不同變形性質(zhì)地圖投影間的直接變換、高斯投影的復(fù)變函數(shù)表示、斜軸墨卡托投影數(shù)學(xué)分析、極區(qū)海圖投影及變換等5個方面，論述了地圖投影計算機代數(shù)分析取得的研究進展，討論了該領(lǐng)域有待進一步解決的主要問題，對推動地圖投影學(xué)的發(fā)展具有積極意義。

地圖投影；投影變換；高斯投影；極區(qū)投影；計算機代數(shù)

地圖作為信息的載體，是描述地球和研究復(fù)雜地理現(xiàn)象的最有效、最直觀的工具，它以科學(xué)的符號系統(tǒng)、地圖投影和綜合方法表達復(fù)雜地理世界的空間結(jié)構(gòu)和空間關(guān)系，是人類文明史上的偉大創(chuàng)想[1-3]。地球曲面和地圖平面之間的矛盾構(gòu)成了地圖最基本的矛盾，解決這一矛盾的數(shù)學(xué)法則構(gòu)成了地圖的數(shù)學(xué)基礎(chǔ)，這是地圖最基本的特性之一，是地圖科學(xué)性和精確性的重要體現(xiàn)。地圖所采用的特殊數(shù)學(xué)法則便是地圖投影。

地圖投影學(xué)也稱為數(shù)學(xué)制圖學(xué)，是研究將地球橢球面(或球面)描寫到地圖平面上，建立地圖數(shù)學(xué)基礎(chǔ)的一門科學(xué)，它在地圖制作和應(yīng)用中起著“基礎(chǔ)”和“骨架”作用，是地圖編制前首要考慮的問題，同時它又是現(xiàn)代地圖學(xué)的重要組成部分[4-7]。地圖投影是地圖的空間數(shù)學(xué)基礎(chǔ)，是地圖學(xué)的理論基礎(chǔ)，是空間信息處理、傳輸和解譯各個階段必不可少的工具，其相關(guān)理論和方法在地圖制圖[8]、地理信息系統(tǒng)[9-10]、空間分析[11]、遙感制圖[12]、大地測量[13]、航海導(dǎo)航[14]等領(lǐng)域中已經(jīng)得到了廣泛的應(yīng)用。

地圖投影的實質(zhì)在于建立地球橢球面(或球面)與地圖平面上的點之間的對應(yīng)關(guān)系。地圖投影計算是地圖數(shù)學(xué)基礎(chǔ)中的重要組成部分，并且是相當(dāng)復(fù)雜的數(shù)學(xué)計算。地圖投影需要處理涉及參考橢球的各類數(shù)學(xué)分析問題，從而不可避免地會遇到大量的橢球偏心率冪級數(shù)展開、隱函數(shù)高階導(dǎo)數(shù)求取、復(fù)變函數(shù)運算等一系列復(fù)雜的符號推導(dǎo)過程。國內(nèi)外許多著名的地圖投影學(xué)者在地圖投影正反解、地圖投影變換、海圖投影及變換方面進行了大量深入系統(tǒng)的研究，取得了一批顯著的成果[15-36]。

值得注意的是，雖然前人對地圖投影中的數(shù)學(xué)分析問題進行了大量的研究，并取得了豐富的研究成果，但由于歷史條件和分析手段的限制，地圖投影數(shù)學(xué)分析問題主要依靠人工推導(dǎo)完成，不僅級數(shù)展開式的次數(shù)不會很高，而且導(dǎo)出的表達式繁瑣冗長，計算效率不高；有時為了計算上的方便采取近似處理，難免會導(dǎo)致一些小的偏差，影響了計算精度；有些算法則包含復(fù)雜的迭代計算，理論分析不甚方便；還有些算法表現(xiàn)為適用于特定參考橢球的數(shù)值形式，不便于推廣應(yīng)用。因此，地圖投影數(shù)學(xué)分析問題解決得并非完美和理想。

計算機代數(shù)(computer algebra)在很多時候又被廣義地理解為“符號計算”，成為與所謂“數(shù)值計算”相對的概念。計算機代數(shù)系統(tǒng)是計算機科學(xué)與數(shù)學(xué)分析和代數(shù)推導(dǎo)相結(jié)合的產(chǎn)物。計算機代數(shù)系統(tǒng)的優(yōu)越性主要在于它能夠進行大規(guī)模的代數(shù)運算，在一定程度上可以使科學(xué)研究和工程技術(shù)人員從枯燥煩瑣的數(shù)學(xué)分析和代數(shù)推理中解脫出來，從而有效地提高工作效率，完成人工推導(dǎo)難以實現(xiàn)的運算過程，同時其程序化設(shè)計可以保證結(jié)果的準(zhǔn)確性，目前流行的計算機代數(shù)系統(tǒng)有Mathematica、Maple、Mathcad等[37]。

近年來，海軍工程大學(xué)地圖投影研究團隊和國內(nèi)外地圖投影研究學(xué)者，利用具有嚴格解析意義的計算機代數(shù)分析方法，借助計算機代數(shù)系統(tǒng)強大的數(shù)學(xué)分析能力，對地圖投影領(lǐng)域中的一些典型數(shù)學(xué)分析過程進行了系統(tǒng)的研究，推導(dǎo)和建立了一系列理論上更為嚴密、形式上更為簡單、精度上更為精確的地圖投影新公式和新算法，實現(xiàn)了地圖投影在一些具體數(shù)學(xué)分析問題上的突破和創(chuàng)新，豐富和完善了地圖投影的理論體系。本文從橢球各緯度間正反解符號表達式、不同變形性質(zhì)地圖投影間的直接變換、高斯投影的復(fù)變函數(shù)表示、斜軸墨卡托投影數(shù)學(xué)分析、極區(qū)海圖投影及變換等5個方面，對地圖投影計算機代數(shù)分析的近期研究成果進行綜述分析，以期為國內(nèi)同行提供參考。

1 橢球各緯度間正反解符號表達式

大地緯度、地心緯度、歸化緯度、等距離緯度、等角緯度、等面積緯度是在測量和地圖學(xué)中常見的6種緯度，實際計算中經(jīng)常會遇到大地緯度和5種輔助緯度間的正反解問題[38]。

地心緯度、歸化緯度與大地緯度之間存在明確的正切關(guān)系式，文獻[39]利用拉格朗日共軛級數(shù)將其正反解展開式表示為正弦函數(shù)倍角項形式，但限于人工推導(dǎo)僅展至6倍。文獻[40]借助計算機代數(shù)系統(tǒng)將倍角項擴展至10倍，并將展開式系數(shù)統(tǒng)一表示為橢球偏心率的冪級數(shù)形式，給出了我國常用大地坐標(biāo)系下展開式的系數(shù)值，精度分析表明，歸化緯度展開式的精度優(yōu)于10-9″，地心緯度展開式的精度優(yōu)于10-8″，均較傳統(tǒng)公式提高3個數(shù)量級。文獻[41]推導(dǎo)出了以歸化緯度、地心緯度解算子午線弧長的展開公式，同時又根據(jù)拉格朗日反演定理，得到了由子午線弧長反解歸化緯度、地心緯度的直接公式，正反解精度均高于傳統(tǒng)基于大地緯度的展開式。

等距離緯度、等角緯度和等面積緯度都是大地緯度的函數(shù)，在等距離投影、等角投影和等面積投影表達式中經(jīng)常用到它們的正解展開式，即將其展開為關(guān)于大地緯度的正弦函數(shù)倍角形式。這一過程涉及非常復(fù)雜的冪級數(shù)展開和復(fù)合函數(shù)高階導(dǎo)數(shù)的求取，以往大都通過人工推導(dǎo)完成，推導(dǎo)過程復(fù)雜冗長，通常只能展至sin 8B(B為大地緯度)，并且為推導(dǎo)方便采取的某些近似往往導(dǎo)致展開式系數(shù)的高階項存在偏差。反解時由于大地緯度一般表現(xiàn)為這些緯度的隱函數(shù)或反函數(shù)形式，有的為非常復(fù)雜的超越函數(shù)，多采用基于正解公式的迭代法，這種方法不但計算效率低，而且最主要的是理論分析不甚方便；另一種方法是直接進行反解變換，即所謂的直接法。文獻[22]經(jīng)過復(fù)雜的Lagrange級數(shù)展開，文獻[42]利用變系數(shù)線性插值方法，給出了以上3種緯度反解的直接展開式，但是由于歷史條件的限制，其間許多推導(dǎo)過程都由人工完成，不僅展開式項數(shù)不高，而且展開式系數(shù)是原正解展開式系數(shù)的多項式形式，不便于記憶，計算也比較復(fù)雜，在實際應(yīng)用中多以具體的數(shù)值形式給出，不便于推廣應(yīng)用。

借助計算機代數(shù)分析方法和計算機代數(shù)系統(tǒng)，文獻[43]對等距離緯度、等角緯度和等面積緯度正解展開式進行了重新推導(dǎo)，發(fā)現(xiàn)并糾正了傳統(tǒng)正解展開式系數(shù)高階項中存在的偏差；文獻[43—46]分別采用冪級數(shù)展開法、Hermite插值法、Lagrange級數(shù)法、符號迭代法推導(dǎo)出了形式一致、展開式系數(shù)完全相同的反解展開式，其中三角函數(shù)倍角項展至sin 8B(B為大地緯度)，從而相互印證了4種方法的正確性，與傳統(tǒng)反解展開式不同的是，反解系數(shù)不再是正解展開式系數(shù)的多項式形式或具體的數(shù)值形式，而是統(tǒng)一表示為橢球偏心率e的冪級數(shù)形式，對于不同的橢球，只需將橢球偏心率代入即可得到該橢球下的相關(guān)系數(shù)，便于推廣使用。在此基礎(chǔ)上，文獻[40]系統(tǒng)推導(dǎo)出了三角函數(shù)倍角項展至sin 10B、系數(shù)展至e10的等距離緯度、等角緯度和等面積緯度正反解展開式，并給出了我國常用大地坐標(biāo)系下展開式的系數(shù)值，精度分析表明，計算機代數(shù)系統(tǒng)下推導(dǎo)出的正反解展開式精度較傳統(tǒng)人工導(dǎo)出的展開式最少提高2個數(shù)量級，最多提高4個數(shù)量級，具體情況如表1所示。

表1等距離緯度、等角緯度和等面積緯度正反解展開式精度比較

Tab.1Theaccuracycomparisonofforwardandinverseexpansionsofrectifying,conformalandauthaliclatitudes

正反解傳統(tǒng)人工導(dǎo)出展開式精度/(″)計算機代數(shù)系統(tǒng)導(dǎo)出展開式精度/(″)精度提高數(shù)量級等距離緯度正解10-510-72等角緯度正解10-410-84等面積緯度正解10-410-84等距離緯度反解10-410-84等角緯度反解10-410-84等面積緯度反解10-410-84

文獻[47]對測量和地圖學(xué)中6種常用緯度進行了系統(tǒng)比較，借助計算機代數(shù)系統(tǒng)推導(dǎo)出了常用緯度間的差異極值點及對應(yīng)差異極值的符號表達式，并將其表示為關(guān)于偏心率e的冪級數(shù)形式；以CGCS2000橢球為例，對各緯度間的差異進行了數(shù)值分析和對比。結(jié)果表明，5種輔助緯度與大地緯度的差異極值點均在π/4右側(cè)；地心緯度與大地緯度差異極值最大，歸化緯度與大地緯度差異極值最小。

文獻[48]利用等距離緯度、等角緯度和等面積緯度這3種緯度和大地緯度間的正反解展開式，導(dǎo)出了它們之間變換的直接展開式，并將式中系數(shù)統(tǒng)一表示為橢球偏心率e的冪級數(shù)形式并展至e10，解決了不同參考橢球下的變換問題，精度分析表明，直接展開式的計算精度優(yōu)于10-8″，滿足地圖投影精密計算的需要。

2 不同變形性質(zhì)地圖投影間的直接變換

地圖投影根據(jù)投影性質(zhì)可分為等距離投影、等角投影和等面積投影，根據(jù)投影后經(jīng)緯線形狀可分為圓柱投影、圓錐投影和方位投影。前人對這些投影間變換問題已進行了不少研究，特別是我國已故著名地圖學(xué)家楊啟和教授在地圖投影變換領(lǐng)域取得了令國內(nèi)外同行矚目的成果，楊啟和教授經(jīng)過幾十年的潛心研究，畢其一生之精力，集其研究之大成，于1989年出版了《地圖投影變換原理與方法》這部巨著，使這個研究領(lǐng)域大大前進一步，該書于2000年被英國Taylor出版公司譯為英文出版，是我國地圖投影方面第一部在國外出版的學(xué)術(shù)著作，標(biāo)志著我國地圖投影變換研究在國際上的領(lǐng)先地位，對于推動地圖學(xué)各分支學(xué)科的現(xiàn)代化水平和提高我國地圖學(xué)在國際上的地位具有重要的意義。

圖1 地圖投影間接變換和直接變換示意圖Fig.1 Sketch of indirect and direct transformations of map projection

值得注意的是，雖然前人對地圖投影變換方法進行了大量的研究，并在一定程度上解決了許多問題，但還不夠全面和充分，特別是對于地球橢球模型下不同變形性質(zhì)間的投影變換，采用的是間接變換法，需要通過中間過渡的方法，反解出原地圖投影點的橢球坐標(biāo)，再代入新投影中求得該點在新投影下的直角坐標(biāo)，而未能建立投影坐標(biāo)間的直接關(guān)系式，由于橢球體下兩種投影表達式較為復(fù)雜，而且變換中有時需要經(jīng)過煩瑣的迭代運算才能反解出大地緯度，使得計算過程變得相當(dāng)復(fù)雜，計算效率不高。因此，需要研究地圖投影直接變換問題，建立不同投影之間精確的對應(yīng)關(guān)系，表達編圖和制圖過程的數(shù)學(xué)實質(zhì)。地圖投影間接變換和直接變換示意圖如圖1所示。

子午線弧長、等量緯度和等面積緯度函數(shù)是等距離投影、等角投影和等面積投影中的重要變量，分別是等距離緯度、等角緯度和等面積緯度的函數(shù)，它們之間的變換是實現(xiàn)不同變形性質(zhì)地圖投影間變換的基礎(chǔ)，而這3種量之間的變換，傳統(tǒng)上是通過解算大地緯度間接實現(xiàn)的，計算公式復(fù)雜冗長，不便于應(yīng)用。

借助計算機代數(shù)分析方法和計算機代數(shù)系統(tǒng)，文獻[49]推導(dǎo)出了子午線弧長正反解公式，并將系數(shù)統(tǒng)一表示為橢球偏心率的冪級數(shù)形式，克服了傳統(tǒng)反解公式大多表示為具體數(shù)值形式的缺陷。文獻[50—52]在等距離緯度、等角緯度和等面積緯度正反解展開式基礎(chǔ)上，推導(dǎo)出了子午線弧長、等量緯度和等面積緯度函數(shù)之間變換的直接展開式，將展開式系數(shù)統(tǒng)一表示為橢球偏心率的冪級數(shù)形式，且展至e8。文獻[40]進一步將展開式系數(shù)擴展至e10，并給出了我國常用大地坐標(biāo)系下展開式的系數(shù)值，精度分析表明直接展開式的計算精度較傳統(tǒng)間接變換公式最少提高2個數(shù)量級，最多提高6個數(shù)量級，可以滿足地圖投影精密計算需要。

在此基礎(chǔ)上，文獻[40，53—54]系統(tǒng)地建立了正軸圓柱投影之間、正軸圓錐投影之間及正軸圓柱投影和正軸圓錐投影間的直接變換模型，避免了過去“圓柱→橢球→圓錐”間接變換導(dǎo)致的誤差，顯著地提高了投影變換的計算精度和計算效率。與傳統(tǒng)變換公式相比，這些直接變換模型形式更為簡單，理論更為嚴密，便于測量和制圖使用。

3 高斯投影的復(fù)變函數(shù)表示

高斯投影被廣泛用作地形圖的數(shù)學(xué)基礎(chǔ)，傳統(tǒng)的高斯投影正(反)解公式表示為經(jīng)差(橫坐標(biāo))的實數(shù)型冪級數(shù)形式，雖然有容易理解和直觀的優(yōu)點，但表達式復(fù)雜冗長，而且對于正解中子午線弧長的計算，給出的是適用于特定橢球的數(shù)值公式，反解中底點緯度則需要迭代求出，較為煩瑣[39]；特別是實用中需要分帶處理，經(jīng)常劃分為3°或6°帶，如圖2所示。針對這一問題，文獻[55—58]采用多種方法導(dǎo)出了適用范圍更廣的高斯投影算法，但由于沒有擺脫實數(shù)表示的限制，仍然沒有解決分帶問題。

圖2 高斯投影分帶示意圖Fig.2 Sketch of zone dividing of Gauss projection

鑒于高斯投影是一種等角投影，而復(fù)變函數(shù)作為一種強有力的數(shù)學(xué)方法，在等角投影中的優(yōu)勢是無可替代的，利用復(fù)變函數(shù)進行等角投影計算和公式推求具有簡單、方便、準(zhǔn)確的優(yōu)點，近年來已有學(xué)者注意到這一問題并進行了研究。文獻[59]采用橢球面在球面的局部描寫，給出了復(fù)變函數(shù)表示的高斯投影近似式，計算精度較低，只能達到0.3 m。文獻[60]討論了橫墨卡托投影的復(fù)變函數(shù)表示，但其給出的反解變換是在子午線弧長正解公式的基礎(chǔ)上迭代得到的，計算過于煩瑣。文獻[61]基于一種有效的遞推公式給出了任意帶寬的高斯投影復(fù)變函數(shù)解法，但所給公式較為復(fù)雜，且遞推過程耗時，計算效率較低。文獻[62]給出了用復(fù)變函數(shù)表示的高斯投影正反解的Fortran程序并進行了計算。

借助計算機代數(shù)分析方法和計算機代數(shù)系統(tǒng)，文獻[63—64]研究了橢球高斯投影的復(fù)變函數(shù)表示，將子午線弧長正反解公式拓展至復(fù)數(shù)域，導(dǎo)出了形式緊湊、結(jié)構(gòu)簡單的正反解公式，但在實際計算中需煩瑣的迭代運算，影響了計算效率。文獻[65]使用計算機代數(shù)系統(tǒng)研究了高斯投影復(fù)變換的數(shù)值計算方法，給出了復(fù)積分計算的積分級數(shù)分析法、橢圓積分函數(shù)法、直接積分法及利用系統(tǒng)函數(shù)求解變換方程。文獻[66]并將等量緯度與子午線弧長變換的直接展開式拓展至復(fù)數(shù)域拓展至復(fù)數(shù)域，導(dǎo)出了復(fù)變函數(shù)表示的高斯投影正反解非迭代公式，不僅提高了計算精度，而且在一定程度上也簡化了計算過程。文獻[67]推導(dǎo)出了球面高斯投影復(fù)變函數(shù)表達式，根據(jù)復(fù)變函數(shù)及常用初等函數(shù)的定義，經(jīng)過一系列數(shù)學(xué)分析過程，給出了球面高斯投影與橫軸墨卡托投影等價性的嚴格證明。文獻[68]在球面高斯投影復(fù)變函數(shù)表示的基礎(chǔ)上，通過一系列數(shù)學(xué)變換，推導(dǎo)出了球面高斯投影正反解、經(jīng)緯線投影方程、投影長度比及子午線收斂角公式的閉合式，與傳統(tǒng)高斯投影冪級數(shù)公式相比，理論上更嚴密、形式上更緊湊。文獻[69]在橢球高斯投影復(fù)變函數(shù)表示的基礎(chǔ)上，通過借助雙曲正弦和公式、雙曲函數(shù)與三角函數(shù)間的函數(shù)關(guān)系，經(jīng)過一系列恒等變換，將高斯投影復(fù)變函數(shù)坐標(biāo)公式及對應(yīng)的長度比、子午線收斂角公式等價變換為實數(shù)形式，具有“不分帶”的優(yōu)點，且表現(xiàn)形式更直觀、清晰。文獻[70]在高斯投影復(fù)變函數(shù)表示的基礎(chǔ)上，給出了基于復(fù)數(shù)等角緯度、復(fù)數(shù)底點緯度的推導(dǎo)過程，得到了高斯投影正反解公式的實數(shù)解，其計算精度達0.1 mm，突破經(jīng)典高斯投影帶寬的限制，適應(yīng)帶寬可達60°。文獻[71]借助計算機代數(shù)系統(tǒng)推導(dǎo)出了等角緯度關(guān)于子午線弧長的展開式，利用復(fù)變函數(shù)及雙曲正切函數(shù)理論，得到了不同中央經(jīng)線高斯投影間的變換公式，與傳統(tǒng)高斯投影換帶公式相比，具有更高的準(zhǔn)確度及更寬的應(yīng)用范圍。文獻[72—73]討論了常用等角投影及其解析變換的復(fù)變函數(shù)表示，給出了高斯投影、墨卡托投影和等角圓錐投影正反解的復(fù)變函數(shù)表示模型，在此基礎(chǔ)上系統(tǒng)地推導(dǎo)出了高斯投影、墨卡托投影和等角圓錐投影間解析變換的復(fù)變函數(shù)表達式。文獻[74]導(dǎo)出了拉格朗日投影正反解的復(fù)變函數(shù)表達式，在此基礎(chǔ)上系統(tǒng)地建立了該投影與高斯投影、墨卡托投影和等角圓錐投影間解析變換的復(fù)變函數(shù)表示模型。這些復(fù)數(shù)變換公式表現(xiàn)為含參考橢球第一偏心率的符號形式，可解決不同參考橢球下的變換問題，與傳統(tǒng)的實數(shù)變換公式相比，其形式更為簡單、理論更為嚴密。

4 斜軸墨卡托投影數(shù)學(xué)分析

墨卡托投影是一類常用的等角投影，因投影面與地球橢球相對位置的不同，可將其分為正軸、橫軸和斜軸墨卡托投影。正軸墨卡托投影適用于東西向赤道附近的測區(qū)，橫墨卡托投影適于南北縱向沿經(jīng)線的測區(qū)。顯然，存在一些地理區(qū)域并不以赤道、經(jīng)線、緯線為中線，而位于某一平面與地球或其他星體的交線附近，如圖3所示。針對這種情況的應(yīng)用需求，人們展開了對斜軸墨卡托投影的研究。

圖3 斜軸墨卡托投影示意圖Fig.3 Sketch of oblique Mercator projection

文獻[75]對旋轉(zhuǎn)橢球進行等角投影，得到橢球斜軸墨卡托投影，這種方法被用于美國阿拉斯加、馬來亞半島和婆羅洲、瑞士、意大利等國家或地區(qū)制圖。文獻[76]曾闡述了Hotine斜軸墨卡托投影是繪制地球資源衛(wèi)星類型數(shù)據(jù)最合適的投影方式，曾用于熱容量測繪衛(wèi)星的投影成像。文獻[77]基于正軸及橫軸墨卡托投影的概念，在確保變換后的赤道投影前后長度不變的前提下，給出了橢球情況下斜軸墨卡托投影的算法。文獻[36]引入復(fù)變函數(shù)對墨卡托投影體系進行了系統(tǒng)的研究，詳細介紹了球面、橢球情況下的正軸、橫軸、斜軸墨卡托投影方式。

隨著國民經(jīng)濟的飛速發(fā)展，我國高速鐵路建設(shè)、高等級公路建設(shè)、長距離的引水和輸油管線建設(shè)等工程日益興起。這些長大線狀工程橫跨東西南北，具有測區(qū)帶狀、路線長等顯著特點，傳統(tǒng)的橫軸高斯投影只適用于南北條帶區(qū)域的投影表示，應(yīng)用于東西跨度較大的長大線狀工程，將導(dǎo)致分帶過于頻繁、換帶計算煩瑣、信息要素表達零散等不足，而采用斜軸墨卡托投影可有效控制投影長度變形，避免分帶計算。文獻[78]介紹了斜圓柱投影的計算方法，使投影圓柱沿工程直伸中心線相切，從而克服了長度變形問題。文獻[79]討論了參考圓球的確定和定位，編制了計算軟件，實現(xiàn)了從常用坐標(biāo)系統(tǒng)向斜軸投影平面坐標(biāo)的轉(zhuǎn)換計算。文獻[80—83]結(jié)合鐵路線路的實測數(shù)據(jù)分析了斜墨卡托投影的應(yīng)用效果，驗證了方法的有效性。然而，上述文獻為計算方便，往往將測區(qū)表面近似為球面，同時未考慮高差對于投影變形的影響，使得算法受距離限制，不適于長線型線路。因此，為更精確地表示測區(qū)，在長線型線路設(shè)計中需考慮地球扁率的影響，而這一問題往往涉及非常復(fù)雜冗長的數(shù)學(xué)分析，人工推導(dǎo)難度極大。

文獻[84—86]提出將接近于長大線狀工程的法截面進行旋轉(zhuǎn)，得到法截面子午線橢球，借助計算機代數(shù)系統(tǒng)對法截面子午線橢球構(gòu)建和變換涉及的復(fù)雜數(shù)學(xué)分析問題進行了深入研究，建立了法截面子午線橢球高斯投影、空間幾何和工程應(yīng)用理論體系，避免了高斯分帶計算，減小了投影變形，在施工允許精度范圍內(nèi)將單個坐標(biāo)系線路長度延伸至1000 km以上，甚至達到7000 km，有效提高了高斯投影在東西跨度較大區(qū)域的使用效率。文獻[87]基于最小二乘、坐標(biāo)系轉(zhuǎn)換和橢球變換理論，借助計算機代數(shù)系統(tǒng)構(gòu)建了斜軸變形橢球高斯投影方法，避免了因分帶帶來的誤差，減小了投影綜合變形，其數(shù)學(xué)模型更為簡單，計算過程更為清晰。

5 極區(qū)海圖投影及其變換

隨著全球氣候變暖，極區(qū)商業(yè)價值及戰(zhàn)略地位日漸凸顯，對其進行資源調(diào)查及航線探索已成為近年來國際上爭相關(guān)注的熱點。海圖對于極區(qū)航行具有不可或缺的支撐作用，而中低緯度采用的墨卡托投影由于長度變形過大而不宜在極區(qū)使用，選擇合適的海圖投影方式對編制極區(qū)航行所需航海圖至關(guān)重要。

文獻[88]提出橫向和斜向?qū)Ш郊夹g(shù)，與之配合的投影方式應(yīng)為蘭伯特等角投影、極球面投影或橫墨卡托投影，但沒有對這些海圖投影方式特性及使用方式進行進一步的詳細分析與說明。文獻[89]提出了利用陀螺作為跟蹤儀的格網(wǎng)導(dǎo)航方法，建議在格網(wǎng)導(dǎo)航中使用斜向墨卡托投影。文獻[90]提出了一種基于橫墨卡托投影導(dǎo)航圖的極區(qū)橫向?qū)Ш郊夹g(shù)，將橫向地理坐標(biāo)系與橫墨卡托投影結(jié)合以支持極區(qū)導(dǎo)航。1989年美軍國防制圖局(Defense Mapping Agency，DMA)推薦北緯84°至北極和南緯80°至南極間的圖幅采用通用極球面投影(universal polar stereographic projection，UPS)。文獻[91]對極球面投影上用直線直接代替大圓航線所造成的航向誤差和位置誤差進行了定性分析，對極區(qū)常用投影方式的導(dǎo)航應(yīng)用研究更深入細致。文獻[92]建議北極區(qū)域使用等距離方位投影和等角方位投影、近北極區(qū)域使用蘭伯特圓錐投影和等距離圓錐投影。文獻[93]認為在尋找出一種新的極地投影之前，墨卡托投影和日晷投影仍將作為編制北極航海圖的首選投影。

以往對極區(qū)常用海圖投影方式的投影特性及可用性分析，只給出了籠統(tǒng)的定性分析初步結(jié)論，對其變形特性、誤差分布特性缺乏系統(tǒng)、深入、細致的理論推導(dǎo)和數(shù)值分析，特別是傳統(tǒng)人工推導(dǎo)的極區(qū)海圖投影算法存在以下有待解決的問題：基于球體模型導(dǎo)出，存在近似、精度不高的缺點；存在奇異問題或分帶現(xiàn)象，表達式不統(tǒng)一；未能建立投影坐標(biāo)直接變換的嚴密解析式，計算煩瑣。

借助計算機代數(shù)分析方法和計算機代數(shù)系統(tǒng)，文獻[94]分析了極球面投影的極區(qū)投影變形、經(jīng)緯線形狀及大圓航線與直線的逼近程度，結(jié)果表明極球面投影能夠克服墨卡托投影不足，可滿足極區(qū)海圖投影要求。文獻[95]針對現(xiàn)有研究中未充分考慮地球扁率對極區(qū)極球面投影誤差影響的問題，研究了雙重極球面投影下極區(qū)投影誤差與各重基準(zhǔn)緯度的關(guān)系。文獻[96]研究了極球面投影上直線代替大圓航線量測距離的可行性，提出了一種基于等距圓的距離量測方法，解決了極球面投影海圖上精確量測距離的問題。文獻[97]建立了等距離正圓柱投影平面上等角航線方程并對其曲率進行了分析，提出了一種可滿足給定精度要求的等角航線展繪算法。

鑒于陸圖采用的高斯投影具有等角優(yōu)勢，在中央子午線附近變形較小，如極區(qū)海圖采用高斯投影則便于陸海圖的統(tǒng)一，但高斯投影傳統(tǒng)實數(shù)公式僅在經(jīng)差較小的條帶內(nèi)適用，傳統(tǒng)復(fù)變函數(shù)表達式雖消除了分帶的限制，但是由于等量緯度在極點存在奇異，使得復(fù)變函數(shù)表達式難以在極區(qū)應(yīng)用。文獻[98—99]推導(dǎo)出了極區(qū)球面高斯投影非奇異公式，滿足了極區(qū)內(nèi)連續(xù)投影的需求，基于該公式推導(dǎo)了極區(qū)經(jīng)緯線投影方程，并結(jié)合日晷投影進行了長度變形及子午線偏移角分析，結(jié)果表明在極圈內(nèi)高斯投影長度變形小于日晷投影。文獻[100]在高斯投影復(fù)變函數(shù)的基礎(chǔ)上，引入復(fù)變等角緯度的概念，避免等量緯度在極點的奇異性，并將極點作為高斯投影的坐標(biāo)原點，建立了極區(qū)非奇異高斯投影復(fù)變函數(shù)表示形式，克服了傳統(tǒng)高斯投影分帶的缺陷，使得高斯投影在極區(qū)有一個統(tǒng)一完整的“一體化表示形式”，但在求得復(fù)數(shù)等角余緯度的過程中，使用了極點附近球近似的方法，一定程度上影響了正解推導(dǎo)過程的嚴密性。文獻[101]建立了等角余緯度和等量緯度間的嚴密關(guān)系式，對其進行復(fù)數(shù)開拓，推導(dǎo)出了理論嚴密的極區(qū)不分帶高斯投影正反解復(fù)數(shù)表達式，利用指數(shù)函數(shù)與三角函數(shù)間的關(guān)系式，進一步導(dǎo)出了極區(qū)高斯投影正反解實數(shù)表達式，在此基礎(chǔ)上繪制了極區(qū)不分帶高斯投影示意圖，如圖4所示。

圖4 北極地區(qū)不分帶高斯投影示意圖Fig.4 Sketch of non-zonal Gauss projection in arctic regions

文獻[102]研究了極球面投影、橫墨卡托投影和日晷投影的直接變換問題，系統(tǒng)地推導(dǎo)出了任意兩種投影間的直接變換公式，可供極區(qū)海圖投影變換和航海導(dǎo)航參數(shù)計算使用。文獻[103]討論了常用極區(qū)海圖投影及其變換，給出了各類球面或橢球面投影及其變換公式，以及長度比或角度變形公式，主要包括等角投影、等面積投影、等距離投影和日晷投影，著重研究了等角投影及等角投影間的解析變換。

6 結(jié)論與展望

借助現(xiàn)代計算機代數(shù)系統(tǒng)強大的數(shù)學(xué)分析功能和符號運算能力，將計算機代數(shù)分析方法應(yīng)用于地圖投影數(shù)學(xué)分析，不僅可以極大地提高地圖投影數(shù)學(xué)分析過程的效率，而且可以推導(dǎo)出過去人工推導(dǎo)難以完成的公式和算法，實現(xiàn)地圖投影學(xué)特定領(lǐng)域某些數(shù)學(xué)分析過程的創(chuàng)新。隨著科學(xué)技術(shù)特別是空間科學(xué)和計算機科學(xué)技術(shù)的發(fā)展，地圖投影的理論越來越科學(xué)化和體系化，地圖投影的方法和種類越來越多樣化，地圖投影學(xué)已突破了傳統(tǒng)的研究內(nèi)容和研究方法，正在向進一步深化和多種專業(yè)知識的結(jié)合方向發(fā)展，新的研究和應(yīng)用領(lǐng)域不斷地開拓。筆者認為，地圖投影計算機代數(shù)分析有待進一步開展的研究工作主要有以下幾個方面：

(1) 空間地圖投影計算機代數(shù)分析。空間地圖投影[104-106]是近40年才出現(xiàn)的一個新的研究領(lǐng)域，是圖像數(shù)學(xué)基礎(chǔ)研究的前沿課題，是衛(wèi)星遙感制圖學(xué)發(fā)展必然產(chǎn)生的問題。空間地圖投影涉及地球形狀、地球自轉(zhuǎn)、衛(wèi)星軌道攝動等非常復(fù)雜的數(shù)學(xué)分析問題，傳統(tǒng)人工推演均采用一定的近似分析方法和數(shù)值積分，借助計算機代數(shù)系統(tǒng)強大的符號運算能力，可以推導(dǎo)和建立理論上更為嚴密、形式上更為簡單、精度上更為精確的空間地圖投影新公式和新算法，揭示各類空間地圖投影復(fù)雜數(shù)學(xué)模型和海量數(shù)據(jù)背后隱藏的規(guī)律，實現(xiàn)空間地圖投影在一些具體數(shù)學(xué)分析問題上的突破和創(chuàng)新。

(2) 中小比例尺地圖投影計算機代數(shù)設(shè)計與分析。中小比例尺地圖投影在過去有過一些研究和分析，但由于受時代的局限，在分析深度、廣度和精細程度上都顯得比較粗淺，絕非完美。許多投影是人工擬合得出的數(shù)值形式，參數(shù)選取有一定的隨機性，缺乏嚴格的數(shù)學(xué)基礎(chǔ)。利用先進的計算機代數(shù)分析工具，可以全面對比各類中小比例尺地圖投影特點，對以前缺乏數(shù)學(xué)基礎(chǔ)的各類數(shù)值擬合投影方法盡可能給出符號化的通用數(shù)學(xué)準(zhǔn)則和表達式，提高我國中小比例尺地圖投影的理論研究和設(shè)計水平。

(3) 海圖投影及航線繪算計算機代數(shù)分析。海圖投影及航線繪算是海圖制圖和航線設(shè)計的理論基礎(chǔ)，在海洋監(jiān)測與調(diào)查、海洋劃界和航海中有著廣泛的應(yīng)用。傳統(tǒng)的海圖投影及航線繪算公式和算法，大多表現(xiàn)為數(shù)值形式，普適性不高，符號形式的算法部分存在高階項誤差，同時在極區(qū)的應(yīng)用和研究也很有限。利用計算機代數(shù)分析方法，借助計算機代數(shù)系統(tǒng)對其進行系統(tǒng)分析，可以將以往算法的數(shù)值形式改進為符號形式，建立更適合的新模型和新算法，完善海圖投影及航線繪算的理論體系。

[1] 王家耀, 成毅. 論地圖學(xué)的屬性和地圖的價值[J]. 測繪學(xué)報, 2015, 44(3): 237-241. DOI: 10.11947/j.AGCS.2015.20140406.

WANG Jiayao, CHENG Yi. Discussions on the Attributes of Cartography and the Value of Map[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(3): 237-241. DOI: 10.11947/j.AGCS.2015.20140406.

[2] 王家耀, 何宗宜, 蒲英霞, 等. 地圖學(xué)[M]. 北京: 測繪出版社, 2016.

WANG Jiayao, HE Zongyi, PU Yingxia, et al. Cartography[M]. Beijing: Surveying and Mapping Press, 2016.

[3] 李志林, 劉啟亮, 高培超. 地圖信息論: 從狹義到廣義的發(fā)展回顧[J]. 測繪學(xué)報, 2016, 45(7): 757-767. DOI: 10.11947/j.AGCS.2016.20160235.

LI Zhilin, LIU Qiliang, GAO Peichao. Entropy-Based Cartographic Communication Models: Evolution from Special to General Cartographic Information Theory[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(7): 757-767. DOI: 10.11947/j.AGCS.2016.20160235.

[4] 吳忠性, 楊啟和. 數(shù)學(xué)制圖學(xué)原理[M]. 北京: 測繪出版社, 1989.

WU Zhongxing, YANG Qihe. Principles of Mathematical Cartography[M]. Beijing: Surveying and Mapping Press, 1989.

[5] 楊啟和. 地圖投影變換原理與方法[M]. 北京: 解放軍出版社, 1989.

YANG Qihe. Principles and Methods of Map Projection Transformation[M]. Beijing: PLA Press, 1989.

[6] YANG Qihe, SNYDER J P, TOBLER W R. Map Projection Transformation: Principles and Applications[M]. London: Taylor & Francis, 2000.

[7] 呂曉華, 李少梅. 地圖投影原理與方法[M]. 北京: 測繪出版社, 2016.

Lü Xiaohua, LI Shaomei. Principles and Methods of Map Projection[M]. Beijing: Surveying and Mapping Press,2016.

[8] 何宗宜. 計算機地圖制圖[M]. 北京: 測繪出版社, 2008.

HE Zongyi. Computer Aided Cartography[M]. Beijing: Surveying and Mapping Press, 2008.

[9] 閭國年, 吳平生, 周曉波. 地理信息科學(xué)導(dǎo)論[M]. 北京: 中國科學(xué)技術(shù)出版社, 1999.

Lü Guonian, WU Pingsheng, ZHOU Xiaobo. Introduction to Geographical Information Science[M]. Beijing: China Science Technology Press, 1999.

[10] 龔健雅. 地理信息系統(tǒng)基礎(chǔ)[M]. 北京: 科學(xué)出版社, 2001.

GONG Jianya. Foundation of Geographical Information System[M]. Beijing: Science Press, 2001.

[11] 郭仁忠. 空間分析[M]. 2版. 北京: 高等教育出版社, 2001.

GUO Renzhong. Spatial Analysis[M]. 2nd ed. Beijing: Higher Education Press, 2001.

[12] 龐小平. 遙感制圖與應(yīng)用[M]. 北京: 測繪出版社, 2016.

PANG Xiaoping. Remote Sensing Mapping and Its Application[M]. Beijing: Surveying and Mapping Press, 2016.

[13] 孔祥元, 郭際明, 劉宗泉. 大地測量學(xué)基礎(chǔ)[M]. 武漢: 武漢大學(xué)出版社, 2005.

KONG Xiangyuan, GUO Jiming, LIU Zongquan. Foundation of Geodesy[M]. Wuhan: Wuhan University Press, 2005.

[14] 邊少鋒, 紀(jì)兵, 李厚樸. 衛(wèi)星導(dǎo)航系統(tǒng)概論[M]. 2版. 北京: 測繪出版社, 2016.

BIAN Shaofeng, JI Bing, LI Houpu. Introduction to Satellite Navigation System[M]. 2nd ed. Beijing: Surveying and Mapping Press, 2016.

[15] 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: U.S. Government Printing Office, 1921.

[16] THOMAS P D. Conformal Projections in Geodesy and Cartography[M]. Washington D C: U.S. Government Printing Office, 1952.

[17] 方俊. 地圖投影學(xué)-第二冊[M]. 北京: 科學(xué)出版社, 1958.

FANG Jun. Map Projection[M]. Beijing: Science Press, 1958.

[18] 方炳炎. 地圖投影學(xué)[M]. 北京: 地圖出版社, 1978.

FANG Bingyan. Map Projection[M]. Beijing: Map Press, 1978.

[19] SNYDER J P. Map Projections: A Working Manual[M]. Washington D C: U.S. Government Printing Office, 1987.

[20] 華棠. 海圖數(shù)學(xué)基礎(chǔ)[M]. 北京: 海潮出版社, 1985.

HUA Tang. Mathematical Foundation of Chart[M]. Beijing: Haichao Press, 1985.

[21] 丁佳波. 關(guān)于等角投影解析變換的補充[J]. 測繪學(xué)報, 1985, 11(1): 46-50.

DING Jiabo. A Complement to the Analytical Transformation Between Conformal Projections[J]. Acta Geodaetica et Cartographica Sinica, 1985, 11(1): 46-50.

[22] 孫群, 楊啟和. 底點緯度解算以及等量緯度和面積函數(shù)反解問題的探討[J]. 解放軍測繪學(xué)院學(xué)報, 1985(2): 64-75.

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.

[23] BOWRING B R. The Transverse Mercator Projection-A Solution by Complex Numbers[J]. Survey Review, 1990, 30(237): 325-342.

[24] 胡毓矩, 龔劍文, 黃偉. 地圖投影[M]. 北京: 測繪出版社, 1981.

HU Yuju, GONG Jianwen, HUANG Wei. Map Projection[M]. Beijing: Surveying and Mapping Press, 1981.

[25] CHENG Yang. The Conformal Space Projection[J]. Cartography and Geographic Information Systems, 1996, 23(1): 37-50.

[26] GRAFAREND E W, SYFFUS R. The Solution of the Korn-Lichtenstein Equations of Conformal Mapping: the Direct Generation of Ellipsoidal Gauβ-Krüger Conformal Coordinates or the Transverse Mercator Projection[J]. Journal of Geodesy, 1998, 72(5): 282-293.

[27] GRAFAREND E W, YOU R J, SYFFUS R. Map Projection[M]. 2nd ed. Berlin: Springer-Verlag, 2014.

[28] KESSLER F C. A Visual Basic Algorithm for the Winkel Tripel Projection[J]. Cartography and Geographic Information Science, 2000, 27(2): 177-183.

[29] OZTAN O C, IPBUKER C, ULUGTEKIN N. A Numerical Approach to Pseudo-Projections on Example Franz Mayr Projection[J]. Journal of General Command of Mapping, 2001, 125: 37-50.

[30] 呂曉華, 劉宏林. 地圖投影數(shù)值變換方法綜合評述[J]. 測繪學(xué)院學(xué)報, 2002, 19(2): 150-153.

Lü Xiaohua, LIU Honglin. A Comprehensive Appraisal of Numerical Transformation Method for Map Projection[J]. Journal of Institute of Surveying and Mapping, 2002, 19(2): 150-153.

[31] 任留成. 空間投影理論及其在遙感技術(shù)中的應(yīng)用[M]. 北京: 科學(xué)出版社, 2003.

REN Liucheng. Theory of Space Map Projection with Its Applications in Remote Sensing[M]. Beijing: Science Press, 2003.

[32] IPBUKER C, BILDIRICI I O. Computer Program for the Inverse Transformation of the Winkel Projection[J]. Journal of Surveying Engineering, 2005, 131(4): 125-129.

[33] 孫達, 蒲英霞. 地圖投影[M]. 南京: 南京大學(xué)出版社, 2005.

SUN Da, PU Yingxia. Map Projection[M]. Nanjing: Nanjing University Press, 2005.

[34] 鐘業(yè)勛. 數(shù)理地圖學(xué)—地圖學(xué)及其數(shù)學(xué)原理[M]. 北京: 測繪出版社, 2007.

ZHONG Yexun. Mathematical Cartography-Cartography and Its Mathematical Elements[M]. Beijing: Surveying and Mapping Press, 2007.

[35] KARNEY C F F. Transverse Mercator with an Accuracy of A Few Nanometers[J]. Journal of Geodesy, 2011, 85(8): 475-485.

[36] OSBORNE P. The Mercator Projection[M]. Edinburgh: Edinburgh University Press, 2013.

[37] 邊少鋒, 許江寧. 計算機代數(shù)系統(tǒng)與大地測量數(shù)學(xué)分析[M]. 北京: 國防工業(yè)出版社, 2004.

BIAN Shaofeng, XU Jiangning. Computer Algebra System and Mathematical Analysis in Geodesy[M]. Beijing: National Defense Industry Press, 2004.

[38] 楊啟和. 測量和地圖學(xué)中應(yīng)用的六種緯度及其變換關(guān)系式[J]. 測繪科技通訊, 1995, 18(3): 14-19.

YANG Qihe. Six Kinds of Latitudes with Their Transformation Relations in Surveying and Mapping[J]. Science and Technology Newsletter, 1995, 18(3): 14-19.

[39] 熊介. 橢球大地測量學(xué)[M]. 北京: 解放軍出版社, 1988.

XIONG Jie. Ellipsoidal Geodesy[M]. Beijing: PLA Press, 1988.

[40] 李厚樸, 邊少鋒, 鐘斌. 地理坐標(biāo)系計算機代數(shù)精密分析理論[M]. 北京: 國防工業(yè)出版社, 2015.

LI Houpu, BIAN Shaofeng, ZHONG Bin. Precise Analysis Theory of Geographic Coordinate System by Computer Algebra[J]. Beijing: National Defense Industry Press, 2015.

[41] 過家春, 李厚樸, 莊云玲, 等. 依不同緯度變量的子午線弧長正反解公式的級數(shù)展開[J]. 測繪學(xué)報, 2016, 45(5): 560-565. DOI: 10.11947/j.AGCS.2016.20140575.

GUO Jiachun, LI Houpu, ZHUANG Yunling, et al. Series Expansion for Direct and Inverse Solutions of Meridian in Terms of Different Latitude Variables[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(5): 560-565. DOI: 10.11947/j.AGCS.2016.20140575.

[42] 楊啟和, 楊曉梅. 測量和地圖學(xué)中應(yīng)用的三種緯度函數(shù)及其反解變換的線性插值方法[J]. 測繪學(xué)報, 1997, 26(1): 92-94, 封3.

YANG Qihe, YANG Xiaomei. Three Kinds of Latitudes in Surveying and Mapping and the Linear Interpolation Method for Inverse Solution[J]. Acta Geodaetica et Cartographica Sinica, 1997, 26(1): 92-94, 封3.

[43] 邊少鋒, 紀(jì)兵. 等距離緯度等量緯度和等面積緯度展開式[J]. 測繪學(xué)報, 2007, 36(2): 218-223. DOI: 10.3321/j.issn:1001-1595.2007.02.018.

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. DOI: 10.3321/j.issn:1001-1595.2007.02.018.

[44] 李厚樸, 邊少鋒. 輔助緯度反解公式的Hermite插值法新解[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2008, 33(6): 623-626.

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.

[45] 王瑞, 李厚樸. 輔助緯度反解公式的Lagrange級數(shù)法推演[J]. 海洋測繪, 2008, 28(3): 18-23.

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.

[46] 李忠美, 邊少鋒, 孔海英. 符號迭代法解算橢球大地測量學(xué)反問題[J]. 海洋測繪, 2013, 33(1): 27-29, 33.

LI Zhongmei, BIAN Shaofeng, KONG Haiying. Symbolic Iterative Method for Solving Inverse Problems in Ellipsoidal Geodesy[J]. Hydrographic Surveying and Charting, 2013, 33(1): 27-29, 33.

[47] 李忠美, 李厚樸, 邊少鋒. 常用緯度差異極值符號表達式[J]. 測繪學(xué)報, 2014, 43(2): 214-220. DOI: 10.13485/j.cnki.11-2089.2014.0031.

LI Zhongmei, LI Houpu, BIAN Shaofeng. Symbolic Expressions of Difference Extrema Between Regular Latitudes[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(2): 214-220. DOI: 10.13485/j.cnki.11-2089.2014.0031.

[48] 李厚樸, 邊少鋒, 劉敏. 地圖投影中三種緯度間變換直接展開式[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2013, 38(2): 217-220.

LI Houpu, BIAN Shaofeng, LIU Min. Direct Expansions of Transformations Between Three Kinds of Latitudes Used in Map Projection[J]. Geomatics and Information Science of Wuhan University, 2013, 38(2): 217-220.

[49] BIAN Saofeng, CHEN Yongbing. Solving An Inverse Problem of a Meridian Arc in Terms of Computer Algebra System[J]. Journal of Surveying Engineering, 2006, 132(1): 7-10.

[50] 李厚樸, 劉敏, 孔海英, 等. 子午線弧長和等面積緯度函數(shù)變換的直接展開式[J]. 海洋測繪, 2011, 31(1): 17-19.

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.

[51] 李厚樸, 邊少鋒, 王克平. 等距離緯度和等量緯度變換的直接展開式[C]∥第二十三屆海洋測繪綜合性學(xué)術(shù)研討會. 呼倫貝爾:[s.n.], 2013: 189-195.

LI Houpu, BIAN Shaofeng, WANG Keping. The Direct Expansions for Transformations Rectifying and Isometric Latitudes[C]∥Hydrographic Surveying and Charting Symposium of 24th General Meeting. Hulunbeier:[s.n.], 2013: 149-153.

[52] 李厚樸, 邊少鋒, 陳良友. 等面積緯度函數(shù)和等量緯度變換的直接解算公式[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2011, 36(7): 843-846.

LI Houpu, BIAN Shaofeng, CHEN Liangyou. The Direct Calculating Formulae for Transformations between Authalic Latitude Function and Isometric Latitude[J]. Geomatics and Information Science of Wuhan University, 2011, 36(7): 843-846.

[53] 李厚樸, 邊少鋒. 不同變形性質(zhì)正軸圓柱投影和正軸圓錐投影間的直接變換[J]. 測繪學(xué)報, 2012, 41(4): 536-542.

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.

[54] BIAN Shaofeng, LI Houpu. Mathematical Analysis in Cartography By Means of Computer Algebra System[C]∥Carlos Bateira. Cartography: A Tool for Spatial Analysis. Croatia: InTech, 2012.

[55] ENGSAGER K E, PODER K. A Highly Accurate World Wide Algorithm for the Transverse Mercator Mapping (Almost)[C]∥Proceedings of XXIII International Cartographic Conference (ICC2007). Moscow: [s.n.], 2007.

[56] BERMEJO-SOLERA M, OTERO J. Global Optimization of the Gauss Conformal Mappings of An Ellipsoid to A Sphere[J]. Journal of Geodesy, 2010, 84(8): 481-489.

[57] KARNEY C F F. Transverse Mercator with an Accuracy of A Few Nanometers[J]. Journal of Geodesy, 2011, 85(8): 475-485.

[58] KAWASE K. A General Formula for Calculating Meridian Arc Length and Its Application to Coordinate Conversion in the Gauss-Krüger Projection[J]. Bulletin of the Geospatial Information Authority of Japan, 2011, 59: 1-13.

[59] 程陽. 復(fù)變函數(shù)與等角投影[J]. 測繪學(xué)報, 1985, 14(1): 51-60.

CHENG Yang. Complex Function and Conformal Projection[J]. Acta Geodaetica et Cartographica Sinica, 1985, 14(1): 51-60.

[60] BOWRING B R. The Transverse Mercator Projection——A Solution by Complex Numbers[J]. Survey Review, 1990, 30(237): 325-342.

[61] KLOTZ J. Eine Analytische L?sung der Gauss-Krüger-Abbildung[J]. Zeitschrift für Versicherungswesen, 1993, 118(3): 106-116.

[62] SCHUHR P. Transformationen Zwischen Ellipsoidischen Geographischen Konformen Gauss-Krüger-bzw. UTM-Koordinaten[J]. Forum, 1995(5): 259-264.

[63] 邊少鋒, 柴洪洲, 金際航. 大地坐標(biāo)系與大地基準(zhǔn)[M]. 北京: 國防工業(yè)出版社, 2005.

BIAN Shaofeng, CHAI Hongzhou, JIN Jihang. Geodetic Coordinate System and Datum[M]. Beijing: National Defence Industry Press, 2005.

[64] 李厚樸, 邊少鋒. 高斯投影的復(fù)變函數(shù)表示[J]. 測繪學(xué)報, 2008, 37(1): 5-9. DOI: 10.3321/j.issn:1001-1595.2008.01.002.

LI Houpu, BIAN Shaofeng. The Expressions of Gauss Projection By Complex Numbers[J]. Acta Geodaetica et Cartographica Ainica, 2008, 37(1): 5-9. DOI: 10.3321/j.issn:1001-1595.2008.01.002.

[65] 劉大海. 高斯投影復(fù)變換的數(shù)值計算方法[J]. 測繪科學(xué)技術(shù)學(xué)報, 2012, 29(1): 9-11.

LIU Dahai. Gauss Projection Numerical Method for Complex Transformation[J]. Journal of Geomatics Science and Technology, 2012, 29(1): 9-11.

[66] 李厚樸, 王瑞, 邊少鋒. 復(fù)變函數(shù)表示的高斯投影非迭代公式[J]. 海洋測繪, 2009, 29(6): 17-20.

LI Houpu, WANG Rui, BIAN Shaofeng. The Non-Iterative Formulae of Gauss Projection By Complex Numbers[J]. Hydrographic Surveying and Charting, 2009, 29(6): 17-20.

[67] 李忠美, 于金星, 李厚樸, 等. 高斯投影與橫墨卡托投影等價性證明[J]. 海洋測繪, 2013, 33(3): 17-20.

LI Zhongmei, YU Jinxing, LI Houpu, et al. The Verification of Equivalence Between Gauss and Transverse Mercator Projections[J]. Hydrographic Surveying and Charting, 2013, 33(3): 17-20.

[68] 劉強, 邊少鋒, 李忠美. 球面高斯投影及其變形的閉合公式[J]. 海軍工程大學(xué)學(xué)報, 2015, 27(1): 45-49, 58.

LIU Qiang, BIAN Shaofeng, LI Zhongmei. Closed Formulae of Gauss Projection with Its Deformation on Sphere[J]. Journal of Naval University of Engineering, 2015, 27(1): 45-49, 58.

[69] 邊少鋒, 劉強, 李忠美. 不分帶的高斯投影實數(shù)公式[J]. 測繪通報, 2016(6): 6-9. DOI: 10.13474/j.cnki.11-2246.2016.0178.

BIAN Shaofeng, LIU Qiang, LI Zhongmei. Expressions of Non-Zonal Gauss Projection in Real Form[J]. Bulletin of Surveying and Mapping, 2016(6): 6-9. DOI: 10.13474/j.cnki.11-2246.2016.0178.

[70] 金立新, 許常文, 魏桂華. 高斯投影復(fù)變函數(shù)表示的實數(shù)解[J]. 海洋測繪, 2017, 37(2): 27-31.

JIN Lixin, XU Changwen, WEI Guihua. Real Solution of Gauss Projection Based on Complex Function[J]. Hydrographic Surveying and Charting, 2017, 37(2): 27-31.

[71] 劉強, 邊少鋒, 李忠美. 利用復(fù)變函數(shù)實現(xiàn)高斯換帶的方法[J]. 海軍工程大學(xué)學(xué)報, 2016, 28(1): 15-19.

LIU Qiang, BIAN Shaofeng, LI Zhongmei. Transformation Between Gauss Projections with Different Central Meridians Based on Complex Function[J]. Journal of Naval University of Engineering, 2016, 28(1): 15-19.

[72] 李厚樸, 邊少鋒, 李海波. 常用等角投影及其解析變換的復(fù)變函數(shù)表示[J]. 測繪科學(xué)技術(shù)學(xué)報, 2012, 29(2): 109-112, 117.

LI Houpu, BIAN Shaofeng, LI Haibo. Expressions of Commonly Used Conformal Projections and Their Analytical Transformations by Complex Numbers[J]. Journal of Geomatics Science and Technology, 2012, 29(2): 109-112, 117.

[73] 李厚樸, 邊少鋒. 高斯投影與墨卡托投影解析變換的復(fù)變函數(shù)表達式[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2009, 34(3): 277-279.

LI Houpu, BIAN Shaofeng. Expressions for Analytical Transformation Between Gauss and Mercator Projections By Complex Numbers[J]. Geomatics and Information Science of Wuhan University, 2009, 34(3): 279-273.

[74] 李勝全, 李厚樸, 邊少鋒. 拉格朗日投影與常用等角投影間解析變換的復(fù)變函數(shù)表示[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2012, 37(11): 1382-1385.

LI Shengquan, LI Houpu, BIAN Shaofeng. Expressions of Analytical Transformations Between Lagrange and the Commonly Used Conformal Projections by Complex Numbers[J]. Geomatics and Information Science of Wuhan University, 2012, 37(11): 1382-1385.

[75] HOTINE M. The Orthomorphic Projection of the Spheroid[J]. Empire Survey Review, 1946, 8(62): 300-311.

[76] SNYDER J P. Calculating Map Projections for the Ellipsoid[J]. The American Cartographer, 1979, 6(1): 67-76.

[78] 李全海. 斜圓柱投影及其對GPS空間坐標(biāo)的轉(zhuǎn)換模型[J]. 同濟大學(xué)學(xué)報, 2002, 30(1): 106-110.

LI Quanhai. Model of Oblique Cylindrical Projection for Transforming GPS Space Coordinates to Plane Coordinates[J]. Journal of Tongji University, 2002, 30(1): 106-110.

[79] 沙月進, 高洪興, 胡伍生. 斜軸圓柱投影方法及其在交通工程中的應(yīng)用研究[J]. 公路交通科技, 2004, 21(11): 20-22.

SHA Yuejin, GAO Hongxing, HU Wusheng. Study on Oblique Cylindrical Projection Method and Its Application in Traffic Engineering[J]. Journal of Highway and Transportation Research and Development, 2004, 21(11): 20-22.

[80] 呂慧玲. 斜軸墨卡托投影方法在鄭西客專中的應(yīng)用研究[J]. 測繪信息與工程, 2009, 34(1): 26-28.

Lü Huiling. Application of Oblique Mercator Projection Method to Zhengxi Special Passenger Transport Line[J]. Journal of Geomatics, 2009, 34(1): 26-28.

[81] 陸鵬程, 林冬偉. 斜軸墨卡托投影模型及其應(yīng)用分析[J]. 鐵道勘察, 2010, 36(4): 26-29.

LU Pengcheng, LIN Dongwei. Model for Oblique Mercator Projections as Well as Analysis on Its Application[J]. Railway Investigation and Surveying, 2010, 36(4): 26-29.

[82] 劉家臣. 高速鐵路精密工程控制測量精度研究[J]. 測繪與空間地理信息, 2011, 34(6): 215-218.

LIU Jiachen. Research on Accuracy of Precise Engineering Control Survey for High-Speed Railway[J]. Geomatics & Spatial Information Technology, 2011, 34(6): 215-218.

[83] 王文慶, 帥明明, 龔俊. 斜軸墨卡托投影在東西走向高速鐵路中的應(yīng)用[J]. 地理空間信息, 2013, 11(4): 152-153.

WANG Wenqing, SHUAI Mingming, GONG Jun. Application of Oblique Mercator Projection in the East-West High-Speed Railway[J]. Geospatial Information, 2013, 11(4): 152-153.

[84] 金立新, 付宏平. 法截面子午線橢球高斯投影理論[M]. 西安: 西安地圖出版社, 2012.

JIN Lixin, FU Hongping. Gaussian Projection Theory Based on Ellipsoid with Normal Section as the Central Meridian[M]. Xi’an: Xi’an Map Publishing House, 2012.

[85] 金立新, 付宏平. 法截面子午線橢球空間幾何理論[M]. 西安: 西安地圖出版社, 2013.

JIN Lixin, FU Hongping. Spatial Geometry Theory Based on Ellipsoid with Normal Section as the Central Meridian[M]. Xi’an: Xi’an Map Publishing House, 2013.

[86] 金立新, 付宏平. 法截面子午線橢球工程應(yīng)用研究[M]. 西安: 西安地圖出版社, 2016.

JIN Lixin, FU Hongping. Engineering Application Based on Ellipsoid with Normal Section as the Central Meridian[M]. Xi’an: Xi’an Map Publishing House, 2016.

[87] 邊少鋒, 劉強, 李忠美, 等. 斜軸變形橢球高斯投影方法[J]. 測繪學(xué)報, 2015, 44(10): 1071-1077. DOI: 10.11947/j.AGCS.2015.20140290.

BIAN Shaofeng, LIU Qiang, LI Zhongmei, et al. An Alteration of Gauss Projection Based on Oblique Deformed Ellipsoid[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(10): 1071-1077. DOI: 10.11947/j.AGCS.2015.20140290.

[88] FOX W A W. Transverse Navigation: An Alternative to the Grid System[J]. The Journal of Navigation, 1949, 2(1): 25-35.

[89] KENDALL E C. Gyro/Grid Navigation[J]. The Journal of Navigation, 1956, 9(4): 429-435.

[90] DYER G C. Polar Navigation-A New Transverse Mercator Technique[J]. The Journal of Navigation, 1971, 24(4): 484-495.

[91] NAUMANN J. Grid Navigation with Polar Stereographic Charts[J]. European Journal of Navigation, 2011, 9(1): 4-8.

[92] SKOPELITI A, TSOULOS L. Choosing A Suitable Projection for Navigation in the Arctic[J]. Marine Geodesy, 2013, 36(2): 234-259.

[93] 李樹軍, 張哲, 李惠雯, 等. 編制北極地區(qū)航海圖有關(guān)問題的探討[J]. 海洋測繪, 2012, 32(1): 58-60.

LI Shujun, ZHANG Zhe, LI Huiwen, et al. Research on Compilation of Nautical Charts of Arctic Regions[J]. Hydrographic Surveying and Charting, 2012, 32(1): 58-60.

[94] 溫朝江, 卞鴻巍, 王榮穎, 等. 極區(qū)極球面投影的可用性及誤差分析[J]. 海軍工程大學(xué)學(xué)報, 2014, 26(3): 42-47.

WEN Chaojiang, BIAN Hongwei, WANG Rongying, et al. Availability and Error Analysis of Polar Stereographic Projection in Polar Navigation[J]. Journal of Naval University of Engineering, 2014, 26(3): 42-47.

[95] 溫朝江, 卞鴻巍, 陳秋, 等. 雙重極球面投影的極區(qū)投影誤差分析[J]. 海洋測繪, 2015, 35(1): 34-37.

WEN Chaojiang, BIAN Hongwei, CHEN Qiu, et al. Error Analysis of Double Polar Stereographic Projection in Polar Regions[J]. Hydrographic Surveying and Charting, 2015, 35(1): 34-37.

[96] 溫朝江, 卞鴻巍, 邊少鋒, 等. 基于等距圓的極球面投影距離量測方法[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2015, 40(11): 1504-1508, 1512.

WEN Chaojiang, BIAN Hongwei, BIAN Shaofeng, et al. A Distance Measuring Method on Polar Stereographic Charts Based on Equidistant Circles[J]. Geomatics and Information Science of Wuhan University, 2015, 40(11): 1504-1508, 1512.

[97] 張志衡, 彭認燦, 董箭, 等. 極地海區(qū)等距離正圓柱投影平面上等角航線的展繪方法[J]. 測繪科學(xué)技術(shù)學(xué)報, 2015, 32(5): 535-538, 544.

ZHANG Zhiheng, PENG Rencan, DONG Jian, et al. The Method for Describing Rhumb Line in Polar Sea on Equidistance Cylindrical Projection Plane[J]. Journal of Geomatics Science and Technology, 2015, 32(5): 535-538, 544.

[98] 李忠美, 李厚樸, 邊少鋒. 極區(qū)橫墨卡托投影非奇異公式及投影變形分析[C]∥第二十四屆海洋測繪綜合性學(xué)術(shù)研討會. 義烏: 中國測繪學(xué)會海洋測繪專業(yè)委員會, 2012: 189-195.

LI Zhongmei, LI Houpu, BIAN Shaofeng. The Nonsingular Formulas and Deformation Analysis for Transverse Mercator Projection in Polar Regions[C]∥Hydrographic Surveying and Charting Symposium of 24th General Meeting. Yiwu: CSGPC, 2012: 189-195.

[99] 張曉平, 邊少鋒, 李忠美. 極區(qū)高斯投影與日晷投影的比較[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2015, 40(5): 667-672.

ZHANG Xiaoping, BIAN Shaofeng, LI Zhongmei. Comparisons Between Gauss and Gnomonic Projections in Polar Regions[J]. Geomatics and Information Science of Wuhan University, 2015, 40(5): 667-672.

[100] 邊少鋒, 李忠美, 李厚樸. 極區(qū)非奇異高斯投影復(fù)變函數(shù)表示[J]. 測繪學(xué)報, 2014, 43(4): 348-352, 359. DOI: 10.13485/j.cnki.11-2089.2014.0052.

BIAN Shaofeng, LI Zhongmei, LI Houpu. The Non-Singular Formula of Gauss Projection in Polar Regions by Complex Numbers[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(4): 348-352, 359. DOI: 10.13485/j.cnki.11-2089.2014.0052.

[101] 李忠美, 邊少鋒, 金立新, 等. 極區(qū)不分帶高斯投影的正反解表達式[J]. 測繪學(xué)報, 2017, 46(6): 780-788. DOI: 10.11947/j.AGCS.2017.20170009.

LI Zhongmei, BIAN Shaofeng, JIN Lixin, et al. Forward and Inverse Expressions of Polar Gauss Projection without Zoning Limitations[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(6): 780-788. DOI: 10.11947/j.AGCS.2017.20170009.

[102] 李厚樸, 邊少鋒, 劉強, 等. 常用極區(qū)海圖投影直接變換的閉合公式[J]. 海洋測繪, 2017, 37(2): 32-34, 38.

LI Houpu, BIAN Shaofeng, LIU Qiang, et al. Closed Formulae of Direct Transformations Between Commonly Used Polar Chart Projections[J]. Hydrographic Surveying and Charting, 2017, 37(2): 32-34, 38.

[103] 陳成. 極區(qū)海圖投影及其變換研究[D]. 武漢: 海軍工程大學(xué), 2015.

CHEN Cheng. The Research on Polar Chart Projection with Its Transformation[D]. Wuhan: Naval University of Engineering, 2015.

[104] COLVOCORESSES A P. Space Oblique Mercator[J]. Photogrammetric Engineering and Remote Sensing, 1974, 40(8): 921-926.

[105] SNYDER J P. The Space Oblique Mercator Projection[J]. Photogrammetric Engineering and Remote Sensing, 1978, 44(5): 585-596.

[106] REN Liucheng, CLARKE K C, ZHOU Chenghu, et al. Geometric Rectification of Satellite Imagery with Minimal Ground Control Using Space Oblique Mercator Projection Theory[J]. Cartography and Geographic Information Science, 2010, 37(4): 261-272.

(責(zé)任編輯：陳品馨)

Research Progress in Mathematical Analysis of Map Projection by Computer Algebra

BIAN Shaofeng,LI Houpu,LI Zhongmei

Department of Navigation, Naval University of Engineering, Wuhan 430033, China

Map projection is an important component of modern cartography, and involves many fussy mathematical analysis processes, such as the power series expansions of elliptical functions, differential of complex and implicit functions, elliptical integral and the operation of complex numbers. The derivation of these problems by hand not only consumes much time and energy but also makes mistake easily, and sometimes can not be realized at all because of the impossible complexity. The research achievements in mathematical analysis of map projection by computer algebra are systematically reviewed in five aspects, i.e., the symbolic expressions of forward and inverse solution of ellipsoidal latitudes, the direct transformations between map projections with different distortion properties, expressions of Gauss projection by complex function, mathematical analysis of oblique Mercator projection, polar chart projection with its transformation. Main problems that need to be further solved in this research field are analyzed. It will be helpful to promote the development of map projection.

map projection; projection transformation; Gauss projection; polar projection; computer algebra

The National Natural Science Foundation of China(Nos. 41631072；41571441；41604010)

BIAN Shaofeng(1961—),male,professor,majors in map projection,geodesy and satellite navigation.

LI Houpu

邊少鋒，李厚樸，李忠美.地圖投影計算機代數(shù)分析研究進展[J].測繪學(xué)報，2017,46(10)：1557-1569.

10.11947/j.AGCS.2017.20170396.

BIAN Shaofeng,LI Houpu,LI Zhongmei.Research Progress in Mathematical Analysis of Map Projection by Computer Algebra[J]. Acta Geodaetica et Cartographica Sinica,2017,46(10):1557-1569. DOI:10.11947/j.AGCS.2017.20170396.

P282

A

1001-1595(2017)10-1557-13

國家自然科學(xué)基金(41631072；41571441；41604010)

2017-07-09

修回日期： 2017-09-05

邊少鋒(1961—)，男，教授，研究方向為地圖投影、大地測量和衛(wèi)星導(dǎo)航。

E-mail： sfbian@sina.com

李厚樸

E-mail： lihoupu 1985@126.com