GAMIT使用手冊

1、GAMIT/GLOBK軟件使用手冊一 軟解介紹GAMIT軟件最初由美國麻省理工學(xué)院研制, 后與美國SCRIPPS海洋研究所共同開發(fā)改良。該軟件是世界上最優(yōu)秀的GSP定位和定軌軟件之一, 采用精密星歷和高精度起算點時, 其解算長基線的相對精度能到達10-9量級, 解算短基線的精度能優(yōu)于1mm, 特點是運算速度快、版木更新周期短以及在精度許可范圍內(nèi)自動化處理程度高等, 因此應(yīng)用相當(dāng)廣泛。GAMIT軟件由許多不同功能的模塊組成, 這些模塊可以獨立地運行。按其功能可分成兩個局部: 數(shù)據(jù)準(zhǔn)備和數(shù)據(jù)處理。此外, 該軟件還帶有功能強大的shell程序。目前，比擬著名的GPS數(shù)據(jù)處理軟件主要有美國麻省理工學(xué)院

2、MIT和海洋研究所SIO聯(lián)合研制的GAMIT/GLOBK軟件、瑞士伯爾尼大學(xué)研制的BERNESE軟件、美國噴氣推進實驗室JPL研制的GIPSY軟件等。GAMIT/GLOBK和BERNESE軟件采用相位雙差數(shù)據(jù)作為根本解算數(shù)據(jù)，GIPSY軟件采用非差相位數(shù)據(jù)作為根本解算數(shù)據(jù)，在精度方面，三個軟件沒有明顯的差異，都可得到厘米級的點位坐標(biāo)精度。相比擬而言，GIPSY軟件為美國軍方研制的軟件，國內(nèi)只能得到它的執(zhí)行程序，在國內(nèi)，它的用戶并不多，BERNESE軟件需要購置，它的用戶稍微多一點，GAMIT/GLOBK軟件接近于自由軟件，在國內(nèi)擁有大量用戶。 GLOBK軟件核心思想是卡爾曼濾波卡爾曼濾波理論是

3、一種對動態(tài)系統(tǒng)進行數(shù)據(jù)處理的有效方法, 它利用觀測向量來估計隨時間不斷變化的狀態(tài)向量，其主要目的是綜合處理多元測量數(shù)據(jù)。GLOBK的主要輸人是經(jīng)GAMIT處理后的h-file和近似坐標(biāo), 當(dāng)然,它亦己成功地應(yīng)用于綜合處理其它的GPS軟件如Bernese和GIPSY產(chǎn)生的數(shù)據(jù)以及其它大地測量和SLR觀測數(shù)據(jù)。GLOBK的主要輸出有測站坐標(biāo)的時間序列、測站平均坐標(biāo)、測站速度和多時段軌道參數(shù)，GLOBK可以有效地檢驗不同約束條件下的影響, 因為單時段分析使用了非常寬松的約束條件，所以在GLOBK中就可以對任一參數(shù)強化約束。GAMIT/GLOBK和BERNESE采用雙差作為數(shù)據(jù)分析的根本觀測量，它們的

4、缺陷是不能直接解算鐘差參數(shù)，只能給出測站的基線結(jié)果，除測站坐標(biāo)參數(shù)之外，這些軟件還可以解算的參數(shù)有：衛(wèi)星軌道參數(shù)、衛(wèi)星天線偏差、光壓參數(shù)、地球自轉(zhuǎn)參數(shù)、地球質(zhì)量中心變化、測站對流層延遲參數(shù)、電離層改正參數(shù)等，這使這些軟件的應(yīng)用從大地測量學(xué)已逐漸延伸到地球動力學(xué)、衛(wèi)星動力學(xué)、氣象學(xué)以及地球物理學(xué)等領(lǐng)域，并取得了很多成果。GAMIT軟件的運行平臺是UNIX操作系統(tǒng)，目前，它可在Sun、HP、IBM/RISC、DEC、LINUX等基于intel處理器的工作站上運行。軟件可處理的最大測站和衛(wèi)星數(shù)目可在編譯時設(shè)定。它的根本輸出文件是H-文件，可作為GLOBK軟件的輸入文件，進而估計測站坐標(biāo)與速度、衛(wèi)星軌

5、道參數(shù)和地球定向參數(shù)。數(shù)據(jù)處理前，用戶需準(zhǔn)備所需要的文件，如測站先驗坐標(biāo)文件L-文件和vg-in文件、播送星歷文件、觀測數(shù)據(jù)文件以及其他輔助文件等。GAMIT每個時段觀測數(shù)據(jù)要求的周期最長為1個UTC天，即從UTC的0點到24點北京時間8:0024:00，原那么上不要跨天作業(yè)。GAMIT軟件的組成結(jié)構(gòu)見圖1.1，它由不同功能模塊組成，主要包括數(shù)據(jù)準(zhǔn)備、生成參考軌道、計算殘差和偏導(dǎo)數(shù)、周跳檢測與修復(fù)、最小二乘平差等模塊，這些模塊即可以單獨運行，也可以用批處理命令聯(lián)在一起運行，最大限度地減少人為操作，提高運算效率。軟件的執(zhí)行程序放在/com、/kf/bin和/gamit/bin三個目錄下。圖1.1

