吳光宇,崔元達(dá),陳璟逸
(蘇州北美國際高級中學(xué),江蘇 蘇州 215104)
隨著經(jīng)濟(jì)全球化和航空、航天、汽車領(lǐng)域的迅猛發(fā)展,大型精密復(fù)雜產(chǎn)品的裝配工作量達(dá)到了50%以上,如何滿足大型復(fù)雜產(chǎn)品的快速高精度裝配成了衡量國家制造業(yè)水平的關(guān)鍵因素之一。由于大型復(fù)雜產(chǎn)品存在測量體積大,測量精度要求高、測量效率要求高等顯著的特點,傳統(tǒng)的測量方法已經(jīng)無法滿足實際的工業(yè)生產(chǎn)需求。
近幾年,非接觸式測量方法在各行各業(yè)得到了廣泛的應(yīng)用,例如雙目視覺、多目視覺、線激光、面激光、光柵投影以及攝影測量等,都是利用工業(yè)相機(jī)對圖像進(jìn)行采集處理得到物體的三維模型。
線激光測量方法因其具有非接觸不損傷物體表面,測量效率高、測量精度高、抗環(huán)境光干擾能力強(qiáng)以及信息量大等優(yōu)點,在工業(yè)三維重建[1-5]、缺陷檢測[6-8]以及目標(biāo)跟蹤[9-12]等領(lǐng)域得到了廣泛的應(yīng)用。
國外對光學(xué)掃描的研究起步較早,三維測量設(shè)備也陸續(xù)投入到了實際的工業(yè)生產(chǎn)需要中,德國GOM公司的ATOS工業(yè)三維光學(xué)掃描技術(shù)能夠快速地提供精準(zhǔn)和細(xì)膩的掃描數(shù)據(jù)[13]。然而在測量過程中,掃描設(shè)備移動時,需要將掃描設(shè)備采集的點云數(shù)據(jù)轉(zhuǎn)換到世界坐標(biāo)系下完成整個物體的拼接過程,因此需要在物體表面粘貼大量的標(biāo)識點便于坐標(biāo)系的統(tǒng)一,給自動化測量帶來了較大的困難?;萜展旧a(chǎn)的ProS3白色結(jié)構(gòu)光三維掃描設(shè)備[14]借助于雙目視覺系統(tǒng)結(jié)合白色結(jié)構(gòu)光組成視覺處理單元對物體進(jìn)行三維掃描,可以捕獲到掃描對象的表面紋理,借助于旋轉(zhuǎn)平臺進(jìn)行全自動拼接,實現(xiàn)360°模型測量。但是該設(shè)備只適合于小型部件的測量,而且需要結(jié)合旋轉(zhuǎn)平臺進(jìn)行全景測量,不適合大型產(chǎn)品部件的掃描。
近幾年,國內(nèi)三維掃描技術(shù)突飛猛進(jìn),先后研發(fā)出多種工業(yè)三維掃描儀,天遠(yuǎn)公司的OKIO產(chǎn)品采用藍(lán)色激光對工業(yè)零部件進(jìn)行三維掃描,抗干擾能力強(qiáng)。但是該設(shè)備仍然需要粘貼標(biāo)志點實現(xiàn)全景測量。天遠(yuǎn)公司之后推出的FreeScanTrak無線跟蹤式激光掃描系統(tǒng)基于動態(tài)光學(xué)跟蹤原理,對掃描頭進(jìn)行跟蹤定位并實時精確測量目標(biāo)的三維形狀,實現(xiàn)了在無需貼點條件下的高精度三維空間點云掃描。但是該設(shè)備需要利用固定的跟蹤儀實時跟蹤掃描頭的位置,因此掃描范圍有限,無法實現(xiàn)大型零部件物體的三維測量[15]。
在實際的三維測量工程應(yīng)用中,測量的數(shù)據(jù)不可避免存在誤差,導(dǎo)致測量結(jié)果呈現(xiàn)鋸齒和柔順性較差的情況,為之后產(chǎn)品加工、制造帶來較大的障礙。因此,在航空航天領(lǐng)域、船舶以及汽車等領(lǐng)域,在數(shù)據(jù)采集完成后都需要對數(shù)據(jù)進(jìn)行光順處理。目前有大量的研究學(xué)者試圖實現(xiàn)曲線的自動光順,然而在工程應(yīng)用中并沒有得到有效應(yīng)用[16-19]。
當(dāng)點云和設(shè)計模型的多邊形進(jìn)行比較時,傳統(tǒng)的方法是使用每個測量點和多邊形的邊進(jìn)行比較,確定測量點在多邊形內(nèi)和多邊形外,然而這種方法運(yùn)算效率低,不能滿足實時性要求高的動態(tài)測量中[20-22]。
本文在分析多傳感器結(jié)構(gòu)特性的基礎(chǔ)上,提出一種基聚類分析的大型復(fù)雜精密零件自動化測量方法,不僅在存在重影以及環(huán)境光影響的情況下精確的提取光條的中心,去除光條上的噪聲,還可以對光條曲線進(jìn)行光順,移動大型精密部件,可以實現(xiàn)高速高精度的大型零部件產(chǎn)品的三維數(shù)據(jù)采集,基于射線法分析大型零部件產(chǎn)品的測量點云和設(shè)計模型之間的位置關(guān)系,反饋給控制系統(tǒng),用于后續(xù)的裝配操作,測量效率高,精度高。
將四組視覺傳感器單元(包括CCD相機(jī)和激光器)分別豎直放置于精密零件兩側(cè),垂直于精密零件中軸線,安裝在精密零件兩側(cè)的立柱上;其中兩組視覺傳感器單元放置于精密零件一側(cè)立柱上,另外兩組視覺傳感器單元放置于精密零件另一側(cè)立柱上。
首先利用激光跟蹤儀建立基坐標(biāo)系統(tǒng),然后對每個視覺傳感器單元進(jìn)行相機(jī)參數(shù)標(biāo)定、光平面參數(shù)方程標(biāo)定以及全局標(biāo)定,最終得到相機(jī)的內(nèi)參矩陣、相機(jī)到基坐標(biāo)系下的全局標(biāo)定矩陣以及激光平面在基坐標(biāo)系下的平面方程,完成系統(tǒng)使用以及測量前的預(yù)處理,系統(tǒng)結(jié)構(gòu)如圖1所示。
圖1 系統(tǒng)結(jié)構(gòu)圖
利用激光跟蹤儀建立測量系統(tǒng)的基坐標(biāo)系,坐標(biāo)系的原點建立在精密零件滑動機(jī)構(gòu)的中點,Z軸為精密零件運(yùn)行的軸線方向,Y軸為垂直于Z軸平面,滑動機(jī)構(gòu)頂端連線方向,X軸可以用右手法則確定。
設(shè)計帶多個圓孔的空間靶標(biāo),如圖2所示,利用三坐標(biāo)測量儀對設(shè)計靶標(biāo)進(jìn)行預(yù)先精確測量,確定所有孔的上表面孔中心的三維坐標(biāo)。
圖2 全局標(biāo)定示意圖
利用多對同名點在相機(jī)坐標(biāo)系和激光跟蹤儀坐標(biāo)系下的坐標(biāo)值,確定無約束最優(yōu)目標(biāo)函數(shù)F,求解第j個相機(jī)坐標(biāo)系和激光跟蹤儀坐標(biāo)系之間的轉(zhuǎn)換矩陣H:
(1)
其中:n為參與標(biāo)定的點個數(shù),Pwi=(xwi,ywi,zwi,l)為第i個點在全局坐標(biāo)系下的齊次坐標(biāo),Pci=(xci,yci,zci,l)為第i個點在相機(jī)坐標(biāo)系下的齊次坐標(biāo)。j?[1,4]。該函數(shù)需要滿足正交約束條件,即RTR=I,其中R為坐標(biāo)系轉(zhuǎn)換矩陣的旋轉(zhuǎn)矩陣,I是3×3的單位向量。
利用激光跟蹤儀采集空間靶標(biāo)上的孔位中心的所有點,相機(jī)采集空間靶標(biāo)圖像,得到對應(yīng)的孔位中心點在相機(jī)坐標(biāo)系下的三維坐標(biāo)點,最小化目標(biāo)函數(shù)F,可以得到第j個相機(jī)坐標(biāo)系和激光跟蹤儀坐標(biāo)系之間的轉(zhuǎn)換矩陣Hj,因此利用上述激光跟蹤儀建立的基坐標(biāo)系,可以將相機(jī)坐標(biāo)系轉(zhuǎn)換到基坐標(biāo)系下。同理,可以計算出所有相機(jī)坐標(biāo)系到基坐標(biāo)系之間的轉(zhuǎn)換矩陣,完成全局標(biāo)定。
嚴(yán)格控制所有相機(jī)進(jìn)行同步采集,將采集到的圖像進(jìn)行圖像處理,剔除采集過程中存在的噪聲以及激光重影,得到線激光在每個相機(jī)坐標(biāo)系的三維坐標(biāo);
將相機(jī)采集到的彩色三通道圖像分解為R(x,y)、G(x,y)、B(x,y) 3個單通道的圖像,然后增強(qiáng)圖像的紅色分量,得到增強(qiáng)后的圖像I(x,y):
I(x,y)=R(x,y)-G(x,y)-B(x,y)
(1)
其中:x,y為像素位置。
對圖像中的每一列進(jìn)行上述處理,得到所有滿足條件的激光中心線作為光條大致的候選中心點,使用漫水填充算法,利用8連通區(qū)域連接和候選種子點相連通的激光條中心點,將斷續(xù)的光條進(jìn)行連接。
然后以長度作為閾值,由于光條的重影亮度比正常的光條亮度暗,環(huán)境光雖然有時候較強(qiáng)烈,但是環(huán)境光影響局部區(qū)域,因此按照光條連通區(qū)域最長的區(qū)域作為最終的光條中心點對光條區(qū)域進(jìn)行篩選,得到最終的精確光條中心點。利用激光三角測量原理,可以將圖像的光條中心點轉(zhuǎn)換為相機(jī)坐標(biāo)下的的三維空間點。
原始掃描數(shù)據(jù)往往包含大量不屬于實際光條的稀疏離群點,這些離群點是由掃描過程中的環(huán)境因素如表面反射率性質(zhì)不好、技術(shù)缺陷或人為錯誤造成的,存在噪聲的激光條會導(dǎo)致點云拼接、點云光順等后處理的失敗。因此,在點云數(shù)據(jù)融合和光順前需要對點云的離群點進(jìn)行剔除。為了更有效地反映物體的真實形貌,本文采用統(tǒng)計濾波方法對三維數(shù)據(jù)點云進(jìn)行雜點的去除[23]。首先對每一個點Pi的K鄰域做一個統(tǒng)計分析處理,計算該點到所有臨近點的平均距離,得到均值μ和標(biāo)準(zhǔn)差σ,判定鄰域點是否為噪點的條件如下所示:
(2)
其中:Pj為Pi的其中一個鄰域點,‖PiPj‖ 表示Pi和Pj兩點之間的距離,當(dāng)Pj為噪點時,將其從點云數(shù)據(jù)中移除。
給定空間三維數(shù)據(jù),解決這一問題最常見的方法是確定一個近似樣條曲線可以最小化一個確定的柔順函數(shù)。根據(jù)給定的空間點集推導(dǎo)出一個基于數(shù)據(jù)直接操作的自動局部整流罩過程。利用曲線的幾何性質(zhì),提出了一種有效的局部迭代方法,其流程如圖3所示。
(3)
圖3 曲線柔順流程圖
將相鄰兩個相機(jī)坐標(biāo)系下的坐標(biāo)點云轉(zhuǎn)換到基坐標(biāo)系下,理論上局部重疊的區(qū)域點云應(yīng)該是重合的,然而由于在全局標(biāo)定過程中不可避免地存在誤差,導(dǎo)致重疊區(qū)域的重疊點云不完全重合。如果只使用一個相機(jī)坐標(biāo)系下的點云數(shù)據(jù),會導(dǎo)致點云出現(xiàn)嚴(yán)重的分層現(xiàn)象,因此在兩組點云融合之前,需要確定相鄰兩個相機(jī)測量數(shù)據(jù)之間的對應(yīng)點。
根據(jù)精密零件自身的特性以及設(shè)備安裝的位置,光條的點云在精密零件輪廓上,在基坐標(biāo)系下呈單調(diào)性分布,因此只需要確定L1點云序列的Y值最小點以及L2點云序列的Y值最大點,在這個區(qū)間內(nèi)的點,認(rèn)為是重疊區(qū)域的重疊點云。
借助于采樣一致性初始配準(zhǔn)算法(SAC-IA)進(jìn)行粗匹配,再利用迭代最近點匹配算法對點云進(jìn)行精細(xì)匹配,得到重疊部分區(qū)域點云之間的轉(zhuǎn)換矩陣,即旋轉(zhuǎn)矩陣R和平移向量T。利用轉(zhuǎn)換矩陣將可以對點云重疊部分的數(shù)據(jù)進(jìn)行融合。
在對大型產(chǎn)品的結(jié)果進(jìn)行分析時,需要確定大型產(chǎn)品的測量結(jié)果和設(shè)計模型的截面輪廓的對應(yīng)關(guān)系。
判斷某個斷面的點云是否超限,相當(dāng)于判斷是否有點云在限界標(biāo)準(zhǔn)輪廓范圍內(nèi),采用射線法通過判斷該點水平線和多邊形交點的奇偶性即可判定是否超界。
假設(shè)給定的設(shè)計模型輪廓為多邊形M,該多邊形是由n個頂點{p1={x1,y1},p2={x2,y2},…,pn={xn,yn}}構(gòu)成的。 為測量模型輪廓點云的其中一個點q的標(biāo)識位,以q為端點向右畫一條射線R,然后對整個輪廓多邊形M的每條邊進(jìn)行遍歷,標(biāo)記射線R與所有邊相交次數(shù)的總和c,依據(jù)總和的奇偶性判定測量點云和設(shè)計模型輪廓的關(guān)系:
(4)
如果相交次數(shù)總和c為奇數(shù),則說明q點位于輪廓內(nèi)部,c為偶數(shù),則說明q點位于輪廓外部。
如果超界進(jìn)行報警,軟件界面進(jìn)行超界位置顯示。
為驗證本方法的有效性,采用basler相機(jī),幀率為30 fps,分辨率為300 w像素,配置computar公司的12 mm鏡頭,激光器采用波長為650 nm的紅色激光,抗干擾能力強(qiáng),功率為0.5 mW,仿真實驗在vs2013軟件下配合opencv算法庫實現(xiàn)。
利用高精度激光跟蹤儀對一個三級階梯的標(biāo)準(zhǔn)塊進(jìn)行三維尺寸計量,標(biāo)準(zhǔn)塊選用的材質(zhì)為1018鋼,三級階梯的高度都為5 mm,精度為0.003 mm。
測量時標(biāo)準(zhǔn)塊放置于相機(jī)和激光器焦距范圍視野中央,發(fā)射線激光到標(biāo)準(zhǔn)塊表面,計算機(jī)對相機(jī)采集的圖像進(jìn)行處理,三維重建點云,獲得階梯標(biāo)準(zhǔn)塊的實測階梯差。D為相機(jī)距標(biāo)準(zhǔn)塊的距離。
表1~3可以看出,測量精度隨著物體距相機(jī)的距離增大而降低,因此實驗中使用D≈600,在該位置處對零部件進(jìn)行三維掃描,可以看出,系統(tǒng)測量精度大于0.25 mm,滿足實際測量需求。
表1 絕對差值(D≈600 mm)
表2 絕對差值(D≈800 mm)
將多個相機(jī)進(jìn)行圖像采集和圖像處理,利用標(biāo)定的線結(jié)構(gòu)光光平面參數(shù)以及全局標(biāo)定矩陣,借助于激光三角原理,
表3 絕對差值(D≈1 000 mm)
重建光條中心點的三維坐標(biāo),將光條中心點的三維坐標(biāo)轉(zhuǎn)換到全局坐標(biāo)系下,對點云進(jìn)行去噪以及點云平滑,可以得到在全局坐標(biāo)系下的精密零件表面三維測量點云數(shù)據(jù),如圖4所示,點云處理可以實現(xiàn)模型的三維重建。
圖4 三維重建模型
針對大型精密部件測量工序復(fù)雜,測量精度要求高的需求,本文提出一種魯棒的線激光檢測方法,由于存在環(huán)境光以及重影的影響,采用聚類和灰度重心法對結(jié)構(gòu)光進(jìn)行中心定位,利用三角測量原理和全局標(biāo)定矩陣,將測量的點云數(shù)據(jù)轉(zhuǎn)換到基坐標(biāo)系下,然后將獲取到的三維點云數(shù)據(jù)進(jìn)行去噪,去除由掃描過程中的技術(shù)缺陷而引起的離散噪聲,提升點云的質(zhì)量。點云融合可以將相鄰相機(jī)重疊區(qū)域的光條進(jìn)行合并,有效地減少數(shù)據(jù)量,并且有效地修正激光線檢測過程中帶來的誤差以及產(chǎn)生的鋸齒問題,利用空間點云光順方法,將點云進(jìn)行光順平滑處理,可以真實反映被測量物體的三維信息。最后采用射線法,將測量物體截面上的點云數(shù)據(jù)進(jìn)行分析處理,給操作者提供預(yù)警信號。實驗表明,本文提出的方法可以快速有效地檢測存在環(huán)境光影響下的激光光條中心,重建大型精密零部件的三維形貌。