利用MATLAB模拟光学简单空间滤波系统 2

1、利用MATLAB模拟光学简单空间滤波系统摘要：阿贝成像原理是第一步在透镜的后焦面上得到物的空间频谱分布，第二步 成像则是合频的过程，实则是两次傅立叶变换。利用阿贝-波特实验装置和空间 滤波系统，可以对一幅光学图像进行光学信息处理。通过MATLAB环境编写程 序完成阿贝-波特实验和空间滤波的物理模型的构建并进行计算机模拟。关键词：MATLAB；阿贝成像原理；空间滤波；计算机模拟引言：早在1873年，阿贝(E. Abbe,18401905)在德国蔡司光学器械公司研究如 何提高显微镜的分辨本领问题时,就认识到相干成像的原理。空间滤波的主要目 的是通过有意识地改变像的频谱，使像实现所希望的变化。光学信

2、息处理是一个 更为广阔的领域，它是基于光学频谱分析，利用傅里叶综合技术，通过空域或频 域调制，借助空间滤波技术对光学信息进行处理的过程。阿贝提出的二次成像理 论和20世纪初的阿贝一波特实验，已经为光学信息处理打下了一定的理论基础。在阿贝成像理论的教学中，单纯依靠数学推演来讲解，效果不好，特别是空间 频率、空间滤波等概念的形成有一定的困难。虽然可以通过演示阿贝-波特实验 来加强教学效果，但由于在普通教室难以完成演示实验，在实验室又受仪器、场地 等方面的限制，实验现象不太理想。为此，我们设计出计算机模拟实验，获得较好 的模拟效果。在学习了解了阿贝成像原理的基础上，我们可以通过MATLAB完 成对阿

3、贝-波特实验和空间滤波系统的计算机模拟，观察各种物体的空间频谱分 布，设计各种不同的空间滤波器。1.阿贝成像原理在相干平行光照明下，显微镜的物镜成像可以分成两步：第一步即分频过 程，由入射光经过物平面P1发生衍射在物镜的后焦面P上形成夫琅禾费衍射图 样;第二步称为合频或频谱综合过程，衍射图样作为新的子波源发出的球面波在 像平面上相十叠加成像。相十光的成像过程本质上是两次傅立叶变换，第一次是 将光场空间分布变成频率分布，第二次则是傅立叶逆变换，即将各频谱分量复合 为像。如下图所示，为阿贝成像原理图。贝成像原理阿贝-波特实验是对阿贝成像理论最好的验证和演示。实验一般做法如下图 所示，用平行相十光束

4、照明一张细丝网格，在成像透镜后焦面上出现周期性网格 的傅里叶频谱，由这些傅里叶频谱分量的在组合，从而在像平面上再现网格得像。 若把各种遮挡物放在频谱面上，就能得到不同的像的频谱，从而得到由改变后的 频谱分量重新组合得到的对应的像。物体焦平面像阿贝波特实验2.空间滤波的傅里叶分析2.1 4f系统滤波的傅里叶分析空间滤波就是利用滤波器（包括振幅滤波器、位相滤波器和复数滤波器） 在相十光学信息处理系统的空间频谱上面滤去一些空间频率成分，从而使像平面 的像按我们的意图来改变。以最典型的4f系统为例，用傅里叶分析的手段来讨 论空间滤波过程。如下图所示，图中：匕是准直透镜；和L3为傅里叶变换透 镜，焦距均

5、为f； P、P和P分别是物面、频谱面和像面，且P采用反演坐标。1233设光栅常数为d,缝宽为a，光栅沿X1方向的宽度为L，则它的透过率为:X1XXt( X )= rect ( )* comb(t ) rect ()1 addL采用单位振幅平面波垂直照明，P面上的光场分布正比于物体的频谱T( f )= aL 党 sin c(岑)sin cL( f - m)X ddX dn=-s=d sin c(Lf ) + sin c(d) sin cL( f-分申近 c(d) sin cL( f + ：)+式中f = f。为了避免各级谱重叠，假定Ld，以便对每一级谱实现单独处理。当在P面上放置不同的屏或孔径，

6、作频域处理，就可以得到不同的输出像。 22.2 4f系统构成的低通滤波器的傅里叶分析根据前面内容所述，在P面上选择适当宽度的狭缝，仅让零级谱通过或仅2让零级和正、负一级谱通过，限制高频成分通过，构成低通滤波器。仅让零级谱 成分通过时，紧靠狭缝后的透射频谱为：T( fx )H( f尸 aL sin c(Lf?式中H( f)为相逢的透过率函数。尸3面上输出光场分布为：g(x )= F -1T( f )H( f )= a rect (：3)3x x d L仅让零级和正、负一级谱通过，透射频谱为: TOC o 1-5 h z aL 一 一 . 一 . a1a1T( f )H( f )= sin c(L

7、f ) + sin c( ) sin cL( f 一 ) + sin c( ) sin cL( f + ) dX d dX dP3面上输出光场分布为: HYPERLINK l bookmark13 o Current Document axaxxxg(X )= F-1T( f )H( f )- rect(亍)+sinc( ) rect(3)exp( j2兀寸)+ rect(矽3x x dLdLdL、.,a、X a X2a、2兀X n)+sin c() rect ( )exp( j2兀)= rect ( )1 + 2sin c()cos d Ld d Ld d3.频域低通滤波的MATLAB模拟3

8、.1模拟方法，)T(u， v)根据傅里叶变换的性质，两个函数卷积的傅里叶变换等于傅里叶变换的乘 积。在频谱面上插入空间滤波器相当于频谱分布函数乘以空间滤波器滤波函数的 复振幅透过率函数。空间滤波的光学处理器的模拟系统简图如上图所示，按图通 过计算机模拟仿真可以完成空间滤波实验。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)

9、; %变换象限并取模 imshow(0.01*n); %以一定比例显示图像3.3低通滤波模拟结果在计算机模拟中，用一幅图像代替物体。对这幅图进行傅立叶变换得到相 应的频谱分布。这一步骤相当于实验中透镜所起的傅立叶变换作用。下图中所示 为原图像及其频谱图分布。将设计的低通滤波器与经过傅立叶变换过的频谱相乘。这一步相当于实验 中在频谱面上设置低通滤波器进行滤波。经过低通滤波，滤掉了物的高频信息, 处理图像显示了物的低频信息即光字。由于经过滤波后的图像能量有所损失导致 输出图像比原图像模糊。下图为经过低通滤波后的图像及其频谱图分布。滤波程序如下：f0=imread(3.bmp);F=fft2(f0)

10、;M,N=size(F);subplot(1,2,1);imshow(f0);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;endendendfigure(s+1);subplot(1,2

11、,2);i=1:1:size;j=1:1:size;mesh(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系统的成像过程的分析，并利用透镜的傅里叶 变换性质，把透镜作为一个频谱分析仪，利用空间滤波的方式改变物的频谱结构， 从而改变像。其实质是改变滤波器的振幅透过率函数，从而达到改变像结构的目 的。通过MATLAB模拟了二维光栅的频谱和低通滤波的效果。借助Matlab构建模 型模拟光学频谱分析系统进行空间滤波实验模拟，能够显示复杂的物理现象，使 抽象的问题形像化,从而加深对空间频率、频谱、空间滤波和等的理解。参考文献：