国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

航空重力測量測線網(wǎng)平差中的粗差處理*

2012-11-14 13:48:26周波陽羅志才
大地測量與地球動力學(xué) 2012年2期
關(guān)鍵詞:抗差交叉點(diǎn)線網(wǎng)

周波陽 羅志才 林 旭 周 浩

(1)武漢大學(xué)測繪學(xué)院,武漢 430079

2)武漢大學(xué)地球空間環(huán)境與大地測量教育部重點(diǎn)實(shí)驗(yàn)室,武漢430079)

航空重力測量測線網(wǎng)平差中的粗差處理*

周波陽1)羅志才1,2)林 旭1)周 浩1)

(1)武漢大學(xué)測繪學(xué)院,武漢 430079

2)武漢大學(xué)地球空間環(huán)境與大地測量教育部重點(diǎn)實(shí)驗(yàn)室,武漢430079)

航空重力測量測線網(wǎng)平差的主要任務(wù)是進(jìn)行系統(tǒng)誤差的補(bǔ)償。為了削弱粗差觀測值的影響,提高補(bǔ)償效果,引入閾值法和抗差估計(jì),并采用模擬數(shù)據(jù)驗(yàn)證了這兩種方法的有效性和可靠性。通過比較發(fā)現(xiàn),抗差估計(jì)在航空重力測線網(wǎng)平差中具有更好的實(shí)用價(jià)值。

航空重力測量;測線網(wǎng)平差;閾值法;抗差估計(jì);權(quán)函數(shù)

1 引言

在航空重力測量中,為了盡可能地構(gòu)成重復(fù)觀測,正確評價(jià)測量精度,測線往往布設(shè)成交叉網(wǎng)狀,在交叉點(diǎn)處兩條測線上的重力值理論上應(yīng)該是相等的,但由于測量平臺在高空作業(yè),處于不斷高速的運(yùn)動狀態(tài),觀測環(huán)境較為復(fù)雜,各類觀測值不可避免地受到“污染”,從而導(dǎo)致觀測值在交叉點(diǎn)處存在不符值。采用低通數(shù)字濾波等技術(shù)消除高頻噪聲后,觀測值中仍存在各種誤差源的干擾和影響。這些干擾和影響按其表現(xiàn)形式可分為偶然誤差、系統(tǒng)誤差以及粗差。偶然誤差的原因十分復(fù)雜,空氣擾動、機(jī)器零件的摩擦、氣壓、氣溫的變化等,以及它們的綜合影響均能產(chǎn)生;系統(tǒng)誤差主要由厄特弗斯效應(yīng)、姿態(tài)誤差、衛(wèi)星星歷誤差等引起;GPS、加速度計(jì)等的工作狀態(tài)反生異常時(shí)觀測值中容易出現(xiàn)粗差。

