数字图像处理与模式识别课件.ppt

《数字图像处理与模式识别课件.ppt》由会员分享，可在线阅读，更多相关《数字图像处理与模式识别课件.ppt（88页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数字图像处理与模式识别第1页，此课件共88页哦第二章 图像处理中的常用数学变换2.1 2.1 引言引言2.2 2.2 空域变换空域变换 2.2.1 2.2.1 代数运算代数运算 2.2.2 2.2.2 几何运算几何运算2.3 2.3 离散傅立叶变换离散傅立叶变换 2.3.1 2.3.1 离散傅立叶变换基本概念离散傅立叶变换基本概念 2.3.2 2.3.2 离散傅立叶变换基本性质离散傅立叶变换基本性质 2.3.3 2.3.3 快速离散傅立叶变换快速离散傅立叶变换2.4 2.4 离散离散GaborGabor变换变换 2.4.1 2.4.1 加窗傅立叶变换加窗傅立叶变换 2.4.2 Gabor2.4

2、.2 Gabor变换的基本概念变换的基本概念 2.4.3 2.4.3 离散离散GaborGabor变换变换2.5 2.5 小波变换小波变换 2.5.1 2.5.1 连续小波变换连续小波变换 2.5.2 2.5.2 二进小波变换二进小波变换 2.5.3 2.5.3 离散小波变换离散小波变换 2.5.4 2.5.4 二维离散小波变换二维离散小波变换 2.5.5 2.5.5 小波变换的应用小波变换的应用2.6 PCA2.6 PCA变换变换 2.6.1 PCA2.6.1 PCA的基本概念及问题描述的基本概念及问题描述 2.6.2 PCA 2.6.2 PCA变换的应用变换的应用2.72.7离散余弦变换离

3、散余弦变换2.82.8其他的正交变换其他的正交变换第2页，此课件共88页哦2.1 引言 图像的数学变换的特点在于其有精确的数学背景，是许多图像的数学变换的特点在于其有精确的数学背景，是许多图像处理技术的基础。在这些变换中，一种是在空间域上进行图像处理技术的基础。在这些变换中，一种是在空间域上进行的，这些变换根据处理操作的特点，可以分为图像的代数运算的，这些变换根据处理操作的特点，可以分为图像的代数运算和几何运算，它们都是利用对输入图像进行加工而得到输出图和几何运算，它们都是利用对输入图像进行加工而得到输出图像。另一种重要的数学变换则是将原定义在图像空间的图像以像。另一种重要的数学变换则是将原定

4、义在图像空间的图像以某种形式转换到另外一些空间，并利用输入图像在这些空间的某种形式转换到另外一些空间，并利用输入图像在这些空间的特有性质有效而快速地对图像进行处理和分析。最典型的变换特有性质有效而快速地对图像进行处理和分析。最典型的变换有离散傅立叶变换，它把空域中的图像信号看作二维时间序列，有离散傅立叶变换，它把空域中的图像信号看作二维时间序列，将其变换到频率域来分析图像的频谱特性。将其变换到频率域来分析图像的频谱特性。除了傅立叶变换外，常用的非空域的变换还有除了傅立叶变换外，常用的非空域的变换还有GaborGabor变换、变换、小波变换、离散余弦变换、小波变换、离散余弦变换、PCAPCA变换

5、等等。无论是在空域中的变换等等。无论是在空域中的数学变换还是频域中的数学变换，它们在图像分析、滤波、增数学变换还是频域中的数学变换，它们在图像分析、滤波、增强、压缩等处理中都有着非常典型而重要的应用。强、压缩等处理中都有着非常典型而重要的应用。第3页，此课件共88页哦2.2 空域变换 2.2.1 2.2.1 代数运算代数运算 图像的代数运算是指对两幅图像进行点对点的四则运算而得到一幅图像的代数运算是指对两幅图像进行点对点的四则运算而得到一幅新的输出图像。图像的代数运算在图像处理中有着广泛的应用，它除了新的输出图像。图像的代数运算在图像处理中有着广泛的应用，它除了可以实现自身所需要的算术操作，还

