下面的实例中,处理的数据为2020年第5天的数据,4个远距离的IGS站(chum、cusv、daej、irkj)、解算其他3个站点。
gamit10.71下载链接(https://cloud.189.cn/t/yM3MVfVryENf (访问码:5l39))
1、数据处理前更新表文件
更新方法在之前的博文里面gamit10.71更新tables的几个方法(使用FileZilla和Ubuntu的gftp软件)-CSDN博客。
表文件的网址如下
ftp://garner.ucsd.edu/pub/gamit/tables/
ftp://everest.mit.edu//pub/GRIDS/
每次用更新(可能有的并不用每次吧)
pole.usno 极移表
ut1.usno 地球自转参数表
svs_exclude.info 坏卫星信息文件
svnav.dat 卫星列表
rcvant.dat 接收机天线对照表
antmod.dat 天线相位中心改正正表
dcb.dat 伪距改正参数
每年更新一次(按对应年份使用)
luntab.2020.J2000 月亮表
soltab.2020.J2000 太阳表
nutabl.2020 章动表
vmf1grd.2020 映射函数模型
leap.sec 跳秒
atmdisp_cm.2020
2、文件准备工作
2.1、创建工作文件夹
在“opt”文件夹内新建“work”文件夹。
2.2、在work内新建数据文件夹
新建“brdc、igs、rinex、met”4个文件夹,分别存放当天的广播星历、精密星历和观测值文件、气象文件。分别把对应的文件放在文件夹中(可不用这个met文件,本次演示没用。广播星历和精密星历的下载方法下次出一期)。
(注意:如果文件采用.z压缩格式,则应当使用进行解压;如果解压出的观测值文件仍采用.d的压缩格式,则应使用crx2rnx软件将其转化为.o,如果不会就百度d转o吧)。
2.3、在work文件夹下右键,点击“在终端中打开”,输入“sh_setup -yr 2020”命令,会发现主文件夹下出现了“tables”文件夹。这个终端别退,后面要用呢。
2.4、打开上一步的“tables”文件夹,右键在终端中打开, 输入“ls -sl”命令查看链接情况(这里有的用不到,红色也无所谓的。比如第一个红色,我在第6步中sestbl.设置了atl.list为no就是用不到的)。
2.5、红色连接处理
假如atmdisp_cm.2020没连接上,是红色的(这里是假如,如果你是红绿色盲,我建议立即放弃)。打开第2.3步的终端(上面说了要用的啊),依次输入下面三个命令。
% 移除红色错误链接
rm -f tables/atml.grid
% 链接正确文件,这个文件是第1步更新的表文件
ln -s ~/gg/tables/atmdisp_cm.2020* tables/atml.grid
% 这个是用来检查是否连接成功,在下图为绿色说明成功
ls -l tables/atml.grid
3、更改配置文件
3.1、更改sestbl.文件。
打开tables下的“sestbl.”,该文件内部有详细的说明,翻译下就知道大概意思。我的修改如下,所在行数已经标出来,如果好兄弟还是找不到,我劝你放弃读研(我主要是为了获取ZTD)。
36行 Elevation Cutoff = 10
33行 Interval zen = 1
44行:Met obs source = GPT 50
46行:Output met = Y
49行:DMap = VMF1
50行:WMap = VMF1
52行:Use map.grid = Y
64行:Use otl.grid = Y
71行:Use atml.grid = Y
3.2、更改sites.defaults文件。
打开“tables”文件夹下的“sites.defaults”文件,它是用来控制需要参与解算的测站。把红圈处改为工作文件夹(我的是“work”),蓝圈处改为测站名字。
3.3、生成测站信息文件station.info文件。
3.3.1、将work/tables下的“station.info”文件拷贝到“rinex”文件夹下,打开并编辑,仅保留以#或*开头的前几行,如下图红框的内容。
3.3.2、删除完成,如下图一样。
3.3.3、右键进入终端,打开终端并进入“rinex”文件夹,运行下面命令(20是年份的意思):
sh_upd_stnfo -files *.20o
3.3.4、命令完成后如下图,将生成的“station.info”文件拷贝到text/tables文件夹,覆盖原文件。
注:输入此命令有时会报错:打不开guess-rcvant.dat和rcvant.dat文件,这时只要从tables把这两个文件复制到当前目录下即可。
3.4、测站近似坐标文件lfile.文件。
igs站可不用生成lfile.文件。进入“work/rinex“文件夹内打开终端。解算的每个站点都要依次生成,我结算三个站点,就需要输入三个命令。
sh_rx2apr -site dxin0050.20o -nav ../brdc/brdc0050.20n
sh_rx2apr -site fjwy0050.20o -nav ../brdc/brdc0050.20n
sh_rx2apr -site gsdh0050.20o -nav ../brdc/brdc0050.20n
注:以dxin站为例sh_rx2apr -site dxin0050.20o -nav ../brdc/brdc0050.20n中红色为o文件名字,蓝色为brdc名字。(igs站可以用这个生成lfile.的文件grep bjfs itrf14.apr > lfile.1 )。
将上面的站点生成“*.apr(或者*.lfile)”文件进行合并数据进行合并,输入下面命令:
cat *.apr > lfile.
cat lfile.* > lfile.
会组合生成“lfile.”文件,将其拷贝到text/tables文件夹,覆盖原文件。
4、解算文件
4.1、解算命令
在前面章节2.3的终端下输入“sh_gamit -expt work -d 2020 005 -orbit IGSF -met -metutil Z -noftp”等待解算完成。(work:文件夹。-d:处理一天数据,多天改为-s。2020 055:处理的年份和年纪日。-noftp:不启用ftp。还有几个忘记了,看到了再补充)。
4.2、结果文件
将存放在名称为年积日的文件夹内,此例中为“005”,可供参考的结果文件为:
sh_gamit_2020_005.summary 解算总结
qworka.005 解算记录
oworka.005 解算记录的简略版
4.3、分析结果文件
打开“qworka.005”,划到最底部。nrms是衡量结果的重要指标,通常postfit nrms值在3~10mm,prefit nrms值应在0.2左右,如果存在很大偏差,就需要查找原因了,重新解算了。
5、个人总结
我个人做gamit解算是为了提取天顶总延迟 ZTD。我是大气水汽方向的,所以需要求ZTD还有可将水量(里面的每个站点m文件就是可以提取需要的ZTD和PWV)。教程或许会有一些瑕疵,还请海涵。如果有不对的地方,欢迎大家指正。本文出于帮助初学者利用gamit软件解算GPS基线的目的撰写,但受笔者水平限制,难免存在疏漏和错误,欢迎大家批评指正。