唐昀超，黃礦裕，陳明猷

（1.仲愷農(nóng)業(yè)工程學院，廣州 510225；2.華南農(nóng)業(yè)大學，廣州 510642）

0 引言

傳統(tǒng)的產(chǎn)品開發(fā)過程無法用統(tǒng)一且嚴密的數(shù)學模型來描述復雜的曲面，設計者很難構(gòu)造出既能滿足產(chǎn)品的功能要求又能吸引消費者注意的具有復雜曲面的產(chǎn)品外觀[1]。因此，傳統(tǒng)的產(chǎn)品開發(fā)方式已逐漸滿足不了新的市場需求。曲面重構(gòu)技術作為一種產(chǎn)品設計技術再現(xiàn)過程，能夠獲取已有產(chǎn)品的復雜外形，便于設計者在此基礎上制作出功能相近，但又不完全一樣的產(chǎn)品。通過激光掃描儀、結(jié)構(gòu)光源轉(zhuǎn)換儀或者立體視覺等3D掃描技術獲取物體的表面數(shù)據(jù)，再進一步生成具有拓撲結(jié)構(gòu)信息的網(wǎng)格或CAD模型。

目前，國內(nèi)外學者對三維重構(gòu)技術做了大量研究。Hugues等提出了一種用于無序點云曲面重構(gòu)的算法[2]，首先確定每個點的切平面并保證切平面的全局方向一致，然后通過估計符號距離函數(shù)確定零等值面，該等值面即為重構(gòu)的曲面，該方法能夠平滑地擬合出目標大致的曲面模型。Jules等提出一種基于先驗模型的曲面重構(gòu)算法[3]，通過事先建立目標的先驗模型，然后與未遮擋的點云進行泊松曲面重構(gòu)而成的曲面模型進行融合，解決了在遮擋情況下封閉點云的曲面重構(gòu)問題，因為該方法是基于先驗模型對曲面的孔洞區(qū)域進行填補，所以不適用于具有復雜曲面的模型。曹詩卉等人提出一種利用凸最小化泊松曲面改進的重建算法[4]，能夠處理通過非均勻采樣或者有噪聲的有序點云，得到光順的曲面。賴喜德等提出了水輪機葉片曲面重構(gòu)的技術[5]，采用基于NURBS的曲面重構(gòu)方法，為具有自由曲面特征的葉片類部件提供了曲面重構(gòu)流程。孫智慧等提出基于點云數(shù)據(jù)的植物葉片曲面重構(gòu)方法[6]，該方法首先對點云數(shù)據(jù)進行精簡，然后采用Delaunay三角剖分方法生成初始網(wǎng)格曲面，再對網(wǎng)格曲面進行優(yōu)化處理。基于上述研究，本文將以齒輪作目標，利用雙目視覺從多個視角獲取圖像，生成目標的三維點云數(shù)據(jù)，并對點云進行去噪、配準和平滑處理，最終反求齒輪的三維曲面模型。

1 試驗裝置

本文以齒輪作為研究對象，搭建了如圖1所示的工作臺，以獲取齒輪各個方位的表面輪廓。該工作臺由雙目相機、支架、旋轉(zhuǎn)臺和工件組成。其中雙目相機為維視工業(yè)相機，型號為MV-EM200C，分辨率為1 600×1 200。支架由導軌搭建而成，高度可調(diào)。旋轉(zhuǎn)臺臺面為白色，方便圖像處理階段分割出齒輪區(qū)域，旋轉(zhuǎn)臺轉(zhuǎn)速由PC端控制。雙目相機之間平行排布，基線距離為60 mm，2臺相機均通過網(wǎng)口線與計算機進行通訊。相機位于旋轉(zhuǎn)臺中心左上方，以便獲取齒輪圓周部分的表面數(shù)據(jù)，相機與旋轉(zhuǎn)臺之間的距離約為400 mm。每獲取完一個視角下的圖像數(shù)據(jù)，控制旋轉(zhuǎn)臺就旋轉(zhuǎn)一定角度（約為60 ）繼續(xù)獲取圖像，直至旋轉(zhuǎn)一周，然后翻轉(zhuǎn)齒輪再獲取一次圖像，以便獲取齒輪背面的表面數(shù)據(jù)。

圖1 試驗裝置圖

2 三維點云的獲取

2.1 圖像矯正

因相機鏡頭存在畸變，致使采集圖像失真，影響到測量精度，需通過單目標定來矯正失真的圖像。雙目相機存在著極線對齊的問題，需要利用立體標定解決，通過約束極線對齊，可降低立體匹配復雜度。本文通過張正友標定法完成相機標定并獲得矯正參數(shù)[7]。圖2顯示的是矯正后的齒輪圖像。

圖2 矯正圖像

2.2 圖像預處理

圖像預處理階段包括了灰度化、截取ROI（感興趣區(qū)域）、高斯濾波與閾值分割四個部分。為了減少計算量，截取齒輪所在的部位作為ROI進行后續(xù)圖像處理。為了進一步減少計算量，使用加權灰度化方法將所得的三通道ROI轉(zhuǎn)化成灰度圖像，如式（1）：

式中 （i，j）表示像素點的位置，R、G和B表示紅、綠和藍三種顏色分量的值。

