袁 偉,鄧盛東,何為聰
(1.中國水電顧問集團昆明勘測設計研究院,云南昆明650041;2.江蘇南京地質工程勘察院蘇州分院,江蘇蘇州215000)
隨著空間定位技術的不斷發展和完善,GPS定位以其精度高、速度快、全天候,且不受天氣的干擾等優點已成為測繪領域中確定地面坐標的主要作業方法之一。由其定位方式可知,GPS在實現定位的過程中不可避免地受到多種誤差的影響,如:電離層延遲,對流層延遲,儀器故障、老化,測站的多路徑效應過高,測站附近電磁場干擾過強和暫時性的信息遮蔽等。在實現高精度定位的過程中,必須最大限度地減弱或消除各種因素對定位結果的影響。因此,有必要在GPS解算前對GPS觀測數據的質量進行檢核,對影響定位結果較大的觀測值應舍棄或者返工,這對觀測誤差的分析和精度的提高是有意義的。目前,評定觀測數據質量的方法主要有兩種,一種是根據商業軟件的基線解算結果分析;一種是使用專門的觀測值質量分析軟件。
TEQC(Translation Editing and Quality Checking)是功能強大且簡單易用的GPS/GLONASS/SBAS數據預處理及質量分析軟件。它的主要功能有格式轉換、編輯及質量檢核。本文中主要使用的是TEQC的質量檢核功能。
對TEQC檢核生成的文件實現可視化主要有兩種方式,一種是曲線視圖;一種是全星空視圖。曲線視圖查看的方式在UNAVCO網站上提供了與TEQC相輔助的工具,如:QCVIEW、gt、TEQCPLOT等。全星空視圖查看的方式目前主要使用QC2SKY。
電離層延遲、多路徑效應、接收機噪聲是影響數據質量的主要因素,其中尤以多路徑效應最為復雜,最為明顯。TEQC中電離層延遲、電離層延遲變化率分別采用式(1)和式(2)來計算,文獻[1]中詳細闡明了多路徑效應采用式(3)和式(4)的組合量來檢核其影響。


式中:P1,P2是雙頻偽距觀測值,L1,L2是雙頻載波相位觀測值,M1,M2是雙頻偽距多路徑效應,m1,m2是雙頻載波相位觀測值的多路徑效應,n1,n2是整周模糊度,f1,f2是頻率,λ1,λ2是波長。
根據有無導航文件,TEQC的質量檢核分為qc-lite和qc-full兩種模式。
(1)qc-lite(無導航文件)。
運行teqc+qc*.**o,TEQC對文件*.**o在qc-lite方式下進行質量檢核。質量檢核的結果會生成報告文件*.**S和數據文件*.ion(L1載波的電離層延遲)、*.iod(電離層延遲的變化率)、*.mp1(L1載波C/A碼或P碼偽距的多路徑效應影響)、*.mp2(L2載波P碼偽距的多路徑效應影響)、*.sn1(L1載波的信噪比)、*.sn2(L2載波的信噪比)。
(2)qc-full(輸入導航文件或同一目錄下有對應的N文件)。
同時輸入O文件和N文件,運行teqc+qc-nav*.**n*.**o;或N文件和O文件在同一目錄下,TEQC會自動搜索N文件,無需用-nav指定,運行teqc+qc*.**o,TEQC對文件*.**o在qc-full方式下進行質量檢核。檢核的結果除qc-lite方式下的報告文件和數據文件外,還增加了衛星和接收機天線的位置信息及兩個數據文件*.azi(衛星方位角),* .ele(衛星高度角)。
需要說明的是使用teqc中的參數++sym可以得到一個所有可能出現的標記符號及其注釋內容說明。其常用命令格式為:Teqc+qc++sym-nav* .**n* .**o。
數據質量評價的主要內容有:L1載波上的C/A碼或P碼的多路徑觀測誤差;L2載波上的P碼多路徑觀測誤差;電離層影響誤差;電離層延遲變化率。此外還包括接收機信號的抖動,L1、L2載波的失鎖(周跳),點位多路徑誤差平均值,衛星高度角和方位角,接收機時鐘漂移,信噪比和其它有關的參數。全星空圖查看使用QC2SKY軟件進行,其命令格式為:qc2sky+filename cutoff linewidth background。其 中:filename是TEQC的輸出文件,不帶擴展名;cutoff是截止高度角,單位是度;linewidth是線寬,數值從1至9;background是背景顏色。曲線圖的查看直接運行TEQCPLOT,選擇相應的文件即可。
實測數據選擇某一滑坡監測GPS控制網的觀測數據進行預處理與后處理,并與未經預處理模式下的結果進行比較。該GPS控制網由12個點組成,網形如圖1所示。

圖1 GPS控制網
S文件中衛星記錄情況見圖2。

圖2 衛星記錄情況
圖2 中:-dn行表示理論上可觀測到的衛星數和實際觀測到的衛星的最小差值;+dn行表示理論上可觀測到的衛星數和實際觀測到的衛星的最大差值;+10行表示高度角在10°以上時理論上可觀測到的衛星數;pos行表示定位情況;clk行表示接收機鐘差設置情況。
控制網中部分質量檢核結果視圖見圖3、圖4。
圖3的下邊緣為0°,每一格為10°。如果圖中出現綠色圓圈,則表示在該歷元有周跳發生。

圖3 JC04的全星空視圖

圖4 JC04的L1載波多路徑效應圖
分析質量檢核的結果,刪除觀測信號不好的時段及衛星。數據處理采用兩種方案進行。一種是對未經相關處理的數據進行基線解算和網平差;一種是對相關處理后的數據進行基線解算和網平差。兩種處理方案的設置完全一致,基線解算及網平差的軟件均使用TGO1.63版。
通過計算,未經相關處理的數據經過網平差后,最弱點的點位中誤差為±3.1 mm,相關處理后的數據經過網平差后,最弱點的點位中誤差為±2.4 mm。
采用兩種不同方案處理后各點的點位誤差變化見圖5~圖7。

圖5 X方向點位精度比較

圖6 Y方向點位精度比較

圖7 各點點位精度比較
由實例分析可見,經過質量檢核和相關處理后,平面點位精度較未經相關處理的點位精度有一定程度的提高,它彌補了使用常規商業軟件解算GPS測量數據時僅能通過基線解算后才能分析觀測數據質量的方式。
TEQC工具包功能強大、應用靈活,其相關的輔助軟件QC2SKY和TEQCPLOT可以有效地幫助用戶對質量檢查后的文件繪圖,從而對數據質量分析的結果進行直觀地查看,最大限度地避免了錯誤的發生,這無疑將會提高工作效率,減少資源浪費。
[1]http://igscb.jpl.nasa.gov/igscb/data/format/rinex210.txt[OL]
[2]http://www.unavcao.ucar.edu/software/teqc[OL]
[3]李洪玉,王耀強.利用TEQC對GPS數據質量的分析[J].測繪與空間地理信息,2009(4):80-42
[4]古偉洪,田鵬波.運用TEQC軟件對GPS數據的預處理與質量評定[J].地理空間信息,2008(12):37-39