去除噪声 -凯发k8网页登录

通过线性滤波去除噪声

您可以使用线性滤波来去除某些类型的噪声。某些滤波器，如平均值滤波器或高斯滤波器，适用于此目的。例如，平均值滤波器可用于从照片中去除颗粒噪声。由于每个像素设置为其邻域中像素的平均值，因此会减少由颗粒引起的局部变化。

有关使用 进行线性滤波的详细信息，请参阅 。

使用平均值滤波器和中位数滤波器去除噪声

此示例说明如何使用平均值滤波器和中位数滤波器从图像中去除椒盐噪声，以便比较结果。这两种类型的滤波都将输出像素的值设置为对应输入像素周围邻域中像素值的平均值。然而，使用中位数滤波时，输出像素的值由邻域像素的中位数（而不是均值）决定。中位数远不如均值对极值（称为离群值）敏感。因此，中位数滤波能够更好地去除这些离群值，而不会降低图像的锐度。

注意：中位数滤波是顺序统计量滤波的一种特殊情况，也称为秩滤波。有关顺序统计量滤波的信息，请参阅 ordfilt2 函数的参考页。

将图像读入工作区并显示它。

i = imread('eight.tif');
figure
imshow(i)

对于此示例，向图像中添加椒盐噪声。这种类型的噪声由设置为黑色或白色（数据范围的极值）的随机像素组成。

j = imnoise(i,'salt & pepper',0.02);
figure
imshow(j)

使用平均值滤波器对含噪图像 j 进行滤波并显示结果。该示例使用 3×3 邻域。

kaverage = filter2(fspecial('average',3),j)/255;
figure
imshow(kaverage)

现在使用中位数滤波器对含噪图像 j 进行滤波。该示例还使用 3×3 邻域。将两个滤波后的图像并排显示，以进行比较。请注意，medfilt2 在去除噪声方面做得更好，硬币边缘的模糊更少。

kmedian = medfilt2(j);
imshowpair(kaverage,kmedian,'montage')

通过自适应滤波去除噪声

此示例说明如何使用 wiener2 函数将 wiener 滤波器（一种线性滤波器）以自适应方式应用于图像。wiener 滤波器可自行适应图像局部方差。当方差较大时，wiener2 几乎不执行平滑处理。当方差较小时，wiener2 执行更多平滑处理。

这种方法通常比线性滤波产生更好的结果。自适应滤波器相比类似的线性滤波器更具选择性，它可保留图像的边缘和其他高频部分。此外，它没有设计任务；wiener2 函数处理所有初步计算，并对输入图像实现滤波器。然而，与线性滤波相比，wiener2 确实需要更多计算时间。

当噪声是恒定功率（“白色”）加性噪声（如高斯噪声）时，wiener2 效果最佳。以下示例将 wiener2 应用于添加了高斯噪声的土星图像。

将图像读入工作区中。

rgb = imread('saturn.png');

将图像从真彩色转换为灰度。

i = im2gray(rgb);

向图像中添加高斯噪声

j = imnoise(i,'gaussian',0,0.025);

显示含噪图像。由于图像相当大，因此只显示图像的一部分。

imshow(j(600:1000,1:600));
title('portion of the image with added gaussian noise');

使用 wiener2 函数去除噪声。

k = wiener2(j,[5 5]);

显示处理后的图像。由于图像相当大，因此只显示图像的一部分。

figure
imshow(k(600:1000,1:600));
title('portion of the image with noise removed by wiener filter');