CATS软件在Ubuntu环境下的下载与使用。

        CATS软件，全称为“The Create and Analyze Time Series”，是Simon Williams在NOC（National Oceanography Centre in Liverpool, UK）工作期间写的一个GNSS时间序列分析软件，具体算法可以参考他的论文“CATS: GPS coordinate time series analysis software”。本人（CSND：Eins07）在这里不做原理和算法的过多介绍，主要介绍这个软件的安装与使用。

        CATS软件可以在Linux和MacOS环境下使用（详见官网：CATS - Mike Floyd (MIT)），本人（CSND：Eins07）的环境是Ubuntu18.04，在VM 16 PRO上以虚拟机的方式运行

        首先你应该先做一些准备工作，在终端运行以下代码以保证之后的安装流程

        sudo apt-get install libblas-dev

        sudo apt-get install liblapack-dev

        sudo apt-get install gfortran

        在这里你可能会遇到Ubuntu报错：

                E：无法获得锁 /var/lib/dpkg/lock - open (11: 资源暂时不可用)

                E：无法锁定管理目录（/var/lib/dpkg），是否有其他进程正在占用它？

        不要慌张，你可以稍等一会，玩两分钟手机，然后再运行，或者重启一下虚拟机。如果你不怕麻烦，或者以上两个方法依然无法解决问题的话，你可以通过以下两行代码来解决这个问题（来源：CSDN：小朱不猪）

        “在命令行输入如下命令：
        sudo rm /var/cache/apt/archives/lock
        sudo rm /var/lib/dpkg/lock”

        后续的软件安装和文件修改内容在上文提到的CATS官网中可以看到（Installation一栏），如果你不方便登上官网，或者看英文手册比较吃力，可以参考”hr11asdxi“先生的博客“CATS软件安装[2022.1.20更新]_cats安装_hr11asdfg的博客-CSDN博客”

        我（CSND：Eins07）在这里讲一下，后续的具体软件使用：

        该软件支持对以下噪声模型生成协方差矩阵，从而进行MLE计算：白噪声（WN），幂律噪声（PL），一阶高斯马尔可夫(FOGM) 噪声（相当于自回归AR(1) 噪声），带通噪声（BP），广义高斯马尔可夫噪声（GGM），可变白噪声（即使用形式错误），步变白噪声（即两个时期之间白噪声尺度的变化，必须与至少一个其他噪声模型一起使用），时变白噪声（即白噪声的指数衰减）。

        对于CATS软件，输入的文件格式应为由两部分组成：标题信息和时间序列。标题信息基本上是与时间序列来源的GPS站点相关的参数列表。标题信息以#符号开头。与 CATS 相关的唯一标题信息是偏移列表。该数据集由七列（无特定格式）组成，依次对应于时间（小数年）北、东、上、北误差、东误差、上误差。默认情况下，假定位置以米为单位，但是可以将其作为选项进行更改。如果你的数据没有误差信息，也可以只有前四列。

下图为格式示例

        之后我们使用软件，在终端打开对应文件夹，输入指令：

cats input.neu -A1y1 -Mwh -Ooutput.txt

      即可完成输入文件的白噪声极大似然估计计算，在指令中“input.neu”是你的输入文件，“output.txt”是你的输出文件（名字与后缀可更改，你也可以将名字改为output.wh，方便辨识这是白噪声的分析结果），“-A1y1”是此次计算分析周年与半周年的周期项，“-Mwh”是分析白噪声。如果你想要分析其他噪声，只需要在对应位置更改“-M”开头的指令即可。

        常见噪声的分析指令如下：

白噪声：-Mwh

幂律噪声：-Mpl

闪烁噪声：-Mpl:k-1

随机游走噪声：-Mpl:k-2

一阶高斯马尔可夫噪声：-Mgm

广义高斯马尔可夫噪声：-Mgg

带通噪声：-Mbp

如果你需要分析多种噪声的组合，你只需要将指令并列即可。例如，假如你想要分析闪烁噪声+白噪声的极大似然估计值，你只需要输入指令：

cats input.neu -A1y1 -Mwh -Mpl:k-1 -Ooutput.txt

结果示例如下，其中“+NORT MLE”即为北方向的MLE计算结果

         

        结果会储存在output.txt中。需要注意的是输入输出文件是可以使用多级地址的，例如你想要将结果文件输出在AAA文件夹里，你只要将指令换为：

cats input.neu -A1y1 -Mwh -Mpl:k-1 -OAAA/output.txt

        在统计完MLE值后，我们不能简单选择MLE值较大的模型作为最优噪声模型，需使用Langbein保守估计准则来推测出最优噪声模型。

        Langbein提出了一种保守估计准则，可判断不同噪声模型的优劣，其思路是：①分别计算闪烁噪声+白噪声及白噪声+随机游走噪声组合模型的MLE值，选取MLE值较大的模型作为零假设；②将非整数谱指数幂律噪声+白噪声与闪烁噪声+白噪声+随机游走噪声模型的MLE值分别与零假设作比较，如果MLE差值大于2.6则拒绝零假设，认为该模型更优，否则接受零假设，认为所选的复杂模型无效。若两模型均优于零假设，则选择MLE值较大者作为最优模型；③将带通噪声（Band-pass Noise，BP）+幂律噪声+白噪声与一阶高斯马尔可夫噪声+白噪声+随机游走噪声模型计算得到的MLE值与前面得到的最优值比较，接受带通噪声+幂律噪声+白噪声模型的阈值设为2.6，接受一阶高斯马尔可夫噪声+白噪声+随机游走噪声模型的阈值设为5.2。

​        另外，做一个友情提醒：CATS软件运算速度极慢，所以建议你在处理多个测站的数据时，写个脚本，然后不管它让他24小时跑（需要的可以找我）。我们以一个730历元的数据为例，以下为各个模型的计算时间（单位为秒）：

    ​FN+WN:17

    ​RW+WN:19

    ​FN+RW+WN:155

    ​PL+WN:1030

    ​FOGM+RW+WN:3959

    ​BP+PL+WN:18188

        730个历元，也就两年的数据，在时间序列分析中不算长，如果你要计算长数据（例如超过10年），那么你在计算PL+WN或FOGM+RW+WN这种比较复杂的模型时，请做好以小时甚至天为单位的计算时间的准备，并且不建议计算BP+PL+WN模型。如果你真要进行长序列的分析，这里建议你使用Hector软件。

        以上便是CATS软件的使用与选取，感谢您的阅读，希望我（CSND：Eins07）的内容能对你有所帮助。

参考文献：

[1]:Williams, S.D.P., CATS: GPS coordinate time series analysis software. GPS Solutions, 2008. 12(2): p. 147-153.

[2]:Langbein, J., Noise in GPS displacement measurements from Southern California and southern Nevada. Journal of Geophysical Research, 2008. 113(B5): p. B05405-n/a.