測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計方法

王樂洋,李志強

(1.東華理工大學(xué) 測繪工程學(xué)院，江西 南昌 330013；2．山東建筑大學(xué) 測繪地理信息學(xué)院，山東 濟南 250101)

大地測量數(shù)據(jù)處理所涉及到的觀測模型一般為非線性模型[1-2]，而非線性模型參數(shù)估計一直是測量平差所研究的重點問題之一。對于非線性模型，常用固有曲率和參數(shù)效應(yīng)曲率來刻畫非線性平差模型的非線性強度，進而評估參數(shù)估計的統(tǒng)計推斷效果[3]。處理非線性模型的傳統(tǒng)方法是線性最小二乘法，利用泰勒公式將非線性模型展開并截取至一次項[4]。對于非線性強度較高的模型，如三角網(wǎng)、導(dǎo)線網(wǎng)，線性近似將引起較大的模型誤差[5-7]。隨著測量技術(shù)的不斷發(fā)展，對測量及平差的精度要求不斷提高，因此，傳統(tǒng)方法并不能滿足當(dāng)今科學(xué)技術(shù)的要求，而如何減弱模型線性化帶來的模型誤差成為了提高成果質(zhì)量的重要內(nèi)容。關(guān)于如何削弱模型誤差影響，文獻[7]通過對非線性函數(shù)求三階偏導(dǎo)數(shù)，考慮了二次項及三次項的影響，提出了非線性模型參數(shù)估計的直接解法。為了有效地避免導(dǎo)數(shù)計算，可以通過函數(shù)差分的方法，利用快速差分迭代解算模型獲取非線性最小二乘參數(shù)估值[8]。文獻[9]研究了基于非線性最小二乘的空間三維直角坐標(biāo)轉(zhuǎn)換算法。針對測距定位方程該非線性平差模型，文獻[10]給出了距離函數(shù)二階導(dǎo)數(shù)的表達式，并推導(dǎo)了非線性平差的封閉牛頓迭代公式以及退化條件。對于非線性平差精度評定問題，可以從近似非線性函數(shù)的概率密度入手，采用無需求導(dǎo)的Sterling插值方法研究非線性函數(shù)的統(tǒng)計信息[11]。

目前，解決非線性最小二乘參數(shù)估計問題的方法主要分三類[12]：第一類是線性最小二乘法，用線性模型的理論與方法近似求解，其弊端在于當(dāng)模型的非線性強度較強時，會產(chǎn)生較大的模型誤差[6]。第二類是直接搜索算法，如模擬退火法、單純形法、遺傳算法等，其優(yōu)勢在于不需要求導(dǎo)計算，但是該類方法無法獲得參數(shù)的解析解，且計算耗時很長[13]。第三類方法是迭代解法，如牛頓法、高斯牛頓法、改進的高斯牛頓法等[14]。牛頓法每次迭代時都需要計算目標(biāo)函數(shù)的海森矩陣，當(dāng)遇到非線性函數(shù)復(fù)雜時往往很難獲得[12]。改進的高斯牛頓法與高斯牛頓法具有相同的迭代速率，但由于改進的高斯牛頓法在每次迭代時都需要確定一個步長因子，因此其計算量仍大于高斯牛頓法[13]。高斯牛頓法不僅具有牛頓法的收斂速率，而且每次迭代計算工作量相對最少[12]，且在實際執(zhí)行該算法前，總是通過某種方法求得一個靠近參數(shù)真值的解作為迭代初值進行解算，因此可以避免迭代算法受初值較差的影響，使其迭代收斂。

考慮到使用高斯牛頓迭代法解算非線性模型得到的參數(shù)估值并不能滿足方差最小性和無偏性[15-17]，考慮將Bootstrap采樣方法引入非線性參數(shù)估計問題中。Bootstrap采樣方法也稱為自助法，該方法通過對觀測數(shù)據(jù)的有放回隨機采樣，能夠充分利用觀測值的先驗信息和數(shù)據(jù)性質(zhì)，降低參數(shù)估值偏差，改善參數(shù)估值的質(zhì)量。文獻[18]首次提出了自助法，并給出了采樣策略和相關(guān)證明，提出通過研究由重復(fù)采樣得到的自助世界的經(jīng)驗分布來逼近真實分布的思想。文獻[19]研究了確定置信區(qū)間的Bootstrap改進方法。針對結(jié)構(gòu)復(fù)雜的數(shù)據(jù)，可以采用更具一般性的Bootstrap假設(shè)檢驗方法對統(tǒng)計量進行推斷分析[20]。文獻[21]綜述了自助法在機器學(xué)習(xí)方面的應(yīng)用和進展。針對自助法在標(biāo)準(zhǔn)差估計、偏差計算、回歸分析、區(qū)間估計、置換檢驗、交叉驗證等諸多方面的研究，也獲得了豐富的成果[22-25]。自助法的研究主要集中在該方法在統(tǒng)計學(xué)上的性質(zhì)，目前尚未發(fā)現(xiàn)將該方法用于非線性參數(shù)估計的論文出現(xiàn)。

