1 引言
在水利工程測量中,多數(shù)情況下工程所處位置地形復雜,交通不便,通視條件較差,采用以經(jīng)緯儀、全站儀測量為代表的常規(guī)測量常常效率低下。隨著GPS-RTK測量系統(tǒng)的使用,由于它具有觀測速度快,定位精度高,經(jīng)濟效益高等特點,現(xiàn)在我院多數(shù)水利工程測量都是采用RTK測量技術(shù)來完成。對于GPS-RTK系統(tǒng)來說,由于它采用的是WGS-84固心坐標系,而在實際工程應(yīng)用中,由于顧及長度變形、高程異常等影響而采用獨立坐標系,這就需要將RTK測量采集的數(shù)據(jù)在兩坐標系中進行轉(zhuǎn)換。
2 國家坐標系及獨立坐標系的建立
2.1 國家坐標系的建立
在我國,由于歷史原因先后采用不同的參考橢球體和大地起算數(shù)據(jù)而形成多個國家坐標系,主要國家坐標系有1954北京坐標系、1980西安坐標系、2000國家坐標系和WGS-84坐標系。前兩個是參心坐標系,后兩個是固心坐標系。由于他們采用不同的橢球體參數(shù),所以地面上同一個點在不同的坐標系中有不同的坐標值。
國家坐標系的主要作用是在全國建立一個統(tǒng)一的平面和高程基準,為發(fā)展國民經(jīng)濟、空間技術(shù)及國防建設(shè)提供技術(shù)支撐,也為防災(zāi)、減災(zāi)、環(huán)境監(jiān)測及當代地球科學研究提供基礎(chǔ)資料。
2.2 獨立坐標系的建立
在工程應(yīng)用中,由于起算數(shù)據(jù)收集困難、測區(qū)遠離中央子午線及滿足特殊要求等諸多原因,如在水利工程測量中,常要測定或放樣水工建筑物的精確位置,要計算料場的土石方貯量和水庫的庫容。規(guī)范要求投影長度變形不大于一定的值(如《工程測量規(guī)范》為2.5cm/km,《水利水電工程測量規(guī)范(規(guī)范設(shè)計階段)》為5.0cm/km)。如果采用國家坐標系統(tǒng)在許多情況下(如高海拔地區(qū)、離中央子午線較遠地方等)不能滿足這一要求,這就要求建立地方獨立坐標系。
在常規(guī)測量中,這種獨立坐標系只是一種高斯平面直角坐標系,而在采用GPS-RTK采集數(shù)據(jù)時,獨立坐標系就是一種不同于國家坐標系的參心坐標系。
跟國家坐標系一樣,建立獨立坐標要確定的主要元素有:坐標系的起算數(shù)據(jù)、中央子午線、參考橢球體參數(shù)及投影面高程等。對于起算數(shù)據(jù),可以采用國家坐標系的坐標和方位角或任意假設(shè)坐標和方位角。在RTK測量中,我們常采用基線的某一端點的單點定位解作為起點,然后以另一點定向,用測距儀測出基線邊長,經(jīng)改正后算出基線端點的坐標;中央子午線常采用測區(qū)中央的子午線;投影面常采用測區(qū)的平均高程面。參考橢球體一般是基于原來的參考橢球體做某種改動,使改變后的參考橢球面與投影面擬合最好,投影變形可以減到最小,也便于與國家坐標系統(tǒng)進行換算。
3 坐標系的轉(zhuǎn)換
GPS-RTK接收機采集的坐標數(shù)據(jù)是基于WGS-84橢球下的大地坐標,而我們經(jīng)常使用的獨立坐標系是基于某種局部橢球體下的平面直角坐標,這兩種坐標是不同坐標基準下的兩種表現(xiàn)形式。利用WGS-84下的大地坐標來推求獨立坐標系中的平面直角坐標,必然要求得兩坐標系之間轉(zhuǎn)換參數(shù)。求取轉(zhuǎn)換參數(shù)的基本思路是利用兩坐標系中必要個數(shù)的公共點,根據(jù)相應(yīng)的橢球參數(shù)及中央子午線采用最小二乘法嚴密平差解算轉(zhuǎn)換參數(shù),具體操作是由轉(zhuǎn)換模型把不同坐標基準下的坐標轉(zhuǎn)換為同基準下的不同坐標形式,再進行同基準下不同坐標形式的轉(zhuǎn)換,從而得到所要的獨立坐標系中的平面直角坐標。轉(zhuǎn)換的難點是WGS-84橢球與獨立坐標系局部橢球的變換。
3.1 常用的坐標轉(zhuǎn)換方法
以一個點為基準點,進行橢球變換的方法稱為單點模式。這個基準點可以是實際觀測的現(xiàn)實存在的點,也可以是使用多個點歸算得到的一個等效的虛擬點。主要方法有橢球膨脹法、橢球平移法和橢球變形法。
以多個點作為橢球變換的參考點稱為多點模式。如武測模型法求定轉(zhuǎn)換參數(shù)。
在各種變換方法中,常采用橢球膨脹法來轉(zhuǎn)換不同基準下的大地坐標,它是橢球變換的一種特殊情況,設(shè)某點在WGS-84中的大地坐標為(),膨脹前后橢球中心、方向保持不變,橢球扁率保持不變,橢球長半軸變化 =。其它各參數(shù)的變化為, ,,,,,,,即平移,旋轉(zhuǎn)及比例變化均為0。代入橢球變換模型的廣義大地微分公式得各變化值:
(1)
將(1)式代入(2)式中即得到轉(zhuǎn)換后局部橢球的大地坐標()為: (2)
具體計算方法見文獻[4]。接著將局部橢球的大地坐標()按選定的中央子午線經(jīng)過高斯正算模型計算可得到局部橢球下的高斯平面直角坐標(),然后結(jié)合兩已知點的坐標進行四參數(shù)(兩個平移參數(shù),一個旋轉(zhuǎn)參數(shù)和一個比例參數(shù))的求取,上面的計算可以通過南方測繪工程之星或其它商用軟件來完成。將求得的四參數(shù)輸入RTK采集器中就可進行數(shù)據(jù)采集了。
3.2 南方測繪工程之星(Engineering Star2.0)的坐標轉(zhuǎn)換方法
ES2.0(Engineering Star2.0)是南方測繪公司開發(fā)的基于WinCE操作系統(tǒng)下一款針對GPS-RTK測量系統(tǒng)的圖形化野外測量數(shù)據(jù)處理和管理軟件。由于它操作簡便、高效、功能強大且符合國內(nèi)測量行業(yè)的野外生產(chǎn)習慣,因此,在測量和相關(guān)領(lǐng)域應(yīng)用廣泛。
使用ES2.0進行RTK作業(yè)時,我們首先要進行控制點轉(zhuǎn)換平面坐標四參數(shù)和高程擬合參數(shù)的確定,在數(shù)據(jù)采集過程中測量點的坐標轉(zhuǎn)換流程如下:
(1)把WGS84大地坐標當作北京54大地坐標直接投影在高斯平面上,得到公共點在高斯平面上的北京54平面直角坐標,得到的成果是近似值,具有一定的系統(tǒng)誤差。
(2)通過四參數(shù)的轉(zhuǎn)換得到獨立坐標系的平面坐標,這里使用的是我們上面介紹的由WGS84坐標和獨立坐標系中坐標求取的四參數(shù)。
(3)通過測區(qū)內(nèi)的高程點擬合參數(shù)得到獨立坐標系的高程。
(4)內(nèi)業(yè)數(shù)據(jù)處理,用PSION手簿與電腦連接后,采集數(shù)據(jù)可以直接轉(zhuǎn)換為南方CASS繪圖軟件DAT文件數(shù)據(jù)格式的獨立坐標系平面直角坐標,也可下載RTK文件數(shù)據(jù)格式的WGS-84大地坐標予以保留。在這兩個文件格式中,我們可以直接查看、修改、刪除不符合采集精度要求的點的坐標、高程。
從它的轉(zhuǎn)換過程可以看出,在使用ES2.0進行RTK作業(yè)采集測量點時,它的坐標轉(zhuǎn)換方式與傳統(tǒng)方式有所不同,這也是ES2.0的獨到之處,它提高了坐標數(shù)據(jù)轉(zhuǎn)換的速度且最后的結(jié)果是一樣的。
3.3 ES2.0在工程中的應(yīng)用
在建水縣洗馬塘水庫測量中,根據(jù)收集的資料現(xiàn)場沒有高等級平面控制點,只有一個黃海高程點。在采用GPS-RTK測量時我們按上述方法建立獨立坐標系,基本思路是在獨立坐標系中選擇該水庫平均高程面1970m為橢球投影面,中央子午線為測區(qū)中央子午線102°50′,四參數(shù)轉(zhuǎn)換的兩公共點采用控制點“BZH”和“K1”點(其中“BZH”為基準點,“K1”為定向點),各控制點采用水準測量聯(lián)測,以便進行高程擬合。具體在RTK采集器中的操作如下:第一步,建立工程名稱,接著在橢球設(shè)置中橢球系名稱項選擇自定義橢球 “User Define”,橢球系長軸值為,橢球系扁率為 =298.257223563;第二步,在橢球參數(shù)設(shè)置中確定投影方式為“高斯投影”,輸入中央子午線為102.50,投影高為1970m;第三步,在四參數(shù)設(shè)置中選擇“啟用四參數(shù)”;第四步,四參數(shù)的求取,兩已知控制點坐標為下表1所示,通過計算后得到的轉(zhuǎn)換參數(shù)為:平移參數(shù) ,,旋轉(zhuǎn)參數(shù),比例參數(shù)。第五步,新建一個文件就可以進行目標點測量了。
表1 公共控制點的坐標表
在測量過程中選取了幾個點同時用RTK測量和全站儀測量,其點位較差如下表2。說明RTK測量技術(shù)能滿足四等以下控制測量、工程測量及地形測量工作。
表2 RTK測量與全站儀測量點位較差表
4 結(jié)語
在采用GPS-RTK技術(shù)測量時,建立地方獨立坐標系或連測國家坐標系確定坐標轉(zhuǎn)換參數(shù)是一項工程一開始就要做的事前。在這里只討論了常用的平面坐標四參數(shù)的求取和南方測繪工程之星中的坐標轉(zhuǎn)換方法,還有GPS-RTK測量中高程擬合方法,點校正的三參數(shù)法和坐標轉(zhuǎn)換的七參數(shù)法沒有討論,其基本原理是一樣的。清楚地認識這些原理,對于我們在分析、處理RTK數(shù)據(jù)過程中剔除粗差、減小誤差及提高效率方面都很有幫助。
參考文獻:
[1] GB 50026-93.工程測量規(guī)范[S].
[2] SL 197-97.水利水電工程測量規(guī)范[S].
[3] 劉大杰,施一民,過靜珺.全球定位系統(tǒng)(GPS)的原理與數(shù)據(jù)處理[M].上海:同濟大學出版社,1999.
[4] 孔祥元,郭際明,劉宗泉.大地測量學基礎(chǔ)[M].武漢:武漢大學出版社,2005.
[5] 董鈞祥,楊德宏.測量坐標轉(zhuǎn)換及其應(yīng)用[J].昆明理工大學學報(理工版),2006,31(3),1-4.
[6] 李世安,劉經(jīng)南,應(yīng)用GPS建立區(qū)域獨立坐標系中橢球變換的研究[J].武漢大學學報(信息科學版),2005,30(10),1-4.