李 立 蘭武輝 吳潔琪

（江西省測繪成果質(zhì)量監(jiān)督檢驗測試中心 江西南昌 330046）

基于TEQC的GPS觀測數(shù)據(jù)預(yù)處理

李 立 蘭武輝 吳潔琪

（江西省測繪成果質(zhì)量監(jiān)督檢驗測試中心 江西南昌 330046）

GPS外業(yè)觀測數(shù)據(jù)的預(yù)處理與質(zhì)量分析對所得的成果起著重要的保證作用。在基于TEQC軟件的基礎(chǔ)上，對GPS外業(yè)觀測數(shù)據(jù)的預(yù)處理與質(zhì)量分析進(jìn)行闡述，以保證GPS數(shù)據(jù)采集成果滿足質(zhì)量要求。

數(shù)據(jù)預(yù)處理；質(zhì)量分析；TEQC軟件

1 引言

GPS數(shù)據(jù)預(yù)處理與質(zhì)量分析的結(jié)果是測量工作者密切關(guān)注的問題，因而采用何種數(shù)據(jù)預(yù)處理與質(zhì)量分析軟件一直是國內(nèi)外學(xué)者爭相研究的課題?？傮w來說，數(shù)據(jù)預(yù)處理與質(zhì)量分析軟件有TGO軟件、Solutions軟件等。UNAVCO組織開發(fā)的TEQC軟件已經(jīng)被大多數(shù)用戶所采用。

2 TEQC軟件功能介紹

TEQC是一個命令行運(yùn)行軟件，運(yùn)行環(huán)境為DOS或Windows Dos。對應(yīng)于不同的接收器和不同類型的數(shù)據(jù)格式，使用的命令是不同的。例如：

（1）將Ashtech觀測的文件轉(zhuǎn)換為RINEX格式文件，使用命令為：

＞teqc–ash d–O.obs L1+L2+CA+P2+ P1 BB003A12.038＞B0030380.120

（2）將Trimble觀測的文件轉(zhuǎn)換為RINEX格式文件，使用命令為：

＞teqc–tr do–week****+nav 00010160.02n 00010160.dat＞00010160.02o

2.1 文件頭部分的編輯

TEQC軟件可對RINEX觀測數(shù)據(jù)文件的頭部分進(jìn)行編輯，命令如下：

＞teqc–O.mo“A001”kunm1360.12o＞A0011360.12o

2.2 RINEX文件的切割

TEQC軟件可對RINEX文件進(jìn)行任意的切割，以滿足用戶對特定觀測時間段內(nèi)的數(shù)據(jù)使用要求。命令有如下幾種：

（1）+d［Y、M、d、h、m、s］，表示從觀測時間開始到設(shè)置的時間間隔結(jié)束這段時間的觀測數(shù)據(jù)文件。例如：＞teqc+dm 20 kunm1360.12o＞temp1360.12o，這個命令表示取RINEX文件從觀測時間開始到20min后的時間段的數(shù)據(jù)；

（2）-d［Y、M、d、h、m、s］，表示以結(jié)束觀測為下限所提取的時間間隔；

（3）-st YYMMddhhmmss［.sss］，表示從此刻起到觀測結(jié)束這段時間段內(nèi)的觀測數(shù)據(jù)；

（4）-e YYMMddhhmmss［.sss］，表示從開始觀測到此刻結(jié)束這段時間內(nèi)的觀測數(shù)據(jù)。

2.3 RINEX文件的合并

TEQC軟件可對多個RINEX文件進(jìn)行合并，以此優(yōu)化存儲空間。命令如下：

＞teqc file1 file2……＞file

2.4 衛(wèi)星系統(tǒng)的選擇和取舍

在對GPS數(shù)據(jù)進(jìn)行預(yù)處理時，會發(fā)現(xiàn)GPS接收機(jī)接受某些衛(wèi)星的信號質(zhì)量不佳，這時就需要用TEQC軟件對衛(wèi)星系統(tǒng)進(jìn)行選擇和取舍，以使成果達(dá)到要求的精度。相關(guān)的TEQC命令如下：

（1）＞teqc–R輸入文件 ＞輸出文件，此命令表示去掉GLONASS衛(wèi)星數(shù)據(jù)；

（2）＞teqc–G#輸入文件 ＞輸出文件，此命令表示禁用#號GPS衛(wèi)星。

2.5 設(shè)置衛(wèi)星高度角

在考慮對流層和電離層對GPS觀測造成的影響時，需要設(shè)置衛(wèi)星高度角。命令為：＞teqc set_mask #輸入文件＞輸出文件

