韓 成，秦貴和，宮 宇，張 超，薛耀紅

（1.吉林大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院，長春 130012；2.長春理工大學(xué) 計(jì)算機(jī)科學(xué)技術(shù)學(xué)院，長春 130022；3.吉林交通職業(yè)技術(shù)學(xué)院 外語系，長春 130021）

0 引 言

一般情況下，結(jié)構(gòu)光三維重構(gòu)方法主要分為系統(tǒng)參數(shù)標(biāo)定、結(jié)構(gòu)光編解碼和特征點(diǎn)三維信息的計(jì)算3個部分。其中，系統(tǒng)參數(shù)標(biāo)定部分主要是通過最小二乘法計(jì)算攝像機(jī)和投影機(jī)的內(nèi)外參數(shù)；特征點(diǎn)三維信息的計(jì)算大多采用針孔模型和三角形相似原理；對于結(jié)構(gòu)光的編碼方法來說，利用各種編碼策略可以使像素點(diǎn)本身帶有先驗(yàn)編碼信息，編碼策略選擇的合適與否影響著匹配算法的復(fù)雜性和三維重構(gòu)的準(zhǔn)確性。對結(jié)構(gòu)光編碼方法的研究已有很多，其中時(shí)間編碼中的線結(jié)構(gòu)光［1-2］方法不能快速實(shí)時(shí)地對物體進(jìn)行重構(gòu)，但它的編碼原理簡單、易于實(shí)現(xiàn)；現(xiàn)有的單幅空間編碼［3］方法雖然能夠?qū)ξ矬w實(shí)現(xiàn)實(shí)時(shí)重構(gòu)，但通過尋找角點(diǎn)［4］或條紋邊界［5］的方法得到的特征點(diǎn)坐標(biāo)信息不夠精確，并且大多數(shù)情況下是使用彩色條紋編碼［6-8］的形式，并不適用于無約束的三維信息解算。對于解碼方法，張超等［9］提出一系列圖像預(yù)處理方法，可以更加準(zhǔn)確地提取較細(xì)格雷碼的條紋；于曉洋等［10］采用一種基于邊緣導(dǎo)向的亞像素定位技術(shù)提取格雷碼條紋邊緣；Fechteler等［6］提出了自適應(yīng)彩色分類方法和捕獲3D模型的方法，解決了環(huán)境光、顏色、人臉特征對重構(gòu)質(zhì)量的影響，提高了系統(tǒng)的魯棒性。Brink等［11］根據(jù)圖中已記錄條紋的潛在關(guān)聯(lián)性和鄰接性，提出了最大生成樹算法，使用該算法可以得到比以往方法更精確、更可靠的重構(gòu)結(jié)果?？琢罡坏龋?2］利用序列的窗口特性和空間鄰域編碼思想對各條紋邊界進(jìn)行標(biāo)識，采用全局優(yōu)化思想對條紋邊界進(jìn)行最優(yōu)鄰域匹配，不需要假定物體全表面單調(diào)，對其連續(xù)性也不做嚴(yán)格要求。Zhang等［13］引入了多通道動態(tài)規(guī)劃思想，克服了圖像中的條紋遮擋問題，但常會出現(xiàn)對局部單調(diào)區(qū)域之間的邊界集的解碼誤差。

在以往的研究中，彩色結(jié)構(gòu)光編碼方法多用于空間編碼和直接編碼，而且大多數(shù)都能用于動態(tài)物體的測量，但由于受到CCD攝像機(jī)的顏色保真度及分辨率的影響而無法得到理想的效果和精度，不容易找到清晰的、無二義性的對應(yīng)，對三維物體的測量精度不高，很難用于物體的精確測量。一般情況下，運(yùn)用的顏色多于3種，容易出現(xiàn)顏色混淆不易識別的問題，而且通常采樣密度和分辨率都不高。彩色結(jié)構(gòu)光的編碼和解碼技術(shù)在結(jié)構(gòu)光重構(gòu)的數(shù)據(jù)獲取中起著至關(guān)重要的作用，進(jìn)而影響著系統(tǒng)的測量精度、速度和可靠性。本文提出了一種具有拓?fù)浣Y(jié)構(gòu)和內(nèi)部結(jié)構(gòu)相關(guān)性的彩色結(jié)構(gòu)光編碼方法及其解碼方法。

1 系統(tǒng)參數(shù)的標(biāo)定