6、能为许多复杂的图像处理提供准备。可以实现自身所需要的算术操作，还能为许多复杂的图像处理提供准备。1.1.加法运算加法运算加法运算加法运算 2.2.减法运算（差分）减法运算（差分）减法运算（差分）减法运算（差分）第4页，此课件共88页哦+=第5页，此课件共88页哦=第6页，此课件共88页哦（a）原图 （b）梯度运算第7页，此课件共88页哦2.2.2 2.2.2 几何运算几何运算 几何运算可以改变图像中物体之间的空间关系。这几何运算可以改变图像中物体之间的空间关系。这种运算可以看成是图像内的各物体在图像内移动的过程。种运算可以看成是图像内的各物体在图像内移动的过程。例如，物体的转动、扭曲、倾斜、拉

7、伸等等，都是几何例如，物体的转动、扭曲、倾斜、拉伸等等，都是几何运算的结果。运算的结果。第8页，此课件共88页哦0,0 xy 旋转第9页，此课件共88页哦0,0 xy 水平镜像第10页，此课件共88页哦0,0 xy 垂直镜像第11页，此课件共88页哦平移放缩第12页，此课件共88页哦旋转第13页，此课件共88页哦复杂变换 右图显示了在失真和相应的校正图像中的四边形区域，四边的右图显示了在失真和相应的校正图像中的四边形区域，四边的顶点是相应的顶点是相应的“控制点控制点”。假设四边形区域中的几何形变过程用双线性。假设四边形区域中的几何形变过程用双线性方程对来建模，即：方程对来建模，即：FDCBAF

8、DCAB第14页，此课件共88页哦几何变换的应用举例几何变换的应用举例几何变换的应用举例几何变换的应用举例 图像在生成过程中，由于系统本身具有非线性或拍摄角度图像在生成过程中，由于系统本身具有非线性或拍摄角度不同，会使生成的图像产生几何失真。几何失真一般分为不同，会使生成的图像产生几何失真。几何失真一般分为系统失真和非系统失真。系统失真是有规律的、能预测的；系统失真和非系统失真。系统失真是有规律的、能预测的；非系统失真则是随机的。非系统失真则是随机的。但对图像作定量分析时，就要对失真的图像进行几何校正（即但对图像作定量分析时，就要对失真的图像进行几何校正（即将存在几何失真的图像校正成无几何失真

9、的图像），以免影响将存在几何失真的图像校正成无几何失真的图像），以免影响分析精度。基本方法是先建立几何校正的数学模型；其次利用分析精度。基本方法是先建立几何校正的数学模型；其次利用已知条件确定模型参数；最后根据模型对图像进行几何校正。已知条件确定模型参数；最后根据模型对图像进行几何校正。通常分为两步：通常分为两步：(1)(1)图像空间的坐标变换；图像空间的坐标变换；(2)(2)确定校正空间各象素的灰度值。确定校正空间各象素的灰度值。第15页，此课件共88页哦灰度级插值 输出象素通常被映射到输入图像中的非整数位置，即位于四输出象素通常被映射到输入图像中的非整数位置，即位于四个输入象素之间。因此，

10、为了决定与该位置相对应的灰度值，个输入象素之间。因此，为了决定与该位置相对应的灰度值，必须进行插值运算。常用的插值方法有必须进行插值运算。常用的插值方法有3 3种：种：1 1）最近邻插值(Nearest Neighbor Interpolation)2 2）双线性插值）双线性插值(Bilinear Interpolation)(Bilinear Interpolation)3）三次立方插值第16页，此课件共88页哦1）最近邻插值(Nearest Neighbor Interpolation)(Nearest Neighbor Interpolation)最简单的插值方法是最近邻插值，即选择离它

