《利用MATLAB模拟光学简单空间滤波系统》由会员分享，可在线阅读，更多相关《利用MATLAB模拟光学简单空间滤波系统（8页珍藏版）》请在装配图网上搜索。
1、利用MATLAB模拟光学简单空间滤波系统摘要：阿贝成像原理是第一步在透镜的后焦面上得到物的空间频谱分布，第二步 成像则是合频的过程，实则是两次傅立叶变换。利用阿贝-波特实验装置和空间 滤波系统，可以对一幅光学图像进行光学信息处理。通过MATLAB环境编写程 序完成阿贝-波特实验和空间滤波的物理模型的构建并进行计算机模拟。关键词：MATLAB；阿贝成像原理；空间滤波；计算机模拟引言：早在1873年，阿贝(E. Abbe,18401905)在德国蔡司光学器械公司研究如 何提高显微镜的分辨本领问题时,就认识到相干成像的原理。空间滤波的主要目 的是通过有意识地改变像的频谱,使像实现所希望的变化。光学信
2、息处理是一个 更为广阔的领域，它是基于光学频谱分析，利用傅里叶综合技术，通过空域或频 域调制，借助空间滤波技术对光学信息进行处理的过程。阿贝提出的二次成像理 论和20世纪初的阿贝波特实验,已经为光学信息处理打下了一定的理论基础。在阿贝成像理论的教学中,单纯依靠数学推演来讲解,效果不好,特别是空间 频率、空间滤波等概念的形成有一定的困难。虽然可以通过演示阿贝- 波特实 验来加强教学效果,但由于在普通教室难以完成演示实验,在实验室又受仪器、场 地等方面的限制,实验现象不太理想。为此,我们设计出计算机模拟实验, 获得 较好的模拟效果。在学习了解了阿贝成像原理的基础上，我们可以通过MATLAB 完成对
3、阿贝-波特实验和空间滤波系统的计算机模拟，观察各种物体的空间频谱 分布，设计各种不同的空间滤波器。1.阿贝成像原理在相干平行光照明下,显微镜的物镜成像可以分成两步：第一步即分频过程，由入射光经过物平面P发生衍射在物镜的后焦面P上形成夫琅禾费衍射图12样;第二步称为合频或频谱综合过程，衍射图样作为新的子波源发出的球面波在 像平面上相干叠加成像。相干光的成像过程本质上是两次傅立叶变换,第一次是 将光场空间分布变成频率分布，第二次则是傅立叶逆变换,即将各频谱分量复合 为像。如下图所示，为阿贝成像原理图。P1wP阿 贝成像原理P3阿贝-波特实验是对阿贝成像理论最好的验证和演示。实验一般做法如下图 所示
4、，用平行相干光束照明一张细丝网格，在成像透镜后焦面上出现周期性网格 的傅里叶频谱，由这些傅里叶频谱分量的在组合，从而在像平面上再现网格得像。 若把各种遮挡物放在频谱面上，就能得到不同的像的频谱，从而得到由改变后的 频谱分量重新组合得到的对应的像。物体焦平面像阿贝波特实验2空间滤波的傅里叶分析2.1 4f系统滤波的傅里叶分析空间滤波就是利用滤波器（包括振幅滤波器、位相滤波器和复数滤波器） 在相干光学信息处理系统的空间频谱上面滤去一些空间频率成分，从而使像平面 的像按我们的意图来改变。以最典型的4f系统为例，用傅里叶分析的手段来讨 论空间滤波过程。如下图所示，图中：L是准直透镜；L和L为傅里叶变换
5、透 123镜，焦距均为f; P、P和P分别是物面、频谱面和像面，且P采用反演坐标。1233设光栅常数为d,缝宽为a,光栅沿x方向的宽度为L,则它的透过率为:1t( x )= rect (2 )* 丄 comb( 1) rect (2)1 a d d L采用单位振幅平面波垂直照明,P 面上的光场分布正比于物体的频谱2aLT( f )= 丁 xd縈. ammL sm c()sm cL(f - )d x dn=g= sin c(Lf ) + sinsin cL(f 1) + sin c(a) sin cL(f + ；)+d x x d d x d式中f =二。为了避免各级谱重叠，假定Ld,以便对每一
6、级谱实现单独处理。 九f2当在P面上放置不同的屏或孔径，作频域处理，就可以得到不同的输出像。22.2 4f系统构成的低通滤波器的傅里叶分析根据前面内容所述，在 P 面上选择适当宽度的狭缝，仅让零级谱通过或仅2让零级和正、负一级谱通过，限制高频成分通过，构成低通滤波器。仅让零级谱 成分通过时，紧靠狭缝后的透射频谱为:aLT( f )H( f )=sin c(Lf )x x d x式中 H( f )为相逢的透过率函数。 P 面上输出光场分布为:x3g( 3)= F -1T(f )H(f )= a rect()3x x dL仅让零级和正、负一级谱通过，透射频谱为:aL11T( f )H( f )=s
7、m c(Lf )+smsm cL(f - ) + smsm cL(f + )x x dxx dx dP3面上输出光场分布为:a x a x x xg(X严 W(fx)H(fx)= d 如()+血C(d)如(用P( j% 方)+ 如(Laxx ax 4./2a、2兀 x)+sinc() rect( 3)exp(- j2兀)= rect( 3 )1 + 2smc()cos 3 dLd d Ld d3频域低通滤波的MATLAB模拟3.1模拟方法T(u, v)根据傅里叶变换的性质，两个函数卷积的傅里叶变换等于傅里叶变换的乘 积。在频谱面上插入空间滤波器相当于频谱分布函数乘以空间滤波器滤波函数的 复振幅
8、透过率函数。空间滤波的光学处理器的模拟系统简图如上图所示，按图通 过计算机模拟仿真可以完成空间滤波实验。3.2二维光栅的频谱 将二维光栅作为物,则可在傅立叶面上观测到如图所示的频谱分布。在MATLAB 中输入以下指令： x=ones(150,150); %创建矩阵 x(1:9:150,:)=0; %得到 1 维光栅 y=x.*(x); %得到 2 维光栅 m=fft2(y,200,200); %傅立叶变换 n=abs(fftshift(m); %变换象限并取模 imshow(0.01*n); %以一定比例显示图像File Edit View. Insert Tools- Desktop Win
9、dow Help33低通滤波模拟结果在计算机模拟中，用一幅图像代替物体。对这幅图进行傅立叶变换得到相 应的频谱分布。这一步骤相当于实验中透镜所起的傅立叶变换作用。下图中所示 为原图像及其频谱图分布。将设计的低通滤波器与经过傅立叶变换过的频谱相乘。这一步相当于实验 中在频谱面上设置低通滤波器进行滤波。经过低通滤波，滤掉了物的高频信息,处理图像显示了物的低频信息即光字。由于经过滤波后的图像能量有所损失导致输出图像比原图像模糊。下图为经过低通滤波后的图像及其频谱图分布。滤波程序如下：fO=imread(3.bmp);F=fft2(f0);M,N=size(F);subplot(l,2,l);imsh
10、ow(fO);title(原图像)； size=256;lenafft=fftshift(fft2(f0);subplot(1,2,2);i=1:1:size;j=1:1:size;mesh(i,j,abs(lenafft(i,j);axis(1 size 1 size 0 40000);title(频谱图);r=88;for s=1for i=1:sizefor j=1:sizeif sqrt(i-128)A2+(j-128)A2)r lenafft(i,j)=0;endendend figure(s+1); subplot(1,2,2);i=1:1:size;j=1:1:size; mes
11、h(i,j,abs(lenafft(i,j); axis(1 size 1 size 0 40000); f0=ifft2(lenafft);figure(s+1);subplot(1,2,1); imshow(abs(f0);endclear结论：本文在阿贝成像原理的基础上介绍了阿贝波特实验，并对其实验结果做 出简要的分析和总结。试验充分的证明了阿贝成像理论的正确性，为后面的空间 滤波系统奠定了基础。通过对4f系统的成像过程的分析，并利用透镜的傅里叶 变换性质，把透镜作为一个频谱分析仪，利用空间滤波的方式改变物的频谱结构， 从而改变像。其实质是改变滤波器的振幅透过率函数，从而达到改变像结构的
12、目 的。通过MATLAB模拟了二维光栅的频谱和低通滤波的效果。借助Matlab构建模 型模拟光学频谱分析系统进行空间滤波实验模拟,能够显示复杂的物理现象,使 抽象的问题形像化,从而加深对空间频率、频谱、空间滤波和等的理解。参考文献：1 苏显渝,李继陶.信息光学.北京:科学出版社,1999.2 何钰西南交大应用物理系，物理与工程 V ol. 16 No. 2 20063 蒲利春等，大学应用物理实验，科学出版社，2011.4 许录平.数字图像处理.北京:科学出版社,2007.5 彭芳麟，计算物理基础，高等教育出版社，2010.6 谢嘉宁,赵建林,光学空间滤波过程的计算机仿真，光子学报 Vol.31NO.7.2002.
利用MATLAB模拟光学简单空间滤波系统