為了計(jì)算得到特征點(diǎn)的三維信息，需要對系統(tǒng)的參數(shù)進(jìn)行標(biāo)定。標(biāo)定的精度和準(zhǔn)確性是至關(guān)重要的，每一個系統(tǒng)參數(shù)的標(biāo)定誤差都將對三維點(diǎn)云計(jì)算的結(jié)果產(chǎn)生極大影響。對系統(tǒng)參數(shù)的標(biāo)定主要分為攝像機(jī)的內(nèi)外參數(shù)標(biāo)定和投影機(jī)的內(nèi)外參數(shù)標(biāo)定兩個部分［14］。在進(jìn)行標(biāo)定的過程中，需要將攝像機(jī)和投影機(jī)二者之間的關(guān)系統(tǒng)一在同一個坐標(biāo)系下，因此，在進(jìn)行系統(tǒng)參數(shù)標(biāo)定之前必須建立一個世界坐標(biāo)系。

首先將具有特定尺寸、已知特征的平面標(biāo)定圖像貼附于參考平面板上，然后將此參考平面板固定于精密移動平臺上，并保證平面標(biāo)定圖像能夠被攝像機(jī)完整地?cái)z取，同時(shí)投影機(jī)也能夠?qū)⑵矫鏄?biāo)定圖像完整地投射到參考平面板上。如圖1所示，在參考平面板沒有移動時(shí)，設(shè)平面標(biāo)定圖像所在的平面為世界坐標(biāo)系的XY平面，平面標(biāo)定圖像中沿著特征點(diǎn)向右為世界坐標(biāo)系的X軸正方向，沿著特征點(diǎn)向下為世界坐標(biāo)系的Y軸正方向，垂直于平面標(biāo)定圖像向內(nèi)為世界坐標(biāo)系的Z軸正方向，平面標(biāo)定圖像中最左上特征點(diǎn)位置為世界坐標(biāo)系的原點(diǎn)。

圖1 世界坐標(biāo)系的確定Fig.1 World coordinate system

1.1 攝像機(jī)參數(shù)標(biāo)定

為了減少標(biāo)定時(shí)間，提高標(biāo)定精度，通常假定攝像機(jī)的成像模型是小孔成像［14］。在對攝像機(jī)內(nèi)外參數(shù)進(jìn)行標(biāo)定時(shí)，對攝像機(jī)的成像畸變進(jìn)行校正。首先通過線性的針孔模型計(jì)算得到攝像機(jī)的初始內(nèi)參，然后使用最大似然法對攝像機(jī)進(jìn)行非線性求精，得到考慮畸變系數(shù)的優(yōu)化內(nèi)參。在多個方向上移動參考平面板，每次移動完成后都由攝像機(jī)攝取一幅標(biāo)定圖像。為了盡可能地提高移動過程中的精確度，使用精密移動平臺以保證每次移動的距離都是預(yù)先設(shè)定的，這樣便可以知道平面標(biāo)定圖像上每一個特征點(diǎn)在世界坐標(biāo)系中的確切坐標(biāo)值。

標(biāo)定圖像攝取完成后，需要對攝取的標(biāo)定圖像進(jìn)行畸變校正。使用如圖2所示的極坐標(biāo)模型，假設(shè)無畸變條件下的一點(diǎn)P0（x0，y0）在畸變圖像中的位置為P1（x1，y1），P0和P1到原點(diǎn)O的距離分別為r0和r1，且在不考慮切向畸變的情況下P0和P1在同一條由原點(diǎn)O發(fā)出的射線上，則P0和P1的關(guān)系為

式中：r0為無畸變點(diǎn)到原點(diǎn)O的距離；r1為畸變點(diǎn)到原點(diǎn)O的距離；ai為借助攝像機(jī)已經(jīng)求出的系數(shù)。

圖2 畸變校正的極坐標(biāo)模型Fig.2 Polar coordinate model of distortion correction

根據(jù)此模型進(jìn)行畸變校正的效果對比圖如圖3所示。

圖3 畸變校正效果Fig.3 Effect of distortion correction

對攝取的標(biāo)定圖像進(jìn)行畸變校正而得到每一幅平面標(biāo)定圖像中各個特征點(diǎn)在標(biāo)定圖像中的圖像坐標(biāo)，然后根據(jù)公式

計(jì)算出特征點(diǎn)在世界坐標(biāo)系中的世界坐標(biāo)和在圖像坐標(biāo)系中的圖像坐標(biāo)之間的單應(yīng)性變換矩陣MC。式（4）中，（ui，vi）為攝像機(jī)圖像坐標(biāo)系中特征點(diǎn)的二維圖像坐標(biāo)；MCij為MC矩陣中第i行第j列的元素值；（Xi，Yi，Zi）為世界坐標(biāo)系中特征點(diǎn)的三維世界坐標(biāo)。

