Bernese 5.0軟件介紹及常見(jiàn)問(wèn)題分析

蔣欠欠,郭英,戚憲東

(山東科技大學(xué)測(cè)繪 科學(xué)與工程學(xué)院,山東 青島 266590)

?

Bernese 5.0軟件介紹及常見(jiàn)問(wèn)題分析

蔣欠欠,郭英,戚憲東

(山東科技大學(xué)測(cè)繪 科學(xué)與工程學(xué)院,山東 青島 266590)

介紹了Bernese5.0軟件數(shù)據(jù)處理的準(zhǔn)備工作、手工處理內(nèi)容、BPE中控制文件的組成,并通過(guò)自動(dòng)批處理PPP.PCF處理GPS數(shù)據(jù),分析遇到的問(wèn)題及解決的方法。

Bernese5.0;BPE;數(shù)據(jù)處理;問(wèn)題分析

0 引 言

目前在GPS高精度定位中,廣泛使用的數(shù)據(jù)處理軟件有美國(guó)麻省理工學(xué)院和美國(guó)加利福尼亞大學(xué)海洋研究所共同研發(fā)的采用雙差方法進(jìn)行數(shù)據(jù)處理的GAMIT軟件,美國(guó)噴氣動(dòng)力實(shí)驗(yàn)室研發(fā)的以非差方法進(jìn)行數(shù)據(jù)處理的GIPSY軟件[1-2],瑞士伯爾尼大學(xué)研發(fā)的Bernese軟件。升級(jí)后的Bernese5.0軟件既能以非差方式進(jìn)行單點(diǎn)定位,又可以用雙差方法進(jìn)行整網(wǎng)平差[3-4],可以手工處理或批處理(BPE)GPS數(shù)據(jù)和GLONASS數(shù)據(jù),以及處理SLR數(shù)據(jù)、永久GPS跟蹤站觀測(cè)網(wǎng)或混合不同類(lèi)型接收機(jī)的觀測(cè)網(wǎng),還可用于處理精密定軌、估計(jì)地球自轉(zhuǎn)和電離層等多種中間產(chǎn)品數(shù)據(jù)。Bernese軟件具有高精度、高效率的特點(diǎn),適用于科學(xué)研究、實(shí)踐教學(xué),也適用于測(cè)繪單位的數(shù)據(jù)產(chǎn)品計(jì)算,學(xué)習(xí)Bernese軟件具有一定的現(xiàn)實(shí)意義[5]。

但是,利用Bernese軟件進(jìn)行數(shù)據(jù)處理時(shí),各部分的步驟需要緊密銜接,若某一操作步驟出現(xiàn)問(wèn)題,軟件將無(wú)法正常運(yùn)行。因此,解決操作過(guò)程出現(xiàn)的問(wèn)題十分關(guān)鍵。

本文通過(guò)PPP.PCF(Precise Point Positioning.Process Control File)處理2015年第121天LHAZ測(cè)站的GPS數(shù)據(jù),分析并解決在處理過(guò)程遇到的問(wèn)題。

1 Bernese 5.0的結(jié)構(gòu)和主要功能模塊

Bernese5.0軟件根據(jù)數(shù)據(jù)處理功能可分為手工處理和BPE兩大部分[6]。手工處理數(shù)據(jù)和BPE處理數(shù)據(jù)之前,都需要進(jìn)行數(shù)據(jù)文件的準(zhǔn)備,其數(shù)據(jù)處理流程如圖1所示。

圖1 Bernese5.0數(shù)據(jù)處理流程圖

1.1 手工處理部分

手工處理部分主要有格式轉(zhuǎn)換部分、軌道部分、數(shù)據(jù)處理部分、模擬部分和常用工具部分五塊內(nèi)容[7]。

1) 格式轉(zhuǎn)換部分

將RINEX格式的原始觀測(cè)文件、原始導(dǎo)航文件、原始?xì)庀笪募D(zhuǎn)換為二進(jìn)制的Bernese格式。

注意檢查數(shù)據(jù)類(lèi)型、RINEX文件頭信息等,觀測(cè)文件轉(zhuǎn)換為*.PZO、*.PZO、*.CZH和*.CZO四種Bernese格式,在自己新建項(xiàng)目下的RAW中檢查格式轉(zhuǎn)換是否正確完成。

