光柵投影三維測(cè)量解相位方法的對(duì)比研究

徐昆鵬，王建華

（1.徐州工業(yè)職業(yè)技術(shù)學(xué)院 機(jī)電工程學(xué)院，江蘇徐州，221140；2.青島理工大學(xué) 自動(dòng)化工程學(xué)院，山東青島，266520）

0 引言

基于光柵投影的光學(xué)三維面形測(cè)量具有非接觸、精度高、速度快、低成本、全場(chǎng)測(cè)量和易于實(shí)現(xiàn)等特點(diǎn)，在逆向工程、工業(yè)檢測(cè)、質(zhì)量控制、虛擬現(xiàn)實(shí)、醫(yī)學(xué)診斷、文物修復(fù)和物體識(shí)別等領(lǐng)域應(yīng)用廣泛[1]。相移輪廓術(shù)(PSM)和傅里葉變換輪廓術(shù)(FTP)是兩種傳統(tǒng)方法。相移法[2]計(jì)算簡(jiǎn)單，測(cè)量精度高，但對(duì)噪聲比較敏感，噪聲會(huì)引起解相位結(jié)果存在大量虛假2π跳變，增加了相位展開(kāi)的難度，且至少需要三幅以上的條紋圖像，適用于靜態(tài)測(cè)量；傅里葉變換法[3]是M.Takeda等人在1983年提出來(lái)的，它只需要一幅條紋圖案，即可恢復(fù)物體三維形貌，可用于動(dòng)態(tài)物體的三維形貌測(cè)量；但由于傅里葉變換不具有局部分析的能力，頻譜混疊等問(wèn)題，影響了其測(cè)量精度。小波變換具有多分辨率分析特點(diǎn)，抑噪能力強(qiáng)，可分為一維小波變換輪廓術(shù)(1D WTP)和二維小波變換輪廓術(shù)(2D WTP)[4-6]。小波脊提取是小波變換輪廓術(shù)的關(guān)鍵步驟，最大模法是常見(jiàn)的小波脊提取方法。

本文分析了相移法、傅里葉變換法和一維、二維小波變換法解相位的原理，對(duì)比了幾種方法的解相位結(jié)果。針對(duì)條紋圖存在較強(qiáng)噪聲干擾，小波變換最大模提取算法會(huì)產(chǎn)生較大誤差，提出了基于價(jià)值函數(shù)的二維小波變換小波脊提取算法。首先，提取二維小波變換系數(shù)模的最大值點(diǎn)和部分局部極值點(diǎn)，共同作為小波脊候選點(diǎn)；其次，根據(jù)噪聲產(chǎn)生的較大系數(shù)模可以使用小波脊連續(xù)性進(jìn)行去除這一原理，在二維小波變換系數(shù)的模上引入尺度因子的梯度，建立價(jià)值函數(shù)，評(píng)估所有候選點(diǎn)的價(jià)值；最后，確定每行條紋的最優(yōu)小波脊。仿真和實(shí)驗(yàn)證明，該方法的抑噪能力明顯優(yōu)于其他幾種解相位方法。

1 原理

1.1 相移輪廓術(shù)

相移法[7]的投影條紋一般為余弦條紋（或正弦），得到條紋的灰度值用下式（1）表示：

式中，a(x,y)為背景光強(qiáng)值，b(x,y)為調(diào)制強(qiáng)度，f0為載波頻率，φ(x,y)為被測(cè)物體高度調(diào)制相位。

沿著垂直于光柵條紋方向在一個(gè)周期內(nèi)移動(dòng)n次，n∈ (3,5)。當(dāng)n= 3時(shí)稱(chēng)為三步相移法，當(dāng)n= 4時(shí)稱(chēng)為四步相移法，當(dāng)n= 5時(shí)稱(chēng)為五步相移法。以四步相移法為例，四幅帶有 /2π相移的灰度值函數(shù)分別用下式（2）～（5）表示：

得到：

求取反正切函數(shù)，得到：

式（6）反正切運(yùn)算得到的相位值 (,)xyψ包裹在(-,]ππ之間，稱(chēng)為相位主值（或包裹相位），需要相位展開(kāi)得到實(shí)際相位。

1.2 傅里葉變換輪廓術(shù)

