利用Anubis 檢核BDS 觀測數(shù)據(jù)質(zhì)量

伊 珣，徐愛功，唐龍江

(遼寧工程技術(shù)大學(xué) 測繪與地理科學(xué)學(xué)院，遼寧 阜新 123000)

0 引言

截止到2020 年7 月，隨著北斗衛(wèi)星導(dǎo)航系統(tǒng)（BeiDou navigation satellite system, BDS）最后一顆衛(wèi)星入網(wǎng)工作，代表著我國自主研發(fā)的導(dǎo)航系統(tǒng)全面建成，很快可以對全球范圍內(nèi)的用戶提供定位服務(wù)。而衛(wèi)星觀測數(shù)據(jù)的質(zhì)量好壞，很大程度上影響了定位和其他服務(wù)的精確性，所以對BDS衛(wèi)星的觀測數(shù)據(jù)進(jìn)行質(zhì)量評估是非常有意義的。

由于國外主流的數(shù)據(jù)格式轉(zhuǎn)換、元數(shù)據(jù)編輯及數(shù)據(jù)質(zhì)量檢查（translation, editing and quality checking, TEQC）軟件只能處理接收機(jī)自主交換格式（receiver independent exchange format, RINEX）2.11 版及其之前版本的全球定位系統(tǒng)（global positioning system, GPS）和格洛納斯衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system, GLONASS）數(shù)據(jù)，無法處理BDS 和高版本的RINEX 數(shù)據(jù)，且軟件不開源，可視化效果不夠理想，不能滿足當(dāng)前質(zhì)量分析的需求。德國聯(lián)邦制圖和大地測量局（Das Bundesamt für Kartographie und Geod?sie, BKG）的恩特里普（Ntrip）克連特（Client）軟件，主要應(yīng)用于單點(diǎn)定位，數(shù)據(jù)檢核方面功能不全面。爾特克利布（RTKLIB）處理結(jié)果輸出不夠詳細(xì)，且檢核質(zhì)量和精度不高。格-努特（G-Nut）/阿努比斯（Anubis）軟件相比于TEQC 軟件，同時支持RINEX2 和RINEX3 版本的BDS 數(shù)據(jù)質(zhì)量檢核，精度高且功能強(qiáng)大，所以實驗利用Anubis 軟件來處理北斗衛(wèi)星導(dǎo)航（區(qū)域）系統(tǒng)即北斗二號（BeiDou navigation satellite (regional) system, BDS-2）及第三代北斗衛(wèi)星導(dǎo)航系統(tǒng)即北斗三號（BeiDou navigation satellite system with global coverage, BDS-3）的觀測數(shù)據(jù)質(zhì)量[1-5]，并編程實現(xiàn)數(shù)據(jù)的批量處理與提取。

許多學(xué)者利用Anubis 等軟件對BDS 的數(shù)據(jù)質(zhì)量進(jìn)行過檢核與分析。文獻(xiàn)[6]利用Anubis 軟件，對江門川島區(qū)域三個測站的全球衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system, GNSS）數(shù)據(jù)進(jìn)行了檢核，GPS/GLONASS/BDS 的數(shù)據(jù)可用率均大于90%，伽利略衛(wèi)星導(dǎo)航系統(tǒng)（Galileo navigation satellite system, Galileo）的數(shù)據(jù)可用率為80.02%，GLONASS 的多路徑誤差值最小，各測站的周跳比均小于1；文獻(xiàn)[7]對BDS-2 和BDS-3 中圓地球軌道（medium Earth orbit, MEO）衛(wèi)星的B3I 頻率數(shù)據(jù)質(zhì)量進(jìn)行了分析，發(fā)現(xiàn)BDS-3 衛(wèi)星數(shù)據(jù)完整率和信噪比高于BDS-2，多路徑誤差小于BDS-2，沒有在BDS-3 中觀測到BDS-2 存在的系統(tǒng)誤差[7]；文獻(xiàn)[8]采用Anubis 軟件和TEQC 軟件對GPS、GLONASS、BDS 和Galileo 觀測數(shù)據(jù)的進(jìn)行質(zhì)量檢查分析，結(jié)果表明這兩款軟件檢核結(jié)果相當(dāng)，數(shù)據(jù)質(zhì)量均符合國際全球衛(wèi)星導(dǎo)航定位系統(tǒng)服務(wù)組織（International global navigation satellite systems（GNSS） Service, IGS）的經(jīng)驗值；文獻(xiàn)[9]利用Anubis 軟件對GPS/BDS 觀測數(shù)據(jù)的質(zhì)量進(jìn)行了分析對比，結(jié)果表明，BDS/GPS 觀測數(shù)據(jù)的質(zhì)量基本處于同一水平。本文將從數(shù)據(jù)完整率、信噪比、周跳和多路徑誤差4 個方面，對BDS-2/BDS-3 B1I/B3I服務(wù)信號的觀測數(shù)據(jù)質(zhì)量進(jìn)行評估。

