光学信息处理讲义

光学信息处理1.引言自六十年代激光出现以来，光学的重要发展之一是形成了一个新的光学分支——傅里叶光学。傅里叶光学是指把数学中的傅里叶分析方法用于波动光学，把通讯理论中关于时间、时域、时间调制、频率、频谱等概念相应地改为空间、空域、空间调制、空间频率、空间频谱，并用傅里叶变换的观点来描述和处理波动光学中学波的传播、干涉、衍射等。傅里叶变换已经成为光信息处理的极为重要的工具。光学信息处理就是对光学图像或光波的振幅分布作进一步的处理。自从阿贝成像理论提出以后，近代光学信息处理通常是在频域中进行。由于光的衍射,图像的夫琅和费衍射分布,即图像的空间频谱分布与图像的空间分布规律不同,这使得在频谱面上对其进行处理可获得一些特殊的图像处理效果。近代光学信息处理具有容量大，速度快，设备简单，可以处理二维图像信息等许多优点，是一门既古老又年青的迅速发展的学科。光学信息存储、遥感、医疗、产品质量检验等方面有着重要的应用。2.实验目的1)通过实验，加强对傅里叶光学中有关空间频率、空间频谱和空间滤波等概念的理解。2)掌握光学滤波技术，观察各种光学滤波器产生的滤波效果，加深对光学信息处理基本思想的认识。3)加深对卷积定理的理解4)了解用光栅滤波实现图像相加减及光学微分的原理和方法。5)了解黑白图像等密度的假彩色编码。3.实验原理1)二维傅里叶变换和空间频谱在信息光学中常用傅里叶变换来表达和处理光的成像过程。设在物屏X-Y平面上光场的复振幅分布为g(x，y)，根据傅里叶变换特性，可以将这样一个空间分布展开成一系列二维基元函数的线性叠加，即)](2exp[yfxfiyx+π∫∫+∞∞−+=yxyxyxdfdfyfxfiffGyxg)](2exp[),(),(π（1）式中fx、fy为x、y方向的空间频率，即单位长度内振幅起伏的次数，G(fx，fy)表示原函数g(x，y)中相应于空间频率为fx、fy的基元函数的权重，亦即各种空间频率的成分占多大的比例，也称为光场（opticalfield）g(x，y)的空间频谱。G(fx、fy)可由g(x，y)的傅里叶变换求得∫∫+∞∞−+−=dxdyyfxfiyxgffGyxyx)](2exp[),(),(π(2)g(x，y)与G(fx，fy)是一对傅里叶变换式，G(fx，fy)称为g(x，y)的傅里叶的变换，g(x，y)是G(fx，fy)的逆变换，它们分别描述了光场的空间分布及光场的频率分布，这两种描述是等1效的。当g(x，y)是空间周期函数时，空间频率是不连续的。例如空间周期为x0的一维函数g(x)，即g(x)=g(x+x0)。描述空间周期为x0的一维光栅时，光栅面上光振幅分布可展成傅里叶级数∑∑==)2exp()2exp()(0xnfiGxfiGxgnnnππ(3)上式中，n＝0，±1，±2，……；f0=1/x0，称为基频；fn=nf0，是基频的整数倍频，称为n次谐波的频率。Gn是g(x)的空间频率，由傅里叶变换得dxxnfixgxGxxn)2exp()(102/2/000π−=∫+−(4)2)透镜的二维傅里叶变换性质在光学上，透镜是一个傅里叶变换器，它具有二维傅里叶变换的本领。理论证明，若在焦距为F的正透镜L的前焦面（X-Y面）上放一光场振幅透过率为g(x，y)的物屏，并以波长为λ的相干平行光照射，则在L的后焦面（X′-Y′面）上就得到g(x，y)的傅里叶变换，即g(x，y)的频谱，此即夫琅禾费衍射情况。其空间频谱∫∫+∞∞−′+′−=′′dxdyyFyxFxiyxgFyFxG)](2exp[),(),(λλπλλ(5)其中空间频率fx、fy与透镜像方焦面（频谱面）上的坐标有如下关系Fxfxλ/′=Fyfyλ/′=（6）显然，),(FyFxGλλ′′就是空间频率为（FyFxλλ′′,）的频谱项的复振幅，是物的复振幅分布的傅里叶变换，这就为函数的傅里叶变换提供了一种光学手段，将抽象的函数演算变成了实实在在的物理过程。由于FyFxλλ′′,分别正比于x′，y′，所以当λ、F一定时，频谱面上远离坐标原点的点对应于物频谱中的高频部分，中心点x′=y′=0，fx=fy=0对应于零频。3)阿贝成像原理现在我们知道，物体应该看成是大量空间信息的集合体，光信息处理涉及的空间信息的频谱不再是一个抽象的数学概念，它是展现在透镜焦平面上的物理实在。1873年，德国人阿贝从波动光学的观点提出了一种成像理论。他把物体或图片看成包含一系列空间频率的衍射屏，在相干光照明下，物体通过透镜成像的过程分为两步：第一步是通过物的夫琅禾费衍射，在物镜后焦面上形成一个衍射图样，第二步是这些衍射图样发出的子波相干涉，在像平面上相干迭加形成物的像，通过目镜可以观察到这个像。ABCLFYHA′B′C′G平行光图1阿贝成像原理图2按频谱分析理论，谱面上的每一点均具有以下四点明确的物理意义：（1）谱面上任一光点对应着物面上的一个空间频率成分。（2）光点离谱面中心的距离，标志着物面上该频率成分的高低，离中心远的点代表物面上的高频成分，反映物的细节部分。靠近中心的点，代表物面的低频成分，反映物的粗轮廓。中心亮点是0级衍射即零频，反映在像面上呈现均匀背景。（3）光点的方向，指出物平面上该频率成分的方向，例如横向的谱点表示物面有纵向栅缝。（4）光点的强弱则显示物面上该频率成分的幅度大小。由以上定性分析可以看出，阿贝的二次成像理论的第一次衍射是透镜对物作空间傅里叶变换，它把物的各种空间频率和相应的振幅一一展现在它的焦平面上。一般情况下，物体透过率的分布不是简单的空间周期函数，它们具有复杂的空间频谱，故透镜焦平面上的衍射图样也是极复杂的。第二次衍射是指空间频谱的衍射波在像平面上的相干迭加。如果在第二次衍射中，物体的全部空间频谱都参与相干迭加成像，则像面与物面完全相似。如果在展现物的空间频谱的透镜焦平面上插入某种光学器件（称之为空间滤波器），使某些空间频率成分被滤掉或被改变，则像平面上的像就会被改变，这就是空间滤波和光学信息处理的基本思想。在实际光学成像系统中，像和物不可能完全一样。这是由于透镜的孔径是有限的，总有一些衍射角比较大的高次光线（高频信息）不能进入物镜而被丢掉。所以像的信息总是比物的少些。由于高频信息主要反映物的细节，因此，无论显微镜有多大的放大倍数，也不可能在像面上分辨出这些细节。这是限制显微镜分辨本领的根本原因。当物镜孔径极其小时，有可能只有零级衍射通过物镜，这时像面上有亮的均匀背景而无像分布。4)空间滤波和光信息处理光信息处理是通过空间滤波器来实现的，所谓空间滤波器是指在图1中透镜的后焦平面上放置某种光学元件来改造或选取所需要的信息，以实现光信息处理。这种光学器件称为空间滤波器。（b）高通（c）带通（d）方向（a）低通图2空间滤波器图2给出了几种常用的空间滤波器。（a）低通滤波：目的是滤去高频成分，保留低频成分。由于低频成分集中在谱面的光轴（中心）附近，高频成分落在远离中心的地方，经低通滤波后图像的精细结构将消失，黑白突变处也变得模糊。（b）高通滤波：目的是滤去低频成分而让高频成分通过，其结果正好与低通滤波相反，使物的细节及边缘清晰。（c）带通滤波：根据需要，有选择的滤掉某些频率成分。（d）方向滤波：只让某一方向，例如纵向的频率成分通过，则像面上将突出了物的横向线条。5)卷积一个二维函数的卷积定义为：ηξηξηξ∫∫−−=∞∞−ddyxhfyxhyxf),(),(),(*),((7)3卷积本身的概念比较抽象，卷积运算比较复杂。但是在理论上可以证明：两个函数乘积的傅里叶变换，等于它们各自的傅里叶变换的卷积。反之，两个函数卷积的傅里叶变换，等于它们各自傅里叶变换的乘积。这个卷积定理指出：傅里叶变换可以化复杂的卷积运算为简单的乘积运算,从而提供了计算卷积的一种简单手段，并利用数字计算机快速的算出(其算法称为快速傅立叶变换算法(FFT))。6)图像相加减利用正弦光栅作为空间滤波器，可在相干光学处理系统中对图像进行实时的相加和相减运算。如图3所示，以正弦光栅作为滤波器，将其置于4F系统的滤波平面上。正弦光栅的复振幅透过率为：)2cos(),(0010ϕπν++=ufGGug(8)式中Ffuλ=,u,ν为场点坐标；F为傅里叶变换透镜的焦距；0f为光栅频率；f是与场点坐标u对应的频率变量；0ϕ表示光栅条纹的初位相，它决定了光栅相对于坐标原点的位置。将图像A和图像B置于物面上。对于每一种物图像U0(x,y)，在4F配置的系统中，它的像函数为物函数的对称函数U0(-x,-y)与),(νug的傅里叶变换的卷积(省略共有的系数)。经过计算，像的函数为00)),((21)),((21),(),('00100100ϕϕλλiieyFfxUGeyFfxUGyxUGyxU−−+−+−−−+−−=(9)即在像面上呈现三个像，分别为0级像，-1级像，+1级像，它们都是原的的再现。因此，在此系统中A，B各有三个像A0，A-1，A+1，B0，B-1，B+1，调节物面上A，B的距离b，使像面上A+1与B-1图像重合，此时A，B距离满足Ffb02λ=(10)当光栅条纹的初位相00ϕ=，即光栅条纹与轴线重合时，即A+1与B-1像之间没有相差，则实现了复振幅相加运算。当光栅条纹的初位相02πϕ=，即光栅条纹偏离轴线14周期时，A+1与B-1像的位置不变，但产生相反的相移，位相差为π，则实现了复振幅相减运算。7)光学微分光学微分不仅是一种重要的光学—数学运算，在光学图像处理中也是突出信息的一种重要方法。在图象识别技术中，突出图象的边缘是一种重要的识别方法。对于一张比较模糊的图象，由于突出了其边缘轮廓而变得易于辨认。为了突出图象的边缘轮廓，我们可以用空间滤波的方法，去掉图象中的低频成分而突出图象的高频成分，从而使轮廓突出。虽然利用高通滤波可使像边缘增强，但由于光能量损失太大，因而使像的能见度大大降低，减弱了信号，利用光学微分法可以得到较满意的结果。复合光栅可以用来作为空间滤波器实现图像微分。所谓的复合光栅是指在同一块全息干板上制作的两个栅线平行但空间频率稍许不同的光栅。为书写简单，设其初始位置时的振幅透过率函数为)2cos()2cos(),(22110ufGufGGugππν++=(11)光学微分的光路系统仍采用4F系统，待微分的图像置于输入面的原点位置，微分滤波器置于频谱面上。作和上面一样的分析，此时像面上将有五幅图像，与G0对应的0级像，与G1，f1相联系的为±1级A±1，与G2，f2相联系的为±级A2±1,两对±1级像中心相对0级中4心的距离分别为2211,FfxFfxλλ=Δ=Δ(12)由于f1,f2之间存在差额2112,)(fffff−=Δ，使得A+1像与A2+1像之间或A-1像与A2-1像之间略有位错，其数值为fFxxxΔ=Δ−Δ=Δλ12(13)这就实现了微分运算的第一步。为实现二者相减，按照图相减同样的做法，将复合光栅平移一合适距离，引起二张图像的相移uΔuffΔ−=)(212πϕ(14)令πϕ=，得fuΔ=Δ21。此时在±1级像所在处的场为010102121)'(21)'(21)','(UGexUGxxUGyxUiΔ=+Δ+=+ϕ(令21GG=)(15)这样就实现了光学微分运算。8）θ调制θ调制实验是对阿贝的二步成像理论的一个巧妙应用。将一个物体用不同的光栅来进行编码，制作成θ片。将θ片置于白光照明中，在频谱面上进行适当的空间滤波处理，便可在输出面上得到一个假彩色的像。我们知道，如果在一个透镜的前面放置一块光栅并用一束单色平行光垂直的照射它，在透镜的后焦面（即频谱面）上就会形成一串的衍射光斑，其方向将垂直于光栅的方向。如果有一个二维的图形，其不同部分由取向不同光栅制成（调制），显而易见，他们的衍射光斑也将有不同的取向，即在透镜的后焦平面（频谱面）上，各部分的频谱分布也将有不同，如果我们挡住某一部分的频谱，在频谱面后的这部分图象将会消失，可见，输入图象中各部分的频谱，只存在于调制光栅的频谱点附近。如果我们用白光照射θ片，则在频谱上可得到彩色的频谱斑（色散作用），每个彩色斑的颜色