8-空间滤波解析

Chapter8空间滤波空间滤波：在光学系统的傅立叶频谱面上放置具有一定透过率的模片，改变光场的频谱结构，按照预定的需要使像质得到改善。•空间滤波技术基于傅立叶变换的理论基础，因而可称为傅立叶综合技术。空间滤波器1873年德国科学家阿贝创建的二次成像理论及其相应的实验是空间滤波与光学信息处理的先导。1935年泽尼克(Zernike)发明了相衬显微术并且第一次成功地观察到微小的相位物体－细菌，极大地推动了光学信息处理技术的发展。(1953年诺贝尔物理学奖)1946年法国科学家杜费（Duffieux）把光学成像看成线性变换，采用傅立叶方法成功地分析成像过程，发表了《傅立叶变换及其在光学中的应用》1953年艾里斯(Elias)《光学处理的傅立叶方法》，《光学和通讯理论》，1956年奥尼尔《光学中的空间滤波》为光学信息处理提供了有力的数学工具。60年代提出复数空间滤波概念，由于激光的出现和全息术的重大发展，光学信息处理进入了蓬勃发展的时期。光学信息处理技术在图像的修正处理、特征识别、孔径天线雷达数据的处理等方面得到了很好的应用。20世纪80年代以后，光学信息处理进入了一个更加活跃的时期。8.1空间滤波的基本原理8.2系统与滤波器8.3空间滤波应用举例8.4傅立叶变换透镜8.1空间滤波的基本原理一、阿贝二次衍射成像理论第一步，入射光场经物平面发生夫琅禾费衍射，在透镜的后焦面上形成一系列衍射斑(频谱)，称为第一次衍射像；第二步，衍射斑作为新的相干波源，由它发出的次波在像平面上相互叠加而形成物体的像，该像称为第二次衍射像。把物体看作是不同空间频率信息的集合，认为相干照明下成像过程可分作两步完成。xh说明：两次傅立叶变换频谱面上的光场分布与物的结构密切相关。原点附近分布着物的低频信息，即傅里叶低频分量，离原点越远，分布着物的较高的频率分量．即傅里叶高频分量。意义：不仅用傅立叶变换阐述了显微镜成像过程的机理，更重要的是首次引入了频谱的概念，启发人们可以通过改造频谱来改造信息。xh二、Abbe-Port实验用平行相干光束照明一张细丝网格，在成像透镜的后焦面上出现周期性网格状的傅里叶频谱，各傅里叶频谱分量在像平面上重新组合，复现网格的像。只呈现网格的垂直结构只呈现网格的水平结构均匀光场网格衬度反转像的性质发生变化结论：对像的垂直结构有贡献的是谱面上的水平方向的分量，而谱面上的垂直方向的分量只对像的水平结构有贡献；零频分量是一个直流分量，它只代表像的本底；挡住零频分量，在一定条件下可使像发生衬度反转；仅允许低频分量通过时，像的边缘锐度降低；仅允许高频分量通过时，像的边缘效应增强；采用选择型滤波器，可望完全改变像的性质。意义充分证明了阿贝成像理论的正确性．像的结构直接依赖于频谱的结构，只要改变频谱的构成，便可改变像的构成。五、空间滤波的傅里叶分析输出平面频谱面物平面物：一维光栅——朗奇（Ronchi）光栅Lxrectdxcombdaxrectxt11111频谱面上x2lfxdLcdacdLcdaucLcdaLdmLcdamcdaLTm1sinsin1sinsinsinsinsinxxxxx(1)滤波器是一个适当宽度的狭缝，只允许零级谱通过像的周期信息丢失，像的结构变成一矩形缝mdmLcdamcdaLTxxsinsinxxxLcdaLHTsinLxrectdaHTFxg313xx(2)狭缝加宽能允许零级和正、负一级频谱通过像场与物的周期相同，但由于高频信息的丢失，像的结构变成余弦振幅光栅dLcdaincdLcdacLcdaLHT1sin1sinsinsinxxxxxdxdacLxrectdadxjLxrectdacdxjLxrectdacLxrectdaHTFxg3333333132cossin212expsin2expsinxx(3)滤波面上放置双缝，只允许正、负二级频谱通过像场复振幅的周期是物的周期的一半，像的结构是余弦振幅光栅dLcdLcdacdaLHT2sin2sinsinxxxxdxdacLxrectdaHTFxg33134cos2sin2xx(4)在频谱面中心放置不透明的小圆屏a=d/2像场的复振幅分布仍为光栅结构，并且周期与物相同，但强度分布是均匀的，即实际上看不见条纹ad/2对应物体上亮的部分变暗，暗的部分变亮，实现了对比度反转xxxxLcdaLTHTsinLxrectdaLxrectdxcombdaxrectLcdaLFTFxg33331131sinxx例1在相干照明4f系统中，在物平面上放置一正弦光栅，其振幅透过率(1)在频谱面的中央设置一小圆屏挡住光栅的零级谱，求像的强度分布及可见度；(2)移动小圆屏，挡住光栅的+1级谱，像面的强度分布和可见度又如何?t(x1)t0t1cos(2x0x1)例2在相干照明4f系统中，在物平面上有两个图像，它们的中心在X轴上，距离坐标原点分别为a和-a，今在频谱面上放置一正弦光栅，其振幅透过率为试证明在像面中心可得到两个图像相加思考：如何在像面中心得到两个图像相减？12H(x,h)1cos(2ax)8.2空间滤波器空间l滤波系统频谱面物平面频谱功率谱11,yxfhx,F11,,yxfFFhxfxlx2fylh22,,hxhxFI功率谱分布特征：x2lfy2lfxh(1)频率特性(2)对称特性(3)方向特性中心的空间频率为零，由中心点向外其空间频率值越来越高，P2面上任一点的空间频率值与该点到中心点的距离成正比，其方向由中心点通过该点的矢径给出频谱面P2上的功率谱呈中心对称分布。若物图像存在线状构造，则其功率谱是沿着与此线状结构正交的直线方向r分布，且物图像中的线状构造越密集，则在P2沿r方向空间频谱分布延伸越远；反之亦然xh由在半圆中不同直径的32个环状PN结硅光二极管元件和另半圆呈辐射状分布的32个楔形PN结硅光二极管元件组成，据此可测出整个频谱面上各处的光强分布频谱分析器，又称为衍射图像采样器楔-环探测器针尖缺陷检查、掩模线宽测量、织物疵病以及纸张印刷质量的检查等。四、空间频率滤波系统１．4f系统傅里叶变换频谱面准直透镜傅里叶变换输出平面33332222133,,,,,yxhyxffyfxHfyfxFFyxgllll11,yxffyfxHll22,2．双透镜滤波系统准直透镜傅里叶变换和成像频谱面输出平面照明和傅里叶变换透镜输出平面傅里叶变换和成像频谱面3．单透镜滤波系统输出平面傅里叶变换和成像透镜频谱面双透镜和单透镜滤波系统优点：系统结构简单缺点：P2面上给出的物频谱都不是物函数的准确的傅立叶变换，而附带有二次相位因子空间滤波器的透过率函数一般是复函数根据透过率函数的性质，空间滤波器可分为振幅型、相位型和复数滤波器hxhxhx,exp,,jAH一、二元振幅滤波器只允许位于频谱面中心及其附近的低频分量通过，可以用来滤掉高频噪声（1）低通滤波器带针孔的不透明模板低通滤波器结构（2）高通滤波器带不透明小圆屏的透明模片阻挡低频分量而允许高频通过，以增强像的边缘，提高对模糊图像的识别能力或实现衬度反转，但由于能量损失较大，所以得到的结果一般较暗。高通滤波器结构（3）带通滤波器只允许特定区间的空间频谱通过，可以去除周期性信号中的随机噪音带通滤波器结构（4）方向滤波器阻挡(或允许)特定方向上的频谱分量通过，可以突出图像中的方向性特征三、位相型滤波器只改变频谱的位相结构，不改变振幅分布光学效率高，通常用真空镀膜的方法得到四、复数型滤波器可同时改变各频谱成分的振幅和位相分布，滤波函数是复函数二、振幅滤波器只改变傅里叶频谱的振幅分布，而不改变它的位相分布8.3空间滤波器应用举例一、策尼克相衬显微镜利用相位滤波器将物体的相位变化转换成可以观察到的光的强弱变化设位相物体的复振幅透过率为：一个普通的显微镜对上述物体所成的像，其强度可以写成用单位振幅的单色平面光波照明,透过的光场复振幅分布1jf(x1,y1)f(x1,y1)1时1111,exp,yxjyxt1111,1,yxjyxt1,1211yxjI像的强度和相移成线性关系滤波后像面复振幅分布像强度分布变相板本底光位相延迟/2，式中取正号，相位值大的部位光强也强，叫做正相衬；本底光位相延迟3/2，式中取负号时，相位值大的部位光强弱，叫做负相衬。可以用于观察抛光表面质量以及透明材料不均匀性检测等其它,10,,hxhxjHhxhxhxhx,,,,jjHF3333,,yxjjyxg3323333,21,1,yxyxjyxI例1：用相衬法观察一位相物体，若人眼可分辨的最小对比度V＝0.03，所用光波波长l＝600nm，试问：（1）当位相板上零级谱的振幅透过率t＝1时，可观察到的最小位相变化是多少？（2）当t＝0.01时，可观察到的最小位相变化又是多少？（3）若检测透明玻璃的不平度，在（1）的情况下，可检测的玻璃的不平度为多少？（设玻璃的折射率nG=1.52）假定成像缺陷是由于成像系统严重离焦引起二、补偿滤波器在频谱平面上放置适当的滤波器，使得滤波器的传递函数补偿原来系统传递函数的缺陷，则两者的乘积产生一个较为满意的频率响应，使照片的质量得到部分改善相应的传递函数arcircarh211adrrrJaarcircaRH00222211~rrrrrardrJraH2002221122hxrrar2高频损失严重，而且在某一中间频率区域，传递函数的符号发生反转吸收板用来衰减很强的低频峰值，以便提高像的对比，突出细节相移板使H的第一个负瓣相移，以纠正对比反转。8.4傅里叶变换透镜无拦阻通过截止频率截止频率、空间带宽积和视场1、截止频率D1/2某些空间频率DfDDu21fDDuulllx2sin1某些空间频率不能通过D1/2fDD21fDDllllx2sin1叶谱，各空间频率成分受到透镜孔径程度不同的阻拦；透镜后焦面上可以得到相应的空间频率成分的物体准确的傅里叶谱；虽然物可能有更高的空间频率成分，但因这些分量全部被透镜的有限孔径所阻拦，在后焦面上完全得不到物的傅里叶谱中的这些高频成分。结论：渐晕效应fDDlx21fDDfDDlxl2211fDDlx21例1设物函数中含有从低频到高频的各种结构信息，物被直径为d=2cm的圆孔所限制，将它放在直径D＝4cm、焦距f=50cm的透镜的前焦面上。今用波长l＝600nm的单色光垂直照射该物并测量透镜后焦面上的光强分布。问：（1）物函数中什么频率范围内的频谱可以通过测量得到准确值？（2）什么频率范围内的信息被完全阻止？2、傅立叶透镜的信息容量——空间带