TBC和LGO在GNSS數(shù)據(jù)處理中的對比和應用研究

胡玉祥，王智，張洪德，孟慶年

(青島市勘察測繪研究院，山東 青島 266032)

1 引 言

在地鐵建設(shè)中，建立高精度的GPS首級控制網(wǎng)必不可少，GPS數(shù)據(jù)觀測往往混合使用不同類型接收機[1]。一般來說，不同GPS接收機廠商有自己專有格式的數(shù)據(jù)存儲形式，需要相應配套的數(shù)據(jù)處理軟件進行處理、轉(zhuǎn)換[2，4，5]。RINEX格式是一種與接收機類型無關(guān)的數(shù)據(jù)交換格式，該格式采用文本文件形式存儲數(shù)據(jù)，方便直觀，這就為不用接收機數(shù)據(jù)混合處理提供了方便。

LGO和TBC作為兩款流行的GPS數(shù)據(jù)處理隨機軟件，由于解算模型和生產(chǎn)工藝面的差異，兩者解算精度有較大區(qū)別[7，9]。考慮到TBC軟件的優(yōu)勢，本文將青島地鐵徠卡接收機RINEX格式數(shù)據(jù)用TBC基線解算后平差，并與LGO基線解算后平差結(jié)果進行對比，結(jié)果表明，TBC基線解算結(jié)果優(yōu)于LGO基線解算結(jié)果。

2 GPS靜態(tài)數(shù)據(jù)處理

從算法角度講，GPS靜態(tài)數(shù)據(jù)處理流程分為數(shù)據(jù)傳輸、格式轉(zhuǎn)換、基線解算和網(wǎng)平差四個階段(如圖1所示)。數(shù)據(jù)傳輸主要利用配套的隨機軟件，以下主要敘述格式轉(zhuǎn)換、基線解算和網(wǎng)平差的相關(guān)知識。

2.1 格式轉(zhuǎn)化

TEQC是由UNAVCO Faclity研制的為地學研究GPS監(jiān)測站數(shù)據(jù)管理服務(wù)的公開免費軟件，主要有格式轉(zhuǎn)換(Translate)、數(shù)據(jù)編輯(Edit)、質(zhì)量檢查(QualityCheck)單點定位(Coordinate)四方面的功能[3]。

為了保證軟件的兼容性，將接收機數(shù)據(jù)轉(zhuǎn)換成標準的RINEX格式，把TRIMBLE接收機所接收到的數(shù)

圖1 GPS靜態(tài)數(shù)據(jù)處理流程

據(jù)文件轉(zhuǎn)換成RINEX觀測數(shù)據(jù)文件、導航數(shù)據(jù)文件和氣象文件，命令格式如下：

teqc -tr do -week 1323 +nav trimbel.05n trimbel.dat > trimbel. 05o

式中，-tr 指明接收機的類型為trimble；do 指明輸入文件為Dat 文件、輸出為RINEX 觀測數(shù)據(jù)文件(o文件)；-week 1323 (可選) 指明觀測日期對應的GPS 周數(shù)，或以年/月/日方式表示(即上述命令也可以表示成teqc -tr do -week 2015/05/18 +nav trimbel.15n trimbel.dat > trimbel. 15o)；+nav 指明同時輸出RINEX 導航數(shù)據(jù)文件；轉(zhuǎn)換結(jié)果文件為觀測數(shù)據(jù)文件trimbel.15o 和導航文件trimbel.15n。

同理，將徠卡接收機數(shù)據(jù)文件轉(zhuǎn)換成RINEX格式，命令如下：

teqc -lei mdb +nav 0319_2508.10n 0319_2508.m00 > 0319_2508.10o

2.2 基線解算

基線解算過程就是由多臺GPS接收機在野外通過同步觀測采集到的觀測數(shù)據(jù)，確定接收機間基線向量及其方差-協(xié)方差陣的過程。

同一時段觀測值的數(shù)據(jù)剔除率宜小于10%。GPS控制測量外業(yè)觀測的全部數(shù)據(jù)應經(jīng)同步環(huán)、異步環(huán)及復測基線檢核，并應滿足下列要求：

(1)同步環(huán)各坐標分量閉合差及環(huán)線全長閉合差，應滿足下列各式要求：

式中：n——同步環(huán)中基線邊的個數(shù)

W——同步環(huán)環(huán)線全長閉合差(mm)

σ——基線測量中誤差，單位為毫米(mm)