2) 軌道部分

運(yùn)用精密星歷文件、地球自轉(zhuǎn)參數(shù)文件生成軌道列表文件,進(jìn)而生成標(biāo)準(zhǔn)軌道。注意廣播星歷與精密星歷生成標(biāo)準(zhǔn)軌道參數(shù)設(shè)置的區(qū)別,在OUT目錄下檢查衛(wèi)星軌道精度。

3) 數(shù)據(jù)處理部分

這是軟件的重要部分,包括數(shù)據(jù)預(yù)處理(接收機(jī)鐘差改正、生成基線、周跳探測(cè)與修復(fù)等)、基于法方程的參數(shù)估計(jì)。注意各參數(shù)的設(shè)置,由于生成的文件較多,在OUT目錄下檢查各文件中的精度。

4) 模擬部分

根據(jù)統(tǒng)計(jì)數(shù)據(jù)信息(觀測(cè)值的RMS、偏差和周跳)產(chǎn)生模擬的GPS和GLONASS觀測(cè)量(碼或相位)。

5) 常用工具部分

進(jìn)行檢查、編輯轉(zhuǎn)換成Bernese格式的數(shù)據(jù)文件,進(jìn)行坐標(biāo)的顯示、比較,殘差的瀏覽、檢核,輸出文件的處理等。

1.2 批處理部分

BPE則由程序、腳本文件、控制文件等組成,其中控制文件是BPE的核心,主要的控制文件有PPP.PCF、RNX2SNX.PCF、CLKDET.PCF和BASTST.PCF。

1) PPP.PCF

非差數(shù)據(jù)批處理,主要運(yùn)用到的程序模塊有:RNXSMT(數(shù)據(jù)預(yù)處理),RXOBV3(RINEX數(shù)據(jù)格式轉(zhuǎn)換為Bernese格式),POLUPD、PRETAB、ORBGEN(軌道計(jì)算),CODSPP、GPSEST、RESRMS、SATMRK(數(shù)據(jù)預(yù)處理及參數(shù)設(shè)置),GPSEST(基線法方程文件),ADDNEQ2(基于法方程得多時(shí)段解)。

2) RNX2SNX.PCF

求差數(shù)據(jù)批處理,與非差數(shù)據(jù)批處理所用到的程序模塊類(lèi)似,在數(shù)據(jù)預(yù)處理階段用到的程序有:CODSPP、GPSEST、RESRMS、SATMRK、MAUPRP、SNGDIF。從而可求出測(cè)站點(diǎn)位速率、對(duì)流層等參數(shù),解算的結(jié)果以SINEX的格式輸出。

3) CLKDET.PCF

根據(jù)網(wǎng)內(nèi)的觀測(cè)數(shù)據(jù),計(jì)算接收機(jī)鐘差和衛(wèi)星鐘差,生成一個(gè)獨(dú)立的鐘文件。若觀測(cè)文件中存在較大的殘差或缺少觀測(cè)數(shù)據(jù)將不會(huì)被計(jì)算。

4) BASTST.PCF

基線解算,它所用到的數(shù)據(jù)從PPP.PCF、RNX2SNX.PCF,CLKDET.PCF生成的文件中獲取。

1.3 數(shù)據(jù)文件準(zhǔn)備

數(shù)據(jù)文件的準(zhǔn)備工作包括下載RINEX格式的原始觀測(cè)數(shù)據(jù)(*.O文件)存放到ORX和RAW目錄下,將對(duì)應(yīng)的精密星歷文件(*.SP3)、鐘差文件(*.CLK)、地球自轉(zhuǎn)參數(shù)文件(*.ERP)、該月的兩個(gè)碼偏差文件(*.DCB)下載并存放在ORB目錄下,下載電離層文件(*.ION)存放到ATM目錄下,將測(cè)站坐標(biāo)文件(*.CRD)、測(cè)站所在板塊文件(*.PLD)、測(cè)站速度場(chǎng)文件(*.VEL)、海潮文件(*.BLQ)下載并復(fù)制到STA目錄下,鐘差文件存放到OUT目錄中。更新GEN文件,主要更新衛(wèi)星問(wèn)題文件(*.CRX),衛(wèi)星參數(shù)文件(*.ION),地球自轉(zhuǎn)參數(shù)文件(*.ERP),相位中心改正文件等。數(shù)據(jù)文件的下載地址為: ftp://ftp.unibe.ch/aiub/BSWUSER50/ ftp://cddis.gsfc.nasa.gov/gps/.