11、所映射到的位置最最简单的插值方法是最近邻插值，即选择离它所映射到的位置最近的输入象素的灰度值为插值结果。数学表示为：近的输入象素的灰度值为插值结果。数学表示为：2）双线性插值(Bilinear Interpolation)双线性插值法是对最近邻法的一种改进，即用线性内插方法，根双线性插值法是对最近邻法的一种改进，即用线性内插方法，根据点的四个相邻点的灰度值，分别在据点的四个相邻点的灰度值，分别在x x和和y y方向上进行两次插值，方向上进行两次插值，计算出的值。最后形成的插值函数为一双曲抛物面方程：计算出的值。最后形成的插值函数为一双曲抛物面方程：第17页，此课件共88页哦(0,0)f(0,0

12、)(x,0)(0,y)(0,1)(x,1)(1,1)(1,0)f(1,0)(x,y)f(x,y)灰度双线性插值示意图yx第18页，此课件共88页哦首先，在x x方向上作线性插值，对上端的两个顶尖进行方向上作线性插值，对上端的两个顶尖进行线性插值得：线性插值得：类似的，对于底端两个顶点进行线性插值有：y方向上作线性插值，以确定：最后得到双线性插值插值公式为：第19页，此课件共88页哦3 3）三次立方插值 该方法利用三次多项式该方法利用三次多项式 来逼近理论上的最佳插值函数来逼近理论上的最佳插值函数 ，其数学表达式为：其数学表达式为：上式中的是周围象素沿方向离原点的距离。待求象素的灰度值由上式中的

13、是周围象素沿方向离原点的距离。待求象素的灰度值由其周围其周围1616个点的灰度值加权内插得到。可推导出待求象素的灰度个点的灰度值加权内插得到。可推导出待求象素的灰度值计算式为：值计算式为：第20页，此课件共88页哦.2 .1 0 1 2S(x)x三次立方插值原理图0uv(x,y)(i,j)(i+1,j)(i+1,j+1)(i,j+1)(i.1,j.1)(i.1,j+2)(i+2,j.1)(i+2,j+2)第21页，此课件共88页哦其中：第22页，此课件共88页哦2.3 离散傅立叶变换2.3.1 傅立叶定义理论基础、连续与离散的傅立叶变换。2.3.2 2.3.2 二维傅立叶变换特性二维傅立叶变换

14、特性可分离性、周期与共轭对称、平移性；旋转特性、线性与相似性、均值性；拉普拉斯、卷积与相关。2.3.3 快速傅立叶变换FFT算法、逆向FFT算法、算法实现。第23页，此课件共88页哦连续与离散的傅立叶变换 一维连续傅立叶变换 二维连续傅立叶变换 离散傅立叶变换 离散傅立叶变换的计算与显示3.1 傅立叶变换理论基础第24页，此课件共88页哦2.3.1 傅立叶变换导言:傅立叶变换离散傅立叶变换的计算与显示离散傅立叶变换的计算举例离散傅立叶变换的显示第25页，此课件共88页哦离散傅立叶变换的计算举例xf(x0)=f(x0+x)01231234第26页，此课件共88页哦2.3.1 傅立叶变换导言:傅立

