Wiener滤波器

Halcom

Wiener滤波器:
SNR = 0.0001;   % 信噪比
resim = Wiener(im, 1.2, 30,SNR);       % 逆滤波

1. function resim = Wiener(ifbl, LEN, THETA, SNR)
   
2. % Wiener滤波器
   
3. %函数输入:
   
4. %         ifbl:  输入的图像矩阵
   
5. %         LEN:   模糊旋转长度，模糊的像素个数
   
6. %         THETA: 模糊旋转角
   
7. %         SNR：信噪比
   
8. %函数输出:
   
9. %         resim: 重构滤波图像
   
10. 
    
11. ifbl = medfilt2(abs(ifbl));%中值滤波
    
12. fbl = fft2(ifbl); % 转化到频域
    
13. % 点扩展函数PSF
    
14. PSF = fspecial('motion',LEN,THETA);
    
15. % 转化PSF函数到期望的维数 光传递函数OTF
    
16. OTF = psf2otf(PSF,size(fbl));
    
17. OTFC = conj(OTF);   % 共轭操作
    
18. modOTF = OTF.*OTFC; % 平方
    
19. % 检测是否存在0值，若为0，则置为0.000001
    
20. for i = 1:size(OTF, 1)
    
21.     for j = 1:size(OTF, 2)
    
22.         if OTF(i, j) == 0
    
23.             OTF(i, j) = 0.000001;
    
24.         end
    
25.     end
    
26. end
    
27. % 使用Wiener滤波器公式，去模糊图像
    
28. debl = ((modOTF./(modOTF+SNR))./(OTF)).*fbl;
    
29. lbed = ifft2(debl);    % 傅里叶反变换
    
30. resim = lbed;          % 重构滤波图像

复制代码