demo4_Lucy_Rich_deblur.m

% Example 4 - Restore fuzzy, noisy image
% --------------------------------------

% start with clean slate
close all, clear all, clc

% Original image
I = imread('board.tif');
I = I(50+[1:256],2+[1:256],:);
figure, imshow(I), title('Original Image')
set(1,'pos',[9 71 330 314])
pause %blur this image

% PSF blur + noise
PSF = fspecial('gaussian',5,5);
Blurred = imfilter(I,PSF,'symmetric','conv');
BlurredNoisy = imnoise(Blurred,'gaussian',0,0.002);
figure, imshow(BlurredNoisy), title('Blurred & Noisy')
set(2,'pos',[349 70 330 314])
pause %try to restore

% Restore image
luc1 = deconvlucy(BlurredNoisy,PSF,5);
figure, imshow(luc1), title('Restored Image, NUMIT = 5')
set(3,'pos',[688 72 330 314])

% Copyright 2003-2010 RBemis The MathWorks, Inc.