(2)異步環(huán)閉合符合下式的規(guī)定：

式中：n——異步環(huán)中基線邊的個數(shù)

σ——基線測量中誤差，單位為毫米(mm)

(3)復測基線較差不超過下式規(guī)定

式中：σ——基線測量中誤差，單位為毫米(mm)

數(shù)據(jù)檢驗中，當重復基線、同步環(huán)、異步環(huán)或附合路線中的基線超限時，應舍棄基線后重新構(gòu)成異步環(huán)，所含的異步基線數(shù)和閉合差應符合規(guī)范要求，否則要進行重新觀測。

2.3 網(wǎng)平差

基線解算得到的三維基線向量僅能確定GPS網(wǎng)的幾何形狀，但無法提供最終網(wǎng)中點的絕對坐標，通過引入起算點坐標可以達到引入絕對基準的目的。

(1)三維無約束平差

應將全部獨立基線構(gòu)成閉合圖形，以三維基線向量及其相應方差協(xié)方差陣作為觀測信息，在WGS-84坐標系中進行三維無約束平差，并提供WGS-84的三維坐標、觀測值的總改正數(shù)、基線邊長以及點位和邊長的精度信息。基線向量改正數(shù)的絕對值應滿足下列各式的要求：

v△x≤3σ

v△y≤3σ

v△z≤3σ

式中：σ——基線測量中誤差，單位為毫米(mm)

(2)二維約束平差

平差前，在青島城市坐標系下，對起算點進行兼容性檢查分析。分別用其中的若干個已知點計算其余點，將計算值和已知值進行比較，如果坐標分量互差不大于 2.0 cm，說明已知點間的兼容性良好，可以作為起算點，如果坐標分量互差大于 2.0 cm，則剔除該點。

將兼容性良好的點作為起算基準，計算出控制網(wǎng)點的青島城市坐標系坐標成果，并將其作為地面施工測量的控制點成果使用。約束平差后，輸出成果包括各點位的坐標、基線向量改正數(shù)、基線邊長和方位角、邊長和方位的精度信息等。

基線向量的改正數(shù)與同名基線無約束平差相應改正數(shù)的較差應滿足下列各式要求：

dv△x≤2σ

dv△y≤2σ

dv△z≤2σ

式中：σ——基線測量中誤差，單位為毫米(mm)

測量成果輸出宜包括相應坐標系中的三維或二維坐標、基線向量改正數(shù)、基線邊長、方位角、轉(zhuǎn)換參數(shù)及精度等信息。

3 青島地鐵實例分析

青島地鐵1號線某車輛段，由于跨度大、地域復雜，周邊僅有少數(shù)高等級城市控制點，需要在此基礎(chǔ)上布設(shè)衛(wèi)星定位控制點與已有高等級控制點進行連續(xù)觀測，如圖2所示。2017年1月23日，利用4臺徠卡GS14雙頻接收機進行觀測，共7個時段。外業(yè)測量過程嚴格執(zhí)行測量規(guī)范的要求，對外業(yè)觀測數(shù)據(jù)進行檢查，剔除不合格數(shù)據(jù)，確保內(nèi)業(yè)解算數(shù)據(jù)質(zhì)量。首先利用LGO進行基線解算，然后將原始數(shù)據(jù)通過TEQC轉(zhuǎn)換成RINEX標準格式用TBC進行基線解算，將兩種方式解算得到的基線向量分別利用武漢大學開發(fā)的CosaGPS進行網(wǎng)平差，平差過程除了使用的基線外，其余的參數(shù)設(shè)置完全相同。本文分析分為以下兩部分：①從不同的指標分析兩種隨機軟件基線解算效果；②對兩種隨機軟件基線解算后網(wǎng)平差精度進行評價。

圖2 GPS觀測網(wǎng)圖

(1)TBC和LGO基線解算結(jié)果對比

分別利用TBC和LGO進行基線解算，基線向量殘差如表1所示，用TBC解算得到的殘差絕對值與LGO解算得到的殘差絕對值的差值如圖3所示。

基線解算向量殘差 表1

圖3 基線向量殘差絕對值之差

圖4 基線向量殘差絕對值之差百分比統(tǒng)計

從表1、圖3、圖4可以看出，①對于同等條件下觀測得到的基線向量，用TBC解算得到的X殘差絕對值有78%小于LGO解算殘差絕對值，Y殘差絕對值有66%小于LGO解算殘差絕對值。②從基線向量殘差分析，總體上用TBC解算基線優(yōu)于LGO解算。

