劉建勇 蓋忠奎

摘要：GPS數(shù)據(jù)處理平差基準(zhǔn)需要通過高等級基準(zhǔn)點進行傳遞。高等級基準(zhǔn)點的質(zhì)量影響GPS網(wǎng)坐標(biāo)解算的準(zhǔn)確性，故在得到基線解算之后，如何提前探測松弛解中基準(zhǔn)站點質(zhì)量優(yōu)劣，以便在后續(xù)平差時放心對合格的基準(zhǔn)站進行約束，剔除或者松弛異常基準(zhǔn)點是本文探討的問題。本文將研究在松弛解中應(yīng)用HELMERT模型對松弛解的基準(zhǔn)進行一致性檢驗，提前探測出有問題的基準(zhǔn)站點，有利于得到更為準(zhǔn)確的平差結(jié)果。

關(guān)鍵詞：HELMERT；松弛解；基準(zhǔn)站；坐標(biāo)一致性

一、概述

在高精度GPS網(wǎng)解算中，基準(zhǔn)站是作為高級控制點來實現(xiàn)坐標(biāo)框架傳遞的，中國現(xiàn)在推薦的坐標(biāo)系統(tǒng)為CGCS2000，其框架歷元為2000.0，距今已經(jīng)有十多年的時間跨度，在這個時間跨度內(nèi)，中國的地殼板塊及其子板塊發(fā)生著不同的運動趨勢，尤其在中國西部地區(qū)。所以GPS解算中，尤其在大跨度GPS網(wǎng)所選用的基準(zhǔn)站的CGCS2000的坐標(biāo)跟現(xiàn)在這些站的瞬時坐標(biāo)是否具有一致性和和較好的附和性，是個十分重要的問題，因為如果不具有一致性，則將誤差引入到待解算點去。所以必須提前對基準(zhǔn)站進行一致性驗證。

HELMERT（赫爾默特）能夠描述兩套坐標(biāo)之間在平移、尺度、定向的相關(guān)變化特性，并通過轉(zhuǎn)換殘差可以探測出不符合相似變化的異常坐標(biāo)。本文將研究使用HELMERT（赫爾默特）法來實現(xiàn)基準(zhǔn)站坐標(biāo)一致性分析的目的，詳細(xì)研究HELMERT模型的理論，及其在轉(zhuǎn)換分析中的方法。

二、HELMERT模型簡介

HELMERT簡介：HELMERT主要是通過對兩套坐標(biāo)進行對比，通過兩套坐標(biāo)公用點求取赫爾默特七參數(shù)（三個平移參數(shù)、三個旋轉(zhuǎn)參數(shù)、一個比例因子），并按照七參數(shù)將檢驗坐標(biāo)系坐標(biāo)轉(zhuǎn)化到參考坐標(biāo)下，同時求取殘差。

參數(shù)平差模型采用最小約束（Minimum Constraint Solution）它的實現(xiàn)主要是對比兩套坐標(biāo)系，前提是承認(rèn)這兩套系統(tǒng)其質(zhì)量重心是一致的。通過殘差可以分析待檢驗坐標(biāo)系中哪些站不符合參考坐標(biāo)，可以探測出變化異常點，以備平差前將其剔除。

HELMERT轉(zhuǎn)化模型如下：在大區(qū)域或者全球的GPS網(wǎng)一般使用大地直角坐標(biāo)系進行轉(zhuǎn)換。其轉(zhuǎn)換參數(shù)解算和坐標(biāo)系轉(zhuǎn)換模型如下：

X′=（1+k）·Rz（rz）·Ry（ry）·Rx（rx）·（X+X0）（2.1）

X0平移參數(shù)；Ri（ri）為X/Y/Z三個方向的旋轉(zhuǎn)參數(shù)；k比例因子。

利用模型2.1可以實現(xiàn)原坐標(biāo)系到參考坐標(biāo)系的轉(zhuǎn)換參數(shù)。同時，利用經(jīng)過最小二乘平差的轉(zhuǎn)換參數(shù)，將原坐標(biāo)系坐標(biāo)轉(zhuǎn)換至參考坐標(biāo)系，然后求取轉(zhuǎn)換后的坐標(biāo)于參考坐標(biāo)的差異，及殘差。通過殘差分析可以發(fā)現(xiàn)比較大的殘差，即為異常點。

三、案例分析

