还剩7页未读,继续阅读
文本内容:
基于小波变换的人脸识别近年来,小波变换在科技界备受重视,不仅形成了一个新的数学分支,而且被广泛地应用于模式识别、信号处理、语音识别与合成、图像处理、计算机视觉等工程技术领域小波变换具有良好的时频域局部化特性,且其可通过对高频成分采取逐步精细的时域取样步长,从而达到聚焦对象任意细节的目的,这一特性被称为小波变换的“变聚焦”特性,小波变换也因此被人们冠以“数学显微镜”的美誉具体到人脸识别方面,小波变换能够将人脸图像分解成具有不同分辨率、频率特征以及不同方向特性的一系列子带信号,从而更好地实现不同分辨率的人脸图像特征提取小波变换的研究背景
4.1法国数学家傅立叶于1807年提出了著名的傅立叶变换,第一次引入“频率”的概念傅立叶变换用信号的频谱特性来研究和表示信号的时频特性,通过将复杂的时间信号转换到频率域中,使很多在时域中模糊不清的问题,在频域中一目了然在早期的信号处理领域,傅立叶变换具有重要的影响和地位定义信号/t为在-8,+8内绝对可积的一个连续函数,则/⑴的傅立叶变换定义如下Fy=4-1\+CQftYi0,dtJ—00傅立叶变换的逆变换为/0=口加204-2171J-8从上面两个式子可以看出,式4T通过无限的时间量来实现对单个频率的频谱计算,该式表明尸
①这一频域过程的任一频率的值都是由整个时间域上的量所决定的可见,式4-1和4-2只是同一能量信号的两种不同表现形式尽管傅立叶变换可以关联信号的时频特征,从而分别从时域和频域对信号进行分析,但却无法将两者有效地结合起来,因此傅立叶变换在信号的局部化分析方面存在严重不足但在许多实际应用中,如地震信号分析、核医学图像信号分析等,研究者们往往需要了解某个局部时段上出现了哪个频率,或是某个频率出现在哪个时段上,即信号的时频局部化特征,傅立叶变换对于此类分析无能为力因此需要一种如下的数学工具可以将信号的时域和频域结合起来构成信号的时频谱,描述和分析其时频联合特征,这就是所谓的时频局部化分析方法,即时频分析法1964年,Gabor等人在傅立叶变换的基础上引入了一个时间局部化“窗函数”gt,改进了傅立叶变换的不足,形成窗口化傅立叶变换,又称“Gabor变换”定义“窗函数”g⑴在有限的区间外恒等于零或很快地趋于零,用函数gt-T乘以/⑴,其⑺=4-3Gf-T e-j0tdt效果等同于在附近打开一个窗口,即:式4-3即为函数ft关于gt的Gabor变换由定义可知,信号/⑴的Gabor变换可以反映j4-4N=dco^ej0tg{t-T^p,TY1Tf该信号在/=T附近的频谱特性其逆变换公式为:可见G/0/的确包含了信号〃t的全部信息,且Gabor窗口位置可以随着T的变化而平移,符合信号时频局部化分析的要求虽然Gabor变换一定程度上克服了傅立叶变换缺乏时频局部分析能力的不足,且兼顾信号的时频分辨率,但其本身仍存在不可克服的局限性,即Gabor窗口不具有自适应性,其大小是固定不变的,因此Gabor变换只能进行单一分辨率的分析而在实际研究中,我们常常希望窗口的大小会随着频率的高低而改变,比如在研究高频信号时,希望窗口开得小一点;反之,在研究低频信号时,则希望窗口开得大一点,这样才更符合实际研究中低频信号分辨率比高频信号分辨率低的特点,因此需要研究更好的解决办法来改善Gabor变换的不足为了克服前面所描述的傅立叶变换和Gabor变换存在的不足,学者们提出了小波变换WaveletTransform小波变换与逆变换
4.
24.
2.1连续小波变换和离散小波变换小波变换是一种窗E1面积即窗口大小固定但窗口形状可变的时频局部化分析方法,其高频部分的时间分辨率较高而频率分辨率较低,而低频部分的频率分辨率较高而时间分辨率较低,因此对信号具有良好的自适应性,被冠以“数学显微镜”的美誉小波变换又可分为连续小波变换和离散小波变换两种连续小波变换的概念是由Morlet等人提出⑺设%为小波变换的核函数,若核函数〃⑴Z C R若满足容许性条件cot4-
5、t-b%/必力=4—6dxa L2则称该函数以幻为基小波一维信号f tCR的连续小波变换可定义为:根据上述对基小波定义,可知j\|°“dxoo,且°xe a0a.b R,在无穷远处趋近于零一般的,记:4-7函数心力x是由基小波函数x经过尺度a的伸缩和b的平移之后所得常用的连续小波包括Morlet小波、Daubechies小波、三次样条小波、Meyer小波和Simoncelli小波等与连续小波变换相对应的是离散小波变换,其一般形式为°W0/X〃2,〃=冲…=|/「[/X/4-8f dx,\-其中以〃⑴=〃也为小波基,劭、%为两个常量且劭〉0离散小波最具I I代表性的为二进小波,即4为2的嘉次2,,%取整由于小波变换在时域和频域同时兼有局部化能力,且能逐步聚焦到对象的任何细节进行分析,因此在人脸识别方面得到众多研究者的关注
4.
2.2几种常见的小波同傅立叶分析不同,小波分析的基小波函数不是唯一存在的,所有满足小波条件的函数都可以作为小波函数,那么小波函数的选取就成了十分重要的问题⑻1Haar小波A.Haar于1990年提出一种正交函数系,定义如下:1/210%1/2X其它这是一种最简单的正交小波,即4-9i//ty/x-ndx=0=±1,±2,n4-102Daubechies dbN小波系该小波是Daubechies从两尺度方程系数%/出发设计出来的离散正交小波一般简写为dbN,N是小波的阶数小波”和尺度函数阈中的支撑区为2NT3的消失矩为N除N=1外Haar小波,dbN不具对称性〔即非线性相位;dbN没有显式表达式除N=1外但仍}的传递函数的模的平N—1方有显式表达式假设公,其中,J为二项式的系数,则有|町⑼/=cos2^Asin2£4-11其中/=万12N-I7Zk=o3Biorthogonal biorNr.Nd小波系Biorthogonal函数系的主要特征体现在具有线性相位性,它主要应用在信号与图像的重构中通常的用法是采用一个函数进行分解,用另外一个小波函数进行重构Biorthogonal函数系通常表示为biorNr.Nd的形式Nr=l Nd=l,3,5Nr=2Nd=2,4,6,8Nr=3Nd=l,3,5,7,9Nr=4Nd=4Nr=5Nd=5Nr=6Nd=8其中,r表示重构,d表示分解4Coiflet coifN小波系coif let也是函数由Daubechies构造的一个小波函数,它具有coifN N=L2,3,4,5这一系列,coiflet具有比dbN更好的对称性从支撑长度的角度看,coifN具有和db3N及syni3N相同的支撑长度;从消失矩的数目来看,coifN具有和db2N及sym2N相同的消失矩数目5SymletsA symN小波系Symlets函数系是由Daubechies提出的近似对称的小波函数,它是对db函数的一种改进Symlets函数系通常表示为symNN=2,3,…,8的形式6Morlet morl小波Morlet函数定义为乎%=0%055%,它的尺度函数不存在,且不具有正交性7Mexican Hatmexh小波2中x=7r-l/4l-x2e-x2/2vr4-12Mexican Hat函数为它是Gauss函数的二阶导数,因为它像墨西哥帽的截面,所以有时称这个函数为墨西哥帽函数/•coi//xdx=Q4-13J-CO墨西哥帽函数在时间域与频率域都有很好的局部化,并且满足由于它的尺度函数不存在,所以不具有正交性8Meyer函数21底/s2u巨@-J27-12hi2兀2r中0=12%一cos40——2〃/205122%4T万,,8万——CD—0兀29T T4-14Meyer小波函数中和尺度函数都是在频率域中进行定义的,是具有紧支撑的正交小波其中,为构造Meyer小波的辅助函数,且有/2〃co——33尸22_4——60——,⑼=卜2»j cos40上-同—13322兀4不24-15069——
4.
2.3二维小波变换与逆变换把对一维的表示推广到二维,考虑二维尺度函数是可分离的情况,可有3个二维小波⑼,则二维尺度函数和小波函数可表示为”x,y=d*y4-16,乂.y=〃4/y4-17设上优,%2}表示一幅离散图像,用低通滤波器人和高通滤波器g分别对的的每一行作滤波,并作隔点抽样,然后再用它们分别对与的每一列滤波并作隔点抽样,得到图像低频概貌J”和图像/高频细节%3%+】,则有如下小波正变换分解算法4-18(〃i,%)=Z E廉2々一%Jg(2%—(匕d+T k)c2k)k\k其小波逆变换(重构算法)如下式c为,%—2n J*2Z百匕—2々%左jk1,k=E2n\n\〃+Z Z2%1—2%g%2—2%“+R,%%i=Z Z以「Xfe-2%%外,%k2n n\n\x+E E虱「2冬解〉-k2n2M+1n,%n\n\x4-19对于NXN像素的图像,小波变换能分解J层,整数,Wlog£o在每一尺度下的.包含前一阶段的低频信息,而
4、明和泪分别包含前一阶段横向、纵向和对角方向的边缘细信息图二维小波正变换框图
4.1图二维小波变换逆变换框图
4.2L表示抽样,为2点取1点的抽样,即只剩下一半样数的分解过程,打表示插样,即得到的样数为原先样数的两倍人脸图像的小波变换
4.3小波在图像处理上的应用思路主要采用将空间或时间域上的图像信号变换到小波域上,成为多层次的小波系数,根据小波的特性,分析小波系数的特点,针对不同需求,结合常规的图像处理方法提出更符合小波分析的新算法来处理小波系数,再对处理后的小波系数进行反变换,将得到所需的目标图像小波变换的优点在于具有良好的时间和频率特性,应用范围较广采用小波分解图像,可降低分解后图像子带的分辨率,大大减少相应的计算复杂度,并可提供更多的空间和频率局部信息图像输入图像输出小波正变换图像处理小波变小波反变换呈图由于一幅图像的信息主要包含在低频部分,而图像细节体现在高频部分,故可通过小波变换得到低频系数,也就是图像的主要信息在本算法中,采用二维离散小波变换,进行二维小波频域分解,实际是进行两次一维小波变换经过一级小波变换之后,一幅图像被分解成为如图
4.4(a)所示的4个子带包括低频区域LL,高频区域LH、HL、HH,每个区域都是一幅图像LL区域表示的是原图像的平滑图像,它包含了原图像的大部分信息,刻画人脸表情和姿态的不变特征;LH区域保持了原图像的垂直边缘细节,体现人脸的轮廓和鼻子的垂直特征;HL区域保持了原图像的水平边缘细节,体现了人脸的眼睛、嘴巴的水平特征;HH区域保持了原图像的对角线细节,受噪声、表情和姿势的影响较大经过一级小波分解后,每个区域的图像维数都是原图像维数的四分之一图中子带LL是低频成分,集中了原始图像的大部分信息,LH和HL反映了沿水平和垂直方向的图像变化,HH子带则代表图像的高频成分还可以运用小波算法继续对LL子带进行分解,图
4.4b为图
4.4aLL LHHLHH的二级小波分解图组实现.其中,低通滤波器起平滑作用,得到图像的缓变成分;带通滤波器起差分作用,得到图像的高频成a-波分解b二级分.使用具有紧支集的规范正交小波基可以构造这类滤波器图
4.
5、
4.6分别是人脸图像的一小波分解级分解、二级分解图像图二维图像的小波分
4.4图一级小波分解图像图二级小波分解图像
4.
54.6通常使用等价的滤波器人本脸文图研像究的的二方维法小是波分首先对原始图像进行进行小波变换与反变换,小波变换后,图像左上角表现了图像的低频分量,它是构成图像信息及其能量的主要的组成部分右上角处表现了图像的水平方向的高频分量左下角表现垂直高频分量而右下角处表现了图像的对角信息而且将图像第一次分解后的高频分量除去后,利用其余的系数仍然可以得到完整的重构图像通过适当的保留或者去除小波变换后的一些特定的分量,达到滤波的目的高通或者低通图图像的小波变换与反变换处理
4.7对二层离散小波重构的图像用PCA方法获取图像的特征矩阵⑵,把训练图像和测试图像投影到该方法所得到的子空间上,利用Fisher线性判别法完成人脸的识别其人脸识别流程图如图
4.8所示。