1 北斗衛(wèi)星導(dǎo)航系統(tǒng)簡介

2020 年7 月，第55 顆BDS 衛(wèi)星成功入網(wǎng)，至此BDS 全面建成。BDS-2 的第1-16 號衛(wèi)星提供B1I/B2I/B3I 三個公開服務(wù)信號；BDS-3 第19-30號衛(wèi)星和第 32-46 號衛(wèi)星提供 B1I/B3I/B1C/B2a/B2b 信號；第59-60 號衛(wèi)星提供B1I/B3I 信號；第31 號衛(wèi)星、第56-58 號衛(wèi)星和第61 號衛(wèi)星是BDS-3 的實驗衛(wèi)星，信號體制與BDS-3 不完全相同[10]。本文對IGS 提供的2020-02-29 的觀測文件進(jìn)行統(tǒng)計。目前全球分別有283 個和285 個測站能夠觀測到BDS-2 和BDS-3 衛(wèi)星，觀測站的數(shù)量呈遞增趨勢，越來越多的國家和地區(qū)重視對BDS衛(wèi)星的觀測。圖1 為BDS 衛(wèi)星觀測站數(shù)量逐年趨勢圖。

圖1 BDS 衛(wèi)星觀測站數(shù)量逐年趨勢

2 Anubis 軟件功能與拓展

2.1 Anubis 軟件簡介

Anubis 軟件是由捷克國家大地測量、地形與地圖制圖研究所開發(fā)的一款可以對所有衛(wèi)星導(dǎo)航系統(tǒng)觀測數(shù)據(jù)都適用的免費(fèi)開源質(zhì)量分析軟件，由C++語言編寫，可以通過命令行窗口進(jìn)行多系統(tǒng)的文件處理，內(nèi)容強(qiáng)大，處理數(shù)據(jù)精準(zhǔn)。用戶可以根據(jù)自己需求調(diào)整配置文件，來獲取想要信息的詳細(xì)程度。檢核質(zhì)量精準(zhǔn)，容易實現(xiàn)可視化，支持視窗（Windows）和利努克斯（Linux）系統(tǒng)，支持 RINEX2.xx 和 RINEX3.0x 格式。Anubis 軟件檢核的結(jié)果文件為.xtr 后綴文件，信噪比、周跳和多路徑誤差均可以精確到每顆衛(wèi)星的每個信號。

2.2 Anubis 軟件數(shù)據(jù)批量處理與提取

Anubis 軟件每次只能處理一天的數(shù)據(jù)，而且每次處理數(shù)據(jù)前都需要手動更改配置文件，修改文件名和觀測時段等各種信息。如果處理數(shù)據(jù)的時間跨度大、測站數(shù)量較多時，效率就很低了。本文利用派森（Python）語言對Anubis 軟件的功能進(jìn)行了拓展，輸入處理數(shù)據(jù)的起始年積日和結(jié)束年積日后，程序就可以自動更改配置文件。由于本文實驗分年份進(jìn)行，所以設(shè)置其可對一年內(nèi)的多個測站的數(shù)據(jù)質(zhì)量檢查，自動提取所需要的質(zhì)量評估結(jié)果，包括各個信號和各衛(wèi)星的多路徑誤差、周跳、信噪比和數(shù)據(jù)完整率，程序主函數(shù)如圖2，提取結(jié)果如表1 所示。

圖2 拓展功能主函數(shù)

表1 提取的結(jié)果文件

程序的處理流程為：1）從布勒德姆（BRDM）文件中獲取文件的年月日，并在配置文件中修改為對應(yīng)日期；2）如果缺少文件則會報錯，并跳過出錯的年積日，繼續(xù)進(jìn)行檢核；3）獲取文件夾下的o 文件和BRDM 文件名，并在配置文件中修改；4）執(zhí)行Anubis 軟件的檢核命令A(yù)nubis -x anub-2.1.2.xml -l process.log；5）讀取結(jié)果文件，將其中需要的部分提取出來，并分類存到指定文件夾的文件中；6）重復(fù)以上步驟直至處理完所有數(shù)據(jù)。

