图像处理和傅里叶变换原理和运用.docVIP

图像处理与傅里叶变换 1背景 傅里叶变换是一个非常复杂的理论，我们在图像处理中集中关注于其傅里叶离散变换离散傅立叶变换(Discrete Fourier Transform) 1.1离散傅立叶变换 图象是由灰度（RGB）组成的二维离散数据矩阵，则对它进行傅立叶变换是离散的傅立叶变换。 对图像数据f(x,y)（x=0,1,… ,Ｍ-1； y=0,1,… ,N-1）。则其离散傅立叶变换定义可表示为： 式中，u=0,1,…, M-1;v= 0,1,…, N-1 其逆变换为 式中，x=0,1,…, M-1;y= 0,1,…, N-1 在图象处理中，一般总是选择方形数据，即Ｍ＝Ｎ 影像f(x,y)的振幅谱或傅立叶频谱： 相位谱： 能量谱（功率谱） 1.2快速傅里叶变化 可分离性的优点是二维的傅立叶变换或逆变换由两个连续的一维傅立叶变换变换来实现，对于一个影像f(x,y)，可以先沿着其每一列求一维傅立叶变换，再对其每一行再求一维变换 正变化 逆变换 由于二维的傅立叶变换具有可分离性，故只讨论一维快速傅立叶变换。 正变换 逆变换 由于计算机进行运算的时间主要取决于所用的乘法的次数。 按照上式进行一维离散由空间域向频率域傅立叶变换时，对于N个F（u)值，中的每一个都要进行N次运算，运算时间与N2成正比。 1965年库里-图基（ Cooly-Tudey)提出将运算操作降到Nlog2N数量级的算法，即N可以分解为一些较小整数的乘积，当N为2的幂（即N=2P，其中P是整数时），效率最高，实现起来也最简单。这就是快速傅立叶变换。 1.3关于基图像（频率矩形） 由二维离散傅里叶反变换式。可知，由于u和v均有0,1,…,N-1的N个可能的取值，所以f(x,y)由N2个频率分量组成，所以每个频率分量都与一个特定的(u,v)值相对应；且对于某个特定的(u,v)值来说，当(x,y)取遍所有可能的值(x=0，1，…，N-1；y=0，1，…，N-1）时，就可得到对应于该特定的(u,v)值的一幅基图像。基图像可表示为。 所以，一幅图像的灰度平均值可由DFT在原点处的值求得 证明 周期性与共轭对称性 F(u,v)=F(u+mM,v+Nn) F(u,v)=F(-u,-v) 也就是说，图3.7的频谱图(a)和(b)实质上是函数的傅里叶频谱图。 1.3运用 1.3.1频率滤波 在频域中，频率越大说明原始信号变化速度越快；频率越小说明原始信号越平缓。当频率为0时，表示直流信号，没有变化。因此，频率的大小反应了信号的变化快慢。高频分量解释信号的突变部分，而低频分量决定信号的整体形象。 在图像处理中，频域反应了图像在空域灰度变化剧烈程度，也就是图像灰度的变化速度，也就是图像的梯度大小。对图像而言，图像的边缘部分是突变部分，变化较快，因此反应在频域上是高频分量；图像的噪声大部分情况下是高频部分；图像平缓变化部分则为低频分量。也就是说，傅立叶变换提供另外一个角度来观察图像，可以将图像从灰度分布转化到频率分布上来观察图像的特征。书面一点说就是，傅里叶变换提供了一条从空域到频率自由转换的途径。对图像处理而言，以下概念非常的重要： 图像高频分量：图像突变部分；在某些情况下指图像边缘信息，某些情况下指噪声，更多是两者的混合； 低频分量：图像变化平缓的部分，也就是图像轮廓信息 高通滤波器：让图像使低频分量抑制，高频分量通过 低通滤波器：与高通相反，让图像使高频分量抑制，低频分量通过 带通滤波器：使图像在某一部分的频率信息通过，其他过低或过高都抑制 模板运算与卷积定理 在时域内做模板运算，实际上就是对图像进行卷积。模板运算是图像处理一个很重要的处理过程，很多图像处理过程，比如增强/去噪（这两个分不清楚），边缘检测中普遍用到。根据卷积定理，时域卷积等价与频域乘积。因此，在时域内对图像做模板运算就等效于在频域内对图像做滤波处理。 比如说一个均值模板，其频域响应为一个低通滤波器；在时域内对图像作均值滤波就等效于在频域内对图像用均值模板的频域响应对图像的频域响应作一个低通滤波。 图像去噪 图像去噪就是压制图像的噪音部分。因此，如果噪音是高频额，从频域的角度来看，就是需要用一个低通滤波器对图像进行处理。通过低通滤波器可以抑制图像的高频分量。但是这种情况下常常会造成边缘信息的抑制。常见的去噪模板有均值模板，高斯模板等。这两种滤波器都是在局部区域抑制图像的高频分量，模糊图像边缘的同时也抑制了噪声。还有一种非线性滤波-中值滤波器。中值滤波器对脉冲型噪声有很好的去掉。因为脉冲点都是突变的点，排序以后输出中值，那么那些最大点和最小点就可以去掉了。中值滤波对高斯噪音效果较差。 椒盐噪声：对于椒盐采用中值滤波可以很好的去除。用均值也可以取得一定的效果，但是会引起边缘的模糊。 高斯白噪声：白噪音在整个