傅里葉變換解相位的基本原理是在頻域中把代表背景光的零頻分量和代表噪聲的高頻分量通過(guò)合適的濾波器濾除掉，剩下的是包含有用相位信息的基頻分量，再對(duì)基頻分量進(jìn)行逆傅里葉變換，通過(guò)正切運(yùn)算即可解出被測(cè)物體高度調(diào)制相位主值。

根據(jù)歐拉公式，式（1）可以（7）表示為：

式中：

沿x方向?qū)κ剑?）采用一維傅里葉變換，得到：

上式（8）由-1級(jí)C*(f+f0)、0級(jí)A(f)、1級(jí)C(f-f0)三種頻譜信號(hào)組成，如圖1所示。

圖1 一維傅里葉變換的基頻分量提取

其中A(f)其代表零頻背景信號(hào)，而C(f-f0)、C*(f+f0)包含了所需的相位信號(hào)。這里選取 0(-)Cff，逆傅里葉變換得到：

其實(shí)部和虛部分別為：

被測(cè)物體高度調(diào)制相位主值可以寫(xiě)成：

同相移法類(lèi)似，傅里葉變換法得到的相位主值也是截?cái)嗟模窗?-,]ππ之間，需要相位展開(kāi)得到實(shí)際相位。

1.3 小波變換輪廓術(shù)

1.3.1 一維小波變換法

圖2(a)為模擬被測(cè)物體1Peaks× ，圖像大小為512pixels×512pixels，被測(cè)物體高度調(diào)制后的變形條紋如圖2(b)所示，選取變形條紋第256行像素的灰度值組成一維數(shù)組g(256,:)(x)。

圖2 被測(cè)物體和變形條紋

對(duì)一維數(shù)組g(256,:)(x)一維小波變換：

式中： (,)Wbs為一維小波變換系數(shù)，b為平移因子，s為伸縮因子（或稱(chēng)尺度因子），母小波 ()xψ選擇一維Morlet小波，表達(dá)式為：

一維小波變換系數(shù)是二維的，其模和相位分布分別用公式（13）、（14）表示：

式中：Modulus是一維小波變換系數(shù)的模，如圖3(a)所示，橫軸為平移因子b，取值范圍為1～512pixles，豎軸為系數(shù)模，白色表示系數(shù)模比較大，黑色表示系數(shù)模比較小；phase_map是一維小波變換系數(shù)的相位，如圖3(b)所示，橫軸為尺度因子s，取值范圍為1～64，豎軸為相位。

圖3 變形條紋的一維小波變換

小波變換系數(shù)模的最大值定義為小波變換的“脊”，“脊線”則代表了信號(hào)和Morlet小波最為相似的位置，提取出“脊線”位置處對(duì)應(yīng)的相位，即為第256行像素的相位主值，條紋所有行像素采用同樣方式得到相位主值，從而獲得整幅條紋圖案的相位主值分布[8-11]。

1.3.2 二維小波變換法

二維小波變換不但含有平移因子和尺度因子，而且含有旋轉(zhuǎn)因子，因此具有方向選擇性，其多分辨率分析能力強(qiáng)于一維小波變換。其原理是把小波不斷平移、縮放和旋轉(zhuǎn)，與原始信號(hào)進(jìn)行比較，計(jì)算得到小波變換系數(shù)矩陣。系數(shù)模代表小波與信號(hào)的相似程度，系數(shù)模值越高，表示小波與信號(hào)越相似。通過(guò)尋找二維小波系數(shù)的最大模，來(lái)恢復(fù)相位主值[12]。

在二維小波變換中，母小波采用二維Morlet復(fù)小波[13-14]，表示為：

式中：θ為小波相對(duì)于x旋轉(zhuǎn)角度，k為小波中心頻率。

變形條紋 (,)gxy二維小波變換表達(dá)式為：

式中：W(a,b,s,θ)為小波變換系數(shù)，(a,b)為兩個(gè)方向上的平移因子，s為伸縮因子（或稱(chēng)尺度因子），θ為旋轉(zhuǎn)角度，rθ為旋轉(zhuǎn)矩陣。

