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

博文

tsjumps---时间序列探测阶跃和异常值

已有 457 次阅读 2019-6-28 11:49 |个人分类:时间序列|系统分类:科研笔记| 去粗差、粗差探测

    感慨一下,GAMIT/GLOBK这软件确实强大,内置了许多很实用的fortran脚本,还需进一步去挖掘。

    这不,在分析时间序列过程中又发现一个好用的时间序列去粗差脚本,虽然我不用了,哇哈哈哈哈哈。

    不过这个工具还是不错的,记录一下吧。(tsjumps官方说明附在最后)

(装逼写英文的,哇哈哈哈哈)

  tsjumps input files:org/prt,"value" file created by "ensum" 

  tsjumps output files : a new "value" file remove jumps and outliers which can be used in "multibase" so that "sh_baseline" can direct them easily to the GMT ploting scripts.


How to get the input "value" file ? 

tsjumps Example:



TSJUMPS: Program to detect jumps and outliers in time series files

TSJUMPS: Programs to detect jumps and outliers in time series

Runstring:
% tsjumps <options>

where <options> are any of the following:
 
   -v[alues] file  -- Name of values file from ensum for analysis
   -p[rt/org] file -- Name of print or org file.  Needed to get experiment
                      start and stop times for writing a rename file
   -o[ut] values file -- Name of file to output new values file to (can
                      be then used in ensum and/or multibase.
   -n[ew]          -- sets that only time series that have jumps in them
                      will be written to the new values file.
   -r[ename] file  -- Name of rename file to create or append to.
   -a[ppend rename] -- Append to new renames to the rename file
   -s[cale criteria] -- Sets the values of the jump detection that depend on the
                      postfit chi**2 of fit to the series.  The checks here
                      are of the form value/sigma compared to
                      scale*sqrt(chi**2/f).
                      Three arguments are specified:
                      DRES -- Change between adjacent data points
                      SIGNIFICANCE -- Scale for significance of final jump
                              estimates.
                      Negative values will cause the default to be used.
                      Example: -s 3.0 2.0 would require jumps be
                        3-sigma changes between values or 3-sigma difference
                        between value and average to points before and after
                        (in this case a pair of jumps is added to allow for
                        and outlier. Any jump estimate that is 
                        less than two sigma after all jumps are estimated 
                        would be removed.
    -m[ax absolute criteria] -- Sets the values for absolute criteria (ie.
                      these do not depend on postfit chi**2/f. Three
                      values are given:
                      ABSDRES -- Absolute value of (residual change/sigma)
                      SLOPE   -- Largest slope that will be allowed; if
                                 exceeded slope set to zero.  (mm/yr)
                      SIGMA   -- Largest data sigma that will be checked
                                 when looking for jumps.
    -f[irst] rename -- Sets the first numeric value for renaming sites ie.,
                      sites names become _GPS -> _GNN -> _GNM .. where
                      NN is the value entered and NM=NM+1 and so.  If this
                      option is not given, sites will not be renamed.

    -d[ebug]       -- Turns on debug given detailed outputs of processing.




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

上一篇:全球igs站分布图
下一篇:ubuntu使用apt、apt-get无法安装和更新软件解决办法

0

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

数据加载中...

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

GMT+8, 2019-8-25 06:04

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部