基于頻率合成的相位解包裹方法

摘 要：相位展開是投射正弦條紋的結(jié)構(gòu)光三維掃描技術(shù)的關(guān)鍵環(huán)節(jié)。本文提出了一種絕對(duì)相位展開方法，通過四步相移法，投射三組不同頻率的正弦光柵條紋，再以低頻的展開相位為基礎(chǔ)，展開中頻光柵和高頻光柵，相位展開不依賴展開路徑、不依賴相鄰像素、不需要任何參數(shù)的絕對(duì)相位。

關(guān)鍵詞：頻率合成;相位展開;相移法;三維掃描技術(shù)

中圖分類號(hào)：TP391.41文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1003-5168（2019）29-0084-04

Phase Unwrapping Method Based on Frequency Synthesis

ZHOU Bo

（Sanya College，Sanya Hainan 572022）

Abstract： Phase unwrapping is a key part of the structured light 3D scanning technique for projecting sinusoidal fringes. This paper proposed an absolute phase expansion method， three sets of sinusoidal grating strips of different frequencies were projected by the four-step phase shift method， and the intermediate frequency grating and the high frequency grating were developed based on the unwrapped phase of the low frequency， phase unwrapping did not depend on the unwrapped path， did not depend on adjacent pixels， and did not require any parameters of the absolute phase.

Keywords： frequency synthesis;phase unwrapping;phase shift method;three-dimensional scanning technology

投射正弦光柵的結(jié)構(gòu)光三維掃描技術(shù)是一種重要的光學(xué)三維測(cè)量方法，通過相移法可以獲得高精度的相位，但其相位為[0，2π]的包裹相位，需要將包裹相位展開才能得到全場(chǎng)相位[1]。目前，相位解包裹方法主要有空間解包裹和時(shí)間解包裹兩大類?？臻g相位解包裹方法從一幅包裹相位圖出發(fā)，通過加減[2nπ]，使相鄰像點(diǎn)的相位差在[-π，π]之間實(shí)現(xiàn)相位展開，空間解包裹方法只適用物體輪廓連續(xù)、無遮擋且圖像條紋無交疊的情況。時(shí)間相位解包方法通過投射不同頻率的光柵條紋，獲取多個(gè)相位圖進(jìn)行相位展開。

自Saldner和Huntley于1993年提出時(shí)間相位展開方法開始，人們就如何減少所需頻率數(shù)目、提高算法效率開展研究，探索出多種算法[2-4]。典型算法包括線性增長(zhǎng)、指數(shù)增長(zhǎng)、負(fù)指數(shù)增長(zhǎng)等時(shí)間相位展開方法，其中，負(fù)指數(shù)增長(zhǎng)算法具有很高的條紋級(jí)次定位精度，相位展開精度高，但是所需的數(shù)據(jù)采集量較大，數(shù)據(jù)處理時(shí)間長(zhǎng)。其改進(jìn)方法包括三頻法、外差法等，精度與負(fù)指數(shù)增長(zhǎng)法相比較低，但縮短了測(cè)量和處理時(shí)間。三頻外差法將三頻展開法與外差法結(jié)合起來，相位精度也達(dá)不到負(fù)指數(shù)增長(zhǎng)法的精度，尤其是存在條紋交疊時(shí)，這兩種方法都無法給出令人滿意的結(jié)果。

本文提出的頻率合成相位解包裹方法也是一種時(shí)間相位解包裹方法，采用四步相移法，向物體表面投射三組頻率光柵條紋，要求其光柵頻率滿足[f1+f3-2f2=1]。基于光拍原理，該方法既具三頻外差法的數(shù)據(jù)處理速度，同時(shí)也能夠達(dá)到負(fù)指數(shù)增長(zhǎng)法的精度。

1 頻率合成相位展開理論

1.1 相移法

相移法有三步相移、四步相移、五步相移以及等間隔相移等方法，以四步相移為例，需向物體表面投射四個(gè)初始相位相差[π2]的同頻率光柵，攝像機(jī)像平面的光強(qiáng)分布[I]可表示為：

[I1x=I0+Acos?x+0·π2]? ? ? ? ? ? ? ? ?（1）

[I2x=I0+Acos?x+π2]? ? ? ? ? ? ? ? ? ? ? ? （2）

[I3x=I0+Acos?x+2π2]? ? ? ? ? ? ? ? ? ?（3）

[I4x=I0+Acos?x+3π2]? ? ? ? ? ? ? ? ?（4）

式中，[I1]、[I2]、[I3]、[I4]為攝像機(jī)獲取的4個(gè)光強(qiáng)分布;[I0]為背景光強(qiáng);[A]為條紋強(qiáng)度幅值;[?x]為相位。

聯(lián)立上述方程，求得相位函數(shù)為：

[?x=arctanI4x-I2xI1x-I3x]? ? ? ? ? ? ? ? ? ? ?（5）

1.2 頻率合成原理

兩個(gè)振動(dòng)方向相同，振幅相等，頻率相差很小的單色光波疊加可以形成拍，把正弦光柵作為光波。頻率分別為[ω1]和[ω2]的兩個(gè)光波，即[E1=A1cosω1x+φ1]、[E2=A2cosω2x+φ2]，其合成振動(dòng)[E=E1+E2]，假設(shè)[A1=A2=A]，[φ1=φ2=φ]，則

[E=2Acosω1-ω22x?cosω1+ω22x+φ]? ? ? ? （6）

當(dāng)[ω1?ω2]時(shí)，且[ω1+ω2?ω1+ω2]，因此合成振動(dòng)可以看作是振幅緩變的簡(jiǎn)諧振動(dòng)。

合成光波的光強(qiáng)為：

