TEQC軟件下的GPS數據質量控制

胡玉坤，劉根友，沙文東

(1. 寧波市測繪設計研究院，浙江 寧波 315042； 2. 中國科學院測量與地球物理研究所，湖北 武漢 430077；3. 武漢市測繪設計研究院，湖北 武漢 430022)

?

TEQC軟件下的GPS數據質量控制

胡玉坤1，劉根友2，沙文東3

(1. 寧波市測繪設計研究院，浙江 寧波 315042； 2. 中國科學院測量與地球物理研究所，湖北 武漢 430077；3. 武漢市測繪設計研究院，湖北 武漢 430022)

利用基于Matlab二次開發的SPEC模塊在Windows下對TEQC質量檢核生成的繪圖文件進行直觀顯示，便于了解測站觀測環境以剔除不良時間段和衛星的觀測數據。對剔除前后的數據利用自編軟件進行基線解算，試驗表明，利用TEQC處理后的基線精度有了一定程度的提高。

Matlab二次開發；SPEC模塊；TEQC質量檢核；基線解算

GPS定位精度取決于原始觀測數據的質量和算法的選擇，因此在進行數據預處理之前，應當對GPS數據進行質量檢查。通常評定觀測數據的質量有兩種方法：一是根據商業軟件的基線解算結果進行分析；二是采用專業的質量評定軟件TEQC[1]。TEQC(Translate/Edit/Quality Check/Coordinate)是由UNAVCO Facility 研制的為地學研究GPS監測站數據管理服務的公開免費軟件，主要功能有格式轉換、數據編輯和質量檢核,可用于檢查雙頻GPS接收機的動態和靜態數據質量。其簡短的幾條命令便可對觀測數據的信噪比、周跳、多路徑呈現出一個直觀的顯示，極大地方便了觀測數據中粗差的剔除[2-5]。本文主要利用基于Matlab開發的SPEC模塊對TEQC處理文件直觀顯示，并利用自編軟件對處理前后的基線進行解算來驗證處理前后數據的質量。

一、TEQC質量檢核原理與可視化分析

1. TEQC質量檢核原理

TEQC的質量檢核是通過偽距和載波觀測值的線性組合來計算L1偽距C/A碼或P碼的多路徑、L2偽距P碼多路徑、電離層延遲和電離層延遲的變化率，以及接收機的鐘漂和信噪比等信息。根據是否使用導航電文，可以分為qc-lite和qc-full模式。在qc-full模式下可以生成1個質量匯總文件和8個繪圖文件。

雙頻測相偽距和碼偽距[6]如下

(1)

L1-L2=(α-1)I+λ1n1-λ2n2+m1-m2

(2)

由式(2)中的第一式與第三式可求得電電離層延遲I

I=(P1-L1+λ1n1+m1-M1)/2

(3)

將式(3)代入式(2)可得

(4)

記式(4)左邊為mp1,右邊忽略相位多路徑可得

(5)

(6)

(7)

對于質量匯總文件，幾個數據質量的關鍵性指標(多路徑、數據完整性、周跳比及信噪比)都在其中。mp1和mp2是隨著時間變化的，與Mi相比mi幾乎可以忽略。可將mp1和mp2直接視為偽距多路徑的影響，習慣將其在一定時間窗口內的均方根作為反映觀測數據質量的指標。IGS數據質量檢測表明，mp1平均值小于0.5，mp2平均值小于0.75。因此本文對于mp1和mp2的閾值分別取0.5和0.75。總的來說，mp1和mp2越小，說明抗多路徑效應的能力越強，數據質量越好。數據完整性為預期歷元數與實際觀測歷元數的比值，O/Slips為觀測值個數與周跳的比值，也可用CSR=1000/O/Slips表示。信噪比(SNR值)是指接收機載波信號強度與噪聲強度的比值。信噪比值越高，觀測數據的質量越好。文件中以5°為區間列出了10°—90°高度角范圍內所有歷元信噪比的平均值。

2. 繪圖工具及其二次開發

