楊　志， 覃　輝， 文鴻雁， 袁明月

(1.廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室， 廣西 桂林　541004；　2.桂林理工大學(xué) 測(cè)繪地理信息學(xué)院， 廣西 桂林　541004；　3.桂林理工大學(xué) 廣西礦冶與環(huán)境科學(xué)實(shí)驗(yàn)中心， 廣西 桂林　541006)

?

基于投影變形大小的投影帶和投影面選擇

楊志1，2，3， 覃輝1，2，3， 文鴻雁1，2，3， 袁明月1，2，3

(1.廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室， 廣西 桂林541004；2.桂林理工大學(xué) 測(cè)繪地理信息學(xué)院， 廣西 桂林541004；3.桂林理工大學(xué) 廣西礦冶與環(huán)境科學(xué)實(shí)驗(yàn)中心， 廣西 桂林541006)

[摘要]我國(guó)國(guó)家坐標(biāo)系統(tǒng)普遍采用高斯-克呂格正射投影( 簡(jiǎn)稱“高斯投影”)，即先由地面投影到參考橢球面，再由參考橢球面投影到高斯平面。經(jīng)過(guò)以上投影變換，定會(huì)引起一定的投影變形，變形的大小隨投影點(diǎn)與大地水準(zhǔn)面的距離不同而變化。針對(duì)該問(wèn)題，工程測(cè)量中通常采用的處理方法是：a.采用任意分帶形式，將中央子午線移至能控制投影變形的范圍。b.根據(jù)實(shí)際需要，重新設(shè)定投影面高程。c.選定一個(gè)抵償面。是先找到最佳中央子午線，然后計(jì)算抵償面高程，得出最佳的投影帶和投影面，通過(guò)實(shí)例計(jì)算對(duì)比，能得到最佳的投影帶和投影面，在一定范圍內(nèi)，橫坐標(biāo)跨度最大，同時(shí)也能夠準(zhǔn)確的確定大地測(cè)量所需要的平面控制網(wǎng)的坐標(biāo)系。

[關(guān)鍵詞]高斯-克呂格正射投影； 投影變形； 中央子午線； 投影面高程； 抵償面

1概述

對(duì)于工程測(cè)量，尤其是在城市測(cè)量，不僅有測(cè)量大比例尺地形圖的任務(wù)，同時(shí)還要滿足各類建筑工程和市政建筑施工放樣等工作的要求。要建立經(jīng)濟(jì)合理的，并且要能夠滿足各種目的的工程平面控制網(wǎng)坐標(biāo)系，在測(cè)繪行業(yè)是一個(gè)重要的研究方向[1]。

隨著測(cè)繪技術(shù)的發(fā)展，特別是GTP定位技術(shù)，其作用不可小覷。在我國(guó)道路建設(shè)中，尤其是在實(shí)際測(cè)量鐵路系統(tǒng)中，GTP技術(shù)的應(yīng)用最為廣泛。隨著最近幾年來(lái)高速鐵路的迅猛發(fā)展，對(duì)其勘測(cè)、設(shè)計(jì)變得越來(lái)越重要，為降低或者解除高差和高斯投影導(dǎo)致的邊長(zhǎng)變形對(duì)工程建設(shè)施工的影響，常常需要建立適宜的獨(dú)立坐標(biāo)系來(lái)滿足工程的要求，所以就需要向鐵路獨(dú)立坐標(biāo)系轉(zhuǎn)換[2]。但因?yàn)榈缆窚y(cè)量有其自身的特點(diǎn)，通常會(huì)出現(xiàn)測(cè)區(qū)遠(yuǎn)離中央子午線或者測(cè)區(qū)平均高程偏大，從而造成長(zhǎng)度變形增大，很難達(dá)到工程實(shí)踐精度的嚴(yán)格要求。為此，需要建立科學(xué)合理的與本地區(qū)相契合的獨(dú)立坐標(biāo)系。通常而言，鐵路獨(dú)立坐標(biāo)系建立的辦法要分兩步走，首先通過(guò)使用國(guó)家坐標(biāo)系已知點(diǎn)束縛GTP網(wǎng)；隨后在國(guó)家坐標(biāo)系橢球基本參數(shù)的基礎(chǔ)上，運(yùn)用變動(dòng)中央子午線和抵償投影面的方法來(lái)建立鐵路獨(dú)立坐標(biāo)系[3]。