航空重力測量測線網(wǎng)平差中的粗差處理通常和系統(tǒng)誤差的補(bǔ)償同步進(jìn)行,忽略粗差直接采用最小二乘來估計(jì)系統(tǒng)誤差模型參數(shù)時(shí),所建立的誤差模型并不能真實(shí)反映系統(tǒng)誤差的變化規(guī)律。平差后,大粗差的影響也會被分配到許多正常觀測值上,從而扭曲平差結(jié)果,參數(shù)估值的效率和可靠性就會受到嚴(yán)重的損害[1]。有效地消除或削弱觀測數(shù)據(jù)中粗差對于參數(shù)估值的不良影響,是提高航空重力測量網(wǎng)平差精度的關(guān)鍵環(huán)節(jié)。一般來說,對粗差的處理方式主要有3類:第一類即計(jì)算機(jī)輔助圖形可視化法,通過人機(jī)交互的方式來實(shí)現(xiàn)粗差的探測和剔除。當(dāng)前數(shù)據(jù)處理都趨于自動化,人為干預(yù)排除粗差的方式將被淘汰[2];第二類是在平差前采用某些方法對粗差定位并剔除,Dixon準(zhǔn)則[3]、Grubbs準(zhǔn)則[3]、χ2檢驗(yàn)[4]、τ標(biāo)準(zhǔn)檢驗(yàn)[4]、F檢驗(yàn)[5]、t檢驗(yàn)[6]、小波分析方法[3,7]以及閾值法[3]可以對粗差定位。其中Dixon準(zhǔn)則無需估算樣本均值、樣本標(biāo)準(zhǔn)差和殘余誤差,在將所有觀測值按照從小到大的順序排列后,Dixon準(zhǔn)則認(rèn)為粗差觀測值僅存在于重排序之后的數(shù)據(jù)兩端。Dixon準(zhǔn)則僅適合于小樣本(3<n<30),采用Dixon準(zhǔn)則進(jìn)行小樣本數(shù)據(jù)的粗差探測能夠得到較嚴(yán)密的結(jié)果[8]。Grubbs檢驗(yàn)、F檢驗(yàn)等統(tǒng)計(jì)檢驗(yàn)法都派生于最小二乘估計(jì),對于偏離正態(tài)分布的數(shù)據(jù)是不可靠的,它們幾乎不可能把剔除粗差與估計(jì)量的計(jì)算分開,利用不可靠的估計(jì)量來進(jìn)行統(tǒng)計(jì)檢驗(yàn)就很難保證這種統(tǒng)計(jì)檢驗(yàn)方法的正確性[2]。小波分析方法實(shí)質(zhì)上是對觀測值信號進(jìn)行濾波,利用小波變換的高低頻分離的特點(diǎn),可在不丟失原始信號重要信息成分的前提下,對信號進(jìn)行放大并對其進(jìn)行特征提取,它也是一種較為有效的粗差探測方法,但算法較為復(fù)雜。閾值法計(jì)算方便簡單,適用于觀測數(shù)據(jù)的快速粗差探測;第三類則是基于多維M估計(jì)建立參數(shù)模型的抗差解式,采用等價(jià)權(quán)思想將其轉(zhuǎn)化為最小二乘解形式[9]。采用迭代的抗差M估計(jì),既保留了最小二乘法的優(yōu)越性,又保證被估參數(shù)既能抵制模型誤差又能抗拒粗差的影響[1,10]。本文將抗差M估計(jì)和閾值法應(yīng)用于航空重力測量測線網(wǎng)平差中的粗差處理,并通過模擬計(jì)算對這兩種方法進(jìn)行評價(jià)。

2 系統(tǒng)誤差補(bǔ)償?shù)膬刹椒?/h2>

航空重力測量的測線網(wǎng)平差實(shí)際上是系統(tǒng)誤差的補(bǔ)償問題,進(jìn)行系統(tǒng)誤差補(bǔ)償?shù)囊话惴椒樽詸z校測線網(wǎng)平差法,但自檢校平差會導(dǎo)致法方程奇異,為了解決這個(gè)問題,將自檢校測線網(wǎng)平差過程進(jìn)行簡化,將誤差理論模型簡化為實(shí)用模型,即為誤差補(bǔ)償?shù)膬刹椒ǎ?1-13]。

2.1 交叉點(diǎn)平差

第一階段使用條件平差法對測線交叉點(diǎn)進(jìn)行平差。在主測線、副測線的交叉點(diǎn)處可建立如下的誤差方程:

式中,Δgik與Δgjk分別為主測線i、副測線j在交叉點(diǎn)k處的重力異常,Δgik-Δgjk表示交叉點(diǎn)不符值,Vik、Vjk為改正數(shù)。對于具有多個(gè)交叉點(diǎn)的測線網(wǎng),可以寫出交叉點(diǎn)條件方程的矩陣形式為

其中,V由每一條測線上交叉點(diǎn)處的改正數(shù)組成,系數(shù)陣B是由1、-1和大量的0構(gòu)成,W為交叉點(diǎn)不符值。依據(jù)最小二乘平差理論可解得

一般來說,測線上各個(gè)測點(diǎn)為獨(dú)立等精度觀測,權(quán)陣P為單位陣。

2.2 測線系統(tǒng)誤差估計(jì)

以觀測時(shí)間t為自變量,以式得到的交叉點(diǎn)不符值的改正數(shù)V為虛擬觀測值,遍歷所有測線對每個(gè)交叉點(diǎn)建立如下形式的誤差方程

vi為V中的第i個(gè)元素,f(t)為系統(tǒng)誤差模型,Δi為噪聲,m為該測線上的交叉點(diǎn)個(gè)數(shù)。f(t)有如下兩種形式:

1)一般多項(xiàng)式模型

2)三角多項(xiàng)式形式

則誤差方程的矩陣形式為

3 粗差處理方法

3.1 閾值法

Huber[2]在1977年采用Tukey污染分布模擬了一些含有不同數(shù)量粗差的正態(tài)分布數(shù)據(jù),其污染正態(tài)分布模式為

其中Φ是標(biāo)準(zhǔn)正態(tài)分布,ε為污染率,是一小量,μ和σ分別為母體均值和標(biāo)準(zhǔn)差。式(9)表明,與均值之差的絕對值大于3σ的觀測值被認(rèn)為是粗差。而實(shí)際觀測中,被測量的真實(shí)值大多是未知,μ和σ不能精確知道,通常以樣本均值和標(biāo)準(zhǔn)差S代替。閾值法的實(shí)用形式為

在實(shí)際操作過程中,把交叉點(diǎn)平差后的不符值的改正數(shù)V作為觀測值向量X,當(dāng)觀測值xi(xi實(shí)際上為交叉點(diǎn)不符值的改正數(shù)vi)滿足式(10)時(shí),被認(rèn)作是粗差。

3.2 抗差M估計(jì)

經(jīng)典最小二乘估計(jì)是利用一組來自于母體為正態(tài)分布的觀測值,求定母體參數(shù)的一種參數(shù)估計(jì),它的最優(yōu)統(tǒng)計(jì)特性(不能脫離觀測值母體)是正態(tài)分布,即觀測值中僅含有偶然誤差這一前提[14]。因此,最小二乘估計(jì)對偏離正態(tài)分布的數(shù)據(jù)是不可靠的。與最小二乘估計(jì)相比,抗差估計(jì)應(yīng)具備兩大特點(diǎn):一是它能消除和削弱粗差對估值的影響;二是它基本上具備經(jīng)典估計(jì)的優(yōu)良特性。

按M估計(jì)原理,取極值函數(shù)為

其中,ai為系數(shù)陣A的第i行。對X求導(dǎo),并令其為0,同時(shí)記φ(Vi)=?ρ/?Vi,則有

令φ(Vi)/Vi=Wi(權(quán)因子)=piiWi為等價(jià)權(quán)元素,則式(12)可寫為

由此,可得到抗差M估值

抗差估計(jì)采用迭代解法,每一次迭代相當(dāng)于進(jìn)行一次最小二乘估計(jì),但應(yīng)以等價(jià)權(quán)代替先驗(yàn)權(quán)。式(7)的第k+1步迭代解為

抗差估計(jì)的結(jié)果取決于所選取的等價(jià)權(quán)函數(shù),權(quán)函數(shù)不一樣,抗差估計(jì)的模型就不同,抗差的效果也不一樣[15]。常見的等價(jià)權(quán)模型有Huber權(quán)函數(shù)、Hampel權(quán)函數(shù)、Tukey權(quán)函數(shù)和 Andrews權(quán)函數(shù)等[9]。本文僅以Huber等價(jià)權(quán)模型為例進(jìn)行了抗差估計(jì),其具體形式為:

式中c一般取2.0。

4 算例及結(jié)果分析

本文模擬了某航空重力測區(qū),它由10條橫向測線和10條縱向測線組成,飛機(jī)起伏飛行,飛行時(shí)離地面高度約為1.8 km。應(yīng)用兩步處理法對觀測數(shù)據(jù)進(jìn)行處理,第一階段交叉點(diǎn)平差后主、副測線驗(yàn)后單位權(quán)中誤差見表1~2。為了更好地檢驗(yàn)閾值法和抗差M估計(jì)的實(shí)際效果,在這些交叉點(diǎn)不符值的改正數(shù)中隨機(jī)挑選10個(gè)加入絕對值分別為(10~18)×10-5ms-2和(6~10)×10-5ms-2的粗差,形成兩組觀測值data1和data2,其交叉點(diǎn)的分布如圖1。在兩步處理法第二階段中分別采用閾值法和抗差估計(jì)來消除粗差的影響,系統(tǒng)誤差模型為2階和3階一般多項(xiàng)式,表3~9給出了相應(yīng)的最小二乘估計(jì)(LSE)和抗差估計(jì)(RE)結(jié)果。