如果尺度因子s有六十四種不同的取值，旋轉(zhuǎn)角度因子s有四種不同的取值，那么將整幅變形條紋二維小波變換后得到一個(gè)四維數(shù)據(jù)為512×512×64×4，取出每個(gè)512×512矩陣第一個(gè)位置點(diǎn)的數(shù)值，組成64×4系數(shù)矩陣，它是第1行、第1列像素點(diǎn)在不同縮放尺度和旋轉(zhuǎn)角度下的二維小波變換系數(shù)矩陣，如圖4所示。

圖4 二維小波變換的四維數(shù)據(jù)

根據(jù)式（13）、（14）求出該像素點(diǎn)的模和相位，找到系數(shù)模最大值，進(jìn)一步得到最大值點(diǎn)的相位值，即為該像素點(diǎn)的相位主值。依次提取出512×512矩陣其他像素點(diǎn)的數(shù)據(jù)，提取最大系數(shù)模位置相對(duì)應(yīng)的相位，從而得到整幅條紋的相位主值。

1.4 幾種方法的仿真對(duì)比

計(jì)算機(jī)生成模擬被測(cè)物體為0.5Peaks，條紋周期為16Pixels，相位展開(kāi)采用Itoh提出的基本相位展開(kāi)方法[15-16]（逐行相位展開(kāi)方法），以均方根誤差（RMSE）作為評(píng)價(jià)標(biāo)準(zhǔn)，對(duì)比各種方法的面形恢復(fù)、誤差如表1所示。通過(guò)對(duì)比幾種方法的解相位結(jié)果可以看出，相移法測(cè)量精度最高，誤差最小。

表1 理想無(wú)噪聲干擾條紋圖，不同方法的三維面形恢復(fù)和誤差

2 基于評(píng)價(jià)函數(shù)二維小波變換脊提取算法

小波變換直接最大模小波脊提取方法的特點(diǎn)是算法簡(jiǎn)單，速度快。但是，當(dāng)條紋圖像中存在噪聲時(shí)，會(huì)出現(xiàn)多個(gè)局部極大值點(diǎn)，小波變換系數(shù)模不是最大值點(diǎn)，也是局部極大值點(diǎn)。如果仍然采用模最大值搜索“脊”，將找到錯(cuò)誤的“脊”，而導(dǎo)致恢復(fù)的三維形貌存在較大誤差。通過(guò)前面的分析可以看出，系數(shù)模的“脊”一般是連續(xù)的，而噪聲產(chǎn)生的局部極值點(diǎn)可以選用系數(shù)模“脊線”的連續(xù)性進(jìn)行去除。根據(jù)以上特性，提出了基于評(píng)價(jià)函數(shù)的二維小波脊提取算法。

2.1 評(píng)價(jià)函數(shù)

評(píng)價(jià)函數(shù)定義為：

式中：φ(b)為尺度因子，b為平移因子，S[φ(b),b]為小波變換系數(shù)((),)bbφ位置的模，代表了尺度因子的梯度，C0和C1為權(quán)值。

根據(jù)前面分析，系數(shù)模越大且路徑越光滑，則越能代表實(shí)際“脊線”。由式（17）可以看出，S[φ(b),b]越大，越小，價(jià)值函數(shù)值也越小，越能反映真實(shí)的“小波脊線”。

假設(shè)C0和C1取值為1，價(jià)值函數(shù)最小的路徑經(jīng)過(guò)點(diǎn)(p,j)，其中尺度因子p=φ(j)，j為平移因子，則后面一點(diǎn)(q,j+ 1)的價(jià)值函數(shù)可表示為：

2.2 算法步驟

Step1: 對(duì)條紋圖進(jìn)行二維小波變換，得到一個(gè)四維的系數(shù)矩陣。

Step2: 0θ= 時(shí)，取出尺度因子為1,2……64時(shí)每個(gè)矩陣的第一行系數(shù)，組成64512× 系數(shù)矩陣，它代表了 0θ= 時(shí)、第一行像素的系數(shù)矩陣。采用相同的方法，取出其他三個(gè)旋轉(zhuǎn)角度中每個(gè)矩陣的第一行系數(shù)，組成各自的第一行像素系數(shù)矩陣，如圖5所示。

圖5 基于價(jià)值函數(shù)的二維小波變換四維系數(shù)矩陣

