一.数据准备
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
GAMIT相关阅读:http://www.linuxidc.com/GAMIT/
GAMIT10.4下载在Linux公社5号FTP服务器,具体下载见 http://www.linuxidc.net/thread-1186-1-1.html
本文永久更新链接地址:http://www.linuxidc.com/Linux/2016-02/128249.htm