高斯模糊公式原理分析

《高斯模糊公式原理分析》由会员分享，可在线阅读，更多相关《高斯模糊公式原理分析（5页珍藏版）》请在装配图网上搜索。

1、最近有感于部分网友对高斯模糊滤镜的研究，现总结如下。高斯模糊是数字图像模板处理法 的一种。其模板是根据二维正态分布（高斯分布）函数计算出来的。正态分布最早由A棣 莫弗在求二项分布的渐近公式中得到，C.F.高斯在研究测量误差时从另一个角度导出了它, P.S.拉普拉斯和高斯研究了它的性质。一维正态分布的函数定义：在这个函数中，第一个参数M是遵从正态分布的随机变量的均值，第二个参数O2是此随机 变量的方差，所以正态分布记作N（p,o2）。遵从正态分布的随机变量其概率规律为：取|j 邻近的值的概率大，而取离j越远的值的概率越小；o越小，分布越集中在j附近，o越大， 分布越分散。正态分布的密度函数的特点

2、是：关于j对称，在j处达到最大值，在正（负） 无穷远处取值为0在jo处有拐点。它的形状是中间高两边低，图像是一条位于x轴上 方的钟形曲线。当j = 0，o2=1时，称为标准正态分布，记为N（0,1）。两个参数的意义：j 期望，o2方差。下面我们解决第一个疑问：高斯模糊滤镜中的半径是什么？答案是高斯半径就是公式中的 高斯曲线的图形和半径的含义如下图（来自Adobe SDK中技术支持专家的文档）所示:由此可见高斯半径（o）对曲线形状的影响：o越小，曲线越高越尖，o越大，曲线越低越 平缓。对二维图像来说，是一个钟形曲面，高斯半径越小，曲面越高越尖越陡峭；高斯半径 越大，曲面越低越平缓。因此高斯半径越

3、小，则模糊越小，高斯半径越大，则模糊程度越大。 我们将看到PS对高斯半径的范围定义是0.1250。当半径为0.1时，高斯模板在计算 后只有中间像素为1,其他像素均为0 （实际上只是趋近0）,即图像不会有变化。第二个疑问，高斯模板大小和高斯半径的关系？这是一个一直困扰我们的误解。因为我们的 思维进入了物理实现的误区。在物理实现中，高斯模板有界，从而使我们忽略了这个问题的 真正答案：高斯模板在逻辑上是无边界的。也就是说高斯模板本质上是逻辑上无穷拓展曲面 的一个近似。因此，模板大小我们应该认为它是无穷大的。只不过在计算的时候，因为在远 处趋近0,因此在某个阈值之下我们不再考虑这些值，这个阈值就是模板

4、边界。下面为二维高斯曲面的公式（x,y代表像素的模板坐标，模板中心位置为原点）：即：g(x,y) = ( 1 / (2*pi*o人2) ) * exp( -(x2+y2)/(2*o人2);根据这个公式，我们可以计算出不同半径下的高斯模板，实际上模板是无穷大的，只是在中 心较远处，他们会趋近0。例如，我们计算出r=0.7时的一个归一化后的高斯模板：/guass radius=0.7000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000000.000000 0.000000 0.0000

5、01 0.000012 0.000033 0.000012 0.000001 0.000000 0.0000000.000000 0.000001 0.000093 0.001976 0.005481 0.001976 0.000093 0.000001 0.0000000.000000 0.000012 0.001976 0.042189 0.117046 0.042189 0.001976 0.000012 0.0000000.000000 0.000033 0.005481 0.117046 0.324724 0.117046 0.005481 0.000033 0.0000000.00

6、0000 0.000012 0.001976 0.042189 0.117046 0.042189 0.001976 0.000012 0.0000000.000000 0.000001 0.000093 0.001976 0.005481 0.001976 0.000093 0.000001 0.0000000.000000 0.000000 0.000001 0.000012 0.000033 0.000012 0.000001 0.000000 0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0

7、.000000 0.000000在网络上众所周知流传的高斯3*3模板实际上是对高斯曲面的一个整数除法形式的近似：1 2 12 4 2 /161 2 1实际验证，我们发现这个3*3模板实际上是对高斯半径约为0.849时的一个近似，当r=0.849 时，其 3*3 归一模板为(在 MATLAB 中，输入h=fspecial(gaussian： 3, 0.849)即可得到这个模板)：(guass radius=0.849000)0.062467 0.125000 0.0624670.125000 0.250131 0.1250000.062467 0.125000 0.062467然后我们可以用Ma

8、tlab中的imfilter来对图像进行高斯模糊的处理：img = imread(c:demo.bmp);h = fspecial(gaussian； 3, 0.849);img2 = imfilter(img, h);subplot(121), imshow(img); title(原 图)subplot(122), imshow(img2); title(高斯模糊后)效果如下：我们可以在Matlab中用如下语句绘制高斯曲面:%绘制高斯模糊曲面！r=0.849; %高斯半径(从0.1到250)x=-3:0.2:3;y=x；X,Y = meshgrid(x,y);Z=exp(-(X.人2+Y.

9、人2)./(2*r*r)/(2*pi*r*r); mesh(X,Y,Z)效果如下图：在数字信号处理中，高斯模糊算法是一种滤波器，它的时域和频域曲线如下所示:1015 Z3 S 30Tim亡 dofYidn020.+口启DBMonndlKd Frequency 声町Frequency 4 and n2DO-2D-+O益B出3亡邑3Leakaae Factor: Relative sidelcbe attenuation: Mainlobe width f-3dBl:从频域曲线看出，高斯模糊本质上一种低通滤波器。体现在图像处理上，图像的边缘等灰度 变化剧烈的地方对应高频信息，将被滤除。最后，我们给

10、出计算高斯模糊模板的C语言代码，请注意，由于高斯模板是对称的，实际 上我们只需要计算出大约1/4模板即可。但这里的代码我们没有做这样的优化。我们输出 的模板是(2*N + 1)*(2N + 1)，高斯半径用r表示。#inelude #inelude #define N 3 /* 模板大小：(2N + 1) * (2N + 1) */void main()double a2*N + 12*N + 1; /* 高斯模板 */double r=0.6; /* 高斯半径：0.1, 250 */double A=1/(2*M_PI*r*r);int i,j;for(i = _1*N;i = N;i+)f

11、or(j = _1*N;j = N;j +)ai+Nj + N=A*exp(-1)*(i*i+j*j)/ (2*r*r);同时，作为比较，我们给出在Matlab中生成高斯模板的代码(fspecial的代码局部，其中 p3是第三个参数即高斯半径)：case gaussian % Gaussian filtersiz = (p2-1 )/2; %注：p2即模板边长，默认值为33std = p3;注：p3即高斯半径，默认为为0.5x,y = meshgrid(-siz 2):siz(2),-siz(1):siz(1);arg = -(x.*x + y.*y)/2*std*std);h = exp(arg);h(heps*max(h(:) = 0;sumh = sum(h(:);注：模板归一化if sumh = 0 /h = h/sumh;end;