%% 学习目标: 约束最小二乘算法对彩色模糊和加性噪声图像进行恢复效果实验
I = imread('01.bmp');
I = I(0+[1:256],0+[1:256],:);
figure;imshow(I);title('原始图像');
%% 生成PSF点扩展函数
PSF = fspecial('gaussian',11,5);
%% 对图像模糊加噪
Blurred = imfilter(I,PSF,'conv');
V = .03;
BlurredNoisy = imnoise(Blurred,'gaussian',0,V);
figure;imshow(BlurredNoisy);title('模糊加噪的图像');
%% 最小二乘方法恢复
NP = V*numel(I);
[reg1 LAGRA] = deconvreg(BlurredNoisy,PSF,NP);
figure,imshow(reg1),title('恢复图像');
imwrite(I,'1.bmp','bmp');
imwrite(Blurred,'2.bmp','bmp');
imwrite(reg1,'3.bmp','bmp');
%% 大仙QQ:1960009019
%% 在线教育微信公众号:大仙一品堂