劉 碩，劉光博，劉尚國

(1.山東科技大學(xué) 測繪與空間信息學(xué)院，山東 青島 266590；2.通用技術(shù)集團(tuán)工程設(shè)計有限公司，山東 濟(jì)南 250031)

在工業(yè)測量領(lǐng)域，常面臨觀測目標(biāo)尺寸大、精度要求高、測量現(xiàn)場環(huán)境復(fù)雜等一系列問題，對高精度控制網(wǎng)的構(gòu)建提出了更高的要求[1-3]。激光跟蹤儀采用激光干涉測距，距離測量的精度極高，且測量范圍大，自動化程度高，操作簡便，已成為目前構(gòu)建精密三維控制網(wǎng)的重要手段[4-5]。然而，為了保證其動態(tài)測量速度的需要，激光跟蹤儀測角系統(tǒng)一般采用光柵度盤，很難消除度盤誤差對其測角誤差的影響[6]，而測角誤差是影響其點位誤差的主要因素[7]。因此，相比于微米級測距精度，激光跟蹤儀的測角精度并沒有顯著提高，鑒于激光跟蹤儀的測距精度對空間點位精度起主導(dǎo)作用[8]，為消除測角誤差的影響，在構(gòu)建精密三維網(wǎng)時通常僅使用激光跟蹤儀的高精度測距結(jié)果。而在具體平差計算時，必須解決自由網(wǎng)平差基準(zhǔn)秩虧的問題。文獻(xiàn)[9]采用混合總體最小二乘迭代算法求解多測站空間坐標(biāo)轉(zhuǎn)換參數(shù)，并分析了測站定位與定向?qū)す飧檭x光束法平差的影響。文獻(xiàn)[10]則針對激光跟蹤三維測邊控制網(wǎng)中的秩虧問題，提出利用矩陣奇異值分解性質(zhì)求解未知參數(shù)，實現(xiàn)對激光跟蹤儀測角誤差的有效約束，方法實用性強(qiáng)且易于編程實現(xiàn)，可解決自由網(wǎng)平差基準(zhǔn)秩虧的問題，但在附加重心基準(zhǔn)約束及其平差結(jié)果校驗等方面，尚需開展進(jìn)一步探討。

本研究以Leica AT930型激光跟蹤儀為研究對象，通過激光跟蹤儀多站位測量，構(gòu)建了高精度三維測邊網(wǎng)。利用公共點轉(zhuǎn)換模型，求取激光跟蹤儀各測站點在基準(zhǔn)坐標(biāo)系下的近似坐標(biāo)，再利用測邊網(wǎng)平差前后重心位置不變的性質(zhì)，引入重心基準(zhǔn)作為附加約束條件；并利用該性質(zhì)作為測邊網(wǎng)平差后的檢核條件，以解決測邊網(wǎng)中秩虧自由網(wǎng)平差的問題。

1 激光跟蹤三維測邊網(wǎng)的構(gòu)建

第i個測站對第j個全局控制點的激光干涉測距值

(1)

對式(1)進(jìn)行線性化，可得誤差方程式：

vij=aijδXi+bijδYi+cijδZi-aijδxj-bijδyj-cijδzj-lij，

(2)

式中，系數(shù)aij、bij、cij與常數(shù)lij分別如下所示。

其中，lij表示的是第i個測站與第j個控制點間激光干涉測距值與近似值的差值。

對測量空間內(nèi)m個測站、n個控制點應(yīng)用式(1)，建立誤差方程并寫成矩陣形式，則有

Vmn×1=Amn×3(m+n)δX3(m+n)×1-lmn×1。

(3)

式中：V為殘差向量，δX為待求未知參數(shù)向量，A為系數(shù)矩陣，l為距離差值常數(shù)向量，分別如下所示。

V=[v11,v12,…,v1n,…,vm1,vm2,…,vmn]T，l=[l11,l12,…,l1n,…,lm1,lm2,…,lmn]T，

δX=[δX1,δY1,δZ1,…,δXm,δYm,δZm,…,δx1,δy1,δz1,…,δxn,δyn,δzn]T，

2 平差前測邊網(wǎng)近似三維坐標(biāo)的解算

在利用重心基準(zhǔn)的激光跟蹤三維測邊網(wǎng)平差中，測站點與全局控制點的近似三維坐標(biāo)計算是平差處理的首要問題。全局控制點的近似三維坐標(biāo)，一般取自基準(zhǔn)坐標(biāo)系(假定第1測站坐標(biāo)系為基準(zhǔn)坐標(biāo)系)下的單測站觀測成果。而關(guān)于各測站近似坐標(biāo)的處理，本研究首先考慮利用空間距離后方交會的方法進(jìn)行解算。如圖1所示，假設(shè)PA、PB、PC三個控制點在基準(zhǔn)坐標(biāo)系下的坐標(biāo)已知，PE點為未知點。在PE點處設(shè)置激光跟蹤儀測站，依次觀測PA、PB、PC控制點，獲取3個激光干涉測距值D1、D2、D3。

