高精度高斯投影正反算公式及其應(yīng)用
關(guān)鍵詞:高斯投影;大地坐標(biāo);高斯平面坐標(biāo);牛頓迭代法
一、引言
高斯投影是正形投影的一種,具有投影長度比與方向無關(guān),投影前后角度不變等優(yōu)點,廣泛應(yīng)用于各種工程項目中。地球橢球面是一個不可展曲面,而地圖是一個無裂縫無重疊的平面,高斯投影采用泰勒級數(shù)展開的方法將地球橢球面轉(zhuǎn)換到地圖平面,但由于展開式為無窮級數(shù),且高次項系數(shù)極為復(fù)雜,所以現(xiàn)在一般選取低次項系數(shù)做近似替換。為了克服投影變形,通常采取分帶投影的方法,我國規(guī)定按經(jīng)差6°和3°進行投影分帶[1]。隨著工程項目跨度逐步增大,東西跨度早已超過6°帶的投影范圍,這就意味著需要新的方法來實現(xiàn)大跨度高精度的高斯投影。
本文根據(jù)高斯投影的一般原理,將投影正算公式展開至第9項。同時,對于高斯投影反算,現(xiàn)在多采用迭代的方法先求得底點緯度Bf,再將Bf代入反算公式求得大地經(jīng)緯度[2]。由于反算公式高階項系數(shù)推導(dǎo)困難,反算精度受限。本文根據(jù)多元非線性函數(shù)的牛頓迭代算法[3],推導(dǎo)出了由高斯平面坐標(biāo)(x、y)直接迭代求得大地坐標(biāo)(B、L)的迭代公式,有效的避開了反算公式高階項系數(shù)缺失的問題。本文分別以12°、24°為帶寬,計算分析投影精度。結(jié)果表明,投影精度有較大提高,投影帶寬度得到相應(yīng)的擴展。
二、高斯投影坐標(biāo)正算
高斯投影正算是已知大地坐標(biāo)(B,L),求高斯平面坐標(biāo)(x,y)。高斯投影正算公式為
(1)
式中,x、y為高斯平面坐標(biāo),x0為子午線弧長,l為經(jīng)差, ( )是緯度的函數(shù)。計算公式為
(2)
(3)
(4)
式中, 為橢球長半徑, 為第一偏心率, 為第二偏心率, 為該投影帶中央子午線經(jīng)度, 為大地緯度, 為投影點到 的大地經(jīng)度差, ,以弧度為單位。 , , 。
三、高斯投影坐標(biāo)反算的牛頓迭代算法
高斯投影坐標(biāo)反算是已知投影帶中央子午線經(jīng)度和高斯平面坐標(biāo)(x,y),反解大地經(jīng)緯度(B,L)。傳統(tǒng)反算方法為先迭代算出底點緯度Bf,然后代入反算公式,進而求得大地坐標(biāo)(B,L)[4-6]。雖然運用了迭代的思想,可是仍然受到反算公式高次項缺失的制約,導(dǎo)致反算的精度不理想。本文采用多元非線性牛頓迭代法,無需解算底點緯度Bf,直接由高斯平面坐標(biāo)(x,y)迭代算出大地坐標(biāo)(B,L)。若計算機性能許可,可以達到任意精度。
由(1)式得
(5)
令
(6)
則
(7)
得如下迭代格式:
(8)
當(dāng)同時滿足
(9)
時終止迭代。其中 根據(jù)具體要求而定。
四、精度分析
根據(jù)以上所推公式,用C#語言編寫12°、24°帶高斯投影正反算程序( 取1×10-16)。在緯度為45°處取點,其大地坐標(biāo)為(B0,L0)。首先將該點歸算到高斯投影平面上,求得它的高斯平面坐標(biāo),再用迭代反算求得該點的大地坐標(biāo)(B1,L1)。(B0,L0)與(B1,L1)相減,即可得到投影誤差(△B,△L)。其結(jié)果見表1和表2。
表1 采用12°帶高斯投影坐標(biāo)正反算結(jié)果對照表
B(°′″) L(°′″)
45 00 00
11 59 00
45 00
11 59
0 00
0 00
由以上結(jié)果可知,當(dāng)帶寬為12°時,投影帶邊沿投影精度可達4×10-9s。當(dāng)帶寬擴展為24°時,在投影帶邊沿投影精度達1×10-6s。
五、結(jié)束語
本文將高斯投影正算公式展開至第9項,同時采用多元非線性函數(shù)的牛頓迭代法進行投影反算。通過實驗分析可知,采用該方法進行投影計算,在24°投影帶的邊沿處,正反算經(jīng)緯度差值小于0.000001s,投影精度較高。
目前我國正在修建的以高速鐵路為代表的大型工程對測量的精度要求很高,所采用的高斯投影獨立坐標(biāo)系每個投影帶可控制范圍很小,給施工測量帶來很多不便。通過本文所述的方法進行投影計算,能在很大程度上拓寬每個投影帶的投影寬度,對這種大跨度、高精度的工程項目的投影計算有一定的借鑒意義。
參考文獻:
[1] 孔祥元,郭際明,劉宗泉.大地測量學(xué)基礎(chǔ)[M].武漢:武漢大學(xué)出版社,2001.
[2] 孔祥元,梅是義.控制測量學(xué)[M].武漢:武漢大學(xué)出版社,1996.
[3] 王尊正.數(shù)值分析基礎(chǔ)教程[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,1992.
[4] 楊建華,楊志強,王騰軍.高斯投影反算中求地點緯度值的牛頓迭代法[J].西安科技大學(xué)學(xué)報,2005,25(1):57-59.
欄目分類
- 新媒體語境下地方文化故事化傳播的策略研究
- 博物館數(shù)字化傳播對文化遺產(chǎn)保護的影響
- 數(shù)字化轉(zhuǎn)型背景下文化創(chuàng)意產(chǎn)業(yè)人才供需匹配研究
- 公共圖書館詩詞閱讀推廣研究
- 基于“工匠精神”培養(yǎng)的高職院校圖書館閱讀推廣策略研究
- 剪紙和文學(xué)在文旅融合中的發(fā)展研究
- “書磚”古籍的修復(fù)研究——以《及第征祥錄》為例
- 基于OBE教學(xué)理念的工業(yè)物聯(lián)網(wǎng)課程人才培養(yǎng)模式研究
- 紙質(zhì)在文化創(chuàng)意產(chǎn)品設(shè)計中的應(yīng)用
- 數(shù)字經(jīng)濟驅(qū)動下的制造業(yè)轉(zhuǎn)型路徑與企業(yè)績效關(guān)系研究
- 2023JCR影響因子正式公布!
- 國內(nèi)核心期刊分級情況概覽及說明!本篇適用人群:需要發(fā)南核、北核、CSCD、科核、AMI、SCD、RCCSE期刊的學(xué)者
- 我用了一個很復(fù)雜的圖,幫你們解釋下“23版最新北大核心目錄有效期問題”。
- 重磅!CSSCI來源期刊(2023-2024版)最新期刊目錄看點分析!全網(wǎng)首發(fā)!
- CSSCI官方早就公布了最新南核目錄,有心的人已經(jīng)拿到并且投入使用!附南核目錄新增期刊!
- 北大核心期刊目錄換屆,我們應(yīng)該熟知的10個知識點。
- 注意,最新期刊論文格式標(biāo)準(zhǔn)已發(fā)布,論文寫作規(guī)則發(fā)生重大變化!文字版GB/T 7713.2—2022 學(xué)術(shù)論文編寫規(guī)則
- 盤點那些評職稱超管用的資源,1,3和5已經(jīng)“絕種”了
- 職稱話題| 為什么黨校更認可省市級黨報?是否有什么說據(jù)?還有哪些機構(gòu)認可黨報?
- 《農(nóng)業(yè)經(jīng)濟》論文投稿解析,難度指數(shù)四顆星,附好發(fā)選題!