织梦CMS - 轻松建站从此开始!

罗索

高斯平滑 高斯模糊 高斯滤波器 ( Gaussian Smoothing, Gaussian

落鹤生 发布于 2012-11-01 09:46 点击:次 
高斯平滑 高斯模糊 高斯滤波器 ( Gaussian Smoothing, Gaussian Blur, Gaussian Filter ) C++ 实现 在之前我们提到过了均值滤波器, 就是说某像素的颜色, 由以其为中心的九宫格的像素平均值来决定. 在这个基础上又发展成了带权的平均滤波器, 这里的高斯平滑或者说滤波器
TAG:

高斯平滑 高斯模糊 高斯滤波器 ( Gaussian Smoothing, Gaussian Blur, Gaussian Filter ) C++ 实现

发展到现在这个平滑算法的时候, 我已经完全不知道如何去命名这篇文章了, 只好罗列出一些关键字来方便搜索了.

在之前我们提到过了均值滤波器, 就是说某像素的颜色, 由以其为中心的九宫格的像素平均值来决定. 在这个基础上又发展成了带权的平均滤波器, 这里的高斯平滑或者说滤波器就是这样一种带权的平均滤波器. 那么这些权重如何分布呢? 我们先来看几个经典的模板例子:

尝试了使用这些滤波器对我们原来的图进行操作, 得到了这样的一组结果:

原图:

raw

3x3 高斯:

3x3

5x5 高斯:

5x5

 

单纯从效果来看, 两个模板都起到了平滑的作用, 只是程度有深浅的区分. 那么从理论上来说为什么能起到平滑的作用呢? 很显然, 像素的颜色不仅由自身决定了, 同时有其周围的像素加权决定, 客观上减小了和周围像素的差异. 同时这些权重的设定满足了越近权重越大的规律. 从理论来讲, 这些权重的分布满足了著名的所谓高斯分布:

  这就是1维的计算公式

这就是2维的计算公式

x, y表示的就是当前点到对应点的距离, 而那些具体的模板就是由这里公式中的一些特例计算而来. 需要说明的是不只有这么一些特例, 从wikipedia可以方便地找到那些复杂的模板比如像:

Sample Gaussian matrix

This is a sample matrix, produced by sampling the Gaussian filter kernel (with σ = 0.84089642) at the midpoints of each pixel and then normalising. Note that the center element (at [4, 4]) has the largest value, decreasing symmetrically as distance from the center increases.

0.00000067 0.00002292 0.00019117 0.00038771 0.00019117 0.00002292 0.00000067
0.00002292 0.00078633 0.00655965 0.01330373 0.00655965 0.00078633 0.00002292
0.00019117 0.00655965 0.05472157 0.11098164 0.05472157 0.00655965 0.00019117
0.00038771 0.01330373 0.11098164 0.22508352 0.11098164 0.01330373 0.00038771
0.00019117 0.00655965 0.05472157 0.11098164 0.05472157 0.00655965 0.00019117
0.00002292 0.00078633 0.00655965 0.01330373 0.00655965 0.00078633 0.00002292
0.00000067 0.00002292 0.00019117 0.00038771 0.00019117 0.00002292 0.00000067

是不是看到就头大了:) 不过没关系, 对于一般的应用来说, 前面的例子已经可以完成任务了.  代码的话我们还是给一份5x5的example:

  1. /** 
  2. ** method to remove noise from the corrupted image by gaussian filter value 
  3. * @param corrupted input grayscale binary array with corrupted info 
  4. * @param smooth output data for smooth result, 
  5. the memory need to be allocated outside of the function 
  6. * @param width width of the input grayscale image 
  7. * @param height height of the input grayscale image 
  8. */ 
  9. void gaussianFilter2 (unsigned char* corrupted, unsigned char* smooth
  10. int width, int height) 
  11.     int templates[25] = { 1, 4, 7, 4, 1,  
  12.                           4, 16, 26, 16, 4,  
  13.                           7, 26, 41, 26, 7, 
  14.                           4, 16, 26, 16, 4,  
  15.                           1, 4, 7, 4, 1 };       
  16.      
  17.     memcpy ( smooth, corrupted, width*height*sizeof(unsigned char) ); 
  18.     for (int j=2;j<height-2;j++) 
  19.     { 
  20.         for (int i=2;i<width-2;i++) 
  21.         { 
  22.             int sum = 0; 
  23.             int index = 0; 
  24.             for ( int m=j-2; m<j+3; m++) 
  25.             { 
  26.                 for (int n=i-2; n<i+3; n++) 
  27.                 { 
  28.                     sum += corrupted [ m*width + n] * templates[index++] ; 
  29.                 } 
  30.             } 
  31.             sum /= 273; 
  32.             if (sum > 255) 
  33.                 sum = 255; 
  34.             smooth [ j*width+i ] = sum; 
  35.         } 
  36.     } 

附带说一些,很明显,和均值滤波器类似, 这个滤波器没有消除校验噪声的作用.

(hhygcy)
本站文章除注明转载外,均为本站原创或编译欢迎任何形式的转载,但请务必注明出处,尊重他人劳动,同学习共成长。转载请注明:文章转载自:罗索实验室 [http://www.rosoo.net/a/201206/16120.html]
本文出处:CSDN博客 作者:hhygcy 原文
顶一下
(4)
100%
踩一下
(0)
0%
------分隔线----------------------------
相关文章
发表评论
请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
评价:
表情:
用户名: 验证码:点击我更换图片
栏目列表
将本文分享到微信
织梦二维码生成器
推荐内容