非模型化相機標(biāo)定方法及測角誤差研究*

李華晉，鄒 劍，王 偉，林嘉睿，邾繼貴，葉聲華

(天津大學(xué)精密測試技術(shù)及儀器國家重點實驗室，天津 300072)

在產(chǎn)品存在復(fù)雜曲面的工程與工業(yè)生產(chǎn)行業(yè)，如汽車、飛機、船舶等，基于圖像的視覺測量技術(shù)以其非接觸式、操作簡單、適應(yīng)性強、工作環(huán)境要求低等特點，越來越多地被應(yīng)用于裝配和質(zhì)量控制中。［1－3］相機標(biāo)定是視覺測量的關(guān)鍵技術(shù)和步驟，目前已產(chǎn)生了很多種標(biāo)定方法，如傳統(tǒng)標(biāo)定方法、自標(biāo)定方法、主動視覺標(biāo)定方法等。［4］這些方法均建立了畸變模型或補償畸變的成像模型，解算相機主距(f)、主點(Cx，Cy)等模型參數(shù)，例如傳統(tǒng)方法的經(jīng)典代 表:DLT 方 法［5、6］、Tsai 方 法［7、8］、Weng 迭 代法［9］等。此類方法中，由于模型各參數(shù)間的相關(guān)性、非線性求解的不確定性、非模型化參數(shù)等因素的影響使模型補償存在局限性，難以進(jìn)一步提高精度以滿足工業(yè)現(xiàn)場越來越高的要求。

非模型化相機標(biāo)定拋棄了模型求解方式，對相機成像進(jìn)行光學(xué)模擬標(biāo)定，逐像元精密標(biāo)定CCD全像面，建立視場角度－像面坐標(biāo)映射表，實現(xiàn)高精度的畸變細(xì)化修正。此方法具有可補償非模型化參數(shù)，無模型求解的不確定性問題等優(yōu)點;但一維靶標(biāo)的加工精度、圖像處理精度、分度機構(gòu)的精度等關(guān)鍵誤差項會影響標(biāo)定精度，故研究其誤差分布特性，對此標(biāo)定技術(shù)的發(fā)展具有重要意義。

1 標(biāo)定原理

非模型化相機標(biāo)定是獲得分度映射表的過程，每對一個方向標(biāo)定可以得到一個分度映射表，一般選擇水平、豎直兩個正交方向?qū)ο鄼C進(jìn)行兩次標(biāo)定，由于標(biāo)定過程相同，本文僅研究水平方向。

根據(jù)垂線法，空間直線在一個視場方位成像為像面上的一條曲線。在某一方向上，直線以固定視場角間隔連續(xù)旋轉(zhuǎn)，則在像面上的對應(yīng)為一系列曲線，并將像面細(xì)分成連續(xù)的角度分度區(qū)域。使用多齒分度臺帶動相機轉(zhuǎn)動，以相對運動的方式實現(xiàn)一維靶標(biāo)在不同視場角下拍攝。如圖1所示，多齒分度臺帶動相機，以固定分度角間隔精密旋轉(zhuǎn)，每個視場拍攝一次靶標(biāo)，處理圖像提取靶標(biāo)像點，即分度點坐標(biāo)，與視場角一起存入分度映射表，經(jīng)過對全像面的拍攝提取，從而實現(xiàn)對整個像面的標(biāo)定。

圖1 相機標(biāo)定裝置原理圖

像面的細(xì)分過程如圖2所示，分度點隨視場角沿u方向增加，逐漸覆蓋整個像面。每列分度點對應(yīng)一個標(biāo)定角度，像面被劃分成多個分度塊(分度塊:在u、v方向上均相鄰的4個分度點圍成的像面區(qū)域)，分度角度間隔和一維靶標(biāo)點間隔分別決定了分度塊u、v方向邊長。在理想針孔成像模型下，u坐標(biāo)值與視場角度值一一對應(yīng)，像點u坐標(biāo)直接代表視場角度。實際上，成像畸變及其他因素的影響，使在同一視場方位成像的一列分度點以曲線，甚至鋸齒線形排列，因此需要同時使用u、v坐標(biāo)，按與像點臨近的多個分度點插值解算角度值。

圖2 像面細(xì)分過程

被測像點處于某個分度塊之中，當(dāng)分度角度間隔足夠小時，便可使用線性插值算法，其僅會引入極微小的分度插值誤差(0.000 1 pixel量級)。一般情況下，分度塊是不規(guī)則的四邊形，如圖3所示，被測A點在1，2，3，4為頂點的分度塊內(nèi)，按照線性加權(quán)插值法，以u坐標(biāo)為主要權(quán)重、v坐標(biāo)為輔助權(quán)重插值，則A點的空間角度θA(uA，vA)見式(1):

式中:UA1為被測點A與分度點1間的u坐標(biāo)間距，其他類同;VA1為被測點A與分度點1間的v坐標(biāo)間距，其他類同;θ1為分度點1處的角度值;θE為分度角度間隔。

圖3 角度插值算法

2 測角誤差分析

