psnr.zip_The Given

I = imread('tissue.png'); I = I(125+(1:256),1:256,:); figure;imshow(I);title('Original Image'); text(size(I,2),size(I,1)+15, ... 'Image courtesy of Alan Partin, Johns Hopkins University', ... 'FontSize',7,'HorizontalAlignment','right'); PSF = fspecial('gaussian',11,5); Blurred = imfilter(I,PSF,'conv'); figure;imshow(Blurred); title('Blurred'); V = .02; BlurredNoisy = imnoise(Blurred,'gaussian',0,V); figure;imshow(BlurredNoisy); title('Blurred & Noisy'); NP = V*numel(I); % noise power [reg1, LAGRA] = deconvreg(BlurredNoisy,PSF,NP); figure,imshow(reg1),title('Restored with NP'); reg2 = deconvreg(BlurredNoisy,PSF,NP*1.3); figure;imshow(reg2); title('Restored with larger NP'); reg3 = deconvreg(BlurredNoisy,PSF,NP/1.3); figure;imshow(reg3); title('Restored with smaller NP'); Edged = edgetaper(BlurredNoisy,PSF); reg4 = deconvreg(Edged,PSF,NP/1.3); figure;imshow(reg4); title('Edgetaper effect');