3 BDS 數(shù)據(jù)質(zhì)量檢核

3.1 信噪比檢核

圖3 RINEX 文件中的信噪比

3.2 多路徑誤差檢核

在GNSS 觀測中，由于接收機(jī)可以接收到周圍環(huán)境反射的衛(wèi)星信號，對直接接收的衛(wèi)星信號產(chǎn)生干擾，導(dǎo)致觀測結(jié)果與真值產(chǎn)生偏差，就是多路徑誤差[13]。多路徑誤差通常反映了周圍環(huán)境和其他因素對信號傳播的影響，多路徑誤差越小，說明抗多路徑能力越強(qiáng)。多路徑誤差通?？梢酝ㄟ^偽距和載波相位觀測值組合求得[14]，即

式中：Mk為偽距多路徑誤差；kP為雙頻偽距觀測值；iL和Lj為雙頻載波相位觀測值；fi、fj和fk為載波頻率；k、i和j為頻率編號。

3.3 數(shù)據(jù)完整性檢核

受到測站周圍環(huán)境的高大建筑物或植物等影響，以及在信號傳輸過程中造成的數(shù)據(jù)缺失，數(shù)據(jù)解析過程中校驗不通過等原因，在用戶接收到衛(wèi)星的觀測數(shù)據(jù)時數(shù)據(jù)不完整，導(dǎo)致定位結(jié)果偏差，嚴(yán)重時甚至無法解算定位。數(shù)據(jù)完整性是可用歷元和理論歷元的比值，反映了數(shù)據(jù)的連續(xù)情況和質(zhì)量水平，數(shù)據(jù)完整性的比值越大，數(shù)據(jù)質(zhì)量就越好。一般來說，數(shù)據(jù)完整性小于85%，就應(yīng)該考慮優(yōu)化測站的位置或者周圍環(huán)境，數(shù)據(jù)完整性定義為

式中：C、T分別為可用歷元和理論歷元，理論歷元就是自觀測時間起，根據(jù)采樣間隔和結(jié)束時間來計算出來的，對評價測站指標(biāo)有重要意義。

3.4 周跳比檢核

周跳是指GNSS 接收機(jī)在跟蹤載波信號時發(fā)生的暫時失鎖現(xiàn)象，導(dǎo)致整周計數(shù)發(fā)生跳變或者中斷。周跳是影響測量精度的一個重要指標(biāo)，當(dāng)數(shù)據(jù)發(fā)生異常、變化超過某一個閾值，就認(rèn)為發(fā)生了周跳，有的學(xué)者利用周跳比來反映周跳的情況，周跳比的值越小表示周跳越少，數(shù)據(jù)質(zhì)量越好，計算公式為

式中：S、O分別為周跳數(shù)與總數(shù)據(jù)個數(shù)；J為周跳比。

4 BDS 數(shù)據(jù)質(zhì)量檢核實驗分析

本文利用2017 年、2018 年和2019 年共3 年，全球26 個測站的觀測數(shù)據(jù)對BDS 數(shù)據(jù)質(zhì)量進(jìn)行檢核。采集數(shù)據(jù)時，衛(wèi)星截止高度角均設(shè)置為15°，采樣間隔為30 s，采用BRDM 混合廣播星歷。利用Anubis 的質(zhì)量檢核結(jié)果來評估分析BDS-2 和BDS-3的數(shù)據(jù)質(zhì)量，表2 為所選測站分布地區(qū)，圖4 為各測站每年的觀測天數(shù)和測站觀測BDS 衛(wèi)星數(shù)。

In regards to the short term clinical outcomes, we studied the procedure time, the time to resume diet,the time to full ambulation, the duration of the total hospital stay and the complication rate.

圖4 測站觀測數(shù)據(jù)情況

表2 測站分布

由圖4 可知，KIR8 站和MAR7 站的年觀測天數(shù)均少于100 d，ARUC 站、BSHM 站、FFMJ 站和HUEG 站的2017 年觀測天數(shù)少于250 d，其余的測站每年的平均觀測天數(shù)均在250 d 以上，觀測時間滿足實驗要求。從圖4（b）可以看出，測站觀測的BDS 衛(wèi)星數(shù)量呈現(xiàn)逐年增加的趨勢，在2019 年呈現(xiàn)爆發(fā)式增長，除TIT2 站僅觀測16 顆以外，其余測站觀測數(shù)量均在20 顆以上，說明最近幾年BDS不僅發(fā)展迅速，而且越來越被各個國家和地區(qū)重視。

