Bernese數(shù)據(jù)處理流程

Bernese處理數(shù)據(jù)前的準(zhǔn)備工作（一）建立項(xiàng)目和定義時(shí)段在Bernese處理數(shù)據(jù)之前，首先要定義一個(gè)項(xiàng)目并激活為它當(dāng)前使用項(xiàng)目，生成項(xiàng)目下的子目錄，然后將數(shù)據(jù)文件拷貝進(jìn)這些子目錄中，還要收集并指定一些與項(xiàng)目有關(guān)的其它基礎(chǔ)信息。從Bernese的5.0版本開始，用戶還不得不為即將處理的數(shù)據(jù)選擇處理時(shí)段。先來定義一個(gè)項(xiàng)目。首先選擇"Campaign>Editlistofcampaigns",出現(xiàn)如下圖所示的界面。這里要求輸入項(xiàng)目存在的路徑和名稱。值得一提的是，上圖表示了兩種路徑的表示方法。第一個(gè)項(xiàng)目Intro很顯然是即將生成在D盤的根目錄下，而第二個(gè)項(xiàng)目Intro_01的路徑則為軟件規(guī)定項(xiàng)目存在的默認(rèn)路徑。也就是說，${P}符號的意義就是將項(xiàng)目生成在軟件規(guī)定的默認(rèn)目錄下，至于默認(rèn)目錄在哪里，下面將給出解答。點(diǎn)擊ASave進(jìn)行保存。然后選擇"Campaign>Selectactivecampaign",選擇一個(gè)當(dāng)前激活的項(xiàng)目。

這里選擇默認(rèn)目錄下的Intro_01項(xiàng)目為激活狀態(tài)，點(diǎn)擊OK接下來就是在剛才激活的項(xiàng)目下生成子目錄。選擇"Campaign>Createnewcampaign",出現(xiàn)下圖所示。這里左側(cè)的一列即為項(xiàng)目下的子目錄。點(diǎn)擊ARun,然后找到Bernese軟件安裝目錄下的"...\Bernese\GPSDATA"文件夾，就可以看到剛才新建的這個(gè)Intro_01項(xiàng)目，并且在項(xiàng)目文件夾內(nèi)同時(shí)生成了和上圖左側(cè)一列相同的共9個(gè)同名文件夾。從這里可以看出，新建項(xiàng)目的默認(rèn)路徑就在"...\Bernese\GPSDATA"。至此，一個(gè)項(xiàng)目的新建過程完成。其中，這9個(gè)文件夾存放的內(nèi)容分別為：ATM存放項(xiàng)目相關(guān)的大氣層文彳（例如電離層文件ION,電離層文件TRPBPE:存放BPE批處理時(shí)每一步驟完成的狀態(tài)信息文件OBS存放Bernese的觀測文件ORB存放跟軌道相關(guān)的文件（軌道文件、地球自轉(zhuǎn)參數(shù)、衛(wèi)星鐘差文件等）ORX存放原始RINEX文件OUT存放輸出文件RAW存放可用于計(jì)算的RINEX文件SOL:存放結(jié)果文件（例如法方程文件，SINEX文件）STA:存放項(xiàng)目相關(guān)的坐標(biāo)和坐標(biāo)信息文件等，項(xiàng)目時(shí)段信息表也在這里接下來設(shè)置時(shí)段。時(shí)段是一段涵蓋了所有即將處理的數(shù)據(jù)所處的時(shí)間段。也就是說,接下來所有處理的數(shù)據(jù)都要在這里所設(shè)置的時(shí)段以內(nèi)。首先選擇"Campaign>Editsessiontable"將會(huì)出現(xiàn)時(shí)段表的編輯界面，如下所示。需要說明的是，時(shí)段表的編輯界面默認(rèn)是只會(huì)顯示上圖的第一行的。IDENTIFIER這列表示的是時(shí)段的四位標(biāo)識符，前三位是年積日，后一位如果為0,則表示時(shí)段為該天全天，時(shí)段表的默認(rèn)狀態(tài)就是0。這一行下面的由A~X共24個(gè)大寫字母則是我手工添加進(jìn)去的，每個(gè)代表一個(gè)小時(shí)。通過后兩列的起始和結(jié)束歷元可以自己指定特定的時(shí)段。這里對時(shí)段表的編輯，是為了下面的步驟選擇時(shí)段做的準(zhǔn)備工作。最后點(diǎn)擊eave保存修改。選擇"Configure>Setsession/computedate",出現(xiàn)如下圖選擇時(shí)段的窗口。

