尹　暉,李小祥,甘喆淵

(1.武漢大學(xué) 測繪學(xué)院，湖北 武漢 430079;2.中國電力科學(xué)研究院高壓研究所,湖北 武漢 430074)

?

橢球變換法建立地方獨立坐標(biāo)系的變形研究

尹暉1,李小祥1,甘喆淵2

(1.武漢大學(xué) 測繪學(xué)院，湖北 武漢 430079;2.中國電力科學(xué)研究院高壓研究所,湖北 武漢 430074)

隨著我國城市基礎(chǔ)設(shè)施建設(shè)的迅速發(fā)展，對大比例地圖測繪的精度要求越來越高。而采用已有橢球參數(shù)來進(jìn)行高斯平面投影卻不能滿足具體項目的高精度要求，必須使用新的橢球參數(shù)來建立地方獨立坐標(biāo)系統(tǒng)。一般的做法是平移中央子午線和改變投影高程面來滿足城市測量中要求邊長變形小于2.5 cm/km的限值。過去大部分的研究主要集中在從理論上探討如何改變橢球參數(shù)來建立獨立坐標(biāo)系，缺少應(yīng)用實際算例對橢球參數(shù)選擇的對比和適用范圍分析。本文針對3種不同方法：橢球膨脹法、橢球平移法和橢球變形法進(jìn)行實例分析，選擇不同橢球參數(shù)進(jìn)行計算得到的新坐標(biāo)系統(tǒng)，并進(jìn)行變形分析，最后給出不同方法的優(yōu)缺點及適用范圍，為不同工程實踐提供適宜的參考方案。

1高斯投影變形分析

由于地球是一個近似的大地橢球體，在地圖投影的過程中就不可避免地產(chǎn)生投影變形。投影變形是指將橢球面上的元素(距離、角度、圖形)投影到平面上產(chǎn)生的長度、角度和面積的變形。高斯投影是正形投影，角度不產(chǎn)生變形[1]，投影后的中央子午線是一條沒有變形的直線。在中央子午線附近區(qū)域變形很小，但自中央子午線向投影帶邊緣，變形逐漸增加，變形最大處在投影帶內(nèi)赤道的兩端[2]。

通常用高斯投影長度比m來表示高斯投影平面上的微分長度與地球橢球面上相應(yīng)長度之比。用m-1表示投影變形的變形量大小。

設(shè)高斯平面上的某點K，用其大地坐標(biāo)(B,L)和平面坐標(biāo)(x,y)來表示的高斯投影長度比的計算公式[1]分別為式(1)和式(2)。

(1)

(2)

其中：u2=e′2cos2B，B,L為K點的經(jīng)緯度，e′為橢球體的第二偏心率，y為K點的橫坐標(biāo)，R為該點處橢球的平均曲率半徑，l為K點的經(jīng)度差，l=L-L0，L0為選擇的中央子午線經(jīng)度。根據(jù)參考文獻(xiàn)[3]可知，在l≤3.5°時，高斯投影的長度比精度優(yōu)于1×10-7。表1給出了不同經(jīng)差l和緯度B下的投影變形量的數(shù)值。

表1　不同經(jīng)差l和緯度B的變形量　mm/km

從表1計算的投影變形量結(jié)果可見：由于我國地理位置位于緯度3°51′N至53°33′N的區(qū)域，在采用3°分帶，最大經(jīng)差1°30′的情況下，靠近投影帶邊緣的變形達(dá)到了200~350 mm/km，這種投影產(chǎn)生的長度變形遠(yuǎn)大于目前城市測量中要求長度變形小于2.5 cm/km的限值規(guī)定；從式(2)計算可以得到在橫坐標(biāo)y大于45 km之后，長度變形將超過1/40 000，大于工程測量中所要求的小于2.5 cm/km的長度變形誤差，需要對橢球參數(shù)和投影參數(shù)進(jìn)行控制。

2橢球變換法建立地方獨立坐標(biāo)系

2.1建立地方獨立坐標(biāo)系的常規(guī)方法

建立地方獨立坐標(biāo)系的目的就是為了控制長度變形，且保持角度不發(fā)生改變。即實際地面測得的距離轉(zhuǎn)換到平面坐標(biāo)系中的兩者距離之差要滿足實際工程的精度要求。

把地面距離觀測值歸算到橢球面需要歸算方向(方位角)改正和距離(長度)改正。由于地面測量距離的觀測值是兩控制點之間的直線長度，空間直線的長度與點的鉛垂線沒有關(guān)系，故可以直接沿端點的法線歸算致橢球面上，產(chǎn)生的變形為

(3)

