查看“︁高斯模糊”︁的源代码

{{noteTA|G1=IT|time=2010-02-28T12:55:10+00:00}}
[[File:Cappadocia Gaussian Blur.svg|thumb|250px|]]
'''高斯模糊'''（英语：Gaussian Blur），也叫高斯平滑，是在{{lang|en|[[Adobe Photoshop]]}}、{{lang|en|[[GIMP]]}}以及{{lang|en|[[Paint.NET]]}}等图像处理软件中广泛使用的处理效果，通常用它来减少图像噪声以及降低细节层次。这种模糊技术生成的图像，其视觉效果就像是经过一个半透明屏幕在观察图像，这与镜头焦外成像效果[[散景]]以及普通照明阴影中的效果都明显不同。高斯平滑也用于[[计算机视觉]]算法中的预先处理阶段，以增强图像在不同比例大小下的图像效果（参见[[尺度空间]]表示以及[[尺度空间]]实现）。
从数学的角度来看，图像的高斯模糊过程就是图像与[[正态分布]]做卷积。由于正态分布又叫作“高斯分布”，所以这项技术就叫作高斯模糊。图像与圆形[[方框模糊]]做卷积将会生成更加精确的焦外成像效果。由于高斯函数的[[傅立叶变换]]是另外一个高斯函数，所以高斯模糊对于图像来说就是一個[[低通滤波器]]。

==机理==

高斯模糊是一种图像模糊[[滤波器]]，它用[[正态分布]]计算图像中每个像素的[[变换]]。N维空间正态分布方程为

:<math>G(r) = \frac{1}{\sqrt{2\pi \sigma^2}^N} e^{-r^2/(2 \sigma^2)}</math>

在二维空间定义为

:<math>G(u,v) = \frac{1}{2\pi \sigma^2} e^{-(u^2 + v^2)/(2 \sigma^2)}</math>

其中''r''是模糊半径 （<math>r^2 = u^2 + v^2</math>），σ是正态分布的[[标准偏差]]。在二维空间中，这个公式生成的曲面的[[等高线]]是从中心开始呈正态分布的[[同心圆]]。分布不为零的像素组成的[[卷积]]矩阵与原始图像做变换。每个像素的值都是周围相邻像素值的[[加权平均]]。原始像素的值有最大的高斯分布值，所以有最大的权重，相邻像素随着距离原始像素越来越远，其权重也越来越小。这样进行模糊处理比其它的均衡模糊滤波器更高地保留了边缘效果，参见[[尺度空间实现]]。

理论上来讲，图像中每点的分布都不为零，这也就是说每个像素的计算都需要包含整幅图像。在实际应用中，在计算[[高斯函数]]的离散近似时，在大概3σ距离之外的像素都可以看作不起作用，这些像素的计算也就可以忽略。通常，图像处理程序只需要计算<math> (6 \sigma + 1) \times (6 \sigma + 1) </math>的矩阵就可以保证相关像素影响。对于边界上的点，通常采用复制周围的点到另一面再进行加权平均运算。

除了圆形对称之外，高斯模糊也可以在二维图像上对两个独立的一维空间分别进行计算，这叫作'''線性可分'''。这也就是说，使用二维矩阵变换得到的效果也可以通过在水平方向进行一维高斯矩阵变换加上竖直方向的一维高斯矩阵变换得到。从计算的角度来看，这是一项有用的特性，因为这样只需要<math>O(n \times M \times N) + O(m \times M \times N)</math>次计算，而不可分的矩阵则需要<math>O(m \times n \times M \times N)</math>次计算，其中<math>M</math>,<math>N</math>是需要进行滤波的图像的维数，<math>m</math>、<math>n</math>是滤波器的维数。

对一幅图像进行多次连续高斯模糊的效果与一次更大的高斯模糊可以产生同样的效果，大的高斯模糊的半径是所用多个高斯模糊半径平方和的平方根。例如，使用半径分别为6和8的两次高斯模糊变换得到的效果等同于一次半径为10的高斯模糊效果，<math>\sqrt{6\times6 + 8\times8} = 10</math>。根据这个关系，使用多个连续较小的高斯模糊处理不会比单个高斯较大处理时间要少。

在减小图像尺寸的场合经常使用高斯模糊。在进行欠采样的时候，通常在采样之前对图像进行低通滤波处理。这样就可以保证在采样图像中不会出现虚假的高频信息。高斯模糊有很好的特性，如没有明显的边界，这样就不会在滤波图像中形成震荡。

==高斯矩阵示例==
这是一个计算σ = 0.84089642的高斯分布生成的示例矩阵。注意中心元素 (4,4)处有最大值，随着距离中心越远数值对称地减小。

{| class="wikitable"
|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
|}

注意中心处的0.22508352比3σ外的0.00019117大1177倍。

==例子==

下面的例子展示了高斯模糊的效果；图2是图1经过高斯模糊滤波器（σ = 2）处理之后的结果。

[[File:H3224.GIF|图1:原始图像]]
[[File:H3224020.GIF|图2:经过高斯模糊处理的图像]]

==参见==

*[[无限脉冲响应]]（IIR）
*[[尺度空间实现]]
*[[高斯滤波器]]
*[[高斯]]

==参考文献==
* https://web.archive.org/web/20061109221710/http://www.cee.hw.ac.uk/hipr/html/gsmooth.html

[[Category:图像处理]]
[[Category:计算机图形学]]
[[Category:高斯函數]]