高斯分布

filter)包含很多种包括低通、高通、带通等在图像上说的高斯滤波通常是指的高斯模糊(Gaussian

Blur)是一种高斯低通滤波。

通常这个算法也可以用来模糊图像提供模糊滤镜。

也可以用来过滤自然界的高斯白噪声。

高斯分布正态分布是一个常见的连续概率分布正态分布的数学期望值或期望值

\mu

σ等于尺度参数决定了分布的幅度。

正态分布的概率密度函数曲线呈钟形所以又被称为钟形曲线。

我们常说的标准正态分布是位置参数

\mu0,方差\sigma^21

高斯滤波器是一种根据高斯函数的形状来选择权值的线性平滑滤波器对于抑制服从正态分布的噪声非常有效一维零均值高斯函数为

(2)

g(x,y)\frac{1}{2\pi\sigma^2}e^{\frac{-x^2y^2}{2\sigma^2}}\tag{3}

2.高斯滤波性质

高斯函数具有五个重要的性质这些性质使得它在早期图像处理中特别有用这些性质表明高斯平滑滤波器无论在空间域还是在频率域都是十分有效的低通滤波器且在实际图像处理中得到了工程人员的有效使用高斯函数具有五个十分重要的性质它们是

1.二维高斯函数具有旋转对称性即滤波器在各个方向上的平滑程度是相同的一般来说一幅图像的边缘方向是事先不知道的因此在滤波前是无法确定一个方向上比另一方向上需要更多的平滑旋转对称性意味着高斯平滑滤波器在后续边缘检测中不会偏向任一方向

2.高斯函数是单值函数这表明高斯滤波器用像素邻域的加权均值来代替该点的像素值而每一邻域像素点权值是随该点与中心点的距离单调增减的这一性质是很重要的因为边缘是一种图像局部特征如果平滑运算对离算子中心很远的像素点仍然有很大作用则平滑运算会使图像失真

3.高斯函数的傅立叶变换频谱是单瓣的正如下面所示这一性质是高斯函数付立叶变换等于高斯函数本身这一事实的直接推论图像常被不希望的高频信号所污染(噪声和细纹理)而所希望的图像特征如边缘既含有低频分量又含有高频分量高斯函数付立叶变换的单瓣意味着平滑图像不会被不需要的高频信号所污染同时保留了大部分所需信号

4.高斯滤波器宽度(决定着平滑程度)是由参数σ表征的而且σ和平滑程度的关系是非常简单的σ越大高斯滤波器的频带就越宽平滑程度就越好通过调节平滑程度参数σ可在图像特征过分模糊(过平滑)与平滑图像中由于噪声和细纹理所引起的过多的不希望突变量(欠平滑)之间取得折衷

5.由于高斯函数的可分离性较大尺寸的高斯滤波器可以得以有效地实现二维高斯函数卷积可以分两步来进行首先将图像与一维高斯函数进行卷积然后将卷积结果与方向垂直的相同一维高斯函数卷积因此二维高斯滤波的计算量随滤波模板宽度成线性增长而不是成平方增长

要模糊一张图像可以直接用均值滤波来做简单的模糊但是这样做显然不大合理因为图像是连续的离卷积核中心的点关系更加密切越远的点关系越疏远这个时候就需要加权平均。

明显的离中心点越近的像素点权重越大。

而正态分布显然是一种可取的权重分配方式又由于图像是二维的所以需要使用二维的高斯函数。

所以高斯滤波的本质是利用高斯函数来生成高斯核高斯卷积模板来对图像进行卷积操作。

理论上高斯分布在所有定义域上都有非负值这就需要一个无限大的卷积核但是实际上仅需要取均值的三倍标准差即

3\sigma

高斯滤波最重要的就是找到高斯模板然后进行卷积以3X3高斯模板为例假设中心点的坐标为00根据二维高斯函数

g(x,y)

σ0.8(这个值不宜过大否则就会变成均值滤波)可以根据坐标值来算出对应的高斯模板。

\sigma

为了防止高斯滤波后的图像偏亮或者偏暗我们还需要对图像进行归一化这九个点的权重和等于0.9125991因此需要分别对这九个数除以0.9125991最终得到的高斯模板为

有了高斯模板就可以对图像进行卷积了但是在FPGA中对于小数的运算不友好于是我们将这个3X3模板扩大了16倍得到了近似的整数模板卷积完成后再除以16来做定浮点数的近似计算。

16倍后的模板高斯卷积模板为

clear;

exp(-(11)/(2*sigma*sigma))/(2*pi*sigma*sigma);

exp(-(10)/(2*sigma*sigma))/(2*pi*sigma*sigma);

exp(-(00)/(2*sigma*sigma))/(2*pi*sigma*sigma);

A*4

floor(gauss_normal/gauss_normal(1,1));生成卷积模板后的FPGA实现这里就不作过多介绍了参考前几章的3X3卷积核以及均值滤波的实现方式改变参数即可。