Gamit課程結(jié)課報(bào)告

1、精選優(yōu)質(zhì)文檔-傾情為你奉上研究生課程學(xué)習(xí)報(bào)告 GNSS數(shù)據(jù)處理軟件應(yīng)用 GAMIT軟件實(shí)踐 2016年12月17日專心-專注-專業(yè)Gamit基線解算一、Gamit介紹Gamit是美國(guó)麻省理工學(xué)院（MIT）和斯克里普斯海洋研究所（SIO）聯(lián)合開發(fā)的高精度GPS數(shù)據(jù)處理軟件，它主要用于分析研究地殼形變、精密單點(diǎn)定位等領(lǐng)域。GAMIT可以解算衛(wèi)星軌道、測(cè)站坐標(biāo)、大氣延遲、整周模糊度等，其最主要的特點(diǎn)是其解算精度高，且免費(fèi)開放源代碼，用戶可以根據(jù)需要對(duì)源程序做相應(yīng)的修改，以便于科研工作。GAMIT軟件所需的數(shù)據(jù)是Rinex格式的，可以處理各種不同型號(hào)的GPS接收機(jī)采集的數(shù)據(jù)。利用高精度的GAMIT軟件

2、可以確定地面站的三維坐標(biāo)，在利用精密星歷和高精度起算點(diǎn)的情況下，基線解的相對(duì)精度能夠達(dá)到10-9左右，解算短基線的精度能優(yōu)于1 mm，是世界上最優(yōu)秀的GPS軟件之一。目前，GAMIT軟件支持在Unix和Linux操作系統(tǒng)上運(yùn)行。二、Gamit安裝首先安裝Linux系統(tǒng)，這里用的是Ubuntu16.04,然后安裝Gmait所需要的一些編譯器：GUN gfortran 4.2、libx11-dev。然后將gamit10.5拷貝到/opt/下，在終端中用sudo chmod命令給予+x權(quán)限。進(jìn)入/opt/gamit10.5/libraries中，用vim編輯Makefile.Config文件，修改其

3、x11libpath和x11incpath路徑和操作系統(tǒng)版本號(hào)（在終端中輸入uname a可以查詢系統(tǒng)版本號(hào)）。最后輸入./install_software，開始安裝，耐心等待安裝完成。最后gamit會(huì)在終端提示：Dont forget to set your: path to include /opt/gamit10.5/gamit/bin and /opt/gamit10.5/kf/bin以及一個(gè)幫助文檔的環(huán)境變量：HELP_DIR environment variable in you shell profile (in .cshrc/.tcshrc add: setenv HELP_D

4、IR /opt/gamit105./help)這里我用的是bshell，所以在/.basrhc中將這些環(huán)境變量添加進(jìn)去，然后輸入source ./bashrc使環(huán)境變量立即生效最后輸入doy命令，測(cè)試是否安裝成功。三、數(shù)據(jù)處理3.1 table文件更新想要進(jìn)行高精度的數(shù)據(jù)處理，table文件是不可或缺的。主要有：地球定向EOF文件、衛(wèi)星鐘文件、潮汐文件、太陽(yáng)、月亮表、固體地球潮汐表、跳秒文件?？梢栽?pub/gamit/tables中更新。3.2數(shù)據(jù)準(zhǔn)備我們小組用北京房山、昆明、拉薩、烏魯木齊、韶山五個(gè)站的數(shù)據(jù)進(jìn)行處理。然后在ftp:/cddis.gsf

5、/pub/gps/data/daily/中下載這五個(gè)站的o文件和n文件，在/中下載igs精密星歷，我用的是2016年123天，第1895周，周一的數(shù)據(jù)。3.3鏈接文件首先需要建立一個(gè)工程，以4個(gè)字母為宜，這里為wang，然后在wang/下建立四個(gè)文件夾，分別是brdc、igs、123、rinex然后將tables拷貝到wang/下，分別把下載的o文件、n文件和sp3文件拷貝到rinex、brdc、igs中，然后在終端中用cd命令進(jìn)入到wang/123/下，用ln s命令將o文件、n文件和sp3文件鏈接到123文件夾下，然后在使用links.d

6、ay2016 123 wang命令，將tables所需的文件鏈接過來，然后輸入ls l檢查鏈接是否斷掉，主要檢查下表這幾個(gè)文件：源文件鏈接文件otl_FES2004.gridotl.gridvmflgrid.2016map.gridluntab.2016.J2000luntab.soltab.2016.J2000soltab.nutabl.2016nutabl.pmu.usnopmu.ANU_grid.atlatl.gridatmdisp_cm.2016atml.grid確認(rèn)無誤后，進(jìn)行下一步3.4 準(zhǔn)備控制文件首先將tables中的itrf08.apr拷貝到rinex中，用grep命令，將五

