劉慕溪

(廣東省核工業(yè)地質(zhì)局測(cè)繪院，廣東 廣州510800)

0 引言

隨著衛(wèi)星、導(dǎo)航定位技術(shù)的發(fā)展，空間定位測(cè)量技術(shù)(尤其是GPS測(cè)量)在測(cè)繪中應(yīng)用越來越廣泛;GNSS系統(tǒng)，如美國的GPS、俄羅斯的GLONASS和北斗衛(wèi)星導(dǎo)航系統(tǒng)，都能為用戶提供海量的數(shù)據(jù)信息［1］。全球已建成上百個(gè)全球連續(xù)運(yùn)行參考站和數(shù)千個(gè)區(qū)域站，不但可以觀測(cè)到地表位移，進(jìn)行形變分析以及板塊運(yùn)動(dòng)研究，而且還能觀測(cè)到震后形變等微小地殼形變，為地殼研究提供基本的原始資料。但測(cè)量不可避免的存在誤差，GPS測(cè)量中的誤差主要體現(xiàn)在與衛(wèi)星有關(guān)的誤差，與接收機(jī)有關(guān)的誤差，與信號(hào)傳播有關(guān)的誤差，如何剔除GPS觀測(cè)值中的誤差，降低噪聲對(duì)觀測(cè)成果的影響成為一個(gè)重要問題。對(duì)流層濕延遲、電離層延遲、天線相位中心誤差和衛(wèi)星軌道誤差等都是噪聲的潛在來源，隨著算法和模型的改進(jìn)，該類噪聲能得到很好的消除或減弱。此外GPS觀測(cè)值中包含著時(shí)空相關(guān)的噪聲，以及非構(gòu)造信號(hào)等共性誤差，對(duì)觀測(cè)結(jié)果帶來不確定性［2］。本文討論了基于主成分分析的噪聲剔除方法，結(jié)果表明主成分分析法可以較好地消除GPS解中的共性誤差，提高GPS測(cè)量定位的精度。

1 主成分分析

1.1 主成分分析概述

主成分分析又叫經(jīng)驗(yàn)正交函數(shù)分析，是現(xiàn)代數(shù)據(jù)分析的一個(gè)有效工具，廣泛應(yīng)用于經(jīng)濟(jì)、金融、測(cè)繪等領(lǐng)域，是一種非參數(shù)數(shù)據(jù)處理方法，計(jì)算簡單、效率高［3］。主成分分析法在測(cè)量控制網(wǎng)的靈敏度分析、粗差檢驗(yàn)等得到了較好應(yīng)用。通過尋找觀測(cè)數(shù)據(jù)的方差－協(xié)方差陣的前k個(gè)主分量，在保持?jǐn)?shù)據(jù)信息損失最小的前提下，經(jīng)過線性轉(zhuǎn)換將原始自變量中相關(guān)的維數(shù)消除，轉(zhuǎn)換到低維向量空間，低維中各主分量是相互正交的，從混雜的數(shù)據(jù)中提取出相關(guān)的信息。由于前幾個(gè)主分量綜合了原始數(shù)據(jù)的最大信息量，因此不會(huì)對(duì)數(shù)據(jù)分析結(jié)果帶來太大的影響，同時(shí)避免了法方程病態(tài)問題的出現(xiàn)，得到參數(shù)的精確估值，提高了解的準(zhǔn)確性和可靠性，同時(shí)可以揭示數(shù)據(jù)中隱含的一些規(guī)律及結(jié)構(gòu)特征［4］。

1.2 主成分分析原理

觀測(cè)P個(gè)變量x1，x2，…，xp，n個(gè)樣本的數(shù)據(jù)資料數(shù)據(jù)矩陣為:

主成分分析就是將p個(gè)觀測(cè)變量綜合成為p個(gè)新的變量(綜合變量)，即

簡寫為:

要求模型滿足以下條件。