Bernese處理數(shù)據(jù)前的準(zhǔn)備工作（二）創(chuàng)建測站文件上一篇介紹了Bernese處理數(shù)據(jù)前的準(zhǔn)備工作之建立項(xiàng)目和定義時(shí)段，接下來要做的就是為創(chuàng)建測站文件了。對一些站文件，Bernese軟件的"Campaign>Editstationfiles"菜單提供掩碼形式的編輯。因?yàn)樗械恼疚募际且訟SCII形式存在的，所以用戶可以用各種工具瀏覽并修改。創(chuàng)建測站文件的過程主要分為三個(gè)步驟：創(chuàng)建參考系/參考速度文件，創(chuàng)建測站信息文件和其他測站文件。首先來說明如何創(chuàng)建參考系和參考速度文件。在將要處理的項(xiàng)目中，如果任何站所處的參考系在ITRF中的精度都沒有到達(dá)分米級別，那么建議在處理過程中至少加入一個(gè)項(xiàng)目所處地域附近的IGS跟蹤站一并參與處理。也就是說，對于參與處理的項(xiàng)目，都必須建立相應(yīng)的參考系和參考速度，對此用戶可以從ITRF或者IGS參考框架選擇與項(xiàng)目對應(yīng)的SINEX文件。Bernese軟件中的SNX2Q0("Conversion>SINEXtonormalequations")這個(gè)程序就是獲取SINEX文件中的參考系和參考速度的，在使用之前必須先將SINEX文件拷貝到項(xiàng)目子目錄SOL中去。需要說明的是，ITRF參考框架最新為ITRF2000,用戶可以到aHYPERLINKftp://lareg.ensg.ign.fr/pub/itrf/itrf2000”下載相應(yīng)的.SNX格式文件。IGS參考框架最新為IGS00b,用戶可以到“http://www.aiub.unibe.ch/download/BSWUSER50/STA”，選擇下載IGS_00_R.CRD^IGS_00_R.VEL文件。這兩個(gè)參考系和參考速度文件已經(jīng)是Bernese標(biāo)準(zhǔn)格式了，可以直接拷貝到項(xiàng)目子目錄STA路徑下。用"Campaign>Editstationfiles>Stationcoordinates/Stationvelocities"可以手動(dòng)編輯修改參考系和參考速度文件。官方建議這些文件僅含有先驗(yàn)參考系和參考速度達(dá)到1cm以上精度的測站，剩下的其他測站的信息則會(huì)在以后的PPP處理過程中再添加進(jìn)去。【實(shí)踐】我在按上述方法下載了IGS_00_R.CRD和IGS_00_R.VEL這兩個(gè)文件并拷貝到“…\Bernese\GPSDATA\Intro_01\STA”路徑后，分別選擇"Campaign>Editstationfiles>Stationcoordinates/Stationvelocities",按提示選擇這兩個(gè)文件，就會(huì)出現(xiàn)如下圖所示的手動(dòng)編輯參考系和參考速度文件的界面，這也就驗(yàn)證了上述操作步驟的正確性。GO■岬WagiMi”■白呻fiWdKX審蟒鼻鹿3rg1翻NiC^MM^r娶愁/第兇呻EDIT皿8工7￥FILEriieit￡dited；/Ftojiu?￡its/BtEtieEe/gpseata/1：itr*_a1it5S_oatvelTitie：|LiS?^!Tr^<'JRJ5I.SS^T7^^?0<^Fi7NsStat!onJi&rnvY(m/y)rE-事IfrlAta1,AJAC100T7B.0BE-o.dmgI-0.0053-1ig-l+l?-1￡二上汕口；寸-G-.CLCfl-Q.QPL4F.WMi.±1|JA.,.'iT5Tj50^021^1?i'.工二￡N-ZJIMT刨401ij4MUij2一電-0^.(5054J件Emic50137W041%.g、,；.I*J5T二|±l-0-.ODF4tun比d強(qiáng)二±|-ff.0228'0.0(?!Tji-1備1卜AMC24041250^4'--0-C1U--0-0110ifIJL!rJL'***I.J-Pl18*產(chǎn)?樹勺3便*151軾-，上；1-：騰用口*；0；WQgM刊FftHCDfTVUL”.接下來是創(chuàng)建測站信息文件。有兩點(diǎn)非常重要，一是確保處理過程中測站的名字正確；二是確保測站上接收機(jī)和天線型號以及天線高的正確。登陸"ftp://ftp.unibe.ch/aiub/BSWUSER50/STA/",可以下載IGS.STA和EUREF.ST腦兩個(gè)文件。它們包含了IGS和EUREF艮蹤站的測站信息文件，分別由IGS和EUREF中心維護(hù)，由igs.snx和euref.snx文件轉(zhuǎn)換而來。用戶也可以自己下載.snx文件(例如"HYPERLINK/pub/station/general/igs.snx"),然后將文件放到項(xiàng)目的SOL目錄下，通過Bernese軟件的SNX2ST理序("Conversion>SINEXtostationinformation")來手動(dòng)完成.snx文件到.STA文件的轉(zhuǎn)換。同樣的，程序RNX2STA"RINEX>RINEXutilities>Extractstationinformation")可以用來提取已經(jīng)存放在項(xiàng)目RAW!錄下的一系列RINEX的頭文件信息，從而獲取測站的信息文件。官方Manual建議用戶用觀測協(xié)議或者任何其他測站設(shè)置的描述來驗(yàn)證測站信息。無論如何，在項(xiàng)目開始處理之前，要檢查和修改好測站信息文件，可以通過"Campaign>Editstationfiles>Stationinformation"來實(shí)現(xiàn)。測站信息文件中的內(nèi)容定義了測站的名稱(在所有測站文件中需要保持統(tǒng)一)，觀測儀器和天線偏心率(一定要在"Phasecentereccentricity"文件中定義)等重要信息，因此一定要準(zhǔn)備充分和正確。【實(shí)踐】這一部分我的實(shí)踐如下。我先下載了IGS.STA和EUREF.ST腦兩個(gè)文件，放在了Intro_01項(xiàng)目的STA目錄下，然后選擇"Campaign>Editstationfiles>Stationinformation",在彈出的對話框中分別選擇IGS.STA和EUREF.STA出現(xiàn)的界面如下。