目前能夠繪制TEQC生成的繪圖文件的軟件大都是基于DOS的，人機交互性差。包括QCVIEW32、QC2SKY和CF2PS等[7-9]，可以顯示單個衛星或分色顯示所有衛星的數據隨時間變化的曲線圖。本文利用Matlab軟件，通過二次開發的SPEC程序模塊，可以實現在Windows界面下直接打開TEQC并直觀地顯示其繪圖文件，根據用戶需求的不同，采用零相位數字濾波對多路徑效應進行濾波，繪制多路徑天空圖[10-13]。

二、GPS數據質量分析

對處于樓頂開闊視野下的一條短基線進行24 h連續觀測(如圖1所示)，為了便于橫向對比，可將兩接收機分別安裝不同的天線(TRIM59800.00和安騰科技公司生產的天線)。用TEQC處理觀測數據，通過上文所設定的標準剔除觀測質量較差的時間段及衛星數據。在剔除的過程中會用到TEQC的數據編輯功能，如RINEX文件的切割、特定衛星的禁用等。為了更好地對比TEQC處理前后GPS數據質量的變化。對TEQC處理前后的數據采用自編軟件解算[14]。

圖1　現場基線布設圖

采用2015年6月5日至6月6日的觀測數據，觀測文件分別為Trim.15o、Anteng.15o。在IGS網站上下載相應的廣播星歷，qc-full模式下即可得到質量匯總文件*.15S和8個繪圖文件*.azi、*.ele、*.iod、*.ion、*.mp1、*.mp2、*.sn1、*.sn2。打開Trim.15S和Anteng.15S可以看到觀測時段、多路徑效應、信噪比、周跳比及衛星高度角的變化等詳細信息。

由表1可以看出，由于觀測時間較長且觀測環境較好，因此多路徑效應和周跳比皆小于選取的閾值。由于是短基線(10m)，故可視A、B兩點處于同等觀測條件，B點的多路徑效應比A點的大，周跳比A點較B點的大，說明了A點觀測數據的質量優于B點。

表1　mp1、mp2、O/Slips統計值

運用基于Matlab開發的SPEC程序模塊繪制多路徑效應隨時間的變化圖及信噪比隨高度角的變化圖，檢測出質量較差的衛星數據及其對應時段。圖2為原始觀測數據每顆衛星L1與L2載波的多路徑效應隨時間變化圖，圖3為原始觀測數據L1與L2載波的信噪比隨高度角的變化圖。

圖2　觀測時段內每顆衛星的多路徑效應

圖3　信噪比隨高度角的變化

從圖2可以看出，PRN04、PRN07的多路徑效應比較顯著，15:24至16：10和1：12至3.30時間段的多路徑效應比較大。從圖3可以看出，L1的信噪比要明顯優于L2的信噪比，并且隨著高度角的增加，信噪比逐漸增大。利用TEQC中的數據編輯功能對上述衛星及觀測時段的數據進行剔除。輸入指令：>teqc-G04,07 Anteng.15o>Anteng1.15o可得到剔除PRN04、PRN07衛星的觀測文件Anteng1.15o。對異常時間段的剔除可以分段處理，輸入指令：>teqc-st 20150605080000-e 20150605152400 Anteng1.15o>Anteng2.15o可得到起始時刻至15：24的觀測文件Anteng2.15o，同理可以獲得16：10 至次日1：12的觀測文件Anteng3.15o和3：30至終止時刻的觀測文件Anteng4.15o。利用TEQC對上述3個文件進行質量檢核，其mp1、mp2及周跳比見表2。

表2　mp1、mp2、O/Slips統計值

從表2可以看出，相對于表1中Anteng.15o文件的mp1、mp2及O/Slips，剔除后3個觀測時段文件皆有提高。其中mp1、mp2平均提高了10%左右，O/Slips提高了120%左右。說明在多路徑效應比較顯著的地方，周跳現象比較嚴重，這為后續聯合多路徑消除周跳工作開拓了思路。

利用自編軟件對原始短基線和處理后形成的3個時段的短基線進行解算，為了進一步分析mp1、mp2對解算精度的影響，可以分別采用L1模式、L2模式及LC模式解算。自編基線解算軟件界面如圖4所示。解算精度見表3。

圖4　自編基線解算軟件界面