(2)兩種隨機軟件基線解算后網(wǎng)平差精度對比

分別利用TBC和LGO進行基線解算后得到的基線，選取相同的獨立基線參與二維約束平差，平差精度統(tǒng)計結(jié)果見表2；最弱點、最弱邊精度信息見表3、表4。

兩種隨即軟件基線解算后二維約束平差總體信息 表2

二維約束平差最弱點點位及其精度 表3

二維約束平差最弱邊及其精度 表4

從表2、表3、表4可以看出，①用TBC基線解算后網(wǎng)平差的驗后單位權(quán)中誤差為 3.06 mm，小于用LGO基線解算后網(wǎng)平差的驗后單位權(quán)中誤差 8.19 mm。②用TBC基線解算后網(wǎng)平差的最弱點點位中誤差為 2.5 mm，小于用LGO基線解算后網(wǎng)平差最弱點中誤差 3.3 mm。③用TBC基線解算后網(wǎng)平差的最弱邊邊長相對中誤差為1/49000，優(yōu)于用LGO基線解算后網(wǎng)平差最弱邊邊長相對中誤差1/46000。④用TBC解算基線后網(wǎng)平差精度整體優(yōu)于用LGO基線解算后網(wǎng)平差精度。

4 結(jié) 語

針對GPS作業(yè)過程中不同型號接收機的混合使用，相應配套的基線解算隨機軟件解算精度相差較大，本文考慮將不同類型接收機原始數(shù)據(jù)轉(zhuǎn)換成標準的RINEX格式后統(tǒng)一使用TBC進行解算。該軟件內(nèi)部算法嚴密，數(shù)據(jù)處理模型綜合考慮了大氣折射、對流層、電離層等外界影響并進行了相應的改正，因而基線解算結(jié)果較好。將本文的方法應用到青島地鐵1號線GPS靜態(tài)數(shù)據(jù)解算中，結(jié)果表明：在軟件參數(shù)設(shè)置、先驗精度等條件完全相同的情況下，用TBC基線解算后二維約束平差單位權(quán)中誤差為 3.06 mm、最弱點中誤差為 2.5 mm、最弱邊邊長相對中誤差1/49000，用LGO基線解算后二維約束平差單位權(quán)中誤差為 8.19 mm、最弱點中誤差為 3.3 mm、最弱邊邊長相對中誤差1/46000。因此，在GPS靜態(tài)數(shù)據(jù)解算中，可以將不同接收機數(shù)據(jù)轉(zhuǎn)換成標準RINEX格式，并用TBC進行基線解算。

由于本文只進行了LGO和TBC兩種隨機軟件之間的對比驗證，沒有與其他軟件進行對比。因此，后續(xù)工作將研究TBC軟件與其他隨機軟件之間解算效果的比較，以便找給出更佳的基線解算方案。

[1] 李征航，黃勁松. GPS測量與數(shù)據(jù)處理[M]. 武漢：武漢大學出版社，2005.

[2] CH2001-92. 全球定位系統(tǒng)(GPS)測量規(guī)范[S].

[3] 吳月霞. 以地鐵車站為核心的地下空間開發(fā)利用研究[D]. 上海：同濟大學，2008.

[4] 徐紹銓，李征航，李振洪等. 隔河巖大壩外觀變形GPS自動化監(jiān)測系統(tǒng)的建立[J]. 武漢測繪科技大學學報，1998，23(S).

[5] 黃丁發(fā)，熊永良，周樂韜等. GPS衛(wèi)星導航定位技術(shù)與方法[M]. 北京：科學出版社，2009.

[6] 黃勁松，李英冰. GPS測量與數(shù)據(jù)處理實習教程[M]. 武漢：武漢大學出版社，2010.

[7] 盧獻建，任超. GPS數(shù)據(jù)處理科研軟件與商業(yè)軟件對比分析[J]. 全球定位系統(tǒng)，2007，32(5)：29～32.

[8] 王潛心. GPS數(shù)據(jù)處理軟件中若干問題的研究[D]. 長沙：中南大學，2005.

[9] 王國祥. GPS基線解算與質(zhì)量控制[J]. 山東科技大學，2007：7～9.

[10] 王明孝，張之孔，馬懷武等. 幾種不同的GPS數(shù)據(jù)軟件的對比[J]. 全球定位系統(tǒng)，2012，37(1)：65～68.