這兩張圖分別顯示的是IGS和EUREF艮蹤站的相關(guān)信息。然后，我到IGS的數(shù)據(jù)庫下載了2008年9月25日，跟蹤站wuhn的30s一天的數(shù)據(jù)文件wuhn2690.08o,將其放在RAV\M錄下，然后選擇"RINEX>RINEXutilities>Extractstationinformation",出現(xiàn)如下圖所示的界面。

我在INPUTFILES一欄里面選擇的是原始RINEX觀測文件，并鍵入文件名，在RESULFILES中為即將提取生成的.STA文件命名，在這里我同樣命名為wuhan2690.STAo點(diǎn)擊*un按鈕后，會(huì)出現(xiàn)提示對話框，忽略它，繼續(xù)點(diǎn)擊主窗口右下角的^Output按鈕，會(huì)出現(xiàn)從RN文件中提取出的STA文件的相關(guān)信息，這表明RNX2ST網(wǎng)序運(yùn)行成功，如下圖所示。E七tqfg;sta曰也工也口3：葉CifSBafEir加工口Y看r盛:0兄5.GPun■曰st1s白壬工stsTisnrnt"r?a*i*nfrmarinfk711?iD"設(shè)等Lg：S<p^\lnx；r0-.01Xt?□itsession;力圖另Q”ar￡003Mt4；17—Mp-2V13￡Q;14'oae:L^JiangF1KEMfilsfor-Mtra匚tingth?jrtttioninformationFil4RINEXfiles1SJP}Mnz；roOl^RA^wy^nS^O,i?0OF3_l?舊t—t：110nlr.ajnef&B-cBi.vwc^wu-A^t?inman4afii笠目iqRtCONUHbtASttXECHa-ltII3A5HTD0936E2.361一完成之后在STA文件夾內(nèi)也會(huì)生成相應(yīng)的wuhn0269.STA這個(gè)文件。有的項(xiàng)目處理中，還需要一些其他的測站文件。創(chuàng)建FIX文件選擇處理過程中的參考站。在"Campaign>Editstationfiles>Stationselectionlists"中，如果選擇之前在STA目錄下的參考系文件IGS_00.CRD也就意味著列表中包含了IGS的參考站。創(chuàng)建PLD文件——NUVE限型中的構(gòu)造板塊分配。在"Campaign>Editstationfiles>Tectonicplateassignment"中可以進(jìn)行操作。禾ij用"Campaign>Editstationfiles>Abbrevationtable"可以來自定義處理過程中的測站名成為表單。有時(shí)還需要建立海洋潮汐負(fù)荷等文件，這里不做詳細(xì)討論。至此，關(guān)于Bernese開始數(shù)據(jù)處理前的文件準(zhǔn)備工作完畢。Bernese數(shù)據(jù)處理實(shí)例本文將對Bernese進(jìn)行非差和雙差處理流程進(jìn)行介紹和總結(jié)。一■、基本概述操作系統(tǒng)：Windows8軟件版本：BerneseGPSSoftware5.0中國境內(nèi)選取10個(gè)IGS跟蹤站：bjfs,chan,guao,kunm,lhaz,shao,wuhn,xian,tnml,twtf數(shù)據(jù)日期：2008年8月8日，對應(yīng)的年積日為221,GPS周為1491+5。二、數(shù)據(jù)準(zhǔn)備至UHYPERLINK/pub/gps/data/daily/2008/221/08o下載以上十個(gè)跟蹤站的RINEX格式原始數(shù)據(jù)文件：bjfs2210.08o,chan2210.08o,guao2210.08o,kunm2210.08o,lhaz2210.08o,shao2210.08o,wuhn2210.08o,xian2210.08o,tnml2210.08o,twtf2210.08o。至U/pub/gps/products/1491/下載對應(yīng)當(dāng)天的igs發(fā)布的精密星歷文件、精密衛(wèi)星鐘差文件，以及對應(yīng)該周的地球自轉(zhuǎn)參數(shù)文件：igs14915.sp3、igs14915.clk、igs14917.erp（注：IGS格式的地球自轉(zhuǎn)參數(shù)文件是以周為單位發(fā)布的，因此文件名的最后一位為7）。之后修改igs14917.erp文件的后綴為Bernese軟件能識別的igs14917.iep格式。到ftp://ftp.unibe.ch/aiub/CODE/2008/下載對應(yīng)當(dāng)天的電離層文件、對應(yīng)該月的兩個(gè)DC/：件：COD14915.IOZP1C10808.DCB和P1P20808.DCR三、建立項(xiàng)目在軟件界面以此點(diǎn)選Menu>Campaign>Editlistofcampaigns,添加項(xiàng)目Olympic作為此次處理的項(xiàng)目名稱，之后點(diǎn)擊Save保存。點(diǎn)選Menu>Campaign>Selectactivecampaign,在彈出的窗口中選擇剛才新建的Olympic項(xiàng)目為當(dāng)前項(xiàng)目。