圖1 交叉點(diǎn)的分布Fig.1 Distribution of crosspoints

由表1可見,用閾值法進(jìn)行粗差的剔除后,數(shù)據(jù)的標(biāo)準(zhǔn)偏差變小,測線網(wǎng)精度有所提高。這說明閾值法能有效對粗差進(jìn)行識別和定位,是一種可靠的粗差探測手段,但它的缺點(diǎn)也很明顯,其對絕對值較小的粗差不敏感,且在剔除粗差觀測值之后,會導(dǎo)致交叉點(diǎn)個(gè)數(shù)變少,這將給后續(xù)研究帶來不利影響。

從表4~7可以看出,當(dāng)部分觀測值含有粗差時(shí),主、副測線的最小二乘估計(jì)驗(yàn)后單位權(quán)中誤差均大于抗差估計(jì)驗(yàn)后單位權(quán)中誤差,這說明當(dāng)觀測值含有粗差時(shí),從整體精度上,采用抗差估計(jì)進(jìn)行系統(tǒng)誤差補(bǔ)償?shù)某晒黠@優(yōu)于采用最小二乘估計(jì)進(jìn)行系統(tǒng)誤差補(bǔ)償?shù)某晒?。但我們也不能單純的以交叉點(diǎn)不符值的均方差大小作為衡量抗差效果唯一標(biāo)準(zhǔn)。當(dāng)信噪比一定時(shí),計(jì)算結(jié)果很大程度上取決于誤差模型的合理選擇,從本算例來看,二階多項(xiàng)式的抗差效果更好。反之,當(dāng)系統(tǒng)誤差模型一定時(shí),信噪比越低,觀測值含粗差的絕對值較大時(shí),抗差效果越明顯。

表1 不含粗差時(shí)主測線系統(tǒng)誤差模型驗(yàn)后單位權(quán)中誤差(單位:10-5ms-2)Tab.1 Posteriori unit weight mean errors of system error model on main line of data without outliers(unit:10-5ms-2)

表2 不含粗差時(shí)副測線系統(tǒng)誤差模型驗(yàn)后單位權(quán)中誤差(單位:10-5ms-2)Tab.2 Posteriori unit weight mean errors of system error model on vice line of data without outliers(unit:10-5ms-2)

表3 閾值法處理結(jié)果(單位:10-5ms-2)Tab.3 Results of Thresholding method(unit:10-5ms-2)

表4 data1主測線系統(tǒng)誤差模型參數(shù)驗(yàn)后單位權(quán)中誤差(單位:10-5ms-2)Tab.4 Posteriori unit weight mean errors of system error model on main line of data1(unit:10-5ms-2)

表5 data1副測線系統(tǒng)誤差模型參數(shù)驗(yàn)后單位權(quán)中誤差(單位:10-5ms-2)Tab.5 Posteriori unit weight mean errors of system error model on vice line of data1(unit:10-5ms-2)

表6 data2主測線系統(tǒng)誤差模型參數(shù)驗(yàn)后單位權(quán)中誤差(單位:10-5ms-2)Tab.6 Posteriori unit weight mean errors of system error model on main line of data2(unit:10-5ms-2)

表7 data2副測線系統(tǒng)誤差模型參數(shù)驗(yàn)后單位權(quán)中誤差(單位:10-5ms-2)Tab.7 Posteriori unit weight mean errors of system error model on deputy line of data2(unit:10-5ms-2)

表8 data1交叉點(diǎn)最小二乘殘差與抗差估計(jì)殘差(單位:10-5ms-2)Tab.8 Residuals of least squares estimation and robust estimation of data1(unit:10-5ms-2)

表9 data2交叉點(diǎn)最小二乘殘差與抗差估計(jì)殘差(單位:10-5ms-2)Tab.9 Residuals of least squares estimation and robust estimation of data2(unit:10-5ms-2)