1.2 投影機(jī)參數(shù)標(biāo)定

對攝像機(jī)參數(shù)標(biāo)定完成后，便可以根據(jù)MC矩陣的值對投影機(jī)的參數(shù)進(jìn)行標(biāo)定。從攝像機(jī)和投影機(jī)成像模型的角度看，兩者都是使用針孔成像模型，但又存在著一個根本的區(qū)別，即一個是攝取圖像，一個是投射圖像。因此，在對投影機(jī)進(jìn)行標(biāo)定時(shí)將需要借助于攝像機(jī)攝取圖像的能力。

首先通過投影儀將平面標(biāo)定圖像投射到參考平面板上，并由攝像機(jī)對投射到參考平面板上的平面標(biāo)定圖像進(jìn)行攝取。通過線性的針孔模型計(jì)算得到投影機(jī)的初始內(nèi)參，然后使用最大似然法對投影機(jī)進(jìn)行非線性的求精，得到考慮畸變系數(shù)的優(yōu)化內(nèi)參。根據(jù)投影機(jī)的畸變系數(shù)對待投影的標(biāo)定圖像進(jìn)行畸變還原，使待投射的標(biāo)定圖像是具有畸變的，以保證投射出的畸變圖像通過投影機(jī)的投射后將原有的畸變相互抵消掉，使最終投射到標(biāo)定板上的標(biāo)定圖像能夠?qū)?shù)的標(biāo)定結(jié)果起到積極作用，使參數(shù)的標(biāo)定結(jié)果更加精確。如果不對標(biāo)定圖像進(jìn)行畸變還原，則將使得單應(yīng)性變換矩陣MP中的外參值出現(xiàn)誤差?；冞€原的過程與畸變校正的過程是相反的，該過程可由式（2）進(jìn)行描述。

將處理后的標(biāo)定圖像再次投射到參考平面板上，并由攝像機(jī)攝取，根據(jù)公式

有效的注視點(diǎn)要求視網(wǎng)膜中央凹對準(zhǔn)注視目標(biāo)超過200ms,考慮到視覺任務(wù)的處理耗時(shí),操作者的視覺注意處于用戶界面局部區(qū)域并進(jìn)行任務(wù)加工時(shí),注意轉(zhuǎn)移速度慢,注視點(diǎn)停留時(shí)間可能會更長,因此本研究認(rèn)為注視點(diǎn)停留時(shí)間大于某一閾值也可以反映操作者正在關(guān)注注視點(diǎn)所在位置區(qū)域.

求解沒有進(jìn)行畸變還原的標(biāo)定圖像中的特征點(diǎn)與攝像機(jī)攝取的標(biāo)定圖像中的特征點(diǎn)的單應(yīng)性變換矩陣MP。式（5）中，MPij為MP矩陣中第i行第j列的元素值；由式（4）計(jì)算得出特征點(diǎn)在世界坐標(biāo)系下的三維坐標(biāo)為 （Xi，Yi，Zi）。

2 彩色結(jié)構(gòu)光編碼方法

結(jié)構(gòu)光的編解碼技術(shù)是一種主動式地給待測三維實(shí)體進(jìn)行特征點(diǎn)標(biāo)記的有效方法。該方法可以使三維實(shí)體具有更易識別和控制的特征點(diǎn)信息。對于結(jié)構(gòu)光編碼來說，主要包含結(jié)構(gòu)光的編碼設(shè)計(jì)，結(jié)構(gòu)光特征點(diǎn)信息的提取，結(jié)構(gòu)光特征點(diǎn)信息的匹配3個重要組成部分。為了求解出指定特征點(diǎn)的三維坐標(biāo)信息，必須準(zhǔn)確地標(biāo)識出每一對特征點(diǎn)之間的匹配關(guān)系，因此結(jié)構(gòu)光的編解碼技術(shù)便成了確定特征點(diǎn)信息的重要手段。