[Ix=2Acosω1-ω22x2=4A2cos2ω1-ω22x=4A2+4A2cosω1x-ω2x]? ? ? （7）

合成光波的強(qiáng)度在[0?4A2]之間變化，這種光強(qiáng)時(shí)大時(shí)小的現(xiàn)象稱為拍。拍的頻率為兩個(gè)疊加光波頻率之差[ω1-ω2]。

1.3 相位展開步驟

光學(xué)三維測(cè)量所用正弦光柵頻率越高，噪聲越小，相位越準(zhǔn)確。為使三維測(cè)量具有較高的精度，三個(gè)正弦光柵都需要具有較高的頻率，光柵頻率為投射裝置投射的條紋數(shù)，從高到低依次為[f1]、[f2]、[f3]，合成光柵的頻率分別為：

[f12=f1-f2]? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （8）

[f23=f2-f3]? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? （9）

[f13=f1-f3]? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（10）

[f123=f12-f23=f1-f2-f2+f3=f1-2f2+f3]? ? ? （11）

一是從三組不同頻率的條紋圖像出發(fā)，計(jì)算其對(duì)應(yīng)的正余弦圖像，再計(jì)算光拍合成頻率對(duì)應(yīng)的正余弦圖像，各項(xiàng)均可由條紋圖像計(jì)算生成。

[sin?1x=I14x-I12xcos?1x=I13x-I11xsin?2x=I24x-I22xcos?2x=I23x-I21xsin?3x=I34x-I32xcos?3x=I33x-I31xsin?12x=sin?1xcos?2x-cos?1xsin?2xcos?12x=cos?1xcos?2x+cos?1xsin?2xsin?23x=sin?2xcos?3x-cos?2xsin?3xcos?23x=cos?2xcos?3x+cos?2xsin?3xsin?123x=sin?12xcos?23x-cos?12xsin?23xcos?123x=cos?12xcos?23x+cos?12xsin?23x]? ? ? （12）

二是將正余弦圖像代入反正切函數(shù)，計(jì)算包裹相位。

[?1x=arctan2sin?1x，cos?1x?2x=arctan2sin?2x，cos?2x?3x=arctan2sin?3x，cos?3x?12x=arctan2sin?1x-?2x，cos?1x-?2x?23x=arctan2sin?2x-?3x，cos?2x-?3x?123x=arctan2sin?12x-?23x，cos?12x-?23x]? ? ? （13）

三是因?yàn)閇f123=f1-2f2+f3]，設(shè)計(jì)一組頻率，使[f123=1]，頻率為1，意味著在視場(chǎng)范圍內(nèi)僅有一個(gè)正弦條紋，故其包裹相位[?123]不存在相位包裹問題，即其本身就是展開相位[φ123]，以[φ123]為依據(jù)，逐次展開其他包裹相位，相位展開公式如下。

[φ12x=?12x+2π·floorf12f123?φ123x?12xφ23x=?23x+2π·floorf23f123?φ123x?23xφ13x=?13x+2π·floorf13f123?φ123x?13xφ1x=?1x+2π·floorf1f12?φ12x?1xφ2x=?2x+2π·floorf2f12?φ12x?2xφ3x=?3x+2π·floorf3f23?φ23x?3x]? ? ?（14）

由上述過程進(jìn)行的相位展開，每個(gè)像素點(diǎn)獨(dú)立的相位可獨(dú)立進(jìn)行，適用于GPU進(jìn)行多線程并行處理。

2 試驗(yàn)

攝像機(jī)采用了Basler102f黑白攝像機(jī)，像素為1 390×1 040，投影機(jī)采用Lightcrafter3400，被測(cè)物體為汽車后視鏡外殼，如圖1所示。

由式（12）和式（13）計(jì)算得到包裹相位，如圖3至圖8所示。

在光照相對(duì)較差的地方，條紋對(duì)比度較小，則該區(qū)域的相位可靠性較差。

3 結(jié)論

基于光拍合成頻率的相位解包裹算法，以光拍理論為基礎(chǔ)，合成條紋頻率為[f1-f2]、[f2-f3]、[f1-f3]和[f1-2f2+f3]的相位，使相位由低頻[f1-2f2+f3]到中頻[f1-f2]、[f2-f3]、[f1-f3]再到高頻[f1]、[f2]、[f3]逐步展開，每個(gè)像素位置上相位的計(jì)算僅僅與該點(diǎn)在低、中、高頻相位圖上對(duì)應(yīng)點(diǎn)的相位值相關(guān)，與相鄰點(diǎn)無關(guān)，不存在誤差向相鄰點(diǎn)傳遞的問題，不存在條紋交疊的影響，不需要進(jìn)行誤差校正。相位精度僅僅與條紋圖像采集過程中的測(cè)量環(huán)境相關(guān)，比如在光照相對(duì)較差的地方，在物體遮擋和不連續(xù)的位置，條紋對(duì)比度較小，則該區(qū)域的相位可靠性較差。采用調(diào)制度技術(shù)進(jìn)行相位可靠性檢測(cè)，可以剔除這部分區(qū)域，為光學(xué)三維測(cè)量提供可靠、高精度的相位。

參考文獻(xiàn)：

[1]牛小兵.光柵投射三維輪廓測(cè)量技術(shù)的研究[D].天津：天津大學(xué)，2003.

[2]岳慧敏，蘇顯渝.時(shí)間位相展開方法研究進(jìn)展[J].激光雜志，2004（3）：9-12.

[3]周波.投射多頻光柵的物體表面三維輪廓的視覺測(cè)量方法：中國，CN101105393[P].2018-01-16.

[4]梁晉，肖振中，臧順來，等.外差式多頻相移技術(shù)的三維光學(xué)點(diǎn)云測(cè)量研究[J].鍛壓技術(shù)，2008（1）：143-147.