中国激光
CHINESEJOURNALOFLASERSVol.38,No.2February,2011
基于解析图像的小波变换光学载频干涉
全息图相位重建方法
李思坤苏显渝陈文静
(四川大学光电科学技术系,四川成都,610064)
摘要将解析图像和实数小波变换(WT)相结合应用到光学载频全息图分析中。首先对载频全息图进行希尔伯特变换构造解析图像,然后对解析图像进行WT,提取小波脊处对应的WT系数的相位信息即可得到载频全息图的相位信息。由于母小波优秀的空域局部化能力,使得因全息图边缘不连续引入的误差被限制在局部区域。同时,由于不需要人工滤波操作,所提方法一定程度上解决了载频全息图的频谱混叠问题。给出了严格的理论公式推导、计算机模拟和实验验证。
关键词傅里叶光学;干涉计量;相位测量;小波变换
中图分类号O438.2文献标识码Adoi:10.3788/CJL201138.0209002
AnalyticImageBasedWaveletTransformMethodforPhaseReconstructionofOpticalInterferogramwithLinearCarrier
LiSikunSuXianyuChenWenjing
(DepartmentofOptoElectronics,SichuanUniversity,Chengdu,Sichuan610064,China)
AbstractAnalyticimagetechniqueandwavelettransform(WT)techniqueareappliedtoanalysisofinterferometryhologramwithalinearcarrier.Inthismethod,Hilberttransformisfirstlyperformedonthehologramtogettheanalyticsignals.AndthenWTiscalculatedonthesignal.Finally,thephaseinformationofthehologramcanbeobtainedfromtheWTcoefficientsatthewaveletridgeposition.Becauseofthegoodspatiallocalizationabilityofthemotherwavelet,theerrorcausedbytheboundaryofthehologramislimitedinlocalareas.Filteringprocessisavoidedinthismethod,soitcanresolvethefrequencyoverlappingproblemtoacertainextent.Mathematicaldemonstrationisgivenindetail.Computersimulationandexperimentsverifythevalidityoftheproposedmethod.KeywordsFourieroptics;interferometry;phasemeasurement;wavelettransformOCIScodes120.3180;120.5050;100.7410
1引言
干涉计量
[1]
是一种常用的非接触相位测量方
出有多种相位重构方法。例如,王超等提出了一种相位梯度积分的方法;Womack[3]提出了一种实时的方法,先把经过被测物体调制的全息图和参考全息图进行乘法运算,将待求相位成分转化为低频成分,然后进行低通滤波提取被测物体相位,由于需要滤波过程,该方法在图像边缘处有较大误差;为了减小这种方法的边缘误差,Mendiola等[4]提出了一
[2]
法,其中载频干涉只需要一幅干涉图就可进行相位重构,因而经常用于对瞬变场的测量。载频干涉计量通过物光和参考光之间大的倾斜在干涉图的某个方向上产生线性载频,然后利用相位重构方法从载频干涉图中重构出所求的相位。目前已有研究者提
收稿日期:20100511;收到修改稿日期:20100917
基金项目:国家自然科学基金(60838002)资助课题。
作者简介:李思坤(1982),男,博士研究生,主要从事三维面形测量和光信息处理等方面的研究。Email:lisk0409@126.com
导师简介:苏显渝(1944),男,教授,博士生导师,主要从事光学三维传感与机器视觉、光电信息技术及应用等方面的研究。Email:xysu@scu.edu.cn(中国光学学会会员号:6100105)
02090021
中国激光
种迭代的实时方法;Ferrari等[5]提出了一种相移重构相位方法,但这种方法从一幅全息图中得到的子图分辨率比较低。Takeda等提出了一种傅里叶变换方法,该方法首先对全息图进行傅里叶变换,通过人工设定一定的带通滤波器滤出包含了待求相位信息的基频部分,然后进行逆傅里叶变换、相位展开重建干涉图相位信息。该方法有较好的噪声抑制能力,可靠性高,实时性强,是一种常用的方法。但是,由于数据处理过程中需要人工设定滤波器进行滤波,会在全息图边缘处存在较大的误差,另一方面当全息图的零频成分、高次频成分和有用基频成分存在混叠时,人工很难设定合适的滤波器进行滤波获得完整的基频成分[7],这种情况下恢复相位结果常存在较大误差甚至错误。为了解决边界和频谱混叠问题,本文将解析图像和小波变换(WT)[8]结合起来应用到相位提取中。由于WT具有良好的空域局部化特性,使得边界误差被限制在局部区域,同时该方法无需人工滤波操作,一定程度上解决了频谱混叠问题。给出了这种方法在载频全息干涉图中应用的详细理论论证,计算机模拟和实验验证了所提方法的有效性。
[6]
直接从系数之中得到所需的相位信息。因此首先通过希尔伯特变换构造被分析信号的解析信号。对全息图信号f(x)的希尔伯特变换可以表示为
Hfx
=1 -
!
!
f!d!,x-!
(4)
然后构造全息图信号f(x)对应的解析信号,表示为I(x)=f(x)+jHf(x)=A+Bcos(x)+
jBsin(x)=A+Bexpj(x),(5)式中∀j#代表虚数单位。在位置b处对相位函数(x)进行泰勒级数展开得到
(x)=(b)+∃(b)(x-b)+
%(b)(x-b)2+&,
2!
(6)
对大多数的载频条纹图而言∃(b)是缓慢变化的,可以认为%(b)∋%∃(b)∋&∋0,所以可以忽略(6)式中(x-b)的高阶项。对解析信号I(x)进行WT,得到
+!
W(a,b)=
+!
-!+!
I(x)a,b(x)dx=
*
a,b(x)dx=
*
2原理
光学条纹图通常具有以下形式
[6]
-!
+!
A+Bexpj(x)
f(x)=A(x)+B(x)cos(x),(1)
式中(x)是全息图的相位函数,A(x)表示背景光强分布,B(x)代表调制度分布,A(x)和B(x)相对于载频信息通常是缓慢变化的,可以分别看作常数A和B。
对信号I(x)的WT可以表示为
+!
-!
*Aa,b(x)dx+
-!
*
Bexpj(x)a,b(x)dx,(7)
采用墨西哥草帽小波作为母小波,根据小波的容许条件知(7)式中第一项等于零,所以(7)式可以进一
步表示为
+!
W(a,b)=
-!
Bexpj(x)a,b(x)dx=
1a∃(b)2
2
*
W(a,b)=
-!
I(x)(x)dx,
*a,b
(2)
2 Bc0a2∃(b)2exp-
expj(b),
(8)
式中W(a,b)为WT系数,a>0是和频率有关的伸
缩因子,b代表和位置x有关的位移因子,∀*#表示共轭。设(x)代表母小波,则母小波经过伸缩和平移可以得子小波a,b(x)=
1a
x-ba
。这里选
因此,WT系数模值W(a,b)对伸缩因子a的一
阶导数可以表示为
dWa,bda
=
2 Bc0a∃(b)2(
2
择空域局部化能力较好的墨西哥草帽小波作为母小波。墨西哥草帽小波是一个实数小波,其定义式为
12
(x)=c01-x2exp-x,(3)
2式中c0=
2为一个常数。当母小波是实值小3 1/4
02090022
exp-1a∃(b)
2当dW(a,b)da得到
ar(b)=
2-a2∃(b)2,
(9)
=0时W(a,b)具有最大值,这时
2=21,2 fb∃(b)(10)
波时,信号的WT系数仍然是一个实数信号,不能式中ar(b)表示位置b处WT脊处的伸缩因子,fb=
李思坤等:基于解析图像的小波变换光学载频干涉全息图相位重建方法
∃(b)/(2 )是所谓的瞬时频率。在位置b处WT系数W(a,b)的模值最大处定义为小波脊,所以I(x)的小波脊处的变换系数可以表示为W(ar,b)=2
2 Bc0exp(-1)expj(b),
(11)
由于b是和x对应的位置因子,可知小波脊处的WT系数的相位和光学载频全息图的相位相等。所以可以通过这种解析信号和WT结合的方法提取全息图的相位信息。
(b)=arctanImWar,bReWar,b,
(12)
小二乘拟合[12],得到其中等效线性相位∀0x,然后把变形全息图相位和线性载频相位相减就可以解得物体调制相位。
3计算机模拟
通过计算机模拟验证理论分析的正确性。模拟的物光相位函数和全息图调制度[分别如图1(a)和(b)所示]的数学表达式分别为
3(x,y)=-x-300/1002
1 2
y-300/100232
-(13)
2
这里得到的相位截断在(- , )之间,以一定的方法对截断相位进行展开[9~11]就可以得到连续的相位信息。
在光学载频全息图中,以∀0表示线性载频。移走被测相位物体,拍摄一幅参考全息图,用同样的方法可以提取出连续相位∀0x,或者对(x,y)进行最
,
+
B(x,y)=0.72-0.032
y/100-3
,
x/100-3
(14)
模拟全息图背景光强A=0.24,线性载频∀0= /3,则模拟干涉载频全息图可以表示为
f(x,y)=0.24+Bcos∀0x+#(x,y).(15)
图1(a)模拟调制相位,(b)调制度,(c)载频干涉全息图Fig.1(a)Simulatedphase,(b)modulation,(c)hologram
全息图像大小为600pixel(600pixel,如图1(c)所示。线性载频成分受到被测物体调制而发生变形。首先对全息图逐行进行希尔伯特变换构造对应的解析图像,然后对得到的解析图像逐行进行WT,这里选择墨西哥草帽小波作为母小波。图2(a)为对解析图像第150行进行WT后的系数模值图像,在每个位置处沿着尺度因子方向模值最大处为小波脊,即图2(a)中所示每个位置的最亮处。图2(b)为对应WT系数的截断相位分布。图2(c)为得到整幅图像的小波脊处WT系数的模值后根据(11)式计算得到的调制度信息。对比图2(c)和图1(b)可见文中所述方法可以重构干涉图调制度成分。提取小波脊位置处的截断相位,以一定的相位展开方法进行相位展开,得到的连续相位分布如图2(d)所示,图2(e)所示为点对点误差分布,可见采用本文提到的方法可以高精度地解调全息图相位。同时,模拟中对全息图调制度信息和调制相位信息的精确重构验证了理论推导部分的正
确性。
4实验
用马赫曾德尔干涉光路,并有意使干涉仪的两个臂发生倾斜,以便在水平方向上产生线性载频。使用HeNe激光器产生波长为∃=633nm的激光作为光源。在干涉仪的其中一个臂上放置一个薄透镜作为待测相位物体。CCD摄像机分辨率为800pixel(600pixel,拍摄的载频干涉全息图经过剪切后如图3(a)所示,全息图的强度分布很不均匀。干涉图第150行的傅里叶频谱如图3(b)所示,零频成分和有用基频成分严重混叠在一起,这使得用傅里叶变换法对全息图进行相位提取的时候,很难设定合适的滤波器滤出完整的基频成分。全息图边缘不连续对傅里叶分析法的影响很大,通常会造成边缘部分存在较大误差。为了减小边缘误差通常选择汉宁滤波窗进行滤波。这里相应选择汉宁滤波窗。对全息图的频谱滤波后进行逆傅里叶变换,得
02090023
中国激光
图2(a)WT系数模值,(b)WT系数相位,(c)重建调制度分布,(d)重建调制相位,(e)相位误差Fig.2(a)Modulusofwavelettransform(WT)coefficients,(b)phaseofWTcoefficients,
(c)reconstructedmodulation,(d)retrievedphase,(e)errorofphase
图3(a)实拍的全息图,(b)第150行的频谱,(c)傅里叶变换方法得到的截断相位和(d)重建的相位,
(e)本文所提方法解得的截断相位和(f)重建的相位
Fig.3(a)Arealhologram,(b)spectrumofthe150throw,(c)wrappedphaseand(d)reconstructedphaseby
usingFFTmethod,(e)wrappedphaseand(f)reconstructedphasebyusingtheproposedmethod
02090024
李思坤等:基于解析图像的小波变换光学载频干涉全息图相位重建方法
到的截断相位如图3(c)所示。在边缘部位仍然存在较大的误差。以菱形相位展开算法进行相位展开后重建的连续相位如图3(d)所示,边缘误差很大。然后以本文所提方法进行相位提取,得到的截断相位如图3(e)所示。对比图3(c)和(e)可见,本文所提方法具有明显的优势。以同样的相位展开算法进行相位展开,重建的连续相位如图3(f)所示。由于WT母小波优秀的空域局部化能力,使得边缘不连续引入的误差被限制在边缘局部区域内。同时本文所提方法不需要人工滤波过程,在一定程度上解决了频谱混叠问题。
[9]
5结论
将解析信号和实数WT相结合可以对光学载频全息干涉图进行高精度解调分析。通过严格的理论分析论证了解析图像小波脊处的WT系数的相位信息和全息图的相位信息相等的关系,并通过对模拟的和实际拍摄的全息图进行相位提取验证了所提方法的有效性和可靠性。由于WT母小波优秀的空域局部化能力,使得由于全息图边缘不连续引入的误差被限制在边缘的局部区域,提高了这些区域相位重建精度;同时由于不需要人工滤波操作,此方法一定程度上解决了载频全息图的频谱混叠问题。本文所述方法无需已知物光和参考光的振幅分布,对条纹图强度变化不敏感,只需一幅或者两幅全息图就可解得到物体调制相位,适用于瞬变场的测量。另外,该方法还可以同时分析全息图的强度分布信息。边缘误差可以采用已有的图像延拓方法处理(例如零延拓或者Gerchberg迭代延拓),将边缘误差外推解决。
致谢感谢四川大学电子信息学院冯国英教授和王超硕士在研究过程中的交流。
参
考
文
献
Press,1998.396~429
金国藩,李景镇主编.激光测量学[M].北京:科学出版社,1998.396~429
2WangChao,FengGuoying.Anewalgorithmforphasereconstructionfromasinglecarrierfrequencyinterferogram[J].ActaOpticaSinica,2008,28(7):1269~1273
王超,冯国英.从一幅载频全息图中实现相位重构的新算法[J].光学学报,2008,28(7):1269~1273
3K.H.Womack.Interferometricphasemeasurementusingspatialsynchronousdetection[J].Opt.Eng.,1984,23(4):391~3954J.YaezMendiola,M.Servn,D.MalacaraHern ndez.Iterativemethodtoobtainthewrappedphaseinaninterferogramwithalinearcarrier[J].Opt.Commun.,2000,178(46):291~296
5J.A.Ferrari,E.M.Frins.Multiplephaseshiftedinterferogramsobtainedfromasingleinterferogramwithlinearcarrier[J].Opt.Commun.,2007,271(1):59~646M.Takeda,H.Ina,S.Kobayashi.Fouriertransformmethodoffringepatternanalysisforcomputerbasedtopographyandinterferometry[J].J.Opt.Soc.Am,1982,72(1):156~1607YangHu,ChenWenjing,LuChengqiangetal..InfluenceofsamplingonFouriertransfromprofilometry[J].ActaOpticaSinica,1999,19(7):929~934
杨虎,陈文静,陆成强等.抽样对傅里叶变换轮廓术的影响[J].光学学报,1999,19(7):929~934
8ZhouXiang,ZhaoHong.ThreedimensionalprofilometrybasedonMexicanhatwavelettransform[J].ActaOpticaSinica,2009,29(1):197~202
周翔,赵宏.基于Mexicanhat小波变换的三维轮廓术[J].光学学报,2009,29(1):197~202
9ZhangQican.TechnicalStudyofThreeDimensionalShapeMeasurementforDynamicProcess[D].Chengdu:SichuanUniversity,2005.20~25
张启灿.动态过程三维面形测量技术研究[D].成都:四川大学,2005.20~25
10LiSikun,ChenWenjing,SuXianyu.Phaseunwrappingguidedbyamplitudeofwaveletridgecoefficientsinwavelettransformprofilometry[J].ActaOpticaSinica,2008,28(4):715~721李思坤,陈文静,苏显渝.小波变换轮廓术中用小波脊系数幅值引导相位展开的研究[J].光学学报,2008,28(4):715~721
11LiYong,SuXianyu.Fastalgorithmforreliabilityguidedphaseunwrapping[J].OptoElectronicEngineering,2005,32(11):76~79
李勇,苏显渝.用于可靠性导向的相位展开的快速算法[J].光电工程,2005,32(11):76~79
12LiuYuankun,SuXianyu,WuQingyang.MulticameracalibrationbyFTPtechnique[J].ActaPhotonicaSinica,2007,36(9):1734~1737
刘元坤,苏显渝,吴庆阳.基于傅里叶条纹分析的多摄像机标定方法[J].光子学报.2007,36(9):1734~1737
1JinGuofan,LiJingzhen.LaserMetrology[M].Beijing:Science
02090025
因篇幅问题不能全部显示,请点此查看更多更全内容