lowpass.m

清华大学出版社出版的《数字图像处理实训教程》一书以提高动手能力为目标

Img=imread('3-22.jpg');
imshow(Img);
Imgf=double(Img);     % 数据类型转换，MATLAB不支持图像的无符号整型的计算
g=fft2(Imgf);        % 傅立叶变换
g=fftshift(g);     % 转换数据矩阵,使低频分量移到中间，可查看fftshift函数的帮助
[M,N]=size(g);
nn=2;           % 二阶巴特沃斯(Butterworth)低通滤波器
d0=50;
m=fix(M/2); n=fix(N/2);
for i=1:M
       for j=1:N
           d=sqrt((i-m)^2+(j-n)^2);
           h=1/(1+0.814*(d/d0)^(2*nn));  % 计算低通滤波器传递函数
           result(i,j)=h*g(i,j);
       end
end
result=ifftshift(result);
J2=ifft2(result);
J3=uint8(real(J2));
figure,imshow(J3);                      % 显示滤波处理后的图像