逆滤波器

Halcom

逆滤波器
resim = Inverse(im, 1.2, 30);       % 逆滤波

1. function resim = Inverse(ifbl, LEN, THETA)
   
2. % 逆滤波器
   
3. %函数输入:
   
4. %         ifbl:  输入的图像矩阵
   
5. %         THETA: 模糊旋转角
   
6. %         LEN:   模糊旋转长度，模糊的像素个数
   
7. %函数输出:
   
8. %         resim: 重构滤波图像
   
9. if ~isa(ifbl,'double')
   
10.     ifbl = double(ifbl)/255;
    
11. end
    
12. % 转化到频域
    
13. fbl = fft2(ifbl);  % 傅里叶变换
    
14. % 点扩展函数PSF
    
15. PSF = fspecial('motion',LEN,THETA);
    
16. % 转化PSF函数到期望的维数 光传递函数OTF
    
17. OTF = psf2otf(PSF, size(fbl));
    
18. % 检测是否存在0值，若为0，则置为0.000001
    
19. for i = 1:size(OTF, 1)
    
20.     for j = 1:size(OTF, 2)
    
21.         if OTF(i, j) == 0
    
22.             OTF(i, j) = 0.000001;
    
23.         end
    
24.     end
    
25. end
    
26. % 使用逆滤波器重构图像
    
27. fdebl = fbl./OTF;
    
28. % 使用逆傅里叶变换得到重构图像 IFFT
    
29. resim = ifft2(fdebl);

复制代码