此時(shí)點(diǎn)OK之后會(huì)提示缺少處理所需的時(shí)段文件，先忽略它，因?yàn)轳R上設(shè)置時(shí)段就行了。點(diǎn)選Menu>Con巾gure>Setsession/computedate,彈出設(shè)置時(shí)段窗口，在這里將處理時(shí)段設(shè)定為2008年8月8日，點(diǎn)擊Compute,會(huì)將日期自動(dòng)計(jì)算為其他幾種常見格式。點(diǎn)擊OK此時(shí)會(huì)提示找不到SESSIONS件，忽略它，馬上生成項(xiàng)目子目錄后就會(huì)生成默認(rèn)的SESSION^件了。點(diǎn)選Menu>Campaign>Createnewcampaign,點(diǎn)擊Run,就會(huì)在Bernese軟件目錄下的GPSDAT&：件下生成一個(gè)名為Olympic的項(xiàng)目文件夾，文件夾內(nèi)會(huì)有一些子目錄，這些子目錄目前都是空的，除了之前提到過在STA文件夾下自動(dòng)生成的SESSIONS件。

i^^TlNEMiMO唱門i^^TlNEMiMO唱門PiE-EMtFiJKliJi期心業(yè)mrt&riGrxl?s5eurwdwuwrya*dil『F'F--■/-Xs.--.d±J.!￡"*，：,m三日二1土J*jI±j*jI±j：，:.：ZJ±]，工d±l.IhMT1MKwa1KRun^i、□=.一.==*a一i=.■■ifcAamn^-a■_■■■--_-_-：■：二1l.f??c，g*」「域的f.Y_QS*qa相TF＞，“1EH：g1，”5TF--至此，項(xiàng)目創(chuàng)建完成。、數(shù)據(jù)拷貝至項(xiàng)目下Olympic目錄的將之前下載的10個(gè)RINEX格式的IGS跟蹤站全天觀測文件拷貝到和RAV\M錄下，將精密星歷文件igs14915.sp3、精密衛(wèi)星鐘差文件igs14915.clk、地球自轉(zhuǎn)參數(shù)文件igs14917.iep,以及兩個(gè)DCB^d41P1C10808.DCBP1P20808.DC瞋五個(gè)文件一起放入ORBOlympic目錄的四、創(chuàng)建和修改測站信息文件點(diǎn)選Menu>RINEX>RINEXutilities>Extractstationinformation,在彈出的窗口中，INPUTFILES一欄中選擇那十個(gè)原始觀測文件；RESULTFILES一欄，這里把輸出的測站信息文件命名為08_221.STA。點(diǎn)擊Run,就會(huì)從原始觀測文件的頭文件中提取測站名稱、接收機(jī)和天線相關(guān)的信息，存儲(chǔ)在Olympic\STA目錄下的08_221.STA文件中。

