GAMIT學習資料整理

GAMIT10.4安裝（基於Ubuntu 10.04）

1、安裝虛擬機器VMware、Ubuntu10.04；

2、進入終端輸入：sudo passwd root 為root使用者建立密碼，並以root使用者登入，或sudo su回車；

3、系統更新、漢化；

4、安裝gcc：# apt-get install gcc ；

5、安裝csh： # apt-get install csh 修改bash為csh，重新啟動；

6、安裝gfortran ： # apt-get install gfortran ；

7、安裝libx11-dev庫支援 # apt-get install libx11-dev ；

8、修改shall為bash，重新啟動，並設定路徑：

回到使用者根目錄，開啟.bashrc ，在最後加上如下程式碼即可

export PATH=$PATH:/opt/gamit/gamit/bin:/opt/gamit/com:/opt/gamit/kf/bin export setenv HELP_DIR=/opt/gamit/help/

9、將gamit安裝包放在目錄opt/gamit/下

進入目錄 ：# cd /opt/gamit

10、修改install_software檔案內容 ： # gedit install_software ， 開啟install_software檔案，在檔案的中下部修改“usr -name libX11.a”為"usr -name libX11.so" 。（動態共用庫）

11、執行install_software ： # ./install_software， 按提示輸入兩次Y後，修改makefile.config ，在/opt/gamit/libraries裡，修改Makefile.config中的一組引數

（1）MAXSIT 55 、MAXSAT 32 、MAXATM 25 、MAXEPC 5760

（2）# Specific to FC5（F6,F7,F8 ）

然後，在輸入兩次Y完成安裝；

12、安裝完後，開啟終端輸入：doy，檢視程式是否已正確安裝。

(註：在root使用者下安裝完後,回到自己使用者不能使用，

原因:需要修改軟體的許可權（chmod 777 軟體名)

、資料準備階段

（個人理解：這一過程可以先在windows下準備好，到時候再拷到linux的工程資料夾）

1.觀測檔案標準化：

1）觀測檔案o檔案的檔名要小寫，觀測檔案裡的點名一般要大寫。

2）認真檢查觀測檔案內接收機型別、天線型別並檢視GAMIT的ant.dat/rcvant.dat中是否有相應的型別，否則需要更新tables中的檔案。

3）天線高及天線測量方式。

2.需要從網上下載的資料檔案

1）http://www.ngs.noaa.gov/CORS/Gpscal.html 查詢觀測 0-files 檔案所在的GPS周天，年積日與 UTC 之間的關係。

2）導航電文檔案brdcdoyn.11n、精密星曆*.sp3檔案以及用到的IGS跟蹤站 O-files資料。(導航電文檔案下載autodoyn.11n；注意：檔名都要改成小寫)

3.一些常見表檔案說明

dcb.dat:碼相關型接收機偽距改正引數統計表；(缺少這個檔案基線處理最後一部會出錯)

antmod.dat：天線相位中心引數檔案；

rcvant.dat：接收機及天線名稱對照表；

guess_rcvant.dat: (如果rcvant.dat中沒有rinex標頭檔案中的接收機或天線型號，將從該檔案中獲取，使用者可以自己修改)（該檔案在處理過程中會自己產生，一般沒用不需要準備） otl_FES2004.grid: GAMIT10.4中的潮汐改正檔案，696M

註：地球固體潮改正和海洋負荷潮改正：(詳見《衛星導航定位新技術及高精度資料處理方法》p7-p8)

地球固體潮改正：攝動天體（月亮、太陽）對彈性地球的引力使地球表面產生周期性的漲落，稱為固體潮現象。在小於100KM的GPS相對定位中，兩個測站的固體潮影響幾乎是相同的，在差分過程中可抵消，因此可不考慮此項改正。

海洋負荷潮改正：由海洋潮汐的周期性漲落所引起的。與固體潮類似，但小一個數量級。 需要每週更新的：

ut1.： 地球自轉引數； ln –s ***

pole.：極移引數； ln –s ***

需要每年更新的：

leap.sec：跳秒表；

soltab.：太陽星曆； ln –s ***

luntab.：月亮星力； ln –s ***

nutabl.：章動參數列； ln –s ***

gdetic.dat：各種大地坐標系引數；

pmu.usno: ln –s ***

grid.otc: ln –s ***（這個好像有問題，GAMIT10.4中不是這個檔案）？ svs_exclude.dat：需要剔除衛星列表；(沒有的話，處理過程中會提示waring，但也可以計算)（註：應該是每次計算都需要更新的，個人理解）

vg_in：測站坐標速度表；

autcln.cmd：（註：如果沒有的話，資料處理的最後一步 csh btest1.bat不能正常進行） 需要重點準備的：

sestbl：資料處理引數設定表；

sittbl：測站約束表；

station.info：測站資訊表；

lfile.：站初始坐標（球坐標）；（自己製作）

、檔案準備階段

1.檔案結構整理

test：工程目錄，在其下面包括：

1）tables/存放各種表檔案；

2）igs/存放精密星曆*.sp3檔案；

3）rinex/存放觀測資料檔案*doyn.11檔案；

4）brdc/存放廣播星歷檔案brdcdoy0.11n檔案；

5）doy1/ doy2/ ??/存放各個時段解算的檔案及連結檔案。

2.需要重點準備的4個檔案

1) station.info:

記錄各測站的站點、站點名稱、開始觀測時間、結束觀測時間、天線高、 天線量測方法、接收機型別、軟硬體版本號、天線型別等資訊。嚴格按照已經給出的跟蹤站格式編輯。(注意：要現在當前目錄下拷入station.info檔案，把內容刪掉留下前面的幾行，否則後面的步驟沒法進行)

具體操作方法有二：

1、手動操作，按照相應檔案格式對齊操作；可以從樣本格式開始編輯。

2、自動提取觀測檔案的有關資訊，生成 station.info 檔案。

在各時段檔案目錄下操作：

sh_upd_stnfo –files site1doy0.yyo site2doy0.yyo –orbt igsf –u –ref ../tables/station.info

執行上述命令後，就自動從觀測檔案 0-files中提取重要資訊，生成 station.info 檔案。

我對station.info檔案的準備（武大郭博峰）：

1、rinex資料夾下新建資料夾all_station(存放所有點)

2、工程觀測檔案(o-files)複製到all_station，框架網（o-files）複製到all_station

3、進入Terminal

cd test/rinex/all_station

sh_upd_stnfo –files *.11o -ref station.info(此檔案為軟連結檔案，從~/不常更新表/station.info 連結至指定資料夾下，且只有一個錶頭)

4、檢視station.info檔案是否正確

gedit station.info

5、將得到的station.info複製到lfiles資料夾下

2） l-file：

記錄各測站的近似坐標: Gamit要求是球面坐標，即l-file。將所有 o 檔案整理在一起，然後分別執行如下命令:

grep POSITION *.11o > lfile.rnx

rx2apr lfile.rnx 2011 250 (註：根據所有的o檔案中的x、y、z取平均值,檔名為lfile.rnx.apr)

gapr_to_l lfile.rnx.apr lfile. “ “ 2011 250

我對l-file檔案的準備（武大郭博峰）：

1、在rinex資料夾下新建資料夾all_data(存放所有點)

2、進入該資料夾：cd test/rinex/all_data

3、提取POSITION

grep POSITION *.11o >lfile.xyz

4、檢視lfile.xyz檔案是否正確

–s ../tables/soltab.2011.J2000 ./soltab.

ln –s ../tables/nutabl.2011 . /nutabl.

ln –s ../tables/pole.usno . /pole.

ln –s ../tables/ut1.usno . /ut1.

ln –s ../tables/pmu.usno . /pmu.

另外，潮汐改正檔案也是很容易出錯的檔案。以10.4版本為例， 預設的潮汐 改正檔名為otl.grid（此檔案為損壞檔案）（連結前是otl._FES2004.grid ）。一定要檢查一下 otl.grid 是否是真實連結，出現死連結的情況需要手動連結。命令為：

ln -s ../tables/otl._FES2004.grid ./otl.grid （我在算的時候沒有用這個檔案，下次試試用）

（注意：連結時，如果該資料夾中已有該檔案，無論連結成功與否，要重新連結時都要先刪除原檔案在進行連結。）

、資料處理階段

一、批次處理

1) 執行makexp 程式，生成輸入檔案

sh_makexp –expt test –orbt igsf –yr 2011 –doy 250 –sess 99 –srin –nav brdc2500.11n –apr lfile. –sinfo 15 00 00 5760

或者：直接執行 makexp，然後按提示輸入

2) 執行sh_sp3fit 指令碼，生成軌道初始根數；

（曾由於ut1.連結的不是最新的，導致錯誤。）

sh_sp3fit –f igs15623.sp3 –o igsf –d 2011 250 –r BERNE –t –u

3) 執行sh_check_sess 指令碼，檢查衛星一致性；

sh_check_sess –sess 250 –type gfile –file gigsf1.250

4) 執行makej 程式，生成衛星鍾差檔案；

makej brdc2500.11n jbrdc1.250

5) 執行sh_check_sess，檢查衛星一致性；

sh_check_sess –sess 250 –type jfile –file jbrdc1.250

6) 執行makex 程式，生成X 檔案；

makex test.makex.batch

7) 執行fixdrv 程式，生成批次檔；

fixdrv dtest1.250 （進行這步可能發生錯誤，提示

FATAL :121202:2337:41.0 FIXDRV/armake: Invalid T-file interval = 10 27 2012 10 27 2012 --need either T-file or one X-file to run FIXDR 修改fixdrv/armake.f 檔案中的2010為2099或9999）並在此目錄下執行 make -f Makefile 進行重新編譯。

8) 執行fixdrv 生成的批次檔

csh btest1.bat 執行此命令前先把GAMIT.FATAL 和 GAMTI.WARNING 刪除

二、 gamit 基線處理結果品質評價指標： (檢驗q檔案)

A）標準均方根 nrms解算結果為 0.25左右，若是 nrms大於 0.5 ，表明處理結果是有問題的（例如：周跳沒有修復；測站的起算點坐標有問題等）。若是nrms小於 0.5 ，則認為成功解算， 否則需要檢查原因，重新處理；

B) 檢視gamit基線解算精度，一般要求相對精度達到10-8~10-9

C) 檢查所有測站點是否都參與計算

說明：如果有問題，需要檢視

1.檢視station.info檔案： REC#、 ANT#、 HtCod

IGS跟蹤站的天線高都是到DHARP的。

外業觀測的需要問清楚當時量取天線的位置。

2.檢視sestbl檔案

1) 高度角 Elevation Cutoff(不能是0°，按規範來為15°)；

2) Choice of Observable 長基線選擇LC_help、短基線選擇L1L2 Independent

3) 其他小錯誤該正項