對所有光束進(jìn)行編碼是結(jié)構(gòu)光檢測技術(shù)的核心，好的編碼方法可以減少解碼的時(shí)間，提高解碼的精度。目前已有多種結(jié)構(gòu)光編碼方法，如二進(jìn)制編碼、格雷編碼、基于隨機(jī)序列的彩色編碼、直接彩色編碼等［1，3，6］。但如果彩色條紋間隔很近，由于彩色條紋間的干擾，很難準(zhǔn)確地提出中心彩色條紋或標(biāo)定條紋的顏色。本文提出一種具有拓?fù)浣Y(jié)構(gòu)和局部結(jié)構(gòu)化的結(jié)構(gòu)光編碼方法。這種編碼方法具有以下4個優(yōu)點(diǎn)：①編碼得到的結(jié)構(gòu)圖像能夠滿足大尺寸、高精度的要求。②編碼的結(jié)構(gòu)化程度更高，使得編碼圖像的解碼正確率更高，編碼圖像的抗干擾性更強(qiáng)。③由單幅編碼圖像便可以得到整個空間的特征點(diǎn)信息，實(shí)時(shí)性強(qiáng)。④該編碼方法對接收結(jié)構(gòu)光編碼圖像的曲面沒有嚴(yán)格的限制。

在此編碼方法中，使用了8種顏色對圖像進(jìn)行編碼，其中的6種彩色顏色便可以有720種不同的5元素組合，但在圖像編碼的過程中卻不能將這些顏色進(jìn)行任意的組合，因?yàn)榛叶戎迪嘟念伾ㄈ缂t色與品紅色）進(jìn)行相鄰組合將對解碼工作造成阻礙而不能正確地進(jìn)行解碼，如圖4所示。其中顏色塊的顏色分別為紅色、綠色、藍(lán)色、黑色、青色、品紅色、黃色、白色。將所述的8種顏色塊分為兩組，一組是由紅色、綠色、藍(lán)色、品紅、黃色、青色6種顏色塊組成，稱為彩色元組；另一組是由黑色、白色兩種顏色塊組成，稱為黑白元組。根據(jù)彩色元組中6種彩色顏色塊的灰度值的不同將彩色元組劃分為兩組，一組是由紅色、藍(lán)色、品紅色3種彩色顏色元素組成，稱為高值彩色元組，另一組是由綠色、青色、黃色3種顏色元素所組成，稱為低值彩色元組。

使用黑白元組中的元素交替排列組成如圖5（a）的3×3方形區(qū)域，稱為白微區(qū)域；使用黑白元組中的元素按照與白微區(qū)域相同位置塊的顏色相反組成如圖5（b）的區(qū)域，稱為黑微區(qū)域；利用彩色元組中的6種彩色顏色生成微區(qū)域的一層彩色邊界，彩色邊界中與低值彩色元組中元素相鄰的顏色塊為在高值彩色元組中隨機(jī)選取的元素，與高值彩色元組中元素相鄰的顏色塊為在低值彩色元組中隨機(jī)選取的元素。如果彩色邊界的左上方第一個顏色塊為高值彩色元組中的元素，稱為高值彩色原子邊界，如圖6（a）所示；如果彩色邊界的左上方第一個顏色塊為低值彩色元組中的元素，稱為低值彩色原子邊界，如圖6（b）所示。使用低值彩色原子邊界和白微區(qū)域構(gòu)成5×5的方形區(qū)域，稱為低值原子區(qū)域，如圖7（a）所示；使用高值彩色原子邊界和黑微區(qū)域構(gòu)成5×5的方形區(qū)域，稱為高值原子區(qū)域，如圖7（b）所示。

圖5 區(qū)域結(jié)構(gòu)化編碼圖像Fig.5 Regional structured coding image

圖6 彩色邊界圖像Fig.6 Color boundary image

圖7 區(qū)域編碼圖像Fig.7 Regional coding image

交替使用高值原子區(qū)域和低值原子區(qū)域，并保證在高值原子區(qū)域的上方、下方、左邊、右邊相鄰的原子區(qū)域只能為低值原子區(qū)域，在低值原子區(qū)域的上方、下方、左邊、右邊相鄰的原子區(qū)域只能為高值原子區(qū)域，任意兩個高值原子區(qū)域的彩色原子邊界都不能具有相同的顏色塊組成結(jié)構(gòu)，同時(shí)任意兩個低值原子區(qū)域的彩色原子邊界都不能具有相同的顏色塊組成結(jié)構(gòu)，即可生成彩色結(jié)構(gòu)光編碼圖像，如圖8所示。從圖中可以看出，該彩色編碼圖像由多個原子區(qū)域組成而形成整體的拓?fù)浣Y(jié)構(gòu)，每個原子區(qū)域中的編碼又形成了局部拓?fù)浣Y(jié)構(gòu)，在每個原子區(qū)域中的微區(qū)域則形成了局部結(jié)構(gòu)化的編碼，因此該編碼圖像具有較強(qiáng)的抗干擾性。

