码迷,mamicode.com
首页 > 其他好文 > 详细

[数字图像处理]图像去噪初步(1)--均值滤波器

时间:2014-07-12 17:37:54      阅读:259      评论:0      收藏:0      [点我收藏+]

标签:图像去噪   图像处理   图像复原   matlab   

1.图像去噪的前言


    上一篇博文中,我对噪声的类型进行了介绍,也使用的Matlab对各种噪声进行了实现。旧话重提,一幅图像,甚至是一个信号的老化,可以使用以下模型来表示表示。

bubuko.com,布布扣

可以使用以下算式来表示

bubuko.com,布布扣

这里,由于退化函数bubuko.com,布布扣的作用,使得原图像bubuko.com,布布扣产生退化(比如,运动模糊),然后在加上一个加性噪声项bubuko.com,布布扣

      本博文,主要对去除加性噪声的滤波器的性能进行了比较。对于退化函数的去除(称为去卷积或者逆滤波),将放在下一篇博文。


      1.1 实验用图像 

bubuko.com,布布扣


      1.2 实验结果的评价

            实验的步骤为,将实验用图像加上加性噪声,然后使用滤波器进行去噪,比较所得到的图像的画质。这里,就涉及到画质的评价方法。一般的,去噪图像的评价一般使用PSNR(峰值信噪比)。

bubuko.com,布布扣

       对于8-bit的图片而言,这里的MAX为255。PSNR越大,其画质就越好。但是,有些时候,使用PSNR来进行评价,也有不太合理的时候。

       请对比如下三张图片,a)是使用平均滤波器进行了处理,使其有些模糊;b)是使用高斯噪声污染原图;c)是使用椒盐噪声污染的图像。

bubuko.com,布布扣

    问题来了,这三张图像哪张画质最好,哪张最差。嘛,按照大众的答案而言,画质从好到差,应该是

a) > c) > b)

