AverageFiltering.zip_图像去噪_均值滤波

close all; clear all; clc; %1.lena图像 %高斯噪声均值去噪 I=imread('lena.png'); imwrite(I,'lena原图.bmp'); figure,imagesc(I),title('原图');colormap('gray'); [m n]=size(I); v=60; I1=imnoise(I,'gaussian',0,v*v/(m*n));%添加高斯白噪声 figure,imagesc(I1),title('高斯加噪图像');colormap('gray'); d1=avg_filter(I1,3); figure,imagesc(d1),title('高斯噪声均值去噪图像');colormap('gray'); imwrite(d1,'lena+gaussian均值去噪图像.bmp'); % %计算psnr的值 x1=imread('lena原图.bmp'); x2=imread('lena+gaussian均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2,lena); %输入数据的mse为:2963.6596 % %输出的psnrvalue的值为:13.4125 psnrvalue=psnr(x1,x2,'lena+gaussian denoising'); % psnrvalue=psnr(x1,I1,'lena+gaussian denoising'); %计算ssim的值 [mssim, ssim_map] = ssim(x1, x2); % [mssim, ssim_map] = ssim(x1, I1); disp(['lena+gaussian SSIM:',num2str(mssim)]); figure,imshow(max(0, ssim_map).^4); %椒盐噪声均值去噪 I2=imnoise(I,'salt & pepper',0.08);%添加椒盐噪声 figure,imagesc(I2),title('椒盐噪声图像');colormap('gray'); d2=avg_filter(I2,3); figure,imagesc(d2),title('椒盐噪声均值去噪图像');colormap('gray'); imwrite(d2,'lena椒盐噪声均值去噪图像.bmp'); % %计算psnr的值 x1=imread('lena原图.bmp'); x2=imread('lena椒盐噪声均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:1906.2472 % %输出的psnrvalue的值为:15.329 psnrvalue=psnr(x1,x2,'lena+salt&pepper denoising'); % psnrvalue=psnr(x1,I2,'lena+salt&pepper denoising'); %计算ssim的值 [mssim, ssim_map] = ssim(x1, x2); % [mssim, ssim_map] = ssim(x1, I2); disp(['lena+salt&pepper SSIM:',num2str(mssim)]); figure,imshow(max(0, ssim_map).^4); % %2.barbara图像 % %高斯噪声均值去噪 % I=imread('barbara.png'); % figure,imagesc(I),title('原图');colormap('gray'); % I1=imnoise(I,'gaussian',0,0.02);%添加高斯白噪声 % figure,imagesc(I1),title('加噪图像');colormap('gray'); % d1=avg_filter(I1,3); % figure,imagesc(d1),title('均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('barbara原图.bmp'); % x2=imread('barbara均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:3913.8613 % %输出的psnrvalue的值为:12.2047 % % %椒盐噪声均值去噪 % I2=imnoise(I,'salt & pepper',0.02);%添加椒盐噪声 % figure,imagesc(I2),title('椒盐噪声图像');colormap('gray'); % d2=avg_filter(I2,3); % figure,imagesc(d2),title('椒盐噪声均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('barbara原图.bmp'); % x2=imread('barbara椒盐噪声均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:3145.4048 % %输出的psnrvalue的值为:13.154 % % % %3.boat图像 % %高斯噪声均值去噪 % I=imread('boat.png'); % figure,imagesc(I),title('原图');colormap('gray'); % I1=imnoise(I,'gaussian',0,0.02);%添加高斯白噪声 % figure,imagesc(I1),title('加噪图像');colormap('gray'); % d1=avg_filter(I1,3); % figure,imagesc(d1),title('均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('boat原图.bmp'); % x2=imread('boat均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:3308.7119 % %输出的psnrvalue的值为:12.9342 % % %椒盐噪声均值去噪 % I2=imnoise(I,'salt & pepper',0.02);%添加椒盐噪声 % figure,imagesc(I2),title('椒盐噪声图像');colormap('gray'); % d2=avg_filter(I2,3); % figure,imagesc(d2),title('椒盐噪声均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('boat原图.bmp'); % x2=imread('boat椒盐噪声均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:2029.0444 % %输出的psnrvalue的值为:15.0579 % % % %4.cameraman图像 % %高斯噪声均值去噪 % I=imread('cameraman.png'); % figure,imagesc(I),title('原图');colormap('gray'); % I1=imnoise(I,'gaussian',0,0.02);%添加高斯白噪声 % figure,imagesc(I1),title('加噪图像');colormap('gray'); % d1=avg_filter(I1,3); % figure,imagesc(d1),title('均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('cameraman原图.bmp'); % x2=imread('cameraman均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:4935.7704 % %输出的psnrvalue的值为:11.1973 % % % %椒盐噪声均值去噪 % I2=imnoise(I,'salt & pepper',0.02);%添加椒盐噪声 % figure,imagesc(I2),title('椒盐噪声图像');colormap('gray'); % d2=avg_filter(I2,3); % figure,imagesc(d2),title('椒盐噪声均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('cameraman原图.bmp'); % x2=imread('cameraman椒盐噪声均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:2390.0792 % %输出的psnrvalue的值为:14.3467 % % % %5.house图像 % %高斯噪声均值去噪 % I=imread('house.png'); % figure,imagesc(I),title('原图');colormap('gray'); % I1=imnoise(I,'gaussian',0,0.02);%添加高斯白噪声 % figure,imagesc(I1),title('加噪图像');colormap('gray'); % d1=avg_filter(I1,3); % figure,imagesc(d1),title('均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('house原图.bmp'); % x2=imread('house均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:3567.2125 % %输出的psnrvalue的值为:12.6075 % % % %椒盐噪声均值去噪 % I2=imnoise(I,'salt & pepper',0.02);%添加椒盐噪声 % figure,imagesc(I2),title('椒盐噪声图像');colormap('gray'); % d2=avg_filter(I2,3); % figure,imagesc(d2),title('椒盐噪声均值去噪图像');colormap('gray'); % % %计算psnr的值 % x1=imread('house原图.bmp'); % x2=imread('house椒盐噪声均值去噪图像.bmp'); % psnrvalue=psnr(x1,x2); %输入数据的mse为:2394.4944 % %输出的psnrvalue的值为:14.3387