圖8 編碼圖像Fig.8 Coded image

3 彩色結(jié)構(gòu)光編解碼方法

在對攝像機(jī)拍攝的編碼圖像（解碼圖像）進(jìn)行解碼時(shí)，首先需要確定該圖像中的特征點(diǎn)的位置，然后，確定每一個特征點(diǎn)的顏色，最后根據(jù)彩色邊框編碼的拓?fù)浣Y(jié)構(gòu)和顏色邊框內(nèi)部的結(jié)構(gòu)確定解碼圖像中的每一個特征點(diǎn)與原始編碼圖像中的特征點(diǎn)的對應(yīng)關(guān)系，以實(shí)現(xiàn)準(zhǔn)確的解碼。

3.1 基于波形分析的特征點(diǎn)提取方法

由于編碼特征點(diǎn)信息存在顏色的變化、特征區(qū)域變形等干擾，因此，在提取特征點(diǎn)信息的過程中，很難準(zhǔn)確地提取出各個編碼特征點(diǎn)在編碼圖像中的位置，造成無法準(zhǔn)確地對特征點(diǎn)進(jìn)行解碼，從而不能正確地對三維實(shí)體進(jìn)行重構(gòu)操作。為了更加準(zhǔn)確地提取出編碼圖像中的特征點(diǎn)信息并對特征點(diǎn)進(jìn)行解碼，需要對攝像機(jī)攝取的編碼圖像進(jìn)行灰度化，并使用一種基于波形分析的特征點(diǎn)提取方法，確定特征點(diǎn)的位置。

該方法對灰度化的解碼圖像的每一行像素點(diǎn)的灰度值進(jìn)行分析，根據(jù)解碼圖像的編碼規(guī)則可知，每一個特征點(diǎn)都將處于解碼圖像灰度圖波形中的波峰或波谷區(qū)域之內(nèi)，當(dāng)確定出波形圖的波峰和波谷區(qū)域后，就可以對每一個特征區(qū)域進(jìn)行特征點(diǎn)的精確位置求取，即求取該特征區(qū)域的質(zhì)心位置。

在實(shí)際的灰度值波形圖中，波形圖一般都是不平滑的，在波峰和波谷區(qū)域之內(nèi)會出現(xiàn)波形的抖動而影響區(qū)域的確定。為了消除波形中抖動因素對圖像分割的影響，需要在波形分析時(shí)對一些次要的小峰進(jìn)行摒棄處理而引入濾波操作，即使用濾波方法對波形進(jìn)行平滑化處理。使用該方法可以更加有效地得到編碼特征點(diǎn)在解碼圖像中的位置。

3.2 基于聚類分析的顏色分類方法

在利用波形分析方法對解碼圖像進(jìn)行處理后，由于相鄰像素點(diǎn)對特征點(diǎn)的干擾，因此，很難準(zhǔn)確地確定特征點(diǎn)的顏色。文獻(xiàn)［15］提出了基于聚類分析的顏色分類方法，該方法可以有效地消除干擾，提高確定特征點(diǎn)顏色的精確性，但不可以準(zhǔn)確識別黑色且識別黑色的時(shí)間長，本文對該方法進(jìn)行了改進(jìn)，改進(jìn)后的方法可以識別任何顏色，并減少了處理時(shí)間。

為了確定特征點(diǎn)的顏色，分別用a1、a2、…、a8表示彩色結(jié)構(gòu)光解碼圖像中像素點(diǎn)的8種顏色，其中顏色ax的RGB值分別用rx、gx、bx表示，x＝1，2，…，8；設(shè)Mq表示屬于第q種顏色的像素點(diǎn)的個數(shù)，并且初始值都為零。基于聚類分析的顏色分類方法的實(shí)施步驟如下：

（1）設(shè)Nq＝0，q＝1，2，…，8。計(jì)算每一個編碼特征點(diǎn)在RGB顏色空間中的坐標(biāo)分別得到8種顏色在RGB顏色空間中的坐標(biāo)間的距離：

式中：D（i，j，s）表示特征點(diǎn)在RGB顏色空間中的坐標(biāo) （R（i，j），G（i，j），B（i，j））到點(diǎn)（rs，gs，bs）的距離。

（2）找到每一個編碼特征點(diǎn)的8個距離中最小的一個，如果 D（i，j，q）是｛（D（i，j，1），D（i，j，2），…，D（i，j，8）｝中的最小值且q值最小，說明該特征點(diǎn)的顏色屬于q類，Nq＝Nq＋1；并根據(jù)式（7）記錄特征點(diǎn)的RGB值。

