孟慶旭，朱曉軍，彭 飛，王 中

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

船臺測量場中公共點的粗差判別方法

孟慶旭，朱曉軍，彭 飛，王 中

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

針對船臺公共點含粗差時對測量場整體精度的影響，提出了一種穩(wěn)健的公共點粗差判別方法。該方法先通過萊以特準則初步篩選得到可疑點，再以羅曼諾夫斯基準則對可疑點進行判別，確定異常點并剔除，最后采用SVD算法計算轉(zhuǎn)站參數(shù)。通過仿真試驗進行驗證，結(jié)果表明，該方法能準確地剔除異常點，得到可靠的轉(zhuǎn)站參數(shù)估值，從而提高測量場整體精度。

ERS點；粗差判別；測量場：轉(zhuǎn)站參數(shù)；SVD算法

船臺高精度測量場用于潛艇耐壓殼體分段的對接合攏與變形測量，以保證潛艇建造質(zhì)量。測量場分測量坐標系和裝配坐標系，裝配坐標系用于構(gòu)建測量場的基準坐標系，測量坐標系則是不同的儀器設(shè)備在不同的站位所構(gòu)建的局部坐標系。測量坐標系之間相互獨立，要通過一組公共點來實現(xiàn)測量坐標系到裝配坐標系的轉(zhuǎn)換，這些公共點被稱為增強系統(tǒng)參考點(enhanced reference system point，ERS point)，通常被布置在裝配現(xiàn)場，是裝配坐標系的參考基準。在進行大空間測量過程中，ERS點的粗差會影響坐標轉(zhuǎn)換參數(shù)的精度，在轉(zhuǎn)站過程中產(chǎn)生誤差，轉(zhuǎn)站誤差是決定高精度測量場重要的因素之一。

在室內(nèi)造船現(xiàn)場等惡劣的環(huán)境下，腳手架的遮擋、基準平臺的振動、人為操作的失誤都會使測量坐標系中ERS點產(chǎn)生較大的誤差，如果將這些誤差直接帶入坐標系融合過程中，將帶來不可預(yù)知的誤差。此外ERS點集一般都存在冗余，為控制融合誤差，有必要在數(shù)據(jù)融合之前，先對測量坐標系中ERS點的有效性進行檢驗，剔除粗差點以控制測量場總體精度。劉東明提出了一種適用于大旋轉(zhuǎn)角的空間三維直角坐標基準的轉(zhuǎn)換方法，并發(fā)現(xiàn)ERS點粗差對轉(zhuǎn)站參數(shù)影響很大，但沒有給出相應(yīng)的解決辦法[1]。高精度測量場在飛機裝配中已成熟運用，它的構(gòu)建基礎(chǔ)是室內(nèi)定位技術(shù)(iGPS)，不存在ERS點的誤差問題[2]。Shamir基于目標點匹配誤差最小原則，提出了對ERS點的選擇優(yōu)化方法，該方法可以通過合理地選擇或剔除相應(yīng)的ERS點來提高測量精度[3]。吳祖海提出了一種坐標轉(zhuǎn)換公共點粗差定位和降低粗差點影響的新方法，即利用最小二乘平差后的方差，采用方差比值檢驗法，實現(xiàn)在布爾莎七參數(shù)計算過程中的粗差定位和剔除[4]。郭英起將粗差歸入隨機模型，選擇適當?shù)臋?quán)函數(shù)消除或削弱對轉(zhuǎn)站參數(shù)的影響，但是合理權(quán)函數(shù)的選取并不容易[5]。以上幾種方法不適用于大空間高精度測量場中ERS點的粗差判別，難以保證測量場的高精度要求。

本文以羅曼諾夫斯基準則和萊以特準則為基礎(chǔ)[6]，利用空間相似性對ERS點集進行判別，篩選并剔除含有粗差的點，將優(yōu)化后的點集用奇異值分解法(SVD)[7]計算可得到較高精度的轉(zhuǎn)站參數(shù)。

1 激光全站儀轉(zhuǎn)站原理