要求投影的長度變形小于等于2.5cm/km，在常規(guī)方法中只控制投影參數(shù)：改變投影面的中央子午線；改變投影面大地高；同時改變投影面的中央子午線與投影面大地高。距離控制時往往并未控制投影方向的距離改化，引起角度畸變，理論上不夠嚴(yán)密。

對于改變投影面的中央子午線和大地高進(jìn)行高斯3°帶投影仍然變形超限的項目測區(qū)(如長距離道路建設(shè)、管道建設(shè)、城市級區(qū)域工程建設(shè)等)，就無法采用這種方法，需要采用橢球變換的方法建立地方獨立坐標(biāo)系。

2.2幾種常見的橢球變換方法

2.2.1橢球膨脹法

橢球膨脹法的基本思想是保持橢球中心和橢球定向在膨脹前后不發(fā)生變動，對橢球進(jìn)行縮放，縮放后的橢球面與獨立坐標(biāo)系所選的平面相切[4-5]。新橢球的長半軸會發(fā)生改變，橢球原點位置不發(fā)生改變(見圖1)。

圖1　橢球膨脹法變換示意圖

在實際的工程使用中一般改變長半軸有3種方法，產(chǎn)生3種不同的結(jié)果，但都是改變橢球長半軸a而不改變橢球第一偏心率e。橢球參數(shù)變化有如下3種：

1)新橢球的長半軸增加的量為項目區(qū)域現(xiàn)選擇的投影面大地高H。

(4)

2)新橢球的長半軸增加的量為原投影面沿法線方向膨脹到現(xiàn)選擇的投影面的增加量。

(5)

3)新橢球的長半軸增加的量為現(xiàn)選擇的投影面的平均曲率半徑。

(6)

式中：H0為現(xiàn)選擇的投影面的大地高，e為橢球第一偏心率(保持不變)，B0為目標(biāo)區(qū)域基準(zhǔn)緯度，L0為目標(biāo)區(qū)域基準(zhǔn)經(jīng)度。

采用橢球膨脹法得到的控制點大地坐標(biāo)換算到獨立坐標(biāo)系的計算公式為

(7)

式中：

(8)

式中：dBx,dLx,dHx為待定控制點的大地坐標(biāo)變化量;Nx為待定點在原橢球參數(shù)下的子午圈曲率半徑;Mx為待定點在原橢球參數(shù)下的卯酉圈曲率半徑;Lx,Bx,Hx為待定控制點在原坐標(biāo)系下的大地坐標(biāo)。

經(jīng)過橢球膨脹法的變換，待定控制點經(jīng)度L不發(fā)生改變，而緯度B和大地高H發(fā)生改變。

2.2.2橢球平移法

相對于橢球膨脹法改變橢球長半軸，橢球平移法是通過移動橢球中心位置來達(dá)到效果的。其基本思想是將原橢球沿著項目區(qū)域獨立坐標(biāo)系的基準(zhǔn)點所在的法線方向進(jìn)行平移，使得移動后橢球面平移到現(xiàn)選擇的投影面(見圖2)。橢球的其余參數(shù)不發(fā)生改變da=0，de2=0。

圖2　橢球平移法變換示意圖

橢球參數(shù)變化相應(yīng)的公式為

(9)

式中：(dX0，dY0，dZ0)表示橢球中心的變化量;B0,L0,M0,N0,H0定義同橢球膨脹法。

采用橢球平移法得到的控制點大地坐標(biāo)換算到獨立坐標(biāo)系的計算公式為

(10)

式中：dBx,dLx,dHx為待定控制點的大地坐標(biāo)變化量;Mx,Nx,Hx,Lx,Bx,Hx定義同橢球膨脹法。

可見經(jīng)過橢球平移法，點的經(jīng)緯度和大地高都發(fā)生改變。

2.2.3橢球變形法

橢球變形的思想是保持橢球中心不發(fā)生改變，改變橢球的扁率α。具體做法是：先將橢球沿基準(zhǔn)點的法線膨脹到所定義的參考面，變化扁率α,使得基準(zhǔn)點處的法線方向與變換前后重合(見圖3)。

圖3　橢球變形法變換示意圖

這樣變換的結(jié)果是保證了基準(zhǔn)點的經(jīng)度和緯度不發(fā)生改變。橢球參數(shù)相應(yīng)的公式為

(11)

式中：e′2為變換后的橢球第一偏心率的平方;N0為基準(zhǔn)點在原橢球參數(shù)下的子午圈曲率半徑;da為橢球長半軸改變量;B0,N0,H0,a,e定義同橢球膨脹法。

采用橢球變形法得到的控制點大地坐標(biāo)換算到獨立坐標(biāo)系的計算公式為

(12)