如果對于每一個q＝1，2，…，8都有Nq＝Mq，那么基于聚類分析的顏色分類方法結(jié)束，否則繼續(xù)執(zhí)行步驟（3）。

（3）根據(jù)式（8）修改a1，a2，…，ak的 RGB值，然后用每一個Mq記錄對應(yīng)的Nq，轉(zhuǎn)而執(zhí)行步驟（1）。

3.3 編碼特征點(diǎn)的匹配方法

編碼特征點(diǎn)的匹配對特征點(diǎn)三維信息的求取至關(guān)重要，每一對特征點(diǎn)之間必須具有唯一性。通過波形分析得到特征點(diǎn)的位置后，根據(jù)該特征點(diǎn)在解碼圖像中的位置攜帶的拓?fù)浣Y(jié)構(gòu)信息和結(jié)構(gòu)化信息便可以確定特征點(diǎn)的編碼信息，從而能夠確定出與原始編碼圖像特征點(diǎn)之間的對應(yīng)關(guān)系。在對特征點(diǎn)進(jìn)行解碼時(shí)，因?yàn)槊總€特征點(diǎn)的顏色序列在編碼圖像中的組合方式是唯一的，所以可以通過顏色序列之間的拓?fù)潢P(guān)系和顏色序列內(nèi)部的拓?fù)潢P(guān)系來確定出該區(qū)域的結(jié)構(gòu)化模式。對特征點(diǎn)進(jìn)行解碼時(shí)，首先需要根據(jù)5×5子區(qū)域的彩色邊框編碼來唯一確定出該區(qū)域的匹配關(guān)系，如圖9所示。如果特征點(diǎn)是彩色編碼，則可根據(jù)邊框的拓?fù)浣Y(jié)構(gòu)確定特征點(diǎn)的編碼；如果特征點(diǎn)是黑色或白色，則需要對區(qū)域內(nèi)部的3×3子區(qū)域進(jìn)行結(jié)構(gòu)化解碼來確定特征點(diǎn)在特征區(qū)域內(nèi)的位置，再結(jié)合5×5子區(qū)域的解碼便可以確定出特征點(diǎn)的編碼，如圖10所示。

圖9 子區(qū)域5×5的解碼Fig.9 Decoding of 5×5sub－region

圖10 子區(qū)域3×3的解碼Fig.10 Decoding of 3×3sub－region

4 特征點(diǎn)三維信息求取

對特征點(diǎn)的三維信息進(jìn)行求取，需要根據(jù)已知的特征點(diǎn)計(jì)算出精度更高的三維信息，使得重構(gòu)出的三維實(shí)體更符合實(shí)際。由于在使用三角法對特征點(diǎn)的三維信息進(jìn)行求解的過程中，存在如攝像機(jī)與投影機(jī)夾角等參數(shù)的計(jì)算誤差，因此本文中使用無約束的求解方法來提高特征點(diǎn)的求解精度。當(dāng)?shù)玫浇獯a圖像中每個特征點(diǎn)的編碼值與調(diào)制前編碼圖像中特征點(diǎn)的編碼值之間的匹配關(guān)系后，便可以使用基于單應(yīng)性變換矩陣（MC和MP）的最大似然法計(jì)算待測三維實(shí)體表面特征點(diǎn)的三維信息，求取特征點(diǎn)三維信息的解算模型如圖11所示。

圖11 三維信息解算模型Fig.11 Three－dimensional information solution model

在計(jì)算得到攝像機(jī)單應(yīng)性變換矩陣MC和投影機(jī)單應(yīng)性變換矩陣MP的基礎(chǔ)上，通過聯(lián)立式（4）（5），便可以得到一個含有4個表達(dá)式的非齊次方程組：

使用最大似然法對公式（9）進(jìn)行求解，便可得到特征點(diǎn)在世界坐標(biāo)系下的坐標(biāo)值（XWp，YWp，ZWp）。

5 三維重構(gòu)實(shí)驗(yàn)

5.1 三維重構(gòu)系統(tǒng)組成

