%程序L5_3.m:同態濾波器%
   clear
   [image_0,map]=bmpread('L5_3.bmp');            %讀取圖像 
   image_1=log(image_0+1);                       %取自然對數
   image_2=fft2(image_1);                        %做二維傅立葉變換
  
   %產生高通Butterworth濾波器%
   n=3;
   D0=0.05*pi;
   rh=0.8;
   rl=0.5;
   [row,col]=size(image_2);
   for k=1:1:row
     for l=1:1:col
       D1(k,l)=sqrt((k^2+l^2));
       H(k,l)=rl+(rh/(1+(D0/D1(k,l))^(2*n)));
     end
   end

   image_3=(image_2.*H);                         %輸入圖像通過濾波器
   image_4=ifft2(image_3);                       %做二維反傅立葉變換
   image_5=(exp(image_4)-1);                     %取指數函數

   figure(1)               
   imshow(image_0,map)                           %顯示原圖像
   figure(2)
   imshow(real(image_5),map)                     %顯示處理后的圖像