GAMIT 解算-詳細步驟

一.資料準備

1.更新相關的tables檔案;

2.精密星曆下載;//精密星歷不同階段產品開頭為:igu ,igr,igs;

3.廣播星曆下載;//下載融合後的星曆;一般是以brdc開頭的

4.IGS觀測站資料下載;//有些站在某些伺服器上沒有,需要從多個伺服器下載資料;若下載得到的是壓縮格式的o檔案,需用gamit軟體提供的crx2rnx命令進行格式轉換;

二.核心檔案設定

1.lfile.

1.1.提出各測站的o檔案頭裡的XYZ位置;

1.2.從網站下載ITRF框架坐標檔案[ITRF2000_GPS.SSC.txt];裡面有其84下XYZ坐標 ;

1.3.在Ubuntu下用tform工具將其轉換為球坐標;

2.station.info //注意天線高的設定

2.1.從網站下載station.info檔案,裡面有各測站某時段內使用的天線的引數;

2.2.從完整的station.info檔案裡面copy出IGS站的條目,並結合當地使用的天線輸入新的條目;

3.sestbl.

Type of analysis =0-iter

;//疊代多少次[計算結果作為初始值再次計算]

Choice of experiment =relax. //要想點位置精度高用relax. 若目的是求基線後面平差則用baseline;

Zenith delay =N ;//對於短基線

4.sittbl.

4.1.下載sittbl.refined 檔案,裡面有最近的各站的可靠性的一個評測;

4.2.製作sittbl檔案,用以上檔案條目+各站初始位置的精度條目;

要想點位置精度高則設定IGS站的約束小一些,即認為IGS站已經很準,強約束,若目的是求基線後面平差則可以將IGS站的約束設定大些;

三.解算過程

以下為解算的bat檔案:

//1.建立工作目錄(呼叫各種檔案到此處)

mkdir workspace

cd workspace

rm -f *

//2.連結o檔案,sp3檔案,copy廣播星歷檔案;

ln -s ../rinex/*.09o .

ln -s ../igs/*.09o .

cp ../brdc/brdc$1* eexample9.$1

ln -s ../igs/igs$2.sp3 .

//3.連結相關tables表檔案

links.day 2009 $1 example

//4.建立解算專案

sh_makexp -expt example -orbt example -yr 2009 -doy $1 -sess 0 -nav eexample9.$1 -sinfo 30 00 00 2880

//5.軌道積分

sh_sp3fit -f igs$2.sp3 -o example

//6.檢查g檔案

sh_check_sess -sess $1 -type gfile -file gexample9.$1

//7.建立j檔案

makej eexample9.$1

jexample9.$1

//8.檢查g檔案

sh_check_sess -sess $1 -type jfile -file jexample9.$1

//9.生成觀測檔案

makex example.makex.batch