CN110781749B - 基于双偏差双空间局部方向模式的人脸识别方法 - Google Patents

基于双偏差双空间局部方向模式的人脸识别方法 Download PDF

Info

Publication number
CN110781749B
CN110781749B CN201910911208.1A CN201910911208A CN110781749B CN 110781749 B CN110781749 B CN 110781749B CN 201910911208 A CN201910911208 A CN 201910911208A CN 110781749 B CN110781749 B CN 110781749B
Authority
CN
China
Prior art keywords
image
double
deviation
space
face
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
Application number
CN201910911208.1A
Other languages
English (en)
Other versions
CN110781749A (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201910911208.1A priority Critical patent/CN110781749B/zh
Publication of CN110781749A publication Critical patent/CN110781749A/zh
Application granted granted Critical
Publication of CN110781749B publication Critical patent/CN110781749B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/172Classification, e.g. identification
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2413Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on distances to training or reference patterns
    • G06F18/24147Distances to closest patterns, e.g. nearest neighbour classification
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/168Feature extraction; Face representation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Computer Interaction (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)
  • Collating Specific Patterns (AREA)

Abstract

本发明公开了一种基于双偏差双空间局部方向模式的人脸识别方法。本发明步骤:1、根据图像的局部纹理特性,通过计算局部图像的绝对偏差和相对偏差得到梯度空间的有效信息,并引入度量函数;2、计算局部图像各方向的灰度值之和,获得灰度空间的信息;3、利用双空间信息,量化特征对图像进行模式编码,得到整个人脸图像的特征图。将特征图分为若干不重叠子块,依据信息熵加权级联各子块的直方图得到人脸的特征向量;4、利用最近邻分类器,计算卡方距离判断向量间的相似程度,完成识别。本发明在光照、姿态、表情、遮挡等干扰下,能保持很好的稳定性和鲁棒性。

Description

基于双偏差双空间局部方向模式的人脸识别方法
技术领域
本发明属于生物特征识别及信息安全的技术领域,特别涉及一种基于双偏差双空间局部方向模式的人脸识别方法。
背景技术
近年来人脸识别因为其稳定、非接触性和易于获取等特点,已经在实际生活中经常应用,例如地铁安检、银行身份验证、安防监控等等。但是实际应用场景环境复杂多变,光照、表情、遮挡、姿态等条件的变化都会显著影响人脸识别的性能。
人脸识别是利用计算机视觉在图像或视频中找到人脸并识别出其真实身份的一种技术。人脸识别主要有以下几个步骤:人脸检测,人脸表征和人脸匹配等等。
基于局部模式的人脸识别算法聚焦于图像的微纹理结构,对原始图像以新的模式编码获得新的特征图像,新的特征图像对于光照、表情等干扰依然保持很好的稳定性。这种特征在提取人脸局部信息的细节部分的同时,通过级联各个图像块的直方图特征,进一步获得整个人脸图像的信息,因此其对光照、表情、遮挡等环境变化相比全局特征都有更好的鲁棒性。
发明内容
本发明的目的是在光照、表情、遮挡等干扰情况下,提出一种基于双偏差双空间局部方向模式的人脸识别方法,以提高人脸识别性能。本方法在考虑梯度空间中包含稳定的边缘信息的同时兼顾了灰度空间的特征信息,通过综合利用多空间的特征信息,进行模式编码后,得到对环境变化有更强抗干扰能力的特征图。然后对特征图分块提取直方图,利用信息熵加权级联得到整个特征图像的直方图特征向量。最后利用卡方距离通过最近邻分类器进行分类,并将其应用于人脸识别领域。
本发明具体包括如下步骤:
步骤1、通过图像采集设备采集人脸数据信息,采集多个人的人脸数据信息作为训练集;
步骤2、通过记录最大灰度值之和对应方向获得灰度空间信息后,然后综合利用梯度和灰度空间下的信息,量化特征对图像进行双偏差双空间模式编码,得到整个人脸图像的特征图;
步骤3、将得到的特征图分为若干不重叠子块,对每一个子块提取统计直方图,依据子块信息熵加权级联所有子块直方图,得到整个人脸的特征向量;
步骤4、用步骤2和步骤3的方法得到待测试样本的特征向量,利用最近邻分类器,通过卡方距离计算待测试人脸样本特征向量和训练样本间的距离,选择具有和测试样本有最小距离的训练样本类别作为测试样本的类别信息,从而实现人脸识别。
本发明有益效果如下:
本发明对人脸图像进行模式编码,得到特征图,特征图对光照和噪声等环境变量相比原图像有更强的抗干扰能力;模式编码过程中,综合利用了梯度和灰度空间的信息,相比传统局部模式方法对人脸有更强的描述能力;对特征图像分块提取直方图,同时依据信息熵来加权级联各子块直方图,既有区分性的提取了人脸特征,又通过子块级联的方式获得了整个图像的全局信息;利用卡方距离和最近邻分类器来匹配人脸,可以快速精确的识别人脸。
本文发明在光照、姿态、表情和遮挡等各种干扰下,都有很好的鲁棒性,在实际中具有不错的应用价值。
附图说明
图1是本发明所设计的人脸识别流程图。
图2是图像邻域的像素分布示意图。
图3是从人脸图像利用直方图得到人脸特征向量的流程图。
图4是Kirsch算子的八个模板。
具体实施方式
下面结合附图对本发明进一步说明。
如图1和2所示,基于双偏差双空间的局部方向模式人脸识别方法,具体包括如下步骤:
步骤1、根据图像的局部纹理特性,通过计算局部图像的绝对偏差和相对偏差得到梯度空间的有效信息,并引入偏差度量函数,以区分不同偏差对结果的影响;
步骤2、从图像的灰度特征出发,计算局部图像各方向的灰度值之和,获得灰度空间的信息;
步骤3、综合利用梯度和灰度空间下的信息,量化特征对图像进行双偏差双空间模式编码,得到整个人脸图像的特征图,将特征图分为若干不重叠子块,依据信息熵加权级联各子块的直方图得到人脸的特征向量;
步骤4、利用最近邻分类器,计算卡方距离判断向量间的相似程度,通过计算待测试人脸样本特征向量和训练样本特征向量的距离,选择具有和测试样本有最小距离的训练样本类别作为测试样本的类别信息,从而实现人脸识别。
步骤1所述的通过计算局部图像的绝对偏差和相对偏差得到梯度空间的有效信息,具体如下:
如图2和4所示,采用Kirsch算子,计算图像的3×3邻域I和8个模板Mi作卷积后,*是卷积符号,得到八个边缘响应值ei,得到局部图像的相对偏差:
ei=I*Mi i=0,1…7
Figure BDA0002214750550000031
其中m代表的在最大边缘响应值的绝对值对应的方向,也就是通过相对偏差记录的第一个量化编码值。ei代表通过采样后的像素与第i个Kirsch模板卷积后得到的边缘响应值,i=0,1,2,3,4,5,6,7分别代表局部图像以gc为中心逆时针的0°、45°、90°、135°、180°、225°、270°和315°的方向;
通过利用局部图像的前向和后向查分获得各方向的绝对偏差值,进一步利用偏差度量函数区分不同偏差:
Figure BDA0002214750550000041
Figure BDA0002214750550000042
Figure BDA0002214750550000043
其中,λ=0.5,dgi代表中心对称各方向上计算得到的绝对差值,i=0,1,2,3分别代表局部图像以gc为中心逆时针的0°、45°、90°和135°的方向。n代表绝对差值最大值对应的方向,也就是通过绝对偏差得到的量化编码值。gi代表局部图像的3×3邻域I中的周围像素值,
Figure BDA0002214750550000049
就是中心像素值,
Figure BDA0002214750550000044
代表偏差度量函数。
如图3所示,步骤2计算局部图像各方向的灰度值之和,获得灰度空间的信息,具体如下:
利用求和公式来求各方向的灰度值之和:
Figure BDA0002214750550000045
式中,Sgi是各个方向的灰度叠加和,i=0,1,2,3分别代表局部图像以gc为中心逆时针的0°、45°、90°和135°的方向。
Figure BDA0002214750550000046
是求和公式。定义如下:
Figure BDA0002214750550000047
然后利用Sgi求最大灰度和对应的方向s进行量化编码:
Figure BDA0002214750550000048
s就是在灰度空间的编码值。
步骤3所述的量化特征对图像进行模式编码,得到整个人脸图像的特征图,将特征图分为若干不重叠子块,依据信息熵加权级联各子块的直方图得到人脸的特征向量,具体如下:
3-1、DVDSLDP=16×m+4×n+s
式中,DVDSLDP是最终的模式编码值。
3-2、通过对原图像进行模式编码后得到整个人脸的特征图。
3-3、将特征图分为若干不重叠子块,依据信息熵加权级联各子块直方图得到整个人脸特征向量。信息熵和加权系数定义如下:
Figure BDA0002214750550000051
Figure BDA0002214750550000052
其中,p(x)是像素值为x的点数占图像所有像素点数的频率,Ri是图像的一个子块区域,H(X)是得到的信息熵大小,ωi是第i个子块的子块系数,Hi(X)是第i个子块的信息熵,N是不重叠子块的个数。
直方图向量定义如下:
Figure BDA0002214750550000053
其中,Hi(c)就是直方图向量,也就是在图像子块中像素值为c的像素点个数。
3-4、依据信息熵加权级联即得到,整个人脸特征向量,定义如下:
Figure BDA0002214750550000054
其中,LH就是级联后的人脸特征向量。
步骤4所述的利用最近邻分类器,计算卡方距离判断向量间的相似程度,通过计算待测试人脸样本特征向量和训练样本特征向量的距离,选择具有和测试样本有最小距离的训练样本类别作为测试样本的类别信息,从而实现人脸识别,具体如下:
4-1、利用卡方距离计算两个向量的距离:
Figure BDA0002214750550000055
其中,χ2(F1,F2)就是两个向量F1和F2的卡方距离,N是向量F1和F2的长度。值越小,代表两向量越相似。
4-2、利用最近邻分类器,求得待测人脸样本的身份,完成人脸识别:
Figure BDA0002214750550000061
其中,y为待测试人脸的特征向量,Fj为训练样本库中人脸身份标签为j的训练样本特征向量。Identity(y)也就是待测人脸的身份标签,即和测试样本有最小距离的训练样本类别标签,得到待测试人脸的身份类别标签,也就完成了人脸识别。

Claims (3)

1.基于双偏差双空间局部方向模式的人脸识别方法,其特征在于包括如下步骤:
步骤1、根据图像的局部纹理特性,通过计算局部图像的绝对偏差和相对偏差得到梯度空间的有效信息,并引入偏差度量函数,以区分不同偏差对结果的影响;
步骤2、通过记录最大灰度值之和对应方向获得灰度空间信息后,然后综合利用梯度和灰度空间下的信息,量化特征对图像进行双偏差双空间模式编码,得到整个人脸图像的特征图;
步骤3、将得到的特征图分为若干不重叠子块,对每一个子块提取统计直方图,依据子块信息熵加权级联所有子块直方图,得到整个人脸的特征向量;
步骤4、利用最近邻分类器,通过卡方距离计算待测试人脸样本特征向量和训练样本间的距离,选择具有和测试样本有最小距离的训练样本类别作为测试样本的类别信息,从而实现人脸识别;
步骤1所述的绝对偏差、相对偏差以及度量函数,具体实现如下:
采用Kirsch算子,计算图像的3×3邻域I和8个模板Mi作卷积后,*是卷积符号,得到八个边缘响应值ei,得到局部图像的相对偏差:
ei=I*Mi i=0,1…7
Figure FDA0003768896150000011
其中m代表的在最大边缘响应值绝对值对应的方向,也就是通过相对偏差记录的第一个量化编码值;ei代表通过采样后的像素与第i个Kirsch模板卷积后得到的边缘响应值,i=0,1,2,3,4,5,6,7分别代表局部图像以gc为中心逆时针的0°、45°、90°、135°、180°、225°、270°和315°的方向;
通过利用局部图像的前向和后向查分获得各方向的绝对偏差值,进一步利用度量函数区分不同偏差:
Figure FDA0003768896150000012
Figure FDA0003768896150000013
Figure FDA0003768896150000021
其中,λ=0.5,dgi代表中心对称各方向上计算得到的绝对差值,i=0,1,2,3分别代表局部图像以gc为中心逆时针的0°、45°、90°和135°的方向;n代表绝对差值最大值对应的方向,也就是通过绝对偏差得到的量化编码值;gi代表局部图像的3×3邻域I中的周围像素值,
Figure FDA0003768896150000022
就是中心像素值,
Figure FDA0003768896150000023
代表度量函数,衡量不同大小的偏差;
步骤2所述的通过记录最大灰度值之和对应方向获得灰度空间信息后,综合利用梯度和灰度空间下的信息,量化特征对图像进行双偏差双空间模式编码,得到整个人脸图像的特征图;具体实现如下:
2-1、记录最大灰度值之和对应方向获得灰度空间的信息,首先求各方向的灰度值之和:
Figure FDA0003768896150000024
式中,Sgi是各个方向的灰度叠加和,i=0,1,2,3分别代表局部图像以gc为中心逆时针的0°、45°、90°和135°的方向;
Figure FDA0003768896150000025
是求和公式;定义如下:
Figure FDA0003768896150000026
然后利用Sgi求最大灰度值之和对应的方向s对灰度空间特征量化编码:
Figure FDA0003768896150000027
s就是在灰度空间的编码值;
2-2、利用双偏差双空间特征模式进行模式编码:
DVDSLDP=16×m+4×n+s
式中,DVDSLDP就是量化双空间特征后,得到的模式编码值;
2-3、通过步骤2-2对图像中的每一点进行模式编码后,就能够得到整个人脸图像的特征图。
2.根据权利要求1所述的基于双偏差双空间的局部方向模式人脸识别方法,其特征在于步骤3具体实现如下:
3-1、将特征图分为若干不重叠子块,依据信息熵加权级联各子块直方图得到整个人脸特征向量;信息熵和加权系数定义如下:
Figure FDA0003768896150000031
Figure FDA0003768896150000032
其中,p(x)是像素值为x的点的个数占图像所有像素点总数的频率,Ri是图像的一个子块区域,H(X)是得到的信息熵大小,ωi是第i个子块的子块系数,Hi(X)是第i个子块的信息熵,N是不重叠子块的个数;直方图向量定义如下:
Figure FDA0003768896150000033
其中,Hi(c)就是直方图向量,也就是在图像子块中像素值为c的像素点的个数;
3-2、依据信息熵对直方图向量加权级联,即可得到整个人脸的特征向量,定义如下:
Figure FDA0003768896150000034
其中,LH就是级联后的人脸特征向量。
3.根据权利要求2所述的基于双偏差双空间的局部方向模式人脸识别方法,其特征在于步骤4具体实现如下:
4-1、利用卡方距离计算两个向量的距离:
Figure FDA0003768896150000035
其中,χ2(F1,F2)就是两个向量F1和F2的卡方距离,N是向量F1和F2的长度;值越小,代表两向量越相似;
4-2、利用最近邻分类器,求得待测人脸样本的身份,完成人脸识别:
Figure FDA0003768896150000036
其中,y为待测试人脸的特征向量,Fj为训练样本库中人脸身份标签为j的训练样本特征向量;Identity(y)也就是待测人脸的身份标签,即和测试样本有最小距离的训练样本类别标签,得到待测试人脸的身份类别标签,也就完成了人脸识别。
CN201910911208.1A 2019-09-25 2019-09-25 基于双偏差双空间局部方向模式的人脸识别方法 Active CN110781749B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910911208.1A CN110781749B (zh) 2019-09-25 2019-09-25 基于双偏差双空间局部方向模式的人脸识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910911208.1A CN110781749B (zh) 2019-09-25 2019-09-25 基于双偏差双空间局部方向模式的人脸识别方法

Publications (2)

Publication Number Publication Date
CN110781749A CN110781749A (zh) 2020-02-11
CN110781749B true CN110781749B (zh) 2022-09-02

Family

ID=69384548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910911208.1A Active CN110781749B (zh) 2019-09-25 2019-09-25 基于双偏差双空间局部方向模式的人脸识别方法

Country Status (1)

Country Link
CN (1) CN110781749B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102667815A (zh) * 2009-10-02 2012-09-12 高通股份有限公司 用于耐遮挡面部识别的方法和系统
CN110010228A (zh) * 2019-03-26 2019-07-12 广州艾颜佳美容美发设备有限公司 一种基于图像分析的脸部皮肤透视算法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102667815A (zh) * 2009-10-02 2012-09-12 高通股份有限公司 用于耐遮挡面部识别的方法和系统
CN110010228A (zh) * 2019-03-26 2019-07-12 广州艾颜佳美容美发设备有限公司 一种基于图像分析的脸部皮肤透视算法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
双空间局部方向模式的人脸识别;杨恢先 等;《中国图象图形学报》;20171130;全文 *

Also Published As

Publication number Publication date
CN110781749A (zh) 2020-02-11

Similar Documents

Publication Publication Date Title
CN109118479B (zh) 基于胶囊网络的绝缘子缺陷识别定位装置及方法
Sirmacek et al. Urban area detection using local feature points and spatial voting
CN101980250B (zh) 基于降维局部特征描述子和隐条件随机场的目标识别方法
CN108921120B (zh) 一种适应广泛零售场景下的香烟识别方法
CN106960176B (zh) 一种基于超限学习机和颜色特征融合的行人性别识别方法
CN109635726B (zh) 一种基于对称式深度网络结合多尺度池化的滑坡识别方法
CN113076802A (zh) 基于缺少断开图像样本的变电站开关分合状态图像识别方法
CN116137036B (zh) 基于机器学习的基因检测数据智能处理系统
CN110766708A (zh) 基于轮廓相似度的图像比较方法
CN105447488B (zh) 基于素描线段拓扑结构的sar图像目标检测方法
CN108256557B (zh) 结合深度学习和邻域集成的高光谱图像分类方法
CN117475236B (zh) 用于矿产资源勘探的数据处理系统及其方法
CN114612315A (zh) 一种基于多任务学习的高分辨率影像缺失区域重建方法
CN110781749B (zh) 基于双偏差双空间局部方向模式的人脸识别方法
Ma et al. Tree-shaped sampling based hybrid multi-scale feature extraction for texture classification
CN114119577B (zh) 高铁隧道漏缆卡扣检测方法
Wan et al. Statistical multiscale image segmentation via alpha-stable modeling
CN115965796A (zh) 一种基于图像识别的金属锈蚀检测方法及系统
CN111259914B (zh) 一种茶树叶特征信息高光谱提取方法
CN114202536A (zh) 一种基于AT-EnFCM-SCNN的灾害环境特征提取方法
Jyothy et al. Texture-based multiresolution steganalytic features for spatial image steganography
CN107292854B (zh) 基于局部奇异性量化分析的灰度图像增强方法
CN112801950A (zh) 一种基于几何失真度量的图像适配质量评价方法
CN116129280B (zh) 一种遥感影像雪检测的方法
CN117173545B (zh) 一种基于计算机图形学的证照原件识别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant