- 数据准备
1.1 数据下载
本次实例拟使用aira bjnm chan cmum hkws mizu twtf ulab等IGS站数据进行解算,计算日本的两个站mizu, aira的位置,解算时间为2019年129天。
在主目录下创建工程文件目录
mkdir test
链接tables文件夹下中的文件
sh_setup -yr 2019 -apr igs14_comb.apr
(注:igs14_comb.apr文件中保存IGS站的先验信息,采用itrf14_comb.apr也可以,两个文件有细小的差别主要是地球参考框架的实现方式不同)
进入test文件夹下创建精密星历igs 广播星历brdc 观测值rinex文件夹 电离层格网文件夹ionex
(注:电离层格网文件是gamit10.71计算新加入的必须的文件)
cd test
mkdir igs brdc rinex ionex
首先下载rinex观测值数据,实例中均为IGS站,观测值可以从武汉大学IGS数据中心下载,将下载好的文件放入rinex文件夹下(注:计算前要确保下载的文件中接收机、天线类型、天线高、文件名等是否正确)。然后准备广播星历数据,同样可以在武汉大学IGS数据中心下载,下载后放入brdc文件下(注:此次示例仅采用GPS计算,使用brdc广播星历即可,若要计算多系统的需要在cddis下载brdm多系统广播星历)。下载精密星历放在igs文件下(网解仅采用精密轨道文件即可,无需精密钟差)。电离层格网文件需要从cddis中心下载,文件名为igsg1290.19i
1.2 表文件更新(注:表文件对于解算精度很重要,这里仅展示如何计算不考虑精度,因此省略)
2. 修改表文件
2.1 sites.defaults
在tables文件夹下找到并打开该文件
cd tables
vi sites.defaults
修改工程名(test) 与测站名(需要解算的测站)
2.2 station.info
vi station.info
首先删除文件中的多余信息
:7,$d (注:7代表开始行,美元符号代表结尾行)
键入回车即可删除
保存后退出
输入命令编辑station.info文件
sh_upd_stnfo -files . ./rinex/*.*o
生成的station.info文件如下:
(注:该文件原本存在IGS站信息,重新生成主要是为了非IGS站)
2.3 lfile.文件
该文件中主要保存各个站的初始坐标,示例中均为IGS基准站,不作更新即可,直接采用igs14_comb.apr文件
2.4 sittbl.文件
起算点坐标分量约束为5cm, 非起算点为30m(本示例中均为IGS站,不改也可以)
2.5 sestbl. 文件
建议仔细学习一下其中每行的意思,这里不做阐述,采用默认值
如果想采用各种模型,如海潮模型,大气模型等可以在MIT的ftp服务器上下载后放到tables中
ftp://everest.mit.edu/pub/GRIDS/
2.6 process.defaults
这个文件主要是用于控制计算的处理进程,大部分配置不需要更改,可以自己研究一下。可以在这个文件里面设置计算时常,开始时间等。
3.批处理
sh_gamit -d 2019 129 -orbit IGSF -expt test -noftp
运行报错:
这个问题是因为配置文件默认改正海潮模型,需要otl.list和otl.grid文件,这两个文件需要手动下载,感兴趣的话可以自己下载一下这两个文件放入tables文件下即可。作为演示,我们这里更改配置文件不考虑海潮模型
进入tables文件下
vi sestbl.
修改如下配置:
重新运行sh_gamit命令即可。
基线计算结束
- 计算结果
解算完成后需要检查计算结果。计算完成会生成以年纪日为目录名称的文件夹,这里是129,在该文件夹下有一个sh_gamit_ddd.summary文件统计计算结果。其中postfit RMS值在3~10mm,prefit nrms值应在0.2左右。还有宽巷和窄巷模糊度固定率。本示例中的站有在日本的,日本的测站由于板块活动明显数据质量较差,同时为考虑海潮模型,对流层模型改正,部分测站如CMUM站较差,模糊度固定率也偏低。
具体基线解算结果在q文件和o文件中,o文件是简略版。
5 GLOBK平差
在test目录下执行
htoglb glbf /dev/null 129/h* a.*
进入gslon目录下,在安装目录~/gg/tables/文件夹下拷贝globk.cmd和glorg.cmd文件到该目录下,打开globk.cmd文件设置平差的起算点,修改apr_site如下。网平差需要进行两次,第一行apr_site all 10 10 10 0 0 0用于第一次网平差,需要给松约束10 10 10,第二次网平差需要给定义的基准点紧约束。
文章来源:https://www.toymoban.com/news/detail-824605.html
这里采用IGS站的坐标作为起算点,计算日本的两个站的坐标,可以给较强的约束。起算点个数看实际情况给出。
修改glorg.cmd文件中的stab_site行添加固定点。固定的个数同样看实际情况几个都可以。
计算多天的数据可以在该文件中固定点约束前面行加上stab_site clear命令
在gsoln目录下执行 ls . ./glbf/*.glx > test.gdl
最后执行 globk 6 globk.org globk.log test.gdl globk.cmd
在globk.org文件中即可看到最终的结果,apr文件是最终的平差结果,位置估计部分可以看到基准点使用个数和调整结果。也可以在在终端使用命令 grep ‘POS S’ globk.org查看使用的基准点等。
文章来源地址https://www.toymoban.com/news/detail-824605.html
到了这里,关于gamit10.71解算教程(入门级别,仅供参考)的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!