频域滤波之巴特沃斯低通滤波

最新推荐文章于 2024-01-30 17:01:01 发布

weixin_30725315

最新推荐文章于 2024-01-30 17:01:01 发布

阅读量1k

点赞数

文章标签： matlab

原文链接：http://www.cnblogs.com/natalie/p/4469361.html

版权

本文介绍了频域滤波中巴特沃斯低通滤波的原理和MATLAB实现。通过计算傅里叶变换，乘以滤波传递函数H(u,v)，再进行傅里叶反变换，实现图像的低通滤波。巴特沃斯滤波器的传递函数公式为H(u,v)=1/1+[D(u,v)/D0]^2n，而指数低通滤波器的传递函数为H(u,v)=exp-(D(u,v)/D0)^n。" 116633911,10552109,在Linux中编译与运行汇编程序的指南,"['Linux操作系统', '汇编语言', '编程工具', '软件开发']

摘要由CSDN通过智能技术生成

频率滤波的基本步骤是：

（1）计算需要增强图像的傅里叶变换。

（2）乘以一个滤波传递函数H(u,v)。

（3）将得到的结果进行傅里叶反变换，得到g(x,y).

理想低通滤波器是通过低频，截至高频。

巴特沃斯低通滤波器的传递函数为H(u,v)=1/1+[D(u,v)/D0]^2n

相关MATLAB程序如下:

I=imread('pic8.jpg');
I=rgb2gray(I);
% 叠加椒盐噪声
J1=imnoise(I,'salt & pepper');  
subplot(1,2,1)
imshow(J1),title('原始图像')
% 数据类型转换，MATLAB不支持图像的无符号整型的计算
f=double(J1);    
% 傅里叶变换
g=fft2(f);        
g=fftshift(g);     
[M,N]=size(g);
% 二阶巴特沃斯(Butterworth)低通滤波器
nn=2;           
d0=25;
m=fix(M/2); n=fix(N/2);
% 计算低通滤波器传递函数
for i=1:M
       for j=1:N
           d=sqrt((i-m)^2+(j-n)^2);
           h=1/(1+0.414*(d/d0)^(2*nn));  
           result(i,j)=h*g(i,j);
       end
end
result=