GPS时间序列分析---剔除跳跃点，拟合时间序列

通常利用GPS时间序列进行数据分析时，会遇到大地震的发生，这个时候会导致GPS的观测结果出现很大的跳跃值，这对后续的数据处理和分析带来了困难。这里分享一个最近了解的，可以用于处理这一问题的工具包---TSAnalyzer。下面主要介绍该工具包的下载、安装和简单的数据处理方法。 

（1）GNSS数据的获取 

GNSS数据下载至内华达大地测量实验室 (NGL) (http://geodesy.unr.edu)，该网站提供了全球范围内的GNSS时间序列。漫游地图，找到我们需要的站点的数据，这里我以CHUM站点为例，进入数据下载的界面。点击界面的tenv3或者tenv格式的数据，进入数据页面，右键选择另存，下载数据成功。 

NGL提供的GNSS数据

（2）TSAnalyzer软件的获取 

今天（2023-3-27）在谷歌冲浪的时候，无意之间发现这个很多GNSS数据处理程序的网站：

GPS_TOOLBOX
里面有各种论文发表后开源的程序包，包括matlab、C++、python的工具，大家可以下载学习。这里我们找到

TSAnalyzer, a GNSS Time Series Analysis Software -- by WU Dingcheng, YAN Haoming, SHEN Yingchun
网站提供了linux版本和用户友好的window版本，此处我本人下载的是window版本。下载对应的软件，点击运行 

（3）GPS数据的处理

直接点击运行，可以进入以下的界面。这个软件比较友好的地方是自带了示例的数据，我们点击Menu，进入Example文件夹，选择chum.neu文件，打开文件。在主界面的左上方，点击Add File中的文件，则可以显示文件加载的结果。

根据网站数据的信息，CHUM站点记录到了在2007年发生了一次阶跃，这表明是地震活动导致的。因此我们需要消除这个数据的跳跃。 

接下来，首先去除离群值，点击右面板的IOR factor，点击clean。然后选择地震信号中断的功能。标记阶跃值【本人是手动标记的，应该有更精确的输入标记的方式，后续得多参照使用说明书】 

点击detrend，选择continuous，得到去阶跃值的结果，然后选择数据导出Save，即可得到去除阶跃值的数据结果。

期待你的点赞收藏