参考了http://gnss.help/2016/09/07/GLOBK-StepbyStep/和http://wap.sciencenet.cn/blog-3391834-1140957.html两篇博文,感谢前辈的无私奉献 Thanks♪(・ω・)ノ
globk各个文件的介绍参见globk的说明书和以上两篇博文。
在解算了22个CORS+5个IGS站三年的数据后,得到单日解,现准备利用glred命令将前三年的单日解h文件进行合并。
1.先将第一年,即2013年的各个单日解文件夹内的h文件提取出来,在工程文件夹ah13下新建文件夹hfile,输入命令如下:
find ./ -name ‘hah13*’ | xargs -i cp {} ./hfile
2.将h文件转化成GLOBK能够识别的二进制文件,在工程文件夹ah13下输入命令:
htoglb ./glbf ./tables/svs_demo.svs ./hfile/h*
3.仍然在ah13文件夹下输入命令:
sh_glred -cmd
调取globk_comb.cmd、glorg_comb.cmd两个命令文件至gsoln文件夹下(该文件夹之前sh_gamit批处理时已经产生)
4.进入gsoln文件夹下,将需要合并的h文件制作成文件列表,为下一步的合并做准备。由于此次将2013年全年的h文件全部合并,所以将所有文件都制作成文件列表。在gsoln文件下输入命令:
ls …/glbf/*.glx > demo_glx.gdl
5.修改gsoln文件夹下globk_comb.cmd文件:
原始设置:
修改成
前两项是注释掉调用glorg命令的,在合并多日解文件中不需要调用glorg。(globk_comb.cmd和glorg_comb.cmd命令文件中,命令行首列空格表示执行命令,首列非空表示此行命令没有被执行)
将原来的注释行首列删除,使命令得以执行。并将合并后的新文件命名为:H2013_COMB.GLX
6.在gsoln文件夹下执行命令:
glred 6 glred.org glred.log demo_glx.gdl globk_comb.cmd
7.按照上述合并2013年h文件的方法,分别在制作2014和2015年的h文件,得到三个合并后的GLX文件:
H2013_COMB.GLX、H2014_COMB.GLX、H2015_COMB.GLX