圖5 為各測站3 年數(shù)據(jù)完整率的平均值。

圖5 B1I、B3I 信號的數(shù)據(jù)完整率

從圖5 可以看出：部分測站沒有觀測B3I 信號，B1I 和B3I 信號在不同測站的數(shù)據(jù)完整率有所差異，BSHM 站、KIR8 站、LPGS 站和MET3 站B3I 信號數(shù)據(jù)完整率優(yōu)于B1I 信號，其余測站B1I信號均優(yōu)于B3I 信號；數(shù)據(jù)完整性最大為ARUC站的B1I 信號，完整率為97.39%，最小的是GAMG站的B3I 信號，僅為56.91%。除ARUC 站、GAMG站、HUEG 站、KIR8 站、TIT2 站和YEL2 站外，其余測站的兩種信號數(shù)據(jù)完整率均大于80%。

圖6 為各測站3 年的信噪比平均值。從圖6 可以看出：LPGS 站BDS-2 B3I 信號的信噪比大于BDS-3 B3I 信號，MET3 站BDS-2 B1I 信號的信噪比大于BDS-3B1I 信號；除LPGS 站及MET3 站之外，其他測站BDS-3 B1I 信號和BDS-3 B3I 信號的信噪比均優(yōu)于BDS-2 衛(wèi)星。ARUC 站BDS-2 B3I信號的信噪比最低，為32.97 dB·Hz，HOB2 站BDS-3 B3I 的信噪比最大，為47.59 dB·Hz；除ARUC 站、BSHM 站和MET3 站外，其他測站BDS-3 B1I 信號和BDS-3 B3I 信號的信噪比均大于40 dB·Hz，由此可以看出，BDS-3 信號的信噪比整體優(yōu)于BDS-2信號的信噪比。

圖6 各測站信噪比

圖7 為各測站3 年的多路徑誤差平均值。

圖7 多路徑誤差

從圖7 可以看出：BDS-3 B1I 信號的質(zhì)量相較于BDS-2 B1I 信號的質(zhì)量有了很大的提升，只有CEDU 站和KRGG 站BDS-3 B1I 信號的多路徑誤差大于BDS-2 B1I 信號；B3I 信號的整體質(zhì)量優(yōu)于B1I 信號，說明測站周圍環(huán)境對B1I 信號的影響較大。BDS-2 B1I 信號的多路徑誤差很差，最多達(dá)到0.71 m。BDS-3 B3I 信號的整體質(zhì)量來說很好，除了ZIM3 站外，其余均在0.3 m 以內(nèi)。

由于BSHM 站的B1I 信號和B3I 信號的周跳異常，都在1 000 以上，所以將其剔除后重新進(jìn)行繪圖，如圖8 所示。

從圖8 可以看出：在CEDU 站、GAMG 站、HARB 站、HOB2 站、PADO 站、RGDG 站、UNB3站和YEL2 站上，B1I 信號及B3I 信號的周跳比均小于10，信號質(zhì)量很穩(wěn)定，其余跟蹤站均有大量周跳產(chǎn)生，其中DYNG 站及HUEG 站的B1I 信號及B3I 信號，F(xiàn)FMJ 站及WTZZ 站的B1I 信號的周跳比超過了100，大約發(fā)生了400 次周跳，穩(wěn)定性很差。

圖8 各測站周跳比

選取上述實驗中結(jié)果較好的 CEDU 站和HOB2 站，檢核2017—2019 年的觀測數(shù)據(jù)，對BDS-2和BDS-3 衛(wèi)星的數(shù)據(jù)質(zhì)量進(jìn)行進(jìn)一步分析。因所選數(shù)據(jù)僅能觀測到BDS-3 的MEO 衛(wèi)星，下文僅對BDS-2/BDS-3 的MEO 衛(wèi)星進(jìn)行區(qū)分，地球靜止軌道（geostationary Earth orbit,GEO）衛(wèi)星和傾斜地球同步軌道（inclined geosynchronous orbits, IGSO）衛(wèi)星均為BDS-2 衛(wèi)星。

圖9 為CEDU 站和HOB2 站3 年的數(shù)據(jù)完整率。