2 數(shù)據(jù)處理所遇問(wèn)題及解決方法

通過(guò)BPE中PPP.PCF運(yùn)行2015年第121天LHAZ測(cè)站的GPS數(shù)據(jù),其中所遇到的問(wèn)題及解決方法如圖所示。

1) 文件類(lèi)型問(wèn)題,如圖2所示。

圖2 文件類(lèi)型問(wèn)題

問(wèn)題一:是由精密星歷文件的格式?jīng)]有被Bernese識(shí)別引起的。將其*.PRE格式改為Bernese能識(shí)別的*.SP3格式。

問(wèn)題二:是由地球自轉(zhuǎn)參數(shù)文件文件沒(méi)有被Bernese識(shí)別引起的,將IGS18427.ERP的后綴名改為Bernese能運(yùn)行的*.IEP格式。

2) 數(shù)據(jù)格式轉(zhuǎn)換問(wèn)題,如圖3所示。

圖3 數(shù)據(jù)格式問(wèn)題

問(wèn)題一:是由RNXGRA程序中*.STA命名不正確引起的。將選項(xiàng)station information file重新命名為15-121.STA即可。一般情況下,將station information file內(nèi)容清空,默認(rèn)STA目錄下的.STA文件名。

問(wèn)題二:是由STA目錄中命名為15-121.STA的文件與RAW目錄下的lhaz1210.SMT中的測(cè)站名沒(méi)保持一致引起的。兩者測(cè)站名需保持一致且 OLD STATION NAME要正確命名,一般為:測(cè)站名+*。

問(wèn)題三:是由觀測(cè)數(shù)據(jù)與其它的衛(wèi)星系統(tǒng)相拒引起的。在程序SMTBV3中,將SATELLITE SYSTEM SELSETION選項(xiàng)由GPS改為ALL即可。

3) 生成軌道文件時(shí)出現(xiàn)的問(wèn)題,如圖4所示。

圖4 生成軌道文件時(shí)出現(xiàn)的問(wèn)題

問(wèn)題一:是由GEN中缺少2015年的壞衛(wèi)星文件引起的。下載2015年的壞衛(wèi)星文件并拷貝到GEN目錄中,這可以通過(guò)及時(shí)更新GEN文件實(shí)現(xiàn)。

問(wèn)題二:是由于衛(wèi)星文件中有太多衛(wèi)星而衛(wèi)星最大數(shù)量為200引起的。通過(guò)刪除部分SATELLIT.I08中的衛(wèi)星使其衛(wèi)星數(shù)小于200來(lái)解決此問(wèn)題。

4) 數(shù)據(jù)預(yù)處理時(shí)出現(xiàn)的問(wèn)題,如圖5所示。

圖5 數(shù)據(jù)預(yù)處理時(shí)出現(xiàn)的問(wèn)題

此問(wèn)題是由缺少I(mǎi)GS151210.CLK文件引起的。打開(kāi)程序CODSPP接收機(jī)時(shí)鐘同步計(jì)算,在INPUT ORBITS選項(xiàng)中Satellite clocks找不到IGS151210.CLK文件,通過(guò)檢查CODSPP之前的程序,在PRETAB生成軌道列表文件中將RESULT FILES的Satellite clock file正確命名并保存即可生成CODSPP運(yùn)行所需的*.CLK文件。

5) 數(shù)據(jù)不完整問(wèn)題,如圖6所示。

圖6 數(shù)據(jù)不完整問(wèn)題

問(wèn)題一:是由OUT目錄中沒(méi)有IGS18425.CLK文件引起的。下載的IGS18425.CLK鐘差文件既要拷貝到ORB中也要復(fù)制到OUT目錄下,數(shù)據(jù)處理用到OUT目錄下的鐘差文件。

