[发明专利]一种构建GNSS区域大地参考框架的方法

权利要求书

1.一种构建GNSS区域大地参考框架的方法，其特征在于GNSS区域大地参考框架的构建具体包括以下步骤：

步骤1：区域子网的松约束解及协方差矩阵解算

利用GAMIT软件处理区域子网持续测站的监测数据，得到区域子网无约束解及其协方差矩阵，利用GLOBK软件并入全球网松约束解；

步骤2：最优约束方案的确定

2-1：内约束台站的选择标准

区域性台站网内约束台站的中心尽量靠近整体网的中心，内约束台站至少有三年的长期连续观测；台站的坐标位置残差小于5mm，速度残差小于3mm/y；

2-2：内约束方案及台站相对权重的确定

以ITRF给出的GNSS台站径向分量为先验值，经3～4次迭代提取每个内约束台站径向分量的方差，当方差大于某一标准时，将从内约束站列表中剔除精度较差的内约束台站；对旋转分量进行等权拟合，计算水平残差偏差指数，当其偏差指数超过某一阈值时，从内约束站点列表中剔除，以台站残差与偏差指数之比作为各内约束台站的相对权重，对旋转分量进行加权拟合，剔除野值，经过3～4次迭代，直至不再出现较大偏差的台站，并确定台站相对权重，得到最优约束方案；

步骤3：Helmert转换及瞬时测站坐标的确定

利用所有自由网络的松约束坐标解进行Helmert转换，并通过最小二乘(LS)平差，在区域参考系实施过程中保持整个相对几何结构不发生内部变形，对单日松弛解坐标解进行参数变换，以保持松散约束解的内部几何结构；

步骤4：区域滤波的实施

利用主分量分析实现区域滤波，所述主分量分析为经验正交函数分析(EmpiricalOrthogonal Function,EOF)或特征向量分析，利用GNSS台站坐标时间序列作正交分解得到时间域的主分量和空间域的特征向量；

步骤5：参考历元坐标、速度及非线性项的确定

对于每个测站而言，其约束坐标时间序列不仅包括线性趋势还包括非线性变化，利用时间序列分析得到每个台站的在参考历元t₀的测站坐标及速度以及非线性项，所述时间序列采用下述14式表示的模型进行分析：

式中：表示季节项振幅；k＝1表示周年项,；k＝2表示半周年项；ω^k和分别表示相应的频率和初始相位；表示振幅的季节项变化；表示t₁开始，t₂历元结束的窗函数的振幅；表示震后形变；t^k为震后形变开始历元；τ^k为特征衰减时刻；为震后指数衰减模型或对数衰减模型；表示其它非线性项。