GLOBK平差+绘制时序记录
看了陈超老师的课,结合自己的处理过程,记录如下:
准备工作
相较于gamit的准备工作,globk的平差准备工作要简单很多。
- 在工程目录(demo)下:
sh_glred -cmd
将tables中的globk.cmd和glred.cmd复制到demo下的gsoln中,当然也可以手动copy
- 在demo中新建文件夹hfile,在demo中查找所有的h文件,具体命令如下:
mkdir hfile *#新建hfile文件夹*
find ./ -name "hdemo*"
- 把所查找到的h文件复制到新建的hfile下:
find ./ -name "hdemo*" | xargs -i cp {} ./hfile/
- 为了防止后面出现错误,先进行:
cp ./tables/svnav.dat ./hfile/
rm glbf/*
- 将h文件转换为二进制文件
htoglb ./glbf/ ./tables/svs_demo.svs ./hfile/h*
- 制作h文件列表(gsoln下)
ls ../glbf/h*.glx > demo.gdl
平差
做重复性
glred 6 globk_rep.prt globk_rep.log demo.gdl globk.cmd
做点位平差
globk 6 globk_demo.prt globk_demo.log demo.gdl globk.cmd
绘制时间序列
要现在linux中安装好GMT
sh_glred -expt demo -s [start year] [start day] [stop year] [stop day] -opt H G E