达卡人技术博客™
        dakaren技术博客,博广天下–源码算法|图像增强|经验总结|设计方案|学习教程！ >>> 本站作为一个技术站点，并不是简单复制转载一些文章而已，所以不可能每天都更新，也没有那么多精力，希望体谅，但我们会努力和大家交流最有用的东西。也欢迎您来稿来信，共同进步。谢谢dajia's支持! 给大家推荐一个不错的CV网站>>计算机视觉论坛：http://cvchina.net

中值滤波（median filter）是一种有效消除椒盐噪声的算法。
基本原理是，对图像中所有点进行遍历，对于每个点，取以其为中心的 n * n 的矩形区域，对矩形区域中点的灰度值进行排序，取其中间值替换当前点。
这样，如果这个点灰度值与周围点相差较大的话，就可以将其平滑化。
当然，这里要求 n 是个大于 1 的奇数。

在这里，我对算法进行了优化，使只有当前点是矩形区域灰度值的最大或者最小值的时候才进行替换，这样效果好了很多。
原始图像：

椒盐噪声：

中值滤波后：

可以看到，图像损失很小，只有边缘处有小的细节损失。同时，由于所加的噪点比较密集，少数地方仍然有噪点，这是因为矩形区域噪点太多，以至于中值本身也是个噪点了。

源代码：

    BmpPixmap &
BmpPixmap::median_filter (int n)
{
    assert (n >= 3 && n % 2);

    int ii, jj, nn;
    Byte model [3][n * n];
    BmpPixmap *temp = new BmpPixmap (*this);
    for (i = n / 2; i < height - n / 2; i++) {
        for (j = n / 2; j < width - n / 2; j++) {
            /*-----------------------------------------------------------------------------
             *  Put n*n pixels around current pixels to the model.
             *-----------------------------------------------------------------------------*/
            for (ii = i - n / 2, nn = 0; ii <= i + n / 2; ii++) {
                for (jj = j - n / 2; jj <= j + n / 2; jj++, nn++) {
                    model [0][nn] = pdata [ii][jj]->get_blue ();
                    model [1][nn] = pdata [ii][jj]->get_green ();
                    model [2][nn] = pdata [ii][jj]->get_red ();
                }
            }
            /*-----------------------------------------------------------------------------
             *  Sort the model.
             *-----------------------------------------------------------------------------*/
            qsort (model [0], nn, sizeof (Byte), cmp_Byte);
            qsort (model [1], nn, sizeof (Byte), cmp_Byte);
            qsort (model [2], nn, sizeof (Byte), cmp_Byte);

            if (pdata [i][j]->get_blue () == model [0][nn - 1] ||
                    pdata [i][j]->get_blue () == model [0][0] ||
                    pdata [i][j]->get_green () == model [1][nn - 1] ||
                    pdata [i][j]->get_green () == model [1][0] ||
                    pdata [i][j]->get_red () == model [2][nn - 1] ||
                    pdata [i][j]->get_red () == model [2][0])
            {
                temp->pdata [i][j]->set (model [0][nn / 2], model [1][nn / 2], model [2][nn / 2]);
            }
        }
    }
    return *temp;
}  /* -----  end of method BmpPixmap::median_filter  ----- */

转载请注明：转自K.I.S.S. – 简单哲学 – 数字图像处理（二）——中值滤波

Related posts:

1. 数字图像处理（四）——均值滤波
2. Wiener2滤波C++代码
3. 常用的数字图像去噪典型算法及matlab实现