一种人脸图像快速特征提取方法

摘要

一种人脸图像快速特征提取方法，第一步，初始化读入训练图像集，将获得的彩色人脸图像转换到灰度空间；第二步，对图像进行线性插值和直方图均衡化；第三步，构造训练集数据矩阵；第四步，训练集应用主成份分析法，求取其低维表述矩阵；第五步，求原训练图像的低维表述矩阵的隶属度矩阵；第六步，计算低维表述矩阵的类内散度矩阵和类间散度矩阵；第七步，求得训练人脸图像集的特征表述矩阵。本发明充分考虑了训练集中各个人的人脸图像之间的关系，利用了其对于某一人脸子集的属于和不属于的信息，得到较好的类中心估计，弱化了某一图像的剧烈扰动对于该人其他图像的扰动，从而更加有效地提取人脸特征。

主权项

1.一种人脸图像快速特征提取方法，其特征在于按照以下步骤进行：步骤1：初始化，采集包含R个人的彩色人脸图像组成含有M个人脸图像的集合，R为人数且为正整数，M为正整数，第r个人的人脸图像数为t_r(r＝1，2，...，R，t_r为自然数且t_r≥4)，r为正整数，每幅人脸图像的数据大小为W_c×H_c×3，其中，c为正整数且c＝1，2，...，M，W_c，H_c为正整数，分别表示第c幅人脸图像的列数和行数，然后将获得的彩色图像转化为灰度图像，得到灰度图像的大小为W_c×H_c；步骤2：对步骤1中处理后的人脸图像进行图像预处理，方法如下：步骤2.1：用公知的线性插值法分别对第g个图像进行尺度调整，g为正整数且g＝1,2,...,M，将第g个图像规格化为W×H的图像，W,H为正整数，分别表示尺度归一化之后的人脸图像统一的列数和行数；步骤2.2：对步骤2.1中处理后的人脸图像进行直方图均衡化；步骤3：构造训练集矩阵：对步骤2中得到的M幅人脸图像的数据格式进行调整，令原来的一幅图像的像素数据矩阵为A，则有A＝[A₁，A₂，…，A_W]                   （1）其中W表示该图像的列数，A_w表示像素数据矩阵A的第w列，w为正整数且w＝1,2,...,W，将像素数据矩阵A中数据按照如下方法进行调整，得到新的W×H行1列的像素数据矩阵A'：$A^{\prime }={[A_1^T,A_W^T,...,A_W^T]}^T---\left(2\right)$将原来的M幅图像均按照上述方法，调整为W×H行1列，再将每一幅图像转换后的数据作为另一矩阵的一列，构造训练集矩阵，即Γ＝[Γ₁,Γ₂,...,Γ_M]                 （3）其中，Γ即为由M个人脸图像构成的训练集矩阵，Γ_l是第l个图像数据经过式（1）和式（2）变换得到的数据，l为正整数，且l＝1,2,...,M；步骤4：对训练集矩阵应用主成份分析法，求取训练集矩阵的低维表述矩阵，其方法如下：步骤4.1：求取训练集矩阵Γ的均值向量，对训练集矩阵Γ每一行对应位置的数据进行数值平均，得到训练集矩阵Γ的均值向量，即$\bar{\Gamma }=\frac{1}{M}\underset{u=1}{\overset{M}{\Sigma }}{\Gamma }_u---\left(4\right)$其中，u为正整数，且u＝1,2,...,M，为训练集矩阵Γ的均值向量；步骤4.2：将训练集矩阵Γ根据其均值向量进行中心化，并求取中心化之后的协方差矩阵的转置矩阵，即$C=\underset{y=1}{\overset{M}{\Sigma }}{({\Gamma }_y-\bar{\Gamma })}^T({\Gamma }_y-\bar{\Gamma })---\left(5\right)$其中，C即为中心化后的训练集矩阵的协方差矩阵的转置矩阵，y为正整数，且y＝1,2,...，M，Γ_y为原训练集矩阵的第y列；步骤4.3：对步骤4.2中求得的矩阵C求取特征值和特征向量，并将特征值按照数值大小降序排列，按照排列后的特征值的顺序重新排列特征值对应特征向量，使矩阵C的特征向量排序后构成新的特征向量矩阵，即：Q＝[Q₁,Q₂,...,Q_M]                          （6）U＝[μ₁,μ₂,...,μ_M],μ₁＞μ₂＞...＞μ_M    （7）其中，Q为矩阵C的特征向量排序后构成的特征向量矩阵，U为矩阵C的特征值组成的行向量，μ_ρ为矩阵C的特征值，Q_ρ为μ_ρ对应的特征向量，ρ为正整数，且ρ＝1，2，...，M；步骤4.4：求取简化特征向量矩阵，对步骤4.3中求的特征向量矩阵Q选取前L列，构成矩阵C的简化特征向量矩阵，L为正整数，且L＝M-R，即S＝[Q₁,Q₂,...,Q_L]                          （8）其中S为矩阵C的简化特征向量矩阵；步骤4.5：将中心化之后的训练集矩阵在简化特征向量矩阵S上投影，得到特征脸矩阵F，特征脸矩阵F的每一列按照下式求取：$F_x=({\Gamma }_x-\bar{\Gamma })S---\left(9\right)$其中，F_x是特征脸矩阵F的第x列，x为正整数，且有x＝1,2,...,L；步骤4.6：将训练集矩阵Γ在特征脸矩阵F上投影，得到训练集矩阵Γ的基于特征脸矩阵F的低维表述矩阵X，即X＝F^TΓ                            （10）其中，X为低维表述矩阵；步骤5：求取步骤4中所得的低维表述矩阵X的隶属度矩阵，方法如下：步骤5.1：设定一个M行M列的矩阵G，求取低维表述矩阵X中第a列与第b列的欧几里德距离，作为矩阵G的第a行第b列的元素G_ab，a,b是正整数且有a,b∈[1,M]，特别地，当a＝b时，设置G_ab为无穷大；步骤5.2：创建矩阵G的副本G'，将矩阵G′每一列中的元素按照数值大小升序排列，再将矩阵G'的元素替换为矩阵G的同一列中与该元素数值相等的元素的在矩阵G中的行标号，特别地，若G'的某列中存在数值相等的元素，则用矩阵G的对应列中与该元素数值相等的元素的行标号按照数值大小以升序的方式替换G′中出现的相等元素；步骤5.3：提取训练集矩阵Γ中每列的属于信息，选取G′的前k行元素构成新的矩阵B，k为正整数且k＝3，则矩阵B的第e行第f列元素B_ef的值Δ_ef表示在低维表述矩阵X中，与X的第f列欧几里德距离最小的k个列按照与第f列欧几里德距离升序排序，则X的第Δ_ef列排在第e位，其中，e,f为正整数，且有e＝1,2,3和f＝1,2,...,M，B_ef为矩阵B的第e行第f列元素，Δef为B_ef的数值且为正整数；步骤5.4：提取训练集矩阵Γ中每列的不属于信息，选取G′的除最后一行之外的后q行元素构成新的矩阵D，q为正整数且q＝10，则矩阵D的第d行第h列元素D_dh的值θ_dh表示在低维表述矩阵X中，与X的第h列欧几里德距离最大的q个列按照与第h列欧几里德距离升序排序，则X的第θ_dh列排在第d位，其中，d,h为正整数，且有d＝1,2,...,10和h＝1,2,...,M，D_dh为矩阵D的第d行第h列元素，θ_dh为D_dh的数据且为正整数；步骤5.5：创建一个R行M列的矩阵E作为低维表述矩阵X的隶属度矩阵，并按照下式计算隶属度矩阵E的元素值：${\xi }_{\mathrm{ij}}=\left\{\begin{array}{c}0.51+({0.49n}_{\mathrm{ij}}-{0.49o}_{\mathrm{ij}})/k,{\Gamma }_j\in V_i\\ (0.{49n}_{\mathrm{ij}}-{0.49o}_{\mathrm{ij}})/k,{\Gamma }_j\notin V_i\end{array}\right.---\left(11\right)$其中，i,j为正整数，且有i＝1,2,...,R和j＝1,2,...,M，ξ_ij是低维表述矩阵X的隶属度矩阵E第i行j列的元素且表示用于训练图片中第j个人脸图像隶属于第i个人的程度，V_i表示训练集中所有表示第i个人的人脸图像的列的集合，Γ_j∈V_i表示训练集中第j幅人脸图像是第i个人的人脸图像，表示训练集中第j幅人脸图像不是第i个人的人脸图像，n_ij表示矩阵B中第j列各元素数值所对应编号的人脸图像中属于第i个人的个数，o_ij表示矩阵D中第j列各元素数值所对应编号的人脸图像中属于第i个人的个数，并且，当0.49n_ij-0.49o_ij＜0时，令0.49n_ij-0.49o_ij＝0；步骤5.6：将步骤5.5中得到的低维表述矩阵X的隶属度矩阵E按列归一化，得到归一化之后的隶属度矩阵E'，即：${\xi }_{\mathrm{IJ}}^{\prime }=\frac{{\xi }_{\mathrm{IJ}}}{\underset{I=1}{\overset{R}{\Sigma }}{\xi }_{\mathrm{IJ}}}---\left(12\right)$其中，I,J为正整数，且有I＝1,2,...R和J＝1,2,...,M，ξ′_IJ是E′的第I行第J列的元素；步骤6：计算低维表述矩阵X的类内散度矩阵和类间散度矩阵，按如下方法进行：步骤6.1：计算低维表述矩阵X的类中心和平均中心，计算公式如下：${\tilde{X}}_p=\frac{\underset{\eta =1}{\overset{M}{\Sigma }}{\xi }_{\mathrm{p\eta }}^{\prime }{X}_{\eta }}{\underset{\alpha =1}{\overset{M}{\Sigma }}{\xi }_{\mathrm{p\alpha }}^{\prime }}---\left(13\right)$$\bar{X}=\frac{1}{M}\underset{ϵ=1}{\overset{M}{\Sigma }}{X}_{ϵ}---\left(14\right)$其中，p为正整数且p＝1,2,...,R，表示第p个人的人脸图像在低维表述矩阵X中对应的所有列数据的类中心，η,α,ε为正整数，X_η表示低维表述矩阵X的第η列，ξ′_pη表示归一化之后的隶属度矩阵E′第p行第η列；步骤6.2：按照如下公式计算低维表述矩阵X的类内散度矩阵和类间散度矩阵：$S_B=\underset{r=1}{\overset{R}{\Sigma }}{t}_r({\tilde{X}}_r-\bar{X}){({\bar{X}}_r-\bar{X})}^T---\left(16\right)$其中，S_W表示低维表述矩阵X的类内散度矩阵，S_B表示低维表述矩阵X的类间散度矩阵，x,γ为正整数，且有x＝1,2,...,R和γ＝1,2,...,M，X_γ表示低维表述矩阵X的第γ列，表示低维表述矩阵X中第x个人对应的所有列数据的类中心，X_γ∈{V_x}表示X_γ属于X中对应于第x个人的所有列组成的集合，t_r为正整数且表示第r个人的人脸图像数目；步骤7：获得训练人脸图像集的特征表述矩阵，方法如下：步骤7.1：求取的特征向量及其对应的特征值，并将特征值按照数值大小降序排列，按照排列后的特征值的顺序重新排列特征值对应特征向量，使矩阵的特征向量排序后构成新的特征向量矩阵，即：Λ＝[Λ₁,Λ₂,...,Λ_M]                        （17）其中，Λ为矩阵的特征向量排序后构成的特征向量矩阵，Ф为矩阵的特征值组成的行向量，φ_z为矩阵的特征值，Λ_z为φ_z对应的特征向量，z为正整数，且z＝1,2,...，M；步骤7.2：求最优投影矩阵，对步骤7.1中求的特征向量矩阵Λ选取前R-1列，构成低维表述矩阵X的最优投影矩阵Y，即：Y＝[Λ₁,Λ₂,...,Λ_R-1]                       （19）步骤7.3：求特征表述矩阵，得到训练人脸图像集的特征表述，按照如下公式：P＝Y^TX                                       （20）其中，P 即为所提取的训练人脸图像集的特征表述矩阵。