6、 GAMIT/GLOBK軟件組成結(jié)構(gòu)圖二 批處理實例確認GAMIT/GLOBK軟件安裝成功之后，為了能正確處理數(shù)據(jù)和繪圖，還需要做兩件事情。第一，更新表文件。在ftp效勞器上下載最新的tables表文件/pub/gamit/tables/，然后復(fù)制到安裝目錄下的tables中，遇到相同文件選擇覆蓋即可；第二，安裝并配置好繪圖工具GMTThe Generic Mapping Tools。下面的實例中，處理的數(shù)據(jù)為2021年第56到60共5天的3個IGS站BJFS、和WUHN、KUNM和2個待求站chdu、pixi的數(shù)據(jù)。一、處理前的準(zhǔn)備1、在主文件夾內(nèi)新建

7、test工程文件夾，工程內(nèi)新建brdc、igs、rinex三個文件夾，分別存放當(dāng)天的播送星歷、精密星歷和觀測值文件。注意：如果文件采用.z壓縮格式，那么應(yīng)當(dāng)使用gunzip命令進行解壓；如果解壓出的觀測值文件仍采用.d的壓縮格式，那么應(yīng)使用crx2rnx命令將其轉(zhuǎn)化為.o的標(biāo)準(zhǔn)RINEX格式sh_crx2rnx -f *.*d可用于.d到.o文件的批量轉(zhuǎn)換。2、終端進入test工程文件夾內(nèi)鏈接tables，運行：   sh_setup -yr 20213、生成文件。將test/tables下的文件拷貝到rinex文件夾下，翻開

8、并編輯，僅保存以#或*開頭的前幾行，保存并關(guān)閉。翻開終端并進入rinex文件夾，運行   sh_upd_stnfo -files *.10o運行成功以后文件便存放了此工程內(nèi)的開始、結(jié)束時間，站名和接收機、天線類型等信息。4、建立lfile.文件。lfile.是測站的先驗坐標(biāo)文件。翻開終端并進入rinex文件夾，用批處理的方式生成lfile.文件可以分為三步：a提取觀測值.o文件的先驗XYZ坐標(biāo)   grep POSITION *.10o > lfile.rnxb將.rnx文件轉(zhuǎn)化為.apr文件   rx2

9、apr lfile.rnx 2021 056c由.apr文件生成lfile.文件   gapr_to_l lfile.rnx.apr lfile. 2021 056將和lfile.兩個文件拷貝到test/tables文件夾，覆蓋原文件。5、sestbl.的配置。test/tables下的sestbl.文件是測段分析策略文件，該文件內(nèi)部有詳細的說明。一般來說采用默認配置即可，通常需要修改的三個地方是：   Choiceof Experiment       選擇處理方式   

10、;Choice of Observable      選擇觀測值類型和模糊度解算   Use otl.grid              選擇是否使用潮汐文件對于Choice of Experiment，選擇BASELINE時將固定軌道并在GAMIT處理中和輸出h-文件時忽略軌道參數(shù)；選擇RELAX時將采用松弛解，合并全球IGS h-文件時需要。要想點位置精度高用RELAX；假設(shè)目的是求基線后面平差那么用BASELINE。在此實例中采用默認的BASELINE。對于Choice

11、 of Observable，選擇LC_AUTCLN為采用寬巷模糊度值并用偽距在autcln中解算；對于小于幾公里的基線，用L1和L2獨立載波相位觀測值L1,L2_INDEPENDENT或者僅用LIL1_ONLY，相比用無電離層組合LC_HELP可以減少噪聲水平。對于Use otl.grid，由于這里已在ftp上更新tables，有了最新的otl.grid文件，所以這里選擇Y。6、sittbl.的配置。test/tables下的sittbl.文件對各個測站的先驗坐標(biāo)或鐘差、大氣模型等進行約束。對高精度的坐標(biāo)采取強約束，而對待求點采用松弛約束。如IGS站的坐標(biāo)分量約束在較小的175px，對未知點

12、的約束可以到510m。7、sites.defaults和process.defaults的配置。sites.defaults文件用來控制需要參與解算的測站。在文檔末尾可以根據(jù)提示編輯，來給定那些不參與解算的測站或是測站的某些天。process.defaults文件用來控制處理過程中的很多細節(jié)，比方sampling interval, number of epochs, start time for processing，default globk .apr file等等，根據(jù)需求和提示進行編輯。在這里均使用它們的默認值。二、利用GAMIT解算基線用終端進入test的工程文件夾，輸入批處理命令進

13、行解算：     sh_gamit -expt test -s 2021 056 060 -noftp -dopt D ao c x >& sh_gamit.log參數(shù)說明：   -expt：指定四個字符的工程名稱   -d：指定需要處理的指定日期，例如-d 2021 56 60，指的是處理2021年第56和60天。   -s：指定需要處理的時間序列，例如-s 2021 56 60，指的是處理2021年第56到60天。   -orbit：衛(wèi)星軌道類型。 

