基于频率域的平滑滤波及其应用-数字图像处理论文

《基于频率域的平滑滤波及其应用-数字图像处理论文》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、数字图像处理论文频率域低通滤波及其应用第六组：201583130145赵明翔摘要：数字图像处理的主要目的是改善图像视觉效果，以便于人们对其更进一步的理解，而图像增强技术是采用某种技术手段改变图像的视觉效果，其中一种技术就是频率域增强技术，本论文重点介绍了用MATLAB对图像分别进行理想低通滤波、布特沃斯低通滤波和高斯低通滤波。关键词：MATLAB滤波理想低通滤波布特沃斯低通滤波高斯低通滤波一．引言滤波是将信号中特定波段频率滤波的操作，是抑制和防止干扰的一项重要措施。根据观察某一随机过程的结果，对另一与之

2、有关的随机过程进行估计的概率理论与方法。而低通滤波则是让低频信号能正常通过，而超过设定临界值的高频信号则会被阻隔、减弱，但是阻隔或减弱的幅度会根据不同的频率以及不同的滤波程序而改变。二．频率域滤波基础频率域滤波有修改一幅图像的傅里叶变换然后计算其反变换得到处理后的结果组成。由此，若给定一副大小M×N的数字图像f(x,y),则基本滤波公式有如下形式-1其中，-1是IDFT，是输入图像的DFT，是滤波函数34数字图像处理论文是滤波后的（输出）图像。函数F，H和g是大小与输出图像相同的阵列。乘积是阵列相乘而得

3、到的，滤波函数通过修改输入的图像数据的变换得到处理后的输出。变换中的低频与图像中缓慢变化的灰度分量有关，而高频由灰度的尖锐过程造成。因此我们所利用的让高频衰减而让低频通过的滤波器（也就是我们下面即将讨论的低通滤波器）将会使图像变得模糊。三．常用的低通滤波器在这门课程中我们主要学习的低通滤波器有三种，分别是：理想滤波器、布特沃斯滤波器和高斯滤波器。这三种滤波器涵盖了从非常尖锐到非常平滑的滤波范围。3.1理想滤波器在以原点为圆心，以D0为半径的圆内，无衰减地通过所有频率，而在该圆外“切断”所有频率的二维低通

4、滤波器，称为理想低通滤波器（ILPF）；它由下面的函数确定：1，00，>0其中，0是一个正常数，是频率域中点与频率矩形中心的距离，即221/2（*）其中，与以前一样，P和Q是和34数字图像处理论文填充后的尺寸。图4.40(a)显示了的透视图，图4.40(b)显示了以图像显示的滤波器。理想这一名称表明在半径0的圆内，所有频率无衰减地通过，而在此圆之外的所有频率则完全被衰减。理想低通滤波器是关于原点径向对称的，这意味着该滤波器完全由一个径向横截面来定义，如图4.40(c)所示。将该横截面旋转360o得到二维

5、滤波器。3.2布特沃斯低通滤波器截止频率位于距原点0处的n阶布特沃斯低通滤波器（BLPF）的传递函数定义为其中，由（*）式给出。图4.44显示了该BLPF函数的透视图、图像显示和径向剖面图。3.3高斯低通滤波器34数字图像处理论文二维高斯滤波器的形式由下式给出：其中，是距频率矩形中心的距离。是关于中心的扩展度的度量。通过令=，我们可以使用其他滤波器的表示法来表示该滤波器：其中，是截止频率。当时，GLPF下降到其最大值的0.607处。四．Matlab实现:4.1理想低通滤波matlab实现:函数程序：理想

6、低通滤波器（以D0=5为例）•I1=imread(‘D:Matlabproject低通、高通滤波实验原图.jpg’);•imshow(I1);title('原图‘);•s=fftshift(fft2(I1));•[M,N]=size(s);%分别返回s的行数到M中，列数到N中•n1=floor(M/2);%对M/2进行取整•n2=floor(N/2);%对N/2进行取整•%ILPF滤波，d0=5，15，30（程序中以d0=5为例）•d0=5;%初始化d0•fori=1:M•forj=1:N•d=s

7、qrt((i-n1)^2+(j-n2)^2);%点（i,j）到傅立叶变换中心的距离•ifd<=d0%点（i,j）在通带内的情况•h=1;%通带变换函数•else%点（i,j）在阻带内的情况•h=0;%阻带变换函数•end•s(i,j)=h*s(i,j);%ILPF滤波后的频域表示•end•end•s=ifftshift(s);•%对s进行反FFT移动%对s进行二维反离散的Fourier34数字图像处理论文变换后，取复数的实部转化为无符号8位整数•s=uint8(real(ifft2(s)));•figu

8、re;%创建图形图像对象•Imshow(s);%显示ILPF滤波后的图像•title(‘ILPF滤波(d0=5)‘)4.2布特沃斯低通滤波matlab实现函数程序:巴特沃斯低通滤波：（以D0=5为例）•I1=imread('D:Matlabproject低通、高通滤波实验原图.jpg');•imshow(I1);title('原图');•s=fftshift(fft2(I1));•[M,N]=size(s);%分别返回s的行数到M中，