gamit10.71解算教程（入门级别，仅供参考）

数据准备
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，第二次网平差需要给定义的基准点紧约束。

这里采用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查看使用的基准点等。

gamit10.71解算教程（入门级别，仅供参考）相关推荐

解决GAMIT10.71解算北斗二号三号观测值问题
发现问题之前在用GAMIT10.71处理北斗数据时发现,利用sh_gamit命令处理的同一批站点的观测文件,全程无报错,但最后结果有些站的RMS为0,即没参与计算. 查看autcln.post.su ...
菜鸟零基础建站入门指引（仅供参考）
菜鸟零基础建站入门指引(仅供参考) 自主建站需要有域名.空间和网站备案综合说明:参考https://help.aliyun.com/knowledge_detail/39784.html 大致流程: ...
最新版WinRAR5.61去广告代码教程分享（仅供学习交流）
最新版WinRAR5.61去广告代码教程分享(仅供学习交流) 第一步:到WinRAR官网www.rarlab.com下载自己需要的版本,选择Chinese Simplified 64bit 安装即可. ...
最新版WinRAR5.70去广告代码教程分享（仅供学习交流）
过段时间会失效,每次都弹框提示40天后购买,替换原绿色版中的rarreg.key即可使用了. 最新版WinRAR5.61去广告代码教程分享(仅供学习交流) 第一步:到WinRAR官网www.rarla ...
CSS入门。（仅供参考）
文档出处部分Coogle和<原创>,如有雷同纯属巧合 <仅供参考> 首先说明以下内容是本人自己整理(仅供参考!!!!!!!!) 首先浏览器我建议使用Google Chrome ...
第十四届蓝桥杯校内模拟赛第二期-Java个人题解（仅供参考）
刚刚结束第十四届蓝桥杯校内模拟赛第二期,在这记录下自己的代码不保证正确! 不保证正确! 不保证正确! 有问题一起改正!! 题解 2048 代码: package _14届模拟2;public cla ...
【学生请假系统——开题报告分享（仅供参考呀）】
目录一.研究目的和意义: 二.课题研究现状: 三.课题研究主要内容.实施方案及创新点: 四.课题进度安排: 五.主要参考文献:(仅供参考呀!) 一.研究目的和意义: 1.研究目的开发高校学生请假系 ...
2019第十届蓝桥杯C++B组题解（赛后重写的，不确保答案正确性，仅供参考）
先说一下这次的感受吧,我们考场比较乱,开始比赛了,还有的电脑有故障,(向这些人发出同情),第一次认真参加比赛,真正比赛的时候感觉没有那么正式,很乱,各种小问题,(例如博主就没找到题目在哪里,找到后又不 ...
项目线上问题应急处理规范（不同公司流程规范不一样，仅供参考）
前言:整理了一下项目线上问题应急处理规范,仅供参考! 1.1线上报障流程 1.2.目的为了明确线上事件的报告.应急响应和处理机制,确保能以快速.高效.准确的应急处理能力来保证各类系统的正常运转,当发 ...

gamit10.71解算教程（入门级别，仅供参考）

gamit10.71解算教程（入门级别，仅供参考）相关推荐

最新文章

热门文章