14-2.m

来自「使用matlab实现图像处理的几个程序」· M 代码 · 共 19 行

I = checkerboard(8);
% 创建一个棋盘图像
PSF = fspecial('gaussian',7,10);
%创建点扩散函数PSF
V = .01;
BlurredNoisy = imnoise(imfilter(I,PSF),'gaussian',0,V);
%对模糊图像添加高斯噪声
NOISEPOWER = V*prod(size(I));
%噪声功率
[J LAGRA] = deconvreg(BlurredNoisy,PSF,NOISEPOWER);
subplot(221); imshow(BlurredNoisy);
title('A = Blurred and Noisy');
subplot(222); imshow(J);
title('[J LAGRA] = deconvreg(A,PSF,NP)');
subplot(223); imshow(deconvreg(BlurredNoisy,PSF,[],LAGRA/10));
title('deconvreg(A,PSF,[],0.1*LAGRA)');
subplot(224); imshow(deconvreg(BlurredNoisy,PSF,[],LAGRA*10));
title('deconvreg(A,PSF,[],10*LAGRA)');