表8~9列出了不含粗差時(shí)最小二乘殘差、加粗差后最小二乘殘差以及加粗差后抗差M估計(jì)的殘差(系統(tǒng)誤差模型為二階多項(xiàng)式)。不難發(fā)現(xiàn),加粗差后抗差M估計(jì)的殘差比最小二乘估計(jì)的殘差更接近于粗差值本身,在異常觀測值上抗差M估計(jì)具有明顯的抵御粗差干擾的能力,這是因?yàn)樵诘?jì)算時(shí),抗差M估計(jì)利用等價(jià)權(quán)模型對異常觀測值實(shí)施了降權(quán)操作。通過對殘差進(jìn)行分析對比,抗差估計(jì)不僅能對粗差的位置進(jìn)行識別,而且還能對粗差的大小進(jìn)行估計(jì)。上述結(jié)果還表明,當(dāng)粗差絕對值較小時(shí),抗差估計(jì)仍然能夠?qū)ζ溥M(jìn)行識別和定位。從這方面來講,抗差估計(jì)比閾值法更加有效和可靠。

5 結(jié)語

對于各測線上的所有采樣點(diǎn)的系統(tǒng)誤差補(bǔ)償來說,閾值法和抗差估計(jì)均具有一定的抗差效果。閾值法是直接剔除異常觀測值,而抗差M估計(jì)能在拒絕和接受一個(gè)觀測值之間起一個(gè)平滑的作用。相比較而言,抗差估計(jì)可靠性更高,更加有效。值得注意的是,任何一種抗差估計(jì)方法的實(shí)際抗差能力都是有一定條件的,當(dāng)粗差個(gè)數(shù)過多、絕對值較小時(shí),閾值法和抗差估計(jì)對粗差的處理結(jié)果往往不夠理想。后續(xù)研究將考慮更加復(fù)雜的系統(tǒng)誤差模型(如三角多項(xiàng)式)和更為豐富的權(quán)函數(shù)模型以及對系統(tǒng)誤差模型的參數(shù)進(jìn)行顯著性檢驗(yàn)。

1 楊元喜.抗差估計(jì)理論及其應(yīng)用[M].北京:八一出版社,1993.(Yang Yuanxi.Theory of robust estimation and its application[M].Beijing:Bayi Press,1993)

2 黃幼才.數(shù)據(jù)探測與抗差估計(jì)[M].北京:測繪出版社,1990.(Huang Youcai.Data detection and robust estimation[M].Beijing:Mapping Press,1993)

3 Kern M,et al.Outlier detection algorithms and their performance in GOCE gravity field Processing[J].Journal of Geodesy,2005,78(9):509-519.

4 Hwang C,Wang C G and Lee L H.Adjustment of relative gravity measurements using weighted and datum-free constraints[J].Computeramp;Geoscience,2002,28:1 005-1 015.

5 李德仁.誤差處理和可靠性理論[M].北京:測繪出版社,1988.(Li Deren.Error processing and theory of reliablity[M].Beijing:Mapping Press,1988)

6 丁士俊,陶本藻.半?yún)?shù)模型的統(tǒng)計(jì)診斷量與粗差檢驗(yàn)[J].大地測量與地球動力學(xué),2005,(3):24-28.(Ding Shijun and Tao Benzao.Measures of statistical diagnosis and gross test[J].Journal of Geodesy and Geodynamics,2005,(3):24-28)

7 Wang Y.Jump and sharp cusp detection by wavelets[J].Biometrika,1995,82(2):385-397.

8 楊茂興.小樣本容量測量數(shù)據(jù)中粗差的剔除[J].計(jì)量與測試技術(shù),2005,(1):27-28.(Yang Maoxing.The coarse error eliminating in small sample space[J].Metrology and Measurement Technique,2005(1):27-28)

9 楊元喜.參數(shù)平差模型的抗差最小二乘解[J].測繪通報(bào),1994,(6):33-35.(Yang Yuanxi.Robust estimation solutions of parameter adjustment model[J].Bulletin of Surveying and Mapping,1994,(6):33-35)

10 周江文.經(jīng)典誤差理論與抗差估計(jì)[J].測繪學(xué)報(bào),1989,18(2):115-120.(Zhou Jiangwen.Classical theory of error and robust estimation[J].Acta Geodaetica et Carto graphica Sinica,1989,18(2):115-120)