2投影變形的基本概念與精度要求

2.1投影變形概念

投影變形是指將地面觀測(cè)的實(shí)際長(zhǎng)度歸算到參考橢球面后，再將橢球面上的實(shí)際長(zhǎng)度歸算到高斯投影面上的變形值。投影面和投影帶的選擇對(duì)于平面控制網(wǎng)測(cè)量中，主要是解決長(zhǎng)度變形的問(wèn)題[4]。

2.2投影變形的原因

投影變形主要由以下原因所引起[5]：

① 實(shí)際測(cè)量的邊長(zhǎng)歸算到參考橢球面上的變形，其值為ΔS1，據(jù)相關(guān)公式計(jì)算，在不同高程面上對(duì)應(yīng)的變形值見(jiàn)表1。從表1可知：ΔS1是負(fù)值，表明將地面實(shí)際測(cè)量的邊長(zhǎng)歸算到參考橢球面上，總是減小的；歸算邊高出參考橢球面的平均高程Hm與ΔS1的絕對(duì)值成正比，ΔS1隨Hm的增大而增大。

② 將參考橢球面上的邊長(zhǎng)歸算到高斯投影面上的變形，其值為ΔS2，根據(jù)相關(guān)公式，可以計(jì)算每公里長(zhǎng)度的投影變形值和相對(duì)投影變形值，見(jiàn)表1。

從表1可見(jiàn)：ΔS2的值總是正數(shù)，由此可知將橢球面上的長(zhǎng)度投影到高斯投影面上[6]，投影變形值總是增大的；ΔS2值隨著歸算邊兩端橫坐標(biāo)的平均值ym的平方增大而增大，由此可知，離中央子午線距離越遠(yuǎn)，其變形值越大。

表1 不同高程面上距離的相對(duì)變形與每公里投影變形及相對(duì)投影變形值Table1 Therelativedeforamtionofdistanceindifferentheightsurface&theprojectiondistortionvaluesandrelativeprojectiondistortionvaluesperkilometer不同高程面上的相對(duì)變形每公里長(zhǎng)度投影變形及相對(duì)投影變形值Hm/mΔS1/mmΔS1/Sy/mΔS2/mmΔS2/S010-1.61/637000101.21/81000020-3.11/318500204.91/12000030-4.71/2120003011.11/9000040-6.31/1590004019.71/15000050-7.81/1274005030.71/3200060-9.41/1060006044.31/2200070-11.01/910007060.31/1650080-12.61/790008078.71/1270090-14.11/700009099.61/10000100-15.71/63700100133.01/8000160-25.11/39000———1000-1571/6370— — —

2.3投影變形的精度要求

在各類工程建設(shè)中，施工放樣工作能否順利進(jìn)行，就必須要求控制點(diǎn)坐標(biāo)反算的邊長(zhǎng)與實(shí)地測(cè)量的邊長(zhǎng)應(yīng)該相等，也就是說(shuō)由上述兩項(xiàng)歸算投影改正而引起的長(zhǎng)度變形或改正數(shù)，應(yīng)該小于施工放樣的精度要求[7]。即每公里變形值應(yīng)小于25 mm。

3長(zhǎng)度投影變形影響因素

3.1實(shí)際測(cè)量真實(shí)長(zhǎng)度歸化至國(guó)家統(tǒng)一的參考橢球面

測(cè)量的實(shí)際長(zhǎng)度歸化到國(guó)家統(tǒng)一的參考橢球面就會(huì)產(chǎn)生高差投影變形[8，9]?？砂词?1)計(jì)算：

