人脸识别中的眼睛定位方法与流程

本发明涉及人脸识别中的眼睛定位方法。
背景技术：
：采用人脸图像来进行身份认证的人脸识别是一种使用非常自然和便利的生物特征识别技术。光照的变化是影响人脸识别系统性能的最主要因素。为了避免光照的影响，人脸识别已广泛采用近红外图像。通常，人脸识别过程包含人脸检测、眼睛定位、预处理、特征提取和比对。鲁棒的眼睛定位在人脸识别系统中起着非常重要的作用。传统的眼睛定位方法可大致分为三类：基于模板的方法、基于外观的方法、基于几何特征的方法。基于几何特征的方法根据眼睛的特性来定位眼睛位置比如边缘和虹膜的强度。这些方法在对比度高的图像上才能获得准确的定位。基于模板的方法，首先基于眼睛形状设计通用眼睛模型，然后用模板匹配方法搜索眼睛。基于外观的方法需要收集大量的训练数据，包含不同个体的不同人脸朝向和不同光照情况下的眼睛。这些训练数据用作训练分类器，并检测训练数据是否通过分类。这种方法能够准确检测眼睛，但是对戴眼镜的人脸的检测效果不稳定。特别地，在人脸识别过程中，反光区域和眼镜边缘是人脸检测和眼睛定位的瓶颈。技术实现要素：本发明提供一种无论图像质量好坏均可以快速、稳定地定位人脸中眼睛的方法及应用该眼睛定位方法的人脸识别系统。本发明的技术方案一提供一种人脸识别中眼睛的定位方法。该眼睛定位方法，首先检测并去除图像中的反光区域，然后进行人脸检测、 检测和去除眼镜的黑框、基于SVM(SupportVectorMachine)的眼睛定位，最后用基于PCA(PrincipalComponentAnalysis)的眼睛定位来定位眼睛位置。在所述检测并去除图像中的反光区域处理中，计算人脸图像的灰度直方图，用Otsu方法确定反射区域的阈值，如果像素满足设定的条件，则被认为在反射区域中。用相邻的不在反射区域的像素平均灰度值替换在反射区域的像素灰度值；在所述人脸检测处理中，在去除反射区域的图像上采用Ada-Boost算法的Viola-Jones人脸检测方法检测人脸；在所述检测和去除眼镜的黑框处理中，对检测到的人脸区域的每个像素，计算归一化的梯度向量，用一个适当的阈值得到一个二值图像，在二值图像上检测眼镜的左边黑框和右边黑框，用相邻的不在黑框区域的像素平均灰度值替换在黑框区域的像素灰度值；在所述的基于SVM的眼睛定位处理中，构造眼睛的训练集，计算以眼睛准确位置为中心的固定大小区域的Gabor系数向量，构造非眼睛的训练集，计算以真实眼睛的像素距离为中心的固定大小区域的Gabor系数向量，对眼睛训练集和非眼睛训练集，训练具有二次核函数的非线性SVM。对于以眼睛周围像素为中心的区域，计算Gabor系数得到SVM评估，以评估的最大值所在的像素作为眼睛位置，并被称为眼睛位置的置信度；在所述的基于PCA的眼睛定位处理中，如果在SVM中眼睛定位的置信度大于预先设定的阈值，则该位置为最终定位结果，不再进行处理，如果在SVM中眼睛定位的置信度不大于预先设定的阈值，采用PCA再次估算眼睛位置，利用手工找到的眼睛位置归一化脸部图像大小，在特定点计算Gabor系数，构造包含这些系数的向量，并作PCA得到线性变换矩阵，对检测到的人脸区域作一系列旋转和缩放变换，对变换后的每个图像，计算Gabor系数，并计算人脸检测的置信度，选取置信度最大的图像，并把该图像中眼睛的平均位置作为原始图像的眼睛位置。本发明的技术方案二提供一种人脸识别系统。该人脸识别系统包括：人脸检测及眼睛定位模块，应用技术方案一的眼睛定位方法，检测并去除反光区域、检测人脸、检测并去除眼镜的黑框、眼睛定位；预处理模块，为避免光照的影响，归一化局部均值和方差；特征提取模块，对检测到的人脸区域提取采样点，对每个采样点计算M个频率和N个方向的Gabor系数绝对值，其中，M、N为大于0的自然数；及比对模块，对于两幅人脸图像，对应的采样点计算M*N维向量的归一化内积，并将内积值相加得到相似度。如果相似度大于预设的阈值，则判断为两幅图像来自同一个人。发明的效果根据本发明的眼睛定位方法及人脸识别系统，无论图像质量好坏均能够快速、稳定地定位眼睛位置。附图说明图1是人脸识别系统的框图。图2是眼睛定位的说明图。具体实施方式下面，结合附图，对本发明的眼睛定位方法及使用该眼睛定位方法的人脸识别系统进行详细说明。眼睛定位方法包括检测并去除反光区域的处理、人脸检测的处理、检测和去除眼镜的黑框的处理、基于SVM的眼睛定位的处理、基于PCA的眼睛定位的处理。<检测并去除反光区域的处理>检测并去除反光区域的处理包括如下(1)～(4)的处理。(1)计算输入图像I的灰度直方图。hist(k)，k＝0，1，…，255为计算出的直方图。(2)用Otsu方法确定反射区域的阈值。Otsu法是公知的方法，这里不再详细叙述。确定出的阈值如下。Th=argmax200sks255(ω1(k)ω1(k)(μ1(k)-μ2(k))2)]]>其中，ω1(k)=Σm=0k-1hist(m)Σm=0255hist(m),ω2(k)=Σm=k255hist(m)Σm=0255hist(m),]]>μ1(k)=Σm=0k-1hist(m)mΣm=0k-1hist(m),μ2(k)=Σm=k255hist(m)Σm=k255hist(m)]]>(3)如果一个像素(i0，j0)满足以下两个条件，则认为该像素在反射区域中，·I(i0，j0)＞Th·rate(i0，j0)＞Th_raterate(i0,j0)=Σ|i-io|≤h|j-jo|≤hI(i,j)/(2h+1)2(Σ|i-io|≤H|j-jo|≤HI(i,j)-Σ|i-io|≤h|j-jo|≤hI(i,j))/((2H+1)2-(2h+1)2)]]>其中h和H是预定义的常量(H＞h)，Th_rate是预定义的阈值。(4)用相邻的不在反射区域的像素的平均灰度值替换在反射区域的像素的灰度值。<人脸检测的处理>采用Ada-Boost算法的Viola-Jones人脸检测方法，根据去除反射区域的图像来检测人脸。<检测和去除眼镜的黑框的处理>检测和去除眼镜的黑框的处理包括如下(1)～(4)的处理。(1)检测对于检测到的人脸区域的每个像素(i，j)，计算归一化的梯度向量grad(i，j)，用一个适当的阈值Th_grad得到一个二值图像B。B(i,j)=1,grad(i,j)≥Th_grad0,otherwise]]>(2)检测眼镜的左边黑框根据图像B计算出如下序列，C(i)=Σj=1w2B(i,j),i=h1,...,h2]]>其中w和h是检测到人脸图像的宽度和高度，h1和h2是预设的常数。在i＝h/2的邻域，找到序列C的局部最大值点，在该点计算B的一个联通区域。然后，用一个抛物线近似该联通区域，如果二次系数大于0则认为该联通区域为眼镜边缘。(3)检测眼镜的右边黑框眼镜右边黑框的检测类似于左边黑框的检测。(4)去除用相邻的不在黑框区域的像素的平均灰度值替换在黑框区域的像素的灰度值<基于SVM的眼睛定位的处理>基于SVM的眼睛定位的处理包括如下(1)～(2)的处理。(1)SVM训练构造眼睛的训练集，计算以眼睛准确位置为中心的固定大小区域的Gabor系数向量。构造非眼睛训练集，计算以真实眼睛的像素距离为中心的固定大小区域的Gabor系数向量。采用上述训练集，训练具有二次核函数的非线性SVM如下：c(x)＝∑iαik(si，x)+b其中si是一个支持向量，αi是权重，b是偏差，k是核函数。(2)眼睛位置估计对于以眼睛周围像素为中心的区域，计算Gabor系数得到SVM评估。以评估的最大值所在的像素作为眼睛位置，这种估算被称为眼睛位置的置信度<基于PCA的眼睛定位的处理>如果在SVM中眼睛定位的置信度大于预先设定的阈值，则该位置为最终定位结果，不再执行基于PCA的眼睛定位。如果置信度不大于预先设定的阈值，采用PCA再次估算眼睛位置。基于PCA的眼睛定位的处理包括(1)～(2)的处理。(1)PCA训练利用手工找到的眼睛位置归一化脸部图像大小，在特定点计算Gabor系数，构造包含这些系数的向量。对向量做PCA。设A为得到的线性变换。(2)估算眼睛位置对检测到的人脸区域作一系列旋转和缩放变换(以质心为中心，左右分别旋转3度、5度、10度、15度等，并分别缩小为0.95、0.90、0.85倍，放大为1.05、1.10、1.15倍等)，对变换后的每个图像，像上面一样计算Gabor系数Gabor_vec。然后按下式计算人脸检测的置信度。score＝||A*Gabor_vec||/||Gabor_vec||选取置信度最大的图像，并把该图像中眼睛的平均位置作为原始图像的眼睛位置。如上所述，能够定位出眼睛的位置。下面，结合图1和图2对使用了本发明的眼睛定位方法的人脸识别系统进行说明。人脸识别系统包括人脸检测及眼睛定位模块、预处理模块、特征提取模块及比对模块。人脸检测及眼睛定位模块执行如下处理。·检测并去除反光区域；·人脸检测；·检测并去除眼镜的黑框；·基于SVM的眼睛定位；·基于PCA的眼睛定位预处理模块执行如下处理。归一化局部均值和方差，归一化为128个7×7大小区域的均值和方差。特征提取模块执行如下处理：在检测到的人脸区域取10×10个采样点，对每个采样点计算5个频率12个方向的Gabor系数绝对值。比对模块执行如下处理。对两幅人脸图像的每个对应采样点计算两个包含60个元素的向量的归一化内积，并将这些值相加得到相似度。如果相似度大于预设的阈值，则认为两张人脸图像来自同一人。本发明人通过实验验证了本发明的效果。在本公司数据库上的测试结果如下：该数据库包含40人的400幅图像，眼睛定位的准确率为97％，等错误率为0.2％。因此，通过本发明的眼睛定位方法，对于现有技术无法定位正确的眼睛的情况，也能够快速、稳定地定位眼睛。如上所述，对本发明的优选的实施方式进行了说明，但上述实施方式仅仅是作为例子来提示的，并不对本发明的保护范围进行限定， 例如，上述涉及的归一化区域大小及个数、所取的采样点的数量、向量的维数等也仅仅是例示，根据实际要求的速度、定位精度等，也可以是其他的数量。当前第1页1 2 3