||
得到的原始时间序列中难免会有粗差,去粗差的工具和方法有很多。
我这里使用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文件重新解算时间序列就可以得到去除粗差的序列了。
去除粗差前后序列
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 06:32
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社