(1)

式中：D為實(shí)測(cè)水平距離；RA為長(zhǎng)度所在方向參考橢球面法截弧的曲率半徑；Hm為觀測(cè)邊的平均大地高。

3.2橢球面上的長(zhǎng)度投影到高斯平面

(2)

式中：S為參考橢球面的長(zhǎng)度；Rm為測(cè)線兩端平均緯度處參考橢球面的平均曲率半徑；ym為測(cè)線在高斯平面上離中央子午線垂距的平均值。

3.3長(zhǎng)度綜合變形[10]

通過(guò)兩次變形后，地面真實(shí)長(zhǎng)度與高斯投影面的長(zhǎng)度之差為長(zhǎng)度綜合變形，按式(3)計(jì)算：

(3)

在實(shí)際計(jì)算中，為了方便起見(jiàn)，又不影響必要的精度，取Rm=RA=6 371 km，采用不同投影面的同一距離近似相等，即S≈D，可以寫成相對(duì)變形的形式為：

(4)

由公式(4)可知：采用國(guó)家統(tǒng)一坐標(biāo)系統(tǒng)所產(chǎn)生的長(zhǎng)度綜合變形與測(cè)區(qū)內(nèi)的平均高程和測(cè)區(qū)所處投影帶內(nèi)的位置有關(guān)，根據(jù)(4)式可以計(jì)算出已知測(cè)區(qū)內(nèi)的長(zhǎng)度相對(duì)變形的大小[11-13]。

可以把長(zhǎng)度綜合變形的容許值代入(4)式得：

H=0.783y2(10-4)±0.159

(5)

式中：以H為縱坐標(biāo)軸，y為橫坐標(biāo)軸，式(5)繪制成圖如圖1所示。

圖1　投影面高程與橫坐標(biāo)投影變形Figure 1　The height varying in the projection surface and 　the projection varying in abscissa

圖1中：適應(yīng)區(qū)是指A區(qū)、B區(qū)適合采用國(guó)家統(tǒng)一坐標(biāo)系，沒(méi)有必要建立獨(dú)立坐標(biāo)系，其他不適應(yīng)區(qū)域則需要建立獨(dú)立坐標(biāo)系才能滿足實(shí)際需要。其中C區(qū)域?yàn)楦吆０蔚貐^(qū)，長(zhǎng)度綜合變形中占主導(dǎo)作用的是高程歸算變形；在D區(qū)域?yàn)殡x中央子午線較遠(yuǎn)，則投影變形占主導(dǎo)作用的是高斯投影所產(chǎn)生的變形；在E區(qū)域內(nèi)則是前面兩者綜合影響長(zhǎng)度變形[14]。

4實(shí)例分析

本文以一范圍在東經(jīng)114°22′～ 114°41′，北緯 37°52′～38°11′的區(qū)域?yàn)槔M(jìn)行分析比較，該地區(qū)東西距114°子午線分別為yW= 32.2 km和yE= 60.0 km。兩邊的平均高程為HW=80 m和HE=61 m，該地區(qū)平均高程為Hu=70 m，高程異常為Hg=40 m。由公式(6)計(jì)算橫坐標(biāo)區(qū)間ym，ym=(0～58.6)km。

(6)

① 當(dāng)投影面高程H=0 m ，中央經(jīng)線L=114°時(shí)，規(guī)劃區(qū)的東西邊投影變形值分別為：

西邊的變形值為：

(7)

東邊的邊線值為：

(8)

從式(8)可知：東邊投影變形值已經(jīng)大于25 mm/km，所以需要選擇合適的方法來(lái)解決投影變形值過(guò)大的影響。

② 當(dāng)變更中央子午線L=114°30′，H=0 m時(shí)，選擇投影帶和投影面：

首先，計(jì)算市中心距統(tǒng)一3°帶中央子午線114°的橫坐標(biāo)近似值y0：

