ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

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

2021-06-16 19:31:30  阅读:206  来源: 互联网

标签:文件 lfile otl 001 sh 0.050 数据处理 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 -doy -update <links/files/all> -apr -series -localeop
-upd_l -tables

>

标签:文件,lfile,otl,001,sh,0.050,数据处理,Gamit
来源: https://blog.csdn.net/qq_40417070/article/details/117965052

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有