[发明专利]一种区域CORS坐标时间序列噪声模型的建立方法

摘要

本发明涉及一种区域CORS坐标时间序列噪声模型的建立方法，其特征在于包括步骤1，获取区域CORS站的观测数据；步骤2，采用GNSS数据预处理软件TEQC；步骤3，采用GAMIT/GLOBK基线解算平差；步骤4，建立高精度坐标时间序列；步骤5，基于奇异谱分析的空间滤波；步骤6，基于MLE估计的噪声模型选择；步骤7：顾及有色噪声模型区域CORS速度场估计。本发明通过采用时域性的频域特征分析，对区域CORS时间序列进行空间滤波，为区域CORS时间序列不同噪声模型选择最大似然值进行判断，从而确定N（北）、E（东）、U（高程）方向最优噪声模型，并顾及最优噪声模型估计出区域CORS站速度场。

主权项

一种区域CORS坐标时间序列噪声模型的建立方法，其特征在于包括如下具体步骤：步骤1，获取区域CORS站的观测数据：包括观测单天数据O文件并将其转为RINEX格式、导航N文件以及从IGS网站下载对应于观测时段的SP3精密星历文件；步骤2，采用GNSS数据预处理软件TEQC：包括由GNSS数据预处理软件TEQC对CORS站观测数据进行数据编辑和质量检核，删减不合格测段；步骤3，采用GAMIT/GLOBK基线解算平差：包括制定合理的基线解算策略，引入区域CORS网周边IGS站点的同步观测数据，用GAMIT软件解算基线；将GAMIT基线解算的松弛解h文件用GLOBK进行平差，调用globk、glorg模块解算，从org约束平差文件中获取区域CORS各站坐标时间序列；步骤4，建立高精度坐标时间序列：包括利用三个坐标值之差通过“拉宜达法则”来剔除序列中野值，采用“最邻近准则”算法进行插值；步骤5，基于奇异谱分析的空间滤波：包括采用奇异谱分析来对区域CORS坐标时间序列进行空间滤波，对区域CORS网坐标时间序列的趋势项进行提取，同时对具有周期变化特点的面积较小的区域CORS网，保留其周期项进项噪声分析建模，得到剔除共性误差后的区域CORS站坐标时间序列误差棒；步骤6，基于MLE估计的噪声模型选择：包括基于MLE最大似然估计对空间滤波后的区域CORS坐标时间序列进行噪声分析；利用极大似然估计准则定量测试噪声组合模型；依据极大似然估计函数的拟合自由度给出蒙提卡罗准则下的选择标准：当两种模型的拟合自由度相同时，MLE差值>0的噪声模型更优；当两种模型的拟合自由度相差1时，MLE差值>2.9作为模型显著区分的阈值；当两种模型的拟合自由度相差2时，MLE差值>4.7作为模型显著区分的阈值；步骤7，顾及有色噪声模型区域CORS速度场估计：包括建立区域CORS速度场估计的基于N、U 、E三分量的最优噪声模型的方法。