ch5_3_1.m

来自「Matlab7.x图像处理的光盘内容 人民邮电出版社的最新」· M 代码 · 共 27 行

%%%%%  5-3-1 Buterworth低通滤波器实现例程如下：
clear all;
I1=imread('E:\MATLAB7\work\5-2\pinyulvboyuantu.bmp');
figure,imshow(I1,[]);
f=double(I1);
g=fft2(f);
%傅立叶变换
g=fftshift(g);
%转换数据矩阵
[N1,N2]=size(g);
n=2;
d0=50;
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
for j=1:N2
d=sqrt((i-n1)^2+(j-n2)^2);
h=1/(1+0.414*(d/d0)^(2*n));
%计算Butterworth低通转换函数
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
X2=ifft2(result);
X3=uint8(real(X2));
figure,imshow(X3)