GPS靜態觀測中數據質量實時檢測

盧獻健 晏紅波 黃 鷹

1 桂林理工大學測繪地理信息學院,桂林市雁山街319號，541006 2 廣西空間信息與測繪重點實驗室，桂林市雁山街319號，541006 3 桂林好測科技信息有限公司，桂林市雁山街319號，541006

?

GPS靜態觀測中數據質量實時檢測

盧獻健1，2晏紅波1，2黃 鷹3

1 桂林理工大學測繪地理信息學院,桂林市雁山街319號，541006 2 廣西空間信息與測繪重點實驗室，桂林市雁山街319號，541006 3 桂林好測科技信息有限公司，桂林市雁山街319號，541006

在分析觀測數據解碼以及觀測數據質量指標的基礎上編制軟件，對采集到的原始GPS觀測數據進行解碼與多路徑效應、周跳等質量指標的分析計算，并對質量指標進行直觀、實時的圖形顯示。介紹所寫軟件在工程實踐中的應用，并與TEQC軟件進行比較，證明了軟件的正確性和可靠性。

靜態觀測；數據質量；實時檢測；數據解碼

衛星靜態觀測數據的質量對提高定位精度和工作效率有重要意義[1]。而在實際工作中有時因某個測站的觀測數據質量太差而不得不重測，影響工作進度和效率。觀測中若能及時進行觀測數據的質量檢測發現異常情況，則可以采取措施避免因觀測數據質量不合格產生的不利影響，提高效率并保證觀測成果的可靠性[2]。

GPS觀測與常規觀測相比，有其自身的特殊性：采用電子記錄、數據量大、受觀測的時間和地點影響較大。雖然有些接收機能顯示觀測時的衛星個數、分布、信噪比等信息，但是僅有這些信息難以準確反映觀測值質量，而且觀測者難以跟蹤記錄這些信息的變化。采用衛星預報軟件雖然可以預報觀測點的衛星情況，但是無法準確反映多路徑效應和周跳個數等信息。TEQC軟件在GNSS觀測數據質量的檢測中得到廣泛的應用[1,3-5]，但其適用于事后質量檢測，無法達到數據質量實時檢測的目的。就目前而言，要進行觀測數據質量的實時檢測，需要完成接收機與計算機間的數據傳輸、數據解碼及數據質量分析計算等工作[6-7]。

本文以使用Trimble BD970 OEM板為主板的接收機作為實驗設備，利用VB.NET語言編寫程序，實現接收機與計算機的實時通訊及二進制數據的解碼，在此基礎上進行周跳探測、多路徑效應及衛星高度角等的計算，并用圖形顯示的方式對數據質量指標進行可視化。最后介紹了本文方法在工程實踐中的應用及與TEQC軟件的比較，驗證了本文方法及程序的正確性。

1 觀測數據傳輸與解碼

GPS接收機一般采用二進制格式記錄原始觀測數據，而衛星觀測數據處理中的標準數據格式為Rinex格式，是ASCⅡ碼文本格式。因此，為了便于數據處理，首先要進行數據的傳輸和解碼。由于各接收機生產商的二進制格式不盡相同，數據解碼的具體操作略有不同，但基本步驟大致相同：1)明確變量的含義；2)獲取變量值；3)數值換算。對于Trimble BD970 OEM板卡，利用VB.NET的MSCOMM控件可以簡便、快速地建立接收機與計算機間的數據通訊，實現穩定、快速的數據傳輸。本文方法中通過定時觸發MSCOMM控件的Timer事件來實現接收機與計算機間的數據實時傳輸。數據解碼方法可以參考文獻[8]。

2 數據質量指標

觀測數據質量的評價指標有多路徑效應大小、周跳個數、數據有效率和信噪比等。

2.1 多路徑效應影響的計算

利用測碼偽距與載波相位觀測量來計算測站多路徑效應的綜合影響[9]：

(1)

(2)

式中，P1、P2分別為L1、L2兩波段上的測碼偽距觀測值；φ1、φ2分別為兩個波段上的載波相位觀測值；κ為L1、L2兩波段頻率f1、f2之比的平方。多路徑效應的評價指標采用均方差(RMS)表示。

2.2 周跳探測

電離層殘差法[10]可依據前幾個歷元的數據對下一個歷元是否發生周跳進行探測，因此適用于本文實時檢測的目的。

設某一歷元t電離層殘差組合為：

(3)