分析標(biāo)定過程，影響標(biāo)定精度的關(guān)鍵因素有:一維靶標(biāo)自身精度、分度機構(gòu)性能、圖像處理精度。他們共同作用于分度點坐標(biāo)、分度間隔、被測點坐標(biāo)，影響測角精度。各參數(shù)項如表1，其中分度點和被測像點坐標(biāo)均是使用質(zhì)心法［10－13］進(jìn)行圖像處理獲得的，但前者作用于標(biāo)定過程，而后者作用于測量過程。前者經(jīng)過一次提取，存入分度映射表，在測量時每次讀取固定不變，屬于系統(tǒng)誤差;后者直接影響測量結(jié)果，每次測量引入的誤差值不同，屬隨機誤差。

表1 影響測角結(jié)果的誤差項

2.1 系統(tǒng)誤差

設(shè)一維靶標(biāo)上有N個靶點，空間標(biāo)定角度值與像面上的分度列/線一一對應(yīng)，分度列包含N個分度點，故標(biāo)定角度值與分度點存在一對多的關(guān)系。一維靶標(biāo)的理想情況是:N→∞，且靶標(biāo)點沿嚴(yán)格的直線排列，那么此時的分度線一般是連續(xù)且光滑的曲線(因受成像畸變影響)，定義為理想分度線。實際系統(tǒng)中，由于靶點直線度、圖像處理精度等影響，分度線與理想分度線間存在偏差，影響著分度線左右以N個分度點為頂點的2N個分度塊，在分度塊內(nèi)插值提取角度時引入測角誤差。

如圖4，分度角 θ1、θ2、θ3對應(yīng)的分度點列的一種局部分布情況，O、P、Q、R是 θ2的鋸齒形分度線的分度點拐點，虛線段l是理想分度線的一部分，陰影部分是實際分度線與理想分度線的偏差，PQ段兩者重合。當(dāng)被測點沿理想分度線軌跡l移動，則引入的誤差項沿V方向分布曲線與分度線形狀吻合，若被測點的v坐標(biāo)保持不變，則引入的θd為定值。θd見式(2)。

圖4 平滑度誤差項的影響過程示意圖

式中:x為分度線與理想分度線的偏差;f為相機的有效焦距。

為分析平滑度對測角誤差的綜合影響特性，以誤差項極值θD作為引入誤差項的評價標(biāo)準(zhǔn)，將平滑度a代入式(2)，則θD:

系統(tǒng)的測角范圍為±10°，若分度角間隔固定，則像面分度間隔f(tanθ3－tanθ2)變化不大，設(shè)為 2 pixel，θE為40″時，a－θD關(guān)系曲線如圖 5。可見降低分度線平滑度可以線性降低引入的誤差項極值，當(dāng)直線度保持在0.2 pixel時，系統(tǒng)誤差可降低到4″以下。由表1知，分度線平滑度的降低，可以通過提高靶點空間直線度、采用更精密的坐標(biāo)提取方法實現(xiàn)。

圖5 a－θD關(guān)系曲線

2.2 隨機誤差

測角的隨機誤差項有像點坐標(biāo)提取誤差和分度角間隔。像點坐標(biāo)的提取是實現(xiàn)測量的基礎(chǔ)及關(guān)鍵步驟，系統(tǒng)選用紅外LED作為特征點，使用質(zhì)心法提取像點中心坐標(biāo)，提取誤差服從正態(tài)分布，σuA=0.05 pixel、σvA=0.05 pixel，σθE=0.6″。

式(1)中參與計算的5個像點坐標(biāo)中，只有被測像點坐標(biāo)(uA，vA)引入的是隨機誤差項。按式(4)合成隨機誤差標(biāo)準(zhǔn)差，當(dāng) θE=40″時 σθ=1.60″［14－15］。

3 實驗

為驗證標(biāo)定方法的可行性和關(guān)鍵誤差研究的正確性，設(shè)計實驗:按照標(biāo)定裝置布置實現(xiàn)，但僅打開一個紅外LED代替整個一維靶標(biāo)作為被測點，如圖6，驗證測角結(jié)果。實驗設(shè)備有:FL2G－50S5M－C工業(yè)相機配合Schneider 23 mm鏡頭;一維靶標(biāo)是29個紅外LED靶點構(gòu)成，標(biāo)定距離3 m，以40″為分度間隔。在水平角－10°～10°內(nèi)標(biāo)定，u方向分布1 876個分度列，每一列有29個分度點。見圖7，分度列間的分度角度間隔相同，但像素間距不同;分度列內(nèi)各分度點像素間距不同。

圖6 實驗現(xiàn)場圖

在－10°～10°范圍內(nèi)，以 0.5°間隔測角 3 個測回，獲得240個共40個角度處的測量數(shù)據(jù)，使用往、返測回的平均值作為測角數(shù)據(jù)，每個測回的測角結(jié)果偏差體現(xiàn)了測角隨機誤差。

圖7 像面分度點分布