圖9 數(shù)據(jù)完整率

由圖9 可知:CEDU 站的B1I 和B3I 在2017年1—3 月數(shù)據(jù)完整性很差，最低只有10%，同時期的HOB2 站沒有觀測數(shù)據(jù)，此時間BDS 剛剛發(fā)射了第22 顆和第23 顆衛(wèi)星，還處于發(fā)展階段，尚未在全球范圍內(nèi)提供高精度定位服務(wù)，而CEDU 站和HOB2 站均處于大洋洲，導(dǎo)致數(shù)據(jù)完整率差。在此之后兩個測站的B1I 信號完整率均在80%以上，B3I 信號在2019 年較差，完整率在60%～80%間浮動。

圖10 和圖11 為CEDU 站和HOB2 站的信噪比和多路徑誤差。

圖10 信噪比

由圖10 可以看出所有統(tǒng)計的衛(wèi)星信噪比均在40 dB·Hz 以上，3 年期間基本保持穩(wěn)定，波動不大，BDS-2 的MEO 衛(wèi)星和BDS-3 的MEO 衛(wèi)星信噪比相當(dāng)，且皆大于GEO 和IGSO 衛(wèi)星。由圖11可知，除了BDS-2 的MEO 衛(wèi)星以外，其余衛(wèi)星多路徑誤差均在0.2 m 以內(nèi)。除2017 年CEDU 站的GEO 衛(wèi)星之外，GEO 和IGSO 的多路徑誤差均優(yōu)于MEO 衛(wèi)星，且BDS-3 優(yōu)于BDS-2 的GEO 衛(wèi)星0.05 m 左右。

圖11 多路徑誤差

結(jié)合實驗結(jié)果，BDS-3 的信噪比和多路徑誤差優(yōu)于BDS-2。BDS-3 與BDS-2 的MEO 衛(wèi)星信噪比差別不大，但多路徑誤差減少了0.05 m。B1I 信號的數(shù)據(jù)完整率整體優(yōu)于B3I 信號，但B3I 信號的多路徑誤差小于B1I 信號。

5 結(jié)束語

本文利用Anubis 軟件對26 個測站2017—2019 共3 年的數(shù)據(jù)進(jìn)行檢核，并對其數(shù)據(jù)完整率、信噪比、多路徑誤差和周跳4 個檢核標(biāo)準(zhǔn)進(jìn)行分析，結(jié)論如下：

1）在同時觀測到兩種信號的測站中，有12 個測站B1I 信號完整率大于B3I 信號，僅有4 個測站B3I 信號完整率大于B1I 信號。其余測站的兩種信號數(shù)據(jù)完整率均大于 80%。除 ARUC 站、GAMG 站、HUEG 站、KIR8 站、TIT2 站和YEL2站外，其余測站的兩種信號數(shù)據(jù)完整率均大于80%。CEDU 站和HOB 站在2017 年4 月前的完整率很差，低于60%，隨著BDS 的發(fā)展，在2017 年中旬之后有大幅提升。

2）除ARUC 站、BSHM 站和MET3 站外，其余測站的BDS-3 兩種信號信噪比均大于40 dB·Hz，總體質(zhì)量能滿足定位需求。僅有LPGS 站和MET3站的BDS-2 信噪比高于BDS-3，可以看出BDS-3的B1I 信號和B3I 信號的信噪比大體上優(yōu)于BDS-2衛(wèi)星。BDS-2 和BDS-3 的MEO 衛(wèi)星信噪比差別不大，皆大于GEO 和IGSO 衛(wèi)星，在選取的時間段內(nèi)，信噪比均大于40 dB·Hz，非常穩(wěn)定。

3）多路徑誤差方面，BDS-3 相較于BDS-2 提升了很多，B3I 信號整體質(zhì)量優(yōu)于B1I 信號，說明測站周圍環(huán)境對B1I 信號的影響較大。大部分測站均在0.3 m 以內(nèi)，質(zhì)量滿足要求。GEO 和IGSO的多路徑誤差均優(yōu)于MEO 衛(wèi)星，BDS-3 的MEO衛(wèi)星優(yōu)于BDS-2 的同類衛(wèi)星0.05 m 左右。

4）有16 個測站的周跳比大于10，并有4 個測站的周跳比超過100，說明BDS 衛(wèi)星的數(shù)據(jù)質(zhì)量不夠穩(wěn)定，B1I 信號的浮動最大。