2.6 質(zhì)量檢核

數(shù)據(jù)質(zhì)量檢核在GPS實際工作中十分重要，電離層、對流層、多路徑效應(yīng)等對GPS觀測數(shù)據(jù)造成的影響不容忽視。TEQC軟件對GPS觀測數(shù)據(jù)的質(zhì)量檢核有兩種模式：qc-lite模式和qc-full模式。

(1)qc-lite模式。此模式只需RINEX觀測數(shù)據(jù)文件， 不需導(dǎo)航文件， 命令為：＞teqc+qc kunm1360.12o。此命令運(yùn)行完畢后，會產(chǎn)生7個文件 ：kunm1360.12S、kunm1360.mp1、kunm1360.mp2、kunm1360.iod、kunm1360.ion、kunm1360.sn1、kunm1360.sn2。

(2)qc-full模式。此模式不僅需要RINEX觀測數(shù)據(jù)文件，還需要RINEX導(dǎo)航數(shù)據(jù)文件，命令為：＞teqc +qc–nav kunm1360.12n kunm1360.12o。此命令運(yùn)行完畢后， 會產(chǎn)生 9個文件：kunm1360.12S、kunm1360.mp1、kunm1360.mp2、kunm1360.iod、kunm1360.ion、kunm1360.sn1、kunm1360.sn2、kunm1360.azi、kunm1360.ele。

3 實例

3.1 操作步驟

本次實例中采用的數(shù)據(jù)源自中國地殼運(yùn)動觀測網(wǎng)絡(luò)中的GPS數(shù)據(jù)，相關(guān)數(shù)據(jù)的下載鏈接為：http:// www.igs.org.cn:8080/，在此采用上海佘山、吉林長春、臺 灣桃 園 的 監(jiān) 測 點 文 件， 即 shao0880.13o、shao0880.13n;chan0880.13o、chan0880.13n;twtf0880.13o、twtf0880.13n。由于下載了觀測數(shù)據(jù)和導(dǎo)航數(shù)據(jù)文件，因此采用qc-full模式進(jìn)行數(shù)據(jù)的質(zhì)量檢核。步驟如下：

（1）在開始程序中的運(yùn)行功能里輸入cmd，進(jìn)入到DOS命令編輯區(qū)，如圖3-1所示：

圖3-1

（2）輸入DOS命令進(jìn)入到C盤根目錄，找到相應(yīng)的 teqc文件，在這個文件夾里面有 teqc.exe、DOS4GW.exe、qcview32.exe以及相關(guān)點的文件，如圖3-2所示：

圖3-2

（3）運(yùn)行 teqc.exe，并啟用 qc-full模式對shao0880.13o的相關(guān)數(shù)據(jù)進(jìn)行質(zhì)量分析，如圖3-3所示：

圖3-3

（4）運(yùn)行完qc-full模式之后，會在本身目錄產(chǎn)生9個新文件，將這9個新文件放在一個bk1文件夾里，如圖3-4所示：

圖3-4

（5）運(yùn)用 QCVIEW32查看各衛(wèi)星的方位角（azi）、高度角（ele）、電離層延遲變化率（iod）、電離層延遲誤差（ion）、L1載波信噪比（sn1）、L2載波信噪比 （sn2）、L1載波C/A碼或P碼偽距的多路徑影響（mp1）、L2載波P碼

偽距的多路徑影響（mp2）。由于其他兩個點相應(yīng)數(shù)據(jù)的處理與此相似，在這里不再重復(fù)敘述步驟。接下來，以17號衛(wèi)星為例來說明數(shù)據(jù)的質(zhì)量和對數(shù)據(jù)進(jìn)行處理。

3.2 視圖文件分析

（1）從圖3-5可以看出，17號衛(wèi)星在第1820個歷元（對應(yīng)時刻15h10min）處和第2550個歷元（對應(yīng)時刻21h15min）處出現(xiàn)了周跳現(xiàn)象；

圖3-5 17號衛(wèi)星處理前電離層延遲誤差

（2）從圖3-6可以看出，17號衛(wèi)星在第2550個歷元（對應(yīng)時刻21h15min）處電離層的變化率出現(xiàn)了突然的變化；

圖3-6 17號衛(wèi)星處理前電離層延遲變化率

（3）從圖3-7可以看出，17號衛(wèi)星在第1820個歷元（對應(yīng)時刻15h10min）處和第2680個歷元（對應(yīng)時刻21h15min）處L1載波的多路徑誤差的絕對值已經(jīng)超過1m；

圖3-7 17號衛(wèi)星處理前L1載波多路徑誤差

