专利名称：一种基于asm算法的人脸特征定位方法
技术领域：
本发明属于图像处理技术领域，主要涉及生物特征鉴别中的人脸识别技术。
背景技术：
在当今信息化吋代，如何准确鉴定一个人的身份，保护信息安全是一个必须解决的关键 社会问题。为此，生物特征鉴别技术悄然新起，并成为目前世界信息安全管理领域的前沿研 究课题。生物特征鉴别技术是指利用人体所固有的生理特征或行为特征来进行个人身份鉴定。 人脸识别技术是生物特征鉴别技术的一个分支，是计算机图像处理技术和模式识别技术在个 人身份识别领域的应用，由于其终身不变性和图像来源不需要人特别配合等特性，近年來己 经成为生物特征鉴别的热门发展方向。人脸自动识别技术，在居民小区、公司单位的门禁系 统、银行、公安、机场、网络等方面应用广泛，具有巨大的经济价值和现实意义。现在它己 经使用在边境检査、通观认证、银行提款、信息管理和建筑物安全管理等，还可以使人们摆 脱记忆信用卡号、银行帐号、身份证号、网络登录号的繁琐。随着数字信号处理技术和图像 处理技术的发展，人脸识别已经越来越受到人们的重视。详见文献John G. Daugman, "High Confidence Recognition of Persons," 77^ PraceW"g 7EE五3J^ /她rw加》wa/ Qir""/zaw Co"，eweSecwr/(v 7fec/z"o/ogy, pp. 254-263, 2000所述。
在人脸识别技术中，人脸图像的检测定位是人脸识别的第一步，它的执行准确度将直接 影响整个人脸识别系统的识别率。在实际中，由于人脸区域常常受到头发，眼睛，胡子和一 些饰物的遮挡和光照的影响，人脸定位算法的准确性和有效性还有待进一步提高。如何在存 在头发，眼睛，胡子和一些饰物的遮挡和光照影响问题的低质量人脸图像中，快速精确地定 位出人脸，并对其边界或位置用数学模型进行描述是我们研究的主要问题。详见文献JohnG. Daugman, "High Confidence Visual Recognition of Persons by a Test of Statistical Independence," 7hmsac"'OM尸"股r"朋d Mac/zz'股/"te〃/gewce, volume 15, no. 11, pp. 1148-1161,
1993。
现在通常使用的人脸定位的方法大致有以下三大类，应用最为广泛的是基于统计学的定 位方法主动形状模型(即active shape model, ASM)算法、主动外观模型(即active appearance model, AAM)算法等。
(1)基于肤色的彩色图定位方法。它通过肤色模型进行粗定位，寻找含有肤色的区域 作为人脸的候选区，然后再在这些候选区利用人脸模板进行匹配或者利用圆检测定位眼睛从 而定位出人脸区。在实际应用中此类方法对光照有一定的要求，光照改变的情况下定位结果会出现不可预料的结果。详见文献Rein-Lien Hau and Anil K.Jain. Face Detection in Color Image [J]. IEEE Trans. Pattern Analysis and Machine Intelligence, may 2002.禾口文献M.H. Yang, N.A. huju. Detecting human faces in color images [J]. /VoceefZ/wgA. Co"y^尸e"ce /附(3g/wg iVocesi7."g, 1998， 1:127-130.
(2)基于变换域的定位方法。基于小波变换的人脸定位方法。它是通过一定的算子， 提取出人脸图像中的不同频率的成分，从而搜索到人脸所在的位置。其缺点是在于背景细节 的干扰，使得人脸定位结果不准确。详见文献Jarmo Ilonen， Joni-Kristian Kamarainen. "Image Feature Localization by Multiple Hypothesis Testing of Gabor Features" 7 awracrio附 /扁ge praceM/wg, w/.77, no.3, March 2008。
基于统计学的定位方法。不管人脸的哪个构成的部分都有轮廓，这种方法通过对感兴趣 的轮廓事先进行训练，后再进行搜索，做到了对定位的目标有针对性，从而能达到比较准确 的效果。详见文献Kwok-Wai Wan, Kin-Man Lam , Kit-Chong Ng. "An accurate active shape model for facial feature extraction" 26 (2005) 2409—2423 。
发明内容
本发明的任务是提供一种基于ASM算法的人脸检测定位方法，它具有在不同光照情况 下定位准确的特点。
为了方便地描述本发明内容，首先对一些术语进行定义。
定义l:轮廓特征点。是指能够表征人脸形状和五官形状的外围轮廓上的离散点。
定义2:配准。由于训练样本中的不同图像的特征点是在各自的坐标下进行标定的，其
绝对位置和图像尺寸存在差异，配准就是把两幅图像的通过某种变换到同一坐标系中，使差
异达到某种程度的统一，各S的形状特点保持不变。
定义3:奇异值分解。wxw阶矩阵^可以写成^=^51/'的形式，t/为m阶正交阵，F为 n阶正交阵，S-c^agfa,,^"."^)^, >0(!' = l_,r),r = M"/t(J) , [/和F中分别是^4的奇异
向量，而S是^的奇异值。奇异值分解和特征值问题是紧密联系的，X4'的正交单位特征向 量组成，特征值组成， Ji的正交单位特征向量组成F ，特征值(与Jl相同)组成， 奇异值分解提供了一些关于^的信息，例如非零奇异值的数目(S的阶数)和J的秩相同， 一旦秩r确定，那么t/的前r列构成了」的列向量空间的正交基。 本发明的技术方案如下
一种基于ASM算法的人脸特征定位方法，如图1所示，包含下列歩骤 歩骤1.选取人脸库中的M张图片作为样本集并对其进行特征点的标定；对人脸库中的每张样本图片的主要器官的轮廓进行w个特征点的手工标定，则一张样本 图片由w个特征点的坐标值1 = (、...^ ,乂,...，;;,,)7'表示；M张图片的样本集表示成
X4 = ",X2,...,ZM);
步骤2.采用Procrustes Analysis算法对样本集x4进行配准，如图2所示，具体包括； 步骤2-1.选取样本集x4中的第A:列向量x,为基准，计算第A:列向量和第m列向量
的平均值&和c",，其中1SA:《M， w^A:且l《m^M;
步骤2-2.计算矩阵5 = ^;,-/^;)^^-其中/为一个和X具有相同维数的单位
向量；并对矩阵5进行奇异值分解，得到^=^/5"，其中t/为2/7阶正交矩阵，F为2"阶正 交矩阵，S是矩阵5的奇异值。
步骤2-3.构造旋转矩阵i = f/f^ ，位移矩阵T = / ^ - icOT 。
步骤2-4. x^经过以j^为基准进行配准后，得到的向量i;二fj^+r。
步骤2-5.对x4中的每一列不同于A的向量做步骤2-2至步骤2-4的变换，得到与x4同 样大小的新的样本集z4 = 其中《=4。
步骤3.计算平均人脸模型P， F-丄f^;
步骤4.采用对三种分辨率下的图片进行轮廓搜索的策略，建立统计灰度搜索模型(图3
所示)，具体方法如下
针对人脸库中的每张样本图片的每一个特征点，以该特征点的坐标为中心，取其法线方 向两侧各w个像素点的灰度值按顺序组成向量g"&,&,…,g^wf ;然后求出向量
g = (gpg2，...，g(2—f的一阶差分向量《=(&-其次构建该特征点的灰度
搜索模型g、g'/l;lg(叫-g,;接着构建该人脸样本图片的搜索模型G'-(g'"g ，...,/ );
最后构建人脸库中整个人脸样本图片集的统计灰度搜索模型浐。
步骤5.对待定位人脸图像i;进行搜索定位，包括以下具体步骤； 步骤5-i.搜索的初始位置为y = m+r ，其中f= +户6为待定位人脸图像j;的 人脸模型，其中i^(p,，/72,…，A) , 6二(4 4..4) ， 6由协方差矩阵￡ = lf;(r, -7"y,的f个最大的特征值从大到小的顺序排列而成，p由6中特征 Af -1*、 7 八" 7
值对应的正交特征向量构成，/的取值为不等式J]A,々a24确定的最小值，权值系数a通
,-=1 i=
常取0.98或0.95， fc向量初始元素值为0;
步骤5-2.根据待定位人脸图像1;的人脸模型5 = ]7+/^在待定位人脸图像1;中提取n个
特征点，为每个特征点坐标为屮心，取其法线方向两侧各/个像素点的灰度值，根据步骤4 的方法计算得到本次搜索的局部灰度模型7，其中/>沐；
步骤5-3.匹配统计灰度搜索模型f和局部灰度模型；'以确定n个特征点的移动方向c/r (如图4所示)，移动方向c^确定后，姿态参数控制的等式(y + dy卜M(l + A,t^)[y] + ^r成
步骤5-4.位置参数t^由等式M(s(i+ife),6)+^)[y+办]+:r+^r"+w确定，等式变 步骤5-5. ;r+办aF+p(6+w)成立，更新形状参数必二p—、办。
步骤5-6.循环5-1到5-5，直至本次循环与上次循环结果的"个特征点的坐标值之和的 差异值占本次循环结果的比例小于阈值Tg时(阈值Tg的取值范围可以是(0，0.001])，搜索
完成，定位结束，得到待定位人脸图像i;的人脸模型f = 7+外。
需要说明的是
1. 步骤2中对训练样本X4采用Procrustes Analysis方法进行配准时选用的基准向量 XJ1S &《M)为样本集X4中的任意一列。
2. 步骤2-5中Z4 = ("，:r2，..":rM)是原始的训练样本集X4经过Procrustes Analysis算法配
准后在同一坐标系下的坐标数据。
3. 歩骤5中模型的搜索模型必须在待识别图像的搜索模型中进行移动搜索，才能找到最 佳匹配，因此模型的灰度搜索模型维数/必须大于w。
本发明首先对一组人脸样本图片进行手工标定提取特征点，然后采用Procrustes Analysis 算法对样本集进行配准，得到T均人脸模型；然后采用对三种分辨率下的图片进行轮廓搜索
的策略，建立统计灰度搜索模型；最后对统计灰度搜索模型和待定位人脸图像j;的局部灰度模型釆取匹配运算，对待定位人脸图像i;进行循环迭代式搜索定位。本发明采用ASM算法
和Procrustes分析方法相结合的方法，通过Procrustes分析方法对坐标系、点模型分布不均和 形状大小存在差异的图像进行配准，从而获得比较准确的姿态参数，能够达到使数据进行更 精确变换的目的。利用Procrustes分析方法和ASM算法结合定位人脸是本发明的一个特色， 与一般的人脸定位方法相比，本发明定位准确率要高的同时速度也在一定程度上有所提高， 在人脸定位方面具有很强的通用性。
图1是整体定位方法流程图。
图2是Procrustes Analysis算法流程图。
图3是灰度搜索模型示意图。
图4是搜索过程的示意图。
具体实施例方式
采用本发明的方法，使用C语言进行程序的编写，然后在matlab的平台上进行仿真实验 并得到结果。对imm—face—db人脸库中的240张拍摄好的、包括40个人在不同光照条件、不 同表情和不同姿势的灰度人脸图像作为源数据，与传统的ASM算法结果的数据分析相比较， 定位的准确率有较大的提高。
综上所述，本发明的方法充分利用人脸的特征信息，结合ASM算法的优势，从而实现快 速准确地从所提供的人脸原始图像中检测定位人脸区域。
权利要求
1、一种基于ASM算法的人脸特征定位方法，包含下列步骤步骤1.选取人脸库中的M张图片作为样本集并对其进行特征点的标定；对人脸库中的每张样本图片的主要器官的轮廓进行n个特征点的手工标定，则一张样本图片由n个特征点的坐标值X＝(x1，...，xn，y1，...，yn)T表示；M张图片的样本集表示成XA＝(X1，X2，...，XM)；步骤2.采用Procrustes Analysis算法对样本集XA进行配准，具体包括；步骤2-1.选取样本集XA中的第k列向量Xk为基准，计算第k列向量Xk和第m列向量Xm的平均值ck和cm，其中1≤k≤M，m≠k且1≤m≤M；步骤2-2.计算矩阵B＝(Xm-I·cm)(Xk-I·ck)T，其中I为一个和X具有相同维数的单位向量；并对矩阵B进行奇异值分解，得到B＝USVT，其中U为2n阶正交矩阵，V为2n阶正交矩阵，S是矩阵B的奇异值；步骤2-3.构造旋转矩阵R＝UVT，位移矩阵T＝I·ck-R·cm；步骤2-4.Xm经过以Xk为基准进行配准后，得到的向量Ym＝R·Xm+T；步骤2-5.对XA中的每一列不同于Xk的向量做步骤2-2至步骤2-4的变换，得到与XA同样大小的新的样本集YA＝(Y1，Y2，...，YM)，其中Yk＝Xk；步骤3.计算平均人脸模型<overscore>Y</overscore>，<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><mover> <mi>Y</mi> <mo>‾</mo></mover><mo>=</mo><mfrac> <mn>1</mn> <mi>M</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>M</mi></munderover><msub> <mi>Y</mi> <mi>i</mi></msub><mo>;</mo> </mrow>]]></math> id="icf0001" file="A2009100596480002C1.tif" wi="22" he="10" top= "174" left = "85" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>步骤4.采用对三种分辨率下的图片进行轮廓搜索的策略，建立统计灰度搜索模型(图3所示)，具体方法如下针对人脸库中的每张样本图片的每一个特征点，以该特征点的坐标为中心，取其法线方向两侧各w个像素点的灰度值按顺序组成向量g＝(g1，g2，...，g(2w+1))T；然后求出向量g＝(g1，g2，...，g(2w+1))T的一阶差分向量g′＝(g2-g1，g3-g2，...，g(2w+1)-g2w)；其次构建该特征点的灰度搜索模型<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><msup> <mi>g</mi> <mo>*</mo></msup><mo>=</mo><msup> <mi>g</mi> <mo>′</mo></msup><mo>/</mo><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mrow><mn>2</mn><mi>w</mi> </mrow></munderover><mo>|</mo><msub> <mi>g</mi> <mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo> </mrow></msub><mo>-</mo><msub> <mi>g</mi> <mi>i</mi></msub><mo>|</mo><mo>;</mo> </mrow>]]></math> id="icf0002" file="A2009100596480002C2.tif" wi="39" he="10" top= "237" left = "37" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>接着构建该人脸样本图片的搜索模型G*＝(g*1，g*2，...，g*n)；最后构建人脸库中整个人脸样本图片集的统计灰度搜索模型<maths id="math0003" num="0003" ><math><![CDATA[ <mrow><msup> <mover><mi>Y</mi><mo>‾</mo> </mover> <mo>*</mo></msup><mo>=</mo><mfrac> <mn>1</mn> <mi>M</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>M</mi></munderover><msub> <msup><mi>G</mi><mo>*</mo> </msup> <mi>i</mi></msub><mo>;</mo> </mrow>]]></math> id="icf0003" file="A2009100596480002C3.tif" wi="26" he="10" top= "254" left = "130" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>步骤5.对待定位人脸图像Ys进行搜索定位，包括以下具体步骤；步骤5-1.搜索的初始位置为y＝M(s，θ)[Ys]+T，其中Ys＝Y+Pb为待定位人脸图像Ys的人脸模型，其中P＝(p1，p2，...，pt)，b＝(λ1 λ2...λt)，b由协方差矩阵<maths id="math0004" num="0004" ><math><![CDATA[ <mrow><mi>E</mi><mo>=</mo><mfrac> <mn>1</mn> <mrow><mi>M</mi><mo>-</mo><mn>1</mn> </mrow></mfrac><munderover> <mi>Σ</mi> <mrow><mi>j</mi><mo>=</mo><mn>1</mn> </mrow> <mi>M</mi></munderover><mrow> <mo>(</mo> <msub><mi>Y</mi><mi>j</mi> </msub> <mo>-</mo> <mover><mi>Y</mi><mo>‾</mo> </mover> <mo>)</mo></mrow><msup> <mrow><mo>(</mo><msub> <mi>Y</mi> <mi>j</mi></msub><mo>-</mo><mover> <mi>Y</mi> <mo>‾</mo></mover><mo>)</mo> </mrow> <mi>T</mi></msup> </mrow>]]></math> id="icf0004" file="A2009100596480003C1.tif" wi="48" he="10" top= "48" left = "24" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>的t个最大的特征值从大到小的顺序排列而成，P由b中特征值对应的正交特征向量构成，t的取值为不等式<maths id="math0005" num="0005" ><math><![CDATA[ <mrow><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>t</mi></munderover><msub> <mi>λ</mi> <mi>i</mi></msub><mo>≥</mo><mi>α</mi><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mrow><mn>2</mn><mi>n</mi> </mrow></munderover><msub> <mi>λ</mi> <mi>i</mi></msub> </mrow>]]></math> id="icf0005" file="A2009100596480003C2.tif" wi="22" he="9" top= "64" left = "109" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>确定的最小值，权值系数α通常取0.98或0.95，b向量初始元素值为0；步骤5-2.根据待定位人脸图像Ys的人脸模型Ys＝Y+Pb在待定位人脸图像Ys中提取n个特征点，为每个特征点坐标为中心，取其法线方向两侧各l个像素点的灰度值，根据步骤4的方法计算得到本次搜索的局部灰度模型y*，其中l＞w；步骤5-3.匹配统计灰度搜索模型Y*和局部灰度模型y*以确定n个特征点的移动方向dY，移动方向dY确定后，姿态参数控制的等式(y+dY)＝M(1+ds，dθ)[y]+dT成立；步骤5-4.位置参数dy由等式M(s(1+ds)，θ+dθ)[Y+dy]+T+dT＝y+dY确定，等式变换得到dy＝M((s(1+ds))-1，-(θ+dθ))[M(s，θ)[Y]+dY-dT]-y；步骤5-5.Y+dy≈Y+P(b+db)成立，更新形状参数db＝P-1·dy；步骤5-6.循环5-1到5-5，直至本次循环与上次循环结果的n个特征点的坐标值之和的差异值占本次循环结果的比例小于阈值Tg时，搜索完成，定位结束，得到待定位人脸图像Ys的人脸模型Ys＝Y+Pb。
2、根据权利要求1所述的基于ASM算法的人脸特征定位方法，其特征在于，步骤5-6 中所述阈值Tg的取值范围是(0,0.001]。
全文摘要
一种基于ASM算法的人脸特征定位方法，属于图像处理技术领域，主要涉及生物特征鉴别中的人脸识别技术。本发明首先对一组人脸样本图片进行手工标定提取特征点，然后采用Procrustes Analysis算法对样本集进行配准，得到平均人脸模型；然后采用对三种分辨率下的图片进行轮廓搜索的策略，建立统计灰度搜索模型；最后对统计灰度搜索模型和待定位人脸图像Y<sub>s</sub>的局部灰度模型采取匹配运算，对待定位人脸图像Y<sub>s</sub>进行循环迭代式搜索定位。本发明采用ASM算法和Procrustes分析方法相结合的方法，可以有效地提高人脸定位速度和精度，具有很强的通用性。
文档编号G06K9/00GK101593272SQ200910059648
公开日2009年12月2日 申请日期2009年6月18日 优先权日2009年6月18日
发明者华 徐, 梅 解 申请人:电子科技大学