平滑的频率域滤波器可以考虑三种滤波器:理想滤波器,巴特沃斯滤波器和高斯滤波器。这里我们介绍的是理想低通滤波器,简单的来说就是,截断傅里叶变换的高频成分,给定初始距离为D0,大于D0的为0,小于或等于D0的为1.如图所示:
代码如下:
clc
clear
A=imread('Fig0441(a)(characters_test_pattern).tif');
A=double(A)/255;
[m,n]=size(A);
for i=1:size(A,1)
for j=1:size(A,2)
B(i,j)=(-1)^(i+j)*A(i,j);
end
end
FB=fft2(B); %傅里叶变换
for i=1:size(A,1)
for j=1:size(A,2)
if sqrt((i-m/2)^2+(j-n/2)^2)>50 %50相当于D0
H(i,j)=0;
else
H(i,j)=1;
end
end
end
FB=FB.*H; %这里是H理想低通滤波器
C=abs(real(ifft2(FB)));
subplot(2,1,1)
imshow(A)
subplot(2,1,2)
imshow(C)
而理想高通滤波器正好和低通相反,大于D0的为1,小于或等于D0的为0.如图所示:
而图像显示为边界比较清晰