應(yīng)用頻率積分相位解調(diào)測量徑向畸變

楊初平，劉建斌，譚穗妍，翁嘉文

（華南農(nóng)業(yè)大學(xué)物理系，廣州510642）

應(yīng)用頻率積分相位解調(diào)測量徑向畸變

楊初平，劉建斌，譚穗妍，翁嘉文

（華南農(nóng)業(yè)大學(xué)物理系，廣州510642）

為了測量光學(xué)成像系統(tǒng)的徑向畸變，采用載頻條紋模板，應(yīng)用瞬時頻率積分法提取因徑向畸變而產(chǎn)生的徑向調(diào)制相位；推導(dǎo)了條紋徑向調(diào)制相位與瞬時頻率的關(guān)系式，并導(dǎo)出徑向調(diào)制相位和徑向畸變位移關(guān)系；采用小波頻率估計提取畸變條紋徑向瞬時頻率，并對其進(jìn)行積分獲得畸變條紋的徑向調(diào)制相位；應(yīng)用徑向調(diào)制相位和立方卷積插值算法對畸變圖像進(jìn)行了校正，得出了詳細(xì)的理論分析和實驗結(jié)果。結(jié)果表明，上述方法是可行的。

測量與計量；畸變測量；載頻條紋；相位解調(diào)；徑向畸變

引 言

畸變普遍存在于光學(xué)成像系統(tǒng)中。徑向畸變測量就是獲得徑向畸變位移函數(shù)即畸變像點相對于無畸變像點的位移，校正畸變圖像。通過設(shè)定特征點、提取特征點的畸變位移的標(biāo)定模板［1-7］是目前測量徑向畸變的主要方法。如果把這些特征點的灰度按正弦函數(shù)分布排列，則測量因徑向畸變導(dǎo)致的特征點的位移將轉(zhuǎn)化成測量灰度條紋徑向調(diào)制相位。本文中以縱向正弦載頻灰度條紋作為模板，把徑向位置畸變轉(zhuǎn)化為徑向調(diào)制相位，采用小波頻率估計提取畸變條紋徑向瞬時頻率并進(jìn)行積分運算得到徑向調(diào)制相位。與目前提取條紋調(diào)制相位都需要進(jìn)行相位解包［8-14］不同，采用瞬時頻率積分相位解調(diào)方法提取變形條紋調(diào)制相位方法沒有產(chǎn)生包裹相位，因而無需經(jīng)相位解包過程，便可直接得到調(diào)制相位。然后換算為徑向位置畸變，實現(xiàn)對圖像的校正。

1 圖像畸變模型

畸變圖像建立如圖1所示的坐標(biāo)系。徑向畸變相對于圖像中心點O具有旋轉(zhuǎn)對稱性，在同一圓周上的點具有相等的位置畸變。設(shè)（x0，y0），（x，y），（xd，yd）分別表示圖像中心點、理想像點和對應(yīng)的畸

Fig.1 Diagram illustrating radial distortion of fringe patterns

式中，徑向畸變位移函數(shù)T（r）是描述畸變像點相對于無畸變像點的位移，表示為：

通過實驗得到Δr，代入（1）式并結(jié)合灰度插值得到校正圖像。

2 徑向畸變相位

2.1 徑向調(diào)制相位與徑向瞬時頻率

以條形正弦灰度模板進(jìn)行成像，采集的畸變條紋和無畸變條紋的關(guān)系如圖1所示，其中徑向OAB經(jīng)過圖像中心O點。采集的無畸變條紋像的灰度可表示如下：

式中，A1（x）是背景灰度，B1（x）是正弦灰度的幅值，相位分布是（2πf0x＋φ），f0是頻率，φ是初始相位。

一方面，由于畸變，采集的條形模板灰度不再是正弦分布，可表示為：

式中，相位分布是［2πf0x＋φ（x）＋φ］，不再是周期性分布，而是出現(xiàn)沿著徑向發(fā)生變化的調(diào)制相位函數(shù)φ（x）?；儣l紋像與無畸變條紋像比較，因徑向畸變產(chǎn)生的相位差即徑向調(diào)制相位函數(shù)為：

