ch5_3_2.m

来自「这是用MATLAB来进行数字图象处理的程序。都是一些很常用的」· M 代码 · 共 28 行

%%%%% 5-3-2  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=5;
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
for j=1:N2
d=sqrt((i-n1)^2+(j-n2)^2);
if d==0
h=0;
else
h=1/(1+(d0/d)^(2*n));
end
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
X2=ifft2(result);
X3=uint8(real(X2));
figure,imshow(X3)