Step3: 取θ= 0時(shí)，第一行像素的64 × 512系數(shù)矩陣的模和如圖6所示。在圖6中依次找出平移因子為b= 1,b= 2， ······b=512時(shí)，模的最大值和最大值90%的局部極值點(diǎn)，共同組成候選脊點(diǎn)，圖7為候選脊點(diǎn)的模，而候選脊點(diǎn)所對(duì)應(yīng)尺度因子如圖8所示。由于噪聲的影響，所有候選點(diǎn)都有可能是最終的“脊”。因此，根據(jù)式（18）計(jì)算所有候選點(diǎn)的價(jià)值；第一列的價(jià)值都為0，第2列價(jià)值的計(jì)算過(guò)程如圖9所示。

圖6 二維小波變換系數(shù)矩陣 0θ= 第一行像素的模

圖7 候選脊點(diǎn)的模

圖8 候選脊點(diǎn)的尺度因子

圖9 候選脊點(diǎn)的價(jià)值計(jì)算

Step4: 找到第512列像素點(diǎn)的所有候選點(diǎn)的最小價(jià)值，并記錄該候選點(diǎn)的位置（即尺度因子和系數(shù)模），然后找出該最小價(jià)值的求解路徑，從而確定第511像素點(diǎn)的小波脊位置，然后按照上述方法反向?qū)ふ业剿凶顑?yōu)路徑，即可找到第512～1列位置，即為 0θ= 時(shí)第一行像素的最優(yōu)小波脊線，如圖10所示。

圖10 最優(yōu)小波脊線

3 計(jì)算機(jī)模擬和實(shí)驗(yàn)

3.1 計(jì)算機(jī)模擬

采用計(jì)算機(jī)模擬驗(yàn)證所提方法的抑噪性能。模擬被測(cè)物體為0.5×Peaks，條紋圖大小為512pixle×512pixel，采用均方根誤差（RMSE）衡量形貌恢復(fù)精度。依次加入標(biāo)準(zhǔn)方差梯度0.1遞增的高斯白噪聲和椒鹽噪聲，對(duì)比二維小波變換的直接最大模的小波脊提取和基于價(jià)值函數(shù)的小波脊提取，如圖11和12所示。信噪比較低時(shí)，基于價(jià)值函數(shù)的小波脊提取算法的均方根誤差明顯小于直接最大模的小波脊提取方法。

圖11 不同椒鹽噪聲等級(jí)時(shí)的RMSE

圖12 不同高斯白噪聲等級(jí)時(shí)的RMSE

3.2 實(shí)驗(yàn)

測(cè)量系統(tǒng)主要由數(shù)字投影儀（三星SP-P410M）、CCD相機(jī)（大恒水星系列MER-050-560U3M）和個(gè)人計(jì)算機(jī)（Intel Core i5-4258U CPU, 4G RAM）組成，如圖13所示。

圖13 測(cè)量系統(tǒng)

被測(cè)物體為三維凸起曲面，條紋圖存在噪聲污染，采集到的條紋圖如圖14(a)所示，分別采用二維小波變換最大模的小波脊提取和基于價(jià)值函數(shù)的小波脊提取算法獲得包裹相位，然后用Itoh逐行逐列相位展開(kāi)方法得到實(shí)際相位分布，結(jié)果如圖14(b)和(c)所示。通過(guò)對(duì)比可以看出，相對(duì)于最大模提取算法，基于價(jià)值函數(shù)的二維小波脊提取形貌恢復(fù)精度更高，結(jié)果更好。

圖14 曲面實(shí)驗(yàn)對(duì)比

4 結(jié)論

光柵條紋投影輪廓術(shù)的解相位方法主要有相移法、傅里葉變換法和小波變換法。本文分析了幾種方法的原理，并對(duì)測(cè)量精度進(jìn)行了對(duì)比分析，研究發(fā)現(xiàn)相移法解相位精度最高，但所需條紋數(shù)量較多，適合靜態(tài)物體測(cè)量；傅里葉變換和小波變換只需一幅條紋即可獲得相位主值，但測(cè)量精度低于相移法。此外，當(dāng)條紋圖存在較強(qiáng)噪聲干擾，小波變換最大模提取算法會(huì)產(chǎn)生較大誤差，提出了一種基于價(jià)值函數(shù)的二維小波變換小波脊提取算法，相對(duì)于二維小波變換直接最大脊提取算法，抑噪能力更好。