数字图像处理 第3章

第三章图像变换3.1概述3.2傅立叶变换和性质3.3其他可分离变换3.4霍特林变换3.1概述

为了有效和快速地对图像进行处理和分析常常需要将原定义在图像空间的图像以某种形式转换到另外一些空间，并利用在这些空间的特有性质方便地进行一定的加工，最后再转换回图像空间以得到所需要的效果。3.1概述图像变换可分离变换统计变换霍特林变换傅里叶变换快速傅里叶变换离散余弦变换沃尔什变换哈达玛变换

3.1概述一、图像变换的引入

1.方法：对图像信息进行变换，使能量保持但重新分配。

2.目的：有利于加工、处理（滤除不必要信息(如噪声)，加强/提取感兴趣的部分或特征）。

3.1概述二、用途

1．提取图像特征（如）：（1）直流分量；（2）目标物边缘：F（u,v）高频分量。2．图像压缩：正交变换能量集中，对集中（小）部分进行编码。3．图像增强：低通滤波，平滑噪声；高通滤波，锐化边缘。3.1概述

图像变换是图像处理和分析技术的基础。在图像处理和分析技术的发展中，傅里叶变换曾经起过并仍起着重要的作用。图像傅立叶变换的物理意义

图像的频率是表征图像中灰度变化剧烈程度的指标，是灰度在平面空间上的梯度。

如：大面积的沙漠在图像中是一片灰度变化缓慢的区域，对应的频率值很低；而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域，对应的频率值较高。

傅立叶变换在实际中有非常明显的物理意义，从物理效果看，傅立叶变换是将图像从空间域转换到频率域，其逆变换是将图像从频率域转换到空间域。换句话说，傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数，傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数

傅立叶变换以前，图像（未压缩的位图）是由对在连续空间（现实空间）上的采样得到一系列点的集合，通常用一个二维矩阵表示空间上各点，记为z=f(x,y)。又因空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就必须由梯度来表示，这样我们才能通过观察图像得知物体在三维空间中的对应关系。为什么要提梯度？因为实际上对图像进行二维傅立叶变换得到频谱图，就是图

像梯度的分布图。当然，频谱图上的各点与图像上各点并不存在一一对应的关系，这一点与是否采取移频处理没有关系。

对频谱移频到原点以后，可以看出图像的频率分布是以原点为圆心，对称分布的。将频谱移频到圆心除了可以清晰地出图像频率分布以外，还有一个好处，它可以分离出有周期性规律的干扰信号，比如正弦干扰。一幅频谱图如果带有正弦干扰，移频到原点上就可以看出，除了中心以外还存在以另一点为中心、对称分布的亮点集合，这个集合就是干扰噪音产生的。这时可以很直观的通过在该位置放置带阻滤波器消除干扰。3.1二维离散傅里叶变换（DFT）尺寸为M×N的离散图像函数的DFT反变换可以通过对F(u,v)求IDFT获得（3.3）（3.4）

DFT变换进行图像处理时有如下特点：（1）直流成分为F(0,0)。（2）幅度谱|F(u,v)|对称于原点。（3）图像f(x,y)平移后，幅度谱不发生变化，仅有相位发生了变化。（3.5）（3.6）

频率域幅值与频率

空间域灰度傅立叶变换举例

傅立叶频谱图上明暗不一的亮点的意义图像上某一点与邻域点差异的强弱，即梯度的大小，也即该点的频率的大小（可以这么理解，图像中的低频部分指低梯度的点，高频部分相反）。傅立叶频谱图上明暗不一的亮点的意义一般来讲，梯度大则该点的亮度强，否则该点亮度弱。这样通过观察傅立叶变换后的频谱图，也叫功率图，我们就可以直观地看出图像的能量分布：如果频谱图中暗的点数更多，那么实际图像是比较柔和的（因为各点与邻域差异都不大，梯度相对较小）；反之，如果频谱图中亮的点数多，那么实际图像一定是尖锐的、边界分明且边界两边像素差异较大的。傅立叶变换在图像处理的重要作用1.图像增强与图像去噪

绝大部分噪音都是图像的高频分量，通过低通滤波器来滤除高频—噪声;边缘也是图像的高频分量，可以通过添加高频分量来增强原始图像的边缘；2.图像分割之边缘检测。提取图像高频分量3.图像特征提取：

形状特征：傅里叶描述子

纹理特征：直接通过傅里叶系数来计算纹理特征

其他特征：将提取的特征值进行傅里叶变换来使特征具有平移、

伸缩、旋转不变性4.图像压缩

可以直接通过傅里叶系数来压缩数据；常用的离散余弦变换是傅立叶变换的实变换；

二维离散傅立叶变换的性质1.线性性质：2.比例性质：3.可分离性：可分离性二维离散傅立叶变换DFT可分离性的基本思想是：二维DFT可分离为两次一维DFT

应用：二维快速傅立叶算法FFT，是通过计算两次一维FFT实现的可分离性傅立变换的可分离性质先进行列变换，然后进行行变换。可分离性可分离性二维的傅立叶变换可以通过两次一维傅立叶变换得到。因此对图像进行傅立叶变换可以先对行（水平方向）进行一维的傅立叶变换，得出的结果再对列（垂直方向）进行一维的傅立叶变换。这样就是使得图像的傅立叶变换实际上把频率分成水平分量和垂直分量，即u分量和v分量。图像是一幅只包括单条水平线的简单图像，显然它只在垂直方向上有灰度跳变（垂直方向是梯度的方向），所以从它的频谱图中只能看到垂直分量。这样就能很好的解释变换前后所出现的垂直现象。4.空间位移：5.频率位移：二维离散傅立叶变换的性质

