Lucy_Richardson滤波器

Halcom

Lucy_Richardson滤波器
im1 = Lucy_Richardson(im, 5, 10, 50);

1. function resim = Lucy_Richardson(ifbl, LEN, THETA, iterations)
   
2. % Lucy_Richardson滤波器
   
3. %函数输入:
   
4. %         ifbl:  输入的图像矩阵
   
5. %         LEN:   模糊旋转长度，模糊的像素个数
   
6. %         THETA: 模糊旋转角
   
7. %         iterations: 迭代次数
   
8. %函数输出:
   
9. %         resim: 重构滤波图像
   
10. if ~isa(ifbl,'double')
    
11.     ifbl = double(ifbl);
    
12. end
    
13. 
    
14. ifbl = medfilt2(abs(ifbl)); %中值滤波
    
15. est = ifbl; % 初始化模糊图像赋值操作
    
16. % 点扩展函数PSF
    
17. PSF = fspecial('motion',LEN,THETA);
    
18. % 转化PSF函数到期望的维数 光传递函数OTF
    
19. OTF = psf2otf(PSF,size(ifbl));
    
20. i = 1;
    
21. while i<=iterations
    
22.     fest = fft2(est);    % 转化到频域
    
23.     fblest = OTF.*fest;  % 频率相乘
    
24.     ifblest = ifft2(fblest); % 傅里叶反变换，转换到时域
    
25.     % 计算模糊前的图像 与 去模糊后的图像 的比值
    
26.     iratio = ifbl./(ifblest+eps);
    
27.     firatio = fft2(iratio);   % 转化到频域
    
28.     corrvec = OTF .* firatio; % 计算相关性向量
    
29.     icorrvec = ifft2(corrvec);% 傅里叶反变换，转换到时域
    
30.     % 计算估计的去模糊图像矩阵
    
31.     aftercorr = icorrvec.*est;
    
32.     est = aftercorr;
    
33.     i = i+1;
    
34. end
    
35. resim = abs(est); % 重构滤波图像

复制代码