Gamit 数据处理,相关的文件配置

一、处理前准备

1、在主文件夹内新建test项目文件,项目内新建brdc、igs和rinex三个文件夹,分别存放广播星历,精密星历几观测值文件,所用的命令分别为sh_get_nav、sh_get_orbits和sh_get_rinex (若文件为.Z,用gunzip命令解压,若仍为d,用命令sh_crx2rnx -f 命令解压为o文件)
2、进入test项目文件夹,链接tables,运行

sh_setup -yr 2015 -apr itrf08_com.apr

3、生成station.info文件。
将test/tables下的station.info文件拷贝到rinex文件夹下,打开并编辑,仅保留以#或*开头的前几行,保存并关闭。
打开终端进入rinex文件夹,运行

sh_upd_stnfo -files *.150

4、建立lfile.文件(先验坐标文件)。(非IGS站时,才需此步)
打开终端进入rinex文件夹,用批处理的方式生成lfile.文件,分为三步:

(1) 提取观测值文件的先验坐标:
   grep POSITION *.15o > lfile.rnx
(2) 将.rnx文件转化为.apr文件
   rx2apr lfile.rnx 2015 001
(3) 由.apr文件生成lfile.文件
   gapr_to_l lfile.rnx.apr lfile. "" 2015 001
将station.info和lfile.两个文件复制到test/tables下,覆盖原文件。

5、sestbl.文件的配置。

  首先配置 sestbl.这里列出主要需要配置的几项内容:
   Choice of Experiment = BASELINE
	如果是解算基线的话建议选择 BASELINE,即固定轨道
	Choice of Observable = LC_AUTCLN
	选择观测值类型,如果是双频接收机并且是大于 5KM 的基线话建
	议选取 LC_AUTCLN,如果是小于 5KM 的基线也可以选择 LC_AUTCLN,
	其中 LC_HELP 和 LC_AUTCLN 的差别不大,LC_AUTCLN 略好些。
	Zenith Delay Estimation = Y
	Interval zen = 2
	对流层参数的设置,如果是短基线(大概小于 40KM)可以设置
	成 N,即不估计对流层天定延迟,采用双差的方法来消除对流层延迟。
	如果是大于 40KM 或者 50KM 的基线要设置成 Y。第二项是如果估计对
	流层延迟的话,设置多长时间估计一个对流层延迟参数,单位是小时。
	Tides applied = 31
	Use otl.list = Y
	Use otl.grid = Y
	设置潮汐改正参数,如果为 31 则考虑海潮改正,为 23 则不加海
	潮改正,如果考虑海潮改正的话,需要 tables 下有两个海潮文件,
	otl.list 和 otl.grid,其中前者是可有可无的,当两者都采用 Y 时,
	程序会根据测站位置首先利用 otl.list(只包括 IGS 站的海潮模型
	参数,以及位于该 IGS 站方圆 10KM 范围内都是一套参数)中的海潮
	参数,如果该测站没有出现在 otl.list 中,则会根据该测站的位置
	和 otl.grid 文件内差出该测站的海潮改正。
	以上为 sestbl.中主要的设置项,其他项对基线解算的影响很小,
	可以采用默认值。

6、sittbl.文件的配置。

下面设置 sittbl.文件,即设置各个站坐标的约束量,对其算点
	坐标采取紧约束,对非起算点的站采用松弛约束。
	执行
	vi sittbl.
	SITE
	打开 sittbl.文件。
	FIX
	--COORD.CONSTR.--
	<< default for regional stations >>
	ALL
	NNN30.
	30.
	30.
	<< IGS core stations >>
	BJFS BJFS_GPS NNN 0.050 0.050 0.05
	LHAZ LHAZ_GPS NNN 0.050 0.050 0.05
	TWTF TWTF_GPS NNN 0.050 0.050 0.05
	USUD USUD_GPS NNN 0.050 0.050 0.05
	对起算点的 XYZ 坐标分量约束为 5cm,非起算点的站采用松弛约束,
	约束量为 30 米。

7、sites.defaults文件的配置。

用来控制需要参与解算的测站,可定义不参与结算的测站。

8、process.default文件的配置。

用来控制处理过程中的细节。

二、理由GAMIT解算基线。

用终端进入test项目文件夹,输入批处理命令:
sh_gamit -expt test -s 2015 001 008 -orbit IGSF -noftp
参数说明:
	-d:处理数据的指定日期,-d 2015 001 008 ,处理001和008天的数据
	-s:处理数据的指定日期,-s 2015 001 008 ,处理001到008天的数据

结果文件:
qtesta.001:解算记录
otesta.001:解算记录的精简版,一般用此文件
htesta.15001:协方差矩阵、参数平均值

postfit nrms 优于0.3时为最佳,大于1时有问题。
sh_setup -yr <yr> -doy <doy> -update <links/files/all> -apr <apr> -series <eop> -localeop  
                        -upd_l -tables <dir> >
  • 3
    点赞
  • 30
    收藏
    觉得还不错? 一键收藏
  • 3
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值