第三章 图像变换.ppt

1、数字图像处理 (Digital Image Processing）,湖北师范学院教育信息与技术学院,第三章 图像变换,3.1 引言 3.2 空域变换 3.3 图像的频域变换 3.4 离散傅立叶变换 3.5 离散余弦变换 3.6 KL变换 3.7 其他正交变换,3.1 引言,图像的数学变换的特点在于其有精确的数学背景，是许多图像处理技 术的基础。在这些变换中，一种是在空间域上进行的，这些变换根据处理 操作的特点，可以分为图像的代数运算和几何运算，它们都是利用对输入 图像进行加工而得到输出图像。另一种重要的数学变换则是将原定义在图 像空间的图像以某种形式转换到另外一些空间，并利用输入图像在这些空

2、间的特有性质有效而快速地对图像进行处理和分析。最典型的变换有离散 傅立叶变换，它把空域中的图像信号看作二维时间序列，将其变换到频率 域来分析图像的频谱特性。 空域变换有如加法、减法等的代数变换，也有如旋转、拉伸等的几何变 换；频域变换除了傅立叶变换外，常用的非空域的变换还有离散余弦变换等。无论是在空域中的数学变换还是频域中的数学变换，它们在图像分析、滤波、增强、压缩等处理中都有着非常典型而重要的应用。,3.2 空域变换,1. 代数变换 图像的代数变换是指对两幅图像进行点对点的四则运算而得到一幅新的 输出图像。图像的代数运算在图像处理中有着广泛的应用，它除了可以实现 自身所需要的算术操作，还能为

3、许多复杂的图像处理提供准备。 1）加法运算 2）减法运算（差分）,+,=,=,3）其他灰度变换 灰度变换的变换函数曲线如下图：,(1) 图象求反,(2) 对比度拉伸,(3) 动态范围压缩,2. 几何变换 图像在生成过程中，由于系统本身具有非线性或拍摄角度不同，会使生 成的图像产生几何失真。几何失真一般分为系统失真和非系统失真。系统失 真是有规律的、能预测的；非系统失真则是随机的。 几何变换可以改变图像中物体之间的空间关系。这种运算可以看成是图 像内的各物体在图像内移动的过程。例如，物体的转动、扭曲、倾斜、拉伸 等等，都是几何运算的结果。如下图：,1）平移,2）放缩,3）旋转,4）水平镜像,5）

4、垂直镜像,6）一般的几何变换 下图显示了在失真和相应的校正图像中的四边形区域，四边的顶点是相 应的“控制点”。假设四边形区域中的几何形变过程用双线性方程对来建模， 即：,7）几何变换中灰度插值 对图像作定量分析时，就要对失真的图像进行几何校正（即将存在几 何失真的图像校正成无几何失真的图像），以免影响分析精度。基本方 法是先建立几何校正的数学模型；其次利用已知条件确定模型参数；最后 根据模型对图像进行几何校正。通常分为两步： (1) 图像空间的坐标变换； (2) 确定校正空间各象素的灰度值。 输出象素通常被映射到输入图像中的非整数位置，即位于四个输入象 素之间。因此，为了决定与该位置相对应的灰

5、度值，必须进行插值运算。 常用的插值方法有3种： (1) 最近邻插值(Nearest Neighbor Interpolation) (2) 双线性插值(Bilinear Interpolation) (3) 三次立方插值,(1) 最近邻插值(Nearest Neighbor Interpolation) 最简单的插值方法是最近邻插值，即选择离它所映射到的位置最近的 输入象素的灰度值为插值结果。数学表示为： (2) 双线性插值(Bilinear Interpolation) 双线性插值法是对最近邻法的一种改进，即用线性内插方法，根据点 的四个相邻点的灰度值，分别在x和y方向上进行两次线性插值插

6、值。如下 图：,首先，在x方向上作线性插值，对上端的两个顶尖进行线性插值得：,类似的，对于底端两个顶点进行线性插值有：,y方向上作线性插值，以确定：,最后得到双线性插值公式为：,(3) 三次立方插值 该方法利用三次多项式 来逼近理论上的最佳插值函数 ，其 数学表达式为： 上式中的是周围象素沿方向离原点的距离。待求象素的灰度值由其周 围16个点的灰度值加权内插得到。可推导出待求象素的灰度值计算式为：,其中：,3.3 图像频域变换,线性系统 系统的定义：接受一个输入，并产生相应输出的任何实体。系统的输入是一个或两个变量的函数，输出是相同变量的另一个函数。,线性系统定义 对于特定的系统，有： X1(

7、t) y1(t) X2(t) y2(t) 该系统是线性的当且仅当：,1.图像的频域变换的理论基础,x1(t) + x2(t) y1(t) + y2(t) 从而有： ax1(t) ay1(t),线性系统移不变性的定义 对于某线性系统，有： x(t) y(t) 当输入信号沿时间轴平移T，有： x(t - T) y(t - T) 则称该线性系统具有移不变性。,卷积的定义 对于一个线性系统的输入f(t)和输出y(t)，其间必定存在关系： ,h(t) 称为线性系统的单位冲激响应函数，其含义为：当线性系统输入f(t)为单位脉冲函数时，线性系统的输出响应。 ,上式称之为卷积积分。,脉冲函数的极限定义：脉冲函

8、数可以看成是一系列函数的极限，这些函数的振幅逐渐增大，持续时间逐渐减少，而保持面积不变。,脉冲函数的定义(也叫函数):,离散一维卷积,二维卷积的定义,离散二维卷积,相关的定义 任意两个信号的相关函数定义为： ,相关与卷积的关系：,2.正交变换,正交变换 连续函数集合的正交性,当C = 1时，称集合为归一化正交函数集合,即每 一个向量为单位向量 其物理意义为多维空间坐标的基轴方向互相正交。,正交函数集合的完备性,若f(x)是定义在t0和t0+T 区间的实值信号，平方可积。可以表示为：,对任意小的0，存在充分大的N，用N个有限展开式估计f(x)时：,可有：,则称函数U 集合是完备的。,正交函数集合

9、完备性的物理意义 任何数量的奇函数累加仍为奇函数 任何数量的偶函数累加仍为偶函数,因此为了能用累加展开式来表示一个任意函数，就要求这个函数集合中既有奇函数又有偶函数,正交函数集合完备性图例,(a) 完备 (b) 不完备,正交函数的离散情况,N 个正交向量,当C1时，成为归一正交化,正交函数的离散情况,N 个正交向量矩阵,必满足：,一维正交变换,对于一维向量f，用上述正交矩阵进行运算：,若要恢复f ，则,以上过程称为正交变换。,一般范式酉变换,若A为复数矩阵，正交的条件为：,其中A*为A的复数共轭矩阵，满足这个条件的矩阵为酉矩阵(unitary matrix)。对于任意向量f的运算称为酉变换(u

10、nitary transform)：,二维酉变换,NN 二维函数可以类似于一维用正交序列展开和恢复。,变换核的可分离性,为一维完备正交基向量的集合。用矩阵表示：,通常选择A =B 。,二维酉变换，A=B时，二维酉变换正变换表示为：,用矩阵表示：,类似的，对于MN 的二维函数f(x,y),基图像,对反变换,可看成是基图像,权因子,酉变换的性质,1. 酉矩阵是正交阵,4.酉变换能量的紧缩 正交酉变换往往趋于将信号能量压缩到相,对少的变换系数中，由于总能量保持不变，因此许多变换系数将包含很少的能量,KL变换可以达到最大的能量紧缩。,5.酉变换去相关 当输入向量元素间高度相关时，变换系数趋向于去相关，

11、这意味着协方差矩阵的非对角项和对角项相比趋于变小。,KL变换可以达到完全的去相关。,则F(u,v)的均值为：,F(u,v)的协方差为：,7. 其他性质： (1) A为酉阵，则其行列式值|A|=1 (2) 若a为向量，则作酉变换后向量模保持不变： b=Aa，则|b|=|a|。,将图像看成是线性叠加系统； 图像在空域上具有很强的相关性； 图像变换是将图像从空域变换到其它域如频域的数变换； 借助于正交变换的特性可使在空域上的复杂计算转换到频域后得到简化； 借助于频域特性的分析，将更有利于获得图像的各种特性和进行特殊处理。,图像变换定义,可进行图像变换的基本条件,满足正交、完备两个条件的函数集合或矩阵

12、才能用于图象的分析。,常用的几种变换：傅里叶变换、WALSH变换、哈达玛变换、K-L变换等，都满足正交性和完备性两个条件。,3.离散图像变换,将离散图象的正交变换为图象信号在一组二维离散完备正交基上的展开,这种正交基展示具有无损重构的性质,以及图象能量的集中和图象信号元素的去相关性能,在图象处理中具有重要的作用。 若离散图象f(m,n)及其在离散完备正交基a(u,v;m,n)上的展开系数为g(u,v)，即：,离散图像的正交变换,1、二维离散完备正交基a(u,v;m,n)的正交性满足,离散图像正交变换的特性,正交性保证变换后图象的紧缩性，图象的去相关性和保证任何被截断的级数展开将使均方误差和为最

13、小。,2、二维离散完备正交基a(u,v;m,n)的完备性满足,完备性保证变换后图象无失真的重构，即保证了当包括了全部系数时，重构误差将为零。,3.4 离散傅立叶变换,调谐信号(欧拉公式)：,傅立叶积分：,其中t代表时间，f代表频率。,基本数学概念,1.连续傅立叶变换,f(x)为连续可积函数，其傅立叶变换定义为：,R(u)，I(u)分别称为傅里叶变换F(u)的实部和虚部。,傅立叶变换的定义(一维),其反变换为：,通常f(x) 的傅里叶变换为复数，可有通用表示式为：,可进一步写为指数形式：,其中：,称之为f(x) 的幅度谱、振幅谱或富里叶谱。,称之为f(x) 的相位谱、相位角。,变换分析的直观说明

14、：,把一个信号的波形分解为许多不同频率正弦波之和。,一维傅立叶变换举例：,方波信号：,经过傅立叶变换后：,几种特殊函数的傅里叶变换：,矩形函数：,矩形函数的傅里叶变换：,sin(x)/x类函数：,sin(x)/x类函数的傅里叶变换：,常数函数：,常数函数的傅里叶变换：,脉冲函数：,脉冲函数的傅里叶变换：,余弦函数：,余弦函数的傅里叶变换：,一维离散傅立叶变换(DFT),一维离散傅立叶变换公式为：,逆变换为：,二维离散傅立叶变换(2DFT),二维傅立叶变换由一维傅立叶变换推广而来：,2）2D傅立叶变换 傅立叶变换可推广到二维函数。如果二维函数f(x,y)满足狄里赫莱条件，那么存在下面的二维傅立叶

15、变换对：,类似于一维傅立叶变换，二维傅立叶变换的幅度谱和相位谱：,2.离散傅立叶变换 如果x(n)为一数字序列，0nN-1,则其离散傅立叶变换定义如下：,其中，u,m均取0,1,M-1；v,n均取0,1,N-1； W1=exp(-j2/M)；W2=exp(-j2/N)）。,二维离散傅立叶变换：,3.离散傅立叶变换的性质 傅立叶变换有许多重其要的性质，这些性质为实际应用提供了诸多便利。下面以二维傅立叶变换为例，介绍几个主要的性质。 1）可分离性 令： 则：,2）线性 傅立叶变换是线性变换，满足线性变换的叠加性：,3）共轭对称性 如果F(u,v)是f(x,y)的傅立叶变换，F*(-u,-v)是傅立

16、叶变换的共轭函数，那么：,4）旋转性 如果空间域函数旋转的角度为0,那么在变换域中此函数的傅立叶变换也旋转同样的角度，即：,5）比例变换性 如果如果是的傅立叶变换，a和b是两个标量，那么：,6）Parseval定理 这个性质也称为能量保持定理。如果F(u,v)是f(x,y)的傅立叶变换，那么有下式成立：,这个性质说明变换前后的能量保持不变。,7）相关定理 两个二维函数f(x,y),g(x,y)的相关函数定义如下：,符号“”表示相关运算。傅立叶变换的一个重要性质是相关定理：,8）卷积定理 两个二维函数f(x,y),g(x,y)的卷积运算定义如下：,符号“*”表示卷积运算。根据上面的定义，傅立叶变

17、换的卷积定理如下：,4. 快速傅立叶变换 1965年，库利图基提出把原始的N点序列依次分解成一系列短序 列，然后求出这些短序列的离散傅立叶变换，以此来减少乘法运算，这 就是快速傅立叶变换（Fast Fourier Transform,FFT）。 对于一个有限长序列x(n)，0nN-1, ,按n的奇偶把x(n) 分解为两个N/2点的子序列: 其傅立叶变换为：,前一半的值,后一半的值,复乘：,复加：,每个 点DFT 分解成两个 点DFT。,第一级,第二级,第三级,上面称为DIT-FFT算法。也可以在频域进行分解，相应算法是DIF-FFT，其蝶形分解图如下：,5. 离散傅立叶变换的显示 下图是一个图

18、的DFT频谱图：,离散傅立叶变换的显示对称平移后,其他几个图像的频谱图：,注意观察对应关系,Jean Buptiste Joseph Fourier和他的付立叶变换,(a)输入图像 (b)幅值谱 (c)相位谱 (d)由幅值谱重构的图象 (e)由相位谱重构的图象,结论：相位谱可能具有更重要的应用,3.5 离散余弦变换,1D离散余弦变换的正变换 余弦变换为： 其中， 反变换为：,2. 2D离散余弦变换的正变换 正变换为： 其中， a(v)与a(u)定义类似。 反变换为：,3. 2D DCT基图像,二维DCT基图像,DCT变换举例：,3.6 KL变换,1. 离散KL变换 KL（Karhunen-Lo

19、eve）或（DKT），也称为Hotelling变换、特征向 量变换(Eigenvector-Based Transform)、主分量（PCA）变换等。 它是一种利用图像的统计性质/统计模型的变换。常用在数据压缩、特 征提取等方面。 1）定义 设图像f NN (x,y), 采样了M次，得到集合 fi(x,y),i=1,2,M。对每个 图像按行或按列依次排列，得采样图像：Xi=(xi1,xi2,xiN2)T, 则有： Mx为均值向量： Mx=EX， CX为协方差矩阵：CX=E(X-Mx)(X-Mx)T,近似表示：,令ei,i, i=1,2,N2分别表示矩阵的特征矢量与特征值，将i减序排列，1 i2

20、N2, 构造变换矩阵：,KL变换：,KL反变换：,变换后均值为0，方差为：,2）性质,2. 应用 1）压缩 CX实对称矩阵，总可以找到标准正交的特征向量集合构成A， A-1=A，由Y利用 X=AY+MX 重建X。 压缩时，取k个大的i，并构造出Ak,，则 离散K-L变换在最小平方误差的意义上最优。,特点：1）比其它方法图像压缩的效率高； 2）图像标准化（旋转）。 其缺点是：1）非分离，需要计算CX ，及其特征值、特征向量； 2）无快速算法。,2）图像旋转 第一基向量与数据中最大变化的方向相对应。若目标已抽出，希望 与某个标准的、或不变的方向对准。需要处理目标中各像素的坐标。如 下图所示：,二维目标的旋转。(a)原始数据的散布指明单位特征向量的方向； (b)利用变换Y=AX作数据旋转；(c)利用变换Y=A(X-MX)作数据旋转和