CN106295609B - 基于块稀疏结构低秩表示的单样本人脸识别方法 - Google Patents
基于块稀疏结构低秩表示的单样本人脸识别方法 Download PDFInfo
- Publication number
- CN106295609B CN106295609B CN201610701068.1A CN201610701068A CN106295609B CN 106295609 B CN106295609 B CN 106295609B CN 201610701068 A CN201610701068 A CN 201610701068A CN 106295609 B CN106295609 B CN 106295609B
- Authority
- CN
- China
- Prior art keywords
- pixel
- block
- rank representation
- low
- sub
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012360 testing method Methods 0.000 claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 35
- 238000012549 training Methods 0.000 claims abstract description 25
- 239000013598 vector Substances 0.000 claims description 8
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims description 6
- 230000000717 retained effect Effects 0.000 claims description 5
- 238000000638 solvent extraction Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000005286 illumination Methods 0.000 abstract description 6
- 230000003416 augmentation Effects 0.000 abstract 1
- 230000011218 segmentation Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000004438 eyesight Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000001149 cognitive effect Effects 0.000 description 2
- 238000003909 pattern recognition Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000016776 visual perception Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/16—Human faces, e.g. facial parts, sketches or expressions
- G06V40/172—Classification, e.g. identification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2136—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on sparsity criteria, e.g. with an overcomplete basis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了基于块稀疏结构低秩表示的单样本人脸识别方法,该方法首先将人脸划分为多个块,每个块又划分为多个重叠的子块,假设同一块内的子块在同一子空间内;然后基于低秩表示模型用全部训练样本对应块内的所有子块构成的局部字典去低秩表示所有测试图像的对应块的中心子块所构成的测试矩阵以实现每个人所对应的子空间的有效分割,并加入块稀疏约束以增强模型的鉴别性,利用非严格的增广拉格朗日乘法求解模型以获得低秩表示系数矩阵;在此基础通过判定表示系数的大小实现对测试图像块的分类;最后对所有测试图像块进行投票最终确定分类结果。本发明不仅对表情、光照变化和遮挡等具有很好的鲁棒性,识别精度高,而且支持高效的并行计算。
Description
技术领域
本发明涉及单样本人脸识别方法,特别是涉及每个待识别对象仅有一幅训练图像的基于块稀疏结构低秩表示的单样本人脸识别方法,属于人脸识别技术领域。
背景技术
自动人脸识别是利用计算机分析人脸图像或视频,实现辨认身份的一种技术。作为一种最直接、最自然、最容易被人接受的生物特征识别技术,自动人脸识别一直是模式识别和计算机视觉领域最热门的研究课题之一,其在理论和实际应用方面都有非常重要的研究意义。在理论方面,自动人脸识别问题的深入研究和有效解决,可以极大地促进模式识别、图像处理与分析、计算机视觉、神经计算和认知心理学等相关学科的发展。在实际应用方面,人脸识别技术在公共安全和日常生活中有着广泛的应用,特别是“9.11”恐怖袭击事件之后,恐怖组织直接威胁到国家安全和平常百姓的生命安全,人脸识别技术作为一种简单、有效的反恐手段而受到前所未有的重视,被广泛应用于商场、地铁和机场等公共场所。
在众多已有的人脸识别方法中,鉴别特征的提取是人脸识别的关键。过去几十年中,在众多的特征提取算法中子空间学习方法受到了广泛的关注。其基本思想是通过一定的准则函数寻找一组基向量,把原始的高维特征投影到这组基向量所张成的低维子空间中,从而使得数据更加紧凑更具有鉴别性,另外也可以降低计算量。子空间学习方法中最具代表性的当属主成分分析(M.A.Turk,A.P.Pentland,“Eigenfaces for recognition,”Journal of cognitive neuroscience,1991,3(1):71-86)和线性鉴别分析(P.N.Belhumeur,J.P.Hespanha.Kriegman,D.J.“Eigenfaces vs.Fisherfaces:Recognition using class specific linear projection,”IEEE Transactions onPattern Analysis and Machine Intelligence,1997,19(7):711-720),其已经被广泛地应用于人脸识别中并取得了很好的识别结果。近年来J.Wright等人提出的基于稀疏表示的人脸识别方法(A.Yang,A.Ganesh,S.Sastry,and Y.Ma,“Robust Face Recognition viaSparse Representation,”IEEE Trans.Pattern Analysis and Machine Intelligence,vol.31,no.2,pp.210-227,2009)使得特征空间的选择不再那么重要,其基本原理是利用所有训练图像构造字典,再通过求解一个欠定方程组来求得测试图像的最稀疏线性组合系数,然后根据这些系数来对图像进行识别分类。
在实际生活中,例如护照验证、安全监控、身份证验证通常一个人只有一个样本,这种情况下的人脸识别问题我们称之为“单样本”问题。在这种情况下,很多传统的人脸识别方法以及基于稀疏或协同表示的人脸识别方法性能下降严重甚至无法工作。这主要是因为这些传统的方法对单样本情况下产生的野值点或遮挡、光照等稀疏噪声不够鲁棒。
发明内容
本发明所要解决的技术问题是:提供基于块稀疏结构低秩表示的单样本人脸识别方法,为单样本人脸识别问题提供了一种简单有效的解决方案。
本发明为解决上述技术问题采用以下技术方案:
基于块稀疏结构低秩表示的单样本人脸识别方法,包括如下步骤:
步骤1,用大小相等的正方形窗口对所有类的训练人脸图像和测试人脸图像进行分块,正方形窗口滑过图像的每个像素点,在图像的每一个像素点处得到对应的一个正方形块,接着在每个正方形块内部,再确定出一个中心子块及其若干邻域子块;
步骤2,假设每个正方形块内的不同子块属于同一线性子空间,利用所有类的训练人脸图像像素i处对应的正方形块内的全部子块构建局部字典Ai,所有类的测试人脸图像像素i处对应的正方形块内的全部子块构建测试矩阵Xi,使用低秩表示模型描述不同测试人脸图像所对应的子空间;
步骤3,在步骤2的低秩表示模型中加入块稀疏约束,得到新的低秩表示模型;
步骤4,通过非严格的增广拉格朗日乘法求解步骤3新的低秩表示模型,得到低秩表示系数矩阵;
步骤5,根据低秩表示系数矩阵,对每个测试人脸图像的各正方形块进行分类;
步骤6,识别每个测试人脸图像的各正方形块的类别后,通过投票的方法确定测试人脸图像所属的类别。
作为本发明的一种优选方案,所述步骤1的具体过程为:
1.1定义图像中任意一个像素i为中心、半径为R的正方形上的像素为该像素的邻域像素,则像素i的领域像素集表示为每个邻域像素对应一个以该邻域像素为中心的S×S子块,S为大于等于3的奇数,子块中的中的S2个像素表示成向量形式为同样地,中心像素i也对应一个S×S子块,表示成向量形式为中心像素i及其邻域像素对应的所有子块形成一个以像素i为中心、大小为(S+2R)×(S+2R)的正方形块,P表示邻域像素的总个数;
1.2采用边缘像素境像的方法来处理图像的边缘像素,得到以图像的边缘像素为中心的正方形块。
作为本发明的一种优选方案,所述步骤2的具体过程为:
2.1局部字典其中表示第k个或类训练人脸图像在像素i处对应的正方形块内的子块构成的字典,k=1,2,…,K,K表示训练人脸图像的总个数或总类别数,P表示像素i的邻域像素的总个数;
2.2测试矩阵其中表示第j个测试人脸图像在像素i处对应的正方形块内的子块构成的矩阵,j=1,2,…,M,M表示测试人脸图像的总个数,P表示像素i的邻域像素的总个数;
2.3使用低秩表示模型描述不同测试人脸图像所对应的子空间:
其中,Zi表示低秩表示系数矩阵,Ei表示误差,λ表示||Zi||*的正则项系数,||·||F表示F范数,||·||*表示核范数。
作为本发明的一种优选方案,步骤3所述新的低秩表示模型为:
其中,Zi表示低秩表示系数矩阵,Ei表示误差,λ表示||Zi||*的正则项系数,β表示∑j∑k||δk,j(Zi)||F的正则项系数,δk,j表示映射函数,δk,j(Zi)表示将Zi中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0,||·||F表示F范数,||·||*表示核范数,Xi表示测试矩阵,Ai表示局部字典。
作为本发明的一种优选方案,所述步骤4的具体过程为:
4.1将步骤3新的低秩表示模型转化为如下形式,忽略像素i:
s.t.X=AZ+E,Z=H,Z=Q
其中,Z表示低秩表示系数矩阵,E表示误差,λ表示||H||*的正则项系数,β表示∑j∑k||δk,j(Q)||F的正则项系数,δk,j(Q)表示将Q中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0,||·||F表示F范数,||·||*表示核范数,X表示测试矩阵,A表示局部字典;
4.2将4.1的模型转换为增广拉格朗日函数L:
其中,Tr(·)表示矩阵的迹,Y1、Y2、Y3表示拉格朗日乘子,T表示转置,μ为大于0的惩罚参数;
4.3Z、H、Q、E、Y1、Y2、Y3初始化为0,μ=10-6,固定其他未知数,更新H,则模型的目标函数可以转化为如下求H的目标函数:
该目标函数可以通过SVT算子求得闭合解;
4.4获得4.3中H的最优解后,固定其他未知数,更新Q,则模型可以转化为如下求Q的目标函数:
令W=Z+Y3/μ,上述求Q的目标函数可以转化为如下形式:
对于任意的k、j,
4.5固定其他未知数,通过如下方式更新Z:
Z=(ATA+2I)-1(ATX-ATE+H+Q+(ATY1-Y2-Y3)/μ)
4.6固定其他未知数,通过如下方式更新E:
E=(μ(X-AZ)+Y1)/(1+μ)
4.7更新拉格朗日乘子:
Y1=Y1+μ(X-AZ-E)
Y2=Y2+μ(Z-H)
Y3=Y3+μ(Z-Q)
4.8更新参数μ:μ=min(ρμ,μmax),其中,μmax=1010,ρ=1;
4.9检查是否满足以下收敛条件:
||X-AZ-E||∞<ε
||Z-H||∞<ε
||Z-Q||∞<ε
其中,ε=10-8,当不满足上述收敛条件时,重复4.3-4.9,直至满足收敛条件,得到低秩表示系数矩阵。
作为本发明的一种优选方案,所述步骤5的分类公式为:
其中,表示第j个测试人脸图像的像素i对应的正方形块所属的类别,δk,j(Zi)表示将Zi中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0,||·||F表示F范数。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1、本发明一方面由于低秩表示对表情、光照变化和遮挡等具有很好的鲁棒性,因而具有较高的识别精度,另一方面源于加入的块稀疏约束增强了低秩模型的判别能力。
2、本发明无需特征提取,分类简单易行;且充分利用了样本的互补性和相关性。
3、本发明判别图像块类别的工作可以并行处理,因而可以大大节省计算时间,提供计算效率。
附图说明
图1是本发明基于块稀疏结构低秩表示的单样本人脸识别方法的算法流程图。
图2是本发明基于块稀疏结构低秩表示的单样本人脸识别方法的像素邻域示意图;其中,(a)P=8,R=1;(b)P=16,R=2;(c)P=24,R=3。
图3是本发明基于块稀疏结构低秩表示的单样本人脸识别方法的子块关系示意图,以P=8,R=1为例;其中,(a)为中心像素i0的子块,(b)为邻域像素i1的子块,(c)为邻域像素i5的子块。
图4是本发明基于块稀疏结构低秩表示的单样本人脸识别方法的块与子块两级划分示意图,以P=16,R=2为例。
图5是本发明基于块稀疏结构低秩表示的单样本人脸识别方法块稀疏结构图。
图6是本发明基于块稀疏结构低秩表示的单样本人脸识别方法在Extended YaleB人脸库上的性能效果对比图。
图7是本发明基于块稀疏结构低秩表示的单样本人脸识别方法在AR人脸库Session1上的性能效果对比图。
图8是本发明基于块稀疏结构低秩表示的单样本人脸识别方法在AR人脸库Session2上的性能效果对比图。
图9是本发明基于块稀疏结构低秩表示的单样本人脸识别方法在LFW人脸库上的性能效果对比图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
近年来,低秩表示LRR(G.Liu,Z.Lin,Robust recovery of subspace structuresby low-rank representation)由于能够实现多个子空间结构的有效分割从而吸引了越来越多的关注。它试图揭示隐含在高维空间的成员关系。由于低秩约束对野值和各种图像变化更加鲁棒,因而,LRR也应该可以用来解决单样本问题。此外,低秩表示模型利用样本间的相关性和互补性从全体数据中捕捉本质特征,并且符合人类视觉感知的稀疏编码策略,因此理论上适合用于提升单样本人脸识别性能。基于这一想法,本发明提出一种基于块稀疏结构低秩表示的单样本人脸识别方法。
结合图1所示,本发明基于块稀疏结构低秩表示的单样本人脸识别方法,包括以下步骤:
1、将所有单训练样本和测试样本的人脸图像划分为多个重叠的块,每个块中心对应图像的一个像素点;接着在每一块内部,再划分出一个中心子块和周边的若干邻域子块,提出人脸图像的每一个块内的不同子块属于同一线性子空间的合理假设,具体如下:
(1)将所有单训练样本和测试样本的人脸图像划分为多个重叠的块,每个块中心对应图像的一个像素点,以该像素点为中心、半径为R的正方形上的像素为该像素的邻域像素,则像素i的邻域像素集可以表示为如图2的(a)、(b)、(c)所示,分别为不同P,R情况下的邻域集。邻域集中的每个像素对应一个以其为中心的S×S子块(S取大于等于3的奇数),子块中的S2个像素表示成向量形式为同样地,中心像素i也对应一个S×S的子块,表示成向量形式为中心像素i及其邻域像素对应的所有小块形成了一个以像素i为中心的大块,大小为(S+2R)×(S+2R),如图3的(a)、(b)、(c)所示。对所有训练人脸图像和测试人脸图像使用上述固定大小的正方形大块进行分块,大块中心对应图像中的像素点,每个大块由图3所示的一个中心子块和周边的若干邻域子块组成。如图4所示,为图像的两级划分示意图。
(2)对于图像边缘像素,采用边缘像素境像的方法来处理因大块超出图像边缘而导致其内部部分像素值缺失的现象,得到以图像边缘像素为中心的大块。
(3)由于大块内的中心子块与邻域子块具有极强的相似性,因此可以假设其属于同一线性子空间。
2、一个单训练样本对应一个类,使用所有类的训练人脸图像的像素i处对应的块内的全部子块构成局部字典Ai,所有类的测试人脸图像的像素i处对应的块内的全部子块构成测试矩阵Xi;使用如下的低秩表示模型来描述不同人脸所对应的子空间:
其中,Zi表示低秩表示系数矩阵,Ei表示误差,λ表示||Zi||*的正则项系数。
3、根据1.2所述的低秩表示模型,每张测试图像的局部块内的子块可以由训练图像在相同位置的子块及其邻域子块线性联合表示,所以在理想情况下利用低秩模型获取低维表征系数后,表征系数矩阵会呈现出块稀疏的结构,如图5所示。为了进一步加强系数矩阵的块稀疏性,我们在低秩表达式模型中显式地加入一个块稀疏约束来进一步强化所得系数矩阵的块稀疏性,即:
其中,β表示∑j∑k||δk,j(Zi)||F的正则项系数,δk,j表示映射函数即Rm×n→Rm×n,δk,j(Zi)表示将Zi中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0。
4、通过非严格的增广拉格朗日乘法求解1.3中的低秩模型获得表示系数矩阵Zi,具体步骤如下:
(1)为了方便,我们忽略1.3中低秩模型的上标i,将基于块稀疏结构的低秩模型转化为如下的形式:
s.t.X=AZ+E,Z=H,Z=Q
(2)将(1)中的低秩模型转换为如下的增广拉格朗日函数L:
其中,Tr(·)表示矩阵的迹,Y1、Y2、Y3表示拉格朗日乘子,T表示转置,μ为大于0的惩罚参数。
(3)Z、H、Q、E、Y1、Y2、Y3初始化为0,μ=10-6,μmax=1010,ρ=1,ε=10-8。
(4)首先固定其他未知数,更新H,则模型的目标函数可以转化为如下求H的目标函数:
该目标函数可以通过SVT(Singular Value Thresholding)算子求得闭合解。
(5)获得(4)中H的最优解后,固定其他未知数,更新Q,则模型可以转化为如下求Q的目标函数:
该目标函数可以一块一块地去解,令W=Z+Y3/μ,那么上述求Q的目标函数可以转化为如下形式:
对于任意的k、j,
(6)固定其他未知数,通过如下方式更新Z:
Z=(ATA+2I)-1(ATX-ATE+H+Q+(ATY1-Y2-Y3)/μ)
(7)固定其他未知数,通过如下方式更新E:
E=(μ(X-AZ)+Y1)/(1+μ)
(8)更新拉格朗日乘子:
Y1=Y1+μ(X-AZ-E)
Y2=Y2+μ(Z-H)
Y3=Y3+μ(Z-Q)
(9)更新参数μ:
μ=min(ρμ,μmax)
(10)检查是否满足以下收敛条件:
||X-AZ-E||∞<ε
||Z-H||∞<ε
||Z-Q||∞<ε
若不满足,则重复(4)-(10)的操作,直至满足收敛条件。
5、根据表示系数矩阵Zi直接进行分类,第j个测试样本的第i个区域块的分类结果由下式获得:
其中,表示第j个测试人脸图像的像素i对应的正方形块所属的类别。
6、对每个测试人脸图像块矩阵Xi的低秩表示分类可以并行进行,即可以同时对多个测试人脸图像块矩阵进行步骤1.2、1.3、1.4、1.5操作。
7、识别每个测试人脸块的类别之后,通过投票的方法最终确定测试人脸图像所属的类别。
本发明基于块稀疏结构低秩表示的单样本人脸识别方法,每块的低秩表示分类可以并行操作,图6为在Extended Yale B人脸库上的性能效果对比图,图7为在AR人脸库Session1上性能效果对比图,图8为在AR人脸库Session2上的性能效果对比图,图9为在LFW人脸库上的性能效果对比图。从图6可以看出,本发明方法取得了最佳的识别率。另外图6也充分说明本发明方法对光照变化具有很强的鲁棒性,因为Extended Yale B人脸库中的人脸图像光照变化非常大。从图7和图8可以发现,本发明方法可以在不依赖于额外训练样本集的条件下仍然获得最佳性能,同时对光照、表情、遮挡、时间变化都具有很好的鲁棒性。从图9可以看出,本发明方法在无约束环境下仍然可以获得最佳的性能。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (6)
1.基于块稀疏结构低秩表示的单样本人脸识别方法,其特征在于,包括如下步骤:
步骤1,用大小相等的正方形窗口对所有类的训练人脸图像和测试人脸图像进行分块,正方形窗口滑过图像的每个像素点,在图像的每一个像素点处得到对应的一个正方形块,接着在每个正方形块内部,再确定出一个中心子块及其若干邻域子块;
步骤2,假设每个正方形块内的不同子块属于同一线性子空间,利用所有类的训练人脸图像像素i处对应的正方形块内的全部子块构建局部字典Ai,所有类的测试人脸图像像素i处对应的正方形块内的全部子块构建测试矩阵Xi,使用低秩表示模型描述不同测试人脸图像所对应的子空间;
步骤3,在步骤2的低秩表示模型中加入块稀疏约束,得到新的低秩表示模型;
步骤4,通过非严格的增广拉格朗日乘法求解步骤3新的低秩表示模型,得到低秩表示系数矩阵;
步骤5,根据低秩表示系数矩阵,对每个测试人脸图像的各正方形块进行分类;
步骤6,识别每个测试人脸图像的各正方形块的类别后,通过投票的方法确定测试人脸图像所属的类别。
2.根据权利要求1所述基于块稀疏结构低秩表示的单样本人脸识别方法,其特征在于,所述步骤1的具体过程为:
1.1定义图像中任意一个像素i为中心、半径为R的正方形上的像素为该像素的邻域像素,则像素i的邻域像素集表示为每个邻域像素对应一个以该邻域像素为中心的S×S子块,S为大于等于3的奇数,子块中的中的S2个像素表示成向量形式为同样地,中心像素i也对应一个S×S子块,表示成向量形式为中心像素i及其邻域像素对应的所有子块形成一个以像素i为中心、大小为(S+2R)×(S+2R)的正方形块,P表示邻域像素的总个数;
1.2采用边缘像素镜像的方法来处理图像的边缘像素,得到以图像的边缘像素为中心的正方形块。
3.根据权利要求1所述基于块稀疏结构低秩表示的单样本人脸识别方法,其特征在于,所述步骤2的具体过程为:
2.1局部字典其中表示第k个或类训练人脸图像在像素i处对应的正方形块内的子块构成的字典,k=1,2,…,K,K表示训练人脸图像的总个数或总类别数,P表示像素i的邻域像素的总个数;
2.2测试矩阵其中表示第j个测试人脸图像在像素i处对应的正方形块内的子块构成的矩阵,j=1,2,…,M,M表示测试人脸图像的总个数,P表示像素i的邻域像素的总个数;
2.3使用低秩表示模型描述不同测试人脸图像所对应的子空间:
s.t.Xi=AiZi+Ei
其中,Zi表示低秩表示系数矩阵,Ei表示误差,λ表示||Zi||*的正则项系数,||·||F表示F范数,||·||*表示核范数。
4.根据权利要求1所述基于块稀疏结构低秩表示的单样本人脸识别方法,其特征在于,步骤3所述新的低秩表示模型为:
s.t.Xi=AiZi+Ei
其中,Zi表示低秩表示系数矩阵,Ei表示误差,λ表示||Zi||*的正则项系数,β表示∑j∑k||δk,j(Zi)||F的正则项系数,δk,j表示映射函数,δk,j(Zi)表示将Zi中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0,||·||F表示F范数,||·||*表示核范数,Xi表示测试矩阵,Ai表示局部字典。
5.根据权利要求1所述基于块稀疏结构低秩表示的单样本人脸识别方法,其特征在于,所述步骤4的具体过程为:
4.1将步骤3新的低秩表示模型转化为如下形式,忽略像素i:
s.t.X=AZ+E,Z=H,Z=Q
其中,Z表示低秩表示系数矩阵,E表示误差,λ表示||H||*的正则项系数,β表示∑j∑k||δk,j(Q)||F的正则项系数,δk,j(Q)表示将Q中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0,||·||F表示F范数,||·||*表示核范数,X表示测试矩阵,A表示局部字典;
4.2将4.1的模型转换为增广拉格朗日函数L:
其中,Tr(·)表示矩阵的迹,Y1、Y2、Y3表示拉格朗日乘子,T表示转置,μ为大于0的惩罚参数;
4.3 Z、H、Q、E、Y1、Y2、Y3初始化为0,μ=10-6,固定其他未知数,更新H,则模型的目标函数可以转化为如下求H的目标函数:
该目标函数可以通过SVT算子求得闭合解;
4.4获得4.3中H的最优解后,固定其他未知数,更新Q,则模型可以转化为如下求Q的目标函数:
令W=Z+Y3/μ,上述求Q的目标函数可以转化为如下形式:
对于任意的k、j,
4.5固定其他未知数,通过如下方式更新Z:
Z=(ATA+2I)-1(ATX-ATE+H+Q+(ATY1-Y2-Y3)/μ)
4.6固定其他未知数,通过如下方式更新E:
E=(μ(X-AZ)+Y1)/(1+μ)
4.7更新拉格朗日乘子:
Y1=Y1+μ(X-AZ-E)
Y2=Y2+μ(Z-H)
Y3=Y3+μ(Z-Q)
4.8更新参数μ:μ=min(ρμ,μmax),其中,μmax=1010,ρ=1;
4.9检查是否满足以下收敛条件:
||X-AZ-E||∞<ε
||Z-H||∞<ε
||Z-Q||∞<ε
其中,ε=10-8,当不满足上述收敛条件时,重复4.3-4.9,直至满足收敛条件,得到低秩表示系数矩阵。
6.根据权利要求1所述基于块稀疏结构低秩表示的单样本人脸识别方法,其特征在于,所述步骤5的分类公式为:
其中,表示第j个测试人脸图像的像素i对应的正方形块所属的类别,δk,j(Zi)表示将Zi中与第k类训练人脸图像和第j个测试人脸图像相关的系数保留下来,其余全部置0,||·||F表示F范数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610701068.1A CN106295609B (zh) | 2016-08-22 | 2016-08-22 | 基于块稀疏结构低秩表示的单样本人脸识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610701068.1A CN106295609B (zh) | 2016-08-22 | 2016-08-22 | 基于块稀疏结构低秩表示的单样本人脸识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106295609A CN106295609A (zh) | 2017-01-04 |
CN106295609B true CN106295609B (zh) | 2019-05-10 |
Family
ID=57662073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610701068.1A Active CN106295609B (zh) | 2016-08-22 | 2016-08-22 | 基于块稀疏结构低秩表示的单样本人脸识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106295609B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107992846A (zh) * | 2017-12-15 | 2018-05-04 | 贵阳宏益房地产开发有限公司 | 遮挡人脸识别方法及装置 |
CN108446589B (zh) * | 2018-02-07 | 2022-03-22 | 杭州电子科技大学 | 复杂环境下基于低秩分解和辅助字典的人脸识别方法 |
CN109063555B (zh) * | 2018-06-26 | 2021-07-02 | 杭州电子科技大学 | 基于低秩分解和稀疏表示残差对比的多姿态人脸识别方法 |
CN108875686A (zh) * | 2018-07-02 | 2018-11-23 | 扬州大学 | 一种基于非负稀疏的鉴别性低秩表示分类方法 |
CN109711283B (zh) * | 2018-12-10 | 2022-11-15 | 广东工业大学 | 一种联合双字典和误差矩阵的遮挡表情识别方法 |
CN109740615B (zh) * | 2018-12-29 | 2023-04-07 | 武汉大学 | 一种对抗攻击样本扰动的去除方法 |
CN112580444B (zh) * | 2020-12-02 | 2022-08-05 | 河海大学 | 基于循环自编码和块稀疏结构表示的单样本人脸识别方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010144259A1 (en) * | 2009-06-09 | 2010-12-16 | Arizona Board Of Regents Acting For And On Behalf Of Arizona State University | Ultra-low dimensional representation for face recognition under varying expressions |
CN103632138A (zh) * | 2013-11-20 | 2014-03-12 | 南京信息工程大学 | 一种低秩分块稀疏表示的人脸识别方法 |
CN103942531A (zh) * | 2014-03-06 | 2014-07-23 | 中南民族大学 | 一种人脸识别系统及其方法 |
CN104392246A (zh) * | 2014-12-03 | 2015-03-04 | 北京理工大学 | 一种基于类间类内面部变化字典的单样本人脸识别方法 |
CN104732186A (zh) * | 2013-12-18 | 2015-06-24 | 南京理工大学 | 基于局部子空间稀疏表示的单样本人脸识别方法 |
CN105046276A (zh) * | 2015-07-14 | 2015-11-11 | 西安电子科技大学 | 基于低秩表示的高光谱图像波段选择方法 |
CN105469359A (zh) * | 2015-12-09 | 2016-04-06 | 武汉工程大学 | 基于局部约束低秩表示的人脸超分辨率重建方法 |
-
2016
- 2016-08-22 CN CN201610701068.1A patent/CN106295609B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010144259A1 (en) * | 2009-06-09 | 2010-12-16 | Arizona Board Of Regents Acting For And On Behalf Of Arizona State University | Ultra-low dimensional representation for face recognition under varying expressions |
CN103632138A (zh) * | 2013-11-20 | 2014-03-12 | 南京信息工程大学 | 一种低秩分块稀疏表示的人脸识别方法 |
CN104732186A (zh) * | 2013-12-18 | 2015-06-24 | 南京理工大学 | 基于局部子空间稀疏表示的单样本人脸识别方法 |
CN103942531A (zh) * | 2014-03-06 | 2014-07-23 | 中南民族大学 | 一种人脸识别系统及其方法 |
CN104392246A (zh) * | 2014-12-03 | 2015-03-04 | 北京理工大学 | 一种基于类间类内面部变化字典的单样本人脸识别方法 |
CN105046276A (zh) * | 2015-07-14 | 2015-11-11 | 西安电子科技大学 | 基于低秩表示的高光谱图像波段选择方法 |
CN105469359A (zh) * | 2015-12-09 | 2016-04-06 | 武汉工程大学 | 基于局部约束低秩表示的人脸超分辨率重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106295609A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106295609B (zh) | 基于块稀疏结构低秩表示的单样本人脸识别方法 | |
Rida et al. | Palmprint identification using an ensemble of sparse representations | |
Zhou et al. | Face recognition based on PCA image reconstruction and LDA | |
Fan et al. | Superpixel guided deep-sparse-representation learning for hyperspectral image classification | |
Wu et al. | Lut-based adaboost for gender classification | |
Liu et al. | Learning expressionlets via universal manifold model for dynamic facial expression recognition | |
Wang et al. | Laplacian LRR on product Grassmann manifolds for human activity clustering in multicamera video surveillance | |
Liong et al. | Face recognition using deep PCA | |
Wang et al. | Pyramid-based multi-scale lbp features for face recognition | |
Thepade et al. | Face gender recognition using multi layer perceptron with OTSU segmentation | |
CN110956113B (zh) | 基于二次协作表示鉴别投影的鲁棒人脸识别方法 | |
CN111325275A (zh) | 基于低秩二维局部鉴别图嵌入的鲁棒图像分类方法及装置 | |
Akbar et al. | Face recognition using hybrid feature space in conjunction with support vector machine | |
Devi et al. | Content based feature combination method for face image retrieval using neural network and SVM classifier for face recognition | |
CN110287973B (zh) | 一种基于低秩鲁棒线性鉴别分析的图像特征提取方法 | |
Siyahjani et al. | A supervised low-rank method for learning invariant subspaces | |
Shiau et al. | A sparse representation method with maximum probability of partial ranking for face recognition | |
Nikan et al. | Partial face recognition based on template matching | |
CN107506694B (zh) | 基于局部中值表示的鲁棒人脸识别方法 | |
Naveen et al. | Human emotions detection using kernel nonlinear collaborative discriminant regression classifier: human emotions detection using KNCDRC | |
Jagadeesh et al. | DBC based Face Recognition using DWT | |
Verschae et al. | A hybrid face detector based on an asymmetrical adaboost cascade detector and a wavelet-bayesian-detector | |
Wang et al. | Kernelized low rank representation on Grassmann manifolds | |
CN112580444B (zh) | 基于循环自编码和块稀疏结构表示的单样本人脸识别方法 | |
Kathavarayan et al. | Identification of untrained facial image in combined global and local preserving feature space |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |