首页 > 代码库 > 每日一练之自适应中值滤波器(基于OpenCV实现)

每日一练之自适应中值滤波器(基于OpenCV实现)

本文主要介绍了自适应的中值滤波器,并基于OpenCV实现了该滤波器,并且将自适应的中值滤波器和常规的中值滤波器对不同概率的椒盐噪声的过滤效果进行了对比。最后,对中值滤波器的优缺点了进行了总结。

空间滤波器

一个空间滤波器包括两个部分:

  • 一个邻域,滤波器进行操作的像素集合,通常是一个矩形区域
  • 对邻域中像素进行的操作

一个滤波器就是在选定的邻域像素上执行预先定义好的操作产生新的像素,并用新的像素替换掉原来像素形成新的图像。
通常,也可以将滤波器称之为核(kernel),模板(template)或者窗口(window)

根据预定义的操作,可以将滤波器分为:

  • 线性滤波器
  • 非线性滤波器

而根据滤波器最终对图像造成的影响,可以将滤波器分为:

  • 平滑滤波器 ,通常用于模糊图像或者去除图像中的噪声
  • 锐化滤波器,突出图像中的边缘细节部分

中值滤波器 Median Filter

中值滤波器是一种常用的非线性滤波器,其基本原理是选择待处理像素的一个邻域中各像素值的中值来代替待处理的像素,其主要功能是像素的灰度值与周围像素比较接近,从而消除孤立的噪声点,所以中值滤波器能够很好的消除椒盐噪声。不仅如此,中值滤波器在消除噪声的同时,还能有效的保护图像的边界信息,不会对图像造成很大的模糊(相比于均值滤波器)。

中值滤波器的效果受滤波窗口尺寸的影响较大,在消除噪声和保护图像的细节存在着矛盾:滤波窗口较小,则能很好的保护图像中的某些细节,但对噪声的过滤效果就不是很好;反之,窗口尺寸较大有较好的噪声过滤效果,但是会对图像造成一定的模糊。另外,根据中值滤波器原理,如果在滤波窗口内的噪声点的个数大于整个窗口内像素的个数,则中值滤波就不能很好的过滤掉噪声。

自适应中值滤波器 Adaptive Median Filter

上面提到常规的中值滤波器,在噪声的密度不是很大的情况下(根据经验,噪声的出现的概率小于0.2),效果不错。但是当概率出现的概率较高时,常规的中值滤波的效果就不是很好了。有一个选择就是增大滤波器的窗口大小,这虽然在一定程度上能解决上述的问题,但是会给图像造成较大的模糊。

常规的中值滤波器的窗口尺寸是固定大小不变的,就不能同时兼顾去噪和保护图像的细节。这时就要寻求一种改变,根据预先设定好的条件,在滤波的过程中,动态的改变滤波器的窗口尺寸大小,这就是自适应中值滤波器 Adaptive Median Filter。在滤波的过程中,自适应中值滤波器会根据预先设定好的条件,改变滤波窗口的尺寸大小,同时还会根据一定的条件判断当前像素是不是噪声,如果是则用邻域中值替换掉当前像素;不是,则不作改变。

自适应中值滤波器有三个目的:

  • 滤除椒盐噪声
  • 平滑其他非脉冲噪声
  • 尽可能的保护图像中细节信息,避免图像边缘的细化或者粗化。

自使用中值滤波算法描述

自适应滤波器不但能够滤除概率较大的椒盐噪声,而且能够更好的保护图像的细节,这是常规的中值滤波器做不到的。自适应的中值滤波器也需要一个矩形的窗口Sxy<script type="math/tex" id="MathJax-Element-1">S_{xy}</script>,和常规中值滤波器不同的是这个窗口的大小会在滤波处理的过程中进行改变(增大)。需要注意的是,滤波器的输出是一个像素值,该值用来替换点(x,y)<script type="math/tex" id="MathJax-Element-2">(x,y)</script>处的像素值,点(x,y)<script type="math/tex" id="MathJax-Element-3">(x,y)</script>是滤波窗口的中心位置。

在描述自适应中值滤波器时需要用到如下的符号:

  • Zmin=Sxy<script type="math/tex" id="MathJax-Element-4">Z_{min}=S_{xy}</script>中的最小灰度值
  • Zmax=Sxy<script type="math/tex" id="MathJax-Element-5">Z_{max}=S_{xy}</script>中的最大灰度值
  • Zmed=Sxy<script type="math/tex" id="MathJax-Element-6">Z_{med}=S_{xy}</script>中的灰度值的中值
  • Zxy<script type="math/tex" id="MathJax-Element-7">Z_{xy}</script>表示坐标(x,y)<script type="math/tex" id="MathJax-Element-8">(x,y)</script>处的灰度值
  • Smax=Sxy<script type="math/tex" id="MathJax-Element-9">S_{max}=S_{xy}</script>允许的最大窗口尺寸

自适应中值滤波器有两个处理过程,分别记为:A和B。
A :
A1 = Zmed?Zmin<script type="math/tex" id="MathJax-Element-10">Z_{med}-Z_{min}</script>
A2 = Zmed?Zmax<script type="math/tex" id="MathJax-Element-11">Z_{med}-Z_{max}</script>
如果A1 > 0 且 A2 < 0,跳转到 B;
否则,增大窗口的尺寸
如果增大后窗口的尺寸 Smax<script type="math/tex" id="MathJax-Element-12">\leq S_{max}</script>,则重复A过程。
否则,输出Zmed<script type="math/tex" id="MathJax-Element-13">Z_{med}</script>

B:
B1 = Zxy?Zmin<script type="math/tex" id="MathJax-Element-14">Z_{xy}-Z_{min}</script>
B2 = Zxy?Zmax<script type="math/tex" id="MathJax-Element-15">Z_{xy}-Z_{max}</script>
如果B1 > 0 且 B2 < 0,则输出Zxy<script type="math/tex" id="MathJax-Element-16">Z_{xy}</script>
否则输出Zmed<script type="math/tex" id="MathJax-Element-17">Z_{med}</script>

自适应中值滤波原理说明

过程A的目的是确定当前窗口内得到中值Zmed<script type="math/tex" id="MathJax-Element-18">Z_{med}</script>是否是噪声。如果Zmin<Zmed<Zmax<script type="math/tex" id="MathJax-Element-19">Z_{min} < Z_{med} < Z_{max}</script>,则中值Zmed<script type="math/tex" id="MathJax-Element-20">Z_{med}</script>不是噪声,这时转到过程B测试,当前窗口的中心位置的像素Zxy<script type="math/tex" id="MathJax-Element-21">Z_{xy}</script>是否是一个噪声点。如果Zmin<Zxy<Zmax<script type="math/tex" id="MathJax-Element-22">Z_{min} < Z_{xy} < Z_{max}</script>,则Zxy<script type="math/tex" id="MathJax-Element-23">Z_{xy}</script>不是一个噪声,此时滤波器输出Zxy<script type="math/tex" id="MathJax-Element-24">Z_{xy}</script>;如果不满足上述条件,则可判定Zxy<script type="math/tex" id="MathJax-Element-25">Z_{xy}</script>是噪声,这是输出中值Zmed<script type="math/tex" id="MathJax-Element-26">Z_{med}</script>(在A中已经判断出Zmed<script type="math/tex" id="MathJax-Element-27">Z_{med}</script>不是噪声)。

如果在过程A中,得到则Zmed<script type="math/tex" id="MathJax-Element-28">Z_{med}</script>不符合条件Zmin<Zmed<Zmax<script type="math/tex" id="MathJax-Element-29">Z_{min} < Z_{med} < Z_{max}</script>,则可判断得到的中值Zmed<script type="math/tex" id="MathJax-Element-30">Z_{med}</script>是一个噪声。在这种情况下,需要增大滤波器的窗口尺寸,在一个更大的范围内寻找一个非噪声点的中值,直到找到一个非噪声的中值,跳转到B;或者,窗口的尺寸达到了最大值,这时返回找到的中值,退出。

从上面分析可知,噪声出现的概率较低,自适应中值滤波器可以较快的得出结果,不需要去增加窗口的尺寸;反之,噪声的出现的概率较高,则需要增大滤波器的窗口尺寸,这也符合种中值滤波器的特点:噪声点比较多时,需要更大的滤波器窗口尺寸。

实现

有了算法的详细描述,借助于OpenCV对图像的读写,自适应中值滤波器实现起来也不是很困难。

    int minSize = 3; // 滤波器窗口的起始尺寸
    int maxSize = 7; // 滤波器窗口的最大尺寸
    Mat im1;
    // 扩展图像的边界
    copyMakeBorder(im, im1, maxSize / 2, maxSize / 2, maxSize / 2, maxSize / 2, BorderTypes::BORDER_REFLECT);
    // 图像循环
    for (int j = maxSize / 2; j < im1.rows - maxSize / 2; j++)
    {
        for (int i = maxSize / 2; i < im1.cols * im1.channels() - maxSize / 2; i++)
        {
            im1.at<uchar>(j, i) = adaptiveProcess(im1, j, i, minSize, maxSize);
        }
    }