基于彩色結(jié)構(gòu)光的三維重構(gòu)系統(tǒng)包括硬件部分和軟件部分。系統(tǒng)的硬件部分如圖12所示，使用的計(jì)算機(jī)型號為DELL DIMENSION 5150。所用的投影機(jī)為可以滿足車系統(tǒng)要求的家用投影機(jī)，型號為EPSON EB－W6。攝像機(jī)是三維重構(gòu)系統(tǒng)中的一個關(guān)鍵組成部分，其工作主要是攝取編碼圖像，本系統(tǒng)使用的是Canon EOS。系統(tǒng)的軟件部分包括攝像機(jī)標(biāo)定模塊、投影機(jī)標(biāo)定模塊、編碼結(jié)構(gòu)光圖像采集模塊、數(shù)字圖像處理模塊、結(jié)構(gòu)光編碼圖像的解碼模塊、結(jié)構(gòu)光特征點(diǎn)匹配模塊、圖像與數(shù)據(jù)的并行處理和計(jì)算模塊以及點(diǎn)云數(shù)據(jù)三維信息計(jì)算模塊。

圖12 實(shí)驗(yàn)系統(tǒng)Fig.12 Experimental system

5.2 實(shí)驗(yàn)結(jié)果與分析

在系統(tǒng)參數(shù)標(biāo)定的過程中，標(biāo)定攝像機(jī)拍攝9幅標(biāo)定圖像，標(biāo)定投影機(jī)也拍攝9幅標(biāo)定圖像，包括世界坐標(biāo)系XY平面上的一張，向前、后、左、右移動的各兩張，標(biāo)定圖像每次移動的距離為20 mm。根據(jù)式（4）（5）進(jìn)行求解，可以計(jì)算出系統(tǒng)的參數(shù)值，如表1和表2所示。

表1 單應(yīng)性變換矩陣元素MCij值Table 1 MCijvalue of homography transformation matrix elements

表2 單應(yīng)性變換矩陣元素MPij值Table 2 MPijvalue of homography transformation matrix elements

圖13 調(diào)制的編碼圖像Fig.13 Modulation of coded image

對攝像機(jī)攝取的編碼圖像進(jìn)行解碼是關(guān)系到整個計(jì)算過程是否正確的關(guān)鍵，攝取的編碼圖像如圖13所示。在對頭像編碼圖像中的8741個編碼特征點(diǎn)進(jìn)行解碼時(shí)，僅有13個編碼特征點(diǎn)的解碼發(fā)生誤匹配，解碼正確率超過99.5%，進(jìn)一步證明了本文編解碼方法的有效性。通過系統(tǒng)參數(shù)的標(biāo)定和特征點(diǎn)編碼的匹配后，根據(jù)計(jì)算得到的特征點(diǎn)的世界坐標(biāo)對其進(jìn)行三維重構(gòu)就可以看到重構(gòu)后的三維實(shí)體，如圖14所示。三維點(diǎn)云的計(jì)算與重構(gòu)過程平均花費(fèi)不到0.1s，平均每秒鐘可以達(dá)到15幀，滿足動態(tài)重構(gòu)的需要；對特定特征點(diǎn)進(jìn)行尺寸重構(gòu)后，最小誤差為0.4mm，最大誤差為4.7mm，平均誤差為1.3mm。從圖14中可以看出，三維重構(gòu)后的實(shí)體模型達(dá)到了預(yù)期的三維重構(gòu)效果。

圖14 三維實(shí)體重構(gòu)Fig.14 Reconstruction of three－dimensional objects

6 結(jié)束語

通過實(shí)驗(yàn)可以看出，本文中的三維重構(gòu)方法不僅能夠很好地重構(gòu)出原始三維實(shí)體，而且對三維實(shí)體表面特征數(shù)據(jù)的提取也非常完整。使用文中的編碼方法，只需一幅編碼圖像便可以得到三維實(shí)體表面的全部信息，并且在一定程度上提高了編碼圖像的抗干擾能力，基本實(shí)現(xiàn)了編碼圖像的實(shí)用性和通用性設(shè)計(jì)。在對編碼圖像進(jìn)行解碼時(shí)，使用了波形分析的方法對特征信息進(jìn)行細(xì)化處理，使用基于聚類分析的顏色分類方法確定特征點(diǎn)的顏色，有效地消除了顏色干擾，利用編碼圖像中存在的拓?fù)浣Y(jié)構(gòu)，提高了特征點(diǎn)的解碼的準(zhǔn)確性。

［1］鄭作勇，姚莉，姚婷婷，等.一種簡易的線結(jié)構(gòu)光投影三維形狀獲取方法［J］.軟件學(xué)報(bào)，2006（17）：176－183.Zheng Zuo－yong.Yao Li，Yao Ting－ting，et al.A simple method for shape recovery based on linear structured light［J］.Journal of Software，2006（17）：176－183.