y0=N·l·cosB=1.857·l′·cosB=43.9 km

(9)

式中：平均緯度B=38°；市中心距中央子午線114°的經(jīng)差l′=30′；卯酉圈曲率半徑N=6 383 km。

由式(9)可知：該區(qū)域中心位于統(tǒng)一3°帶中央子午線(114°)以東 43.9 km，對(duì)新建的獨(dú)立坐標(biāo)系，則規(guī)劃區(qū)東邊位于中央子午線以東16.1 km(經(jīng)差11′，yE=16.1 km)，規(guī)劃區(qū)西邊位于中央子午線(114 °30′)以西11.7 km(經(jīng)差-8′，yW=11.7 km)，由此可計(jì)算出投影變形值分別為：

西邊變形值：

(10)

東邊變形值：

(11)

可見(jiàn)，中央子午線兩邊投影變形值都小于25 mm/km，選擇過(guò)該區(qū)域中心點(diǎn)的子午線作為中央子午線，然后投影到H=0 m的高程面上能基本解決投影變形過(guò)大的影響。

③ 當(dāng)改變投影面高程，投影到該區(qū)平均高程H=70 m的高程面上，仍選用114°為中央子午線建立獨(dú)立坐標(biāo)系，由式(8)可知：東邊的高斯投影變形值較大，高程歸化改正不能徹底的抵消其變形，假如中央子午線不改變，就只能人為的改變歸化高程來(lái)抵償部分高斯投影變形值，使其長(zhǎng)度變形仍然小于 1/40 000。

西邊的投影變形為：

(12)

東邊的投影變形為：

(13)

由式(13)可知：這種情況，東邊的投影變形值大于25 mm/km，不能滿足要求。

④ 選用114°30′為中央子午線，投影到市區(qū)平均高程H=70 m的高程面上建立標(biāo)系。

西邊的投影變形為：

(14)

東邊的投影變形為：

(15)

可見(jiàn)，東西兩邊投影變形值都小于25 mm/km，將城區(qū)投影到H=70 m高程面上，選擇過(guò)該區(qū)中心點(diǎn)的子午線作為中央子午線，也能夠解決投影變形值過(guò)大的影響。

⑤ 計(jì)算抵償面的高程，投影到抵償面高程上，中央經(jīng)線L=114°時(shí)。

(17)

假設(shè)該區(qū)中心點(diǎn)y0= 43.9 km處的投影變形的變形值為0，可求出高程修正值ΔH得：ΔH=41 m，可知抵償面的高程為：

Hd=Hu-(H+ΔH)=-81 m

(18)

這時(shí)規(guī)劃區(qū)西邊的歸化高度為：

(19)

東面的歸化高度為：

(20)

計(jì)算投影變形得：

西面的投影變形為：

(21)

東面的投影變形為：

(22)

計(jì)算抵償帶寬度，將式(19)、式(20)代入式(6)得：yW=4.59 kmyE=61.98 km

東西的橫坐標(biāo)跨度為：

Δy=yE-yW=57.39 km

(23)

根據(jù)式(18)、式(21)、式(22)，在不改變中央子午線(L=114°)的情況下，可以選擇Hd=-81 m的抵償高程面作為投影面，也能解決投影變形值過(guò)大的問(wèn)題，但變形值還是相對(duì)較大，由式(23)也可得知，東西跨度也比較小。

⑥ 在選擇Hd=-81 m低償高程面作為投影面建立獨(dú)立坐標(biāo)系，同時(shí)選擇過(guò)該區(qū)中心點(diǎn)y0= 43.9 km的子午線(經(jīng)度為114°30′)作為中央子午線，此時(shí)可以計(jì)算出投影變形值為：

西邊的投影變形為：

(24)

東邊的投影變形為：

(25)

計(jì)算抵償帶寬度：令式(24)、式(25)的右邊等于25 mm/km，計(jì)算東西邊距中央子午線的最大橫坐標(biāo)得：

