基于聯合平差的船臺測量場標定方法研究

豐寅帥，彭 飛，王 中，孟慶旭

(海軍工程大學 艦船與海洋學院，湖北 武漢 430033)

0 引 言

隨著船舶制造行業的不斷發展，船舶制造所能達到的精度越來越高，對船舶制造過程中所能達到的測量精度提出了更高的要求，而傳統的以捶球和標尺為主的接觸式測量方法已經逐漸難以滿足當前船舶制造行業的精度需求，因此，以全站儀為主的非接觸式高精度測量方法已經在造船現場得到了廣泛應用。在目前船舶建造現場的測量作業中，由于船體分段尺寸較大，測量條件復雜[1]，在某一固定站位下無法完整獲取整個船體分段的尺寸和形狀數據[2]，以往常見的方法是利用全站儀進行多級轉站的形式，通過不斷的轉移測量站位來獲取船體分段多個角度的測量數據，但這種方法速度較慢，并且轉站過程本身是帶有無法消除的系統誤差，會造成測量精度下降等問題。因此，為了減少轉站次數，實現船體分段測量數據的統一管理，就需要在船臺區域內構建能夠覆蓋全域的具有統一基準的船臺測量場[3]。

圖1 船臺測量場示意圖Fig. 1 Slipway measurement field

如圖1所示，船臺測量場是以多級轉站測量為基礎，在船臺區域內布置相當數量的公共點，通過多次坐標轉換，將不同站位下所測得的公共點坐標統一到一個基準坐標系下。目前，常見的坐標轉換算法主要有奇異值分解法[4]、四元數法[5]和正交矩陣法[6]，由于坐標轉換算法本身就帶有一定的誤差，且在不同站位所測得的公共點數量可能不同，因此，在將公共點坐標轉換至基準坐標系時會引入一定的參數誤差，同時，受限于全站儀等測量設備的加工精度，在某一站位下對公共點進行測量時也會引入一定的測量誤差，進一步加大了船臺測量場的整體誤差。為了解決此類問題，就需要對船臺測量場進行標定，目前常見的方法是基于最小二乘原理[7-9]，對轉站測量得到的公共點坐標進行平差計算，從而得到滿足誤差最小時的坐標轉換參數。近年來，國內外學者對采用平差計算方法進行了大量的研究，Schaffrin[10-12]利用拉格朗日函數構建了平差模型，并進一步推到得到了加權整體最小二平差算法；Calikins[13]和Meid[14]則是通過對測量現場的不確定度進行研究，給出了一種基于光束平差法的誤差處理方法；何榮基基于間接平差原理，給出了自由設站測量時含約束條件的平差模型。但上述研究大多是針對單次測量或單級轉站所開展的，而對于船臺測量場構建時的多級轉站過程，相關的研究還相對較少。

為了對利用多級轉站方法構建船臺測量場時公共點的最優坐標轉換參數進行求解，本文采用一種基于最小二乘原理和公共點經過多級轉站后返回原站位坐標不變原則的聯合平差算法，并通過構建拉格朗日函數得到平差值的求解方法，實現了船臺測量場公共點位置誤差的最小化處理。通過模擬實驗，對此方法的有效性和精度進行了驗證。

1 誤差方程及其求解方法

1.1 誤差方程

根據誤差理論和測量平差原理[15-16]，在利用多級轉站方法構建船臺測量場時，公共點的殘余誤差方程可以寫為：

在轉站測量中，V為轉站后公共點坐標的誤差矩陣；A為 誤差方程的系數矩陣；為轉站參數誤差矩陣；L為公共點的經過轉站后的坐標和測量坐標的偏差。

當轉站測量時公共點誤差最小時，V取最小值，同時，滿足最小值，式中，為公共點坐標的權值矩陣。

1.2 誤差方程的求解方法

由坐標轉換模型可知，在公共點由第j站轉至第j+1站時，有

式中：j和j+1分別公共點的坐標數據在第j站和第j+1站測量得到；kj，Rj和Tj分別代表公共點坐標從第j站轉至第j+1站時的尺度參數、角度參數和距離參數。其中，公共點在轉站前后相對位置不變，故kj取值為1，坐標轉換原理和歐拉角理論，角度參數可表示為Rj=RXj·RYj·RZj，對應角度分別為 αj，βj， γj，距離參數可表示為：

在船臺測量場中，對于某一個公共點在j+1站的坐標，其位置誤差可以表示為其真實坐標和測量值的差值，即

對式（4）作泰勒展開，并與式（1）聯立，可以得到：

式中，為第i個公共點在第j站的測量坐標初值，分別為第i個公共點從第j站轉至第j+1站的平移參數，分別代表第i個公共點在第j+1站的真實坐標初值。

2 約束條件方程及其求解方法