［2］陶立，孫長庫，何麗，等.基于結(jié)構(gòu)光掃描的彩色三維信息測量技術(shù)［J］.光電子·激光，2006，17（1）：111－114.Tao Li，Sun Chang－ku，He Li，et al.A color 3－D acquisition method based on structured－light scanning［J］.Journal of Optoelectronics·laser，2006，17（1）：111－114.

［3］Jordi Pages，Joaquim Salvi，Carles Matabosch.Implementation of a robust coded structured light technique four dynamic 3Dmeasurements［C］∥Proc of ICIP 2003，Barcelona，Spain，2003.

［4］廖素引，金施群，吳先良，等.單幅圖像的三維重構(gòu)視覺系統(tǒng)特征點(diǎn)提取算法的研究與實(shí)現(xiàn)［J］.宇航計(jì)測技術(shù)，2009，29（6）：12－15.Liao Su－yin，Jin Shi－qun，Wu Xian－liang，et al.Research and realization of the feature point extraction algorithm in 3Dreconstruction system using single image［J］.Journal of Astronautic Metrology and Measurement，2009，29 （6）：12－15.

［5］黃俊春，章炯民.基于彩色結(jié)構(gòu)光的實(shí)時(shí)三維重建［J］.計(jì)算機(jī)應(yīng)用與軟件，2009，26（12）：235－237.Huang Jun－chun，Zhang Jiong－min.Real－time 3D reconstruction based on colour－coded structured light［J］.Computer Applications and Software，2009，26（12）：235－237.

［6］Fechteler P，Eisert P，Rurainsky J.Fast and high resolution 3Dface scanning［C］∥Proc of ICIP 2007，San Antonio，Texas，USA，2007.

［7］Fechteler P，Eisert P.Adaptive color classification for structured light systems［C］∥Proc of CVPR 2008，Anchorage，AK，USA，2008.

［8］Koninckx Thomas P，van Gool Luc.Real－time range acquisition by adaptive structured Light［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2006，28（3）：432－445.

［9］張超，楊華民，韓成，等.基于格雷碼結(jié)構(gòu)光的編解碼研究［J］.長春理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，32（4）：365－368.Zhang Chao，Yang Hua－min，Han Cheng，et al.Study of coding and decoding based on gray－code structure light［J］.Journal of Changchun University of Science and Technology（Natural Science Edition），2009，32（4）：365－368.

［10］于曉洋，吳海濱.基于條紋邊緣解碼的結(jié)構(gòu)光三維測量技術(shù)［J］.電子器件，2008，31（2）：389－392.Yu Xiao－yang，Wu Hai－bin.3Dmeasurement tech－nology by structured light based on fringe edge decoding［J］.Chinese Journal of Electron Devices，2008，31 （2）：389－392.

［11］Brink Willie，Robinson Alan，Rodrigues Marcos.Indexing uncoded stripe patterns in structured light systems by maximum spanning trees［C］∥British Machine Vision Conference，Leeds，UK，2008.

［12］孔令富，陳淑平，郭福田.一種基于編碼結(jié)構(gòu)光的三維重建匹配算法［C］∥計(jì)算機(jī)技術(shù)與應(yīng)用進(jìn)展2007全國第18屆計(jì)算機(jī)技術(shù)與應(yīng)用（CACIS）學(xué)術(shù)會議論文集，浙江，中國，2007.

［13］Zhang Li，Curless Brian，Seitz Steven M.Rapid shape acquisition using color structured light and multi－pass dynamic programming［C］∥Proc of the 1st Int＇l Symp on 3DData Processing Visualization and Transmission（3DPVT）.Padova：IEEE Com－puter Society Press，2002.

［14］韋爭亮，鐘約先，袁朝龍，等.單攝像機(jī)單投影儀三維測量系統(tǒng)標(biāo)定技術(shù)［J］.清華大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，49（2）：202－205.Wei Zheng－liang，Zhong Yue－xian，Yuan Chaolong，et al.Calibration of a 3Dmeasurement system having one camera and one projector［J］.Journal of Tsinghua University（Science and Technology），2009，49（2）：202－205.

［15］范靜濤，韓成，張超，等.一種新的De Bruijn彩色結(jié)構(gòu)光解碼技術(shù)研究［J］.電子學(xué)報(bào)，2012，40（3）：483－488.Fan Jing－tao，Han Cheng，Zhang Chao，et al.Study of a new decoding technology for De Bruijn structured light［J］.ACTA Electronica Sinica，2012，40（3）：483－488.