PurpleSky的个人博客分享 http://blog.sciencenet.cn/u/PurpleSky

博文

使用GAMIT软件分步进行基线解算

已有 9941 次阅读 2015-4-2 12:00 |个人分类:GAMIT/GLOBK|系统分类:科研笔记|关键词:学者| gamit, 基线解算

1.文件准备

  • *RINEX观测数据(O,N,G

  • 每个测站接受机和天线信息(station.info 文件

  • 卫星列表和解算方案(session.info

  • IGS观测站数据;若下载得到的是压缩格式的rinex d文件,可用GAMIT软件提供的sh_crx2rnx命令进行格式转换

  • *SOPACFTP服务器下载以下共同表文件(如果你的数据观测日期晚于所使用的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/igsSP3文件);

  • /$work/tables(公用表);

  • /$work/rinex(观测文件)


3.参数文件的修改和设置

主要是解算控制文件(sestbl.sittbl.),可以据模板设置。

如果打算在解算时使用VMF1映射函数(sestbl.文件中设置),则需要下载VMF1的格网数据文件(mitftp服务器)。并设置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 值,小于这个值的 svpospoint-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等。读者可以在终端中输入这些命令参看帮助文档。




https://m.sciencenet.cn/blog-1365943-879244.html

上一篇:配置CentOS系统下的GAMIT自动发邮件
下一篇:使用GAMIT进行批处理基线解算

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-10 18:04

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部