另一方面，徑向畸變表現(xiàn)為畸變像點相對于無畸變像點的位置移動，如圖1所示，徑向畸變使B點移到A點，產(chǎn)生的徑向畸變位移函數(shù)表示為Δr（x）。則徑向畸變相位函數(shù)與徑向畸變位移函數(shù)的關(guān)系為：

因此，只要得到徑向調(diào)制相位就可以計算徑向畸變位移分布。

畸變條形模板灰度分布不再是固定的頻率f0，而是隨著徑向發(fā)生變化?；儣l形徑向瞬時頻率定義為：

把f0移到等號左邊，兩邊積分得：φ（x）－φ（0）＝在畸變像的中心點，畸變?yōu)?，則瞬時頻率f（0）＝f0，而且φ（0）＝0。因此徑向調(diào)制相位函數(shù)為：

因此，若得到瞬時頻率分布f（x），則可以計算徑向各個畸變點的徑向調(diào)制相位。在以往的條紋調(diào)制相位解調(diào)過程中，由于采用反正切函數(shù)計算條紋每個點的調(diào)制相位，大于2π的調(diào)制相位將被折斷到（－π，π）范圍內(nèi)，為了獲得真實調(diào)制相位，需要對解調(diào)相位進(jìn)行解包處理。采用瞬時頻率積分不會出現(xiàn)相位折斷，因而無需相位解包過程就能夠直接計算真實的調(diào)制相位。

2.2 頻率估計

條紋上各點瞬時頻率的計算可以采用連續(xù)小波變換進(jìn)行估計。對1維函數(shù)f（x），其連續(xù)小波變換系數(shù)Wf（a，b）定義為［11-13］：

母函數(shù)ψ（x）采用Gabor小波函數(shù)［12］：

式中，ψa，b（x）為母函數(shù)經(jīng)過尺度因子a伸縮和平移因子b平移后得到的小波序列，ψa，b*（x）是ψa，b（x）的共軛函數(shù)。對實信號f（x）的各個b點采用不同a值進(jìn)行連續(xù)小波變換，得到一系列的小波系數(shù)Wf（a，b）。小波系數(shù)的幅值A(chǔ)（a，b）和相位φ（a，b）可以通過Wf（a，b）的實部Re［Wf（a，b）］和虛部Im［Wf（a，b）］分別表示為［12］：

對信號各個位置b，通過檢測不同a值的小波變換系數(shù)幅值A(chǔ)（a，b），其中最大值所對應(yīng)的尺度因子a就是最佳尺度，記為a（b），信號在b點的瞬時頻率可以估計為f（b）＝1／a（b）。把變形條紋I2（x，y）過中心點O沿OAB方向的灰度分布作為f（x），按照上述過程可以提取畸變條紋徑向OAB的瞬時頻率分布f（x）＝1／a（x），經(jīng)過瞬時頻率積分得到徑向調(diào)制相位函數(shù)。

3 實驗結(jié)果

調(diào)整攝像機(jī)光軸垂直參考平面，然后采集圖像。圖2是成像系統(tǒng)采集的變形條紋（640pixel× 480pixel），中心點O即零畸變點位于（320，240）。圖3是OAB方向的瞬時頻率分布，以O(shè)為中心呈左右對稱。圖4是對O點右側(cè)瞬時頻率積分得到徑向調(diào)制相位的原始曲線和多項式擬合曲線。按照（4）式把圖4中的徑向調(diào)制相位轉(zhuǎn)化為徑向畸變位移，采用立方卷積插值法［15］對畸變圖像（見圖5）進(jìn)行校正，校正圖像如圖6所示。

Fig.2 Distorted fringe pattern

Fig.3 Local frequency along OAB

Fig.4 Modulated phase along OAB

Fig.5 Distorted image

Fig.6 Calibrated image

4 小 結(jié)

為了測量成像系統(tǒng)的徑向畸變，以正弦灰度條紋模板作為成像對象，把徑向位置畸變轉(zhuǎn)化為變形條紋徑向調(diào)制相位；該方法只需要采集1幅畸變條紋圖像，該變形條紋不僅含有徑向畸變信息，而且在畸變條紋中心點含有無畸變條紋像信息，采用小波頻率估計提取徑向頻率并進(jìn)行瞬時頻率積分計算徑向調(diào)制相位，實現(xiàn)畸變校正。

