|||
处理过程中,为了自己的方便使用,需要编写自己的脚本(用红色标出)
在计算速度时候,我们可以直接利用我们的gps测站求解速度,但是速度具有参考性,我们得到的速度并不是在文献中常说的相对于ITRF框架或者欧亚框架的速度,所以我们此次介绍得到ITRF框架下的速度
在检查完自己的时间序列后,删除粗差后,我们来求解速度
1、首先我们得到全球解,也就是全球的h文件,
下载全球数据利用sh_get_hfiles得到mit的H文件或者sopac的h文件
· sh_get_hfiles-yr 2013 -doy 246 -ndays 5 -net igsall eura
· 或者直接利用sh_glred进行下载
· sh_glred-s 2013 024 2013 025 -expt sina -opt R LA H G T >sh_gred.log 合并全球解与自己的解单天解
这样其实就是相当于你用gamit结算了全球的解,然后一起得到一个参考全球的速度,为了将自己的解与全球的解融合,这里在gamit中结算需要有几个igs站点的结算。
2、将tables的globkvel的cmd与glorg的cmd文件拷贝过来,编辑两个vel的cmd用新的apr和新的eq文件,在org_opt中加入VSUM 使用松约束加入转换旋转对地球eop参数。在glorg中加入rate_org的估计参考站为处理的几个站,globk结算速度时,gdl后面不用+1,默认为结算速度,
3、运行globk的命令
globk 6 sina.prt sina.log sina.gdl globk_velNaNd
4、 得到chi~比例因子 加入到gdl中,删除前面生成的prt文件与org文件重新运行
· sed‘s/$/0.034 /g' sina.gdl
· globk6 sina.prt sina.log sina.gdl globk_velNaNd
5、画出速度
sh_plotvel -fcoco1213_rep.vel -maxsigma 4 -R225/310/0/50 -page L -ps coco1213_rep
注意:
查看globk_vel.org中确信所有的chi2增量在范围0.2-1.5中,如果一些偏高,则应该查找时间序列的outlier,然后分别对每一个站点进行sig_neu,或者对于每一个h文件因子在gdl中进行调整因子,最后查看总的chi增量在chi2 increments自己的解的因子用1左右,并且参考框架的站点的rms小于2mm/yr
当然:我们也可以利用tsfit进行结算速度,这也是有效的方法之一
1、 仅仅准备globkcomb的cmd与glorg的cmd就可以,然后编辑文件,加入PBOB的命令在org_prt中,然后首先运行
glred6 sina.prt sina.log sina.gdl globk_velNaNd
2、得到pos文件
tssum . ex –R *.org
. 是当前目录,ex是随意命名,然后得到pos文件
3、编辑tsfitNaNd文件,运行tsfit
tsfit tsfitNaNd tsfit.sum *pos
命令编辑如下:
eq_file ../tables/ITRF14.eq
*Annual signal
per 365.25
*Remove bad outliers (5-sigma)
nsigma 5
* Savethe velocity file (and an apriori coordinate file)
velfile @.vel
* Savethe edited data list
rep_edits @.rns
* Usethe realistic sigma algorithm for the velocity sigmas
real_sigma
4、直接可以得到velfile,然后直接画图即可。实现的算法就是使用的real_sigma,也就是利用tsview以及velview的工具实现的算法,这里也就是tsfit是它们的最新替代品。
下一节会讲到 如何实现参考框架,以及转化到欧亚板块的速度。
以及如何定义自己的板块,得到自己板块区域的速度
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 21:47
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社