在艦船建造現(xiàn)場，由于船體分段較大和現(xiàn)場設(shè)備的遮擋，需要采用多臺激光全站儀或單臺激光全站儀多站位相結(jié)合的方法完成對船體分段的測量[8-9]。每臺全站儀在各站位形成獨立的測量坐標系，整個室內(nèi)船臺構(gòu)成裝配坐標系，通過ERS點將測量坐標系和裝配坐標系統(tǒng)一。激光全站儀轉(zhuǎn)站的實質(zhì)就是ERS點集的三維剛體匹配問題。為實現(xiàn)轉(zhuǎn)站和保證測量場得的精度，ERS點需合理的分布在整個船臺測量場，對裝配坐標系下的ERS點重復(fù)標定，標定結(jié)果可作為精確理論值。假設(shè)P={p1,p2,…pN}是裝配坐標系中的ERS點集，Q={q1,q2,…qN}是測量坐標系中的ERS點集，N為點集的數(shù)量，轉(zhuǎn)站參數(shù)用T和R來表示，分別為全站儀測量坐標系相對裝配坐標系的平移參數(shù)和旋轉(zhuǎn)參數(shù)。三維剛體相似變換公式為

P=RQ+T

式中

式中，Δx、Δy、Δz分別表示坐標系沿各坐標軸的平移參數(shù)參量；f、w、k分別表示坐標軸的旋轉(zhuǎn)角參數(shù)。

全站儀多站位轉(zhuǎn)站如圖1所示，O1-X1Y1Z1和O2-X2Y2Z2為測量坐標系，方形點表示ERS點，所有的ERS點組成一個空間體，兩個坐標系通過測量ERS點的位置進行轉(zhuǎn)站，虛線表示激光測量光束。由于激光全站儀每次轉(zhuǎn)站測量前需要進行水平調(diào)整，不考慮誤差因素，假設(shè)兩個坐標系轉(zhuǎn)站過程中Z軸的方向始終保持豎直不變，則X軸和Y軸的旋轉(zhuǎn)角f和w為零。若要求解到最優(yōu)的轉(zhuǎn)站參數(shù)T和R，構(gòu)建最小二乘誤差函數(shù)∑2并取最小值

(1)

針對最小二乘誤差函數(shù)可通過奇異值分解法(SVD)求解得到轉(zhuǎn)站參數(shù)R和T。

2 ERS點的粗差判別模型

2.1 ERS點偏差估計

檢驗測量坐標系中的ERS點集是否滿足要求，實際上就是檢驗其是否與裝配坐標系中的ERS點集匹配，點匹配問題起源于多維尺度分析中兩相似點組之間的空間位置關(guān)系的求解問題，因此可以采用相似度評估的方式來對其進行檢驗。相似度的評估以點集相互距離的匹配性和點集向量空間方向的相似度(空間布局的相似度)為條件對ERS點進行誤差估計。

圖2 ERS點集分布示意圖

(2)

式中，pi和qi分別是點集P和Q的ERS點，N為ERS點的數(shù)量。點集相互距離的匹配性可以用兩個點集中對應(yīng)點到質(zhì)心的距離差Δdi來表示。

(3)

點集向量空間方向的相似度Δci則可以采用向量空間余弦相似度的方式來衡量，即

Δci=|sim(pi,pz)-sim(qi,qz)|

(4)

式中，pz是裝配坐標系下與Z軸平行的單位向量；qz是測量坐標系與Z軸平行的單位向量；pi和qi分別表示在各自坐標系下與坐標原點形成的三維空間向量。

(5)

通過距離相似度和向量空間方向余弦相似度的計算，可以得到兩個相似度集合。距離相似度集合ΔD={Δd1，Δd2，…，ΔdN}和向量空間方向相似度集合ΔC={Δc1，Δc2，…，ΔcN}。

2.2 粗差點的剔除

對于相似度集合ΔD={Δd1，Δd2，…，ΔdN}和ΔC={Δc1，Δc2，…，ΔcN}，需要對其進行粗大誤差的判斷，以消除人為因素的影響。粗差誤差點的判別可以分兩步進行，首先計算殘差，并采用3σ準則對測量坐標系中的ERS點集進行判斷，找出異常點。vdi和vci分別是兩點相互距離和空間向量余弦的殘差方程。

(6)

σd和σc分別是兩點相互距離和空間向量余弦的標準差方程。