(1)各主分量不相關(guān)，即Fi，F(xiàn)j互不相關(guān)(i≠j，i，j=1，2，…，p)。

(2)F1的方差大于F2的方差大于F3的方差，依次類推。

于是，稱F1為第1主成分，F(xiàn)2為第2主成分，依此類推，有第p個(gè)主成分。主成分又叫主分量。這里aij稱為主成分系數(shù)。通過前k個(gè)分量能夠反應(yīng)隨機(jī)變量X的主要特征，該數(shù)據(jù)處理方法稱之為主成分分析。

主成分分析法一般分3個(gè)步驟進(jìn)行求解: (1)首先對(duì)分析數(shù)據(jù)組成數(shù)據(jù)矩陣X(m，n)，m為數(shù)據(jù)類型，n為數(shù)據(jù)樣本個(gè)數(shù);(2)數(shù)據(jù)中心化，即將各個(gè)觀測(cè)數(shù)據(jù)減去均值得到新的數(shù)據(jù)矩陣; (3)計(jì)算新矩陣的方差、協(xié)方差，獲取特征向量，實(shí)現(xiàn)在重構(gòu)誤差最小的條件下對(duì)數(shù)據(jù)維數(shù)的有效簡化。

1.3 GPS噪聲剔除中應(yīng)用

以GPS臺(tái)站網(wǎng)的坐標(biāo)時(shí)間序列為原始數(shù)據(jù)，每個(gè)臺(tái)站的坐標(biāo)分量殘差時(shí)間序列排列起來形成一個(gè)(m＞n)的矩陣X，m表示觀測(cè)歷元數(shù)，n表示觀測(cè)類型，數(shù)據(jù)矩陣如下:

做如下定義:

即V構(gòu)成X的正交基底，矩陣X按照KLE展開可得

ak(ti)可由下式求出:

式中，ak(t)是第k個(gè)主成分，vk(x)是對(duì)應(yīng)主成分的響應(yīng)特征矩陣，分別代表時(shí)間特征和空間響應(yīng)［5］。

2 實(shí)例分析

本文對(duì)連續(xù)運(yùn)行的IGS基準(zhǔn)站2001－2007的觀測(cè)數(shù)據(jù)進(jìn)行分析。首先對(duì)GPS觀測(cè)數(shù)據(jù)進(jìn)行預(yù)處理，GPS數(shù)據(jù)預(yù)處理是指及時(shí)對(duì)外業(yè)觀測(cè)數(shù)據(jù)進(jìn)行檢驗(yàn)，探測(cè)數(shù)據(jù)中的粗差，剔除那些質(zhì)量較差的數(shù)據(jù)，以提高觀測(cè)數(shù)據(jù)的質(zhì)量。在數(shù)據(jù)預(yù)處理的過程中，使用TEQC軟件，對(duì)GPS數(shù)據(jù)進(jìn)行質(zhì)量分析［6］。數(shù)據(jù)預(yù)處理后按照如下步驟進(jìn)行數(shù)據(jù)處理［7］。

(1)采用GAMIT軟件，通過無基準(zhǔn)方法解算出各GPS站點(diǎn)的三維坐標(biāo)及其方差－協(xié)方差陣，以每天24 h的GPS觀測(cè)數(shù)據(jù)為基本單位形成單天解。無基準(zhǔn)算法可以有效地避免基準(zhǔn)誤差的干擾，便于將內(nèi)部獨(dú)立的坐標(biāo)基準(zhǔn)轉(zhuǎn)換到外部統(tǒng)一的參考框架。

(2)將獲取的單天解進(jìn)行平差，獲得單天解坐標(biāo)序列。對(duì)坐標(biāo)序列去均值和趨勢(shì)項(xiàng)處理后(此處減去均值是為后續(xù)應(yīng)用主成分分析法進(jìn)行噪聲分析做準(zhǔn)備)，得到GPS基準(zhǔn)站單天解連續(xù)坐標(biāo)時(shí)間序列，限于篇幅以CAND站為例，圖1為CAND站垂向分量坐標(biāo)序列。

