QOCA（Quasi-observation Combination Analysis）

合并松束解获得地壳形变信息

接收gipsy gamit fonda globk bernese 格式的数据 还有接收sinex和blue book的数据

可以得到站点位置，速度，网络参数，同震形变和震后形变，应变率和剪切律，检测误差

代替globk，直接用gamit生成的单天解进行运算sepa格式的qoca，直接得到时间序列和基线重复性。

1、mk_net 生成先验的站点坐标和速度文件，从其他格式转化到qoca格式

  mk_net infile outfile style datum format

  infile 必备的输入文件名字，outfile 是输出文件名字，style是gipsy globk bbook ipgp usgs bernese stainfo sinex格式

  datum 是WGS84 WGS72 NAD27 NAD83 AGD84 BESSEL 默认值为WGS84

  format 可选的 用默认的格式

  例如：mk_net globk_style.apr qoca_style.net globk

  qoca 格式为u v  w 代表东 北 上 单位m/yr

2、mk_qob 转换qoca格式的先验文件为qob格式

  mk_qob [driver file]

  将合并的解编辑号mk_qob.list运行，生成qob文件

  glf为globk/gamit输出的二进制h文件，gamith为gamit输出的ascii文件

3、mrg_dat 得到多少文件在qoca中分析中被用到

  mrg_dat [infile] [outfile]

  infile 为输入文件名字 是qob个数以及qob列表

4、cov_eigen 检验qob协方差矩阵的状况

  cov_eigen [infile] [outfile] [option] [unuse file]

  infile 为mrg_dat输出的文件 outfile为命名 option为1:输出最大和最小的特征值 2:输出全部的特征值。

5、get_obs_tbl

  get_obs_tbl [infile] [outfile]

  infile为mrg_dat的输出文件，得到每一个站点在qob中

6、get_span 得到每一个站点的时间间隔

  get_span [infile] [outfile]

  infile为mrg_dat得到的qob文件，间隔少于一年在文件名有*，间隔少于两年多于一年，则有%

7、准备qoca.drv分两种格式，一种是sepa格式，一种global格式

  先验坐标文件直接用itrf08_comb.apr

  输入qob文件为mrg_dat的生成文件

  utpm tables文件为pmu.bull_a

  leap sec文件leap.sec

  输出文件 自定义名字.out

  map文件 .map

  重新加权文件reweighting file:

   reference frame:                   WGS84

   reference coordinate, rtime:       geodetic 2000.1995

   solution combination mode:         separately

   resolution :                      resolution   0.0  0.0

   network adjust parameter (net_adj): xyzuvws

   filtering mode:                     forward

  c 同震形变文件

  c 震后形变文件

  global格式：估计速度，同震形变，震后形变

  必备的 跟sepa一样net文件，输入qob文件，输出文件，utpm文件，跳秒文件

  djump 仅仅在sepa格式中用到，

  典型的松约束为坐标为10m，速度为1m/yr。典型的紧约束为1cm坐标，1mm/yr速度，参考框架用紧约束，其他用松约束。

  在global模式下运行需要考虑三个部分：

  检查outliers，检查

8、mk_pos_ps画时间序列

  mk_pos mapfile.list site.list ./

  利用qoca得到的map文件的一个列表和想要得到时间序列站点的站点site的列表，直接列表不需要number

  然后最后跟一个存放ps文件的目录。

9、 mk_bsl_ps生成基线重复性

   mk_bsl_ps baseline.output baseline.list ./

  在qoca的drv中设置baseline file 得到baseline.out文件 将想得到的baseline.out的第一个域内的基线名字输入到baseline.list

  然后存放到一个指定的目录中

10、net_update 更新net文件

   net_update <infile> <outfile> <modfile> <option> <mode>

   输入文件原始的net文件，输出更新的net文件 modfile为qoca输出的修改的netfile option 为

   1 三维坐标加速度 2 二维坐标加速度 3 三维坐标 4 二维坐标 5 三维速度 6 二维速度

   mode 为1 更新所有的站点  2 排除紧约束的站点

   在qoca的global模式下 要输入 modified network file

11、compare_cor 比较net文件的区别

  compare_coor <infile> <outfile> <modfile> <option>

   infile为原始的net文件 outfile为输出文件，modfile在qoca生成修改的文件option同上

12、get_v_rel 得到相对速度解

   get_v_rel <infile> <outfile> <option> <ref site>

    infile为输入的map文件，output为输出的map文件，option为

    1、fonda map文件 2、fonda先验坐标文件 3、globk的速度文件 4、genrelref文件 5、qoca速度文件

   这个也可以输入fonda的map文件，fonda先验坐标文件，globk的速度文件，globk的genrelref的输出文件

    可以生成以上文件的相对速度场，ref_site 为参考站点的名字 对于qoca的map文件 选择5后则在参考站点后加以数字2

13、pps2smap 转换gipsy的输出stacov文件为qoca的sepa的map文件

   pps2smap <aprfil> <infile> <outfil> <mode> <fmt> <rtime>

   aprfile 为qoca格式的坐标速度先验文件，输入文件为单点定位的stacov文件的列表，输出map文件 mode为1 先验速度除去2为速度保持

   fmt 为1输出qocamap文件格式 2为输出qoca site_xyz格式 3为输出site_enu格式 rtime 为参考时间单位为年

