先說兩句
前久由於工作需要,花了兩個月時間來弄gamit解算;找了幾個同學都是半瓢水,後來得高人指點,總算是弄通了,其實在我看來,會使用gamit本身沒什麼,能寫一個出來才nb,所以很鄙視某些人把這個東西藏著捂著,或是說一半留一半的做法; ;所以這裡總結一下各步驟貼出來給大家參考,希望有用到的朋友少走彎路.
GAMIT相關閱讀：http://www.linuxidc.com/GAMIT/
GAMIT10.4下載在Linux公社5號FTP服務器，具體下載見 http://www.linuxidc.net/thread-1186-1-1.html
一.數據准備
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