11 黃謨濤,等.海洋重力測量系統(tǒng)誤差補(bǔ)償?shù)膬刹教幚矸ǎ跩].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2002,27(3):251 -255.(Huang Motao,et al.Two-step processing for compensating the systemic errors in marine gravity measurements[J].Geomatics and Information Science of Wuhan University,2002,27(3):251-255)

12 李海.航空重力測量測線網(wǎng)平差的理論與方法[D].解放軍信息工程大學(xué),2002.(Li Hai.Theory and methods of network adjustment of airborne gravimetry[D].The PLA Information Engineering University,2002)

13 蔡劭琨.航空重力測量網(wǎng)絡(luò)平差方法研究[D].國防科學(xué)技術(shù)大學(xué),2009.(Cai Shaokun.Research on the network adjustment of airborne gravimetry[D].National University of Defense Technology,2009)

14 余學(xué)詳,呂偉才.抗差估計(jì)在粗差探測和平差計(jì)算中的運(yùn)用[J].測繪工程,1998,8(3):40-44.(Yu Xuexiang and Lv Weicai.The application of robust estimation to the gross error detection and adjustment caculation[J].Engineering of Surveying and Mapping,1998,8(3):40-44)

15 余學(xué)祥,呂偉才.基于標(biāo)準(zhǔn)化的相關(guān)觀測抗差估計(jì)模型[J].武漢測繪科技大學(xué)學(xué)報(bào),1999,24(1):72-78.(Yu Xuexiang and Lü Weicai.Robust estimation model for correlated observations based on standardized residuals[J].Journal of Wuhan Technical University of Surveying and Mapping,1999,24(1):72-78)

PROCESSING METHODS FOR OUTLIERS IN NETWORK ADJUSTMENT OF AIRBORNE GRAVIMETRY

Zhou Boyang1),Luo Zhicai1,2),Lin Xu1)and Zhou Hao1)

(1)School of Geodesy and Geomatics,Wuhan University,Wuhan 430079
2)Key Laboratory of Geospace Environment and Geodesy,Wuhan University,Wuhan 430079)

The main task of network adjustment of airborne gravimetry is to compensate systemic errors.In order to weaken the influence of outliers and improve the effect of compensation,this paper brings in thresholding and robust estimation,and confirms these two methods’validity and reliability in simulated data processing.By means of comparison,a conclusion is reached that the robust estimation has more practical value in the network adjustment of airborne gravimetry.

airborne gravimetry;network adjustment;thresholding;robust estimation;weight function

1671-5942(2012)02-0110-05

2011-12-25

國家自然科學(xué)基金(41174062);中國博士后科學(xué)基金(20110491189);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金(111110)

周波陽,1983年生,博士研究生,目前從事航空重力測量數(shù)據(jù)處理研究.E-mail:byzhou@whu.edu.cn

P207

A

猜你喜歡
抗差交叉點(diǎn)線網(wǎng)
新型線網(wǎng)城軌乘客信息系統(tǒng)的研究與分析
圍棋棋盤的交叉點(diǎn)
軌道交通COCC線網(wǎng)信號系統(tǒng)設(shè)計(jì)
改善單頻PPP參數(shù)收斂速度的抗差估計(jì)方法
基于高中生命科學(xué)知識交叉點(diǎn)的教學(xué)方法研究
地形簡化對DEM不確定性的抗差性研究
基于抗差最小均方估計(jì)的輸電線路參數(shù)辨識
區(qū)域重力異常值的交叉點(diǎn)平差實(shí)例分析
緊湊型大都市區(qū)軌道線網(wǎng)形態(tài)配置研究
自動售檢票線網(wǎng)化維修管理系統(tǒng)的構(gòu)建
咸丰县| 桐庐县| 治县。| 克拉玛依市| 乌拉特前旗| 辽阳县| 色达县| 九龙坡区| 晋中市| 广丰县| 前郭尔| 崇仁县| 遵义市| 湄潭县| 大方县| 毕节市| 额敏县| 高清| 葫芦岛市| 枝江市| 咸宁市| 太康县| 兴文县| 常宁市| 安庆市| 中牟县| 洪洞县| 潼关县| 漠河县| 图片| 云和县| 宁夏| 丁青县| 云阳县| 额济纳旗| 明光市| 澄江县| 湖州市| 嘉鱼县| 池州市| 伊春市|