15、叶变换F(0)=1/4f(x)exp0=1/4f(0)+f1(1)+f(2)+f(3)=1/4(2+3+4+4)=3.25F(1)=1/4f(x)exp-j2x/4)=1/4(2e0+3e j21 1/4+4e j22/4+4e j23/4)=1/4(-2+j)F(2)=-1/4(1+j0)F(3)=-1/4(2+j)第27页，此课件共88页哦离散傅立叶变换的显示 通过对傅立叶变换模，来显示傅立叶变换图象。由于模的值域大于显示的值域，因此要进行动态值域的压缩D(u,v)=c log(1+|F(u,v)|)其中：c=255/k;k=max(log(1+|F(u,v)|)值域0,k的上限（最大值）

16、第28页，此课件共88页哦离散傅立叶变换的显示第29页，此课件共88页哦第30页，此课件共88页哦离散傅立叶变换的显示对称平移后第31页，此课件共88页哦第32页，此课件共88页哦2.3.2 二维傅立叶变换特性可分离性周期与共轭对称平移性旋转特性 线性与相似性 均值性 拉普拉斯 卷积与相关第33页，此课件共88页哦2.3.2 二维傅立叶变换特性:可分离性先对行做变换：然后对列进行变换f(x,y)（0，0）（N-1，M-1）xyF(x,v)（0，0）（N-1，M-1）xvF(x,v)（0，0）（N-1，M-1）xvF(u,v)（0，0）（N-1，M-1）uv第34页，此课件共88页哦2.3.3

17、快速傅立叶变换:FFT算法思想分析这些表达式得到如下的特性：分析这些表达式得到如下的特性：（1 1）一个）一个N N个点的变换，能够通过将原始表达 式分成两个部分来计算式分成两个部分来计算（2 2）通过计算两个（）通过计算两个（N/2N/2）个点的变换。得到）个点的变换。得到 F Feveneven(u)(u)和和 F Fodd(u)（3 3）奇部与偶部之和得到）奇部与偶部之和得到F(u)F(u)的前(N/2)(N/2)个值。个值。（4 4）奇部与偶部之差得到）奇部与偶部之差得到F(u)F(u)的后的后(N/2)个值。个值。且不需要额外的变换计算且不需要额外的变换计算。第35页，此课件共88页

18、哦2.3.3 快速傅立叶变换:FFT算法思想快速傅立叶变换的思想：快速傅立叶变换的思想：1 1）通过计算两个单点的）通过计算两个单点的DFTDFT，来计算两个点的，来计算两个点的DFTDFT2 2）通过计算两个双点的）通过计算两个双点的DFTDFT，来计算四个点的，来计算四个点的 DFT DFT，以此类推，以此类推3 3）对于任何）对于任何N=2N=2mm的DFT的计算，通过计算两个N/2 点的DFT，来计算N个点的DFT第36页，此课件共88页哦第37页，此课件共88页哦2.4 离散Gabor变换2.4.1 2.4.1 加窗傅立叶变换2.4.2 Gabor2.4.2 Gabor变换的基本概念

19、变换的基本概念2.4.3 2.4.3 离散GaborGabor变换变换第38页，此课件共88页哦引引 言言连续小波变换（连续小波变换（连续小波变换（连续小波变换（CWTCWT）小波变换的性质小波变换的性质离散小波变换（离散小波变换（DWT）二维小波二维小波多分辨率分析多分辨率分析多分辨率分析多分辨率分析快速小波变换（快速小波变换（快速小波变换（快速小波变换（FWTFWT）2.5 小波变换小波变换第39页，此课件共88页哦1引言引言付利叶等变换的局限付利叶等变换的局限小波的提出、发展和应用小波的提出、发展和应用波和小波波和小波第40页，此课件共88页哦第41页，此课件共88页哦 应应用用：将小波

20、用于地震信号的分析与处理；将二进小波变换用于图像的边缘检测、图像压缩与重构；将连续小波变换用于涡流的研究；将小波变换用于噪声中的未知瞬态信号；将小波变换用于语音信号的分析、变换和综合;将正交小波变换用于算子及拟微分算子的化简；将小波变换的自适应性用于解微分方程；将小波变换用于电磁场领域的若干问题研究等，都取得了初步成果。第42页，此课件共88页哦波和小波（波和小波（波和小波（波和小波（WaveletWavelet）第43页，此课件共88页哦第44页，此课件共88页哦2连续小波变换（连续小波变换（CWT）1)小波变换的定义小波变换的定义 设函数f(t)L2(R)，则小波变换的定义如下：其中，积分

21、核为 的函数族。a0为尺度参数（伸缩参数）,b为定位参数（平移参数），该函数称为小波。若a1函数(t)具有伸展作用，若a1函数(t)具有收缩作用。伸缩参数a对(t)的影响如下图：第45页，此课件共88页哦 随着参数a的减小，(t)的支撑区也随之变窄，反之亦然。(t)随伸缩参数a和平移参数b而变化如下图：大大大大a a小小小小a a第46页，此课件共88页哦图中小波函数为 。当a=2,b=15时，2,15(t)的波形从原点向右移至t=15，且波形展宽。当a=0.5，b=-10时，1/2,-10(t)的波形从原点向左移至t=-10,且波形收缩。第47页，此课件共88页哦2）小波函数要满足的条件(1

22、)紧支撑性紧支撑性（Compact support），即在一个很小 的区域之外函数均为零，函数具有速降特性。(2)平均值为零平均值为零，即：而且其高阶矩也为零：第48页，此课件共88页哦小波应是一个具有振荡性和迅速衰减的波。因为：第49页，此课件共88页哦3）小波反变换 对于所有对于所有f(t),(t)f(t),(t)L L2 2(R)(R)，连续小波逆变换定义为：，连续小波逆变换定义为：变换能量守恒：变换能量守恒：a,ba,b第50页，此课件共88页哦4）几种小波（1）Haar小波第51页，此课件共88页哦（2）Mexico Hat 小波 Mexico Hat Mexico Hat 小波是小