之后找到08_221.STA文件，用記事本打開，修改TYPE001:RENAMINGOFSTATIONS下的內(nèi)容，具體做法是將TYPE002:STATIONINFORMATION^白勺STATIONNAMEFLGFRO和TO內(nèi)容不做修改拷貝上去，之后的其他內(nèi)容留空。如圖所示。編輯完之后保存關(guān)閉。t.?和E鄭SjJ際⑴立0士配三配T/巨Ml:MHJUZfCGFVTkTICfflSSTATICKNAME『LG跟IIS3眥腰濫Esffi1#1ABglx11aluo&oooooo-uMrBLBo就陰srsoslBCsnElSTATICKNAME『LG跟IIS3眥腰濫Esffi1#1ABglx11aluo&oooooo-uMrBLBo就陰srsoslBCsnEl-rmhs-6.Vbit_RPQ-0.-M-HM-g^1.!■VL*41w^-11-LHl時(shí)MfliF即勒此尸MEHooooooooooHln心aQ。aa□u?u一■THTd.?i"I..hv--h-ajDooooooooooYYYYIT即HHWSF**=+++***+*;2W12Wl閱如舊20991231(W00OO2M91231?M002。則123100Eg20931231弼。12212O9&1231000Q(0(W00000000DCEOT912310000OO20991231(MJ00002M51231EKi00OdTOOLDSTATIOMHMEKSTA7PJUNME&*T*4*Hm***kk&JFSr*€UAOK3iLEMSHAD?LG*44H￥TWI1制娜LfflO001L319SC001■9B0001IO鼬010101010101.訓(xùn)QQ0(i