本文從處理測角網(wǎng)坐標(biāo)平差模型的傳統(tǒng)方法出發(fā)，針對該方法在處理非線性模型時引起較大模型誤差的局限性，采用高斯牛頓迭代方法求取未知點坐標(biāo)；為減小迭代計算產(chǎn)生的偏差，引入Bootstrap采樣方法，從而進一步改善參數(shù)估值質(zhì)量。針對等精度和不等精度觀測數(shù)據(jù)，本文通過兩個案例研究并與傳統(tǒng)方法進行對比分析，來驗證本文方法的有效性和優(yōu)勢性。

1 測角網(wǎng)坐標(biāo)平差模型及其高斯牛頓迭代解法

以角度為觀測值、未知點坐標(biāo)為參數(shù)的平差模型，稱為測角網(wǎng)坐標(biāo)平差模型。為直觀表示，假設(shè)角度測量的示意圖如圖1所示，其中點j為頂點，∠hjk為待測角。

圖1 測角示意圖

假設(shè)測角觀測方程的向量表達式為：

L=f(X)+ε.

(1)

式中：L∈Rn×1為相互獨立的角度觀測值向量；X∈Rt×1為未知點坐標(biāo)向量；f(·)表示非線性映射關(guān)系，f(X)為n個非線性函數(shù)組成的向量；ε∈Rt×1表示角度觀測誤差向量。

用未知點坐標(biāo)估值和測角真誤差的估值分別代替其真值，得測角網(wǎng)模型的誤差方程：

(2)

(3)

(4)

矩陣表達形式為：

v=Bx-l.

(5)

式中：B為誤差方程的系數(shù)矩陣，x為未知點坐標(biāo)改正數(shù)向量，l為誤差方程常數(shù)項向量。

根據(jù)最小二乘原理，加權(quán)最小二乘解為[26-27]：

x=(BTPB)-1BTPl.

(6)

式中：P是角度觀測值的權(quán)陣。

最后可得測角網(wǎng)模型未知點坐標(biāo)估值為：

(7)

式中：X0為未知點坐標(biāo)初值。

(8)

(9)

將高斯牛頓解法用于求解測角網(wǎng)坐標(biāo)平差模型的關(guān)鍵，是每步迭代后都需要及時更新系數(shù)矩陣B及常數(shù)項向量l。因此，分析歸納誤差矩陣及常數(shù)項向量的特點和性質(zhì)就顯得尤為重要。

假設(shè)式(4)的系數(shù)為：

(10)

等價于

(11)

式中：ai，bi，ci，di由式(10)求得，其數(shù)值的正負(fù)取決于坐標(biāo)增量ΔX,ΔY的正負(fù)；(+ai),(-ai),(+bi),(-bi),(+ci),(-ci),(+di),(-di)的正負(fù)由坐標(biāo)增量與式(6)中的正負(fù)號共同決定。

由圖1可知，對于所有觀測角度L，其誤差方程總共存在7種表達形式，見表1。

表1 誤差方程的不同形式

根據(jù)高斯牛頓迭代原理及系數(shù)矩陣及常數(shù)項向量的更新規(guī)律，測角網(wǎng)坐標(biāo)平差的高斯牛頓迭代解法總結(jié)為算法1。

算法1：高斯牛頓迭代解法。

1)利用余切公式獲取迭代初值X0；

2)計算各基線邊的坐標(biāo)增量ΔX，ΔY及邊長S；

3)由式(10)計算系數(shù)ai，bi，ci，di，并獲取Bk；

4)計算各基線邊的近似坐標(biāo)方位角α0，得到lk；

2 測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計方法

