Gentle''的個人博客分享 http://www.ueservicedoffices.com/u/tao164411096 九分耕耘,一份收獲。

博文

時間序列去粗差

已有 533 次閱讀 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://www.ueservicedoffices.com/blog-3391834-1175812.html

上一篇:時間序列高程方向斷崖式變化
下一篇:時間序列PCA分析記錄

0

該博文允許注冊用戶評論 請點擊登錄 評論 (0 個評論)

數據加載中...
掃一掃,分享此博文

Archiver|手機版|科學網 ( 京ICP備14006957 )

GMT+8, 2019-6-27 20:06

Powered by ScienceNet.cn

Copyright © 2007- 中國科學報社

返回頂部
时时彩平台