14、;  -yrext：給日目錄前添加年前綴，例如2021_006。   -noftp：處理過程中不連接ftp下載數(shù)據(jù)。   -copt：數(shù)據(jù)處理完成后待壓縮的文件類型，例如-copt o q m k x。   -dopt：數(shù)據(jù)處理完成后待刪除的文件類型，例如-dopt D ao c x。結(jié)果文件將存放在名稱為年積日的文件夾內(nèi)，此例中，可供參考的結(jié)果文件為：   sh_gamit_2021_006.summary        解算總結(jié)   qtes

15、ta.006                       解算記錄   otesta.006                       解算記錄的簡略版，一般關(guān)注此文件   htesta.10006                

16、60;   協(xié)防差矩陣、參數(shù)平差值在上述文件中，可根據(jù)描述來判斷解算結(jié)果是否符合相應(yīng)的需求。其中，基線解算結(jié)果O文件，即這里的otesta.006中的postfit_nrms項優(yōu)于0.3左右時最正確；如果大于1.0,那么表示此解存在問題。三、利用GLOBK進行平差處理用終端進入test的工程文件夾內(nèi)，運行：   sh_glred -expt test -s 2021 050 2021 065 -yrext -opt H G E >& sh_glred.log參數(shù)說明：   H：運行htoglb，把文件轉(zhuǎn)換為二進制文件&#

17、160;  G：運行g(shù)lred，合并文件   E：繪圖運行成功后，進入gslon文件夾，以psbase開頭的文件即為各個站三維坐標(biāo)時間序列的圖形表示。用終端進入gsoln目錄下，依次運行：   ls./glbf/h*glx > test.gdl   glred6 globk_comb.prt globk_comb.log test.gdl globk_comb d生成的globk_文件中就包含了解算點的三維坐標(biāo)和相關(guān)參數(shù)。至此，利用GAMIT-GLOBK解算GPS基線并進行平差的實例表達完畢。三 分布處

18、理實例1. 在桌面上新建一個test的文件夾，在test中新建 名稱為 brdc播送星歷、igs存放精密星歷、rinex下載的觀測值文件和自己要解算的觀測值文件。2. 用doy在test的終端里測試，比方輸入：doy 2021 02 01 查看要解算的是GPS第幾周，第幾天，以便下載數(shù)據(jù)用。3. 在終端輸入gftp，下載相關(guān)的文件，一定要選對自己解算需要的數(shù)據(jù)。一般下載*.d.z文件比擬小，好下載。下載完所需要的數(shù)據(jù)后，要解壓 ：語句：gunzip*.Z4. 新建一個文件夾，比方：名字為032，把要解算的測站o文件和n文件，以及精密星歷文件放進來。把032和test進行如下鏈接。5. 翻開ri

19、nex文件，輸入下面的語句，會產(chǎn)生一個lfile.rnx文件： 查看lfile.rnx文件中各點的坐標(biāo)是否正確：不能出現(xiàn)坐標(biāo)值為0。6. 輸入rx2apr lfile.rnx 2021 032,會多出來一個名為：lfile.rnx.apr的文件，接著輸入：gapr_l lfile.rnx.apr lfile. “ 2021 032 生成一個用緯度和精度來表示點的坐標(biāo)的文件lfile.,并把它復(fù)制到tables文件夾。7. 在tables文件夾中新建一個文件 sitelist，在里面輸入?yún)⒖颊军c的名稱，每行的第一列是空格，從第二行開始輸。 8. 在tabie文件夾啟動一個終端，輸入：sh_upd

20、_stnfo -l sitelist 。9. 會生成一個.new文件，把它給成 保存后，再輸入：sh_upd_stnfo -files ./rinex/*13o10. 在032文件中查看的鏈接是否正確，然后做一個l文件的鏈接，翻開終端輸入： ln -s ./tables/lfile ./ 到此準(zhǔn)備工作完成。11. 在032文件翻開終端，輸入：makexp 然后按照以下輸入。12. 按照提示輸入：歷元間隔這里是30s，開始時間，如果從00 00開始一天正好是2880個歷元。精密星歷13. 接下來按照提示輸入：sh_sp3fi

21、t -f igs 17255.sp3 o igsf t 利用 sp3 精密星歷生成 g 文件上面已經(jīng)下載好了 ，并且同時生產(chǎn) t-和 u-文件14. sh_check_sess sess 32 type gfile -file gigsf3.032下載的導(dǎo)航文件的名稱3指的是年的最后一位數(shù)字，032是新建的文件名稱檢查 g-文件15. makej auto0320.13n jauto3.032 利用導(dǎo)航文件生成 j-文件16. sh_check_sess sess 32 type jfile file jauto3.032檢查 j-文件17. makex test.makex.batch 生產(chǎn) x-文件建立的文件名稱18. fixdrv dtest3.032 生成批處理文件19. csh btest3.bat 執(zhí)行批處理腳本，得到 q-文件，即解算結(jié)果。四 精度評估1GPS常用評價站坐標(biāo)精度的指標(biāo)是多時段基線重復(fù)性和多時段坐