黃偉 馮濤



摘 要:垂線偏差在大地測量中是一項(xiàng)重要的數(shù)據(jù),傳統(tǒng)的測量方法受天氣狀況影響較大。該文介紹了一種利用GPS和精密水準(zhǔn)儀來測定垂線偏差的計(jì)算方法,實(shí)驗(yàn)結(jié)果表明,在一定條件下,該方法的精度符合要求。
關(guān)鍵詞:垂線偏差;GPS;嵌入式
在大地測量工程中,為了便于計(jì)算,人們將地球假想成一個(gè)理想的橢球體,建立地心直角坐標(biāo)系,從而把測量觀測數(shù)據(jù)投射到橢球體上進(jìn)行計(jì)算,著名的WGS-84就是國際普遍采用的以GPS為基礎(chǔ)建立的坐標(biāo)系。實(shí)際上,地球并不是一個(gè)質(zhì)量均勻分布的橢球體,等重力勢(shì)面的大地水準(zhǔn)面和橢球面就有了差別,因此,大地水準(zhǔn)面的垂線和橢球面的法線就有了偏差,我們稱之為垂線偏差。
垂線偏差的測量方法有很多,傳統(tǒng)的有天文經(jīng)緯度測量法,重力測量法等,而GPS水準(zhǔn)測量法受天氣影響較小,比較有效的改善了測量的效率和精度。
1 GPS測定垂線偏差算法
1.1 垂線偏差定義
垂線偏差是指大地水準(zhǔn)面的鉛垂線和橢球面上法線方向的夾角,在大地測量中是一項(xiàng)重要的工作。測繪工程中,觀測數(shù)據(jù)是以橢球面法線為基準(zhǔn)進(jìn)行計(jì)算。橢球面與大地水準(zhǔn)面的關(guān)系如圖1所示(這里所指的橢球面是以WGS-84定義的基準(zhǔn)橢球體)。
1.2 垂線偏差測量原理
橢球面和大地水準(zhǔn)面存在著高度差, 我們假設(shè):H為正常高,即海拔高,h為大地高,即離橢球高,N為高程異常(橢球面和大地水準(zhǔn)面的距離)。若在實(shí)際地面做一條基線AB,由此可以推斷出:
做兩式差,
其中△hAB可以通過GPS數(shù)據(jù)得出,△HAB由精密水準(zhǔn)儀測量得出。
在橢球面上建立以A點(diǎn)為原點(diǎn),法線方向?yàn)閆軸,經(jīng)度方向?yàn)閅軸,維度方向?yàn)閆軸的站心直角坐標(biāo)系,如圖2所示,u即是所求的垂線偏差角。假設(shè)u在經(jīng)線方向和為向方向的分量分別為ε和η,而在基線方向的分量角為θ,可以得出:
其中φ為垂線偏差基于經(jīng)線方向角,ψ為基線基于經(jīng)線方向角。而對(duì)于基線方向來說,如圖3所示。
在θ角度較小的時(shí)候可以推出近似公式
其中△NAB可以通過公式解出,而AB可以通過GPS數(shù)據(jù)計(jì)算得出。所以當(dāng)測量兩條基線,或者做多條基線采用最小二乘法便可以通過公式5聯(lián)立方程組解出垂線偏差角。
2 GPS數(shù)據(jù)采集
在linux系統(tǒng)下,通過串口對(duì)GPS接收機(jī)模塊進(jìn)行讀寫,解析GPS發(fā)送的數(shù)據(jù)幀,分離GPS數(shù)據(jù)電文,獲取經(jīng)緯度,以及高差異常數(shù)據(jù),再根據(jù)水準(zhǔn)儀得到的海拔高差信息解出垂線偏差。
2.1 串口通信
GPS接收機(jī)工作后,會(huì)把信息通過串口源源不斷的傳給主機(jī),串口只要負(fù)責(zé)將傳送來的數(shù)據(jù)進(jìn)行讀操作并緩存。在linux系統(tǒng)下具體實(shí)現(xiàn)的步驟是:首先對(duì)串口進(jìn)行配置,打開串口,讀串口數(shù)據(jù),將接收到的GPS串口數(shù)據(jù)緩存,分離GPS數(shù)據(jù)幀,關(guān)閉串口。
2.2 NMEA0183電文分離
NMEA-0183協(xié)議是GPS接收機(jī)應(yīng)當(dāng)遵守的標(biāo)準(zhǔn)協(xié)議,采用ASCII碼,其串行通信默認(rèn)參數(shù)為:波特率為4800bps,數(shù)據(jù)位為8bit,開始位為1bit,停止位為1bit,無奇偶校驗(yàn)。NMEA-0183協(xié)議定義的語句非常多,常用的有$GPGGA、$GPGSA、$GPGSV、$GPRMC、$GPVTG、$GPGLL等。而對(duì)于我們來說,從上一節(jié)來看,需要得到的信息為測量點(diǎn)的經(jīng)緯度信息,高度信息,而$GPGGA語句正好包含我們所需要的信息,所以我們只需要將其提取,并進(jìn)行分割即可。
其標(biāo)準(zhǔn)格式如下:
$GPGGA,<1>,<2>,<3>,<4>,<5>,<6>,<7>,<8>,<9>,M,<10>,M,<11>,<12>*hh
這里我們可以通過提取經(jīng)緯度以及高度信息,并通過坐標(biāo)系變換解出之前提出的空間距離AB:設(shè)U點(diǎn)的大地坐標(biāo)系為U(λ,ψ,h),則通過以下公式變換可以解出其在地心直角坐標(biāo)系的坐標(biāo)U(X,Y,Z),進(jìn)而解出基線空間距離。
至此,所需要的測量數(shù)據(jù)都直接或間接計(jì)算得到,只要根據(jù)第一章節(jié)所討論的方法進(jìn)行解算,就可以解出所要求的垂線偏差。
3 誤差分析
在垂線偏差角分量角公式中: ,兩邊求微分得到: 。因?yàn)閐l對(duì)垂線偏差的影響很小所以省略后一項(xiàng),所以其誤差公式:
可以看出當(dāng)m△H和m△h一定時(shí),mθ和l有關(guān)。
Matlab仿真圖:
如圖4、圖5所示,一定范圍內(nèi),隨著基線長度的增加,在相同偏差分量角誤差的要求下高差誤差變寬,相同高差誤差下,基線長度越長偏差分量角誤差越小(在考慮基線在水平面的情況下)。
如圖6所示,可以看出一定范圍內(nèi)隨著基線長度的增長,偏差分量角的誤差也逐漸變小。
4 結(jié)語
利用GPS測量垂線偏差是一種時(shí)效性比較高的測量技術(shù),地形在一定范圍內(nèi)隨線性變化的場合下,GPS測量垂線偏差是一種精度比較高的理想途徑。
[參考文獻(xiàn)]
[1]肖榮健,鄒強(qiáng).垂線偏差的確定方法[J].大眾科技,2009(8):58-59.
[2]馬利敏.GPS技術(shù)在地質(zhì)測量中的應(yīng)用[J].科學(xué)與財(cái)富,2012(8):212.
[3]王愛生.利用GPS和水準(zhǔn)測量解算垂線偏差[J].測繪通報(bào),2002(2):17-20.
[4]吳紅蘭,黃圣國.嵌入式GPS數(shù)據(jù)采集與分析系統(tǒng)設(shè)計(jì)[J].實(shí)驗(yàn)室研究與探索,2011(1):54-57.
[5]熊波如,孫金坤,李培君,等.高層建筑GPS精密變形監(jiān)測設(shè)計(jì)[J].世界家苑,2011(9):155.