matlab中图像双边滤波_MATLAB 理想滤波器-图像处理

平滑的频率域滤波器可以考虑三种滤波器：理想滤波器，巴特沃斯滤波器和高斯滤波器。这里我们介绍的是理想低通滤波器，简单的来说就是，截断傅里叶变换的高频成分，给定初始距离为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.如图所示：

而图像显示为边界比较清晰