（4）從圖3-8可以看出，17號衛(wèi)星在第1820個歷元（對應(yīng)時刻15h10min）處和第2680個歷元（對應(yīng)時刻21h15min）處L2載波的多路徑誤差的絕對值已經(jīng)超過1m；

圖3-8 17號衛(wèi)星處理前L2載波多路徑誤差

（5）從圖3-9可以看出，17號衛(wèi)星L1載波的信噪比偏高；

圖3-9 17號衛(wèi)星處理前L1載波信噪比

（6）從圖3-10可以看出，17號衛(wèi)星L2載波的信噪比偏高。

圖3-10 17號衛(wèi)星處理前L2載波信噪比

3.3 數(shù)據(jù)預(yù)處理

通過視圖的分析，進(jìn)一步了解了觀測數(shù)據(jù)中存在的問題，現(xiàn)利用TEQC命令對數(shù)據(jù)進(jìn)行編輯，步驟如下：

（1）刪除shao0880.13o文件中第1850歷元以前的數(shù)據(jù)，命令為：teqc-st20130329152500 shao0880.13o＞shao0881.13o；

（2）刪除shao0881.13o文件中第2550歷元以后的數(shù)據(jù)，命令為：teqc-e20130329211500 shao0881.13o＞shao0882.13o；

（3）刪除shao0882.13o文件中質(zhì)量不好的3號、26號衛(wèi)星， 命令為：teqc-G3 shao0882.13o＞shao0883.13oteqc-G26shao0883.13o＞shao0884.13o；

（4）用TEQC軟件對shao0884.13o進(jìn)行質(zhì)量分析，并用QXCIEW32查看處理后的文件。

3.4 處理前后的數(shù)據(jù)在SPSO軟件中基線解算質(zhì)量的對比

圖4-11 數(shù)據(jù)處理前基線解算報告

圖4-12 數(shù)據(jù)處理后基線解算報告

4 結(jié)束語

通過基于TEQC的GPS觀測數(shù)據(jù)預(yù)處理，可以看出：

（1）采用TEQC軟件對GPS觀測數(shù)據(jù)預(yù)處理，可以發(fā)現(xiàn)觀測值中存在的誤差，有利于我們及時糾正誤差，保證后期觀測數(shù)據(jù)平差時的精度。

（2）各種GPS觀測數(shù)據(jù)的誤差與衛(wèi)星高度角有很大的聯(lián)系。隨著衛(wèi)星高度角的降低，信噪比也會降低，多路徑效應(yīng)、電離層影響都會加劇。因此，應(yīng)當(dāng)取適當(dāng)?shù)母叨冉刂菇?，盡量接收高緯度的衛(wèi)星，從而保證GPS觀測數(shù)據(jù)的質(zhì)量。

（3）剛開始GPS的數(shù)據(jù)觀測時，此時接收機(jī)所跟蹤的衛(wèi)星信號受各種干擾嚴(yán)重，因此最好剔除最初的歷元數(shù)據(jù)。

（4）根據(jù)檢查結(jié)果對GPS數(shù)據(jù)進(jìn)行相應(yīng)的預(yù)處理，通過數(shù)據(jù)處理前后的對比分析可以看出經(jīng)TEQC預(yù)處理后的GPS觀測數(shù)據(jù)，有效的提高了GPS基線解算的質(zhì)量。

［1］劉基余.全球定位系統(tǒng)原理與應(yīng)用［M］.北京：測繪出版社，1993.

［2］周忠謨，易杰軍，周琪.GPS衛(wèi)星測量原理與應(yīng)用［M］.北京：測繪出版社，2004.

［3］魏子卿，葛茂榮.GPS相對定位的數(shù)學(xué)模型［M］.北京：測繪出版社，1997.

［4］李征航，黃勁松.GPS測量與數(shù)據(jù)處理［M］.武漢：武漢大學(xué)出版社，2005.

［5］范世杰，郭際明，彭秀英.TEQC在GPS數(shù)據(jù)預(yù)處理中的應(yīng)用與分析［J］.測繪信息與工程，2004，4(2):33-35.

［6］李輝，李建文，高海榮.三級GPS大地控制網(wǎng)觀測數(shù)據(jù)質(zhì)量檢測［J］.測繪科學(xué)技術(shù)學(xué)報，2007，4(2):139-141.

［7］李軍，王繼業(yè)，熊熊等.東北亞地區(qū) GPS觀測數(shù)據(jù)質(zhì)量檢測與分析［J］.武漢大學(xué)學(xué)報，2006，3(3):209-212.