[發明專利]一種區域CORS坐標時間序列噪聲模型的建立方法

在審

• 申請號: 201410614749.5
• 申請日: 2014-11-04
• 公開（公告）號: CN104392414A
• 公開（公告）日: 2015-03-04
• 發明（設計）人: 曹奇;張迎燕;劉志強;岳東杰;王可偉;王性猛;王海;李成仁
• 申請（專利權）人: 河海大學
• 主分類號: G06T5/00
• 分類號: G06T5/00;G06F17/50;G06T7/00;G01S19/03
• 代理公司: 南京經緯專利商標代理有限公司 32200
• 代理人: 吳樹山

1.一種區域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三分量的最優噪聲模型的方法。

2.?根據權利要求1所述的一種區域CORS坐標時間序列噪聲模型的建立方法，其特征在于步驟3所述用GAMIT解算基線，是指其結果相對精度達到10^-9，對步驟1所述單天數據O文件提取的Postfit?nrms均小于0.25；檢驗GLOBK平差得到平均NRMS小于3.0的N、E、U各站坐標時間序列和坐標重復性WRMS。

3.根據權利要求1所述的一種區域CORS坐標時間序列噪聲模型的建立方法，其特征在于步驟5所述區域CORS網各坐標時間序列，是指其各基準站水平誤差區間為[-15，15]mm，高度坐標分量誤差區間為[-30，30]mm，高程坐標分量誤差為水平坐標分量誤差的兩倍。

4.根據權利要求1所述的一種區域CORS坐標時間序列噪聲模型的建立方法，其特征在于步驟6所述對空間濾波后的區域CORS坐標時間序列進行噪聲分析，是指濾波后的區域CORS坐標時間序列選擇可變白噪聲或者純白噪聲模型，根據蒙提卡羅準則來判斷由MLE估計算出其最大似然估值差值，對N、E、U三分量坐標時間序列譜指數與FN、RWN、WH參數估值得到最優有色噪聲模型。

5.根據權利要求1所述的一種區域CORS坐標時間序列噪聲模型的建立方法，其特征在于步驟7所述建立區域CORS速度場估計的基于N、U、?E三分量的最優噪聲模型的方法，是指在ITRF框架下，對區域CORS網各基準站點N、U方向坐標時間序列采用“WH+FN”模型估計速度與精度，E方向坐標時間序列采用“WH+FN+RWN”模型估計速度與精度；與之對比的是對N、U、E方向坐標時間序列采用“WH”模型進行速度精度估計。