式中：dBx,dLx,dHx為待定控制點的大地坐標(biāo)變化量;Mx,Nx,Hx,Lx,Bx定義同橢球膨脹法;α為原橢球扁率;dα為橢球扁率變化量;α′為新橢球扁率。

橢球平移法的變換前后可以發(fā)現(xiàn)待定控制點的經(jīng)度L不發(fā)生改變，而緯度B和大地高H發(fā)生改變。

2.3算法流程

首先，不管采用哪種橢球變換法，需要確定目標(biāo)區(qū)域基準(zhǔn)量：基準(zhǔn)緯度B0、基準(zhǔn)經(jīng)度L0和現(xiàn)選擇的投影面大地高H0，B0、L0和H0可以采用項目區(qū)域已有資料中的基準(zhǔn)點經(jīng)緯度和目標(biāo)投影面大地高來確定，也采用待定控制點經(jīng)緯度和大地高平均值的方法計算得到。

其次是確定橢球變換參數(shù)：①橢球膨脹法是確定橢球長半軸的變化量da；②橢球平移法是確定橢球中心變化量(dX0，dY0，dZ0)；③橢球變形是確定橢球的長半軸的變化量da、新橢球偏心率e′和新橢球扁率α′。

最后是根據(jù)不同的橢球變換方法待定控制點的大地坐標(biāo)改變量dBx,dLx,dHx，從而獲得新橢球下的控制點大地坐標(biāo)，如需要再投影成高斯平面坐標(biāo)。

HDS2003和TGO是GPS數(shù)據(jù)處理的常用軟件，功能強(qiáng)且應(yīng)用面廣。但這兩種軟件只采用了橢球膨脹法來確定橢球變換參數(shù)[6]，缺乏3種方法的效果對比。本文基于C#重新設(shè)計了建立地方獨立坐標(biāo)系計算軟件，如圖4所示。

圖4　橢球變換法建立獨立坐標(biāo)系軟件界面

3數(shù)據(jù)分析

測區(qū)為南方某市級區(qū)域，面積約1萬km2，平均緯度23°30′，平均經(jīng)度114°30′?，F(xiàn)有WGS-84坐標(biāo)系起算點68個。橢球變換法建立地方獨立坐標(biāo)系可以采用平均經(jīng)緯度和大地高作為基準(zhǔn)量，也可以采用控制點的平均經(jīng)緯度和大地高作為基準(zhǔn)量[7]。在軟件設(shè)計中既設(shè)計了以區(qū)域平均經(jīng)緯度和大地高作為基準(zhǔn)量也設(shè)計了使用控制點的平均經(jīng)緯度和大地高作為基準(zhǔn)量進(jìn)行計算。GPS控制點一般在山頂或者樓頂，其大地高比項目區(qū)域平均大地高偏高，故采用區(qū)域平均經(jīng)緯度和大地高作為基準(zhǔn)量，來進(jìn)行數(shù)據(jù)分析[8]。

根據(jù)不同橢球變換法計算得到的橢球參數(shù)成果如表2所示。

表2　不同方法橢球參數(shù)　m

限于篇幅選擇其中3個控制點(其中3號控制點為偏離中央子午線最遠(yuǎn)的點)，計算結(jié)果如表3所示(橢球膨脹法展示第一種參數(shù)下的結(jié)果)。

表3　不同方法計算得到的大地坐標(biāo)

計算得到1、2、3號控制點和全部控制點處平均的微分長度投影變形如表4所示。

表4　3種方法每千米邊長變形對比　cm/km

數(shù)據(jù)分析顯示：①3種橢球變換法都對長度投影變形進(jìn)行了很好的控制，滿足在整個市級區(qū)域2.5 cm/km邊長變形要求。橢球平移法變形較為平均，且最大變形比另外兩種方法小。②在程序設(shè)計過程中因橢球平移法變換后橢球形狀參數(shù)不發(fā)生改變，大地坐標(biāo)投影成高斯平面坐標(biāo)可以采用現(xiàn)有電算公式或軟件，復(fù)雜性大為降低。③橢球變換法與常規(guī)方法相比，在滿足城市測量2.5 cm/km的變形要求下，其控制面積達(dá)到了1萬km2，大大拓展了獨立坐標(biāo)系的使用范圍。

4結(jié)論

高斯投影變形是高斯投影的關(guān)注熱點，尤其在大比例尺地形圖應(yīng)用中，由于精度要求較高，采用原有橢球參數(shù)進(jìn)行投影就達(dá)不到實際工程的高精度要求。

1)常規(guī)方法控制投影變形控制的面積較小，計算復(fù)雜，且不具有理論嚴(yán)密性。橢球變換法建立獨立坐標(biāo)系控制面積大為增加，計算簡單，具有理論嚴(yán)密性。

