||||
1.文件准备
*RINEX观测数据(O,N,G)
每个测站接受机和天线信息(station.info 文件)
卫星列表和解算方案(session.info)
IGS观测站数据;若下载得到的是压缩格式的rinex d文件,可用GAMIT软件提供的sh_crx2rnx命令进行格式转换
*从SOPAC的FTP服务器下载以下共同表文件(如果你的数据观测日期晚于所使用的GAMIT/GLOBK软件的发布日期,这些文件则通常需要更新。下载这些文件后可以直接放入GAMIT/GLOBK安装目录的tables文件夹内,一劳永逸,避免以后重复下载):
antmod.dat 天线相位中心改正表
autcln.cmd 周跳自动探测和修改命令表
dcb.dat 差分码偏差表
gdetic.dat 地球形状参数表
hi.dat 仪器高参数表
leap.sec 跳秒表
luntal.2014.J2000 月亮历表
nutabl.2014章动表
ot1.list/grid 全球格网或各测站的海洋潮汐参数表
pole.usno 极移表
rcvant.dat 接收机及其天线型号对照表
soltab.2014.J2000 太阳历表
svnav.dat 星号对照表
svs_exclude.dat 坏卫星信息表
ut1.usno 国际时间系统表
2.创建工作目录
应包括以下目录($work代表你的要进行GAMIT解算的目录):
/\$work/$doy(具体的年积日);
/$work/brdc(广播星历);
/$work/igs(SP3文件);
/$work/tables(公用表);
/$work/rinex(观测文件)
3.参数文件的修改和设置
主要是解算控制文件(sestbl.和sittbl.),可以据模板设置。
如果打算在解算时使用VMF1映射函数(sestbl.文件中设置),则需要下载VMF1的格网数据文件(mit的ftp服务器)。并设置tables文件夹中的map.grid指向这个文件。
4.初始文件生成
4.1 station.info文件
cd 命令进入tables目录中,station.info文件中的文件头复制到刚新建的文档中,保存。输入命令:
sh_upd_stnfo -files ../rinex/*o
提取测站信息(这里的../rinex/*o是指rinex观测数据所在路径。)
4.2 l-file台站坐标文件
(方法一)批量提取rinex观测数据头文件中的测站位置作为先验坐标(屡试不爽)。
命令示例(在rinex目录中执行):
grep POSITION *o>lfile.rnx
rx2apr lfile.rnx $year $doy($year为具体的年份(4位),$doy为年积日,下同)
gapr_to_l lfile.rnx.apr lfile. \$year $doy
至此,lfile.文件已准备好了(可打开查看)
(方法二)生成L-file命令示例:
sh_rx2apr -site <site> -nav <nav> -ref <ref> -apr <apr> -chi<val>"
这里<site>是一个 RINEX 格式观测文件的文件名,即需要得到坐标的测站, <nav>是一个 RINEX 格式的导航文件,<ref>是一个已知坐标测站的 RINEX 文件名,作为参考站(靠近未知台站好,但没有必要必须靠近),<apr> 是一个 apr 文件的文件名,这个文件包括参考测站的坐标。<val> 是 chi-square 值,小于这个值的 svpos(point-position)解被认为是收敛的。
4.3 链接各项文件
在$doy目录下,输入命令(示例):
ln -s../rinex/*o ./(链接观测数据)
ln -s../igs/*sp3 ./(链接精密星历)
ln -s../brdc/*n ./(链接广播星历)
sh_links.tables -frame J2000 -year $year(链接tables下的公用表文件)(链接后,需检查是否成功,尤其是luntab. soltab. nutabl. pole. ut1.这几个文件,如果链接失败,可以使用ln -s 命令手动链接)。
5.开始解算
依次执行命令:
makexp #数据预处理,按提示键入解算任务名、广播星历名,采样间隔等
sh_sp3fit -f *sp3 -o igsf -d $year $doy #链接精密星历
makej brdc*n jbrdc\$y.$doy #生成卫星钟文件($y代表一位数字表示的观测数据年,即年的最后一位)
makex $expt.makex.batch #创建观测数据($expt代表你的解算任务名)
fixdrv d$expt\$y.$doy #启动数据驱动程序
csh b*bat #执行数据解算
6.备注
GAMIT/GLOBK软件提供了很多下载准备数据的脚本,如下载广播星历的sh_get_nav,下载精密星历的sh_get_orbits以及下载IGS参考站观测数据的sh_get_rinex等。读者可以在终端中输入这些命令参看帮助文档。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-5-10 18:04
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社