圖1 距離后方交會

由空間距離的計算公式，可得

(4)

通過解算上述三元二次方程，即可求得未知測站點PE的三維坐標(biāo)。但使用后方交會方法獲取測站點近似坐標(biāo)時，解的精確與否與三點交會的網(wǎng)形有很大關(guān)系，若網(wǎng)形不佳如PA、PB、PC三點接近共線或共圓，會造成解算結(jié)果的失真，影響下一步測邊網(wǎng)的平差工作。

因此，考慮利用公共點轉(zhuǎn)換模型求取各測站的近似坐標(biāo)。即在某測量空間選取不少于3個目標(biāo)點，并保證所選目標(biāo)點在每一個測站下都得到了觀測。選擇某測站坐標(biāo)系為基準(zhǔn)坐標(biāo)系，將每一測站下所選目標(biāo)點轉(zhuǎn)換至基準(zhǔn)坐標(biāo)系，求取坐標(biāo)轉(zhuǎn)換參數(shù)，然后導(dǎo)入相應(yīng)測站坐標(biāo)，即可求得其在基準(zhǔn)坐標(biāo)系下的坐標(biāo)。將求得的各測站點近似坐標(biāo)、基準(zhǔn)坐標(biāo)系下的單測站觀測成果作為三維測邊網(wǎng)平差前的初始值。

3 基于重心基準(zhǔn)的三維測邊網(wǎng)平差模型

對式(3)應(yīng)用最小二乘原理，可得法方程式：

NδX-W=0，

(5)

式中，法方程系數(shù)矩陣N=ATPA，常數(shù)陣W=ATPl，P為權(quán)矩陣，參考文獻(xiàn)[11]，P可由式(6)確定

Pij=(1/mSSij)2。

(6)

式中mS為儀器出廠時的標(biāo)稱測距精度。

在激光跟蹤三維測邊控制網(wǎng)中，所有測站點、控制點三維坐標(biāo)作為未知參數(shù)參與平差，且基于激光干涉測距值構(gòu)建的激光跟蹤三維測邊網(wǎng)，僅有尺度基準(zhǔn)而缺少位置、方位基準(zhǔn)，因此該三維測邊網(wǎng)的平差解算屬于典型的秩虧自由網(wǎng)平差問題[12]。系數(shù)矩陣A為秩虧矩陣，因此法方程系數(shù)矩陣N為奇異矩陣，式(5)的解不唯一。為了獲得唯一解，需加入新的約束條件。依據(jù)秩虧網(wǎng)平差原理，三維控制網(wǎng)中重心位置是固定不變的。因此，可利用重心基準(zhǔn)條件進(jìn)行約束，其平差模型為：

GTδX=0，

(7)

式中G為附加約束矩陣。激光跟蹤三維測邊網(wǎng)的基準(zhǔn)數(shù)為6[13]，因此，該附加約束矩陣為：

由式(7)可得到：

(8)

(9)

聯(lián)立式(8)與式(9)可得平差后重心坐標(biāo)為：

(10)

從式(10)可以看出，平差后測站點與全局控制點重心坐標(biāo)等于平差前各點的重心坐標(biāo)，該特性可作為基于重心基準(zhǔn)的激光跟蹤三維測邊網(wǎng)平差后的檢核條件。

由秩虧自由網(wǎng)平差公式，可分別得到未知參數(shù)δX及其協(xié)因數(shù)矩陣QδXδX：

δX=(N+GGT)-1ATPl，

(11)

QδXδX=(N+GGT)-1N(N+GGT)-1。

(12)

單位權(quán)中誤差μ和未知參數(shù)精度mδX分別為：

(13)

(14)

4 實驗與算例分析

4.1 實驗測試方案

如圖2所示，在山東科技大學(xué)3D試驗場內(nèi)，選取約10 m×10 m的測量范圍，均勻布設(shè)7個全局控制點，編號為Ph(h=1,2,…,7)，其中P1、P2、P4、P5位于地面上，P3、P6、P7位于墻面上，點間最大高差約為0.9 m。

圖2 測邊控制網(wǎng)及實驗場景

采用Leica AT930型激光跟蹤儀作為測量設(shè)備，共設(shè)置10個測站，編號依次為Stk(k=1,2,…,10)，依次觀測7個全局控制點。觀測前用熱熔膠固定控制點基座，將球棱鏡放置于基座上并調(diào)整好激光入射角度即可進(jìn)行觀測。