14、get_ve_eigen 计算速度矢量误差椭圆的特征值

   get_ve_eigen qoca_global.map > velo_err_ellip.eigen

   得到站点名字  最大轴  最小轴 方位角

15、rot_frame 从NNR框架下转换到相对板块的框架下

   rot_frame <infile> <outfile> <style><plate> <pole_lat><pole_lon><pole_vel>

   框架名字有NAM,EUR,PAC,AFR,IND,ARA,AUS,ANT,NAZ,SAM,SNA

    输入文件 style有两种 1为qoca的net格式，2为qoca的map格式 plate为命名 lat lon为欧拉极的度数 vel为角速度 后面三者任选项

16、calc_bsl_azm计算基线长度和方位角

    calc_bsl_azm <netfile> <infile><outfile> <rtime> <option>

    netfile先验net文件，rtime参考时间，模式1=长度加方位 2=dx，dy，dz 3=de，dn，du

17、syn_bsl_rat 合成两个站点间的基线长度变化率和方位角变化率

   syn_bsl_rat <mapfile> <infile><outfile> <option>

   模式1 为长度和方位角变化率

18、tsf_xyz_neu 从xyz转化为neu

   tsf_xyz_neu <infile><outfile> <option>

   模式1 为xyz到neu 2为neu 到xyz 3 aor到neu 4 ne数字画ne 5ne到ne数字化 6 xyz到neu（globk格式）

19、map_2_apr 将map文件转化为apr文件

   map_2_apr <infile> <outfile> <option>

   模式为1 apr文件转化为小数点坐标十进位坐标  2 apr文件转化为map文件 3 map文件转化为apr文件 4 map文件转化为xyz apr（globk格式）

20、neu_xyz

21、est_pole

  达到一个最佳解需要

      1、设置合适的约束 2、分配合适的相对加权对于每一个qob文件

      3、设置合适的随机扰动处理未知的系统不协调 4、执行合适的矫正最小化系统偏差

      5、选择合适的同震模式和震后模式  6、执行realistic statistics获得合理的不连续估计

    在每一个survey的进行运行sepa的mode以至于生成时间序列分析检测误差，用四个基本的方法检测误差：

   1、检查chi-square的increments

    grep -i subf qoca_sepa.out 得到最后的一列，绝对值大于1，一个典型大于其他的qob文件

    经常错误的站点名字在qob文件中 相同的站点名字在先验坐标文件却不是相同的站点或者错误的先验坐标和速度

    若在out输出有站点的调整为*****，为了矫正这个文件，可以利用运行qoca的sepa的mode单独运行这个qob文件

    在这次运行中，让每一个参数保持松约束，不要估计任何net的参数，不要设置任何link，再去查找坐标矫正，会找到最大的矫正的站

    若有几个站点大于则，不用最大矫正的点，重新运行qocasepa的mode。

    若有些站点没有先验坐标，则qoca默认为用临近的站点的坐标，代替它的先验坐标，这样，则用qoca的运行后更新这些点的坐标即可，

   2、用screening 文件

    在删除大的误差后，处理中等或者小的误差，若有episodic 则用djump 为yes，可以删除同震位移或者相位中心的改变引起。用

    grep -i 站点 sepa.screen 得到站点e n u x y z sigmax sigmay sigmaz 的调整可以明显得到误差

   3、设置加权因子

    利用每一个qob得到的比例因子的平方根，加到qob列表的inlist中，重新运行qob的sepa模式

   4、更新先验坐标，有些新的站点的加入，水平方向5m，垂直方向10m作为检验标准

brief deconstraint：

   为了保证合并的是松束解，在mk_qob的drv文件中 利用 deconstraint file : decon.list

   decon.list里第一行为先验文件名字，第二行前三列为e n u的紧约束单位为米，接三列为想得到的松约束，接三列为xyz想要得到的松约束，接着跟站点名字。

   完成后继续mk_qob的drv 生成qob文件

brief internal constraint

   在sepa的模式下用，可以利用net_adj进行替换实施七参数转换完成

   对于全球的解，直接利用上调整，对于区域的network，得到区域的时间序列，则绑定全球站点挨着这个network比较近的进行约束

   并且分配一定的约束在几个当地站点的水平坐标上

brief resolution

   resolution=1为参数估计全部来于数据，=0 为参数估计来于指定的约束，人为的估计，不能出现在科学文章中

注：在运行得到missing site后，将这些站点放到unuse.site中，编辑mk_qob.drv的文件 重新得到qob文件

接着运行mrg_dat得到qoca的drv的输入文件运行，得到map文件以及out文件。

  sepa的模式的qoca会删除由站点运动引起的先验速度，这有利于研究稳态的蠕变，常规速度比如是震后形变或者季节性变化

  若要保持站点运动在输出的map文件中，应该在map的输出文件的一项中 后加 v

  ITRF 框架是CM参照框架 就是整个地球加大气海洋的质量

模块分析：

   详细看附件的内容，感谢华东师范大学的董大南老师的讲授以及文档

以及文档详细信息   http://qoca.jpl.nasa.gov/

Analyze_strain.docx

Analyze_tseri.docx

mload1.doc

mstokes.docx

pca.docx

st_filter.docx