信噪比

video=mmreader('che2.AVI');
 I1=read(video,590);
I1=rgb2gray(I1);
[m,n]=size(I1);
J=imnoise(I1,'gaussian',0,0.002); %加入高斯噪声
%进行均值滤波
h=fspecial('average',3);    %fspecial函数用于产生预定义滤波器
I2=uint8(round(filter2(h,I1)));   %filter2函数用于图像滤波,此处h是滤波参数（均值），I1是要处理的图像
%进行中值滤波
I3=medfilt2(J,[3,3]);   %medfilt2函数用于图像的中值滤波
%进行维纳滤波
I4=wiener2(J,[3,3]);%进行一次维纳滤波
I5=wiener2(I4,[3,3]);%进行二次维纳滤波
subplot(2,3,1),imshow(I1),title('原图象')
subplot(2,3,2),imshow(J),title('加噪声图象')
subplot(2,3,3),imshow(I2),title('均值滤波后图象')
subplot(2,3,4),imshow(I3),title('中值滤波后图象')
subplot(2,3,5),imshow(I4),title('维纳滤波后图象')
I1=double(I1);
I2=double(I2);
I3=double(I3);
I4=double(I4);
I5=double(I5);
s2=0;
s2=double(s2);
s3=0;
s3=double(s3);
s4=0;
s4=double(s4);
s5=0;
s5=double(s5);
for i=1:m
    for j=1:n

      s2=s2+(I1(i,j)-I2(i,j))^2;
   s3=s3+(I1(i,j)-I3(i,j))^2;
   s4=s4+(I1(i,j)-I4(i,j))^2;
   s5=s5+(I1(i,j)-I5(i,j))^2;
    end
end
MSE2=s2/(m*n)
MSE3=s3/(m*n)
MSE4=s4/(m*n)
MSE5=s5/(m*n)
 PSNR2=10*log10(double(255^2/MSE2))
  PSNR3=10*log10(double(255^2/MSE3))
   PSNR4=10*log10(double(255^2/MSE4))
    PSNR5=10*log10(double(255^2/MSE5))