平滑的频率域滤波器可以考虑三种滤波器:理想滤波器,巴特沃斯滤波器和高斯滤波器。这里我们介绍的是理想低通滤波器,简单的来说就是,截断傅里叶变换的高频成分,给定初始距离为D0,大于D0的为0,小于或等于D0的为1.如图所示:
![c074665b44cfb291e342e847be22825c.png](https://i-blog.csdnimg.cn/blog_migrate/b72827bfc8d936f0b1973cb2ea2e6d9f.png)
代码如下:
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)
![4d3966a35b74547b60a6d2b2e30ab7a6.png](https://i-blog.csdnimg.cn/blog_migrate/690bb624d61c19ab98e943eaac6bdb4f.jpeg)
而理想高通滤波器正好和低通相反,大于D0的为1,小于或等于D0的为0.如图所示:
![92a18b3276329135edc68f651942c649.png](https://i-blog.csdnimg.cn/blog_migrate/07e7178785f85f13fc42cb9b2b595900.png)
而图像显示为边界比较清晰
![1c27c1dce718c912284e1abf3207b2c3.png](https://i-blog.csdnimg.cn/blog_migrate/df59799991c6f25bc88da6f2977b81f5.jpeg)