Bootstrap采樣方法通過重采樣給定觀測數(shù)據(jù)獲取自助樣本來統(tǒng)計推斷未知參數(shù)，該方法通過對已知樣本(假設(shè)樣本容量為n)做有放回隨機循環(huán)抽樣，得到M個樣本容量仍為n的自助樣本，進而估計原始樣本的有效信息。自助法能充分利用總體包含在原始樣本中的所有信息，使得自助樣本的分布能夠不斷逼近原始樣本統(tǒng)計量的分布。該方法的優(yōu)勢性在于它不需要對未知模型及分布做任何假設(shè)，也無需推導(dǎo)估計量的精確表達式，只需通過檢驗樣本內(nèi)統(tǒng)計量的變化來估計未知參數(shù)的整個抽樣分布，即對觀測值進行重采樣并計算估計值[21]。

(12)

根據(jù)Bootstrap方法采樣策略、高斯牛頓迭代原理及系數(shù)矩陣和常數(shù)項向量的更新規(guī)律，測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計方法總結(jié)為算法2。

算法2：Bootstrap重采樣方法。

1)假設(shè)原始觀測值樣本N=(L1,L2,…，Ln)為觀測總體，其權(quán)陣為P=(P1,P2,…,P3)；

2)產(chǎn)生n個隨機數(shù)(i1,i2,…，in)，其中1≤iu≤n，(u=1,2,3,…，n)；

3)對N中的觀測值進行采樣，獲取自助樣本Nr=(Li1,Li2,Li3,…,Lin)，其中1≤r≤1 000；

5)根據(jù)自助樣本中的角度數(shù)據(jù)，利用余切公式獲取坐標(biāo)迭代初值X0；

6)計算采樣到的基線邊的坐標(biāo)增量ΔY,ΔX及邊長S；

7)由式(10)計算系數(shù)ai，bi，ci，di，并獲取系數(shù)矩陣Bk；

8)計算采樣到的基線邊的近似坐標(biāo)方位角α0，并獲取常數(shù)項向量lk；

12)通過式(12)計算未知點坐標(biāo)的自助估計值。

根據(jù)以上計算步驟，可以得到測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計算法迭代流程圖，如圖2所示。

圖2 測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計算法迭代流程圖

從上述迭代步驟可以看出，所有自助樣本中包含的樣本數(shù)據(jù)均來源于原始樣本，并未根據(jù)更多的觀測信息進行計算，所以測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計方法與高斯牛頓算法均是利用相同的角度觀測數(shù)據(jù)最終得到未知點坐標(biāo)估值。

3 算例及分析

3.1 案例1：獨立等精度情況

模擬一個加密控制網(wǎng)，如圖3所示。網(wǎng)中A，B，C，D是已知三角點，P是待定點，獨立等精度獲取12個觀測角度，已知點的坐標(biāo)數(shù)據(jù)見表2。

圖3 測角網(wǎng)示意圖

表2 已知點數(shù)據(jù) m

表3 自助法計算參數(shù)估值的過程

從表3可以看出，自助法采樣過程將原始觀測數(shù)據(jù)采樣成了多個自助樣本，雖然每個自助樣本的樣本容量與原始樣本相同，但并不是改變觀測值的先后排列順序。有放回隨機抽樣過程使得每個自助樣本中可能存在重復(fù)的原始數(shù)據(jù)點，而另外一些原始樣本點沒有出現(xiàn)。因此，每個自助樣本將隨機地異于原始樣本，導(dǎo)致每個自助樣本獲得的參數(shù)估值存在細(xì)微差異。

采用最小二乘法(LS)、基于最小二乘的Bootstrap方法(LSB)、算法1(GNI)及算法2(GNIB)計算該測角網(wǎng)中未知點坐標(biāo)，各方法的參數(shù)估值、范數(shù)結(jié)果見表4。LSB法的具體表現(xiàn)為在獲取自助樣本后、通過LS方法來解算自助樣本獲取未知點坐標(biāo)參數(shù)。

表4 參數(shù)估計結(jié)果