采集圖像的過程中存在著噪聲，因此利用高斯濾波對圖像進行平滑。高斯濾波使用高斯卷積核遍歷圖像，并在卷積核所在的窗口對像素值進行卷積操作，最后把卷積結(jié)果賦給錨點，從而完成濾波操作。若采用的卷積核大小為(2k+1)×(2k+1)，則該核上第i行，第j列的值為：

由于目標周圍的背景大致為白色，為了在立體匹配過程中減少匹配計算量，將左圖像的背景變?yōu)楹谏?，這樣左右圖像的背景部分由于差異太大將無法匹配。首先采用最大類間方差法進行閾值分割[8]，再根據(jù)式（3）將左圖的白色背景部分變?yōu)楹谏Ｗ罱K處理效果如圖3所示。

圖3 圖像處理效果圖

2.3 圖像重構(gòu)

該階段包括立體匹配和獲取三維點云兩個部分，首先通過立體匹配得到左右圖像中匹配點的視差，再利用三角測量原理計算得到三維點的坐標，從而獲取三維點云。立體匹配采用SGBM算法，該算法包含了匹配代價計算、匹配代價聚合、視差的優(yōu)化計算和視差精化四個步驟[9]。通過立體匹配，得到如圖4所示的視差圖，再由式（4）獲取三維點云。

圖4 視差圖

圖5 初始三維點云

3 曲面的獲取

3.1 點云去噪

得到的初始三維點云，不可避免地會出現(xiàn)噪點，如圖5所示。若處理不當，會直接影響最終曲面重構(gòu)的效果。為了獲取較好的去噪效果，本文采用直通濾波、半徑濾波、體素柵格濾波和統(tǒng)計濾波五種濾波方法進行組合去噪。直通濾波算法是通過設定坐標范圍來過濾點云，能夠快速粗略地去除離群點，達到粗處理的目的[10]。半徑濾波算法以搜索點為中心拓展出半徑為R的球體，統(tǒng)計落在球體內(nèi)的點云個數(shù)，當個數(shù)小于給定的閾值K時，則判斷該搜索點為離群點并刪除。體素柵格濾波算法是根據(jù)點云形成三維體素柵格，然后利用每個體素柵格中所有三維點的重心來近似代替體素柵格內(nèi)所有的三維點，通過這樣的采樣方式來得到濾波的效果。統(tǒng)計濾波算法是利用對搜索點與鄰域點集之間的距離統(tǒng)計來判斷該點是否為噪點。濾波后的點云如圖6所示。

圖6 濾波效果圖

3.2 點云配準與光滑

上述所獲取的點云只是在一個相機視角下得到的，無法描述完整物體表面，因此，需要在多個視角下獲取物體不同位置的表面，然后進行配準，拼接出完整的目標。本文將使用ICP（迭代最近點）算法對點云進行配準。ICP是基于自由形態(tài)曲面的配準算法[11]，由于具有較高的準確度和效率，是點云配準的主流算法。配準的效果如圖7所示。

圖7 配準效果圖

配準而成的點云不可避免的存在大量的不規(guī)則重疊，若直接進行曲面重建，得到的曲面將會不光滑，呈鋸齒狀。為了得到光順的表面，采用移動最小二乘法對點云進行平滑化處理[12,13]。最終平滑后的效果如圖8所示。

圖8 平滑效果圖

3.3 曲面重構(gòu)

經(jīng)上述步驟處理后的數(shù)據(jù)仍是離散點云，缺乏拓撲關系，因此本文采用貪婪投影三角化算法對點云數(shù)據(jù)進行曲面重構(gòu)。該方法將三維數(shù)據(jù)點按法線方向投影到平面，然后對平面上的點作平面三角化，得到各點之間的連接關系。三角化過程采用的是基于Delaunay的三角化算法[14]，先在平面上進行區(qū)域生長，通過選擇一個三角片作為初始曲面，不斷地拓展邊界，最后形成一張完整的曲面。

具體過程為：①找出任意P點的k鄰域，當鄰域足夠小的時候，曲面具有局平性，鄰域內(nèi)的所有三維數(shù)據(jù)點和P點在過P點的切平面上必定存在與它們一一對應的唯一投影點；②求出過P點的切平面以及鄰域內(nèi)所有點和P點在該切平面上的投影點。③對投影點進行平面三角剖分。最終曲面重構(gòu)的結(jié)果如圖9所示。

圖9 曲面重構(gòu)效果圖

4 結(jié)論

本文提出一種基于雙目視覺的齒輪曲面重構(gòu)算法。首先分別對不同視角下的齒輪雙目圖像進行處理，包括矯正圖像、閾值分割和立體匹配等步驟，完成圖像重構(gòu)獲得初始三維點云。然后對三維點云進行組合濾波去噪，達到去除噪點和精簡點云數(shù)據(jù)量的目的。之后將多片點云進行配準，拼接成完整的齒輪點云模型，然后利用移動最小二乘法進行平滑，去除重疊的點云數(shù)據(jù)，最終采用貪婪三角化算法將點云轉(zhuǎn)換成帶有拓撲關系的曲面模型。上述算法能夠較好地重構(gòu)出目標的曲面模型，但仍存在著不足。重構(gòu)得到的曲面模型不具備水密性，對于點云數(shù)據(jù)缺失部分，對應的曲面模型也存在孔洞。針對這一問題，在后續(xù)研究中將結(jié)合泊松曲面算法進行完善，得到封閉光滑的曲面模型。