[QOCA学习笔记]利用QOCA软件analyze_tseri模块进行时间序列分析

操作环境：win10+虚拟机ubuntu下QOCA（ver 1.33）
数据：通过解算得到的安徽省CORS参考站坐标时间序列（neu格式）

对于QOCA软件的介绍，在其官网https://qoca.jpl.nasa.gov/上可以查看相关的描述。其analyze_tseri模块可以用来分析GNSS坐标时间序列。对该模块的详细介绍可以在这里下载：analyze_tseri analysis.pdf

与之前pca模块相似，在调用analyze_tseri模块之前亦需要仔细查阅其驱动文件的设置。下面简单介绍一次简单的时间序列分析过程。

1. QOCA的各模块调用均需要事先准备好驱动文件，在其中做好各种参数的设置工作，在运行时只需要在终端输入：analyze_tseri analyze_tseri.drv 即可完成分析时间序列。以驱动文件analyze_tseri.drv为例，按照介绍文件的描述对其进行标注：

=========================================================================
*          << key word controlled driving file format >>                *
* symbol ":" must exist in command lines as index pointor               *
* any non-blank character at first column means comment line            *
* 第一列存在任何非空的的字符，均认为该行为注释行，不执行                   *                             
* empty after ":" means comment line too, but line apears in out file   *
=========================================================================
---------- part 1 -- general i/o informations
  apriori value file:                qoca.apr   
%（站点位置和速度先验值文件）与之前QOCA各模块的先验文件格式一样
%注：对于GAMIT/GLOBK用户，要想在利用QOCA解算前得到这样的先验文件，可以通过mk_net模块得到，本文后面会介绍
  input list file:                   input.list
%QOCA独立模式输出映射文件的列表。第一行有两个输入项。第一项是文件数。第二个项是文件格式类型。
%type=1表示QOCA map格式，type=2表示SOPAC neu格式。其他更多输入文件类型请查阅介绍文件
  sit_list file :                    analyze_tseri.site
%（站点列表文件名）在站点列表文件中，第一行是站点总数。其余的行是站点名称。
  est_parameter file :               analyze_tseri.para
%（参数估计文件）核心文件。第一行为台站总数，下面每个台站第一行台站名和估计参数的数目，接下来是每个估计参数的具体要求。
%每个参数行的第一个数为索引指数，第二个数为参数先验值，第三个值为参数的约束值，第四、五个数为参数的有效区间的起始和结束时间
  output file:                       analyze_tseri.out1
%输出文件名称
  residual file :                    analyze_tseri.resi
%输出的第一个残差文件命名
  res_option :                       4127
%控制第一个残差文件的输出选项
c  omitted_span file:              omit_span.list
%需要删除历元
c  cme_correction file :           pca.seign pca.cpt
%经过pca模块解算后得到的空间向量文件以及分解出来各分量的文件，打开这个选项相当于对滤波后的时间序列进行参数估计
c  specific term_out file :           season_out.term
c  specific term_option :            8160
  resi_file2:                     analyze_tseri.resi2
%第二个残差文件的名称
   res2_option:                      31 4115 4113
 %第二个残差文件的输出选项，本次最多可以输出三个不同选项的残差文件
c  res2_option:                      4113 4115 4099
c  mdl2_out:                       analyze_tseri.mdl2
c  mdl2_option:                      4078 4076 4092
  enu_correlation usage :           yes 
  jump_day removal :               yes
c  cutoff criterion (year):              2.5
  oneline_netformat:                 yes
c  color_noise analysis model :          2
c  atmos_load correction file :      atmos_load.list
c  moisture_load correction file :   mois_load.list
c  snow_load correction file :       snow_load.list
c  ocean_load correction file :       ocean_load.list
  span to est jump apr (est_jump_span, year) :   0.1 0.1
  weak_obs (big sigma) criteria mm:     60.0 60.0 80.0
  outlier (big o-c) criteria mm :     100.0 100.0 200.0
  very bad_obs criteria mm :     1000.0 1000.0 1000.0
  t_interval :                2013.0 2018.5
  end: 
---------- part 2 -- apriori informations
  exit:

2. 介绍完驱动文件，其中有几个选项比较重要：
   （1）参数设置文件样例如下：
   
   不同的索引对应着不同的参数，analyze_tseri可以设置的参数有57个，具体介绍可以下载相关介绍文件仔细查看：QOCA指南.zip
   
   （2）残差控制文件：
   
   
   （3）先验文件的制作：
   用gamit/globk的apr文件，通过mk_net将其转换成QOCA需要用到的apr先验文件即可：在gamit解文件夹/rinex/下
   cat *.apr > globk.apr
   mk_net globk.apr qoca.apr globk
   即可生成名为qoca.apr的qoca需要用的先验坐标文件。
   在执行analyze_tseri操作后，生成的analyze_tseri.out文件中也包含了更新后的台站坐标解，而且是以先验坐标文件的格式展示的，我们可以从第一次analyze_tseri分析后的analyze_tseri.out文件中提取更新过后的先验文件，作为后面pca分析的先验文件，在终端输入：
   
   即可得到新的先验文件。文件格式之前已在介绍pca模块时给出。
3. 运行后终端中显示正常结束后才算解算正常：（如果出现报错应该好好检查之前提到的所有准备文件，包括列表文件中的数据）
   
4. 正常解算后可以得到analyze_tseri.out、analyze_tseri.resi和analyze_tseri.resi1文件。可以对out文件中各参数项进行提取以获得相应的目标信息。也可以直接对残差序列进行分析。