一种利用对数变换和Kirsch算子的人脸光照不变特征提取方法与流程

本发明涉及模式识别技术领域，具体是一种利用对数变换和kirsch算子的人脸光照不变特征提取方法。

背景技术：

基于生物特征之一的人脸作为人体与生俱来的最普通的生理结构，具有其他生物特征所不具有的直接友好性，是生物特征识别研究领域的重要对象。因此，基于人脸的识别技术应运而生，它作为一种最能平衡各方利弊及具多功能的生物特征识别技术，在身份鉴定、人机交互和视频监控等商业应用领域中有着广阔的发展前景。经过近40年的研究和发展，取得了长足的进步，许多优秀的人脸识别算法在理想情况下已经能取得非常好的识别性能。但是人脸识别仍是一个没有彻底解决的难题，光照变化问题便是其未解决的关键问题之一。解决人脸识别中光照问题的方法主要分为三类：提取光照不变特征、光照变化的建模、光照条件标准化。在提取光照不变特征的方法中，又主要分为两种：第一种是在对数域，通过低通滤波消除光照分量，例如：在对数域中，利用离散余弦变换(dct)或小波变换进行滤波处理等方法。第二种是构造除法运算，通过除法运算来消除慢变化的光照分量，例如：梯度脸(gradient-face)、韦伯脸(weber-face)、局部二值模式(localbinarypattern，lbp)等方法。因此能够有效地提取在光照中保持鲁棒性的人脸特征以增强人脸识别在实际中的应用成为了一个至关重要的问题。

技术实现要素：

本发明的目的在于提供一种能够提高人脸识别率的利用对数变换和kirsch算子的人脸光照不变特征提取方法，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：

一种利用对数变换和kirsch算子的人脸光照不变特征提取方法，包括以下步骤：

(1)将原图像变换到对数域；

(2)在对数域中，利用kirsch算子对人脸锐化处理，提取人脸的水平方向细节特征，同时消除了慢变化的人脸光照部分，所得结果就是需要提取的人脸光照不变特征。

作为本发明进一步的方案：根据光照反射模型，任意灰度人脸图像i(x,y)可以看成是反射分量r(x,y)和光照分量l(x,y)的乘积，即i(x,y)＝r(x,y)l(x,y)；所述步骤(1)中，对人脸图像进行对数变换，使其反射分量和光照分量由相乘转换为相加，即lni(x,y)＝lnr(x,y)+lnl(x,y)。

作为本发明进一步的方案：以kirsch算子模板进行锐化处理后，图像i'某像素点(x,y)的灰度值i'(x,y)如公式(1)所示：

i'(x,y)＝-3lni(x-1,y-1)-3lni(x,y-1)-3lni(x+1,y-1)-3lni(x-1,y)

-3lni(x+1,y)+5lni(x-1,y+1)+5lni(x,y+1)+5lni(x+1,y+1)(1)；

为了表述方便，i(x+m,y+n)简记为im,n，同理，r(x+m,y+n)简记为rm,n，l(x+m,y+n)简记为lm,n，则公式(1)简记为公式(2)，如下所示：

i'0,0＝-3lni-1,-1-3lni0,-1-3lni1,-1-3lni-1,0-3lni1,0+5lni-1,1+5lni0,1+5lni1,1(2)；

根据光照反射模型进一步得到公式(3)：

i'0,0＝-3lnr-1,-1l-1,-1-3lnr0,-1l0,-1-3lnr1,-1l1,-1-3lnr-1,0i-1,0

-3lnr1,0i1,0+5lnr-1,1r-1,1+5lnr0,1r0,1+5lnr1,1l1,1(3)；

作为本发明进一步的方案：光照分量l(x,y)变化缓慢，所以有公式(4)：

l-1,-1≈l0,-1≈l1,-1≈l-1,0≈l1,0≈l-1,1≈l0,1≈l1,1(4)；

对于公式(4)的所有变量均用lt替代，得到公式(5)：

因此消除了光照分量。

与现有技术相比，本发明的有益效果是：

本发明对不同角度的光照具有一定的抑制作用，处理后人脸图像的识别率较现有的算法处理后的人脸识别率高；本发明采用对数变换和kirsch算子对人脸图像做光照处理，方法简单，识别速度快。

附图说明

图1为本发明提取人脸光照不变特征的基本流程图。

图2为本发明中以(x,y)为中心的八邻域灰度值矩阵及其对应的简化表达示意图。

图3为提取cmupie人脸数据库光照不变特征前后的人脸示例图。

图4为提取扩展yaleb正面人脸库光照不变特征前后的人脸示例图。

具体实施方式

下面结合具体实施方式对本专利的技术方案作进一步详细地说明。

请参阅图1-4，一种利用对数变换和kirsch算子的人脸光照不变特征提取方法，包括以下步骤：

(1)将原图像变换到对数域；

(2)在对数域中，利用kirsch算子对人脸锐化处理，提取人脸的水平方向细节特征，同时消除了慢变化的人脸光照部分，所得结果就是需要提取的人脸光照不变特征。

