%例程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)