yWmax=45.124 kmyEmax=45.126 km

東西的橫坐標(biāo)跨度為：

Δy=yEmax+yWmax=90.25 km

(26)

從表2中可以清晰的看到：第一種方法與第三種方法變形值都超過(guò)了25 mm/km的精度要求，第二種方法東邊、西邊變形值都比較大，第四種方法也不是最優(yōu)解，第五種改變抵償面高程方法在實(shí)際應(yīng)用中比較多，但所建立的坐標(biāo)系東西變形相對(duì)較大，橫坐標(biāo)跨度也很小，不利于以后的城區(qū)擴(kuò)建工作。而利用第六種方法，先分別計(jì)算東西歸化高度，然后得到抵償面高程，還采用最佳中央子午線建立坐標(biāo)系，不但能得到投影變形值最小，并且其橫坐標(biāo)跨度也能夠達(dá)到最大值，能充分滿足城市今后規(guī)劃布局的需要。

表2 投影變形值及橫坐標(biāo)跨度對(duì)比Table2 Theprojectiondistortionvalueconstrastthespanofabscissa投影面高程/m中央子午西邊投影變形/(mm·km-1)東邊投影變形/(mm·km-1)橫坐標(biāo)跨度/km0114° 6.128.5—0114°30'-17.13-12.65—70114° 11.245.76—70114°30'0.124.6—-81114° -12.52257.39-81114°30'-1.61.490.25

5結(jié)論

針對(duì)投影變形，通常兩種手段來(lái)實(shí)現(xiàn)投影帶和投影面的選擇，首先是在已知的數(shù)據(jù)基礎(chǔ)上，在excel表格里面進(jìn)行初步計(jì)算，看投影在統(tǒng)一坐標(biāo)系里面是否符合標(biāo)準(zhǔn)，如果不符合標(biāo)準(zhǔn)，就用兩種方法解決這個(gè)問(wèn)題，第一是改變投影面高程來(lái)控制變形，第二是改變投影帶的中央子午線來(lái)控制投影變形。

所有測(cè)量工作的第一要?jiǎng)?wù)是建立準(zhǔn)確的坐標(biāo)系，這項(xiàng)要?jiǎng)?wù)是以后測(cè)量工作順利進(jìn)行的關(guān)鍵。上文中的六種處理的情況有幾種能滿足要求，但建立科學(xué)合理的城市坐標(biāo)系還要依據(jù)現(xiàn)有的成果和城市發(fā)展規(guī)劃的要求，既要考慮當(dāng)前的又要著眼未來(lái)，所以為了變形最小，還有投影的有效東西跨度的最大化，在選用常用的抵償面的基礎(chǔ)上，采用第六種情況是最佳的選擇。所以，本文的重點(diǎn)是討論投影變形大小的投影帶和投影面的最佳選擇，是基于現(xiàn)實(shí)情況下的一種最佳考慮選擇。

[參考文獻(xiàn)]

[1]劉連旺，覃輝.6種曲線元上的道路中樁坐標(biāo)解算[J].中南公路工程，2006.31(4)：72-76.

[2]潘國(guó)強(qiáng).高速公路改擴(kuò)建工程中路基加寬方式及其特點(diǎn)綜述[J].中南公路工程，2007，32(5)：34-38.

[3]葛婷，王曉飛，符鋅砂.改擴(kuò)建高速公路工程交通組織方法探討[J].中南公路工程，2010，35(6)：118-123.

[4]楊國(guó)清，王軍德.假3度帶坐標(biāo)系及關(guān)于大城市坐標(biāo)系的建議[J].測(cè)繪科學(xué)，2011，36(5)：22-24.

[5]吳云孫.淺談廣州市平面控制測(cè)量坐標(biāo)系統(tǒng)的選取[J].測(cè)繪通報(bào)，2005(8)：42-43.

[6]覃輝.測(cè)量程序與新型全站儀的應(yīng)用[M].北京：機(jī)械工業(yè)出版社，2006.