［1］ TANGW，YE D，YUAN F，et al.Independent distortion correction algorithm in vision measurement system［J］.Journal of Optoelectronics·Laser，2013，24（2）：308-315（in Chinese）.

［2］ JIN Y，MENG JB，WANG K，et al.The research of plane array camera geometric distortion correction method［J］.Microelectronics＆Computer，2011，28（10）：36-43（in Chinese）.

［3］ SU C Z，WANG E G，HAO JT，et al.Distortion correction for images in planarmetrology［J］.Optics and Precision Engineering，2011，19（1）：161-167（in Chinese）.

［4］ WU K X，DUAN M L.New method of correcting barrel distortion on latticemodel［J］.Journal of Computer Applications，2012，32（4）：1113-1115（in Chinese）.

［5］ KONG B，F(xiàn)ANG T J.A simple and precisemethod for radial distortion calibration［J］.Journal of Image and Graphics，2004，9（4）：429-434（in Chinese）.

［6］ YU J，LIN JM，YANG JY，et al.Distortion measurement and calibration technology for CCD camera system on new target board［J］.Acta Optica Sinica，2007，27（8）：1440-1442（in Chinese）.

［7］ XUW X.The study of lens distortionmeasurementof NFOV based on star pointmethod［J］.Laser Technology，2011，35（5）：593-595（in Chinese）.

［8］ YANG Ch P，WENG JW，YANG L L，et al.Grid fringe pattern Fourier transform profilometry［J］.Laser Technology，2010，34（4）：493-496（in Chinese）.

［9］ YANG Ch P，WENG JW，WANG JW.Distortion measurement and calibration technique based on phase analysis for carrier-fringe pattern［J］.Acta Photonica Sinica，2010，39（2）：316-319（in Chinese）.

［10］ LIB，QIAN X F，LIX H，et al.Phase-unwrapping algorithm based on radial shearing principle［J］.Laser Technology，2013，37（1）：44-47（in Chinese）.

［11］ SUN J，CHENW J，SU X Y，etal.Study themeasurement range ofwavelet transform profilometry［J］.Acta Optica Sinica，2007，27（4）：647-653（in Chinese）.

［12］ ZHONG JG，WENG JW.Phase retrieval of optical fringe patterns from the ridge of a wavelet transform［J］.Optics Letters，2005，30（19）：2560-2562.

［13］ YANG Ch P，LIU Y，LIH.Radial distortion distribution in optical imagingmeasured by wavelet phase analysis［J］.Acta Optica Sinica，2012，32（9）：0912001（in Chinese）.

［14］ WANW B，SU JH，YANG L H，etal.Phase unwrapping algorism for image processing of interferogram［J］.Journal of Applied Optics，2011，32（1）：70-74（in Chinese）.

［15］ YOU Y H，ZHOU X K.Research on optimal interpolation algorithm for digital image［J］.Chinese Space Science and Technology，2005，25（3）：14-18（in Chinese）.

M easurement of radial distortion based on frequency integration phase demodulation

YANG Chuping，LIU Jianbin，TAN Suiyan，WENG Jiawen
（Department of Physics，South China Agricultural University，Guangzhou 510642，China）

Tomeasure the radial distortion in optical imaging systems，adopting the straight sinusoidal carrier-fringe pattern，the radial modulated phase resulting from radial distortion was extracted by means of instantaneous frequency integration.The formula between the instantaneous frequency and the radial modulated phase was deduced，and the conversion formula between the radial modulated phase and the radial distortion displacement was obtained.The instantaneous frequency of the distorted fringe pattern in the radial direction was calculated by using frequency estimation of Gabor wavelet transform，and integrated to obtain the radialmodulated phase.The radialmodulated phase and the cubic convolution interpolation algorithm were used to calibrate the distorted image.Experimental results demonstrate that the method is available.

measurement and metrology；distortion measurement；carrier frequency fringe pattern；phase demodulation；radial distortion

O438

A

10.7510／jgjs.issn.1001-3806.2014.03.026

1001-3806（2014）03-0402-04

國家自然科學(xué)基金資助項目（61307011）；廣東省自然科學(xué)基金資助項目（9151064201000035）

楊初平（1970-），男，副教授，研究方向為光信息處理。

E-mail：yangchp＠scau.edu.cn

2013-07-09；

2013-10-08