5.6 频域滤波(附matlab滤波器代码)

5.6.1 低通滤波

低通滤波器的功能是让低频率通过而滤掉或衰减高频，其作用是过滤掉包含在高频中的噪声。所以低通滤波的效果是图像去噪声平滑增强，但同时也抑制了图像的边界，造成图像不同程度上的模糊。对于大小为M×N的图像，频率点（u,v）与频域中心的距离为D（u，v），其表达式为：

                 D(u,v)=[(u-M/2)^2+(v-N/2)^2]^(1/2)

理想低通滤波器

产生公式为：

其中D0为理想低通滤波器的截止频率。理想低通滤波器的形状如图所示。

理想低通滤波器在半径为D0的范围内，所以频率都可以没有衰减地通过滤波器，该半径之外的所以频率都完全被衰减掉。理想低通滤波器具有平滑图像的作用，但是有很严重的振铃现象。

close all 
clear all
clc

I=imread('coins.png');

I=im2double(I); 

M=2*size(I,1);                 %滤波器的行数

N=2*size(I,2);                 %滤波器的列数

u=-M/2:(M/2-1);                

v=-N/2:(N/2-1);

[U,V]=meshgrid(u, v);

D=sqrt(U.^2+V.^2);

D0=80;                         %截至频率

H=double(D<=D0);               %理想低通滤波器

J=fftshift(fft2(I,size(H, 1),size(H,2)));   %时域转换到频域

K=J.*H;                        %傅里叶反变换

L=ifft2(ifftshift(K)); 

L=L(1:size(I,1),1:size(I, 2));

figure;

subplot(121);imshow(I);

subplot(122);imshow(L);

在程序中，设计了理想低通滤波器，截止频率为80。通过二维离散傅里叶变换将图像转换为频域，频域图像乘以滤波器的系数，然后进行二维傅里叶反变换转换到时与图像。程序运行后，结果如下图所示。通过低通滤波器去掉了图像中的高频部分，图像的边缘变得模糊。

巴特沃斯低通滤波器

产生公式为

其中D0为巴特沃斯低通滤波器的截止频率，参数n为巴特沃斯滤波器的阶数，n越大则滤波器的形状越陡峭。

close all
clear all
clc

I=imread('liftingbody.png');   
 
I=im2double(I);
 
M=2*size(I,1);                               %滤波器的行数

N=2*size(I,2);                               %滤波器的列数

u=-M/2:(M/2-1);

v=-N/2:(N/2-1);

[U,V]=meshgrid(u,v);

D=sqrt(U.^2+V.^2);

D0=50;                                        %截止频率

n=6;                                          %滤波器的阶数

H=1./(1+D./D0).^(2*n);                        %设计巴特沃斯滤波器

J=fftshift(fft2(I,size(H,1),size(H,2)));      %转换到频域

K=J.*H;                                       %滤波处理

L