[7]范一中，王繼剛，趙麗華.抵償投影面的最佳選取問(wèn)題[J].測(cè)繪通報(bào)，2000(2)：20-21.

[8]牛卓立，賴昌意.建立區(qū)域坐標(biāo)系問(wèn)題的討論[J].測(cè)繪工程，1998，7(4)：13-16.

[9]許婭婭，張碧琴.公路施工測(cè)量百問(wèn)[M].北京：人民交通出版社，2006.

[10]田青文.測(cè)量學(xué)(下)[M].西安：西安地質(zhì)學(xué)院出版社，1994，90-115.

[11]武瑞宏.工程獨(dú)立坐標(biāo)系的建立方法研究[J].工程勘察，2009(3)：68-71.

[12]王繼剛，王堅(jiān).具有抵償面得任意待高斯投影直角坐標(biāo)系的選取方法[J].測(cè)繪通報(bào)，2002(11)：31-32.

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

[14]張華海，王寶山，趙長(zhǎng)勝，等.應(yīng)用大地測(cè)量學(xué)[M].徐州：中國(guó)礦業(yè)大學(xué)出版社，2007.

The Selection of Projection Zone and Plane is Based on the Size of Projection Distortion

YANG Zhi1，2，3， QIN Hui1，2，3， WEN Hongyan1，2，3， YUAN Mingyue1，2，3

(1.Guangxi Key Laboratory of Spatial Information and Geomatics， Guilin， Guangxi， 541004， China；2.College of Geomatics and Geoinformation， Guilin University of Technology， Guilin， Guangxi， 541004， China；3.Guangxi Scientific Experiment Center of Mining， Metallurgy and Environment， Guilin University of Technology， Guilin， Guangxi 541006， China)

[Abstract]China’s national coordinate system commonly used Gauss-Kruger orthographic projections( abbreviation “Gauss Projections” ).That is to say the projections go first to the ground，and then to the reference elliTPoid.Next，the projections will finally projected onto the Gaussian flat.After the transform of above projections，it will cause some sort of deforms.The scale of the deforms will be varied with the different distance between the projection point and the geoid.Usually.Aimed at the problem，The approaches be used in engineering survey are：a.Use the form of zonation randomly to put the central meridian to the range that the projection distortion can be able to control；b.According to the actual needs，reset the elevation of projection plane；c.Select the compensation plane.In this paper，we should find the best central meridian at first，and then calculate the elevation of compensating surface，we will get the optimal projection zone and projection plane.By analysis of calculated result obtained the optimal projection zone and projection plane and the span of abscissa is the peak of maximum within a certain range，it also be able to accurately determine the coordinate system of the plane geodetic control network at the same time.

[Key words]gauss-kruger orthographic projections； projection distortion； the central meridian； the elevation of projection plane； compensation plane

[中圖分類號(hào)]U 412.24+1

[文獻(xiàn)標(biāo)識(shí)碼]A

[文章編號(hào)]1674-0610(2016)01-0059-05

[作者簡(jiǎn)介]楊志(1989-)，男，湖南桃江人，研究生，變形監(jiān)測(cè)與數(shù)據(jù)處理。[通訊作者] 文鴻雁(1963-)，男，湖南祁陽(yáng)人，教授，博士，研究方向?yàn)榫芄こ虦y(cè)量與專題信息系統(tǒng)。

[基金項(xiàng)目]國(guó)家自然科學(xué)基金項(xiàng)目(41461089)；廣西“八桂學(xué)者”崗位專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室資助課題(桂科能130511402，1207115-06)；廣西自然科學(xué)基金基金項(xiàng)目(2014GXNSFAA118288)；廣西礦冶與環(huán)境科學(xué)實(shí)驗(yàn)中心資助課題(KH2012ZD004)；廣西研究生教育創(chuàng)新計(jì)劃項(xiàng)目(YCSZ2014151，YCSZ2012083)

[收稿日期]2014-10-23