专利名称：人脸图像的光照处理方法
技术领域：
本发明属于图像处理领域，涉及一种图像的光照处理算法。
背景技术：
在过去的十几年里，人脸识别受到了极大的关注，已经在商业和执法部门得到了广泛应用，如刑事鉴定、信用卡识别、安全系统、现场监控等。在这些实际应用中，人脸识别的应用环境是多样的，由于光照条件、面部表情、姿势、角度、发型等外界条件的影响，人脸识别的过程中要求这些外界条件相对稳定才能达到较高的识别率。由于人脸的表面是一个三维结构，点光源容易形成一些很强的投影，从而导致一些人脸细节的丢失，许多常用的识别方法对光照很敏感，故光照变化对人脸识别的过程影响巨大。在过去的几年中，人脸识别对光照的处理已经有了很大的提高，但是在国际上举行的人脸识别供应商评测(Face Recognition Vendor Test, FRVT)和人脸识别挑战赛 (Face Recognition Grand Challenge)中的测试表明，光照变化仍然是人脸识别性能提高的一个很重要的阻碍。目前有许多方法用来解决光照问题。如Retinex方法，该方法采用高斯核平滑图像来估计光照部分，通过在对数域对原图与光照估计求差得到光照处理图像。其中的高斯滤波核F(x，y)为环绕中心点的高斯函数F(x,y) = Ke-(x^yW其中，σ为高斯标准差，K为归一化因子，K = 1/ Σ F(x，y)。由Retinex方法处理得到的图像可以较好的保留全局信息，但是由于由Retinex 方法缺乏多尺度分析，故不能更好的提取细节信息。后来，在Retinex 的基础上又发展了 LTV(Logarithmic Total Variation)方法。然而，虽然这两种方法不需要对光照进行假设，也不需要多个图像参与训练，但处理后的图像主要表现为全局性信息，相应的细节性信息有所不足。近年来，小波分析开始应用于人脸的光照处理° 如〈〈Multiscale facialstructure representation for face recognition under varyingillumination)). Zhang Taiping, Fang Bin, Yuan Yuan, etc. Pattern Recognition, 2009,42 (2) :251-258 中多尺度人脸结构表示(Multiscale Facial StructureRepresentation, MFSR)方法，该方法用基于软阈值法的小波去噪模型估计光照部分，并认为人脸的不变部分等于原图与去噪后的图像之差。虽然多尺度人脸结构表示 (MFSR)在光照数据库上取得了非常好的效果，但经过该方法光照处理所得的图像突出了人脸的细节信息，而图像中人脸的边缘、轮廓等大尺度特征信息不够清晰，满足不了人脸识别的需要。另外，因为人脸的特征主要表现为水平方向，所以水平边缘信息相对于垂直边缘信息对人脸识别有着更重要的作用。但Retinex采用的高斯滤波各方向上的性质相同，不能有效地同时区分水平边缘信息和垂直边缘信息。
发明内容
本发明所要解决的技术问题是提供一种人脸图像的光照处理方法。本方法采用水平边缘滤波核的单尺度Ret inex光照处理方法对人脸图片进行处理，突出了人脸图像的水平边缘，并将经不同光照处理方法所得的图像按各自权重系数进行融合，融合后的图像同时保留了融合前各方法的优点，从而提升了人脸识别的性能，加大了人脸图像的识别率。本发明公开了一种人脸图像的光照处理方法，包括如下步骤步骤1 输入人脸图像，并将人脸图像归一化；步骤2 将归一化后的人脸图像用不同光照处理方法分别进行光照处理；步骤3 将光照处理后的各图像进行融合，并输出融合图像。所述步骤2中人脸图像进行光照处理时，利用单尺度Retinex方法进行光照处理， 并采用水平边缘滤波核来突出人脸图像的水平边缘。所述水平边缘滤波核为F'(x,y) = K-e-x'na2其中，X，y分别为滤波核的坐标，μ为既定参数，σ为高斯标准差， r _1_
Σ Σ (β"2/2σ2 .一丨“彻2)为归一化因子，i、j分别为滤波核中的坐标。所述既定参数为0至2之间的浮点数。所述步骤3中各图像进行融合时，按光照处理时采用的各光照处理的方法权重系数进行融合。所述各图像进行融合前单位化。本发明提出了一种人脸图像的光照处理方法，采用水平边缘滤波核的单尺度Ret inex光照处理方法对人脸图像进行处理，突出了人脸的水平边缘细节信息，从而提高识别性能，并对不同光照预处理方法得到的图像进行融合，保留了图像的大尺度信息和细节信息，能够有效地提高识别性能。
图1为本发明人脸图像的光照处理方法的流程图；图加为单尺度Retinex光照处理方法的高斯滤波核图像；图2b为本发明人脸图像的光照处理方法的水平边缘滤波核图像；图3为本发明人脸图像的光照处理方法的滤波过程示意图；图如为本发明人脸图像的光照处理方法中采用水平边缘滤波核进行光照处理的实施例原图；图4b为采用标准的单尺度Retinex处理方法(SSR)得到的结果图；图如为本发明人脸图像的光照处理方法的实施例中采用水平边缘滤波核的单尺度Retinex处理方法(H-SSR)得到的结果图像；图fe为用标准Ret inex方法对人脸图像的处理效果；图恥为本发明人脸图像的光照处理方法中采用水平边缘滤波核的单尺度 Retinex方法(H-SSR)对人脸图像的处理效果图；图6a为本发明人脸图像的光照处理方法中各光照处理方法进行融合实施例的原
4图；图6b为本发明人脸图像的光照处理方法中采用水平边缘滤波核的单尺度 Retinex方法(H-SSR)对实施例的原图进行光照处理所得的结果图像；图6c为采用LTV方法对实施例原图进行光照处理所得的结果图像；图6d为采用人脸多尺度结构表示方法(MFSR)对实施例原图进行光照处理所得的结果图像；图7为本发明人脸图像的光照处理方法中实施例的光照处理方法处理后的图像进行融合的示意图；图8为本发明人脸图像的光照处理方法中实施例的采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)和人脸多尺度结构表示方法(MFSR)处理后的图像进行融合后得到的结果图像；图9为本发明人脸图像的光照处理方法中经不同光照处理方法所得图像的识别率曲线图。
具体实施例方式下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。本发明公开了一种人脸图像的光照处理方法，如图1所示，包括如下步骤步骤1 输入人脸图像，并将人脸图像归一化为固定大小；对输入的人脸图像进行处理，根据各人脸图像中双眼的位置截取人脸区域，将得到的人脸区域进行归一化，缩放为固定大小。步骤2 将归一化后的人脸图像用不同光照处理方法分别进行光照处理；采用不同的光照处理方法，将归一化后的各人脸图像分别进行处理。本实施例中， 用SSR表示标准的单尺度Retinex光照处理方法，并以LTV (LogarithmicTotal Variation) 方法和多尺度人脸结构表示(Multiscale Facial StructureRepresentation, MFSR)方法为参考分别对归一化后的人脸图像进行处理。标准的单尺度Retinex光照处理方法(SSR):根据朗伯反射模型，图像I (x，y)等于反射分量R与光照分量L的乘积I(x, y) = R(x, y) · L(x, y)(1)由于光照分量L相对于反射分量R变化缓慢，因此可以通过对人脸图像进行平滑滤波，计算卷积从而获得L。但反射分量R求解涉时及到除法，计算比较耗时，因此重新对图像点(χ，y)定义为R'(x, y) = log I(x, y) — \og[F(x, y) I(x, y)](2)R' (x，y)是标准的单尺度Retinex光照处理方法(SSR)的输出，I (x, y)为灰度图像，(8)代表卷积。通过log变换后，除法运算变为减法运算，从而可以极大的提高运算速度。另外，由于log变换本身在一定程度上也能消减光照效应，所以变换后的式子可以获得更好的光照处理效果。在O)式中，F(x，y)为滤波算子，表示环绕中心点的高斯函数=(3)其中，σ为一常量，代表高斯标准差。
本方法利用单尺度Retinex方法进行光照处理，并采用水平边缘滤波核代替高斯滤波核估计光照部分，以突出人脸图像的水平边缘。水平边缘滤波核定义为F'{x,y) = K-ex2'2^ -e^'2"2(4)其中，x，y为图像中当前点的坐标，μ为既定参数，σ为高斯标准差， κ__1_
Σ Σ (6~"2σ2 .一們为归一化因子，整张图像对应的常量，ι J表示整张图像中的各坐
标点，既定参数μ的值为ο至2之间的浮点数。当μ =1时，所得水平边缘滤波核如图2b 所示。如图2a、图2b所示，相对于图加中的高斯滤波核，图2b中的水平边缘滤波核在垂直方向上分配的权重更大，故图2b所示的图像中垂直方向上的灰度值更大。当如式O)中对处于水平边缘上的点计算卷积时，水平边缘滤波核对水平边缘外的像素点分配较大的权重。如图3所示，图像背景为白色，中间包含一条黑色水平边缘，分别用高斯滤波核和水平边缘滤波核对图像中心的白点计算卷积。虚框表示滤波核的作用范围，水平边缘滤波核在黑点的位置分配的权重较大，则虚框中得到的卷积与黑点处的灰度值更相近，而与虚框中心白点的灰度值差异较大。如果采用高斯滤波核，则虚框中得到的卷积接近白点位置的灰度值，而与黑点的灰度值差异较大。标准的单尺度Retinex光照处理方法(SSR)的输出为当前点灰度值与卷积结果的差，而采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)可以增大当前点灰度值与卷积结果的差，进而突出水平边缘。如图^、4b、k所示，图如为原图，带有一条黑色的水平边缘，图4b为标准的单尺度Retinex处理方法(SSR)处理得到的结果图像，图如为采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)得到的结果图像，标准的单尺度Retinex光照处理方法中处理后的水平边缘中含有一条较亮的线，而基于水平边缘滤波核处理得到的图像效果更接近于原图，识别率更高。如图5b、图5c所示，为对图如进行光照处理的结果图像，相比标准单尺度 Retinex方法(SSR)，采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)在处理后的图像上水平边缘表现更清晰，如左眉毛和嘴巴的左半部分。以公开数据集卡内基梅隆大学人脸表情数据库(CMU PIE)中的(iallery集中的68 张无光照图像作为参考对象，illumination集中的14 张有光照图像作为测试对象。用 SSR表示基于高斯核的单尺度Retinex光照处理方法，σ = 14 ；H-SSR代表基于水平边缘滤波核的单尺度Retinex光照处理方法，既定参数μ = 0，高斯标准差σ = 12。未经处理的原图识别率为19. 65%，经过光照处理后识别率有较大的提高，标准的单尺度Retinex光照处理方法(SSR)的识别率达到了 74. 23%。而采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)的识别率达到了 78. 43%，相比标准的单尺度Retinex 光照处理方法(SSR)提高了 4. 2个百分点，提高了人脸识别的性能。采用LTV方法和人脸多尺度结构表示(MFSR)方法分别对输入的人脸图像进行光照处理。步骤3 将光照处理后的各图像进行融合，并输出融合图像。如图6b、6c、6d所示，为不同光照处理方法的对图6a进行光照处理的效果图。由
6图中可以看出，通过人脸多尺度结构表示方法(MFSR)进行光照处理可以突出细节信息，但对人脸图像的大尺度特征保留效果较差；采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)对人脸图像的大尺度特征保留较好，但对细节信息保留较差；LTV方法主要提取的是小尺度特征，但与人脸多尺度结构表示方法(MFSR)所保留的细节信息存在较大的差异。因此，可以把采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)和LTV 方法处理得到的图像分别与人脸多尺度结构表示方法(MFSR)处理得到的图像进行融合， 融合后的识别效果更好。如图7所示，各图像进行融合时，按光照处理时采用的各光照处理方法的权重系数进行融合，融合的方式如下
R1, R2、Rnew分别为两种光照处理方法以及两种光照处理方法进行融合后输出图像的灰度，β为权重系数。其中，R1为采用水平边缘滤波核的单尺度Retinex光照处理方法(H-SSR)或者LTV光照处理方法的输出图像的灰度值，&为人脸多尺度结构表示方法 (MFSR)光照处理方法的输出图像的灰度值。β表示采用水平边缘滤波核的单尺度Retinex 光照处理方法(H-SSR)或者LTV光照处理方法对应的权重系数，显然，当β =0时，融合后输出图像的灰度值等同于人脸多尺度结构表示方法(MFSR)进行光照处理输出图像的灰度值，当β = 1时，融合后输出图像的灰度值等同于采用水平边缘滤波核的单尺度Retinex 光照处理方法(H-SSR)或者LTV光照处理方法输出图像的灰度值。R1, R2在进行融合前需要归一化为单位向量，即
权利要求
1.一种人脸图像的光照处理方法，其特征在于，包括如下步骤 步骤1 输入人脸图像，并将人脸图像归一化；步骤2 将归一化后的人脸图像用不同光照处理方法分别进行光照处理； 步骤3 将光照处理后的各图像进行融合，并输出融合图像。
2.根据权利要求1所述的方法，其特征在于所述步骤2中人脸图像进行光照处理时， 利用单尺度Retinex方法进行光照处理，并采用水平边缘滤波核，突出人脸图像的水平边缘。
3.根据权利要求2所述的方法，其特征在于所述水平边缘滤波核为其中，x，y为当前点的坐标，μ为既定参数，σ为高斯标准差，尤=J Σ -e^"'2^)为归一化因子，i、j为整张图像中各点的坐标。
4.根据权利要求3所述的方法，其特征在于所述既定参数为0至2之间的浮点数。
5.根据权利要求1所述的方法，其特征在于所述步骤3中各图像进行融合时，按光照处理时采用的各光照处理方法权重系数进行融合。
6.根据权利要求5所述的方法，其特征在于所述各图像进行融合前单位化。
全文摘要
本发明一种人脸图像的光照处理方法，属于图像处理领域。包括如下步骤步骤1输入人脸图像，并将人脸图像归一化为固定大小；步骤2将归一化后的人脸图像用不同光照处理方法分别进行光照处理；步骤3将光照处理后的各图像进行融合，并输出融合图像。本方法采用水平边缘滤波核的单尺度Retinex光照处理方法对人脸图片进行处理，突出了人脸图像的水平边缘，并将经不同光照处理方法所得的图像按各自权重系数进行融合，融合后的图像同时保留了融合前各方法的优点，从而提升了人脸识别的性能，加大了人脸的识别率。
文档编号G06K9/54GK102214292SQ201010146028
公开日2011年10月12日 申请日期2010年4月9日 优先权日2010年4月9日
发明者刘昌平, 谭怒涛, 黄磊 申请人:汉王科技股份有限公司