23、波是GaussGauss函数的二阶导数：它是实值小波，一般形式为：它是实值小波，一般形式为：第52页，此课件共88页哦（3）Morlet 小波 Morlet Morlet 小波是最常用的复值小波，它可由下式给出：小波是最常用的复值小波，它可由下式给出：第53页，此课件共88页哦3.小波变换的性质小波变换的性质（1）线性（2）平移和伸缩的共变性（3）小波变换还有微分运算、局部正则、能量守恒、空间-尺度局部化等特性。第54页，此课件共88页哦4离散小波变换（离散小波变换（DWT）离散小波函数、离散小波变换、反变换分别定义如下：离散小波函数、离散小波变换、反变换分别定义如下：第55页，此课件共88页

24、哦5.二维小波二维小波连续的二维小波函数、小波变换和反变换分别如下：第56页，此课件共88页哦二维离散小波变换和反变换分别为：第57页，此课件共88页哦6多分辨率分析多分辨率分析金字塔算法金字塔算法拉普拉斯金字塔编码拉普拉斯金字塔编码子带编码和解码子带编码和解码第58页，此课件共88页哦v拉普拉斯金字塔编码拉普拉斯金字塔编码第59页，此课件共88页哦子带编码和解码子带编码和解码第60页，此课件共88页哦双通道子带编码第61页，此课件共88页哦双通道子带解码第62页，此课件共88页哦7快速小波变换（快速小波变换（FWT）快速小波变换（快速小波变换（FWT）：）：鱼骨算法鱼骨算法第63页，此课件共

25、88页哦快速小波反变换：快速小波反变换：第64页，此课件共88页哦第65页，此课件共88页哦第66页，此课件共88页哦1PCA(主分量分析主分量分析/K-L)变换变换2DCT与与PCA的关系的关系2.6 PCA变换变换第67页，此课件共88页哦1PCA(主分量分析主分量分析/K-L)变换变换均值：均值：均值：均值：偏差：偏差：偏差：偏差：协方差矩阵：协方差矩阵：协方差矩阵：协方差矩阵：第68页，此课件共88页哦PCAPCA变换：变换：变换：变换：PCAPCA反变换：反变换：反变换：反变换：变换后均值为变换后均值为变换后均值为变换后均值为0 0，方差为：，方差为：，方差为：，方差为：第69页，此

26、课件共88页哦作用：解除相关性；可以用于降维处理 也称为主分量分析(K-L)，用于人 脸识别。如果降到M维的均方误差为：第70页，此课件共88页哦2DCT与与PCA的关系的关系其特征值为：其特征值为：其特征值为：其特征值为：第71页，此课件共88页哦其特征向量为：其特征向量为：其特征向量为：其特征向量为：其根为：其根为：其根为：其根为：第72页，此课件共88页哦2.7 离散余弦变换离散余弦变换(DCT)第73页，此课件共88页哦第74页，此课件共88页哦1Walsh-Hadamard 变换变换2Haar 变换变换3斜变换斜变换4.DST 变换变换5.Hartley 变换变换2.8 其它正交变换

27、其它正交变换第75页，此课件共88页哦1Walsh-Hadamard 变换变换第76页，此课件共88页哦第77页，此课件共88页哦2Haar 变换变换第78页，此课件共88页哦一般的，Haar函数定义为：第79页，此课件共88页哦44的Haar变换矩阵为：第80页，此课件共88页哦88的Haar变换矩阵为：第81页，此课件共88页哦3斜（Slant）变换第82页，此课件共88页哦由22矩阵，通过下面的方式产生NN矩阵：第83页，此课件共88页哦第84页，此课件共88页哦44的slant变换矩阵S为：第85页，此课件共88页哦4DST 变换变换第86页，此课件共88页哦5Hartley 变换变换第87页，此课件共88页哦第88页，此课件共88页哦