2.1 約束條件方程

在轉站過程中，同時滿足公共點在第1站的坐標在經過多次轉站后，當其轉回第1站時，坐標不變，即

若式（7）成立，則有

式（8）和式（9）即為約束條件方程。

2.2 約束條件方程的求解方法

在平差計算時，可將約束條件方程寫成如下形式：

對式（10）作泰勒展開，可得：

令

則式（12）可寫成如下形式：

文獻[17]對式（13）的線性化處理方法進行了詳細闡述，并得到了對應參數的具體表達形式。

3 平差值的求解方法

根據約束條件方程構建拉格朗日函數：

的參數矩陣。

求偏導，設此時 Φ的偏導數為零，則 Φ有極小值，即此時誤差最小，則

經整理，可得：

將式（1）代入，可得：

令

則式（18）可以寫成

令

則式（20）可寫成

則有

將式（23）代入式（18），可得：

同時，未加約束條件的平差值經整理，可以表示為：

4 實驗驗證

本次實驗在船臺測量場的基礎上，對確定公共點在基準坐標系中位置的轉站測量過程進行了模擬，對利用聯合平差算法對所測得的公共點坐標數據進行處理，并通過與未進行聯合平差的結果進行對比，對基于聯合平差的船臺測量場標定方法進行驗證。

4.1 實驗過程

1）在實驗場地內的不同位置共布設22個公共點，在第1測量站位架設全站儀，為了對船臺測量場的實際情況進行模擬，各站位所測得的公共點數量均不相同，每個公共點測量10次。

2）將全站儀轉移至下一站位，重復過程1，直至測完6站，并對測量數據進行記錄。

4.2 實驗結果

在本次實驗中，采用每一測站下每個公共點10次測量數據的平均值作為當前站位的測量坐標，同時，將相鄰的兩站兩兩組合，利用坐標轉換算法對角度參數R和距離參數T進行計算，共得到6組數據。

利用聯合平差算法對測量數據進行處理，其中，平差初值為通過計算得到的6組坐標轉換參數，約束條件方程根據公共點在第1站的坐標進行一輪坐標轉換后返回第1站時與原坐標的差值最小列出，并設置迭代計算時角度閾值為1″，距離閾值為0.1mm。最后，將經過聯合平差處理得到的中誤差和坐標轉換參數與未經過聯合平差處理的數據進行對比。

表1為在經過聯合平差計算前后，所用站位所測得的公共點經過坐標轉換統一至第1站的坐標數據。表2和表3分別表示不同轉站過程的測量數據，在經過聯合平差計算前后坐標轉換參數誤差的對比。從表中數據可以看出，公共點的測量數據經過聯合平差計算后，其坐標轉換參數誤差整體小于聯合平差計算前的坐標轉換誤差。還可以看出，在進行聯合平差前后，其坐標轉換參數誤差的差值較小，這是因為在本次實驗中，每一站為所測的公共點數量較多，并且每個公共點的坐標值均是10次重復測量取平均得到，公

表1 聯合平差前后的公共點坐標Tab. 1 Coordinates of common points before and after joint adjustment

表2 未經聯合平差的轉站參數誤差Tab. 2 Transformation parameter errors without joint adjustment

表3 聯合平差后的轉站參數誤差Tab. 3 Transformation parameter errors after joint adjustment

共點測量坐標本身所帶有的測量誤差較小，因此，坐標轉換參數誤差較小，其差值也較小。表4為在經過聯合平差前后，第6站的公共點坐標在經過一輪坐標轉換后，返回至第6站時的坐標與原坐標的差值，由表中數據可以看出，在經過聯合平差之后，第6站公共的坐標經過坐標轉換返回至第6站后與原坐標的差值更小，進一步說明了利用聯合平差算法對公共點坐標進行處理可以有效提升船臺測量的整體精度。

表4 聯合平差前后第6站的坐標差Tab. 4 Coordinate difference of 6th station before and after joint adjustment

5 結 語

在利用多級轉站的測量方法構建船測量場的過程中，由于公共點測量誤差和坐標轉換算法誤差的存在，會導致不同站位所測得的公共點坐標統一至基準坐標系后與其真實空間位置產生一定的偏差，從而影響船臺測量場的精度。本文采用聯合平差算法實現了多級轉站時最優坐標轉換參數的求解，該方法基于最小二乘原理和誤差最小原則，給出了船臺測量場公共點經過多級轉站測量的誤差方程及其求解方法，并根據公共點經過坐標轉換返回原站位坐標不變的原則，給出了對應的約束條件方程及其求解方法。最后，根據拉格朗日函數，得到了平差值的求解方法。實驗證明，本文所提出的船臺測量場聯合平差標定算法有效減少了船臺測量場的誤差，提高了利用船臺測量場進行測量作業的精度。