7、個(gè)igs站的坐標(biāo)拷貝到lfill.中，然后將lfile.鏈接到tables和123中。接下來更新測(cè)站列表，首先在sites.defaults中將五個(gè)站添加進(jìn)去，然后新建一個(gè)文本文件，輸入五個(gè)站的站名，用sh_upd_stnfo l命令生成，然后用sh_upd_stnfo files./rinex/*.16o生成新的。緊接著在sittbl.中將文件中沒有的五個(gè)igs站添加進(jìn)去，在sestbl.中設(shè)置如下：Choice of Experiment = BaseLine,Choice of Observable = LC_AUTCLN,Use otl.

8、grid = Y, Use met.grid = N, Use map.grid = Y, Use otl.grid=Y, Use atl.grid=Y。3.5分布處理在123目錄下執(zhí)行makexp，按照要求輸入項(xiàng)目名，軌道類型等參數(shù)，執(zhí)行完成后，處理精密星歷：sh_sp3fit f igs18951.sp3 o igsf t，執(zhí)行完成后進(jìn)行sh_check_sess sess 123 type gfile file gigs6.123文件檢查，然后處理廣播星歷：makej bjfs1230.16n jbjfs6.123，處理完成后用上述命令再檢查一次j文件，最后用makex，然后輸入1，生成

9、batch文件，在用fixdrv dwang6.123生成bat批處理文件，最后用csh bwang6.bat進(jìn)行基線解算，最后生成的q-、o-文件就是解算結(jié)果。四、基線解算4.1 解算結(jié)果首先看o-文件中的Postfit nrms為0.，小于0.25，解算質(zhì)量合格。然后接著查看基線解算結(jié)果，如下表（單位：m）：基線x中誤差y中誤差z中誤差LL中誤差BJFS_CHAN0.00520.00660.0058.94340.0041BJFS_LHAZ0.00640.01000.0066.00740.0062BJFS_SHAO0.00640.00790.0059.32830.0037BJFS_URUM0

10、.00660.00860.0069.78240.0062CHAN_LHAZ0.00720.01080.0077.39440.0071CHAN_SHAO0.00660.00840.0066.50910.0037CHAN_URUM0.00740.00950.0079.87770.0069LHAZ_SHAO0.00790.01130.0072.78770.0074LHAZ_URUM0.00610.01010.0066.72260.0038SHAO_URUM0.00820.01020.0076.09840.0072從表中可以看出x，y，z三個(gè)方向的基線解算誤差基本都小于1cm，y方向的誤差在1cm左右

11、，本次基線解算合格。它們的協(xié)方差如下表（單位：m）：基線xyzBJFS_CHAN-0.40903-0.602080.78241BJFS_LHAZ-0.14529-0.357270.81118BJFS_SHAO-0.48006-0.559810.78616BJFS_URUM-0.26825-0.307950.86264CHAN_LHAZ-0.16709-0.465070.77643CHAN_SHAO-0.45949-0.581290.74635CHAN_URUM-0.2703-0.407530.85069LHAZ_SHAO-0.26097-0.404570.83994LHAZ_URUM-0.03

12、105-0.064890.78802SHAO_URUM-0.3564-0.333070.838384.2 結(jié)果對(duì)比與第118天的小組解算結(jié)果比較，結(jié)果如下表（單位：m）：基線BJFS_CHAN-0.00350.01020.0051-0.0035BJFS_LHAZ-0.0003-0.0018-0.00500.0008BJFS_SHAO-0.00740.01940.01230.0004BJFS_URUM-0.0011-0.0036-0.0061-0.0023CHAN_LHAZ0.0032-0.0120-0.0101-0.0002CHAN_SHAO-0.00400.00910.00720.0007C

13、HAN_URUM0.0024-0.0138-0.0113-0.0017LHAZ_SHAO-0.00720.02100.01730.0011LHAZ_URUM-0.0008-0.0018-0.00120.0001SHAO_URUM0.0064-0.0230-0.01840.0001從上表可以看出，同組之間的解算記過最大在1.5cm左右，基本都在毫米級(jí)，因此解算結(jié)果合格。五 總結(jié)因?yàn)镚amit是在Linux下運(yùn)行的，以前沒有接觸過Linux系統(tǒng)，因此剛開始學(xué)習(xí)Linux和安裝Gamit時(shí)遇到了很多問題，后來經(jīng)過向老師請(qǐng)教和百度，終于安裝成功，但是一執(zhí)行Gamit命令總會(huì)提示命令找不到，后來在課上老師講到是因?yàn)榄h(huán)境變量沒有配置對(duì)。然后在剛開始解算時(shí)沒有頭緒，不知道怎么做，就帶這上課時(shí)的印象再去看Gamit下的例子，跟著跑