clc
clear all
close all
I=imread('1.png');
subplot(211);
imshow(I);
title('原图');
I=rgb2gray(I);
subplot(212);
imshow(I);
title('灰度图像');
PSF=fspecial('gaussian',7,10);
V=.0001;
BlurredNoisy=imnoise(imfilter(I,PSF),'gaussian',0,V);
BlurredNoisy=double(BlurredNoisy);
WT=zeros(size(I));
WT(5:end-4,5:end-4)=1;
J1=deconvlucy(BlurredNoisy,PSF);
J2=deconvlucy(BlurredNoisy,PSF,20,sqrt(V));
J3=deconvlucy(BlurredNoisy,PSF,20,sqrt(V),WT);
figure(2)
subplot(211);
imshow(mat2gray(BlurredNoisy));
title('A=Blurred and Noisy');
subplot(212);
imshow(mat2gray(J1));
title('deconvlucy(A,PSF)');
figure(3)
subplot(211);imshow(mat2gray(J2));
title('deconvlucy(A,PSF,NI,DP)');
subplot(212);imshow(mat2gray(J3));
title('deconvlucy(A,PSF,NI,DP,WT)');