從表4結(jié)果可以看出，經(jīng)典LS法得到的坐標(biāo)估值與參數(shù)真值偏離較大?？紤]到參數(shù)估值的真誤差一般由觀測誤差以及線性近似所引起的模型誤差共同影響。而針對非線性強度很強的測角網(wǎng)模型，在具有相同觀測誤差的角度觀測值情況下，LS法得到的參數(shù)估值最差，表明LS法在線性近似的過程中，引入了較大的模型誤差，從而嚴(yán)重影響坐標(biāo)估計值的質(zhì)量。相比于LS法及LSB法，基于逐步線性化的GNI法得到的參數(shù)估值更接近于參數(shù)真值，表明GNI法通過反復(fù)迭代的確能夠削弱線性近似所帶來的模型誤差的影響，使其參數(shù)估值精確度比傳統(tǒng)方法有較為明顯的提高。

同時，相比于傳統(tǒng)的LS法，LSB法能夠改善參數(shù)估值的質(zhì)量；而相比于GNI法，GNIB法在估計值的精確度方面也有較明顯的提升，且均優(yōu)于LS法及LSB法。表明將自助法與傳統(tǒng)的最小二乘法及高斯牛頓迭代法結(jié)合是可行且有效的。自助法的循環(huán)有放回隨機采樣過程獲得的自助樣本能夠充分利用角度觀測值先驗信息，而且該采樣過程不會產(chǎn)生額外的模型誤差，也不會改變模型態(tài)性。尤其是GNIB法，其大量的自助樣本能夠減小GNI法計算坐標(biāo)估值產(chǎn)生的偏差，從而使得GNIB法獲得的坐標(biāo)估值最接近真值。從實驗結(jié)果可以看出，算法2除具備算法1的優(yōu)點外，其采樣角度觀測值獲得的大量自助樣本可以減小GNI法計算產(chǎn)生的偏差，進而改善參數(shù)估值的質(zhì)量。

3.2 案例2：獨立不等精度情況

模擬一個更復(fù)雜的測角網(wǎng)，如圖4所示。網(wǎng)中A,B,C,D是已知點，P1,P2是待定點，獨立不等精度觀測了18個角度，已知點坐標(biāo)值見表5。

圖4 測角網(wǎng)示意圖

表5 已知點數(shù)據(jù) m

表6 角度觀測數(shù)據(jù)

針對不等精度的測角網(wǎng)觀測數(shù)據(jù)，自助法獲取自助樣本及參數(shù)估值的具體實施過程見表7。

表7 Bootstrap參數(shù)估計方法重采樣過程

XLFXL從表7的自助采樣過程可以看出，針對不等精度觀測數(shù)據(jù)，自助法重采樣獲取自助樣本的過程中，需要對觀測值及其權(quán)值同步采樣，保持觀測值與其權(quán)值相對應(yīng)。自助法通過有放回隨機抽樣將樣本容量為18的原始角度觀測值樣本采樣成M個自助樣本，該過程并不是改變觀測值的先后排列順序，其采樣過程使得某些原始數(shù)據(jù)點重復(fù)出現(xiàn)，而另外一些原始樣本點沒有被采樣到。因此，利用高斯牛頓解算每個自助樣本得到的未知點坐標(biāo)估值存在細(xì)微差異。

分別采用線性加權(quán)最小二乘法(WLS)、基于加權(quán)最小二乘的Bootstrap方法(WLSB)、算法1(GNI)以及算法2(GNIB)對該測角網(wǎng)進行參數(shù)估計，各方法的參數(shù)估值、范數(shù)結(jié)果見表8。為更好地實現(xiàn)實驗對照，WLSB法是WLS法與Bootstrap方法結(jié)合后的方法，具體表現(xiàn)為通過WLS法來解算自助樣本。

表8 參數(shù)估計結(jié)果

由表8的結(jié)果可以看出，傳統(tǒng)WLS法獲取的坐標(biāo)估值質(zhì)量最差，是由于該方法需要利用泰勒級數(shù)對測角網(wǎng)模型進行線性化處理，在截斷取至一次項的過程中不可避免地引入了線性化模型誤差。對于非線性強度很強的測角網(wǎng)模型，線性近似將產(chǎn)生大于觀測誤差的模型誤差，該誤差極易引起參數(shù)估值與真值偏離較大；當(dāng)坐標(biāo)初值計算不準(zhǔn)確時，將嚴(yán)重影響坐標(biāo)估計值的質(zhì)量。而相比于WLS法和WLSB法，GNI法得到的參數(shù)估值最優(yōu)，說明GNI法能夠減弱測角網(wǎng)函數(shù)模型線性近似所帶來的模型誤差的影響，其循環(huán)迭代過程可以不斷修正參數(shù)估值，在一定程度上能夠有效改善傳統(tǒng)WLS法對測角網(wǎng)模型線性化過程中所引起的精度損失，使其參數(shù)估值比傳統(tǒng)的線性近似方法更接近參數(shù)真值。

