专利名称：：基于灰度信息的人眼精确定位方法
技术领域：
：本发明涉及一种基于灰度信息的人眼精确定位方法，属信息
技术领域：
。
背景技术：
：自动人脸识别系统的构建已经成为一个模式识别和计算机视觉中的热门研究领域。一般来说，自动人脸识别系统主要包括两个部分，即人脸检测和人脸识别，目前来看，两个部分都有了一些效果较好的方法，虽然人脸检测已经取得了比较好的效果，但是由于人脸存在各种大小和不同的旋转角度，使得我们在进行人脸识别前，必须对检测到的人脸图像进行归一化处理，在人脸各种面部特征中，眼睛是进行人脸归一化的首选特征，在人脸识别中，人眼定位的准确度对人脸识别的效果有着重要的影响，现有的人眼检测方法不能够精确的对眼睛进行定位，特别是在人脸图像分辨率较低，噪音干扰较大的情况下往往定位的误差较大，从而影响了后续的人脸识别系统的性能。
发明内容本发明的目的在于克服现有技术之不足，而提供一种基于灰度信息的人眼精确定位的方法。本发明的技术方案是基于灰度信息的人眼精确定位方法，包括人脸检测、灰度直方图分析、图像增强、Gabor小波滤波、聚类分析及邻域操作步骤；具体如下(l)人脸检测用Adaboost算法对输入的数字图像进行人脸检测操作，得到人脸的位置和尺寸；(2)灰度直方图分析:将检测到的人脸图像统一縮放到130X150像素大小进行灰度化操作，选择人脸图像中央位置的一个大小为31X13像素的矩形区域作为肤色样本分析的对象，设/(力J)是人脸图像中G'，J)处的灰度值，历是矩形区域的灰度平均值，如果"180,则在该矩形区域的灰度直方图中，找到[yzr10，抓10]的范围内出现次数最多的一个灰度值作为肤色的灰度值W，对人脸图像进行灰度直方图分析，确定人脸中肤肤色区域的灰度范围；(3)图像增强确定的肤色的灰度值W'，使用下式对图像进行增强操作，新的灰度图像在(力J)处的灰度值记为/'(力J)。0,/(/，7')</關，/咖=wVa255,，/(,.,力>/簡，0/"<formula>formulaseeoriginaldocumentpage5</formula>d=2c2,c2=2-35如果m》180，则不对图像进行增强操作，即尸U/U(4)Gabor小波滤波对于增强后的人脸图像，使用Gabor小波对图像进行滤波操作，使用公式为<formula>formulaseeoriginaldocumentpage5</formula>如婦。2式中，F(Z，力表示像素点的空域坐标，IIll表示范数运算，、=V，、=C、//'，A=w/8,^{0,1,...,4}和/^{0，1,..,7}分别表示Gabor滤波器的方向和尺度，在Gabor滤波后产生的40个实部图像中，选择^=0,^{2.3，4}，即水平方向上尺度最大的三个实部图像依照式<formula>formulaseeoriginaldocumentpage5</formula>进行合成得到最后的参考图像；(5)聚类分析采用K-Means聚类的方法对合成的参考图像进行分析，得到人眼候选窗口，二值化人眼候选窗口，对二值化后的每个人眼候选窗口中的白色区域进行分析，得到人眼窗口；(6)邻域操作使用两个邻域算子对人眼窗口进行扫描，确定瞳孔中心位置。本发明与现有技术相比，具有检测方法简单、眼睛定位精确的优点。图l为本发明流程图示意图。图2为人脸图像进行灰度直方图分析过程示意图。图3为图像增强过程示意图。图4为Gabor小波滤波后的参考图像。图5为采用K-Means聚类的方法对合成的参考图像进行分析，得到的人眼窗口。图6为用两个邻域算子对人眼窗口进行扫描，确定的瞳孔中心位置。具体实施例方式首先对输入的数字图像使用Adaboost算法进行人脸检测操作，得到人脸的位置和尺寸，把检测到的人脸统一縮放到130X150像素大小，对人脸图像进行灰度化操作，选择人脸图像中央位置的一个大小为31X13像素的矩形区域作为肤色样本分析的对象,设/(/,J)是人脸图像中(/，j)处的灰度值，历是矩形区域的灰度平均值，如果/^180，则在该矩形区域的灰度直方图中，找到[/zr10，/^10]的范围内出现次数最多的一个灰度值作为肤色的灰度值w',使用下式对图像进行增强操作，新的灰度图像在(i，j)处的灰度值记为尸"i>。/'('，y'》o,/(u')</m,n,/m,n=5〃^255,，/(/，/)>/max,/max=w7c2(/(/，;')-/mm)'255如果m^l80,则不对图像进行增强操作，即尸(J'，》=JG'，力，具体过程见图2。对于增强后的人脸图像，使用Gabor小波对图像进行滤波操作，使用公式为e、A'z-一式中，F(;r，》表示像素点的空域坐标，IIll表示范数运算，、,V，J,、ax/尸，^-^W8，^{0,1,..,4}和/^{0,1,..,7}分别表示Gabor滤波器的方向和尺度。在Gabor滤波后产生的40个实部图像中，v={2,3,4},即水平方向上尺度最大的三个实部图像依照下式进行合成得到最后的参考图像，见图3。A+《2"3=1《3=2《2=4仏根据人脸面部器官的分布规律，能粗略估计人眼区域的分布，设(T:A^和力^^分别代表参考图像的宽度和高度，左眼和右眼的坐标分别为e,0r,，和&(义2,/2),《=1.3/4，y产1.6力/5,jf2=2.8/4,y2=1.6///5，以这两点为中心，得到大小为0.32rX0.32fr的两个人眼预处理窗口El和E2,见图4。对每个人眼预处理窗口E，把它划分为三个子集，初始化每一个子集的聚类中心为ce"te/^min(》，ce^er2=max(》，ce/3ter3=(ce"ter,+ce/7ter2)/2。通过K—means算法找至lj新的聚类中心ce"tar,'，ce77t印'，和ce^e/:3'，把ce"t印'所表示的子集设定为白色，其它两个子集设定为黑色，如果一组白色像素组成的块满足如下条件之一，则该块被置为黑色。(1):块的宽度小于块的高度。(2):块的像素点总数小于7(7-10)。计算白块的数量/7，同时计算每个白块的中心位置，记为c,Or,》，QOr,》，…，y)。对这"个中心位置按纵坐标排序得到新的中心序列cZ(义，》，Cr,y)，…，0r,/)，这里c,'(》《…《c/(》，如果/=1,则e(x'力=Cl'Or,j)。如果">=2，则eCr,y)=<2'Or,/)。设得到的两只眼睛的中心位置分别为。Cr,y)和^0r,/)，其中e,0r,/)表示左眼中心位置，e"x,》表示右眼中心位置。/7,表示左眼预处理窗口中白块的数目，仏表示右眼预处理窗口中白块的数目，如果e,Cr,》和&"y)两点所形成的水平夹角大于30度，且/7,〉仏时，眼睛中心位置e,0r,y)被移除，/7,〈A时，眼睛中心位置er(义，力被移除，重新计算新的眼睛中心的位置，最后，以得到的眼睛位置eU,y)为中心，在人脸灰度图像上切取31X13像素大小的矩形窗口作为人眼窗口￡F，具体过程见图5。对眼睛窗口^T中的每一个像素，求该像素周围的3^3邻域的像素值的和，然后用和值替代该像素的灰度值，见图6。用该领域算子对^F扫描一遍，把新得到的人眼窗口图像记为AT(:,j,在)V/(3,3)中，使用同样的方法获得该图像的5X5邻域图像，记为7^(5,5)。找到#/(5.5)中具有最小灰度值的位置，记为A,、"力，则瞳孔中心位置/U》=权利要求1、一种基于灰度信息的人眼精确定位方法，其特征在于该定位方法包括人脸检测、灰度直方图分析、图像增强、Gabor小波滤波、聚类分析及邻域操作步骤；具体如下(1)人脸检测用Adaboost算法对输入的数字图像进行人脸检测操作，得到人脸的位置和尺寸；(2)灰度直方图分析将检测到的人脸图像统一缩放到130×150像素大小进行灰度化操作，选择人脸图像中央位置的一个大小为31×13像素的矩形区域作为肤色样本分析的对象，设I(i，j)是人脸图像中(i，j)处的灰度值，m是矩形区域的灰度平均值，如果m＜180，则在该矩形区域的灰度直方图中，找到[m-10，m+10]的范围内出现次数最多的一个灰度值作为肤色的灰度值si，对人脸图像进行灰度直方图分析，确定人脸中肤色区域的灰度范围；(3)图像增强确定的肤色的灰度值si，使用下式对图像进行增强操作，新的灰度图像在(i，j)处的灰度值记为I”(i，j)。<mathsid="math0001"num="0001"><math><![CDATA[<mrow><msup><mi>I</mi><mo>′</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='{'close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo><</mo><msub><mi>I</mi><mi>min</mi></msub><mo>,</mo><msub><mi>I</mi><mi>min</mi></msub><mo>=</mo><mi>si</mi><mo>/</mo><msub><mi>c</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>255</mn><mo>,</mo></mtd><mtd><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>></mo><msub><mi>I</mi><mi>max</mi></msub><mo>,</mo><msub><mi>I</mi><mi>max</mi></msub><mo>=</mo><mi>si</mi><mo>/</mo><msub><mi>c</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mfrac><mrow><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>I</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>·</mo><mn>255</mn></mrow><mrow><mo>(</mo><msub><mi>I</mi><mi>max</mi></msub><mo>-</mo><msub><mi>I</mi><mi>min</mi></msub><mo>)</mo></mrow></mfrac><mo>,</mo></mtd><mtd><mi>others</mi></mtd></mtr></mtable></mfenced></mrow>]]></math></maths>c1＝2c2，c2＝2.35如果m≥180，则不对图像进行增强操作，即I’(i，j)＝I(i，j)；(4)Gabor小波滤波对于增强后的人脸图像，使用Gabor小波对图像进行滤波操作，使用公式为<mathsid="math0002"num="0002"><math><![CDATA[<mrow><msub><mi>ψ</mi><mrow><mi>μ</mi><mo>,</mo><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mrow><mi>μ</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><msup><mi>σ</mi><mn>2</mn></msup></mfrac><msup><mi>e</mi><mrow><mo>(</mo><mo>-</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mrow><mi>μ</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><msup><mrow><mo>|</mo><mo>|</mo><mi>z</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>/</mo><msup><mrow><mn>2</mn><mi>σ</mi></mrow><mn>2</mn></msup><mo>)</mo></mrow></msup><mo>[</mo><msup><mi>e</mi><mrow><msub><mi>ik</mi><mrow><mi>μ</mi><mo>,</mo><mi>v</mi></mrow></msub><mi>z</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mi>σ</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mrow></msup><mo>]</mo></mrow>]]></math></maths>式中，z＝(x，y)表示像素点的空域坐标，||·||表示范数运算，id="icf0003"file="A2009100947640002C3.tif"wi="15"he="4"top="219"left="24"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>kv＝kmax/fv，id="icf0004"file="A2009100947640002C4.tif"wi="15"he="3"top="220"left="59"img-content="drawing"img-format="tif"orientation="portrait"inline="yes"/>v∈{0，1，…，4}和μ∈{0，1，…，7}分别表示Gabor滤波器的方向和尺度，在Gabor滤波后产生的40个实部图像中，选择μ＝0，v＝{2，3，4}，即水平方向上尺度最大的三个实部图像依照式G(x，y)＝q1f(x，y)*ψ2，0+q2f(x，y)*ψ3，0+q3f(x，y)*ψ4，0<mathsid="math0003"num="0003"><math><![CDATA[<mfencedopen='{'close=''><mtable><mtr><mtd><msub><mi>q</mi><mn>1</mn></msub><mo>+</mo><msub><mi>q</mi><mn>2</mn></msub><mo>+</mo><msub><mi>q</mi><mn>3</mn></msub><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><msub><mi>q</mi><mn>3</mn></msub><mo>=</mo><msub><mrow><mn>2</mn><mi>q</mi></mrow><mn>2</mn></msub><mo>=</mo><msub><mrow><mn>4</mn><mi>q</mi></mrow><mn>1</mn></msub></mtd></mtr></mtable></mfenced>]]></math></maths>进行合成得到最后的参考图像；(5)聚类分析采用K-Means聚类的方法对合成的参考图像进行分析，得到人眼候选窗口，二值化人眼候选窗口，对二值化后的每个人眼候选窗口中的白色区域进行分析，得到人眼窗口；(6)邻域操作使用两个邻域算子对人眼窗口进行扫描，确定瞳孔中心位置。全文摘要本发明涉及一种基于灰度信息的人眼精确定位方法，属信息
技术领域：
。方法为(1)人脸检测从输入图像中定位人脸区域；(2)灰度直方图分析对人脸图像进行灰度直方图分析，确定人脸中肤色区域的灰度范围；(3)图像增强对图像进行灰度调整，使眼部特征更加明显；(4)Gabor小波滤波对增强后的图像进行Gabor小波滤波操作，利用Gabor滤波后的实部图像进行合成，得到参考图像；(5)聚类分析采用K-Means聚类的方法对合成的参考图像进行分析，得到二值化人眼候选窗口，对二值化后的每个人眼候选窗口中的白色区域进行分析，确定人眼的大致位置，得到人眼窗口；(6)邻域操作使用两个邻域算子对人眼窗口进行扫描，确定瞳孔中心位置。本发明具有检测方法简单、眼睛定位精确的优点。文档编号G06T7/00GK101615292SQ20091009476公开日2009年12月30日申请日期2009年7月24日优先权日2009年7月24日发明者丹徐,钱志明申请人:云南大学