这样的排名吧(更多实际例子,请参考https://ece.uwaterloo.ca/~z70wang/research/ssim/)。那么,我们求其的PSNR是这样的。

bubuko.com,布布扣

这明显不科学,三幅图像的PSNR是一样的。反观PSNR的计算式,PSNR计算的时候,使用了MSE这个量。而MSE仅仅表现了两幅图像的灰度值的差,而对于图像的结构,却没有进行任何分析。

       这里使用一种比较好的图像画质评价的方法:SSIM(念做:艾斯-希姆)。这是一种由两张图像的灰度差异,构造差异和对比度去判断两张图的接近程度的方法。详情请参考[文献1],这里只做简单的介绍一下啦。

bubuko.com,布布扣

SSIM从图像亮度(Luminance),图像对比度(Contrast)和图像构造(Structure)去判断处理过的图像与原图的差异。这里,使用了某个区域的内的平均值作为亮度度量,使用方差作为对比度度量,使用协方差作为构造度量,来进行判断。这样,SSIM就比仅使用灰度去判断的PSNR更加准确。一样的,使用SSIM求取上面三幅图象的类似度。

bubuko.com,布布扣

从上表可以看出来,通过使用SSIM进行判断的结果,更加符合主观感受。本文余下的实验,全部使用SSIM去判断画质。


2.几个均值滤波器---线性处理


    2.1 算术均值滤波器

    算术均值滤波器很简单,就是求某个区域内的所有像素的均值,可用下式子表示。
    bubuko.com,布布扣
从式子上可以看出来,这就是一个低通滤波器,会使得画面模糊,有些许去噪能力。稍微做个实验看看。

    将实验用图像加噪,噪声均值为0、方差为0.0298的噪声。

bubuko.com,布布扣

    下面,使用算术均值滤波器,看看去噪效果。

bubuko.com,布布扣

   被去掉了些许,只是些许。再看频率域内的图像,果然是一个低通滤波器,我们都可以脑补出这个滤波器的振幅特性了,对吧?


    2.2 几何均值滤波器

    接下来是几何均值滤波器,很简单,求某个区域内的集合平均值。
    bubuko.com,布布扣
对于这个滤波器,书(《Digital Image Processing》 Rafael C. Gonzalez / Richard E. Woods)上说了,这个滤波器产生的去噪效果可以与算术平均值滤波器基本一样,但是可以更少的丢失细节。
    同样的,将实验用图像加噪,噪声均值为0、方差为0.0298的噪声。
bubuko.com,布布扣
       有结果可见,其去噪效果也不是太理想,但是原本芯片的pin脚什么的,还是比算术均值滤波器要一些的。嘛,当然,几何滤波器有一个致命的缺点,一旦有0值出现,那么这个像素的值立即被决定为0,这也就意味着,几何滤波器不可以去除胡椒噪声。

    2.3 算术均值滤波器与几何均值滤波器的比较

        为了对比几何均值滤波器与算术均值滤波器,我们进行了如下几组实验。由于篇幅问题,我就不贴出图来了。看数据就可以了。
        1.  噪声:高斯噪声,均值0,方差为0.17
bubuko.com,布布扣
         2.噪声:高斯噪声,均值0.2,方差为0.17
bubuko.com,布布扣
         3.噪声:椒盐噪声,胡椒密度0,盐粒密度0.1
bubuko.com,布布扣
这个还是想把结果贴上,几何平均滤波器的实验结果还是具有一定的观赏性的。
bubuko.com,布布扣
         4.噪声:椒盐噪声,胡椒密度0.1,盐粒密度0
bubuko.com,布布扣
       就如同实验结果一样,由于包含了大量了胡椒噪声,几何滤波器坏掉了。得到的结果很糟糕。

       实验结论:实验4的数据说明,由于包含了0值,几何滤波器去噪效果并不太好。但是实验3,几何滤波器的去噪效果真的是很不错。简单而言,算术均值滤波器泛用性比较好,而几何滤波器则擅长于去除盐粒噪声。
================
       拓展运用:既然几何滤波器对于去除盐粒噪声,那么,对于仅仅函数胡椒噪声的图像取反,将胡椒噪声转换为盐粒噪声去处理,所得结果再返回来,那么,几何均值滤波器还是可以去除胡椒噪声的。
================

    2.4 谐波均值滤波器

    其表达式如下所示,
bubuko.com,布布扣
注意式子的分母,这个滤波器不但不能去除椒盐噪声,对于灰度过黑的图像而言,也是要坏事的。书(《Digital Image Processing》 Rafael C. Gonzalez / Richard E. Woods)上又说了,这个滤波器,擅长于高斯噪声的去噪。实验一下看看。
    把实验用图加噪,类型为高斯噪声,均值为0,方差为0.15.
bubuko.com,布布扣
从实验结果看,其实这个滤波器的去噪效果还不如算术平均滤波器和几何平均滤波器。

    2.5 逆谐波均值滤波器

    逆谐波滤波器可以对应多种噪声,式子如下。
bubuko.com,布布扣
这个滤波器,可以通过Q值的变化,来获得一定的效果。当Q为正,这个滤波器可以去除胡椒噪声,当Q为负,这个滤波器可以去除盐粒噪声。贴个结果收活吧。
    bubuko.com,布布扣

3.总结

   本文介绍了几种图像去噪的线性滤波器,并对他们进行了比较。本次实验所用到的代码,我贴在下面。
   另外,评价函数ssim_index(),请自行网上查找。


参考文献

[1] Z. Wang, A. C. Bovik, H. R. Sheikh and E. P. Simoncelli, "Image quality assessment: From error visibility to structural similarity," IEEE Transactions on Image Processing, vol. 13, no. 4, pp. 600-612, Apr. 2004.


原文发于博客:http://blog.csdn.net/thnh169/ 

[数字图像处理]图像去噪初步(1)--均值滤波器,布布扣,bubuko.com

[数字图像处理]图像去噪初步(1)--均值滤波器

标签:图像去噪   图像处理   图像复原   matlab   

原文地址:http://blog.csdn.net/zhoufan900428/article/details/37698107

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!