相比于傳統(tǒng)的WLS法，WLSB能夠提高參數(shù)估值的質(zhì)量；而相比于GNI法，GNIB法在參數(shù)估值的精確度方面也具有一定的提升。結(jié)果表明自助法同樣適用于不等精度的角度觀測數(shù)據(jù)。由于自助法重采樣不需要對模型及分布做任何假設(shè)，也無需推導(dǎo)參數(shù)估值的精確表達式，而有放回隨機采樣過程獲得的自助樣本能夠更好地包含原始觀測值的先驗信息和數(shù)據(jù)性質(zhì)。尤其是GNIB法，其大量的自助樣本能夠減小GNI法計算產(chǎn)生的偏差，提高參數(shù)估值的精確度。實驗驗證了自助法在重采樣不等精度觀測數(shù)據(jù)時需要對觀測值的權(quán)值同步采樣的正確性，也再次驗證了將自助法與高斯牛頓迭代法相結(jié)合、并用于測角網(wǎng)坐標(biāo)平差模型的有效性與優(yōu)勢。

4 結(jié)束語

本文介紹了測角網(wǎng)坐標(biāo)平差模型和Bootstrap重采樣理論，將高斯牛頓迭代解法應(yīng)用于測角網(wǎng)坐標(biāo)平差模型，并在此基礎(chǔ)上結(jié)合自助法，提出了測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計算法，并給出了詳細(xì)的采樣步驟及計算流程圖。針對等精度與不等精度的獨立觀測數(shù)據(jù)，將本文方法用于兩個測角網(wǎng)案例，并與傳統(tǒng)的加權(quán)最小二乘法和基于加權(quán)最小二乘的Bootstrap方法求取的坐標(biāo)估值進行比較。由于測角網(wǎng)模型的非線性強度較強，經(jīng)典的加權(quán)最小二乘法在對測角網(wǎng)模型進行線性化處理過程中引入了較大的模型誤差，導(dǎo)致其坐標(biāo)估值與真值偏離較大。而高斯牛頓法通過迭代削弱了因線性近似帶來的模型誤差的影響，有效改善傳統(tǒng)最小二乘法對測角網(wǎng)模型線性化所引起的參數(shù)精度損失，使其坐標(biāo)估值比傳統(tǒng)方法獲得的結(jié)果更加準(zhǔn)確。測角網(wǎng)坐標(biāo)平差模型的Bootstrap參數(shù)估計算法具備了除高斯牛頓迭代法的所有優(yōu)點外，由于其重采樣過程不僅不會產(chǎn)生額外的模型誤差，而且能夠充分利用角度觀測值性質(zhì)，其大量的自助樣本能夠減小高斯牛頓迭代計算產(chǎn)生的偏差，進一步提高參數(shù)估值的精確度。算例驗證了將高斯牛頓迭代解法應(yīng)用于測角網(wǎng)坐標(biāo)平差模型的必要性與實用性，也驗證了將Bootstrap參數(shù)估計方法用于測角網(wǎng)坐標(biāo)平差模型的有效性與優(yōu)勢性。

在計算效率方面，自助法通過重采樣觀測值得到M個自助樣本，其解算過程盡管會增加計算耗時，但是其原理簡單、適用性強，易于理解和編碼，且重采樣次數(shù)M可調(diào)，它僅需要通過計算機重復(fù)解算自助樣本便可改善未知參數(shù)的質(zhì)量，不失為一種較好的參數(shù)估計方法。本文方法同樣也適用于觀測值相互獨立的非線性模型參數(shù)估計問題，本文僅以測角網(wǎng)坐標(biāo)平差為例。隨著自助法的首次引入測角網(wǎng)坐標(biāo)平差模型中，將拓展自助法在參數(shù)估計中的應(yīng)用，也將是非線性平差理論研究的一個重要補充。

Bootstrap采樣方法依賴于原始觀測數(shù)據(jù)的樣本容量，足夠的觀測數(shù)據(jù)對提高參數(shù)估值的自助估計十分重要。此外，對于相關(guān)觀測數(shù)據(jù)的非線性平差及提高計算效率的問題還需要進一步研究。