1.数据准备
参考博文:http://gnss.help/2016/09/01/sh-gamit/index.html
2.表文件链接
工程目录下链接命令:
sh_setup -yr <year> -doy <doy>
例如:sh_setup -yr 2017 -doy 130
3.格网文件链接检查
表文件链接后工程目录会生成tables文件夹,如果使用海潮格网改正文件输入以下命令检查格网文件是否链接成功:
ls -l tables/otl.grid
出现以上表链接成功
3.配置文件修改
①tables的sestbl. 文件中默认:
Use otl.list = N ; Ocean tidal loading list file from OSO
Use otl.grid = Y ; Ocean tidal loading grid file, GAMIT-format converted from OSO
如果已经有了海潮格网文件,可不用更改,没有的话可到ftp://garner.ucsd.edu/pub/gamit/tables/中下载(如图),放到安装目录下的tables重新在工程目录链接即可;
如果是igs测站,改为
Use otl.list = Y ; Ocean tidal loading list file from OSO
Use otl.grid = N ; Ocean tidal loading grid file, GAMIT-format converted from OSO
表示使用list不使用girb文件,其中list中是从grid中提取出来的igs的海潮改正值,这是没有grid文件也没有影响
②如果设置了
Output met = Y
那么你还可以得到包含对流层参数解算结果的 Z-文件。可以由这个文件进一步处理得到可降水汽 PWV 的值
3.初始文件生成
①初始文件为station.info 和 lfile文件,station.info文件必须制作,制作方法参见博文:http://gnss.help/2017/03/05/create-station-info/index.html,如果是igs测站tables目录里一般会已经含有station.info文件,如已包含不用制作,如果没有可到ftp://garner.ucsd.edu/pub/gamit/setup/下载(如下图)
②Ifile文件批处理时Gamit会自动根据导航电文自动生成
4.基线解算
基线解算命令格式有两种:
sh_gamit -expt <expt> -d <year> <doys> -orbit <IGST>和sh_gamit -expt <expt> -s <year> <d1> <d2> -orbit <IGST>
前者适用于少量数据,后者用于大量观测日期上连续的数据需要处理,<d1> <d2>分别为起止年积日。
例如要处理 2016 年前 3 天的观测数据,你可以这样输入:
sh_gamit -expt demo -d 2016 001 002 003 -orbit IGSF # IGSF 表示使用最终星历
或sh_gamit -expt demo -s 2016 001 003 -orbit IGSF
5.批处理解算结果
q文件 解算记录
o文件 解算记录的简略版,一般关注此文件,文件中postfit_nrms项优于0.3左右时最佳;如果大于1.0,则表示此解存在问题。
h文件 协防差矩阵、参数平差值
6.平差
工程目录下运行命令:
sh_glred -expt <expt> -d <year> <doys> -opt H G E或sh_glred -expt <expt> -s <year> <d1> <d2> -opt H G E
例如:sh_glred -d 2017 130 -expt demo -opt H G E
平差成功会生成gslon文件夹,文件夹里含一个org文件(包含解算点的三维坐标和相关参数)和psbase开头的文件(测站三维坐标时间序列的图),如果仅是单天的文件不会绘制图形。