(7)

如果測量坐標系中的ERS點i的殘差|vdi|>3σd或|vci|>3σc，則認為該點為粗差可疑點。由于ERS點的樣本空間一般不大，3σ準則對于小樣本空間存在一定的局限性，若不存在點i滿足殘差|vdi|>3σd或|vci|>3σc，則選取|vdi|或|vci|的最大值作為異常點。

然后在3σ準則的基礎(chǔ)上，按t分布的實際誤差分布范圍來判斷粗大誤差，即采用羅曼諾夫斯基準則，也稱為t檢驗準則來進行判別。假設(shè)j點被認為是粗差點，則將其從數(shù)列中剔除，剔除后再計算其平均值。

(8)

數(shù)列的標準差則變化為

(9)

根據(jù)ERS點的數(shù)量N和選取的顯著度α，以及查t分布表得到檢驗系數(shù)K(N,α)[6]。如果|vdi|>Kσd或|vci|>Kσc，則認為i點為粗大誤差，并且將其從ERS點集中剔除，否則予以保留。在多個粗差同時存在的情況下，為防止各個粗差相互干擾，特別是大粗差掩蓋小粗差，每輪只剔除一個粗差點，剔除后再開始第二輪檢驗，直至所有的粗差被剔除。

3 數(shù)值仿真試驗

為驗證ERS點誤差模型的準確性和可靠性，在Matlab軟件平臺上編寫模型計算程序和SVD算法求解程序并進行仿真試驗。用激光全站儀測量一組數(shù)據(jù)作為裝配坐標系下點集P的數(shù)據(jù)，ERS點數(shù)量N=9，假設(shè)轉(zhuǎn)站參數(shù)Δx=3000，Δy=4000，Δz=5000，f=0，w=0，k=-π/3。通過Matlab計算出點集Q的理論坐標值，假設(shè)9號點為含有粗差的異常點，1—8號點加入服從正態(tài)分布Δξi～N(0,5)的高斯噪聲，9號點加入服從正態(tài)分布Δεi～N(5,10)的高斯噪聲，得到含有誤差測量坐標系ERS坐標值(見表1)。

將裝配坐標系和測量坐標系下的ERS坐標進行計算得到距離殘差和角度殘差,如圖3、圖4所示，標準差為σd=4.145 647和σc=0.001 242，通過3σ準則發(fā)現(xiàn)無滿足|vd|>3σd或|vc|>3σc的異常點，因此選取|vd|或|vc|最大的點為異常點，即3號點或9號點。

表1 ERS點坐標值 mm

圖3 距離殘差

圖4 角度殘差

為了驗證ERS點的優(yōu)化模型，將剔除前后的ERS點坐標值用SVD算法分別求解轉(zhuǎn)站參數(shù)，計算結(jié)果如表2所示。結(jié)果表明，將含有粗差的點剔除后計算得到的轉(zhuǎn)站參數(shù)與設(shè)定的真值相比偏差較小，具有較高精度。將含有粗差的ERS點參與到轉(zhuǎn)站參數(shù)的計算中會使平移參數(shù)產(chǎn)生較大偏差，并使坐標軸X和Y發(fā)生微小旋轉(zhuǎn)，轉(zhuǎn)站參數(shù)的誤差會產(chǎn)生轉(zhuǎn)站誤差，并且影響到整個測量場的精度。

表2 計算結(jié)果對比

4 結(jié) 語

在船臺大尺寸空間測量場內(nèi)進行轉(zhuǎn)站測量，由于環(huán)境因素和人為因素會導致部分ERS點產(chǎn)生粗差，帶有粗差的ERS點若參與到轉(zhuǎn)站參數(shù)的計算中，會影響轉(zhuǎn)站精度。本文通過萊以特準則和羅曼諾夫斯基準則對測量坐標系下的ERS點進行優(yōu)化，將含有粗差的ERS點剔除，不參與轉(zhuǎn)站參數(shù)的計算，保證了在復(fù)雜環(huán)境下轉(zhuǎn)站測量的精度，進一步提高了測量場的精度和適用性。該方法還有進一步優(yōu)化空間，考慮用不確定度估計的方法來判別粗差點，與測量場的不確定度估計形成統(tǒng)一，有利于船臺測量場的實際應(yīng)用。