4.2 三維坐標(biāo)初始值的計算

假設(shè)各測站點坐標(biāo)為(0,0,0)，選擇St1測站坐標(biāo)系作為基準(zhǔn)坐標(biāo)系，利用公共點轉(zhuǎn)換模型，求取其他各測站點在St1測站下的坐標(biāo)值，將其作為激光跟蹤三維測邊網(wǎng)平差前的初始值，各測站點坐標(biāo)初值與全局控制點坐標(biāo)初值如表1所示。

表1 測站點與全局控制點坐標(biāo)初值

4.3 重心基準(zhǔn)約束下的算法及程序設(shè)計

該激光跟蹤三維測邊網(wǎng)中，有70個觀測方程，51個未知參數(shù)，6個基準(zhǔn)數(shù)，25個多余觀測數(shù)，誤差方程有解。通過軟件編程，實現(xiàn)了附加重心基準(zhǔn)約束條件的秩虧自由網(wǎng)平差解算(圖3)。算法流程如圖4所示，導(dǎo)入測站點、控制點坐標(biāo)初值，并將激光干涉測距值作為真實值參與平差。計算平差前測邊網(wǎng)的重心坐標(biāo)O1與平差后的重心坐標(biāo)O2，用以檢驗該測邊網(wǎng)平差后的準(zhǔn)確性。平差后，單位權(quán)中誤差為5 μm，平差結(jié)果與各未知參數(shù)精度如表2、表3所示。

圖3 平差模型解算界面Fig. 3 Solution interface of adjustment model

圖4 三維測邊網(wǎng)平差算法流程圖

由表2、表3測站點與全局控制點平差結(jié)果，可以反算各測站到全局控制點的距離，其與激光干涉測距值的差值可用以檢驗平差結(jié)果的準(zhǔn)確性。圖5表示平差前與平差后各點坐標(biāo)反算距離與激光干涉測距值的差值。

圖5 平差前后坐標(biāo)反算距離與激光干涉測距值的差值

4.4 實驗分析

1) 對表2、表3依據(jù)重心公式計算可得，平差前該控制網(wǎng)重心坐標(biāo)O1等于平差后點坐標(biāo)計算的重心坐標(biāo)O2，驗證了激光跟蹤三維測邊網(wǎng)平差前后，網(wǎng)重心位置不變的性質(zhì)。

表2 測站點坐標(biāo)平差結(jié)果

2) 測站點與全局控制點平差后，平面精度平均值約為14 μm，而高程精度平均值為62 μm，即平差后點的高程精度略差于平面精度。

3) 由表3可以看出，地面控制點精度平均值為43 μm，墻面控制點精度平均值為115 μm，即地面點精度優(yōu)于墻面點精度，分析原因是部分布設(shè)在墻面上的點位穩(wěn)定性較差，對最終點位測量精度產(chǎn)生了較大影響。

表3 全局控制點坐標(biāo)平差結(jié)果

4) 經(jīng)平差處理后，激光跟蹤儀測站點位誤差最大值為79 μm，最小值為42 μm，平均值為57 μm；而全局控制點點位誤差最大值為160 μm，最小值為35 μm，平均值為74 μm，即測站點精度略高于控制點精度。

5) 由圖5可以看出，利用重心基準(zhǔn)的激光跟蹤三維測邊網(wǎng)平差后，點坐標(biāo)反算距離與激光干涉測距值的差值較平差前更為均勻。平差前該差值最大為26.5 μm，最小為-30.1 μm，平均值為-1.8 μm；平差后該差值最大為10.3 μm，最小為-9.6 μm，平均值僅為0.2 μm，證明該平差模型可對激光跟蹤儀測角誤差進(jìn)行有效的約束。

5 結(jié)論

本研究利用激光跟蹤儀微米級測距精度的優(yōu)勢，在試驗場內(nèi)構(gòu)建了高精度三維測邊網(wǎng)，基于控制網(wǎng)平差前后重心基準(zhǔn)不變的性質(zhì)，較好地解決了三維測邊網(wǎng)中秩虧自由網(wǎng)平差的問題。實驗表明，利用重心基準(zhǔn)的激光跟蹤三維測邊網(wǎng)平差模型，可有效地消除激光跟蹤儀測角誤差對整體測量精度的影響，在10 m×10 m測量范圍內(nèi)，平差后的測站、控制點點位精度平均值優(yōu)于0.1 mm。平差后，點坐標(biāo)反算距離與激光干涉測距值的差值可控制在10 μm以內(nèi)。