Gentle''的个人博客分享 http://blog.sciencenet.cn/u/tao164411096 九分耕耘,一份收获。

博文

[转载]GLOBK求解时间序列精度评定并绘制时序图

已有 214 次阅读 2018-11-19 13:56 |个人分类:GAMIT/GLOBK|系统分类:科研笔记|文章来源:转载

globk执行glred命令,绘制时间序列图。

    glred 6 globk_fjxm.prt globk_fjxm.log fjxm.gdl globk_comb.cmd
    sh_plotcrd –f globk_fjxm.org –e fjxm_globk –s long –c 1 –o 1 –b site.list

    site.list是点名列表,每一行首为空格。


在时间序列曲线上除了标注估计的线性速率及其误差外(下图头部的rate( mm/yr)标识线性速度及误差),还有如下表征时间序列精度的指标:WRMS和NRMS

时间序列精度的指标:

1)wrms( weighted root mean square) 加权均方根误差


2)nrms( normalized root mean square) 归一化均方根误差,无纲量;


wrms越小表示时间序列的重复性越好,(平面1-2 mm ,垂直:3-10mm)
站点序列中散点的分布近似为高斯分布,nrms ~1.0



如上图,DZAX时间序列中出现异常值,可以在globk_comb.cmd中的地震文件选项中进行设置(eq_file ~/gg/tables/itrf08_comb.eq)
即,现在org文件中找到异常值发生的准确时间,然后在~/gg/tables/itrf08_comb.eq中按照如下方法剔除该点:
    rename came_gps came_xcl h0112171200_adri.glx 01 12 17 0 0 01 12 17 24 0
   (rename sitename h-file name start time end time)
设置完成之后,在重新执行glred命令,重新绘制时间序列图;

-xcl:表示从文件中删除




http://blog.sciencenet.cn/blog-3391834-1147061.html

上一篇:GAMIT10.7分布处理流程
下一篇:[转载]globk通用讲解(很实用)

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2019-10-23 18:10

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部