1.數(shù)據(jù)概述：本文數(shù)據(jù)采用2011年288（doy）～294一周的某市CORS網(wǎng)數(shù)據(jù)。待解CORS站14個，基準(zhǔn)站選用周邊5個連續(xù)運行站為SHAO（上海站）、WUHN（武漢站）、BJFS（北京站）、ULAB、GUAO（烏魯木齊站）。每個時段周邊基準(zhǔn)站均齊全，且數(shù)據(jù)時段長度為23.5小時。解算使用精密星歷。

2.解算步驟

1）得到松弛約束解

使用GAMIT軟件解算基線，在參數(shù)設(shè)置時，對固定站采用N/E/U方向分別給與：0.050 0.050 0.100的約束，得到各個固定站的松弛約束解。

基線解算和平差采用軟件：美國麻省理工學(xué)院（MIT）和海洋研究所（SIO）聯(lián)合開發(fā)的GAMIT/GLOBK軟件10.40版本。

2）HELMERT模型驗證坐標(biāo)一致性

基線解算可以得到單天解，單天解中有每個固定站的單天坐標(biāo)。提取各個固定站的單天解坐標(biāo)然后與各個站的ITRF2008坐標(biāo)2005.0歷元進行HELMERT轉(zhuǎn)換。查驗各個站轉(zhuǎn)換殘差。

3）進行單時段全約束平差，查驗平差開方檢驗指標(biāo)。

4）根據(jù)HELMERT的殘差指標(biāo)和平差的開方檢驗指標(biāo)定位異常點。

5）對異常點或者時段進行合理處理后綜合平差。

3.殘差分析

采用HELMERT模型得到基準(zhǔn)站單天松弛解坐標(biāo)與ITRF2008框架的2005.0歷元的坐標(biāo)的轉(zhuǎn)換殘差如下：

根據(jù)經(jīng)驗將轉(zhuǎn)換殘差闋值設(shè)置為：N/E方向4cm；U方向為8個cm。

如上表分析可以看到轉(zhuǎn)換殘差2014.290當(dāng)天BJFS站E方向殘差大于4cm。故BJFS站在2014.290時段可能存在問題。

然后，本文對7個時段分別分進行了單天平差，約束所有基站站，其單天解開方檢驗如下表：

表4.1可以看出2014.290開方檢驗比其他時段都大，也驗證了2014.290時段的基站站不一致，正好與HELMERT檢驗結(jié)果一致。所有該時段BJFS站有問題，故本次測試將2014.290時段剔除。

其他時段，坐標(biāo)轉(zhuǎn)換殘差均符合標(biāo)準(zhǔn)。所以，剔除2014.290天之后，在ITRF2008坐標(biāo)系下平差可以放心使用基準(zhǔn)站的ITRF2008坐標(biāo)，且可以施加較強約束。

四、結(jié)論

在GPS數(shù)據(jù)處理，尤其是對大區(qū)域或者全國區(qū)域的CORS站處理時，必須嚴(yán)格驗證所選基站站的坐標(biāo)一致性。以及在長期數(shù)據(jù)時間序列解算時候，驗證單時段選用基站站的一致性是保證時間序列真實的前提。本文也通過實例分析可以得到結(jié)論：平差的開方檢驗指標(biāo)對基準(zhǔn)坐標(biāo)一致性的敏感程度不高，但通過HELMERT分析可以看出其對異常坐標(biāo)的敏感度比較高，可以及時準(zhǔn)確發(fā)現(xiàn)不一直的異常點。所以赫爾模特模型可以作為高精度數(shù)據(jù)處理的很有效的坐標(biāo)一致性檢驗?zāi)Ｐ图右詰?yīng)用。

參考文獻(xiàn)：

[1]張志偉，童軍.基于赫爾模特法港口坐標(biāo)系統(tǒng)投影變形處理[J].港工技術(shù)，2005

[2]茹樹青，吉長東.RTK坐標(biāo)轉(zhuǎn)換中四參數(shù)法與七參數(shù)法精度比較[J].東北水利水電， 2006（5）：67～68.

[3]曹陽，張家根.城市規(guī)劃空間數(shù)據(jù)Helmert大地坐標(biāo)轉(zhuǎn)換模型研究[J].東北水利水電，2006

[4]Rolf Dach，Urs Huge iversity of Bern， January 2007

[5]T.A.Herring， R.W.King， S.C.McClusky，GLOBK Reference Manual， Department of Earth， Atmospheric， and Planetary Sciences Massachusetts Institute of Technology， December 2010

作者簡介：劉建勇（1980.06-），男，河北石家莊人。大學(xué)，研究方向：航測與制圖。

蓋忠奎（1981.08-），男，河北平山人。大學(xué)，研究方向：測繪及監(jiān)測。