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

博文

时间序列去粗差

已有 511 次阅读 2019-4-27 16:48 |个人分类:时间序列|系统分类:科研笔记

    得到的原始时间序列中难免会有粗差,去粗差的工具和方法有很多。

    我这里使用globk自带的tsfit;tsfit是一个时间序列拟合工具,可以自己设置周年、半周年等参数。。。。

    关于tsfit的详细解释参考说明书。

tsfit去粗差的步骤:

1.首先要使用tssum从org文件中提取PBO格式的pos时间序列文件

2.准备tsfit的控制文件tsfit.cmd,我这里只是为了去除粗差,所以不设置其他参数。

        设置大于5倍中误差的坐标予以剔除

*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

3. 运行tsfit后会得到粗差文件tsfit.rns

  tsfit tsfit.cmd <sum_file> <PBO_files>

4.将tsfit.rns作为eq文件重新解算时间序列就可以得到去除粗差的序列了。



去除粗差前后序列




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

上一篇:时间序列高程方向断崖式变化
下一篇:时间序列PCA分析记录

0

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

数据加载中...
扫一扫,分享此博文

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

GMT+8, 2019-6-21 00:18

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部