CN102915436A - 基于类内变化字典和训练图像的稀疏表示人脸识别方法 - Google Patents

基于类内变化字典和训练图像的稀疏表示人脸识别方法 Download PDF

Info

Publication number
CN102915436A
CN102915436A CN2012104117710A CN201210411771A CN102915436A CN 102915436 A CN102915436 A CN 102915436A CN 2012104117710 A CN2012104117710 A CN 2012104117710A CN 201210411771 A CN201210411771 A CN 201210411771A CN 102915436 A CN102915436 A CN 102915436A
Authority
CN
China
Prior art keywords
centerdot
image
face
class
dictionary
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012104117710A
Other languages
English (en)
Other versions
CN102915436B (zh
Inventor
邓伟洪
胡佳妮
郭军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN201210411771.0A priority Critical patent/CN102915436B/zh
Publication of CN102915436A publication Critical patent/CN102915436A/zh
Application granted granted Critical
Publication of CN102915436B publication Critical patent/CN102915436B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Collating Specific Patterns (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于类内变化字典和训练图像的稀疏表示人脸识别方法,用于解决现有方法在小样本、光照不均、人脸遮挡和表情变化下的局限性问题,提高了人脸识别精度。该方法的实现步骤为:(1)对训练图像集合和测试人脸图像进行图像特征提取,形成训练图像矩阵和测试图像向量,并分别进行归一化处理;(2)从人脸数据库中统计相同人脸在不同外部环境影响下呈现出来的图像纹理差异,形成人脸的类内变化字典;(3)将测试图像表示成为训练图像矩阵和类内变化字典的线性组合,通过L1范数最小化准则求得最优稀疏表示系数;(4)求得原测试图像与由各类别训练图像和类内变化字典重构得到的重构图像的残差,代入类别判定公式得到识别结果。

Description

基于类内变化字典和训练图像的稀疏表示人脸识别方法
技术领域
本发明涉及一种人脸识别方法,确切地说,涉及一种基于类内变化字典和训练图像的稀疏表示人脸识别方法,属于计算机的图像处理和模式识别技术领域。
背景技术
人脸识别是利用计算机对人脸图像进行分析,从中提取有效的鉴别特征,实现身份认证的一种技术。与指纹、虹膜、掌纹等生物特征相比,利用人脸进行身份认证的优势是使用方便、成本低廉、操作隐蔽等。其应用场景十分广泛,例如公安监控、司法认证、民航安检、口岸出入控制、智能门禁等。
人脸识别系统通常包括人脸检测、特征提取和识别算法三个部分。传统的人脸识别技术的研究侧重于特征提取和识别算法,并已经形成了一些经典方法,例如,主成分分析、线性鉴别分析、弹性匹配、神经网络等。有一些识别算法在实验数据中取得了较好的效果,然而,在现实应用中,人脸识别方法仍然面临着下述诸多挑战:
(1)小样本问题:训练集合中的每个类别(人脸)仅有一个或几个训练图像,训练图像的单一性和测试图像的多态性就会导致识别精度较低。
(2)光照问题:光照的角度、强度和传感器的不同,都会导致同一人脸在不同环境下采集的人脸图像存在巨大差异。
(3)人脸遮挡和/或表情问题:墨镜、胡须、发型、化装等因素和/或喜怒哀乐悲恐惊等不同表情都会导致人脸图像模式的多样性。
上述多重因素的相互干扰和叠加,使得人脸识别技术成为一项极富挑战性的课题,因此,面向现实环境的人脸识别技术已经成为当前计算机图像处理和模式识别领域的研究热点。
2009年John Wright等人将压缩感知理论引入人脸识别技术的研究中,提出了基于稀疏表示的人脸识别方法,该方法开创性地指出:测试人脸图像一旦能够表示为训练图像的稀疏线性组合,则特征提取和识别算法将不再重要。基于稀疏表示的人脸识别方法为解决噪声、光照、遮挡、表情等人脸识别难题提供了全新的思路。但是,其高精度识别结果是在严格限定的条件下获得的:对于每个候选人脸图像,至少需要32个不同光照下的训练图像,才足以线性表示复杂环境下的光照。显然,在目前的大部分现实应用中,无法提供如此充分的训练数据,所以传统的基于稀疏表示的人脸识别方法在现实环境中的应用至今非常有限。因此,如何寻找更加简便、实用的基于稀疏表示的人脸识别方法,已经成为业内科技人员关注的焦点课题。
发明内容
有鉴于此,本发明的目的是提供一种基于类内变化字典和训练图像的稀疏表示人脸识别方法,利用类内变化字典和训练图像的稀疏表示来解决现有技术在人脸图像处于小样本、光照不均、人脸遮挡和表情变化情况时差异描述上的不足,并在降低对人脸训练图像数量需求的同时,提高光照、遮挡或表情等复杂环境下的人脸识别精度。
为了达到上述目的,本发明提供了一种基于类内变化字典和训练图像的稀疏表示人脸识别方法,其特征在于:所述方法包括下列操作步骤:
(1)对训练图像集合中的所有人脸图像进行特征提取,得到训练图像矩阵A:A=[A1,A2,...,Ai,...,Ak],且
Figure BDA00002305250400021
并对该训练图像矩阵A中的每一列进行归一化处理,使每一列向量的模都为1;式中,自然数下标i为人脸类别序号,其最大值为k;d1和n分别为训练图像的维数和个数,且第i类训练人脸图像矩阵
Figure BDA00002305250400022
其中的ni为第i类训练人脸图像的个数,且R为实数集;
(2)对测试图像进行特征提取,得到测试图像向量
Figure BDA00002305250400024
并对该测试图像向量y进行归一化处理;式中,d2为测试图像的维数;
(3)从训练图像集合或其他开放人脸图像集合中获取用于描述相同类别人脸在不同环境因素下的类内变化,构成类内变化字典,即矩阵式中,d3为类内变化字典的行数,且d3与d2、d1的数值都相等;p为类内变化字典中描述设定的人脸图像变化的列数;
(4)用训练图像矩阵A和类内变化字典D线性表示测试图像向量
Figure BDA00002305250400031
其表达式为:y=Ax+Dβ+z,式中,x和β分别为线性的稀疏系数,z为噪声项,且x∈Rn×1,β∈Rp×1
Figure BDA00002305250400032
(5)因稀疏系数x和β具有多值性,故采用基于L1范数最小化方法求解得到其最优稀疏系数
Figure BDA00002305250400033
Figure BDA00002305250400034
其目标函数为: x ^ β ^ = arg min | | x β | | l , 且满足 | | [ A , D ] x β - y | | 2 ≤ ϵ ; 式中,ε为预设的误差,且最优稀疏系数 x ^ = [ x 1 , x 2 , · · · x t · · · , x n ] ∈ R n × l
Figure BDA00002305250400038
其中,自然数下标t为最优稀疏系数序号,其最大值为n,该
Figure BDA00002305250400039
中的第t个元素xt是重构测试图像向量y时第t个训练图像的系数;
(6)分别保留最优稀疏系数
Figure BDA000023052504000310
中与第i类训练图像相对应的系数,并将
Figure BDA000023052504000311
中与第i类训练图像无关的系数全部置零,得到下述带通系数
Figure BDA000023052504000312
δ 1 ( x ^ ) = x 1 , x 2 , · · · , x n 1 , 0 , · · · , 0 T ,
δ i ( x ^ ) = 0 , · · · , 0 , x n 1 + n 2 + · · · + n i - 1 + 1 , x n 1 + n 2 + · · · + n i - 1 + 2 , · · · , x n 1 + n 2 + · · · + n i - 1 + n i , 0 , · · · , 0 T ,
……
δ k ( x ^ ) = 0 , · · · , 0 , x n - n k + 1 , x n - n k + 2 , · · · , x n T ;
(7)根据步骤(6)获得的带通系数,分别用每个类别的训练图像和类内变化字典对测试图像y进行重构,得到相应的重构图像向量qi:
Figure BDA000023052504000316
式中i=1,2,...,k;
(8)分别计算原测试图像向量y与重构图像向量qi的残差ri(y)、即两者差值的模:ri(y)=‖y-qi2,式中i=1,2,...,k;残差越小表示重构图像与测试图像越近似或越相像;
(9)将该残差ri(y)代入测试图像向量y的类别判定公式:Identity(y)=argminiri(y),求解得到其中所有k个残差中的最小值,并将该最小残差所对应的类别序号作为最终的人脸识别结果Identity(y)。
本发明与现有方法相比较,具有以下优点:
本发明首创提出一种将测试图像表示为训练图像矩阵和类内变化字典的稀疏线性组合的方法。该类内变化字典用于描述同一人脸在光照、遮挡和/或表情变化时的图像纹理差异。通过补偿测试图像与其相同类别的训练图像之间的差异,类内变化字典可以提高人脸识别算法在表情、遮挡和/或光照变化环境下的识别精度。
另外,在每类训练图像较少或者每类只有单个训练图像的情况下,本发明首创提出的类内变化字典可从独立于训练图像集合的其他开放人脸图像集合中获取,为解决小样本和单样本情况下的人脸识别问题提供了一种全新思路。因为传统的稀疏表示方法只适用于样本充足的情况,基于该性能,本发明使得训练图像稀疏表示方法在小样本和单样本情况下同样能够获得较高的识别精度,因此该方法的推广应用前景看好。
附图说明
图1是本发明基于类内变化字典和训练图像的稀疏表示人脸识别方法示意图。
图2是本发明基于类内变化字典和训练图像的稀疏表示人脸识别方法操作步骤流程图。
图3是本发明实施例中采用的训练图像和测试图像样例的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步的详细描述。
众所周知,人脸识别的核心问题是寻找出同一人脸的共性特征和差异特征,一旦能够有效刻画同一人脸的差异,就能实现复杂、真实环境下的高精度人脸识别。本发明旨在克服现有技术在人脸差异刻画上的不足,提供一种基于类内变化字典和训练图像的稀疏表示人脸识别方法。
类内变化是指同一人脸在不同外部环境下所呈现的差异。例如,某一人脸在佩戴墨镜时图像与其正常光照和表情下的图像的差异表现为墨镜。本发明首创的类内变化字典D用于描述每类人脸在包括不同光照和/或遮挡、和/或表情的复杂外部环境下所呈现的各种图像纹理差异的集合,以便能够设法找出同一人脸的差异特征。类内变化字典中的每一列描述人脸的一种设定的差异变化,也就是将某一个人脸在异常光照和/或遮挡、和/或表情的图像与其在正常光照和表情、且无遮挡的图像相减得到的差值(例如,将某一人脸在佩戴墨镜时的图像和其正常光照和表情下的图像做减法所得的结果即为墨镜;将某一人脸在异常光照下的图像和其正常光照和表情下的图像做减法时的结果就是某种特殊光照)。不同人脸的差异往往存在共性,因此一旦能够有效刻画同一人脸图像的差异,就能实现复杂、真实环境下的高精度人脸识别。分别用每个类别的训练图像矩阵和类内变化字典的线性组合得到重构图像,根据原测试图像与重构图像的残差来实现高精度的人脸识别。图1展示了本发明人脸识别方法的工作原理:即基于类内变化字典和训练图像的稀疏表示人脸识别方法。
参见图2,介绍本发明方法的具体操作步骤:
步骤1,对训练图像集合中的所有人脸图像进行特征提取,得到训练图像矩阵A:A=[A1,A2,...,Ai,...,Ak],且
Figure BDA00002305250400051
并对该训练图像矩阵A中的每一列进行归一化处理,使每一列向量的模都为1;式中,自然数下标i为人脸类别序号,其最大值为k;d1和n分别为训练图像的维数和个数,且第i类训练人脸图像集合
Figure BDA00002305250400052
其中的ni为第i类训练人脸图像的个数,且R为实数集。
步骤2,对测试图像进行特征提取,得到测试图像向量并对该测试图像向量y进行归一化处理;式中,d2为测试图像的维数。
步骤3,从训练图像集合或其他开放人脸图像集合中获取用于描述相同类别的人脸在不同环境因素下的类内变化,构成类内变化字典,即矩阵
Figure BDA00002305250400055
式中,d3为类内变化字典的行数,且d3与d2、d1的数值都相等;p为类内变化字典中描述设定的人脸图像变化的列数。
该步骤使用其他开放人脸图像集合构建或获取类内变化字典时,若发生该人脸图像集合中的图像维数与训练图像维数d1不相等的情况,则采用包括插值法的图像缩放方法,使得其他开放人脸图像集合中图像的维数等于训练图像的维数d1,以保证类内变化字典的行数d3与训练图像的维数d1相等。
构建类内变化字典有下述两种方法:
第一种方法是求解得到用不同光照和/或遮挡、和/或表情的某个人脸图像减去其在所有环境中的平均图像所得到的差值向量的集合,其操作步骤如下:
(A1)从人脸数据库中随机选择与训练人脸不同的m个人,再对每个人选取w幅人脸图像,用于计算类内变化字典;还对每幅人脸图像进行特征提取,将每幅人脸图像表示成一个d3维的列向量;
(A2)分别将每个人的w幅人脸图像标记为矩阵
Figure BDA00002305250400061
再求解每个人的w幅图像的平均人脸图像
Figure BDA00002305250400062
式中,自然数下标j为计算类内变化字典时的人脸类别序号,其最大值为m;矩阵Dj的第b列向量自然数下标b为矩阵Dj的列序号,其最大值为w;
(A3)从m个人的所有人脸图像中,计算得到人脸的类内变化字典D: D = [ D 1 - c 1 × e 1 , D 2 - c 2 × e 2 , · · · , D j - c j × e j , · · · , D m - c m × e m ] ∈ R d 3 × ( w × m ) ; 式中,cj为列向量,ej为全1的行向量,且e1=e2=…=ej=…=em=[1,...,1]∈R1×w,cj×ej是每列均为cj且维数为d3×w的矩阵,(Dj-cj×ej)表示矩阵Dj的每一列都减去其对应的平均人脸图像cj
(A4)采用上述方法得到类内变化字典的列数p=w×m;
(A5)对类内变化字典D中的每列元素进行归一化处理,使得每列元素的模值都为1。
第二种方法是求解得到同一人脸在包括不同光照和/或遮挡、和/或表情的各种环境下获取的不同图像之间每两个图像的差值向量集合,其操作步骤为:
(B1)从人脸数据库中随机选择与训练人脸不同的m个人,再对每个人选取w幅人脸图像,用于计算类内变化字典;还对每幅人脸图像进行特征提取,将每幅人脸图像表示成一个d3维的列向量;
(B2)从每个人的w幅人脸图像中,任意选取其中两幅人脸图像,并计算这两幅人脸图像的差值向量,共完成
Figure BDA00002305250400065
种选取和计算后,将所有差值向量组合成差异矩阵Qj: Q j = [ a j 1 - a j 2 , a j 1 - a j 3 , · · · , a js - a j ( s + 1 ) , · · · , a j ( w - 1 ) - a jw ] ∈ R d 3 × [ w ( w - 1 ) / 2 ] ; 式中,自然数下标j为计算类内变化字典时的人脸类别序号,其最大值为m;第j个人脸图像集合中的第s幅图像
Figure BDA00002305250400067
且自然数序号s的最大值为w;
(B3)从m个人的所有人脸图像中,计算得到人脸的类内变化字典D:
D = [ Q 1 , Q 2 , . . . , Q j , . . . , Q m ] ∈ R d 3 × { m × [ w ( w - 1 ) / 2 ] } ;
(B4)采用上述方法得到类内变化字典的列数
(B5)对类内变化字典D中的每列元素进行归一化处理,使得每列元素的模值都为1。
步骤4,用训练图像矩阵A和类内变化字典D的线性组合表示测试图像向量
Figure BDA00002305250400071
得到的表达式为:y=Ax+Dβ+z,式中,x和β分别为线性的稀疏系数,z为噪声项,且x∈Rn×1,β∈Rp×1
Figure BDA00002305250400072
该步骤中,当测试图像的维数d2与训练图像的维数d1不相等时,使用包括插值法的图像缩放方法,使得测试图像的维数d2与训练图像的维数d1相等。
步骤5,因稀疏系数x和β具有多值性,故采用基于L1范数最小化方法求解得到其最优稀疏系数
Figure BDA00002305250400073
Figure BDA00002305250400074
其目标函数为: x ^ β ^ = arg min | | x β | | l , 且满足 | | [ A , D ] x β - y | | 2 ≤ ϵ ; 式中,ε为预设的误差,且最优稀疏系数 x ^ = [ x 1 , x 2 , · · · x t · · · , x n ] ∈ R n × l
Figure BDA00002305250400078
其中,自然数下标t为最优稀疏系数序号,其最大值为n,该
Figure BDA00002305250400079
中的第t个元素xt是重构测试图像向量y时第t个训练图像的系数。
步骤6,分别保留最优稀疏系数
Figure BDA000023052504000710
中与第i类训练图像相对应的系数,并将
Figure BDA000023052504000711
中与第i类训练图像无关的系数全部置零,得到下述带通系数
δ 1 ( x ^ ) = x 1 , x 2 , · · · , x n 1 , 0 , · · · , 0 T ,
δ i ( x ^ ) = 0 , · · · , 0 , x n 1 + n 2 + · · · + n i - 1 + 1 , x n 1 + n 2 + · · · + n i - 1 + 2 , · · · , x n 1 + n 2 + · · · + n i - 1 + n i , 0 , · · · , 0 T ,
……
δ k ( x ^ ) = 0 , · · · , 0 , x n - n k + 1 , x n - n k + 2 , · · · , x n T ;
步骤7,根据步骤6获得的带通系数,分别用每个类别的训练图像和类内变化字典对测试图像y进行重构,得到相应的的重构图像向量qi:
Figure BDA000023052504000716
式中i=1,2,...,k;
步骤8,分别计算原测试图像向量y与重构图像向量qi的残差ri(y)、即两者差值的模:ri(y)=‖y-qi2,式中i=1,2,...,k;残差越小表示重构图像与测试图像越近似或越相像;
步骤9,将该残差ri(y)代入测试图像向量y的类别判定公式:Identity(y)=arg miniri(y),求解得到其中所有k个残差中的最小值,并将该最小残差所对应的类别序号作为最终的人脸识别结果Identity(y)。
本发明已经进行了多次实施试验,下面介绍实施例的仿真试验情况与效果:
参见图3,介绍本发明方法的仿真实施例的试验条件与内容:
仿真实施例的实验使用的人脸图像来自AR人脸数据库,该数据库包含126个人的4000多幅彩色图像,其中每个人的图像都涵盖表情、光照和遮挡的变化,该数据库分为两个部分(session 1和session 2),两者区别是session 2中的图像是在拍摄session 1之后的14天拍摄的。本发明实施例从AR人脸数据库(session1)中随机选择80个人,每个人选择其13幅图像,每幅图像的像素点阵为165×120(参见图3所示)。以每个人在正常表情和光照下的1幅图像作为训练图像,其他12幅图像全部作为测试图像。
仿真实验是首先提取每幅训练图像的灰度特征(19800维)和Gabor特征(10240维)。在获取图像的Gabor特征时,先将原始的165×120图像缩放成128×128像素,而后通过5个尺度加8个方向的40个Gabor滤波器获取响应特征,再将每个滤波器得到的响应特征进行比例为64的下采样,最后得到10240维(即40×16×16)的Gabor特征。然后,从训练图像以外随机选取20个人,每人13幅图像,以此计算类内变化字典。根据具体实施方式步骤3中的第一种方法构建类内变化字典。根据L1范数最小化方法求解测试图像向量y的最优稀疏表示系数
Figure BDA00002305250400081
Figure BDA00002305250400082
最后,分别计算原测试图像向量y和由不同类别训练样本和类内变化字典重构得到的重构图像向量qi之间的残差ri(y)=‖y-qi2,并根据残差大小判定测试图像所归属的类别,即残差越小表示重构图像与测试图像越接近或越相像。本仿真实验的软件平台为MATLAB7.0。
仿真实施例的实验比较了在常用的灰度特征和Gabor特征下,本发明方法和经典的基于稀疏表示SRC(Sparse Representation Classification)的人脸识别方法之间的错误率,实验结果如下表1所示。
表1两种方法在不同图像特征下的识别错误率对比
特征 图像维数  SRC方法识别错误率   本发明方法的识别错误率
灰度 19800  40.31   10.63
Gabor 10240  22.71   5.00
从表1可以看出:本发明方法在两种常用图像特征下的识别错误率均显著低于现有方法。在灰度特征下,本发明方法可将经典SRC方法的错误率降低73.63%,即[(40.31-10.63)/40.31]×100%;在Gabor特征下,本发明方法可将经典SRC方法的错误率降低77.98%,即[(22.71-5.00)/22.71]×100%。可见,本发明方法的改进效果是非常显著的,具有很好的推广应用前景。

Claims (6)

1.一种基于类内变化字典和训练图像的稀疏表示人脸识别方法,其特征在于:所述方法包括下列操作步骤:
(1)对训练图像集合中的所有人脸图像进行特征提取,得到训练图像矩阵A:A=[A1,A2,...,Ai,...,Ak],且
Figure FDA00002305250300011
并对该训练图像矩阵A中的每一列进行归一化处理,使每一列向量的模都为1;式中,自然数下标i为人脸类别序号,其最大值为k;d1和n分别为训练图像的维数和个数,且第i类训练人脸图像矩阵
Figure FDA00002305250300012
其中的ni为第i类训练人脸图像的个数,且
Figure FDA00002305250300013
R为实数集;
(2)对测试图像进行特征提取,得到测试图像向量
Figure FDA00002305250300014
并对该测试图像向量y进行归一化处理;式中,d2为测试图像的维数;
(3)从训练图像集合或其他开放人脸图像集合中获取用于描述相同类别人脸在不同环境因素下的类内变化,构成类内变化字典,即矩阵
Figure FDA00002305250300015
式中,d3为类内变化字典的行数,且d3与d2、d1的数值都相等;p为类内变化字典中描述设定的人脸图像变化的列数;
(4)用训练图像矩阵A和类内变化字典D线性表示测试图像向量
Figure FDA00002305250300016
其表达式为:y=Ax+Dβ+z,式中,x和β分别为线性的稀疏系数,z为噪声项,且x∈Rn×1,β∈Rp×1
Figure FDA00002305250300017
(5)因稀疏系数x和β具有多值性,故采用基于L1范数最小化方法求解得到其最优稀疏系数
Figure FDA00002305250300018
Figure FDA00002305250300019
其目标函数为: x ^ β ^ = arg min | | x β | | l , 且满足 | | [ A , D ] x β - y | | 2 ≤ ϵ ; 式中,ε为预设的误差,且最优稀疏系数 x ^ = [ x 1 , x 2 , · · · x t · · · , x n ] ∈ R n × l
Figure FDA000023052503000113
其中,自然数下标t为最优稀疏系数序号,其最大值为n,该
Figure FDA000023052503000114
中的第t个元素xt是重构测试图像向量y时第t个训练图像的系数;
(6)分别保留最优稀疏系数
Figure FDA000023052503000115
中与第i类训练图像相对应的系数,并将
Figure FDA000023052503000116
中与第i类训练图像无关的系数全部置零,得到下述带通系数
Figure FDA000023052503000117
δ 1 ( x ^ ) = x 1 , x 2 , · · · , x n 1 , 0 , · · · , 0 T ,
δ i ( x ^ ) = 0 , · · · , 0 , x n 1 + n 2 + · · · + n i - 1 + 1 , x n 1 + n 2 + · · · + n i - 1 + 2 , · · · , x n 1 + n 2 + · · · + n i - 1 + n i , 0 , · · · , 0 T ,
……
δ k ( x ^ ) = 0 , · · · , 0 , x n - n k + 1 , x n - n k + 2 , · · · , x n T ;
(7)根据步骤(6)获得的带通系数,分别用每个类别的训练图像和类内变化字典对测试图像y进行重构,得到相应的重构图像向量qi:
Figure FDA00002305250300022
式中i=1,2,...,k;
(8)分别计算原测试图像向量y与重构图像向量qi的残差ri(y)、即两者差值的模:ri(y)=‖y-qi2,式中i=1,2,...,k;残差越小表示重构图像与测试图像越近似或越相像;
(9)将该残差ri(y)代入测试图像向量y的类别判定公式:Identity(y)=arg miniri(y),求解得到其中所有k个残差中的最小值,并将该最小残差所对应的类别序号作为最终的人脸识别结果Identity(y)。
2.根据权利要求1所述的方法,其特征在于:所述类内变化是指同一人脸在不同外部环境下所呈现的差异;所述类内变化字典D是由描述同一人脸在包括不同光照和/或遮挡、和/或表情的干扰因素下的各种图像纹理差异的集合,类内变化字典中的每一列描述一种设定的人脸差异变化,也就是将某一个人脸在异常光照和/或遮挡、和/或表情的图像与其在正常光照和表情、且无遮挡的图像相减得到的差值。
3.根据权利要求1所述的方法,其特征在于:所述类内变化字典的第一种获取方法是求解得到用不同光照和/或遮挡、和/或表情的某个人脸图像减去其在所有环境中的平均图像所得到的差值向量的集合,其操作步骤如下:
(A1)从人脸数据库中随机选择与训练人脸不同的m个人,再对每个人选取w幅人脸图像,用于计算类内变化字典;还对每幅人脸图像进行特征提取,将每幅人脸图像表示为一个d3维的列向量;
(A2)分别将每个人的w幅人脸图像标记为矩阵
Figure FDA00002305250300023
再求解每个人的w幅图像的平均人脸图像式中,自然数下标j为计算类内变化字典时的人脸类别序号,其最大值为m;矩阵Dj的第b列向量自然数下标b为矩阵Dj的列序号,其最大值为w;
(A3)从m个人的所有人脸图像中,计算得到人脸的类内变化字典D: D = [ D 1 - c 1 × e 1 , D 2 - c 2 × e 2 , · · · , D j - c j × e j , · · · , D m - c m × e m ] ∈ R d 3 × ( w × m ) ; 式中,cj为列向量,ej为全1的行向量,且e1=e2=…=ej=…=em=[1,...,1]∈R1×w,cj×ej是每列均为cj且维数为d3×w的矩阵,(Dj-cj×ej)表示矩阵Dj的每一列都减去其对应的平均人脸图像cj
(A4)采用上述方法得到类内变化字典的列数p=w×m;
(A5)对类内变化字典D中的每列元素进行归一化处理,使得每列元素的模值都为1。
4.根据权利要求1所述的方法,其特征在于:所述类内变化字典的第二种获取方法是求解得到同一人脸在包括不同光照和/或遮挡、和/或表情的各种环境下获取的不同图像之间每两个图像的差值向量的集合,其操作步骤如下:
(B1)从人脸数据库中随机选择与训练人脸不同的m个人,再对每个人选取w幅人脸图像,用于计算类内变化字典;还对每幅人脸图像进行特征提取,将每幅人脸图像表示成一个d3维的列向量;
(B2)从每个人的w幅人脸图像中,任意选取其中两幅人脸图像,并计算这两幅人脸图像的差值向量,共完成
Figure FDA00002305250300031
种选取和计算后,将所有差值向量组合成差异矩阵Qj: Q j = [ a j 1 - a j 2 , a j 1 - a j 3 , · · · , a js - a j ( s + 1 ) , · · · , a j ( w - 1 ) - a jw ] ∈ R d 3 × [ w ( w - 1 ) / 2 ] ; 式中,自然数下标j为计算类内变化字典时的人脸类别序号,其最大值为m;第j个人脸图像集合中的第s幅图像
Figure FDA00002305250300033
且自然数序号s的最大值为w;
(B3)从m个人的所有人脸图像中,计算得到人脸的类内变化字典D:
D = [ Q 1 , Q 2 , . . . , Q j , . . . , Q m ] ∈ R d 3 × { m × [ w ( w - 1 ) / 2 ] } ;
(B4)采用上述方法得到类内变化字典的列数
Figure FDA00002305250300035
(B5)对类内变化字典D中的每列元素进行归一化处理,使得每列元素的模值都为1。
5.根据权利要求1所述的方法,其特征在于:所述步骤(3)中,使用其他开放人脸图像集合获取类内变化字典时,若发生该人脸图像集合中的图像维数与训练图像维数d1不相等的情况,则采用包括插值法的图像缩放方法,使得其他开放人脸图像集合中图像的维数等于训练图像的维数d1,以保证类内变化字典的行数d3与训练图像的维数d1相等。
6.根据权利要求1所述的方法,其特征在于:所述步骤(4)中,当测试图像的维数d2与训练图像的维数d1不相等时,使用包括插值法的图像缩放方法,使得测试图像的维数d2与训练图像的维数d1相等。
CN201210411771.0A 2012-10-25 2012-10-25 基于类内变化字典和训练图像的稀疏表示人脸识别方法 Expired - Fee Related CN102915436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210411771.0A CN102915436B (zh) 2012-10-25 2012-10-25 基于类内变化字典和训练图像的稀疏表示人脸识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210411771.0A CN102915436B (zh) 2012-10-25 2012-10-25 基于类内变化字典和训练图像的稀疏表示人脸识别方法

Publications (2)

Publication Number Publication Date
CN102915436A true CN102915436A (zh) 2013-02-06
CN102915436B CN102915436B (zh) 2015-04-15

Family

ID=47613797

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210411771.0A Expired - Fee Related CN102915436B (zh) 2012-10-25 2012-10-25 基于类内变化字典和训练图像的稀疏表示人脸识别方法

Country Status (1)

Country Link
CN (1) CN102915436B (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310415A (zh) * 2013-03-15 2013-09-18 清华大学 基于人脸的图像缺损补绘方法及系统
CN103324914A (zh) * 2013-05-31 2013-09-25 长安大学 一种基于稀疏系数的人脸图像多表情相互转换方法
CN104166860A (zh) * 2014-07-25 2014-11-26 哈尔滨工业大学深圳研究生院 基于约束的面向单个测试样本的人脸识别方法
CN103984918B (zh) * 2014-04-21 2015-06-10 郑州轻工业学院 一种基于类内及类间变化的人脸图像识别方法
CN104732186A (zh) * 2013-12-18 2015-06-24 南京理工大学 基于局部子空间稀疏表示的单样本人脸识别方法
CN104063714B (zh) * 2014-07-20 2016-05-18 詹曙 一种用于视频监控的、基于cuda并行计算与稀疏表示的快速人脸识别算法
CN105760845A (zh) * 2016-02-29 2016-07-13 南京航空航天大学 一种基于联合表示分类的集体人脸识别方法
CN105844223A (zh) * 2016-03-18 2016-08-10 常州大学 一种类别特色字典和共享字典学习的人脸表示算法
CN105957049A (zh) * 2016-02-03 2016-09-21 北京化工大学 一种基于稀疏表示分类的遥感图像变化检测方法
CN106127185A (zh) * 2016-07-01 2016-11-16 南京邮电大学 一种非限制环境下基于低秩协同的人脸性别识别方法
CN106372595A (zh) * 2016-08-31 2017-02-01 重庆大学 遮挡人脸识别方法及装置
CN106462773A (zh) * 2014-04-16 2017-02-22 工业研究与信息服务处 使用gabor函数的模式识别系统及方法
CN106845362A (zh) * 2016-12-27 2017-06-13 湖南长城信息金融设备有限责任公司 一种基于多尺度变换的稀疏表示的人脸识别方法
CN106980825A (zh) * 2017-03-15 2017-07-25 广东顺德中山大学卡内基梅隆大学国际联合研究院 一种基于归一化像素差特征的人脸姿势分类方法
CN107194378A (zh) * 2017-06-28 2017-09-22 深圳大学 一种基于混合字典学习的人脸识别方法及装置
CN107368803A (zh) * 2017-07-14 2017-11-21 广州智慧城市发展研究院 一种基于类别稀疏表示的人脸识别方法及系统
CN107578029A (zh) * 2017-09-21 2018-01-12 北京邮电大学 计算机辅助图片认证的方法、装置、电子设备及存储介质
CN108052932A (zh) * 2018-01-10 2018-05-18 重庆邮电大学 一种遮挡自适应人脸识别方法
CN108828402A (zh) * 2018-04-26 2018-11-16 广东电网有限责任公司 一种高压电缆故障诊断的方法
CN110210336A (zh) * 2019-05-16 2019-09-06 赣南师范大学 一种低分辨率单样本人脸识别方法
CN111931665A (zh) * 2020-08-13 2020-11-13 重庆邮电大学 一种基于类内变化字典建模的欠采样人脸识别方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030169908A1 (en) * 2002-03-04 2003-09-11 Samsung Electronics Co., Ltd. Method and apparatus of recognizing face using component-based 2nd-order principal component analysis (PCA)/independent component analysis (ICA)
CN101976352A (zh) * 2010-10-29 2011-02-16 上海交通大学 基于小样本学习和稀疏表示的多光照人脸识别方法
CN101976360A (zh) * 2010-10-27 2011-02-16 西安电子科技大学 基于多级分类的稀疏表征人脸识别方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030169908A1 (en) * 2002-03-04 2003-09-11 Samsung Electronics Co., Ltd. Method and apparatus of recognizing face using component-based 2nd-order principal component analysis (PCA)/independent component analysis (ICA)
CN101976360A (zh) * 2010-10-27 2011-02-16 西安电子科技大学 基于多级分类的稀疏表征人脸识别方法
CN101976352A (zh) * 2010-10-29 2011-02-16 上海交通大学 基于小样本学习和稀疏表示的多光照人脸识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
WEIHONG DENG 等: "Extended SRC:Undersampled Face Recognition via Intraclass Variant Dictionary", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》, 30 September 2012 (2012-09-30), pages 1864 - 1870, XP011490713, DOI: doi:10.1109/TPAMI.2012.30 *

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310415A (zh) * 2013-03-15 2013-09-18 清华大学 基于人脸的图像缺损补绘方法及系统
CN103324914A (zh) * 2013-05-31 2013-09-25 长安大学 一种基于稀疏系数的人脸图像多表情相互转换方法
CN103324914B (zh) * 2013-05-31 2016-01-06 长安大学 一种基于稀疏系数的人脸图像多表情相互转换方法
CN104732186A (zh) * 2013-12-18 2015-06-24 南京理工大学 基于局部子空间稀疏表示的单样本人脸识别方法
CN104732186B (zh) * 2013-12-18 2018-03-09 南京理工大学 基于局部子空间稀疏表示的单样本人脸识别方法
CN106462773B (zh) * 2014-04-16 2019-07-12 工业研究与信息服务处 使用gabor函数的模式识别系统及方法
CN106462773A (zh) * 2014-04-16 2017-02-22 工业研究与信息服务处 使用gabor函数的模式识别系统及方法
CN103984918B (zh) * 2014-04-21 2015-06-10 郑州轻工业学院 一种基于类内及类间变化的人脸图像识别方法
CN104063714B (zh) * 2014-07-20 2016-05-18 詹曙 一种用于视频监控的、基于cuda并行计算与稀疏表示的快速人脸识别算法
CN104166860B (zh) * 2014-07-25 2017-06-13 哈尔滨工业大学深圳研究生院 基于约束的面向单个测试样本的人脸识别方法
CN104166860A (zh) * 2014-07-25 2014-11-26 哈尔滨工业大学深圳研究生院 基于约束的面向单个测试样本的人脸识别方法
CN105957049A (zh) * 2016-02-03 2016-09-21 北京化工大学 一种基于稀疏表示分类的遥感图像变化检测方法
CN105957049B (zh) * 2016-02-03 2018-10-23 北京化工大学 一种基于稀疏表示分类的遥感图像变化检测方法
CN105760845A (zh) * 2016-02-29 2016-07-13 南京航空航天大学 一种基于联合表示分类的集体人脸识别方法
CN105844223A (zh) * 2016-03-18 2016-08-10 常州大学 一种类别特色字典和共享字典学习的人脸表示算法
CN106127185A (zh) * 2016-07-01 2016-11-16 南京邮电大学 一种非限制环境下基于低秩协同的人脸性别识别方法
CN106372595A (zh) * 2016-08-31 2017-02-01 重庆大学 遮挡人脸识别方法及装置
CN106845362A (zh) * 2016-12-27 2017-06-13 湖南长城信息金融设备有限责任公司 一种基于多尺度变换的稀疏表示的人脸识别方法
CN106845362B (zh) * 2016-12-27 2020-03-20 湖南长城信息金融设备有限责任公司 一种基于多尺度变换的稀疏表示的人脸识别方法
CN106980825A (zh) * 2017-03-15 2017-07-25 广东顺德中山大学卡内基梅隆大学国际联合研究院 一种基于归一化像素差特征的人脸姿势分类方法
CN107194378A (zh) * 2017-06-28 2017-09-22 深圳大学 一种基于混合字典学习的人脸识别方法及装置
CN107194378B (zh) * 2017-06-28 2020-11-17 深圳大学 一种基于混合字典学习的人脸识别方法及装置
CN107368803A (zh) * 2017-07-14 2017-11-21 广州智慧城市发展研究院 一种基于类别稀疏表示的人脸识别方法及系统
CN107578029B (zh) * 2017-09-21 2020-03-27 北京邮电大学 计算机辅助图片认证的方法、装置、电子设备及存储介质
CN107578029A (zh) * 2017-09-21 2018-01-12 北京邮电大学 计算机辅助图片认证的方法、装置、电子设备及存储介质
CN108052932A (zh) * 2018-01-10 2018-05-18 重庆邮电大学 一种遮挡自适应人脸识别方法
CN108828402A (zh) * 2018-04-26 2018-11-16 广东电网有限责任公司 一种高压电缆故障诊断的方法
CN110210336A (zh) * 2019-05-16 2019-09-06 赣南师范大学 一种低分辨率单样本人脸识别方法
CN110210336B (zh) * 2019-05-16 2021-03-16 赣南师范大学 一种低分辨率单样本人脸识别方法
CN111931665A (zh) * 2020-08-13 2020-11-13 重庆邮电大学 一种基于类内变化字典建模的欠采样人脸识别方法
CN111931665B (zh) * 2020-08-13 2023-02-21 重庆邮电大学 一种基于类内变化字典建模的欠采样人脸识别方法

Also Published As

Publication number Publication date
CN102915436B (zh) 2015-04-15

Similar Documents

Publication Publication Date Title
CN102915436B (zh) 基于类内变化字典和训练图像的稀疏表示人脸识别方法
Han et al. Autoencoder inspired unsupervised feature selection
William et al. Face recognition using facenet (survey, performance test, and comparison)
CN108491849B (zh) 基于三维稠密连接卷积神经网络的高光谱图像分类方法
Svoboda et al. Palmprint recognition via discriminative index learning
Singh et al. Gender and ethnicity classification of iris images using deep class-encoder
Dong et al. Finger vein recognition based on multi-orientation weighted symmetric local graph structure
Su et al. Nonnegative sparse autoencoder for robust endmember extraction from remotely sensed hyperspectral images
Kekre et al. Performance Comparison for Face Recognition using PCA, DCT &WalshTransform of Row Mean and Column Mean
CN112836671B (zh) 一种基于最大化比率和线性判别分析的数据降维方法
CN107341505B (zh) 一种基于图像显著性与Object Bank的场景分类方法
Wang et al. Multi-scale feature extraction algorithm of ear image
Ouarda et al. MLP Neural Network for face recognition based on Gabor Features and Dimensionality Reduction techniques
CN104239856A (zh) 基于Gabor特征和自适应线性回归的人脸识别方法
CN108664941B (zh) 基于测地映射分析的核稀疏描述人脸识别方法
Nimbarte et al. Age Invariant Face Recognition using Convolutional Neural Network.
Jalled Face recognition machine vision system using eigenfaces
CN103745242A (zh) 一种跨设备生物特征识别方法
CN113095305B (zh) 一种医药异物高光谱分类检测方法
Wu et al. Age-invariant face recognition using coupled similarity reference coding
CN103246877A (zh) 基于图像轮廓的人脸识别新方法
Lu et al. Fusing shape and texture information for facial age estimation
Guo et al. Palmprint Recognition Based on Local Fisher Discriminant Analysis.
CN110287973B (zh) 一种基于低秩鲁棒线性鉴别分析的图像特征提取方法
CN101482917B (zh) 一种基于二阶二维主分量分析的人脸识别系统和方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150415

Termination date: 20211025

CF01 Termination of patent right due to non-payment of annual fee