巴特沃斯(butterworth)低通滤波器对受噪声干扰的图像进行平滑处理

《巴特沃斯(butterworth)低通滤波器对受噪声干扰的图像进行平滑处理》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、%频域增强（巴特沃斯原型）%利用巴特沃斯（Butterworth）低通滤波器对受噪声干扰的图像进行平滑处理clc;clear;J1=imread('lena.gif');subplot(3,2,1);imshow(J1);title('原图');f=double(J1);g=fft2(f);%傅立叶变换g=fftshift(g);%转换数据矩阵subplot(3,2,2);x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);z=log(abs(g));%取幅度mesh(z);%以三维坐标显示该图像频谱图title('原图频谱');

2、[M,N]=size(g);nn=2;%二阶巴特沃斯(Butterworth)低通滤波器d0=20;m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt((i-m)^2+(j-n)^2);h=1/(1+0.414*(d/d0)^(2*nn));%计算低通滤波器传递函数result(i,j)=h*g(i,j);endendsubplot(3,2,4);x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);z=log(abs(result));%取幅度mesh(z);%以三维坐标显示该图像频谱图tit

3、le('低通滤波后的频谱');result=ifftshift(result);J2=ifft2(result);J3=uint8(abs(J2));subplot(3,2,3);imshow(J3);title('低通滤波后的图像');%利用巴特沃斯（Butterworth）高通滤波器对受噪声干扰的图像进行平滑处理nn=2;%二阶巴特沃斯(Butterworth)高通滤波器d0=5;m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt((i-m)^2+(j-n)^2);if(d==0)h=0;elseh=1/(1+0.4

4、14*(d0/d)^(2*nn));%计算传递函数endresult(i,j)=h*g(i,j);endendsubplot(3,2,6);x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);z=log(abs(result));%取幅度mesh(z);%以三维坐标显示该图像频谱图title('高通滤波后的频谱');result=ifftshift(result);J2=ifft2(result);J3=uint8(abs(J2));subplot(3,2,5);imshow(J3);title('高通滤波后的图像');