基于自然樣條補償最小二乘的高程誤差研究

王守春,汪長余,王松慶,米景躍,安麗麗

(天津市水利科學研究院,天津 300061)

基于自然樣條補償最小二乘的高程誤差研究

王守春,汪長余,王松慶,米景躍,安麗麗

(天津市水利科學研究院,天津 300061)

基于非參數自然樣條補償最小二乘估計方法,把 GPS高程數據與對應的水準數據的差值作為觀測值,按照非參數自然樣條補償最小二乘的估計方法,得到消弱偶然誤差v之后的平差值,并以此去改正 GPS高程,可以取得較好的效果。并以某河道測量采集的數據作為事例,證明該方法的有效性。

自然樣條內插;補償最小二乘;非參數;GPS高程

目前測量數據處理中,主要用到的有參數估計、半參數估計以及非參數估計。對于參數估計和半參數估計,其中參數部分的形式是需要已知的,因此它們具有較好的外延性。而對于一組不知道各參數之間關系的測量數據,參數估計的方法就顯得無能為力,這恰好是非參數估計所能解決的問題。非參數自然樣條補償估計是非參數估計方法中較為常用的一種。本文基于非參數自然樣條補償最小二乘的方法來分析高程點矯正誤差。

1 非參數自然樣條補償最小二乘

1.1 光滑自然樣條內插

對于一組測量數據(Li,ti),i=1,…,n,希望得到一光滑曲線s(t),其內插數據點為(Li,ti),且滿足s(ti)=Li,當前最為常用的方法是自然樣條函數法。自然樣條函數是一分段多項式函數,也即是:對某一區間[a,b]作劃分,Δ:a

而對于區間[a,t1]、[tn,b],根據兩對觀測值[a, L0],[t1,L1]和[tn,Ln],[tn+1,Ln+1],都為線性

函數。

1.2 非參數自然樣條補償最小二乘準則

已知非參數模型

式中:s(t)是光滑曲線,Δi是白噪聲,根據混雜有白噪聲的觀測數據Li,i=1,…,n來構造函數s(t),也成為非參數回歸。如果對曲線s(t)不做任何限制而選擇直線內插或樣條內插,可以使得殘差的平方和為0,但是這樣的擬合并沒有去掉或消弱噪聲,所以這樣的內插不是最佳的。

其中,前一項為殘差平方和,后一項為補償項。α是大于0的平滑參數,在擬合程度和光滑程度之間起平衡作用。

2 數據處理及結果分析

2.1 未加補償項的自然樣條內插

現有樣本數為74的一組觀測值,按區間長為1分為73個子區間,依據1.1中介紹的自然樣條內插的方法對該組數據進行自然樣條內插,得到的s(t)如圖1所示,圖2為直線內插的s(t)。

圖1 自然樣條內插的曲線

圖2 直線內插結果

由于直線內插和自然樣條內插都不能去除噪聲,所以這2種內插方法對分析誤差來說沒有太大意義。

2.2 非參數補償最小二乘

根據1.2中介紹的回歸準則,在Matlab中進行迭代計算選取平衡參數α=0.05,s(t)及v=s(t)-L如圖3、圖4所示。

圖4中的v也即是1.2中提到的白噪聲,這樣就得到了剔除白噪聲后的差值s(t),并且以此去改正每公里斷面上的所有高程數據,從而可以得到較為準確的高程數據。圖 5、圖 6分別是當α=1,α=10時所對應的s(t)。

3 結 論

通過自然樣條補償最小二乘的估計方法,可以對非參數(觀測值)、白噪聲進行估計,以此來改正各觀測值,提高精度。在河道斷面測量中,用該方法可以消弱由點矯正和基準轉換所帶來的誤差。具體實施步驟如下:

1)根據已知水準點,進行三等水準測量,并在河道兩岸每公里處留一個水準點;

2)用任意設站的RTK對每公里處所留的三等水準點進行測量,同時對每一個需要測量的河道斷面所對應的兩堤特征點(堤腳,堤頂,防浪墻等)進行測量,并用測深儀加RTK的方式測量該斷面的水下高程;

3)對RTK所測量的數據進行高程點矯正,比對每公里留點處的三等水準高程與矯正過的 GPS高程,若整體差值較大,則以留點處的2種高程之差作為觀測值,進行非參數自然樣條補償最小二乘估計,剔除或減弱偶然誤差(白噪聲),并以剔除或消弱偶然誤差的觀測值來改正對應斷面上的所有 GPS高程;

4)生成斷面成果。非參數自然樣條補償最小二乘估計方法中,確定光滑因子α仍然是一難點,還有待于進一步研究。

[1]丁士俊.測量數據的建模與半參數估計[D].武漢:武漢大學,2005.

[2]柴根象,洪圣巖.半參數回歸模型[M].合肥:安徽教育出版社,1995.

[3]孫海燕.測量平差函數模型的若干討論[J].武漢大學學報,2003(S):62-64.

[4]丁士俊,陶本藻.自然樣條半參數模型與系統誤差估計[J].武漢大學學報:信息科學版,2004,29(11):964-967.

[5]潘雄.半參數測量模型在 GPS高程擬合中的應用[J].武漢工業學院學報,2004,23(2):99-101.

[6]胡宏昌.半參數模型的估計方法及其應用[D].武漢:武漢大學,2004.

[7]丁士俊,孫振冰,劉星.廣義補償最小二乘估計及其統計性質[J].大地測量與地球動力學,2006,26(4):22-26.

[8]孫海燕,吳云.半參數回歸與模型精化[J].武漢大學學報:信息科學版,2002,27(2):172-174.

Study on the height errors based on the penalized least square of natural spline

WANG Shou-chun,WANG Chang-yu,WANG Song-qing,M IJing-yue,AN Li-li

(Tianjin Hydraulic Research Institute,Tianjin 300061,China)

Based in the penalized least square estimation of nonparametric natural sp line,take the difference of GPS Height and Leveling Height as observed values,and get the adjusted value w hose accidental errorvhave been reduced.Then co rrect the data of GPS Height w ith the adjusted value,thus a better result could be obtained.Finally,this theo ry is p roved effective by the data of a river surveying.

natural sp line interpolation;penalized least square;non-parameter;GPS Height

P224

A

1006-7949(2010)04-0027-03

2009-10-19

王守春(1984-),男,助理工程師,碩士研究生.

[責任編輯:張德福]