|||
处理过程中,为了自己的方便使用,需要编写自己的脚本(用红色标出)
第一种方法是,得到站点位置后:
1、 grep Unc. sina.org >sina.apr 得到运算点的坐标
作为apr坐标,然后输入到globkcomb的cmd的apr选项中,然后重新进行画图
2、 然后运行在gsoln目录下:
glred6 glred_20150811.prt glred_20150811.log 201407_NSFBay.gdl cmd
sh_plotcrd -f *.org-expt emed00
上述仍然需要编辑时间序列,遇到不好的天数进行重新编辑,所以我们采用:
第二种方法,当然我们介绍一种新的globk处理时间序列的算法tsfit
1、在globk的cmd中org_opt 加入PBOB 命令变为,运行
glred 6 glred_20150811.prt glred_20150811.log 201407_NSFBay.gdl
2、得到pos文件
tssum . ex –R *.org
. 是当前目录,ex是随意命名,然后得到pos文件
3、编辑tsfitNaNd文件,运行tsfit
tsfit tsfitNaNd tsfit.sum *pos
命令编辑如下:
eq_file ../tables/RESET.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、将得到的eq文件 上面输出的eq文件应为tsfit.rns 输入到globkcomd的cmd文件中,重新运行glred,删除原先的org文件以及prt文件
glred6 glred.prt glred.log 201407_NSFBay.gdl globkcmd
5、得到新的org文件,进行画图
sh_plotcrd -f *.org-expt emed00
可以得到不同的图以及误差重命名的GPS站点图
也可以用
sh_plot_pos -f globk_rep.org 生成raw时间序列
备注:当然在精准的时间序列时,可以选择利用globk得到的站点坐标返回gamit进行重新再次处理,然后继续返回这样的方法进行求解时间序列可以得到更佳精确的时间序列。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 01:30
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社