CN102968636A - 一种人脸轮廓的提取方法 - Google Patents
一种人脸轮廓的提取方法 Download PDFInfo
- Publication number
- CN102968636A CN102968636A CN2012105294233A CN201210529423A CN102968636A CN 102968636 A CN102968636 A CN 102968636A CN 2012105294233 A CN2012105294233 A CN 2012105294233A CN 201210529423 A CN201210529423 A CN 201210529423A CN 102968636 A CN102968636 A CN 102968636A
- Authority
- CN
- China
- Prior art keywords
- ray
- facial contour
- point
- extracting method
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种人脸轮廓的提取方法,通过标准化要处理的图像尺寸,以降低后续处理的难度,并根据三庭五眼的人脸自然属性,依赖于人脸的平均轮廓,降低光照对人脸轮廓提取的影响。进而通过过鼻尖的射线与人脸轮廓的交点,然后查找分界点,进而拟合人脸轮廓,精确性更高。
Description
技术领域
本发明涉及一种人脸轮廓的提取方法。
背景技术
人脸轮廓线提取是人脸特征识别和人头三维模型重建的基础。在大多数应用中,人脸,尤其是脸型识别都是以正面人脸轮廓为基础进行的识别。更具体地,本文中提取的人脸轮廓是指的是两耳到下巴的外轮廓线,该轮廓可以用来进行人脸脸型的识别,是人脸识别的一个重要特征依据。
人脸轮廓的提取方法,通常采用几种方法:
方法一、使用Canny算子的边缘检测算法,该方法使用全局阈值对图像进行二值化,进而提取人脸轮廓。当所输入的原始图像光照不均匀时,图像上各像素点亮度差异较大,全局阈值显然难以胜任所有位置的边缘检测。由此可知,该方法极易受光照影响,人脸光照不均匀时,不能准确的找到轮廓。
方法二、使用基于肤色模型和梯度算子的方法,提取人脸轮廓,该类方法利用了肤色信息能够很好的提取人脸轮廓,但不能处理灰度图像。
方法三、基于几何活动轮廓模型的人脸轮廓提取方法,利用人脸形状的椭圆性约束来寻找人脸轮廓,这种方法,能较精确地提取出图象中的人脸轮廓。但模型的初始位置直接关系到结果的准确性,容易陷入局部最优解。
发明内容
本发明的目的在于提出一种新的人脸轮廓的提取方法,利用人脸平均轮廓,有效的降低光照的影响,并准确的提取人脸轮廓。
本发明采用以下技术方案:
一种人脸轮廓的提取方法,包括以下步骤:
1)对要处理的图像进行尺寸标准化:规范两瞳孔间距以及两瞳孔的平均中心点坐标,并统一图像尺寸;
2)根据三庭五眼法确定鼻尖位置:依据两瞳孔间距及两瞳孔的平均中心点坐标确定鼻尖A的平均坐标(x0,y0);
3)过所述鼻尖A向两侧做水平射线K0,进而以K0为基轴,每隔角度m作一条射线Ki,对应有N = 180/m +1条射线,从而射线Ki的方程为y = (x-x0)tan(i*m) + y0,其中0 ≤i ≤N;
4)根据三庭五眼法粗确定人脸轮廓:依据两瞳孔间距及两瞳孔的平均中心点坐标确定左耳垂B的平均坐标为(x2,y2),右耳垂D的平均坐标为(x3,y3),下巴中心点C的平均坐标为(x1,y1),从而过点B、C、D的人脸平均轮廓所对应的方程为y = -( (x-x1) * (x-x1) ) / (2*P) + y1,其中P = (x3 - x2)*( x3 - x2) / ( 8 * (y1-y3) );
5)依据每一射线Ki的方程与粗获得的人脸轮廓的交点集合M,寻找每一射线Ki上人脸轮廓分界点;
6)集合所有分界点,拟合成人脸轮廓。
从以上的技术方案可以看出,依据本发明的上述人脸轮廓的提取方法通过标准化要处理的图像尺寸,以降低后续处理的难度,并根据三庭五眼的人脸自然属性,依赖于人脸的平均轮廓,降低光照对人脸轮廓提取的影响。进而通过过鼻尖的射线与人脸轮廓的交点,然后查找分界点,进而拟合人脸轮廓,精确性更高。
上述人脸轮廓的提取方法,进一步保证的方案是,步骤5)中确定每一射线Ki上人脸轮廓分界点的方法是:
5.1)构建射线Ki上第j个点的向量Vi(j) = [ Ki (j-5)- Ki (j),Ki (j-4)- Ki (j),Ki (j-3)- Ki (j),Ki (j-2)- Ki (j),Ki (j-1)- Ki (j),Ki (j+1)- Ki (j),Ki (j+2)- Ki (j),Ki (j+3)- Ki (j),Ki (j+4)- Ki (j),Ki (j+5)- Ki (j) ],其中Ki (j) (0 < j < R)是射线Ki上第j个点的灰度值,R取A点到射线Ki上mi点的像素长度加上微调参数L;
5.2)计算向量Vi(j)的方差D i(j),D i(j)是第i条射线上第j个点的向量方差;
5.3)将射线Ki上方差D i(j)最大的前3个点,与集合M对应的交点mi求欧式距离,距离最小的点即为对应射线Ki的边界点。
这样,通过灰度方差,不仅可以准确的射线界定Ki的边界点,而且能够更有效地降低光照对人脸轮廓提取的影响。
上述人脸轮廓的提取方法,在步骤2)之前还包括对要处理的图像进行灰度化的步骤,以减小后续步骤处理的计算量。
上述人脸轮廓的提取方法,对图像进行灰度化的方法是取像素点灰度值Y = R×0.299 + G×0.587 + B×0.114,其中R、G、B为彩色图像三个通道的亮度值。
上述人脸轮廓的提取方法,角度m为5o,使得拟合出的人脸轮廓在相对平滑的同时,计算量也不会很大。
上述人脸轮廓的提取方法,使用最小二乘法多项式进行曲线拟合,多项式方程为f(x)=a0+a1x+ a2x2 + … +anxn,其中,多项式的阶数n为9,a0,a1 ,a2… an是拟合后的多项式系数。
具体实施方式
依据本发明的实施例表述如下,为更清楚的进行表达,以树形结构进行描述:
1、对要处理的图像进行尺寸标准化处理;以降低后续处理的难度,尤其是涉及相关计算的处理;注:后续内容中分号为相对应步骤的说明。
1.1、根据1,对图片进行缩放、旋转、平移、剪裁,使得照片尺寸统一为240×320,两瞳孔距离为64像素,两瞳孔的平均中心点坐标为(120 ,160);在该技术条件下,可以获得相对清晰的脸部轮廓,且计算量不大。
2、图像灰度化:判定图像为非灰度图像,则将图像灰度化。
2.1、根据2,可用公式Y = R×0.299 + G×0.587 + B×0.114,把RGB 彩色图像转换为灰度图像。
3、根据人眼间距及人脸五官比例,推测鼻尖A的平均坐标为 (x0,y0),过鼻尖A向脸颊两侧做水平射线K0,过A点,距K0间隔角度m--做射线Ki,共有N = 180/m +1条射线,射线Ki的方程为y = (x-x0)tan(i*5) + y0,其中0 ≤i ≤N。
人脸五官比例体现在人脸的相似性,表现在三庭五眼法,在计算机图形学中称为三庭五眼算法,也成三庭五眼黄金分割法。其原理在于人的脸长与脸宽的一般标准比例,不符合此比例,就会与理想的脸型产生距离。眼睛的宽度,应为同一水平脸部宽度的3/10;下巴长度应为脸长的1/5;眼球中心到眉毛底部的距离,应为脸长的1/10;眼球应为脸长的1/14;鼻子的表面积,要小于脸部总面积的5/100;理想嘴巴宽度应为同一水平脸部宽度的1/2。通过一些算法可以在知晓其中的一些量之后可以粗估其他部分的位置和尺寸。三庭和五眼具有位置的相对固定性,并且一些人具有严格的三庭五眼特征,符合黄金分割比例,大多数人在此基础上差别并不大。
3.1、根据3,鼻尖A的平均坐标可以是(120,205),角度m可以是5,则射线Ki的方程为y = (x-120)tan(i*m) + 205,其中0 ≤ i ≤37,tan表示正切函数。
4、根据人眼间距及人脸五官比例,推测左耳垂B的平均坐标为(x2,y2),右耳垂D的平均坐标为(x3,y3),下巴中心点C的平均坐标为(x1,y1),过点B、C、D的人脸平均轮廓方程为y = -( (x-x1) * (x-x1) ) / (2*P) + y1,其中P = (x3 - x2)*( x3 - x2) / ( 8 * (y1-y3) );这里的平均就主要表现为三庭五眼黄金分割的情形,是人脸的一般属性;那么依据该轮廓与上述射线的交点就,以及边缘点的属性就可以找出准确的人脸轮廓。
4.1、根据4、左耳垂B的平均坐标可以是(47,210),右耳垂D的平均坐标可以是(193,210),下巴中心点C的平均坐标可以是(120,277),P = (x3 - x2)*( x3 - x2) / ( 8 * (y1-y3) ) = (193 – 47)*( 193 – 47) / ( 8 * (277-210) ) = 39.77。人脸平均轮廓方程y = -( (x-x1) * (x-x1) ) / (2*P) + y1 = -( (x-120) * (x-120) ) / (2*39.77) + 277 = -( (x-120) * (x-120) ) /79.54 + 277。
5、求人脸平均轮廓方程和射线Ki的交点集合M ={m0,m1,m2,…,mi}是射线Ki上的人脸平均轮廓边界点集合,其中0 ≤ i ≤ (180/m +1)。
6、寻找射线Ki上的人脸轮廓边界点:以A为起点向脸颊两侧方向,依次遍历每条射线寻找人脸轮廓边界点。
6.1、构建射线Ki上第j个点的向量Vi(j) = [ Ki (j-5)- Ki (j),Ki (j-4)- Ki (j),Ki (j-3)- Ki (j),Ki (j-2)- Ki (j),Ki (j-1)- Ki (j),Ki (j+1)- Ki (j),Ki (j+2)- Ki (j),Ki (j+3)- Ki (j),Ki (j+4)- Ki (j),Ki (j+5)- Ki (j) ],其中Ki (j) (0 < j < R)是射线Ki上第j个点的灰度值,R取A点到射线Ki上mi点的像素长度加上微调参数L。
6.1.1、根据6.1,R的调整微调L = 10。
6.2、计算向量Vi(j)的方差D i(j),D i(j)是第i条射线上第j个点的向量方差。
6.3、将射线Ki上方差D i(j)最大的前3个点,与集合M对应的交点mi求欧式距离,距离最小的即为要寻找的边界点。
7、对每条射线Ki上找到的人脸轮廓边界点,使用最小二乘法多项式进行曲线拟合,得到人脸轮廓曲线方程。
7.1、根据7,拟合的多项式是f(x)=a0+a1x+ a2x2 + … +anxn,其中,n是多项式的阶数,这里取n = 9;a0,a1 ,a2… an是拟合后的多项式系数。
Claims (6)
1.一种人脸轮廓的提取方法,其特征在于,包括以下步骤:
1)对要处理的图像进行尺寸标准化:规范两瞳孔间距以及两瞳孔的平均中心点坐标,并统一图像尺寸;
2)根据三庭五眼法确定鼻尖位置:依据两瞳孔间距及两瞳孔的平均中心点坐标确定鼻尖A的平均坐标(x0,y0);
3)过所述鼻尖A向两侧做水平射线K0,进而以K0为基轴,每隔角度m作一条射线Ki,对应有N = 180/m +1条射线,从而射线Ki的方程为y = (x-x0)tan(i*m) + y0,其中0 ≤i ≤N;
4)根据三庭五眼法粗确定人脸轮廓:依据两瞳孔间距及两瞳孔的平均中心点坐标确定左耳垂B的平均坐标为(x2,y2),右耳垂D的平均坐标为(x3,y3),下巴中心点C的平均坐标为(x1,y1),从而过点B、C、D的人脸平均轮廓所对应的方程为y = -( (x-x1) * (x-x1) ) / (2*P) + y1,其中P = (x3 - x2)*( x3 - x2) / ( 8 * (y1-y3) );
5)依据每一射线Ki的方程与粗获得的人脸轮廓的交点集合M,寻找每一射线Ki上人脸轮廓分界点;
6)集合所有分界点,拟合成人脸轮廓。
2.根据权利要求1所述的人脸轮廓的提取方法,其特征在于,步骤5)中确定每一射线Ki上人脸轮廓分界点的方法是:
5.1)构建射线Ki上第j个点的向量Vi(j) = [ Ki (j-5)- Ki (j),Ki (j-4)- Ki (j),Ki (j-3)- Ki (j),Ki (j-2)- Ki (j),Ki (j-1)- Ki (j),Ki (j+1)- Ki (j),Ki (j+2)- Ki (j),Ki (j+3)- Ki (j),Ki (j+4)- Ki (j),Ki (j+5)- Ki (j) ],其中Ki (j) (0 < j < R)是射线Ki上第j个点的灰度值,R取A点到射线Ki上mi点的像素长度加上微调参数L;
5.2)计算向量Vi(j)的方差D i(j),D i(j)是第i条射线上第j个点的向量方差;
5.3)将射线Ki上方差D i(j)最大的前3个点,与集合M对应的交点mi求欧式距离,距离最小的点即为对应射线Ki的边界点。
3.根据权利要求1或2所述的人脸轮廓的提取方法,其特征在于,在步骤2)之前还包括对要处理的图像进行灰度化的步骤。
4.根据权利要求3所述的人脸轮廓的提取方法,其特征在于,对图像进行灰度化的方法是取像素点灰度值Y = R×0.299 + G×0.587 + B×0.114,其中R、G、B为彩色图像三个通道的亮度值。
5.根据权利要求1或2所述的人脸轮廓的提取方法,其特征在于,角度m为5o。
6.根据权利要求1所述的人脸轮廓的提取方法,其特征在于,使用最小二乘法多项式进行曲线拟合,多项式方程为f(x)=a0+a1x+ a2x2 + … +anxn,其中,多项式的阶数n为9,a0,a1 ,a2… an是拟合后的多项式系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210529423.3A CN102968636B (zh) | 2012-12-11 | 2012-12-11 | 一种人脸轮廓的提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210529423.3A CN102968636B (zh) | 2012-12-11 | 2012-12-11 | 一种人脸轮廓的提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102968636A true CN102968636A (zh) | 2013-03-13 |
CN102968636B CN102968636B (zh) | 2015-09-09 |
Family
ID=47798767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210529423.3A Active CN102968636B (zh) | 2012-12-11 | 2012-12-11 | 一种人脸轮廓的提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102968636B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016074248A1 (zh) * | 2014-11-15 | 2016-05-19 | 深圳市三木通信技术有限公司 | 基于人脸识别的验证使用方法及装置 |
CN110348286A (zh) * | 2019-05-24 | 2019-10-18 | 广东工业大学 | 一种基于最小二乘法的脸型拟合与匹配方法 |
CN110909571A (zh) * | 2018-09-17 | 2020-03-24 | 武汉资联虹康科技股份有限公司 | 一种高精度面部识别空间定位方法 |
CN112257721A (zh) * | 2020-12-08 | 2021-01-22 | 常州微亿智造科技有限公司 | 一种基于Fast ICP的图像目标区域匹配方法 |
CN113591602A (zh) * | 2021-07-08 | 2021-11-02 | 娄浩哲 | 一种基于单视角的人脸三维轮廓特征重建装置及重建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1731418A (zh) * | 2005-08-19 | 2006-02-08 | 清华大学 | 复杂背景图像中鲁棒的眼睛精确定位方法 |
CN101295405A (zh) * | 2008-06-13 | 2008-10-29 | 西北工业大学 | 人像与车辆识别报警跟踪方法 |
CN101464948A (zh) * | 2009-01-14 | 2009-06-24 | 北京航空航天大学 | 一种基于关键点的仿射不变矩的目标识别方法 |
CN101650777A (zh) * | 2009-09-07 | 2010-02-17 | 东南大学 | 一种基于密集点对应的快速三维人脸识别方法 |
-
2012
- 2012-12-11 CN CN201210529423.3A patent/CN102968636B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1731418A (zh) * | 2005-08-19 | 2006-02-08 | 清华大学 | 复杂背景图像中鲁棒的眼睛精确定位方法 |
CN101295405A (zh) * | 2008-06-13 | 2008-10-29 | 西北工业大学 | 人像与车辆识别报警跟踪方法 |
CN101464948A (zh) * | 2009-01-14 | 2009-06-24 | 北京航空航天大学 | 一种基于关键点的仿射不变矩的目标识别方法 |
CN101650777A (zh) * | 2009-09-07 | 2010-02-17 | 东南大学 | 一种基于密集点对应的快速三维人脸识别方法 |
Non-Patent Citations (1)
Title |
---|
陈其杰: "《人脸外轮廓线的提取方法》", 《红外与激光工程》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016074248A1 (zh) * | 2014-11-15 | 2016-05-19 | 深圳市三木通信技术有限公司 | 基于人脸识别的验证使用方法及装置 |
CN110909571A (zh) * | 2018-09-17 | 2020-03-24 | 武汉资联虹康科技股份有限公司 | 一种高精度面部识别空间定位方法 |
CN110909571B (zh) * | 2018-09-17 | 2022-05-03 | 武汉资联虹康科技股份有限公司 | 一种高精度面部识别空间定位方法 |
CN110348286A (zh) * | 2019-05-24 | 2019-10-18 | 广东工业大学 | 一种基于最小二乘法的脸型拟合与匹配方法 |
CN112257721A (zh) * | 2020-12-08 | 2021-01-22 | 常州微亿智造科技有限公司 | 一种基于Fast ICP的图像目标区域匹配方法 |
CN113591602A (zh) * | 2021-07-08 | 2021-11-02 | 娄浩哲 | 一种基于单视角的人脸三维轮廓特征重建装置及重建方法 |
CN113591602B (zh) * | 2021-07-08 | 2024-04-30 | 娄浩哲 | 一种基于单视角的人脸三维轮廓特征重建装置及重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102968636B (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103456010B (zh) | 一种基于特征点定位的人脸卡通画生成方法 | |
CN107403168B (zh) | 一种面部识别系统 | |
CN106250867B (zh) | 一种基于深度数据的人体骨骼跟踪系统的实现方法 | |
CN105139438B (zh) | 视频人脸卡通动画生成方法 | |
CN103761519B (zh) | 一种基于自适应校准的非接触式视线追踪方法 | |
CN104573634A (zh) | 一种三维人脸识别方法 | |
Gu et al. | Feature points extraction from faces | |
CN107122789B (zh) | 基于深度摄像头的多模态信息融合的学习专注度分析方法 | |
CN103839250B (zh) | 用于面部图像处理的方法和设备 | |
CN102737235B (zh) | 基于深度信息和彩色图像的头部姿势估计方法 | |
CN106504751A (zh) | 自适应唇语交互方法以及交互装置 | |
CN102968636B (zh) | 一种人脸轮廓的提取方法 | |
CN104408462B (zh) | 面部特征点快速定位方法 | |
CN104598915A (zh) | 一种手势识别方法与装置 | |
CN106774856A (zh) | 基于唇语的交互方法以及交互装置 | |
CN103218605A (zh) | 一种基于积分投影与边缘检测的快速人眼定位方法 | |
CN105913416A (zh) | 一种自动分割三维人脸模型区域的方法 | |
CN103679175A (zh) | 一种基于深度摄像机的快速3d骨骼模型检测方法 | |
CN103870808A (zh) | 一种手指静脉识别方法 | |
CN103810491A (zh) | 融合深度和灰度图像特征点的头姿估计兴趣点检测方法 | |
CN106600640A (zh) | 一种基于rgb‑d相机的人脸识别辅助眼镜 | |
TW200719871A (en) | A real-time face detection under complex backgrounds | |
CN106874830A (zh) | 一种基于rgb‑d相机和人脸识别的视障人士辅助方法 | |
CN102024156A (zh) | 彩色人脸图像中的唇部区域定位方法 | |
CN103544478A (zh) | 一种全方位人脸检测的方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C56 | Change in the name or address of the patentee | ||
CP01 | Change in the name or title of a patent holder |
Address after: Shun high tech Zone of Ji'nan City, Shandong province 250101 China West Road No. 699 Patentee after: SYNTHESIS ELECTRONIC TECHNOLOGY CO., LTD. Address before: Shun high tech Zone of Ji'nan City, Shandong province 250101 China West Road No. 699 Patentee before: Shandong Synthesis Electronic Technology Co., Ltd. |