理想低通滤波器

Halcom

1. function im5 = freqfilter_ideal_lp(im,D0)
   
2.     if ~isa(im,'double')
   
3.         im1 = double(im)/255;
   
4.     end
   
5. im2 = fft2(im1);      % 傅里叶变换
   
6. im3 = fftshift(im2);  % 中心化
   
7. [N1, N2] = size(im3);
   
8. for i = 1:N1
   
9.     for j = 2:N2
   
10.         if(im3(i,j) < D0)              % 进行理想低通滤波
    
11.             result(i,j) = 0;
    
12.         else
    
13.             result(i,j) = im3(i,j);
    
14.         end
    
15.     end
    
16. end
    
17. result = ifftshift(result); % 反中心化
    
18. im4 = ifft2(result);           % 反变换
    
19. im5 = im2uint8(real(im4));
    
20. 
    
21. end

复制代码