频域低通滤波器（数字图像处理与机器视觉P210-P220）

一、理想低通滤波
%%理想低通滤波
function out=imidealflpf(I,freq)
% imidealflpf 函数 构造理想的频域低通滤波器
% I参数 输入的灰度图像
% freq 参数 低通滤波器的截止频率
% 返回值 out -指定的理想低通滤波器
[M,N]=size(I);
out = ones(M,N);
for i=1:M
for j=1:N
if (sqrt(((i-M/2)^2+(j-N/2)^2))>freq)
out(i,j)=0;
end
end
end
end

二、高斯低通滤波
%% 高斯低通滤波
function out=imgaussflpf(I,sigma)
%imgaussflpf函数 构造频域高斯低通滤波器
%I参数 输入的灰度图像
%sigma参数 高斯函数的sigma参数

[M,N]=size(I);
out=ones(M,N);
for i=1:M
for j=1:N
out(i,j)=exp((-(i-M/2)^2-(j-N/2)^2)/(2*sigma^2));
end
end
end

三、对灰度图像进行频域滤波
function out=imfreqfilt(I,ff)
%imfreqfilt 函数 对灰度图像进行频域滤波
%参数I 输入的空域图像
%参数ff 应用的与原图等大的频域滤镜

if (ndims(I)==3)&&(size(I,3)==3)%RGB图像
I=rgb2gray(I);
end

if (size(I)~=size(ff))
msg1=sprintf(‘%s:滤镜与原图像不等大，检查输入’,mfilename);
msg2=sprintf(‘%s:滤波操作已经取消’,mfilename);
eid=sprintf(‘Images:%s:ImageSizeNotEquqal’,mfilename);
error(eid,’%s %s’,msg1,msg2);
end

%快速傅里叶变换
f=fft2(I);

%移动原点
s=fftshift(f);

%应用滤镜及反变换
out=s.*ff;%对应元素相乘实现频域滤波
out=ifftshift(out);
out=ifft2(out);

%求模值
out=abs(out);

%归一化以便显示
out=out/max(out(:));
end