首先定义滤波器最小的窗口尺寸以及最大的窗口尺寸。
要进行滤波处理,首先要扩展图像的边界,以便对图像的边界像素进行处理。copyMakeBorder根据选择的BorderTypes使用不同的值扩充图像的边界像素,具体可参考OpenCV的文档信息。
下面就是遍历图像的像素,对每个像素进行滤波处理。需要注意一点,不论滤波器多么的复杂,其每次的滤波过程,都是值返回一个值,来替换掉当前窗口的中心的像素值。函数adpativeProcess就是对当前像素的滤波过程,其代码如下:

uchar adaptiveProcess(const Mat &im, int row,int col,int kernelSize,int maxSize)
{
    vector<uchar> pixels;
    for (int a = -kernelSize / 2; a <= kernelSize / 2; a++)
        for (int b = -kernelSize / 2; b <= kernelSize / 2; b++)
        {
            pixels.push_back(im.at<uchar>(row + a, col + b));
        }
    sort(pixels.begin(), pixels.end());
    auto min = pixels[0];
    auto max = pixels[kernelSize * kernelSize - 1];
    auto med = pixels[kernelSize * kernelSize / 2];
    auto zxy = im.at<uchar>(row, col);
    if (med > min && med < max)
    {
        // to B
        if (zxy > min && zxy < max)
            return zxy;
        else
            return med;
    }
    else
    {
        kernelSize += 2;
        if (kernelSize <= maxSize)
            return adpativeProcess(im, row, col, kernelSize, maxSize); // 增大窗口尺寸,继续A过程。
        else
            return med;
    }
}

首先,根据当前窗口的大小,取得所有像素值存放到vector中,然后对vector进行排序,取得像素的最小值、最大值和中值。然后测试当前取得的中值是否在(min,max)之间,如果是,则中值不是噪声点,则开始对当前像素值进行处理,判断其是否是噪声点。如果,测试当前已取得的中值是噪声点,则扩大窗口的尺寸,在更大的空间中重新寻找中值。

上面自适应中值滤波器实现起来比较简单,所以问题就来了:效率及其的低下。这里,这是对自适应中值滤波器的原理的学习,可以忽略这个不必要的细节。

结果对比

技术分享
左边是添加概率为0.2的椒盐噪声,右边是原图。下面是使用常规的中值滤波和本文实现的自适应中值滤波器后的处理结果

技术分享

左边是自适应中值滤波器(最小窗口为3,最大窗口为7)的结果,右图是常规中值滤波器(窗口大小为5)的结果。可以看出,无论是中值滤波还是自适应的中值滤波,都能过滤掉图像中的噪声,自适应中值滤波器的效果要好些,常规的还有一些噪声没有过滤掉。而且,常规的中值滤波器对图像造成的模糊较明显,而自适应中值滤波器很好的的保存了图像中的细节。

下面测试更大概率噪声下,两种滤波器的工作情况。噪声概率为0.4时,
技术分享
可以看出,常规的中值滤波器已经不能很好的过滤掉噪声,而自适应的中值滤波还可以胜任。

中值滤波器总结

  • 中值滤波器能够很好的滤除“椒盐”噪声。椒盐噪声是在图像上随机出现的孤立点,根据中值滤波器的原理,使用邻域像素的中值代替原像素,能够有效的消除这些孤立的噪声点。
  • 和均值滤波器相比,中值滤波在消除噪声的同时,还能在很大程度保护图像的细节,不会造成很大的模糊。
  • 和常规的中值滤波器相比,自适应中值滤波器能够更好的保护图像中的边缘细节部分。
<script type="text/javascript"> $(function () { $(‘pre.prettyprint code‘).each(function () { var lines = $(this).text().split(‘\n‘).length; var $numbering = $(‘
    ‘).addClass(‘pre-numbering‘).hide(); $(this).addClass(‘has-numbering‘).parent().append($numbering); for (i = 1; i <= lines; i++) { $numbering.append($(‘
  • ‘).text(i)); }; $numbering.fadeIn(1700); }); }); </script>

    每日一练之自适应中值滤波器(基于OpenCV实现)