式中，λ、φ、N和f分別為載波的波長、載波相位觀測值、整周模糊度和載波頻率；下標1、2為L1、L2載波；A(t)=-40.3∫Neds，Ne為信號路徑上的電子密度，s為傳播路徑。將式(3)兩邊同時除以λ1并在相鄰歷元間求差，得到電離層殘差法的檢驗量：

(4)

式中，ΔNI(t)表示歷元間電離層殘差的變化值。

當L1、L2載波相位觀測值沒有發生周跳時，ΔN1、ΔN2等于0，則相鄰歷元的電離層殘差值ΔNI(t)不會發生跳變。因此，可以通過ΔNI(t)的變化情況判斷是否發生周跳。本文使用多項式擬合及電離層殘差法共同探測周跳。

實際工作中，常采用總觀測值數量和周跳發生數量的比值來反映整個觀測過程的數據質量，比值越大說明數據中發生周跳的次數越少。其表示為：

(5)

2.3 數據有效率

數據采集中，當接收機(雙頻、單頻)對衛星完成鎖定后，則正常情況下某一歷元所能獲得的觀測值(L1、L2、P1、P2、C1、C2等)個數是已知的。若測站周圍環境變化、接收機出現異常等原因引起衛星信號失鎖，則影響某時刻接收機對某一類型觀測值的采集，造成實際采集數據量小于應采集數據量，或者所采集的觀測數據誤差較大。因此，對于某一時段來說，實際采集數據量與應采集數據量的比值亦可反映觀測數據的質量和測站觀測環境的好壞。該比值稱為觀測數據有效率，表示為：

(6)

3 系統設計與靜態數據實時檢測

3.1 實時檢測功能模塊的設計與實現

在上述分析的基礎上，利用VB.NET編寫數據質量的實時檢測軟件(real-time QC，RTQC)，對使用了Trimble BD970 OEM主板的某國產GPS接收機的原始二進制觀測數據進行實時解碼。對解碼后的衛星觀測數據，利用電離層殘差法進行周跳探測并計算多路徑效應均方差MP1、MP2。最后將周跳探測結果(周跳發生個數與時刻)和MP1、MP2值可視化，利用圖形進行實時動態顯示。同時還可以對觀測時段內已采集數據進行實時的質量統計分析，及時掌握數據質量的變化。實時檢測的流程見圖1。

圖1 數據質量實時檢測的流程Fig.1 Process of real-time data quality detection

為了提高數據存儲、訪問等操作的效率，數據質量檢測的結果采用二維數據表DataTable來存放，周跳信息、MP1與MP2分別存放于單獨的文件，其格式見表1。表中第一列表示數據采樣間隔，其他各列表示對應衛星的數據檢測質量結果。

表1 MP1檢測結果二維數據表格式

數據質量實時檢測的界面見圖2。

圖2 實時質量檢測的界面Fig.2 The program interface of real-time data quality detection

RTQC軟件的主要功能有數據傳輸、質量分析與質量信息統計。數據傳輸采用串口通訊的方式實現接收機與計算機的實時數據通訊；質量分析主要包括觀測數據的多路徑效應(MP1、MP2)與周跳兩項指標的圖形可視化顯示；質量信息統計可以計算觀測時段內任意時刻的數據有效率、o/slips指標及平均多路徑效應等數據質量統計信息；圖形顯示區可以顯示全部衛星和單顆衛星的質量指標，方便進行數據質量分析，有關圖形顯示的參數可以在菜單“視圖”中設置，如時段長，采樣間隔等繪圖參數。

3.2 靜態數據的實時檢測

為驗證本文RTQC軟件對觀測數據質量實時檢測的效果，設計如下實驗：選擇具有不同觀測環境的A、B兩站點，利用RTQC軟件及同一臺接收機分別進行數據質量的實時檢測，其中A站點為某大樓樓頂，具有良好的凈空，觀測條件良好；B站點為大樓下的地面點，距大樓較近，并且其周圍有較高的樹木，觀測條件較差。實驗觀測時段為2 h，數據采集間隔為15 s。觀測開始前，根據時段長和采樣間隔設置繪圖參數，觀測過程中RTQC軟件可以實時進行數據質量檢測。觀測完成后，還可利用TEQC軟件分別對所采集的數據進行質量檢測，以驗證RTQC軟件的結果。

圖3為A、B兩站點前10 min靜態觀測的數據質量檢測結果。

圖3 A、B兩站短時間內檢測結果Fig.3 Real-time detection results on station A and B

