網(wǎng)絡(luò)同步RTK定位原理與精度分析

賀俊偉 李春華 劉 俊 鄧 芳

(1.成都市勘察測繪研究院，四川成都 610081;2.四川省地質(zhì)測繪院，四川成都 610017)

圖1 汶川地震災(zāi)區(qū)CORS基站分布

圖2 天津市CORS系統(tǒng)分布示意

CORS NRTK(Continuously Operating Reference Stations Network Real－Time Kinematic positioning)是集衛(wèi)星大地測量技術(shù)、通信技術(shù)和計算機(jī)網(wǎng)絡(luò)技術(shù)于一體的空間信息實(shí)時服務(wù)技術(shù)，整個系統(tǒng)包括連續(xù)運(yùn)行參考站、數(shù)據(jù)處理與控制中心以及流動用戶部分［1］。

目前，我國在多個地區(qū)和大中城市已經(jīng)建成400個左右的CORS(Continuous Operation Reference Stations縮寫為 CORS)基站［2］，如5.12汶川地震災(zāi)區(qū)(如圖1)、廣東省、地震局和氣象局的專用CORS系統(tǒng)，以及深圳、上海、北京、天津(如圖2)、蘇州、成都等城市GPS綜合服務(wù)系統(tǒng)。

CORSNRTK測量方法(即網(wǎng)絡(luò)RTK)極大地克服了傳統(tǒng)RTK技術(shù)作業(yè)距離方面的限制，其主要優(yōu)勢表現(xiàn)為:(1)可以實(shí)時獲取測站精確的WGS84坐標(biāo)和用戶坐標(biāo);(2)單臺GPS接收機(jī)作業(yè)模式，減少作業(yè)成本，提高工作效率;(3)縮短觀測時間(幾秒)，提高了RTK測量的可靠性和定位精度;(4)其作業(yè)范圍不僅適合于基準(zhǔn)站網(wǎng)所包圍的范圍，而且可以擴(kuò)大到基準(zhǔn)站網(wǎng)外的一定距離，(5)基準(zhǔn)網(wǎng)內(nèi)和網(wǎng)外一定范圍內(nèi)定位精度均勻［3］。

CORSNRTK測試表明:兩次獨(dú)立初始化的坐標(biāo)較差、CORS網(wǎng)內(nèi)與CORS網(wǎng)外(約60 km范圍內(nèi))網(wǎng)絡(luò)RTK的絕對定位精度基本一致，平面定位精度(中誤差)約 ±3.0 cm，大地高定位精度(中誤差)約 ±5 cm［3，4］。

由于每次CORS NRTK測量均為獨(dú)立觀測，每點(diǎn)的精度基本一致，故相鄰點(diǎn)之間的相對定位精度不高，難以滿足城市日常測量(300 m以內(nèi)的邊長)工作的需要。

在CORS NRTK定位精度的基礎(chǔ)上，提出了改善相鄰點(diǎn)之間相對定位精度的作業(yè)方法和數(shù)據(jù)處理模型:即同步RTK測量和測邊網(wǎng)平差或三維坐標(biāo)差平差。

1 RTK測量邊長相對精度分析

下面就CORSNRTK測量的二維平面精度進(jìn)行詳細(xì)的分析。

兩點(diǎn)的邊長計算公式為

其中，Sij相鄰點(diǎn)之間的平面距離;

(Xi，Yi)和(Xj，Yj):RTK 測量并轉(zhuǎn)換后的平面坐標(biāo)。

線形化后

假設(shè) Xi，Xj，Yi，Yj的中誤差相等，均為 ± m，按照誤差傳播定律，則邊長中誤差為

取mX=mY=±2 cm，則

如需保證邊長相對誤差不大于1/1萬，則相鄰點(diǎn)的距離應(yīng)不小于280 m，而實(shí)際作業(yè)中，邊長長度(尤其是在城市建成區(qū)域內(nèi))難以完全滿足此項(xiàng)要求。

2 同步RTK測量原理與實(shí)驗(yàn)方法

在CORS基站網(wǎng)內(nèi)，選取八個GPS點(diǎn)，分別采用靜態(tài)GPS測量、同步RTK測量、非同步RTK測量獲取各點(diǎn)的觀測信息(三維坐標(biāo) B、L、H 或 X、Y、Z)。

觀測方法1:靜態(tài)GPS測量，采用Trimble 5700雙頻GPS接收機(jī)，同步觀測45～60 min，同步環(huán)之間采用邊連接模式，靜態(tài)GPS觀測網(wǎng)見圖3。該GPS網(wǎng)中，最長邊約為746 m，最短邊約為140 m，小于280 m的邊長有5條。進(jìn)行三維無約束平差和二維約束平差計算，獲取各點(diǎn)的靜態(tài)GPS點(diǎn)二維坐標(biāo)(x，y)，并將此坐標(biāo)作為真值。