基線L1模式L2模式LC模式原始0．00410．00470．0087時段10．00370．00350．0068時段20．00390．00410．0076時段30．00380．00330．0069

從表3可以看出，無論是L1模式、L2模式還是LC模式，基線的解算精度都有所提高。說明利用TEQC軟件對觀測數據進行質量分析，并剔除不良觀測時段與衛星，可在一定程度上提高基線結算的精度。

三、結束語

利用TEQC軟件對GPS觀測數據進行質量檢核，操作簡單，能夠快速生成反映GPS數據質量的一系列文件，通過基于Matlab二次開發的SPEC模塊可以在Windows下直觀地顯示GPS觀測數據的質量，以便于用戶了解周邊的觀測環境并剔除不良觀測時段及衛星數據。試驗表明，對剔除前后的基線利用自編軟件解算后，其精度有了一定程度的提高，這對利用TEQC軟件處理GPS數據具有重要的實踐意義。

[1]ESTEY L H, MEERTENS C M. TEQC: The Multi-purpose Toolkit for GPS/GLONASS Data[J]. GPS Solutions, 1999, 3(1): 42-49.[2]李軍, 王繼業, 熊熊, 等. 東北亞地區 GPS 觀測數據質量檢測和分析[J]. 武漢大學學報(信息科學版), 2006, 31(3): 209-212.

[3]古偉洪, 田鵬波, 王振輝. 運用 TEQC 軟件對 GPS 數據的預處理與質量評定[J]. 地理空間信息, 2008, 6(6): 37-39.

[4]魏二虎, 王中平, 龔真春, 等. TEQC 軟件用于 GPS 控制網數據質量檢測的研究[J]. 測繪通報, 2008 (9): 6-9.

[5]YUAN N I E G W. GPS Data Assessment by TEQC Program and Software Menu Development [J]. Gnss World of China, 2003(4)：8.

[6]張勤, 李家權. 全球定位系統 (GPS) 測量原理及其數據處理基礎[M]. 西安：西安地圖出版社, 2001.

[7]趙紅強, 郭際明, 李宗華, 等. 繪圖工具 QCVIEW 在 GPS 數據預處理中的應用[J]. 測繪通報, 2008 (11): 44-46.

[8]王建軍, 吳北平, 馮健. TEQC 和 QCVIEW32 在 GPS 觀測數據預處理中的應用[J]. 工程地球物理學報, 2010, 7(5): 634-638.

[9]劉剛, 李征航, 于曉歆. TEQC 與 CF2PS 在 GPS 數據預處理中的應用[J]. 地理空間信息, 2010 (5): 146-148.

[10]賈瑩媛, 黃張裕, 楊富春, 等. 基于Matlab 的 TEQC 在 GPS 數據預處理中的應用[J]. 海洋測繪, 2012, 32(2): 52-54.

[11]余文坤, 戴吾蛟, 楊哲. 基于 TEQC 的 GNSS 數據質量分析及預處理軟件的設計與實現[J]. 大地測量與地球動力學, 2010, 30(5): 81-85.

[12]聶桂根, 王院. TEQC 軟件對 GPS 數據質量的評定及軟件界面的開發[J]. 全球定位系統, 2003, 28(4): 32-38.

[13]HILLA S. A New Plotting Program for Windows-based TEQC Users[J]. GPS Solutions, 2002, 6(3): 196-200.

[14]陳凱華, 祿占磊. TGO 對基于 TEQC 處理的 GPS 數據質量評定[J]. 河南工程學院學報(自然科學版), 2011, 23(1): 53-55.

GPS Data Quality Control Based on TEQC

HU Yukun,LIU Genyou,SHA Wendong

胡玉坤，劉根友，沙文東.TEQC軟件下的GPS數據質量控制[J].測繪通報，2016(10)：51-53.DOI:10.13474/j.cnki.11-2246.2016.0327.

2015-11-13；

2016-05-18

國家自然科學基金(41321063)；大地測量與地球動力學國家重點實驗室基金(SKLGED2013-4-1-Z)

胡玉坤(1990—)，男，碩士，研究方向為GNSS衛星導航與定位。E-mail：huyukun0608@163.com

P228.4

B

0494-0911(2016)10-0051-03