圖3(a)、圖3(b)分別為A站點約10 min所采集數據的MP1、MP2檢測結果。不難看出，由于A點的觀測條件良好，因此其MP1、MP2的檢測結果較小，MP1、MP2的最大值均沒有超過1 m，MP2的數值較MP1稍大。而圖3(c)、圖3(d)為觀測條件相對較差的B點的MP1、MP2檢測結果。由于其附近有大樓和樹木等因素的影響，圖中所示的多路徑效應的檢測結果較大，最大的接近5.0 m和10 m。另外，通過軟件所顯示的單顆衛星多路徑效應圖形可知，未發生數值突變的部分衛星的MP1、MP2數值也超過了2.0 m，而且還有數據中斷的情況，說明B點的觀測確實受到了其周圍外界因素的影響。因此，相對A點的檢測結果，B點有多顆衛星的MP1、MP2的數值偏大，其受周圍觀測環境的影響較大。

由上述的實驗分析可知，利用RTQC軟件對具有不同觀測條件的測站進行數據質量檢測，其檢測結果符合事前的預期，說明RTQC軟件可以在較短的時間內(幾分鐘)通過對所采集的數據進行實時質量分析，及時反映觀測條件的好壞。

圖4為完成2 h觀測后，A、B兩點整時段的MP1的檢測結果。在完成整時段的數據質量檢測后，RTQC軟件還可以對該時段的數據質量作統計分析，如圖5為A點靜態觀測數據質量的統計信息。

圖4 A、B兩點整時段的MP1檢測結果Fig.4 Detection results of MP1 in observation interval

圖5 A站點觀測質量統計Fig.5 Observation quality statistics results of station A

從圖5可知，A站點各顆衛星的多路徑效應數值穩定，數值變化基本在-0.75～0.75。借助單顆衛星的MP1變化圖不難得出，G12號衛星的MP1值在270～350歷元之間有大于1的間斷變化，但總體上其數值處于較小的浮動區間內。相對而言，B站點的MP1結果有明顯的跳變和數據中斷，發生在多顆衛星及不同的時間段上。

綜合B站點的MP2與周跳檢測結果可知，B站點周圍的環境對觀測造成了影響，觀測數據質量明顯不如A站點，這符合觀測前對站點觀測質量的預測，同時說明本文RTQC軟件的正確性。

3.3 與TEQC軟件的分析比較

為了進一步驗證RTQC軟件的正確性和上述實時數據質量檢測的可靠性，利用該接收機的隨機軟件將觀測數據的格式轉為Rinex格式，并用TEQC軟件對A、B測站點所采集的數據進行事后數據質量檢測。TEQC與RTQC軟件數據質量檢測結果的對比見表2。

表2 兩種軟件數據檢測結果的比較

表2中列舉了兩軟件多路徑效應、o/slips及采集數據量的比較，各項指標基本相同，說明了RTQC軟件檢測結果的可靠性。

4 結 語

在分析數據解碼和數據質量指標計算方法的基礎上，利用VB.NET面向對象的程序語言設計了RTQC軟件，實現衛星觀測數據質量的實時檢測。結合現場實驗以及TEQC軟件驗證了RTQC軟件的正確性和可靠性。RTQC軟件可以實時對靜態觀測數據進行質量檢測，掌握數據質量的變化以便及時作出各種觀測決策，減少不必要的事后返工重測。例如在CORS基站、工程控制網等的布網選點中，可以利用RTQC軟件對站點觀測環境進行現場檢測，以便選定更優的站點位置，保證觀測數據的質量。

[1] 王世君.運用TEQC軟件進行GPS數據質量分析及其在高速鐵路框架網選址中的應用[J]. 鐵道勘察，2011，37(2)：37-40(Wang Shijun. Applications of TEQC Software in Analysis on Quality of GPS Data and Selection for Location of Frame Network for High Speed Railways[J].Railway Investigation and Surveying,2011,37(2):37-40)

[2] Mitelman A,Phelts R E,Akos D,et al.A Real-time Signal Quality Monitor for GPS Augmentation Systems[C].International Technical Meeting of the Satellite Division of the Institute of Navigation,Salt Lake,2000

[3] Rui D.TEQC Software to Make Reference to the Data Analysis about GNSS Continuously Operating Reference Station’s Choice[J].Urban Geotechnical Investigation & Surveying, 2009,1:24

