[发明专利]一种顾及板块运动异向性与局部稳定性的区域参考框架维持方法

权利要求书

1.一种顾及板块运动异向性与局部稳定性的区域参考框架维持方法，其特征在于，包括如下步骤：

步骤一、获取单日区域松弛解；

步骤二、通过对单日区域松弛解进行处理，生成单天时间序列，进行区域基准站的稳定性分析，均匀选取稳定的特征点；

步骤三、通过对单日区域松弛解进行处理，获得CGCS2000结果、ITRF框架瞬时结果，然后进行赫尔默特转换，均匀选取特征点；

步骤四、通过对单日区域松弛解进行解析与恢复，构造区域整体法方程；

步骤五、综合步骤二得到的稳定的特征点、步骤三得到的特征点，选取拟稳点；

步骤六、引入拟稳点构造的基准条件与法方程联合求解，得到待定点CGCS2000的最小约束解。

2.如权利要求1所述的一种顾及板块运动异向性与局部稳定性的区域参考框架维持方法，其特征在于：所述步骤二中进行区域基准站的稳定性分析，均匀选取稳定的特征点的具体过程是：

其中，a为初始位置；b为速率；c,d,e和f分别为年、半年周期项系数；g为偏移；h为震后速率变化；k为震后速率衰减模型；v为误差；t_i为时间；H为阶梯函数，y(t_i)为单站、单分量坐标序列；

起算区域基准站的选取应满足四个原则：①连续性：测站在近1年内有连续观测；②稳定性：站点的坐标时间序列具有稳定“可知”的点位变化速度，稳定性好；③高精度：速度场精度应优于3mm/a；④平衡性：所选站点尽量均匀分布；

按照上述原则，选取区域内稳定的、残差变化最小的点作为特征点。

3.如权利要求1所述的一种顾及板块运动异向性与局部稳定性的区域参考框架维持方法，其特征在于：所述步骤三、通过对单日区域松弛解进行处理，获得CGCS2000结果、ITRF框架瞬时结果，然后进行赫尔默特转换，均匀选取特征点中的转换模型为：

X_ref＝(1+D)·R_x(ε_x)·R_y(ε_y)·R_z(ε_z)·X_in+T (2)

式中：X_ref为参考坐标矩阵；X_in为待转换坐标矩阵；R_i(ε_i)为沿三个坐标轴的欧拉角；T为平移向量矩阵。

4.如权利要求1所述的一种顾及板块运动异向性与局部稳定性的区域参考框架维持方法，其特征在于：所述赫尔默特转换的过程为：

第一步：整体转换异常点剔除：将获取的所有待定点的松弛解与解算的ITRF框架瞬时结果进行赫尔默特转换，分析转换残差，剔除异常站点；

第二步：局部转换异常点分析：分子区再次进行赫尔默特转换，从各区域内选出均匀分布的，转换残差较小且稳定的点作为特征点。

5.如权利要求1所述的一种顾及板块运动异向性与局部稳定性的区域参考框架维持方法，其特征在于：所述步骤四、通过对单日区域松弛解进行解析与恢复，构造区域整体法方程的过程是：

步骤1、对单日区域松弛解搜索公共跟踪站，对参数进行重排，生成转换矩阵；

步骤2、恢复法方程流程：

计算法方程系数：

计算法方程常数：

计算二次型：

其中为输入基线文件中提取的参数估值，为参数的方差协方差矩阵，b为法方程常数，为单位权方差因子；

步骤3、获取无约束的参数解及法方程：从估计解方差阵中减去约束矩阵，具体表示为：

N^unc＝D^est-D^const (6)

其中N^unc为去除先验约束后的法方程，D^est为估计解得方差阵，D^const约束矩阵；

步骤4、采用参数重排生成的转换矩阵改造得到的解，采用参数预消除算法对非公共参数进行消除，得到所有天公共测站的无约束解与法方程。