GAMIT-GLOBK數據處理報告
一．處理任務
利用GAMIT-GLOBK軟件對2011年年積日為94天的shao, lhaz, xian, kunm, bjfs, urum共6個IGS測站的GPS測量數據進行處理，並對處理結果進行評估。
二．處理步驟
-
安裝虛擬機和Linux系統
在win7系統下安裝VMware Workstation v9.0.1虛擬機，並在其上安裝ubuntu 12.10。
安裝完並登陸ubuntu后，在terminal中運用以下指令安裝csh:
sudo apt –get install csh
安裝完成后在terminal執行：
chsh
按提示輸入：
/bin/csh
修改bash為csh。重啟ubuntu系統，打開terminal，輸入：
sudo –sH
輸入密碼后進入root權限，然后輸入：
echo $SHELL，
返回
/bin/csh，說明SHELL修改成功。
安裝庫文件libX11 –dev和編譯器gfortran:
sudo apt –get install libX11 –dev
sudo apt –get install gfortran
-
安裝GAMIT-GLOBK軟件
在/opt目錄下建立gamit目錄作為軟件的安裝目錄，將下載的GAMIT-GLOBK+10.40.zip安裝包文件復制到/opt/gamit目錄下。打開terminal，進入gamit文件夾下，輸入命令：
unzip GAMIT-GLOBK+10.40.zip
解壓得到source文件夾，在terminal中進入source文件夾，運行：
./install_software
如果此時提示權限不足，則輸入
chmod 777 *
為所有用戶和文件設置可讀寫權利，然后運行./install_software命令，當到達提示
時，先進入source/libraries文件夾下用gedit打開Makefile.config，找到第283行：
# ----- for Linux from 0.0.1 to 3.0.0 -- #
OS_ID Linux 0001 3000
將其改為
# ----- for Linux from 0.0.1 to 4.0.0 -- #
OS_ID Linux 0001 4000
這是由於Linux內核版本已更新至3以上；另外在此文件的105行將
MAXSIT 60
MAXSAT 32
MAXATM 13
MAXEPC 2880
的MAXATM修改為25，MAXEPC修改為6000，關閉並保存文件后在terminal里選y繼續安裝，即可完成GAMIT-GLOBK軟件的安裝。
安裝完成后，需要修改csh.cshrc文件，如下圖提示所示：
進入/home/etc文件夾下，用gedit csh.cshrc打開，在endif前加入
stenv PATH /opt/gg/gamit/source/gamit/bin: /opt/gg/gamit/source/gamit/com:
/opt/gg/gamit/source/gamit/kf/bin:$PATH
stenv HELP_DIR /opt/gg/gamit/source/help/
stenv INSTITUTE GNSS CENTRER
保存並關閉文件后退出，在terminal中輸入doy指令，若出現結果則表明GAMIT-GLOBK已成功安裝。
-
工程文件和待處理數據准備
在任意目錄下建立工程文件夾x094，文件名為四字小寫，里面建立brdc，igs，rinex和年積日94四個文件夾：
需要下載的數據有rinex觀測值文件，廣播星歷和sp3星歷，分別放在rinex，brdc和igs文件夾下，可以用sh_get_rinex，sh_get_nav和sh_get_orbits自動下載，也可以通過網站http://cddis.nasa.gov/gnss_datasum.html#daily網站手動下載，本例數據從網上下載，准備好的工程文件內容如下（tree指令可通過sudo apt-get install tree安裝）：
-
tables文件夾更新和准備文件制作
需要每年更新的tables文件有：
章動表nutabl.、太陽表soltal.、月亮表luntab.、跳秒表leap.sec
需要每月更新的文件有：
編碼文件dcb.dat
需要每日更新的文件有：
地球自轉參數ut1.、pole
在有新的接收機或者天線出現時需要更新rcvant.dat，在有新的衛星發射時需要更新svnav.dat。
這些數據可以在ftp://garner.ucsd.edu/pub/gamit/tables/下載，在所有更新完成后，在工程目錄x094下運行sh_setup -yr 2011將tables表鏈接到工程目錄中。
需要准備的文件有：
測站信息文件station.info
將tables文件夾里的station.info文件復制到rinex目錄下，打開並編輯該文件，只保留以*號開始的行信息（文件頭），保存並關閉該文件。使用terminal進入rinex目錄，輸入
sh_upd_stnfo –files *.11o
就可以生成station.info文件：
近似坐標文件lfiles.
近似坐標從觀測o文件中提取，terminal進入rinex目錄，輸入
grep POSITION *.11o >lfiles.rnx
抓取o文件中的近似坐標值到lfiles.rnx文件中，得到測站近似坐標后，采用命令
rx2apr lfiles.rnx 2011 094
得到apr格式的文件lfiles.rnx.apr，然后通過指令
gapr_to_l lfiles.rnx.apr lfiles. 2011 094
將apr文件轉換為lfiles.：
將station.info和lfiles.文件復制到工程文件夾x094下的tables目錄下。
進入工程文件夾x094下的tables目錄，目錄控制文件process.defaults、測站列表文件sites.defaults、解算控制文件sestbl.、精度控制文件sittbl.均保持默認值。在本次處理過程中，不考慮海洋潮汐的影響，可以禁用海洋潮汐模型otl.grid文件，在sestbl.文件中設置為：Tides applied = 23; use otl.grid = N，保存並關閉sestbl.文件
-
數據處理
進入工程目錄x094下，采用批處理命令
sh_gamit –d 2011 094 –expt x094 –noftp
處理過程截圖如下：
-
處理結果及精度評定
最終處理結果保存在qx094a.094文件中。
運用globk進行點位平差處理，首先terminal進入工程目錄x094下，運行
sh_glred –d 2011 94 –expt x094 scal –opt H G
對基線解算數據進行預處理，然后進入gsoln目錄，運行
ls ../glbf/h*glx>x094.gdl
globk 6 globk_comb.prt globk_com.log x094.gdl globk_comb.cmd
處理截圖如下：
