博文
时间序列去粗差
||
得到的原始时间序列中难免会有粗差，去粗差的工具和方法有很多。
我这里使用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文件重新解算时间序列就可以得到去除粗差的序列了。
去除粗差前后序列
https://blog.sciencenet.cn/blog-3391834-1175812.html
上一篇：时间序列高程方向断崖式变化
下一篇：时间序列PCA分析记录