[1] 劉東明，郭際明，施展，等．大角度三維基準轉(zhuǎn)換及公共點的確定[J].測繪通報，2010(1)：11-14.

[2] 喻世臣，康曉峰，翟南，等．基于激光跟蹤儀、iGPS的飛機水平測量技術(shù)研究[J].航空制造技術(shù)，2015，491(21)：119-121.

[3] SHAMIR R R, JOSKOWICZ L, SHOSHAN Y. Optimal Landmarks Selection and Fiducial Marker Placement for Minimal Target Registration Error in Image-guided Neurosurgery[C]∥Proceedings of the SPIE. [S.l.]: SPIE, 2009.

[4] 吳祖海，羅偉釗，李軍．坐標轉(zhuǎn)換中公共點粗差定位與降低粗差點影響方法研究[J].大地測量與地球動力學，2014，34(1)：118-122.

[5] 郭英起，康彬，張秋江，等.基于空間直角坐標系的高精度坐標轉(zhuǎn)換方法研究[J].大地測量與地球動力學，2012，32(3)：125-128.

[6] 費業(yè)泰．誤差理論與數(shù)據(jù)處理[M]．北京：機械工業(yè)出版社，2015.

[7] EGGERT D W, LORUSSO A, FISHER R B. Estimating 3-D Rigid Body Transformations: A Comparison of Four Major Algorithms[J].Machine Vision and Applications,1997,9(5-6):272-290.

[8] 朱志潔，趙耀，王中，等．船體外板點位信息轉(zhuǎn)站測量方法研究[J]．計算機工程，2015，41(10)：286-289.

[9] 潘國榮，張鵬，孔寧．造船精度控制系統(tǒng)中用移站測量獲取點位信息的一種方法[J]．大地測量與地球動力學，2010，30(5)：121-124．

[10] 金漲軍．飛機裝配中大尺寸測量場的建立與優(yōu)化技術(shù)[D].杭州：浙江大學，2016.

[11] 屈昊．基于對偶四元數(shù)的三維物體匹配方法的研究[D].哈爾濱：哈爾濱工程大學，2004.

[12] 楊福芹，戴華陽，鄒定輝，等．一種公共點坐標系轉(zhuǎn)換中的粗差探測方法及其應(yīng)用[J]. 測繪通報,2015(8):38-41.

[13] 郭英起，曹先革，劉江．公共點坐標粗差檢驗方法及其限差的研究[J]. 全球定位系統(tǒng),2012(4):30-32.

[14] 陳義，陸環(huán)．以三維坐標轉(zhuǎn)換為例解算穩(wěn)健總體最小二乘方法[J].測繪學報，2012, 41(5):715-722.

[15] 王建民，張錦，蘇巧梅．觀測數(shù)據(jù)中粗差定位與定值算法[J].武漢大學學報(信息科學版)，2013, 38(10):1225-1228.

TheGrossErrorDiscriminationMethodofERSPointinSlipwayMeasurementField

MENG Qingxu，ZHU Xiaojun,PENG Fei，WANG Zhong

(Department of Naval Architecture and Ocean Engineering，Naval University of Engineering，Wuhan 430033，China)

In view of the influence of the gross error of the ERS point on the accuracy of the measurement field, this paper presents a robust method for judging the gross error of ERS point. In this method, the suspicious points are obtained by the initial screening of the Lai criteria, and then the faulty points are discriminated by the Romanovsky criterion, the anomaly points are determined and removed, and finally the SVD algorithm is used to calculate the station parameters. The results show that the method can accurately eliminate and improve the overall accuracy of the measurement field.

ERS point; gross error discrimination; measurement field; station parameters; SVD algorithm

2017-03-31

國家自然科學基金(51609253)

孟慶旭(1993—)，男，碩士生，研究方向為艦船精度建造技術(shù)研究。E-mail:978536456@qq.com

彭 飛

孟慶旭，朱曉軍，彭飛,等.船臺測量場中公共點的粗差判別方法[J].測繪通報，2017(12)：103-106.

10.13474/j.cnki.11-2246.2017.0389.

P25

A

0494-0911(2017)12-0103-04