根据光照反射模型，任意灰度人脸图像i(x,y)可以看成是反射分量r(x,y)和光照分量l(x,y)的乘积，即i(x,y)＝r(x,y)l(x,y)；所述步骤(1)中，对人脸图像进行对数变换，使其反射分量和光照分量由相乘转换为相加，即lni(x,y)＝lnr(x,y)+lnl(x,y)。

以kirsch算子模板进行锐化处理后，图像i'某像素点(x,y)的灰度值为i'(x,y)。

具体来说，所述利用对数变换和kirsch算子的人脸光照不变特征提取方法，具体步骤如下：

(1)以原图像某一像素点(x,y)为中心的8邻域如图2(a)所示，为简化表述，中心像素点(x,y)处的灰度值i(x,y)表述为i0,0，如图2(b)所示，其邻域点的灰度值i(x+m,y+n)表述为im,n。

根据光照反射模型，任意灰度人脸图像i(x,y)可以看成是反射分量r(x,y)和光照分量l(x,y)的乘积，即i(x,y)＝r(x,y)l(x,y)，其中l(x,y)取决于光源且只代表面部光照的分量，r(x,y)取决于物体的表面特征并包含了面部的关键信息，在本发明中i(x+m,y+n)可表述为rm,nlm,n。将人脸图像每个像素值变换到对数域，在本发明中lni(x+m,y+n)表述为lnrm，n+lnlm，n。

(2)kirsch微分算子可以对人脸图像进行八个特定方向上的锐化处理，加强图像的细节边缘和轮廓特征，这八个方向上的模板分别为

提取图像水平方向细节特征的算子主要是h0和h4算子，本发明在下面的推导过程中，仅以h4模板为例，对人脸图像进行锐化处理，得到i'。具体的原理推导过程如下：

设以kirsch算子的h4矩阵模板进行锐化处理后，人脸图像i'某像素点(x,y)的灰度值i'(x,y)如公式(1)所示：

i'(x,y)＝-3lni(x-1,y-1)-3lni(x,y-1)-3lni(x+1,y-1)-3lni(x-1,y)

-3lni(x+1,y)+5lni(x-1,y+1)+5lni(x,y+1)+5lni(x+1,y+1)(1)

为了表述方便，i(x+m,y+n)简记为im,n，同理，r(x+m,y+n)简记为rm,n，l(x+m,y+n)简记为lm,n，则公式(1)简记为公式(2)，如下所示：

i'0,0＝-3lni-1,-1-3lni0,-1-3lni1,-1-3lni-1,0-3lni1,0+5lni-1,1+5lni0,1+5lni1,1(2)

根据光照反射模型进一步得到公式(3)：

i'0,0＝-3lnr-1,-1l-1,-1-3lnr0,-1l0,-1-3lnr1,-1l1,-1-3lnr-1,0i-1,0

-3lnr1,0i1,0+5lnr-1,1r-1,1+5lnr0,1r0,1+5lnr1,1l1,1(3)

因为光照分量l(x,y)变化缓慢，在局部近似相等，所以有公式(4)：

l-1,-1≈l0,-1≈l1,-1≈l-1,0≈l1,0≈l-1,1≈l0,1≈l1,1(4)

对于公式(4)的所有变量均用lt替代，得到公式(5)：

由公式(5)可知已消除了光照分量。

本发明在cmupie和扩展yaleb人脸库中进行实验：

1、人脸库介绍：cmupie人脸数据库共有68个人的41368张人脸图像，本发明选择的是在正面人脸光照组(c27)中进行实验，c27光照组共包含1428个人脸库。实验时，使用的人脸尺寸为64×64，从每个人脸的21种不同光照中选取一张图像作为训练样本，其余人脸作为测试样本。

扩展yaleb人脸库共包含38个人的9种不同姿态，每种姿态又包含64种不同的光照情况。本发明只在扩展yaleb正面人脸库(共2432张人脸)中进行实验，所有人脸图像按照不同角度可分为5个子集：子集1(θ＜12°)有266个样本、子集2(13°＜θ＜25°)有456个样本、子集3(26°＜θ＜50°)有456个样本、子集4(51°＜θ＜77°)有532个样本、子集5(θ＞77°)有722个样本。实验时，使用的人脸尺寸为192×168，用子集1作为训练样本，其余子集作为测试样本。

2、实验结果：

实验1，将提取的人脸光照不变特征与原图像进行了对比发现经本发明方法处理后的人脸图像特征更为明显，有效地克服了光照对人脸特征的影响。在cmupie正面人脸库和扩展yaleb正面人脸库中提取光照不变特征前后的人脸示例分别如图3、图4所示。

实验2，本发明方法采用l1范数测量距离后利用最近邻分类器分类判别。在cmupie人脸库上的正确识别率是97.78％，在扩展yaleb人脸库的人脸正确识别率是98.80％，对人脸不变特征提取有较好的鲁棒性。

本发明对不同角度的光照具有一定的抑制作用，处理后人脸图像的识别率较现有的算法处理后的人脸识别率高；本发明采用对数运算和kirsch算子对人脸图像做光照处理，方法简单，识别速度快。

上面对本专利的较佳实施方式作了详细说明，但是本专利并不限于上述实施方式，在本领域的普通技术人员所具备的知识范围内，还可以在不脱离本专利宗旨的前提下做出各种变化。