CN102737237B - 基于局部关联保持的人脸图像降维方法 - Google Patents
基于局部关联保持的人脸图像降维方法 Download PDFInfo
- Publication number
- CN102737237B CN102737237B CN201210248646.2A CN201210248646A CN102737237B CN 102737237 B CN102737237 B CN 102737237B CN 201210248646 A CN201210248646 A CN 201210248646A CN 102737237 B CN102737237 B CN 102737237B
- Authority
- CN
- China
- Prior art keywords
- vector
- matrix
- facial image
- row
- calculate
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于局部关联保持的人脸图像降维方法,首先用多维向量表示人脸图像,依据两向量差向量的范数得到每个向量的k个近邻,并按照径向基函数计算每个向量的k个近邻的归一化的权重。计算每个向量与其k个近邻加权之和的差向量,通过每个差向量的转置与其本身相乘得到矩阵,将所有向量对应的矩阵相加,得到局部关联保持矩阵。通过计算局部关联保持矩阵的特征值及特征向量,并选择部分大的特征值对应的特征向量作为基向量组成投影矩阵,实现降维。降维后的人脸图像很好地保持了数据局部关联,有利于图像识别,通过本发明方法提取特征后的分类效果优于PCA及LPP;降低了计算复杂性,揭示了新方法与PCA及LPP间的关系。
Description
技术领域
本发明涉及一种人脸图像降维方法,尤其涉及一种基于局部关联保持的人脸图像降维方法。
背景技术
人脸图像由大量像素点值组成,通过高维向量或高阶矩阵表示,人脸图像识别需要大量的计算与存储代价,导致维数灾难,因此在对人脸图像操作前,需要对人脸图像进行降维处理,即将原始人脸图像映射到一个低维空间,得到低维空间表示人脸图像的最主要特征,降低计算与存储代价,实现人脸图像的自动识别。
目前经典的不考虑数据类别标记的降维方法为PCA(Primal Component Analysis:主成分分析),在计算及理论分析方面简单,特征脸即为PCA在人脸图像识别中最著名的应用。PCA降维只考虑数据的全局分布特征,忽略了数据的非线性结构特征。KPCA扩展了PCA算法,通过借助核变换将数据映射到更高维空间,并在新空间对数据进行线性特征提取,实现基于核的非线性特征提取。方法的优点是不需要明确非线性映射函数,特征提取只需计算高维空间向量内积核函数,提取的特征更加有效地描述了人脸图像的非线性结构。人脸图像数据具有流形特征,应用流形学习技术实现人脸图像降维,并保持源图像的流形特征。主要算法有LLE(local linear Embedding:局部线性嵌入)、NPE(Neighborhood Preserving Embedding:局部保持嵌入)、LPP(Locality Preserving Projection:局部保持投影)。基于流形学习的特征提取本质上在保持数据间局部关系的基础上,学习数据内在的非线性特征,更加有利于人脸图像的识别。上述流形学习方法的缺点是计算复杂,没有给出明确的非线性变换矩阵,不能直接得到未标记人脸图像的变换特征。
发明内容
本发明的目的就是为了解决上述问题,提供一种基于局部关联保持的人脸图像降维方法,它具有在保持图像数据局部关联的前提下,通过简单计算实现人脸特征降维优点。
为了实现上述目的,本发明采用如下技术方案:
一种基于人脸图像局部关联保持的特征降维方法,首先用多维向量表示人脸图像,依据两向量差向量的范数得到每个向量的k个近邻,并按照径向基函数计算每个向量的k个近邻的归一化的权重。计算每个向量与其k个近邻加权之和的差向量,通过每个差向量的转置与其本身相乘得到矩阵,将所有向量对应的矩阵相加,得到局部关联保持矩阵。通过计算局部关联保持矩阵的特征值及特征向量,并选择部分大的特征值对应的特征向量作为基向量组成投影矩阵,实现降维。人脸图像通过投影矩阵影射到低维空间,在低维空间实现人脸图像识别。
本发明的具体步骤为:
步骤一:将m幅大小为s×t个像素点的人脸图像表示为s×t维的行向量x1,x2,…,xi,…,xm,其中m为人脸图像个数,s为图像行像素点数,t为列图像列像素点数,xi表示第i幅人脸图像对应的s×t维的行向量,这m幅人脸图像包含p个人,每人幅图像;
步骤二:对于任意一个行向量xi(i∈{1,2,…,m}),计算dij=‖xi-xj‖(j∈{1,2,…,m}且j≠i),从中选择k个(k=9)使得dij最小的行向量,组成集合记为Ne(xi)。其中‖·‖表示向量的范数,d表示xi-xj所得差向量的范数;
步骤三:计算权重矩阵W,矩阵第i行第j列成员记为wij;
步骤四:计算每个行向量对应的差向量,并计算局部关联保持矩阵V;
步骤五:求解矩阵V的特征值及特征向量,选择d个最大的特征值对应的特征向量,并由该d个特征向量作为列组成一个矩阵,称为投影矩阵M;d=min(m,n),其中n=r(V)为矩阵V的秩;
步骤六:将p个人的任意一个人脸图像行向量通过步骤5得到的投影矩阵M映射到低维空间。
所述步骤三中权重矩阵W的第i行第j列成员wij定义如下:
其中参数σ=2。
所述步骤四中差向量计算方法为:对于任意一个行向量xi,计算其与k(k=9)个近邻行向量的加权之和的差向量ri
所述步骤四中局部关联保持矩阵V的计算如下:
所述步骤六中,对于任意一个行向量xi,计算xiM得到一个d维行向量,其中d=min(m,n)□s×t,从而将人脸图像降到d维,实现人脸图像降维。
本发明的有益效果:降维后的人脸图像很好地保持了数据局部关联,有利于图像识别,Yale人脸数据库上的分类误差率比较,通过本发明方法提取特征后的分类效果优于PCA及LPP;本发明降低了计算复杂性,同时揭示了新方法与PCA及LPP间的关系,当k=m,并且参数σ足够大时,新方法退化为PCA;分析表明,新方法特征=LPP特征+复杂非线性隐性特征。
附图说明
图1为本发明的降维方法流程图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
如图1所示为本发明的降维方法流程图,一种基于局部关联保持的人脸图像降维方法,
具体步骤为:
步骤一:将m幅大小为s×t个像素点的人脸图像表示为s×t维的行向量x1,x2,…,xi,…,xm,其中m为人脸图像个数,s为图像行像素点数,t为列图像列像素点数,xi表示第i幅人脸图像对应的s×t维的行向量,这m幅人脸图像包含p个人,每人幅图像;
步骤二:对于任意一个行向量xi(i∈{1,2,…,m}),计算dij=‖xi-xj‖(j∈{1,2,…,m}且j≠i),从中选择k个(k=9)使得dij最小的行向量,组成集合记为Ne(xi)。其中‖·‖表示向量的范数,dij表示xi-xj所得差向量的范数;
步骤三:计算权重矩阵W,矩阵第i行第j列成员记为wij;
步骤四:计算每个行向量对应的差向量,并计算局部关联保持矩阵V;
步骤五:求解矩阵V的特征值及特征向量,选择d个最大的特征值对应的特征向量,并由该d个特征向量作为列组成一个矩阵,称为投影矩阵M;d=min(m,n),其中n=r(V)为矩阵V的秩;
步骤六:将p个人的任意一个人脸图像行向量通过步骤5得到的投影矩阵M映射到低维空间。
所述步骤三中权重矩阵W的第i行第j列成员wij定义如下:
其中参数σ=2。
所述步骤四中差向量计算方法为:对于任意一个行向量xi,计算其与k(k=9)个近邻行向量的加权之和的差向量ri
所述步骤四中局部关联保持矩阵V的计算如下:
本发明以Yale人脸数据库作为实施实例,该数据库包含15个人的总计165幅人脸图像,每人11幅图像分别在不同光照条件及不同表情下获得。将人脸图像变换为32x32像素的图像,灰度级255。降维前对数据进行归一化处理。具体实验步骤如下:
(1)选择每个人的前4幅图像,将每幅图像变换为1024维的行向量,总计60个向量;
(2)计算向量两两间的差向量的范数,对每个向量选择与其距离最近的9个向量组成其近邻向量集合;
(3)按照公式(1)计算权重矩阵,该矩阵为60X60的矩阵,并按照公式(2)计算每个向量对应的差向量,该差向量为1024维的行向量;
(4)按照公式(3)计算局部关联保持矩阵V,该矩阵为1024X1024维矩阵;
(5)计算V的特征值及特征向量,该矩阵的秩为n=1024,大于60,所以选择降维后的维数为60;取V的最小的60个特征值对应的特征向量,按列组成1024X60维的映射矩阵M;
(6)在Yale人脸库中任选一副人脸图像,转换为1024维行向量,按照发明步骤6得到对应的变换后的60维特征向量,实现降维。
在Yale人脸数据库上的分类误差率比较(采用线性带惩罚项SVM支持向量机分类算法,SVM中的参数C=1);数据集Yale,165幅图降维后采用五折交叉验证分类结果如表1所示:
降维方法 | 分类精度及方差 |
本发明方法 | 65.7303±1.0744 |
PCA降维 | 63.6364±0.5051 |
LPP降维 | 60±1.5335 |
表1 在Yale人脸数据库上的分类误差率比较
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (3)
1.一种基于局部关联保持的人脸图像降维方法,其特征是,首先用多维向量表示人脸图像,依据两向量差向量的范数得到每个向量的k个近邻,并按照径向基函数计算每个向量的k个近邻的归一化的权重;计算每个向量与其k个近邻加权之和的差向量,通过每个差向量的转置与其本身相乘得到矩阵,将所有向量对应的矩阵相加,得到局部关联保持矩阵;通过计算局部关联保持矩阵的特征值及特征向量,并选择部分大的特征值对应的特征向量作为基向量组成投影矩阵,实现降维;人脸图像通过投影矩阵影射到低维空间,在低维空间实现人脸图像识别。
2.一种基于局部关联保持的人脸图像降维方法,其特征是,具体步骤为:
步骤一:将m幅大小为s×t个像素点的人脸图像表示为s×t维的行向量x1,x2,…,xi,…,xm;其中m为人脸图像个数,s为图像行像素点数,t为图像列像素点数,xi表示第i幅人脸图像对应的s×t维的行向量,这m幅人脸图像包含p个人,每人幅图像;
步骤二:对于任意一个行向量xi,i∈{1,2,…,m},计算dij=||xi-xj||j∈{1,2,…,m}但j≠i,从中选择k,k=9个使得dij最小的行向量,组成集合记为Ne(xi);其中||·||表示向量的范数,dij表示xi-xj所得差向量的范数;
步骤三:计算权重矩阵W,矩阵第i行第j列成员记为wij;
步骤四:计算每个行向量对应的差向量,并计算局部关联保持矩阵V;
步骤五:求解矩阵V的特征值及特征向量,选择d个最大的特征值对应的特征向量,并由该d个特征向量作为列组成一个矩阵,称为投影矩阵M;
d=min(m,n),其中n=r(V)为矩阵V的秩;
步骤六:将P个人的任意一个人脸图像行向量通过步骤五得到的投影矩阵M映射到低维空间;
所述步骤三中权重矩阵W的第i行第j列成员wij定义如下:
所述步骤四中差向量计算方法为:对于任意一个行向量xi,计算其与k个近邻行向量的加权之和的差向量ri,k=9;
所述局部关联保持矩阵V的计算如下:
3.如权利要求2所述一种基于局部关联保持的人脸图像降维方法,其特征是,所述步骤六中,对于任意一个行向量xi,计算xiM得到一个d维行向量,其中d=min(m,n)<<s×t,从而将人脸图像降到d维,实现人脸图像降维。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210248646.2A CN102737237B (zh) | 2012-07-18 | 2012-07-18 | 基于局部关联保持的人脸图像降维方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210248646.2A CN102737237B (zh) | 2012-07-18 | 2012-07-18 | 基于局部关联保持的人脸图像降维方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102737237A CN102737237A (zh) | 2012-10-17 |
CN102737237B true CN102737237B (zh) | 2014-06-18 |
Family
ID=46992706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210248646.2A Expired - Fee Related CN102737237B (zh) | 2012-07-18 | 2012-07-18 | 基于局部关联保持的人脸图像降维方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102737237B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514443B (zh) * | 2013-10-15 | 2016-12-07 | 中国矿业大学 | 一种基于lpp特征提取的单样本人脸识别迁移学习方法 |
CN103679709B (zh) * | 2013-11-28 | 2017-08-01 | 深圳职业技术学院 | 光谱图像处理方法 |
CN104063717B (zh) * | 2014-07-09 | 2017-03-22 | 哈尔滨工业大学 | 局部结构保持的全局监督图嵌入复杂***监测数据可视化方法 |
CN104881682B (zh) * | 2015-05-26 | 2019-03-05 | 东南大学 | 一种基于保局映射与主成分分析的图像分类方法 |
CN106407982B (zh) * | 2016-09-23 | 2019-05-14 | 厦门中控智慧信息技术有限公司 | 一种数据处理方法以及设备 |
CN106338977A (zh) * | 2016-10-17 | 2017-01-18 | 沈阳化工大学 | 基于差分局部保持投影的非线性过程故障检测方法 |
CN107704887B (zh) * | 2017-10-20 | 2021-04-02 | 北京工业大学 | 一种基于f范数的局部保持投影的图像识别方法 |
CN108121965B (zh) * | 2017-12-21 | 2020-04-21 | 深圳大学 | 基于鲁棒联合稀疏特征提取的图像识别方法 |
CN109885578B (zh) * | 2019-03-12 | 2021-08-13 | 西北工业大学 | 数据处理方法、装置、设备及存储介质 |
CN111553417B (zh) * | 2020-04-28 | 2022-04-22 | 厦门大学 | 基于判别正则化局部保留投影的图像数据降维方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079103A (zh) * | 2007-06-14 | 2007-11-28 | 上海交通大学 | 基于稀疏贝叶斯回归的人脸姿势识别方法 |
CN101187986A (zh) * | 2007-11-27 | 2008-05-28 | 海信集团有限公司 | 基于有监督的近邻保持嵌入和支持向量机的人脸识别方法 |
-
2012
- 2012-07-18 CN CN201210248646.2A patent/CN102737237B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101079103A (zh) * | 2007-06-14 | 2007-11-28 | 上海交通大学 | 基于稀疏贝叶斯回归的人脸姿势识别方法 |
CN101187986A (zh) * | 2007-11-27 | 2008-05-28 | 海信集团有限公司 | 基于有监督的近邻保持嵌入和支持向量机的人脸识别方法 |
Non-Patent Citations (4)
Title |
---|
Feature dimension reduction for microarray data analysis using locally linear embedding;SHI CHAO et al.;《The Asia Pacific Bioinformatics Conference》;20040916;211-217 * |
SHI CHAO et al..Feature dimension reduction for microarray data analysis using locally linear embedding.《The Asia Pacific Bioinformatics Conference》.2004, |
刘晓宁等.流形学习在三维人脸特征降维中的应用.《计算机应用研究》.2010,第27卷(第10期),3982-3984. |
流形学习在三维人脸特征降维中的应用;刘晓宁等;《计算机应用研究》;20101031;第27卷(第10期);3982-3984 * |
Also Published As
Publication number | Publication date |
---|---|
CN102737237A (zh) | 2012-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102737237B (zh) | 基于局部关联保持的人脸图像降维方法 | |
Qiao | New SVD based initialization strategy for non-negative matrix factorization | |
CN111192292B (zh) | 基于注意力机制与孪生网络的目标跟踪方法及相关设备 | |
US11429903B2 (en) | Privacy-preserving asynchronous federated learning for vertical partitioned data | |
CN110163239B (zh) | 一种基于超像素和条件随机场的弱监督图像语义分割方法 | |
Yun et al. | PNN based crop disease recognition with leaf image features and meteorological data | |
CN106570183B (zh) | 一种彩色图像检索和分类方法 | |
CN105095902B (zh) | 图片特征提取方法及装置 | |
CN104268593A (zh) | 一种小样本情况下多稀疏表示的人脸识别方法 | |
CN108921061A (zh) | 一种表情识别方法、装置和设备 | |
CN103456013B (zh) | 一种表示超像素以及度量超像素之间相似性的方法 | |
Lu | Hybridizing nonlinear independent component analysis and support vector regression with particle swarm optimization for stock index forecasting | |
CN108428055B (zh) | 一种考虑负荷纵向特性的负荷聚类方法 | |
CN106599227B (zh) | 用于获取基于属性值的对象之间的相似度的方法与装置 | |
CN102592129B (zh) | 基于场景驱动的智能手机图像特征点选取方法 | |
Ye et al. | Fast orthogonal linear discriminant analysis with application to image classification | |
Peluffo-Ordónez et al. | Short review of dimensionality reduction methods based on stochastic neighbour embedding | |
CN108629405B (zh) | 提高卷积神经网络计算效率的方法和装置 | |
CN111242183A (zh) | 基于注意力机制的图像识别分类方法及装置 | |
CN108491883B (zh) | 一种基于条件随机场的显著性检测优化方法 | |
CN105023239A (zh) | 基于超像素和最大边界分布的高光谱数据降维方法 | |
Rathie et al. | Stable and generalized-t distributions and applications | |
Zhou et al. | Safe subspace screening for nuclear norm regularized least squares problems | |
CN101504723B (zh) | 一种建立投影空间的方法和装置 | |
Ezoji et al. | Use of matrix polar decomposition for illumination-tolerant face recognition in discrete cosine transform domain |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140618 Termination date: 20200718 |