但實驗存在調(diào)整誤差，如圖8所示由于光心與轉(zhuǎn)軸存在重合度偏差，即標(biāo)定坐標(biāo)系與驗證坐標(biāo)系不重合，偏心距R導(dǎo)致驗證誤差(θ'－θ0)隨θ0呈三角函數(shù)分布，可通過將誤差數(shù)據(jù)進(jìn)行三角函數(shù)擬合去除。

圖8 校準(zhǔn)偏心模型

剔除校準(zhǔn)誤差后的水平角測角誤差分布如圖9。整體誤差分布趨勢如圖10，出現(xiàn)的拐點對應(yīng)圖(4)中分度線的拐點，說明驗證機構(gòu)在與水平面成一定角度的平面內(nèi)旋轉(zhuǎn)，被測像點穿越了分度線的不同區(qū)間段，引入了平滑度系統(tǒng)誤差項。將實驗參數(shù)代入式(3)有:

圖9 測角誤差分布

圖10是測角誤差項中的θd分布，在±6″以內(nèi)。測回間的數(shù)據(jù)偏差體現(xiàn)了隨機誤差項，如圖11所示，隨機誤差分布在±4″以內(nèi)，測角誤差標(biāo)準(zhǔn)差為1.68″與2.2 節(jié)的理論值 σθ=1.60″基本吻合，測角不確定度為5.04″(99.73%的置信概率)。

圖10 θd的分布

圖11 隨機誤差分布

4 結(jié)論

本文首先介紹了非模型化相機標(biāo)定方法的原理、標(biāo)定過程、測角模型。其次，分析了關(guān)鍵系統(tǒng)誤差－分度線平滑度的作用過程和分布特性;對隨機誤差項進(jìn)行合成。最后，設(shè)計實驗驗證了測角結(jié)果及誤差分析的正確性，得到了系統(tǒng)目前條件下，系統(tǒng)誤差－平滑度誤差項的分布特性(20″/pixel)及測角不確定度5.04″(99.73%的置信概率)。標(biāo)定方法的誤差分析，對標(biāo)定參數(shù)的選取有重要的指導(dǎo)意義，系統(tǒng)誤差在目前配置下是在±6″內(nèi)分段線性分布，由于直接受一維靶標(biāo)的靶點直線度影響，故可通過提高一維靶標(biāo)的制造精度進(jìn)一步提高標(biāo)定精度，有很好的應(yīng)用前景。

［1］張德海，梁晉，唐正宗，等.大型復(fù)雜曲面產(chǎn)品近景工業(yè)攝影測量系統(tǒng)開發(fā)［J］.光電工程，2009，36(5):122－128.

［2］李小春，杜學(xué)飛，秦子建，等.模型機航空近景攝影測量系統(tǒng)的開發(fā)和應(yīng)用［J］.巖土力學(xué)，2009，30(2):318－322.

［3］杜歆，徐鋼梅，朱云芳，等.一種基于單相機三維重建的簡單方法［J］.傳感技術(shù)學(xué)報，2007，20(8):1917－1920.

［4］張丹.基于圖像的模擬相機標(biāo)定方法的研究［D］.長春理工大學(xué)，2008.

［5］Abdel-Aziz Y I，Karara H M.Direct Linear Transformation from Comparator Coordinates into Object Space Coordinates［C］//ASP Symposium on Close-Range Photogrammetry，1971，1:1－18.

［6］Faugeras O D，Toscani G.The Calibration Problem for Stereo［C］//Proceedings of IEEE Conference on Computer Vision and Pattern Recognition，1986，1:15－20.

［7］Tsai R，Lenz R K.A Technique for Full Autonomous and Efficient 3D Robotics Hand/Eye Calibration［J］.IEEE Trans Robotics and Automation，1989，5(3):345－358.

［8］Tsai R.An Efficient and Accurate Camera Calibration Technique for 3D Machine Vision［C］//Proc IEEE Conference on Computer Vision and Pattern Recognition，1986，1:364－374.

［9］Weng J Y，Cohen P，Hernion M.Camera Calibration with Distortion Model and Accuracy Evaluation［J］.IEEE Trans Pattern Analysis and Machine Itelligence，1992，14(10):965980.

［10］魏振忠，高明，張廣軍，等.一種光斑圖像中心的亞像素提取方法［J］.光電工程，2009，36(4):7－12.

［11］陳運錦，馮瑩，魏立安，等，光斑質(zhì)心亞像素定位誤差的實驗研究［J］.光電工程，2010，37(2):80－84.

［12］楊君，張濤，宋靖雁，等.星點質(zhì)心亞像元定位的高精度誤差補償法［J］.光學(xué)精密工程，2010，18(4):1002－1010.

［13］Yang J，Zhang T，Shong J Y，et al.High Accuracy Error Compensation Algorithm Forstar Image Sub-Pixel Subdivision Location［J］.Opt.Precision Eng.，2010，18(4):1002－1010

［14］朱曉輝，郭文佳，王向軍.超大范圍坐標(biāo)測量系統(tǒng)的現(xiàn)場測試基準(zhǔn)及其誤差分析［J］.傳感技術(shù)學(xué)報，2005，18(2):265－268.

［15］費業(yè)泰.誤差理論與數(shù)據(jù)處理［M］.北京:機械工業(yè)出版社，2004，79－93.