当图像在频率域时移动时需要用到频率位移性质。频率位移性质图像中心化把图像进行傅立叶变换后，往往要把中心移到u0=v0=N/2的位置上频率位移即将f（x,y）的图像频谱（图像能量集中在低频的4个角，见下图(a)）从原点（0，0）移到中心（N/2,N/2），得到一个完整的频谱，称为频谱中心化(见下图(b))将f(x,y)与一个指数项相乘就相当于把其变换后的频域中心移动到新的位置。F(u,v)与一个指数项相乘就相当于把其反变换后的空域中心移动到新的位置。图像经过二维傅立叶变换后，其变换系数矩阵表明：1、若变换矩阵Fn原点设在中心，其频谱能量集中分布在变换系数短阵的中心附近。若所用的二维傅立叶变换矩阵Fn的原点设在左上角，那么图像信号能量将集中在系数矩阵的四个角上。这是由二维傅立叶变换本身性质决定的。同时也表明一股图像能量集中低频区域。2、变换之后的图像在原点平移之前四角是低频，最亮，平移之后中间部分是低频，最亮，亮度大说明低频的能量大（幅角比较大）

（a）原始图像(b)中心化前的频谱图(c)中心化后的频谱图图3.5图像频谱的中心化变换之后的图像在原点平移之前四角是低频，最亮;平移之后中间部分是低频，最亮，亮度大说明低频的能量大。平均值平均值定义：由傅立叶变换定义：因此，f(x,y)的平均值与傅立叶变换系数的关系为：6.平均值：

如果是一幅图像，在原点的傅里叶变换F(0,0)等于图像的平均灰度级，也称作频率谱的直流成分。

7.离散卷积定理设f(x,y)和g(x,y)是大小分别为A×B和C×D的两个数组，则它们的离散卷积定义为卷积定理

（3.12）（3.13）卷积定理：f(x,y)*g(x,y)<=>F(u,v)G(u,v)f(x,y)g(x,y)<=>F(u,v)*G(u,v)【例3.2】用MATLAB实现图像的傅里叶变换。解：MATLAB程序如下：

A=imread('pout.tif'); %读入图像imshow(A);%显示图像A2=fft2(A); %计算二维傅里叶变换A2=fftshift(A2); %将直流分量移到频谱图的中心 figure,imshow(log(abs(A2)+1),[010]);%显示变换后的频谱图

（a）原始图像（b）图像频谱图3.7傅里叶变换3.2二维离散余弦变换（DCT）

任何实对称函数的傅里叶变换中只含余弦项，余弦变换是傅里叶变换的特例，余弦变换是简化DFT的重要方法。DCT逆变换为：

（3.15）（3.15）二维离散余弦变换离散余弦变换，经常被信号处理和图像处理使用，用于对信号和图像(包括静止图像和运动图像)进行有损数据压缩。这是由于离散余弦变换具有很强的"能量集中"特性:大多数的自然信号(包括声音和图像)的能量都集中在离散余弦变换后的低频部分。例如，在静止图像编码标准JPEG中，在运动图像编码标准MJPEG和MPEG的各个标准中都使用了离散余弦变换。在这些标准制中都使用了二维的离散余弦变换，并将结果进行量化之后进行熵编码。这时对应离散余弦变换中的n是8，并用公式对每个8x8块的每行进行变换，然后每列进行变换。得到的是一个8x8的变换系数矩阵。其中(0,0)位置的元素就是直流分量，矩阵中的其他元素根据其位置表示不同频率的交流分量。

【例3.3】应用MATLAB实现图像的DCT变换。解：MATLAB程序如下：A=imread('pout.tif'); %读入图像I=dct2(A); %对图像作DCT变换 subplot(1,2,1),imshow(A);%显示原图像subplot(1,2,2),imshow(log(abs(I)),[05]);（a）原图（b）DCT系数图3.10离散余弦变换3.3其他可分离变换1923年，美国数学系J.LWalsh提出walsh函数。[1]函数展开有三种：Walsh序的Walsh函数，佩利序的Walsh函数，哈达玛序的Walsh函数。沃尔什变换主要用于图像变换，属于正交变换。这种变换压缩效率低，所以实际使用并不多。但它快速，因为计算只需加减和偶尔的右移操作。沃尔什变换的定义如下：给定一个NXN像素块Pxy（N必须是2的幂），二维WHT定义为：1、沃尔什变换3.3其他可分离变换

沃尔什（Walsh）变换是一种可分离变换。当时，变换核为：

3.3其他可分离变换离散沃尔什变换W(u)为：是z的二进制表达中的第k位。例如n=3，则对z=6（1102），有b0(z)=0，b1(z)=1，b2(z)=1。3.3其他可分离变换由沃尔什变换核组成的矩阵是一个对称矩阵并且其行和列正交（即各行向量与各列向量的内积为0，互相独立），即：

所以离散沃尔什反变换为：

3.3其他可分离变换2-D的沃尔什正变换核和反变换核由以下2式给出：

这2个核完全相同，所以下面2式给出的2-D沃尔什正变换和