問(wèn)題二:是由天線類(lèi)型找不到引起的。通過(guò)編輯.CRX文件,按照提示添加天線類(lèi)型解決此類(lèi)型問(wèn)題。

3 結(jié)束語(yǔ)

Bernese軟件能高效的處理GPS高精度數(shù)據(jù),尤其使用BPE過(guò)程,軟件運(yùn)行出現(xiàn)錯(cuò)誤中斷時(shí),一般通過(guò)BPE目錄查找出現(xiàn)問(wèn)題的文件,找到并分析錯(cuò)誤,進(jìn)而解決問(wèn)題。本文PPP.PCF處理GPS數(shù)據(jù)時(shí)根據(jù)錯(cuò)誤所在的程序模塊將錯(cuò)誤分為以下5種類(lèi)型:

1) 文件類(lèi)型問(wèn)題,2) 數(shù)據(jù)格式轉(zhuǎn)換問(wèn)題,3) 生成軌道文件時(shí)出現(xiàn)的問(wèn)題,4) 數(shù)據(jù)預(yù)處理時(shí)出現(xiàn)的問(wèn)題,5) 數(shù)據(jù)不完整問(wèn)題，并給出相應(yīng)解決錯(cuò)誤的方法。

[1] 李征航,張小紅.衛(wèi)星導(dǎo)航定位新技術(shù)及高精度數(shù)據(jù)處理方法[M].武昌:武漢大學(xué)出版社,2009:172-174.

[2] CHENG Mengfei，WANG Hua CHENG Yingyan ,etal.The trend of the APRGP velocity Field and plates movement derived from GPS data[J].Springer Journal, 2010,53(4):767-772.

[3] DACH R,HUGENTOBLER U,FRIDEZ P,etal.Bernese GPS Software Version 5.0 [M].Astronomical Institute of Bern,2007.

[4] 曲超.GPS數(shù)據(jù)處理軟件Bernese5.0應(yīng)用[J].中國(guó)科技信息,2015(9):56-58.

[5] 呂浩,張友陽(yáng),呂志平，等.基于Bernese非差的GNSS大網(wǎng)分布式解算[J].測(cè)繪科學(xué)技術(shù)學(xué)報(bào),2015,32(5):69-72.

[6] 孟黎,徐杰,任超.Bernese GPS4.2版本數(shù)據(jù)處理軟件的使用與常見(jiàn)錯(cuò)誤分析[J].全球定位系統(tǒng),2007,33(3):42-46.

[7] 邰賀,付研.應(yīng)用Bernese5.0軟件進(jìn)行精密單點(diǎn)定位計(jì)算[J].中國(guó)學(xué)術(shù)期刊電子出版社,2011(2):74-75.

The Introduction of Bernese 5.0 Software and the Analysis of Common Problems

JIANG Qianqian,GUO Ying,QI Xiandong

(CollegeofGeodesyandGeomatics,ShandongUniversityofScienceandTechnology,Qingdao266590,China)

Introduces the data processing preparation, manual processing content, the control file composition of BPE in Bernese5.0 software and through the PPP process of automatic BPE deal with the GPS data, analyzing problems and Putting forward the method to solve the problem.

Bernese5.0; BPE; data processing; problem analysis

10.13442/j.gnss.1008-9268.2017.02.027

2016-09-19

地理空間信息工程國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)(編號(hào)：201305)； 海島(礁)測(cè)繪技術(shù)國(guó)家測(cè)繪地理局重點(diǎn)實(shí)驗(yàn)室(編號(hào)：2014B01)

P228.4

A

1008-9268(2017)02-0117-04

蔣欠欠 (1992-),女，山東濟(jì)寧人,碩士生,主要從事GNSS定位研究。

郭英 (1973-),女，山東臨沂人,博士、副教授、碩士生導(dǎo)師,主要從事GNSS定位技術(shù)及其完備性、電離層層析理論與應(yīng)用、基于GNSS和多傳感器的室內(nèi)高精度定位研究。

聯(lián)系人: 蔣欠欠E-mail:jiangqian0326@163.com