2)根據(jù)《建立相對獨立的平面坐標(biāo)系統(tǒng)管理辦法》的規(guī)定：為避免重復(fù)建設(shè)，一個城市只能建立一個相對獨立的平面坐標(biāo)系統(tǒng)。因此對于基準(zhǔn)經(jīng)緯度和投影面的選擇顯得尤為重要。對基準(zhǔn)的確定可以采用平均的方法，利用GPS技術(shù)快速地獲取項目區(qū)域均勻分布的多個控制點，加以平均來獲得基準(zhǔn)參數(shù)，也可以根據(jù)已有資料來確定。計算得到新的橢球參數(shù)，再把新的控制點大地坐標(biāo)投影到平面坐標(biāo)應(yīng)使用橢球變換的中央子午線，才能更好地控制投影變形。

3)在算法設(shè)計正對3種橢球變換的方法進(jìn)行了程序設(shè)計。經(jīng)過數(shù)據(jù)分析橢球平移法對投影變形的控制效果最好，也在算法上最簡單。特別是進(jìn)行高斯投影時，無需改變橢球參數(shù)，實現(xiàn)了快速計算。

4)根據(jù)以上分析建議使用橢球變換法中的橢球平移法來建立地方獨立坐標(biāo)系。

參考文獻(xiàn)：

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

[2]孔祥元，郭際明，劉宗泉.大地測量學(xué)基礎(chǔ)[M].武漢：武漢大學(xué)出版社，2010.

[3]趙俊生，劉雁春，王克平，等.關(guān)于高斯投影長度變形的探討[J].海洋測繪，2007，27(3)：9-11.

[4]李世安，劉經(jīng)南，施闖.應(yīng)用GPS建立區(qū)域獨立坐標(biāo)系中橢球變換的研究[J].武漢大學(xué)學(xué)報(信息科學(xué)版)，2005，30(10)：888-891.

[5]魏長明，刑建國，劉小剛.利用橢球變換建立獨立坐標(biāo)系[J].礦山測量，2014(1)：30-32.

[6]張宏偉，康世英.基于TGO的獨立坐標(biāo)系與國家坐標(biāo)系轉(zhuǎn)換方法[J].地礦測繪，2010(2)：12-15.

[7]邱云峰，倪津，楊映泉，等.基于CGCS2000大地坐標(biāo)系基準(zhǔn)的建立與應(yīng)用[J].測繪工程，2013，22(4)：42-44.

[8]張顯云，杜宇，張俊.GPS技術(shù)建立地方獨立坐標(biāo)系的方法分析[J].礦冶工程，2012(10)：23-25.

[責(zé)任編輯：劉文霞]

摘要：在分析高斯投影變形的基礎(chǔ)上，針對工程測量中的投影變形精度要求，闡述了3種橢球變換法建立地方獨立坐標(biāo)系的方法，給出具體的計算公式和參數(shù)選取原則。對3種橢球變換法進(jìn)行軟件設(shè)計，簡化計算的工作量。并詳細(xì)介紹如何求解這些參數(shù)和控制點坐標(biāo)值，結(jié)合實例數(shù)據(jù)進(jìn)行分析，得到不同方法建立地方獨立坐標(biāo)系的優(yōu)缺點和適用范圍。

關(guān)鍵詞：高斯投影；橢球變換；橢球大地測量；獨立坐標(biāo)系

Deformation of local independent coordinate established with ellipsoid transform method YIN Hui1,LI Xiaoxiang1,GAN Zheyuan2

(1.School of Geodesy and Geomatics,Wuhan University,Wuhan 430079,China;2.High Voltage Department,China Electric Power Research Institute,Wuhan 430074,China)

Abstract：Based on the deformation of Gauss-Kruger projection and the requirement of projection deformation,this paper presents three different ellipsoid-transform metods to establish the independent coordinate system,specific formula and parameter selection principles.A new software is designed to avoid the shortages of current program and the calculation has been simplified.The details to solve these parameters and control point coordinates value date are analyzed with examples,and some meaningful conclusions based on different methods are drawn to show the strengths and weaknesses as well as the scope of local independent coordinate system.

Key words：Gauss-Kruger projection;ellipsoid transform;geodetic ellipsoid;independent coordinate system

作者簡介：尹暉(1962-)，女，教授，博士.

基金項目：國家自然科學(xué)基金資助項目(51077105)；國家電網(wǎng)公司科技攻關(guān)團(tuán)隊資助項目(SG11013)；大地測量與地球動力學(xué)國家重點實驗室開放基金資助項目(SKLGED2013-3-6-E)

收稿日期：2014-12-01

中圖分類號：P226+.3

文獻(xiàn)標(biāo)識碼：A

文章編號：1006-7949(2016)02-0001-05