(3)將上一步生成的殘余坐標(biāo)時(shí)間序列進(jìn)行主成分分析，把時(shí)間序列分解成時(shí)間域的主分量和空間域的特征分量，按其貢獻(xiàn)率排列，見表1。

表1 主成分分析后得到的前3個(gè)主分量的貢獻(xiàn)值

由表1可知在3個(gè)坐標(biāo)方向主成分分析后前3個(gè)分量的累計(jì)貢獻(xiàn)率達(dá)到90%，綜合了絕大部分的原始信息，故采用前3個(gè)分量作為主分量，對(duì)GPS單天解形成的連續(xù)坐標(biāo)時(shí)間序列進(jìn)行處理，得到站點(diǎn)單天解的殘余序列(圖2)。

圖1 CAND垂向分量坐標(biāo)序列

圖2 經(jīng)過主成分分析剔除噪聲后的坐標(biāo)序列

通過主成分分析法處理后，GPS單天解的坐標(biāo)序列可看出，GPS坐標(biāo)序列的波動(dòng)趨于平穩(wěn)，波動(dòng)振幅有所減小，對(duì)非構(gòu)造信號(hào)、時(shí)空相關(guān)的誤差源引起的站點(diǎn)位移進(jìn)行了剔除。對(duì)數(shù)據(jù)進(jìn)行分析可知，濾波之前垂向分量的均方根誤差為1.751 8，濾波之后為1.556 0，即通過主成分分析濾波方法分離出的誤差約減少了站點(diǎn)11%的不確定度，提高了站點(diǎn)坐標(biāo)的精度、信噪比。此外，濾波前坐標(biāo)序列呈現(xiàn)出明顯的周期性變化，濾波后坐標(biāo)序列的周期性明顯減弱，表明坐標(biāo)序列中存在的誤差呈現(xiàn)出季節(jié)性變化趨勢(shì)，有待于進(jìn)一步分析。

3 結(jié)論

本文基于主成分分析的噪聲剔除方法，針對(duì)GPS信號(hào)中存在的非構(gòu)造信號(hào)誤差，時(shí)空相關(guān)誤差等，能有效的剔除或減弱上述誤差的影響，可提高GPS測(cè)量、定位的精度與可靠性，對(duì)GPS的應(yīng)用具有一定的的參考價(jià)值。

［1］ 夏金偉，趙東保.TEQC數(shù)據(jù)質(zhì)量分析［J］.地理空間信息，2012，(4):60－65.

［2］ 田云鋒，沈正康.GPS坐標(biāo)時(shí)間序列中非構(gòu)造噪聲的剔除方法研究進(jìn)展［J］.地震學(xué)報(bào)，2009，(1):68－81，117.

［3］ Shlens J.A Tutorial on Principal Component Analysis，Center for Neural Science［J］.New York University New York City，NY 10003－6603 and Systems Neuro blogy Laboratory，Salk Insitute for Biological Studies La Jolla，CA 92037，April 22，2009.

［4］ 馮光財(cái)，陳正陽.基于主成分回歸的GPS高程曲面擬合［J］.測(cè)繪科學(xué)，2007，32(1):51－52.

［5］ 殷海濤，甘衛(wèi)軍，熊永良，等.PCA空間濾波在高頻GPS定位中的應(yīng)用研究［J］.武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2011，(7):825－829.

［6］ 賀小星，周世健.某城區(qū)D級(jí)GPS靜態(tài)控制網(wǎng)設(shè)計(jì)與實(shí)施［J］.江西科學(xué)，2012，30(3):295－298.

［7］ 王 敏，張祖勝，許明元，等.2000國家GPS大地控制網(wǎng)的數(shù)據(jù)處理和精度評(píng)估［J］.地球物理學(xué)報(bào)，2005，(4):817－823.