实验四 图像复原
程序：
（1）卷积
I = imread('66.jpg');
PSF = fspecial('gaussian',7,10);
V = .0001;
IF1= imfilter(I,PSF);
BlurredNoisy = imnoise(IF1,'gaussian',0,V);
WT = zeros(size(I));
WT(5:end-4,5:end-4) = 1;
INITPSF = ones(size(PSF));
[J P] = deconvblind(BlurredNoisy,INITPSF);
figure(1);imshow(I,[]);
title('原始图像');
figure(2);
subplot(221);imshow(BlurredNoisy,[]); title('模糊加噪声');
subplot(222);imshow(PSF,[]); title('真实 PSF');
subplot(223);imshow(J,[]); title('模糊图像');
subplot(224);imshow(P,[]); title('复原 PSF');
（2）维纳
I = imread('66.jpg');
I = I(60+[1:256],222+[1:256],:);
figure; subplot(131);imshow(I,[]); title('原始图像');
LEN = 31;
THETA = 11;
PSF = fspecial('motion',LEN,THETA);
Blurred = imfilter(I,PSF,'circular','conv');
subplot(132);imshow(Blurred,[]);
title('模糊后的图像');
len=9;
theta=0;
PSF=fspecial('motion',len,theta);
wnr1=deconvwnr(Blurred,PSF);
subplot(133);imshow(wnr1,[]);title('维纳复原');
运行结果：
评论1