[4] 張亦梅, 特木其勒, 劉可,等. 應用TEQC對湖北省陸態網絡連續站的觀測數據進行質量分析[J]. 大地測量與地球動力學, 2011, 31(S1):94-97(Zhang Yimei,Temuqile, Liu Ke,et al.Quality Analysis of Observation Data of Continuous Station of CMONOC in Hubei by Using TEQC[J].Journal of Geodesy and Geodynamics, 2011,31(Supp1):94-97)

[5] 張顯云，姚黔貴，張鵬飛，等．基于TEQC的GNSS數據質量可視化表達與分析[J]．測繪通報，2014(4)：26-29(Zhang Xianyun,Yao Qiangui,Zhang Pengfei,et al.The Visual Expression and Analysis of GNSS Data Quality Based on TEQC[J].Bulletin of Surveying and Mapping, 2014(4):26-29)

[6] 劉巖，汪劍云，吳北平，等.基于VB.NET的GPS接收機串口通信的實現[J].地理空間信息，2012，10(2)：47-49(Liu Yan,Wang Jianyun,Wu Beiping,et al.Achievement of Serial Communication of GPS Receiver Based on VB.NET[J].Geospatial Information,2012,10(2): 47-49)

[7] 朱靜然，潘樹國.GNSS原始數據解碼為RINEX標準數據格式及數據質量分析[J].現代測繪，2014，37(4)：8-10(Zhu Jingran,Pan Shuguo.Raw Data of GNSS Decoded to Standard RINEX Format Data and Analysis of Data Quality[J].Modern Surveying and Mapping. 2014,37(4):8-10)

[8] 鄒進貴，李勇.GNSS天線陣列接收機數據采集與解碼的實現及應用[J].測繪通報，2015(8)：1-3(Zou Jingui,Li Yong.Implementation and Application of Data Collecting and Decoding Based on GNSS Antenna Array Receiver[J].Bulletin of Surveying and Mapping, 2015(8):1-3)

[9] 魏二虎，王中平，龔真春，等.TEQC軟件用于GPS控制網數據質量檢測的研究[J].測繪通報，2008(9)：6-9(Wei Erhu,Wang Zhongping,Gong Zhenchun,et al.On the Data Quality Detection of GPS Control Network with TEQC Software[J].Bulletin of Surveying and Mapping,2008(9):6-9)

[10]袁鵬，姜衛平，張彥芬.基于電離層殘差法探測和修復周跳多值性的方法研究[J].測繪信息與工程，2011，36(6)：1-3(Yuan Peng,Jiang Weiping,Zhang Yanfen.Multi-Value of Ionosphere Residual Error Method for Cycle Slips Detection and Reparation[J]. Journal of Geomatics,2011,36(6):1-3)

Real-Time Quality Checking of GPS Static Measurement Data

LUXianjian1,2YANHongbo1,2HUANGYing3

1 College of Geomatics and Geoinformation, Guilin University of Technology, 319 Yanshan Street,Guilin, 541006,China 2 Guangxi Key Laboratory of Spatial Information and Geomatics, 319 Yanshan Street,Guilin,541006,China 3 Guilin Good Surveying Information Technology Co Ltd, 319 Yanshan Street,Guilin,541006,China

Based on analysis of GPS observation data decoding and data quality indices, we compile a new program for decoding the original observation data and calculating and analyzing the quality indices, such as multipath effect, cycle slips and so on. The quality indices can be displayed intuitively by graphics in real-time. Moreover, the experimental verification is carried out through engineering practice, and the results prove the correctness and reliability of the designed software, providing guidance and help for practical work.

static measurement; data quality; real-time checking; data decoding

National Natural Science Foundation of China,No.41161072; Natural Science Foundation of Guangxi,No.2014GXNSFAA118288;Guangxi Key Laboratory Spatial Information and Geomatics Foundation,No.130511409,130511415.

YAN Hongbo,PhD candidate,majors in application of GNSS and remote sensing for soil moisture retrieval, E-mail：2009019@glut.edu.cn.

2015-11-03

項目來源：國家自然科學基金(41161072)；廣西自然科學基金(2014GXNSFAA118288)；廣西空間信息與測繪重點實驗室基金(130511409，130511415)。

盧獻健，講師，主要從事GNSS數據處理及其應用研究，E-mail：285922956@qq.com。

晏紅波，博士生，講師，主要從事GNSS應用與遙感土壤濕度反演研究，E-mail：2009019@glut.edu.cn。

10.14075/j.jgg.2016.11.002

1671-5942(2016)011-0946-05

P228

A

About the first author:LU Xianjian,lecturer,majors in GNSS surveying data processin and its application，E-mail：285922956@qq.com.