小波分析在光学信息处理中的应用.docx

小波分析在光学信息处理中的应用马晶 谭立英 冉启文用小波变换分析惠更斯-菲涅耳原理，为小波光学理论描述理论基础.运用该理论，对光学系统的空域滤 波、空频域滤波现象进行分析，其中小波在空域滤波以空间可变处理为例，小波在空频域滤波以匹配滤波为 例，波前滤波则以夫琅和费单缝、圆孔等为例，分析了小波变换在光学领域处理问题的可行性，初步建立了 小波光学的理论框架.PACC：4200；4230APPLICATION OF WAVELET ANALYSIS TO OPTICAL INFORMATION PROCESSING*MA JING TAN LI-YING(State Key Laboratory of Tunable Laser Technology,Harbin Institute of Technology,Harbin 150001) RAN QI-WEN(Department of Mathematics,Harbin Institute of Technology,Harbin 150001)ABSTRACT In this paper,the Huygens-Fresnel principle is restudied by means of wavelet analysis firstly,this builds the basis of the theory for wavelet optics.Then the theory of wavelet analysis is used to analyze the filtering in spatial domain and spatial frequency domain.For example,it gives the space variable processing for the filtering with wavelet in spatial domain,the match filtering for the filtering with wavelet in spatial frequency domain and a single slit and a round aperture for the wave-front filtering.The feasibility for treating the problems in optical field with wavelet transform is discussed,and the theoretical frame of wavelet optics is built.PACC：4200;42301 引言小波变换是一种新的有效的分析手段.它能有效地从光学图像(信号)中提取信息，通过伸缩和平移等运 算功能对光学图像进行多尺度的细化(multiscale analysis) .近年来，人们尝试用光学系统实现小波变换.Szu等人提出了利用二维相关器实现一维小波变换的方案⑵3] .Mendlovic等人提出了一个利用Dammann 光栅的多通道系统[4],其中一个通道实现一个尺度的小波变换.Wang等人采用体全息技术，利用角度复用实 现了光学Harr小波变换，并用于图像特征提取⑸.Bloch等人讨论了利用计算全息方法实现光学小波变换的 方案[6].此外还有利用干涉计量仪实现光学小波变换的方案[7].以上这些工作的目的都是用光学系统实现小波变换，而本文则侧重用小波分析方法和理论来建立新的光 学理论一一小波光学•传统的光学信息处理理论，在数学手段上借助Fourier变换，且认为光学系统中光信 号的传递是平稳过程.而实际上由于外界条件干扰，图像在光学系统中的传递是非平稳的，因此用小波变换 分析更符合实际的光学系统.本文用小波变换分析惠更斯-菲涅耳原理，提出了小波波前滤波思想，为小波光 学理论描述了理论基础.同时运用该理论，对光学系统的空域滤波、空频域滤波现象进行分析，其中小波在 空域滤波以空间可变处理为例，小波在空频域滤波以匹配滤波为例，波前滤波则以夫琅和费单缝、圆孔等为 例，分析了小波变换在光学领域处理问题的可行性，初步建立了小波光学的理论框架.2 光波传播过程的小波波前滤波惠更斯-菲涅耳原理是光的衍射理论的基础，但也可以从另一角度研究光的衍射现象，即认为光波在传 播过程中其波阵面受到限制，衍射孔径对波前进行滤波(波前滤波)，滤波后的传播情况可以看成受限孔径处 光波在各点的权重重新分配，这里引入权重函数h (x,y)，并选其为母小波a,b;c,d(1)其中a,c是尺度因子，b,d是平移因子.此时受限衍射孔径处(x,y)点的光场分布为该光波各点权重的叠加结 果，即(3)其中U'(x,y)为衍射孔径上点(x,y)处未加权时的光场分布函数，U (x,y)可看成U'(x,y)与h (x,y)的卷积a,b;c,d a,b;c,d则在空频域中(4)其中'■ ■ ■ ■ ' ' 1. ' 1 分别是 U (x,y),U(x,y),h (x,y)的 Fourier 变换，亦a,b;c,d a,b;c,d即屏上的光场分布是受限孔径处光场分布与受限孔径处权重因子(母小波)的Fourier变换的乘积，而母小波 的选择应根据具体情况而定. 因此，从广义滤波系统而言，衍射现象也可看作是衍射孔径对波前的一种滤波. 下面以衍射现象为例，说明小波波前滤波思想的正确性.2.1 矩形孔径夫琅和费衍射 考虑一矩形孔径，其振幅透射率为[8]其中l ,l分别是x,y方向上的孔径宽度.xy用一单位振幅的单色平面波垂直照明，则孔径上的场分布即为由⑷式有： ' ' ' ：，■■■ 1 \所以只要找到一种母小波，使其在有效屏上门'•厂I，其余部分为零，就与实际情况完全一致.2.2 圆孔、正弦型位相光栅(夫琅和费衍射)此时仍用一单位振幅的单色平面波垂直照明，则孔径处的场分布U'(x,y)仍等于透射率函数t(x,y), 只是t(x,y)的表达式有所不同.对圆形孔径，有(7)而对正弦型位相光栅，则有(8)与前面对矩形孔径的分析一样，在这里仍需找到一种母小波，使其在有效屏上(P,q)=l,其余部分为零，就 与实际情况完全一致.从上面的几种特殊情况可知，对夫琅和费衍射，只需找到一种母小波，使其Fourier变换(p,q)(频率响应)在有效屏上使(p,q)=l,其余部分为零，由此要求选如下形式小波⑼：(9)其在频域中(10)上面几个例子都是对称光学系统，而实际上(2)式并不要求光学系统必须是对称的.从以上具体例子可以 看出，小波波前滤波理论是正确的. 由于惠更斯-菲涅耳原理是光学理论的基础，这样可以预计，光学系统的 各种信号处理和滤波都可以用小波分析的方法来处理. 因此可认为该理论奠定了小波光学理论基础. 下面以 此理论为基础，分别对光学系统的空域滤波、空频域滤波进行分析.3 空域小波滤波若在输入平面P放置一复透明滤波片h (x,y)作为空域滤波(如图1),于是在紧靠P平面的复光场为1 a,b;c,d 1以.厂7)= /.； I、.门『•工也.“•工hh飞(11)其中用二维母小波h (x,y)描述透明滤波片，f(x,y)是输入光信号，k是系数.a,b;c,d图 1 空域滤波系统在光学系统中，若把输入处理信号f(x,y)的空间坐标反转过来，并在输入空间域横向平移至(X0,y°)点, 即由 f(x -x,y -y) 代替 f(x,y) ，则00(12)这就是 f(x,y) 和 h (x,y) 的卷积积分.a,b;c,d 在空频域中，有(13)其中 1 ■'' ' 1' ' ■ . - ■ ■ 1 i1 ' 1 1 ■.分别是 E(x,y),f(x,y),h (x,y)的 Fourier 变换.a,b;c,d这样根据小波波前滤波理论，用小波变换实现了光学系统空域滤波， 并可通过扫描机构实现各种尺度下的滤波. 同时，由于小波变换的局域化 特征，可根据需要对任一点的局部特征进行提取，而且在设计空域滤波 器时，可在理论上根据需要对各个局部进行设计，并通过智能化控制使 得空域滤波灵活方便，这是Fourier变换所做不到的.4 空域小波滤波与空间可变处理 现在通过举例说明小波空域滤波理论的正确性.例如对空域的空间可变处理，由(11)式，输入-输出关系可由下式描述：以 < '.V> = \\ ./ Cl T V 也―..；(J' TV k.M 八(14)其中空间脉冲响应用小波母函数h (x,y)表示，f(x,y)是输入函数，g(x,y)是输出函数.--'■■ ■■■ - - ' ■■ 1 ■ ' (⑸其中空间可变处理反映在小波函数的a,c伸缩因子和b,d平移因子上.我们可制作一光栅，它的光栅频率是空间坐标的函数，换言之，对一给定的物点(x,y)，光栅某一局部的角空间频率(p ,q )将会使该物点变换到 xy频谱面上的点(a ,0 )，且~ ' - (16)其中入是光源的波长，f是透镜的焦距.5 空频域小波滤波在空频域内(图2P平面)也可以合成一个所需要的线性滤波运算， 其空间滤波器是由置于相干光学系统空频域的透明滤波片组成.如图2， 若待处理信号f(x,y)置于空域P，则输入信号的Fourier变换分布在空 频域P平面上.若有一透明滤波片h (p,q)置于P，于是紧靠P的背2 a,b;c,d 2 2后复光场分布为二声5 q ”5. w(办:p d ：(17)在输出平面p3上的光场分布为(18)(19)这里用小波母函数表示透明滤波片，空频域小波滤波器可由任意形状的 孔径和狭缝组成低通、高通、带通小波滤波器.并可根据小波的特性，对 输入信息的局部进行提取、识别，还可根据其空频域局部的变换情况， 找出输入信号奇异点的变化情况，而通过调整和移动空频滤波器，可实 现各种尺度下的空频域小波滤波.尤其当输入光信号是非平稳信号时，只 有用小波变换才有可能进行滤波处理.图 2 空频域滤波系统6 空频域小波滤波与匹配滤波器下面基于空频域小波滤波思想，以空频域匹配滤波说明光学系统空 频域小波滤波理论的正确性.将用母小波表示的匹配滤波器放在空频域P平面上，假定输入待处2理信号混杂一些相加性的高斯噪声，同时还假定2该小波匹配滤波器单独对信号S（x,y）的响应为So（x,y），而对随机噪声n（x，y）的响应为n0（x，y），有(20)f(x,y)=S(x,y)+n(x,y).滤波器在 x=y=0 时的输出信噪比为其中b 2是输出噪声的均方值.由(14)式(22)(23)由母小波表示的匹配滤波器，目的是要确定一个滤波函数，使其信噪比为最大，而由母小波函数的局域性及紧支撑性，可实现此目的.7 结论通过上面的分析可得到如下结论：1. 小波波前滤波思想是正确的，它与传统理论并不矛盾，由此也奠 定了小波光学信息处理的基础.2. 小波空域滤波可实现光学系统中的各种空域滤波，并由此分析了 空域可变空间处理，而用小波变换实现空间可变处理将比 Fourier 变换 更具灵活性.3. 小波空频域滤波在光学系统中可进行各种空频域滤波，据此讨论 了空频域小波匹配滤波器，由小波的局部特性实现匹配滤波.从上面的讨论初步看出，小波变换适合于光学信息处理系统，并将 给光学信息处理。