圖3 靜態(tài)GPS觀測網(wǎng)

觀測方法2:同步RTK測量，每相鄰點(diǎn)之間進(jìn)行同步RTK觀測，即RTK觀測時間相同，每個同步RTK觀測歷元約為60個，共獲得15條同步觀測結(jié)果(X、Y、Z)，同步RTK觀測示意見圖4。根據(jù)兩點(diǎn)之間的同步觀測值計算坐標(biāo)差和邊長(改化為平距)，并分別按照測邊網(wǎng)(觀測值為Si)和三維坐標(biāo)差(觀測值為ΔXij，ΔYij，ΔZij)進(jìn)行平差計算(二維約束條件與靜態(tài)GPS網(wǎng)一致)，獲取各點(diǎn)的平面坐標(biāo)，與靜態(tài)測量真值坐標(biāo)進(jìn)行比較，并統(tǒng)計其精度指標(biāo)。

圖4 同步RTK觀測示意

方法3:采用非同步RTK測量方法直接獲取每點(diǎn)的坐標(biāo)，每點(diǎn)觀測歷元數(shù)為60個。

3 試驗(yàn)結(jié)果

以靜態(tài)GPS數(shù)據(jù)為基準(zhǔn)，按照式(5)分別計算非同步RTK、同步RTK、測邊網(wǎng)平差、三維(二維)坐標(biāo)差平差的相對誤差。

共統(tǒng)計10條GPS邊長，各邊長度和非同步RTK測量各邊長相對誤差見表1。

從表1可以發(fā)現(xiàn)，最弱邊相對誤差為92.2×10－6，最強(qiáng)邊相對誤差為5.5 ×10－6，平均相對誤差為34.3 ×10－6。

非同步RTK測量相對誤差分布見圖5。

采用不同數(shù)據(jù)處理方法計算邊長(共10條邊)相對誤差，各種數(shù)據(jù)處理方法得到的各條邊相對中誤差分布見圖5(橫軸表示邊長，縱軸表示相對誤差)。統(tǒng)計結(jié)果見表2，統(tǒng)計圖見圖6。

表1 非同步RTK邊長相對誤差

表2 RTK測量相對誤差統(tǒng)計 ×10－6

圖5 邊長相對中誤差分布

從圖5、圖6和表2不難發(fā)現(xiàn)，相對誤差平均值和最大值大小順序依次為:非同步RTK→同步RTK(不平差)→測邊網(wǎng)平差→坐標(biāo)差平差。其中，三維坐標(biāo)差數(shù)據(jù)處理模式的邊長相對精度約為非同步觀測的1/3，相鄰點(diǎn)之間(短邊，邊長范圍在100～1 000 m之間)相對精度可以達(dá)到50×10－6(即《衛(wèi)星定位城市測量技術(shù)規(guī)范》中二級GNSSRTK)的技術(shù)要求。

在實(shí)際控制網(wǎng)加密作業(yè)過程中，可采用同步RTK測量模式，或同步RTK測邊網(wǎng)平差和三維坐標(biāo)差平差，以提高控制點(diǎn)間相鄰邊長的相對精度。

圖6 不同數(shù)據(jù)處理方法的相對誤差統(tǒng)計

4 結(jié)論及建議

同步RTK測量、或同步RTK測量的間接觀測值(邊長或坐標(biāo)差)進(jìn)行數(shù)據(jù)處理，可以極大地提高相鄰點(diǎn)間的相對精度，可以將《衛(wèi)星定位城市測量技術(shù)規(guī)范》(CJJ/T73—2010)中二級GNSS RTK平面測量中相鄰點(diǎn)之間的最短距離從300 m縮短為100 m，并可完全取代城市一、二、三級導(dǎo)線測量，極大地拓寬了CORSNRTK的適用范圍。

［1］周樂韜．連續(xù)運(yùn)行參考站網(wǎng)絡(luò)實(shí)時動態(tài)定位理論、算法和系統(tǒng)實(shí)現(xiàn)［D］．成都:西南交通大學(xué)，2007

［2］中國科學(xué)技術(shù)協(xié)會，中國測繪學(xué)會．測繪科學(xué)與技術(shù)學(xué)科發(fā)展報告［M］．北京:中國科學(xué)技術(shù)出版社，2010

［3］李春華．基于網(wǎng)絡(luò)GPS和精化大地水準(zhǔn)面的區(qū)域?qū)崟r三維定位理論與應(yīng)用［D］．成都:西南交通大學(xué)，2010

［4］韓國超，劉莉萍，李春華，等．NRTK定位原理與性能測試［J］．測繪通報，2010(6):30－33．