图像复原的Matlab实现(源代码)

附录：源代码：
1、高斯噪声的添加以及滤波处理
I= imread('E: \lena.bmp','bmp');
J=imnoise(I,'gaussian',0,0.01);
figure;
subplot(1,2,1);
imshow(I);
title('源图像lena.bmp');
subplot(1,2,2);
imshow(J);
title('加入gaussian噪声后的lena.bmp');
n1=7;sigma1=1.5;n2=3;sigma2=1.5;theta=0;
r=[cos(theta) -sin(theta); sin(theta) cos(theta)];
for i = 1 : n2
for j = 1 : n1
u = r*[j-(n1+1)/2 i-(n2+1)/2]';
h(i,j)=exp(-u(1)^2/(2*sigma1^2))/(sigma1*sqrt(2*pi))*exp(-u(2)^2/
(2*sigma2^2))/(sigma2*sqrt(2*pi)); end
end
h = h / sqrt(sum(sum(h.*h)));
f1=conv2(J,h,'same');
subplot(1,2,2);
figure;
imagesc(f1);
title('高斯平滑后的lena.bmp(7x7)');
colormap(gray);
2、图像模糊及添加噪声
①图像的运动模糊
I= imread('E: \lena.bmp','bmp');
figure;
subplot(1,2,1);
imshow(I);
title('源图像lena.bmp');
f=double(I); % 数据类型转换，MATLAB不支持图像的无符号整型的计算
g=fft2(f); % 傅立叶变换
g=fftshift(g); % 转换数据矩阵
[M,N]=size(g);
a=0.5;b=0.5;T=0.1;
m=fix(M/2);
n=fix(N/2);
j=sqrt(-1);
for i=1:M
for k=1:N
h=(T/(pi*(i*a+k*b)))*sin(pi*(i*a+k*b))*exp(-j*pi*(i*a+k*b)); end
result(i,k)=h*g(i,k);
end
result=ifftshift(result);
J1=ifft(result);
J2=uint8(real(J1));
subplot(1,2,2);
imshow(J2);
title('模糊化lena.bmp');
②图像加噪
figure;subplot(1,2,1);
imshow(J2);
title('运动模糊后的lena.bmp(角度为45)');
J3=imnoise(J2,'gaussian',0,0.01);
subplot(1,2,2);
imshow(J3);
title('加噪并模糊的lena.bmp');
3、图像恢复
①维纳滤波恢复图像
I= imread('E: \lena.bmp','bmp');
H=fspecial('motion',50,45);
J=imfilter(I,H,'circular','conv');
figure;
subplot(1,2,1);
imshow(J);
title('运动模糊后的lena.bmp(角度为45)');
J1=imnoise(J,'gaussian',0,0.01);
subplot(1,2,2);
imshow(J1);
title('加噪并模糊的lena.bmp');
%figure;
J2=deconvwnr(J1,H),[]);
imshow(J2);
title('模糊噪声图像的维纳滤波复原');
figure;
noise=imnoise(zeros(size(I)),'gaussian',0,0.01);
NSR=sum(noise(:).^2)/sum(im2double(I(:)).^2);
J3=deconvwnr(J1,H,NSR);
imshow(J3);
title('引入SNR的维纳滤波复原');
②约束最小二乘方法恢复图像
I= imread('E:\lena.bmp','bmp');
I1=checkerboard(8);
PSF=fspecial('motion',50,45);
V=0.0001;
J=imfilter(I,PSF,'circular','conv'); J1=imnoise(J,'gaussian',0,0.01); figure;subplot(1,2,1);
imshow(J1);
title('模糊加噪图像');
NoisePower=V*prod(size(I));
[G,LAGRA]=deconvreg(J,PSF,NoisePower); subplot(1,2,2);
imshow(G);
title('约束最小二乘滤波复原');。