均值中值.zip

clear all I=imread('11187.bmp'); %读入JPG彩色图像文件 imshow(I) %显示出来 title('输入的彩色JPG图像') imwrite(rgb2gray(I),'原始图像灰度图1.bmp'); %将彩色图片灰度化并保存 I=rgb2gray(I); %灰度化后的数据存入数组 X = im2double(I); % 转为双精度类型 % x = noise(X, 'gaussian', 0.005); % 添加噪声 高斯噪声 方差 0.005 x = noise(X, 'gaussian', 0.01); % 添加噪声 高斯噪声 方差 0.01 % x = noise(X, 'guassian', 0.02); % 高斯噪声 方差 0.02 x = noise(x, 'salt & pepper', 0.02); % 椒盐噪声 密度 0.02 % x = noise(X, 'salt & pepper', 0.05); % 椒盐噪声 密度 0.05 f = ones(3); x_mean1 = meanfilter(x,f); % 均值滤波 x_med1 = medfilter(x_mean1, f); % 中值滤波 % 计算峰值信噪比 %psnr_mean1 = PSNR(x_mean1, X) %psnr_med1 = PSNR(x_med1, X) figure; imshow(x_mean1);title('均值滤波1'); figure; imshow(x_med1); title('中值滤波1');