4-28.m

matlab数字信号与应用源代码4

%例程4－28实现Butterworth高通滤波器

D = imread('rice.png');
figure, imshow(D),
D=double(D);
F=fft2(D);        %傅立叶变换
F=fftshift(F);    %转换数据矩阵
[N1,N2]=size(F);
n=2;
d0=4;
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+(d0/d)^(2*n));     %计算高通转换函数
       FD(i,j)=h*F(i,j);        %高通滤波
end
end
FD=ifftshift(FD);
FD=ifft2(FD);
FD=real(FD);
FD=100*FD;                        %拉升图像灰度范围
FD=uint8(FD);
figure,imshow(FD)