CN110110637A - 一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法 - Google Patents

一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法 Download PDF

Info

Publication number
CN110110637A
CN110110637A CN201910352167.7A CN201910352167A CN110110637A CN 110110637 A CN110110637 A CN 110110637A CN 201910352167 A CN201910352167 A CN 201910352167A CN 110110637 A CN110110637 A CN 110110637A
Authority
CN
China
Prior art keywords
wrinkle
image
face
value
texture
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.)
Pending
Application number
CN201910352167.7A
Other languages
English (en)
Inventor
李嘉路
华芮
李骁
张游龙
王润苗
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Huajia Biological Intelligence Technology Co Ltd
Original Assignee
Shenzhen Huajia Biological Intelligence Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Huajia Biological Intelligence Technology Co Ltd filed Critical Shenzhen Huajia Biological Intelligence Technology Co Ltd
Priority to CN201910352167.7A priority Critical patent/CN110110637A/zh
Publication of CN110110637A publication Critical patent/CN110110637A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/161Detection; Localisation; Normalisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/161Detection; Localisation; Normalisation
    • G06V40/165Detection; Localisation; Normalisation using facial parts and geometric relationships
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/168Feature extraction; Face representation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/10Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
    • G06V40/16Human faces, e.g. facial parts, sketches or expressions
    • G06V40/172Classification, e.g. identification
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Computer Interaction (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开了一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法。本发明涉及人工智能的计算机视觉领域,可基于一张清晰的人脸照片,依次采用人脸识别、图像分割、皱纹识别、特征提取和预测分级的方法,自动识别人脸皮肤图像中皱纹并对严重程度自动分级,自动分级的准确度可高达0.85。该方法通过计算机全自动快速实现,具备专业人员对皱纹评估的准确度,可广泛应用于医美、美容和化妆品等行业,也可辅助与皮肤皱纹相关的研究。

Description

一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法
技术领域
本发明采用的方法涉及人脸识别、图像分割、皱纹识别、特征提取与皱纹自动分级,涉及人工智能的计算机视觉领域和生物统计学领域。
背景技术
人类可以通过肉眼自动去寻找人脸或是人脸图像中的皱纹区域,并对观察到的皱纹进行主观的严重程度评价。但受主观情感的影响和专业知识的限制,普通人对皱纹严重程度的评价往往不够准确;专业医师的评价往往较为合理,但受时间、地点等环境因素的限制,邀请专业医师进行评判并不现实可行。因此,亟需一种方便而准确的皱纹严重程度自动化评估方法,可以让普通群众快速便捷的得到合理的评估结果,同时也可辅助皮肤相关的研究或研发。
发明内容
针对以上难题,我们发明了一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,该方法对人脸图像的进行人脸识别、图像分割、皱纹识别和特征提取,再结合机器学习(machine learning)的方法对皱纹进行分级预测。本发明方法可以基于一张眼部或额头皱纹区域清晰无遮挡的人脸照片,自动对图像中人脸皮肤皱纹严重程度进行分级。自动分级结果的准确度经训练外数据的验证,可达0.85左右。
本发明采用的技术方案包括以下几个步骤:
1.输入清晰的人脸图像,分割人脸图像中的皱纹区域;
2.使用步骤1.1中分割的人脸皱纹区域图像进行皱纹识别,并计算皱纹纹理相关物理值(以下称纹理值);
3.在分割的人脸皱纹区域图像中,提取皱纹相关特征和其他非皱纹相关特征;
4.利用人脸皱纹区域图像的特征进行皱纹严重程度自动分级。
所述的步骤1具体包括:
1.1.使用dlib工具包(Python调用)进行人脸识别并检测68个面部特征点;
1.2.将面部特征点按面部轮廓、左侧眉毛、右侧眉毛、鼻子、左侧眼睛、右侧眼睛和嘴巴进行分类;
1.3.计算左侧眉毛特征点的平均水平间距与右侧眉毛特征点的平均水平间距的差值,根据图像的尺寸设置该差值的阈值,比较差值与阈值判断人脸侧向;
1.4.根据人脸侧向选择需要分割的人脸皱纹区域:
侧向朝前时可选择分割额头区域、眼袋区域;
侧向朝右时可选择分割额头区域、左侧眼角区域、左侧眼袋区域;
侧向朝左时可选择分割额头区域、右侧眼角区域、右侧眼袋区域;
1.5.结合五官特征点的位置、图像尺寸和人脸侧向,分割皱纹区域:
侧向朝前时,额头区域的大小随图像尺寸变化,中心位置根据眉毛特征点计算;眼袋区域大小随图像尺寸变化,中心位置根据眼睛特征点计算;
侧向朝右时,额头区域的大小随图像尺寸变化,比朝前时稍小,中心位置根据眉毛特征点计算;左侧眼角区域的大小随图像尺寸变化,中心位置根据眼睛特征点和眉毛特征点计算;
侧向朝左时,额头区域的大小随图像尺寸变化,比朝前时稍小,中心位置根据眉毛特征点计算;左侧眼角区域的大小随图像尺寸变化,中心位置根据眼睛特征点和眉毛特征点计算。
所述的步骤2具体包括:
2.1.首先将整张待分析皱纹图像归一化处理,利用毛发区域的局部均值相对较小的特点设置阈值剔除毛发的影响,再利用皮肤皱纹区域局部方差比平滑区域较大的特点,设置局部方差的阈值进行筛选,从而得到整幅图像皱纹区域的掩码,最后对皱纹区域归一化处理得到归一化的皱纹图像,后续处理都基于该归一化皱纹图像;
2.2.利用高斯滤波器构建x和y两个方向的梯度算子,使用梯度算子计算归一化图像在x和y方向上的微分,利用微分结果计算图像中每个像素点处的二倍角正余弦,再利用反正切函数计算每个像素处的纹理方向值,扩展到整幅图像上,则得到皱纹纹理方向图像,它反映了在每个像素点处的纹理方向;
2.3.对图像中每个相同大小的区块,利用区块的方向值旋转区块,使得区块纹理呈现竖直状态,使用求和的方式,将所有竖直方向上的信息投影到同一水平方向上,从而使整个区块上的纹理信息投影为一个类正弦波形,计算该类正弦波形的频率作为该区块的皱纹频率参考值,设置类正弦波形的最大波长阈值和最小波长阈值对频率参考值进行筛选作为最终的区块皱纹频率,扩展到整幅图像上每个区块的皱纹频率和整幅图像的皱纹频率平均值;
2.4.使用对局部皱纹纹理的频率和方向都有分析作用的Gabor filter对图像纹理进行分析。首先根据整幅图像的皱纹平均频率计算Gabor filter的尺寸,以该尺寸矩形的中心为原点,计算得到相位角为0°的Gabor filter,当对当前像素进行滤波时,根据该像素点的纹理方向计算相位角,将Gabor filter旋转相位角后,再对当前像素点对应的滤波区域进行滤波,以得到当前像素点的纹理相关物理值,扩展到整幅图像上即可得到皮肤皱纹纹理值图像;
2.5.对皮肤皱纹纹理值图像设置适当的阈值(纹理值<-2.5)进行分割得到初始的皱纹区域,计算初始皱纹区域的8连通区域面积,筛选面积最大的前15%的皱纹作为皱纹的识别结果。
所述的步骤3具体包括:
3.1.根据皮肤皱纹纹理值图像和识别的皱纹,提取皱纹深度、长度、面积和密度等皱纹相关特征:
1)根据皱纹识别的结果,统计每条皱纹或每个局部区域的相关特征作为皱纹特征分布;
2)根据获取的皱纹特征分布,计算单个分布的不同特征值;
3)根据获取的皱纹特征分布,计算皱纹特征分布之间的结合特征值;
3.2.根据原始皮肤皱纹图像,直接获取非皱纹相关特征:
1)计算原始皮肤皱纹图像的局部二值图像(Local Binary Pattern)的Uniform模式,统计LBP的概率密度直方图;
2)采用mahotas工具包(Python调用)计算原始皮肤皱纹图像的Haralick texturefeatures,它包含了灰度共生矩阵(GLCM)的13个相关特征值。
每张图像收集至少200个不同的皱纹特征。
所述的步骤4具体包括:
将提取的图像特征输入至多分类随机森林模型中,使用包外数据(Out of baggingdata)评估此预测模型的准确度。重复此过程30次,得到平均预测准确度为0.854,中位预测准确度为0.875。
本发明的有益效果:
本发明通过运用人脸识别、图像特征提取、机器学习的方法,提出了一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,自动分级结果的准确度为0.85左右;
本发明可对任意一张皱纹区域清晰无遮挡的图像进行分析,获取人脸皱纹的严重程度评判结果,可广泛应用于医美、皮肤衰老研究等领域,也可方便快捷地为普通人群提供合理的皮肤皱纹分级分析。
附图说明
图1是本发明的人脸侧向朝前示例图。
图2是本发明的人脸侧向朝前时,皱纹区域分割示例图。
图3是本发明的人脸侧向朝左示例图。
图4是本发明的人脸侧向朝左时,皱纹区域分割示例图。
图5是本发明的人脸侧向朝右示例图。
图6是本发明的人脸侧向朝右时,皱纹区域分割示例图。
图7是一张人脸眼袋区域皮肤皱纹图。
图8是本发明的方法针对该眼袋区域皮肤皱纹图获取的局部纹理方向示例图。
图9是本发明的方法针对该眼袋区域皮肤皱纹图获取的人脸皱纹纹理示例图。
图10是本发明的方法针对该眼袋区域皮肤皱纹图获取的人脸皱纹识别结果示例图。
图11是本发明的人工分级轻度皱纹示例图。
图12是本发明的人工分级中度皱纹示例图。
图13是本发明的人工分级重度皱纹示例图。
图14是本发明的设计流程图。
具体实施方式
人脸皮肤图像的捕获与接受条件:
本发明中,实验使用的人脸皮肤图像是使用VISIA皮肤检测仪拍照获取,但本发明提出的方法也同样适用于通过其他的数字摄像机或者光学仪器采集的人脸皮肤图像。人脸图像应当清晰,拍摄光照良好,眼袋、眼角和额头皱纹区域无阴影或头发的遮挡,此外图像不得经过滤镜或PS等处理。
人脸图像识别和人脸特征点分类:
本发明中,使用dlib工具包(Python调用)对人脸图像进行检测识别,同时获取包含68个人脸特征点的集合,按照以下方式对人脸特征点进行划分:
特征点集合中0~16特征点为人脸轮廓特征点;
特征点集合中17~21特征点为左侧眉毛特征点;
特征点集合中22~26特征点为右侧眉毛特征点;
特征点集合中27~35特征点为鼻子特征点;
特征点集合中36~41特征点为左侧眼睛特征点;
特征点集合中42~47特征点为右侧眼睛特征点;
特征点集合中48~67特征点为嘴巴特征点。
人脸侧向判断:
本发明中,计算左侧眉毛特征点的平均间距与右侧眉毛的平均间距的差值,同时根据图像的尺寸对差值设定适当的阈值(本实验图像尺寸为3456x5184,设置阈值50):
当差值大于设定的阈值时,人脸侧向朝右(示例图见图5);
当差值小于设定阈值的相反数时,人脸侧向朝左(示例图见图3);
当差值绝对值小于设定阈值时,人脸侧向朝前(示例图见图1)。
分割皱纹区域:
1.当人脸侧向朝右时,可选择分割额头区域、左侧眼角区域、左侧眼袋区域(示例图见图6),相同尺寸的不同图像同一类皱纹区域面积固定,由图像尺寸决定,分布位置由以下方式决定:
1.1.额头区域的中心由左侧眉毛特征点中点和右侧眉毛特征点中点决定:
X(forehead)=[X(left_brow_midp)+X(right_brow_leftp)]/2
Y(forehead)=min{Y(left_brow),Y(right_brow)}-A
其中X(forehead)和Y(forehead)分别是额头区域的中心x坐标和y坐标,X(left_brow_midp)是左侧眉毛特征点中点的x坐标,X(right_brow_leftp)是右侧眉毛特征点最左点的x坐标,Y(left_brow)是所有左侧眉毛特征点的y坐标,Y(right_brow)}是所有右侧眉毛特征点的y坐标,A是随图像尺寸变化的值;
1.2.左侧眼角区域的中心由左侧眼睛特征点的左眼角点和左侧眉毛特征点的左眉毛点决定:
X(left_canthus)=X(left_brow_leftp)-B
Y(left_canthus)=Y(left_eye_leftp)
其中X(left_canthus)和Y(left_canthus)分别是左侧眼角区域中心的x坐标和y坐标,X(left_brow_leftp)是左侧眉毛特征点的最左点的x坐标,B是随图像尺寸变化的值,Y(left_eye_leftp)是左侧眼睛特征点的最左点的y坐标;
1.3.左侧眼袋区域的中心由左侧眼睛特征点决定:
X(left_pouch)=mean{X(left_eye)}
Y(left_pouch)=max{Y(left_eye)}+C
其中X(left_pouch)和Y(left_pouch)分别是左侧眼袋区域中心的x坐标和y坐标,X(left_eye)和Y(left_eye)分别是所有左侧眼睛特征点的x坐标和y坐标,C是随图像尺寸变化的值。
2.当人脸侧向朝左时,可选择分割额头区域、右侧眼角区域、右侧眼袋区域(示例图见图4),相同尺寸的不同图像同一皱纹区域面积固定,由图像尺寸决定,分布位置由以下方式决定:
2.1.额头区域的中心由左侧眉毛特征点中点和右侧眉毛特征点中点决定:
X(forehead)=[X(right_brow_midp)-X(left_brow_rightp)]/2
Y(forehead)=min{Y(left_brow),Y(right_brow)}-A
其中X(forehead)和Y(forehead)分别是额头区域的中心x坐标和y坐标,X(right_brow_midp)是右侧眉毛特征点中点的x坐标,X(left_brow_rightp)是左侧眉毛特征点最右点的x坐标,Y(left_brow)是所有左侧眉毛特征点的y坐标,Y(right_brow)}是所有右侧眉毛特征点的y坐标,A是随图像尺寸变化的值;
2.2.右侧眼角区域的中心由右侧眼睛特征点的右眼角点和右侧眉毛特征点的右眉毛点决定:
X(right_canthus)=X(right_brow_rightp)+B
Y(right_canthus)=Y(right_eye_rightp)
其中X(right_canthus)和Y(right_canthus)分别是右侧眼角区域中心的x坐标和y坐标,X(right_brow_rightp)是右侧眉毛特征点的最右点的x坐标,Y(right_eye_rightp)是右侧眼睛特征点的最右点的y坐标,B是随图像尺寸变化的值;
2.3.右侧眼袋区域的中心由右侧眼睛特征点决定:
X(right_pouch)=mean{X(right_eye)}
Y(right_pouch)=max{Y(right_eye)}+C
其中X(right_pouch)和Y(right_pouch)分别是右侧眼袋区域中心的x坐标和y坐标,X(right_eye)和Y(right_eye)分别是所有右侧眼睛特征点的x坐标和y坐标,C是随图像尺寸变化的值。
3.当人脸侧向朝前时,可选择分割额头区域、左侧眼袋区域、右侧眼袋区域(示例图见图2),相同尺寸的不同图像同一皱纹区域面积固定,由图像尺寸决定,分布位置由以下方式决定:
3.1.额头区域的中心由左侧眉毛特征点中点和右侧眉毛特征点中点决定:
X(forehead)=[X(right_brow_midp)+X(left_brow_midp)]/2
Y(forehead)=min{Y(left_brow),Y(right_brow)}-A
其中X(forehead)和Y(forehead)分别是额头区域的中心x坐标和y坐标,X(left_brow_midp)是左侧眉毛特征点中点的x坐标,X(right_brow_midp)是右侧眉毛特征点中点的x坐标,Y(left_brow)是所有左侧眉毛特征点的y坐标,Y(right_brow)}是所有右侧眉毛特征点的y坐标,A是随图像尺寸变化的值;
3.2.左侧眼袋区域的中心由左侧眼睛特征点决定:
X(left_pouch)=mean{X(left_eye)}
Y(left_pouch)=max{Y(left_eye)}+C
其中X(left_pouch)和Y(left_pouch)分别是左侧眼袋区域中心的x坐标和y坐标,X(left_eye)和Y(left_eye)分别是所有左侧眼睛特征点的x坐标和y坐标,C是随图像尺寸变化的值;
3.3.右侧眼袋区域的中心由右侧眼睛特征点决定:
X(right_pouch)=mean{X(right_eye)}
Y(right_pouch)=max{Y(right_eye)}+C
其中X(right_pouch)和Y(right_pouch)分别是右侧眼袋区域中心的x坐标和y坐标,X(right_eye)和Y(right_eye)分别是所有右侧眼睛特征点的x坐标和y坐标,C是随图像尺寸变化的值。
归一化皱纹区域灰度值:
1.将图像由RGB三通道彩色图像转换为单通道灰度图像,使用整幅图像的平均值和标准差对图像进行归一化处理:
P(x,y)=[I(x,y)-mean(I)]/std(I)
其中I使输入图像的灰度图像,mean(I)是I的平均灰度值,std(I)是I的灰度值标准差,I(x,y)是I位于每个像素处的灰度值,P(x,y)是归一化图像中对应像素点的灰度值;
2.计算归一化处理后的图像的局部标准差和平均值,设置标准差阈值和平均值阈值,利用毛发区域局部平均值较小的特点剔除掉毛发区域,再利用皱纹区域标准差较平滑皮肤区域大的特点筛选出大致的皱纹区域;
3.针对大致的皱纹区域计算平均值与标准差,再次归一化处理皱纹区域,得到最终的皱纹区域归一化图像。
计算皱纹纹理方向(示例图见图8,示例图原图见图7):
1.计算高斯滤波器在水平和竖直两个方向上的梯度作为梯度算子,利用梯度算子与预处理结果图像进行卷积,获得图像中每个在水平和竖直两个方向上的偏微分;
2.使用高斯滤波器对偏微分进行平滑滤波,以减小偏微分矩阵的噪声;
3.根据偏微分矩阵计算出每个点处纹理方向的二倍角正余弦:
sin2θ=2dx(u,v)dy(u,v)/[dx2(u,v)+dy2(u,v)]
cos2θ=[dx2(u,v)+dy2(u,v)]/[dx2(u,v)+dy2(u,v)]
其中sin2θ和cos2θ皱纹纹理方向角的二倍角正弦和二倍角余弦,dx(u,v)和dy(u,v)是位于点(u,v)的x方向和y方向上的微分;
4.通过二倍角正余弦计算出纹理方向的二倍角正切值,再利用反正切函数计算得到图像上每个像素点的皱纹纹理方向角。
计算皱纹频率:
1.将图像分为若干同样大小的区块,计算区块上的皱纹纹理方向平均值作为该区块的皱纹纹理方向,根据皱纹纹理方向旋转区块使得皱纹纹理呈竖直状态;
2.皱纹的纹理在三维空间内呈现出类正弦平面波的形状,使用将同一竖直方向上的灰度值累加的方式将三维空间内的类正弦平面波投射到同一平面上,形成和皱纹纹理形状相关的类正弦波形;
3.计算类正弦波形的波形总长(第一个波峰或波谷和最后一个波峰或波谷之间的像素数)和数量(波峰或波谷数),设置波形数量阈值,过滤掉波形数量小于阈值的区块,计算保留区块上的皱纹纹理平均波长λ(λ=波形总长/数量);
4.利用阈值分割的方法对λ进行筛选:λmin≤λ≤λmax,以获取合理的皱纹纹理波长,当波长过大或过小时,均设置皱纹纹理频率为0;
5.利用筛选后的皱纹纹理波长计算区块图像的皱纹纹理频率,从而估计出整幅图像上所有局部区域的皱纹纹理频率和整幅图像的皱纹纹理平均频率。
计算皱纹纹理值(示例图见图9):
Gabor filter滤波器是用于分析纹理的滤波器,与纹理频率与方向相互关联,使用该滤波器对图像进行滤波可以获取得到皱纹纹理值,具体包括以下步骤:
1.使用皱纹纹理平均频率计算Gabor filter的尺寸,使得Gabor filter可以尽可能包含一条皱纹纹理在法线方向上的所有信息;
2.计算最初始相位角为0°的Gabor filter:
Gabor filter=exp[-(x2+y2)/σ2]cos(2πfx)
其中exp为自然对数的底数,x和y表示Gabor filter内的每个离散点的横纵坐标,σ是Gauss filter的标准差,为Gabor filter大小的1/6,i表示虚数,f表示整幅图像的皱纹平均频率;
3.将Gabor filter按特定相位角步长(5°)进行旋转,直到旋转相位角度达到180°,存储每个旋转后的Gabor filter以方便滤波时直接调用;
4.根据Gabor filter的大小分割出图像的滤波边界,对边界内的每个像素进行滤波,具体滤波步骤如下:
4.1.以当前像素为中心,获取一个和Gabor filter同大小的滤波区域;
4.2.根据当前像素的纹理方向值,将小于0°的纹理方向值增加180°,将大于180°的纹理值减去180°,使得所有纹理方向值都映射到(0°,180°)的范围内;
4.3.根据当前像素的纹理方向值,获取对应相位角的Gabor filter;
4.4.将当前像素对应的滤波区域与相应的Gabor filter相乘并求和,所获值即为当前像素点的纹理值,再遍历滤波边界内所有像素点,则可获取原始皱纹图像的纹理值图像。
识别皱纹(示例图见图10):
所获的纹理值图像,皱纹的皮脊越凸或皮沟越凹纹理值绝对值越大。采用阈值分割的方法,设定皱纹纹理值最大值阈值条件(纹理值<-2.5)分割出皱纹的皮沟区域,再对分割出的皮沟区域筛选其中面积最大的前15%条皮沟,这些皮沟即是最终识别的皱纹。
获取皱纹相关特征分布:
选取皱纹的面积、密度、深度和长度作为特征提取的对象进行分析:
1.通过计算每条识别皱纹的像素数量作为每条识别皱纹的面积,统计每条识别皱纹的面积作为皱纹面积分布;
2.将皱纹图像分为若干区块,计算区块内皱纹面积与区块面积的比例作为该区域内的皱纹密度,统计每个区块的皱纹密度作为皱纹密度分布,同时统计每个区块的皱纹纹理值分布相关特征值信息作为局部皱纹深度相关分布;
3.将每条皱纹的纹理值作为深度分布,计算该条皱纹的若干深度特征值,统计所有皱纹的同一深度特征值作为皱纹深度相关分布;
4.提取每条识别皱纹的骨架,骨架的像素数量作为每条识别皱纹的长度,统计每条识别皱纹的长度作为皱纹长度分布。
提取皱纹特征:
根据整幅图像的皱纹相关特征分布,计算相应特征值用以表征图像:
1.根据皱纹相关特征分布,计算每个特征分布的相关特征值:
1.1.计算分布的最小值、最大值、平均值、10分位点、25分位点、50分位点、75分位点、90分位点等;
1.2.计算子分布的极差,如从最小值到最大值的分布极差,从10分位点到90分位点的子分布极差,从25分位点到75分位点的子分布极差;
1.3.计算子分布的平均绝对离差、标准差,如从最小值到最大值的分布平均绝对离差和标准差,从10分位点到90分位点的平均绝对离差和标准差;
1.4.计算分布的平方和及相关伸缩值;
1.5.计算分布的熵;
1.6.计算分布的统计直方图,计算落入每个直方图的概率密度平方和;
1.7.计算分布的三阶矩(偏度)、四阶矩(峰度);
2.根据皱纹相关特征分布,计算皱纹特征分布的结合特征值:
2.1.计算皱纹长度分布与皱纹深度特征值(如最大深度、平均深度等)分布之间的点积;
2.2.计算皱纹密度分布的密度特征值(如最大密度、50分位点密度)与该区域对应深度特征值(如最大深度、平均深度等)乘积;
2.3.计算皱纹密度分布的密度特征值分布与对应区域的深度特征值(如最大深度、平均深度等)分布的点积。
提取原始图像特征:
根据原始图像直接提取特征,不参考皱纹识别的结果。具体包括以下方案:
1.计算原始皱纹图像的局部二值图像(Local Binary Pattern)的Uniform模式,计算局部二值图像的统计直方图,用以描述图像局部模式和像素排列规则;
2.使用mahotas工具包(python调用)计算原始皱纹图像的Haralick texturefeatures,它基于图像的灰度共生矩阵(GLCM),用以描述图像灰度的方向、间隔、变化幅度等信息特征。
人脸皱纹严重程度人工分级:
人脸皱纹严重程度的人工分级分为三种程度,每种程度的分级标准具体包括以下:
1.轻度(Mild,示例图见图11),满足如下任一标准:
1.1.皱纹完全不可见;
1.2.皱纹隐约可见;
1.3.可见皱纹非常浅,且可见皱纹只覆盖了局部区域(覆盖区域面积<50%);
2.中度(Moderate,示例图见图12),满足如下任一标准:
2.1.部分皱纹可见且深度较浅;
2.2.皱纹较浅,可见皱纹覆盖面积较广(覆盖区域面积>50%);
3.重度(Severe,示例图见图13),满足如下任一标准:
3.1.部分皱纹明显可见且深度较深;
3.2.皱纹较深,可见皱纹覆盖面积较广(覆盖区域面积>50%)。
皱纹严重程度分级:
使用多分类的随机森林方法,结合上述提取的图像特征和人工分级结果,训练所得的随机森林模型用于皱纹严重程度的自动分级。

Claims (5)

1.一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,其特征在于:
1.1.输入清晰的人脸图像,分割人脸图像中的皱纹区域;
1.2.使用步骤1.1中分割的人脸皱纹区域图像进行皱纹识别,并计算皱纹纹理相关物理值(以下称纹理值);
1.3.在分割的人脸皱纹区域图像中,提取皱纹相关特征和非皱纹相关特征;
1.4.利用人脸皱纹区域图像的特征进行皱纹严重程度自动分级。
2.根据权利要求1所述的一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,其特征在于:步骤1.1具体包括:
1)对输入的人脸图像进行人脸识别,并检测五官特征点;
2)根据五官特征点的位置和图像的尺寸,判断人脸图像中人脸的侧向;
3)根据人脸的侧向选择分割人脸皱纹区域;
4)结合五官特征点的位置、图像的尺寸和人脸的侧向,计算分割人脸皱纹区域的大小与位置。
3.根据权利要求1所述的一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,其特征在于:步骤1.2具体包括:
1)使用图像的均值和标准差归一化图像,设置局部均值阈值剔除毛发区域,设置局部标准差阈值得到大致的皱纹区域,最后针对皱纹区域单独计算均值和标准差进行归一化,得到归一化皱纹图像;
2)构造梯度算子,计算归一化皱纹图像的皱纹纹理在水平和竖直两个方向上的微分,再根据微分计算得到皱纹纹理方向角的二倍角正余弦,最后利用反正切函数得到皱纹区域的皱纹纹理方向;
3)利用局部皱纹纹理方向,旋转归一化皱纹图像的局部皱纹纹理使其竖直,将该局部区域的像素值沿竖直方向叠加,即将三维空间内的皱纹纹理平面波变换到二维空间内,形成一个类正弦波形,计算正弦波形的平均波长和波形的个数,并设置平均波长的阈值和波形个数的阈值,筛选出适合的局部纹理平均波长,最后再求出局部皱纹纹理频率,扩展到整幅图像即求得每个区域的皱纹纹理频率;
4)利用皱纹纹理平均构造相位角为0°的皱纹纹理纹理分析滤波器Gabor filter,当对归一化皱纹图像当前像素点进行滤波时,根据滤波器的大小获取滤波区域,利用当前像素点的皱纹纹理方向角旋转滤波器,计算当前点的纹理值,扩散到整幅图像中则可获取纹理值图像;
5)对纹理值图像使用阈值分割的方法分割出皱纹皮沟,作为最终识别的皱纹区域。
4.根据权利要求1所述的一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,其特征在于:步骤1.3具体包括:
1)根据皱纹识别的结果获取皱纹区域的掩码,结合纹理值图像,获取皱纹特征的分布;
2)根据皱纹特征的分布,计算单个分布的特征值和分布之间的结合特征值;
3)直接根据步骤1.中分割的人脸皱纹图像,提取图像特征。
5.根据权利要求1.1所述的一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法,其特征在于:步骤1.4具体包括:
以提取的图像特征作为输入值,使用随机森林方法自动预测该图像的分类,即将该图像中皱纹的严重程度分级。
CN201910352167.7A 2019-04-25 2019-04-25 一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法 Pending CN110110637A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910352167.7A CN110110637A (zh) 2019-04-25 2019-04-25 一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910352167.7A CN110110637A (zh) 2019-04-25 2019-04-25 一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法

Publications (1)

Publication Number Publication Date
CN110110637A true CN110110637A (zh) 2019-08-09

Family

ID=67487413

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910352167.7A Pending CN110110637A (zh) 2019-04-25 2019-04-25 一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法

Country Status (1)

Country Link
CN (1) CN110110637A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110738678A (zh) * 2019-10-18 2020-01-31 厦门美图之家科技有限公司 脸部细纹检测方法、装置、电子设备及可读存储介质
CN110956623A (zh) * 2019-11-29 2020-04-03 深圳和而泰家居在线网络科技有限公司 皱纹检测方法、装置、设备及计算机可读存储介质
CN111767846A (zh) * 2020-06-29 2020-10-13 北京百度网讯科技有限公司 图像识别方法、装置、设备和计算机存储介质
CN111899271A (zh) * 2019-12-16 2020-11-06 西北工业大学 一种图像分割和超声波探伤仪的自动检定的方法及系统
CN112613459A (zh) * 2020-12-30 2021-04-06 深圳艾摩米智能科技有限公司 一种人脸敏感区域检测方法
CN112712054A (zh) * 2021-01-14 2021-04-27 深圳艾摩米智能科技有限公司 脸部皱纹检测方法
CN113499034A (zh) * 2021-06-29 2021-10-15 普希斯(广州)科技股份有限公司 一种皮肤检测方法、系统及美容装置
CN113907717A (zh) * 2021-11-01 2022-01-11 南京工程学院 基于妊娠纹严重性客观评价指标的妊娠纹严重性评价方法
CN114612994A (zh) * 2022-03-23 2022-06-10 深圳伯德睿捷健康科技有限公司 训练皱纹检测模型的方法、装置及检测皱纹的方法、装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006107288A (ja) * 2004-10-07 2006-04-20 Toshiba Corp 個人認証方法、装置及びプログラム
CN104732200A (zh) * 2015-01-28 2015-06-24 广州远信网络科技发展有限公司 一种皮肤类型和皮肤问题的识别方法
US20150356344A1 (en) * 2014-06-09 2015-12-10 Panasonic Intellectual Property Management Co., Ltd. Wrinkle detection apparatus and wrinkle detection method
US20170076146A1 (en) * 2015-09-11 2017-03-16 EyeVerify Inc. Fusing ocular-vascular with facial and/or sub-facial information for biometric systems
CN107392866A (zh) * 2017-07-07 2017-11-24 武汉科技大学 一种光照鲁棒的人脸图像局部纹理增强方法
CN108369644A (zh) * 2017-07-17 2018-08-03 深圳和而泰智能控制股份有限公司 一种定量检测人脸抬头纹的方法和智能终端
CN108932493A (zh) * 2018-06-29 2018-12-04 东北大学 一种面部皮肤质量评价方法
CN109086688A (zh) * 2018-07-13 2018-12-25 北京科莱普云技术有限公司 脸部皱纹检测方法、装置、计算机设备和存储介质

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006107288A (ja) * 2004-10-07 2006-04-20 Toshiba Corp 個人認証方法、装置及びプログラム
US20150356344A1 (en) * 2014-06-09 2015-12-10 Panasonic Intellectual Property Management Co., Ltd. Wrinkle detection apparatus and wrinkle detection method
CN104732200A (zh) * 2015-01-28 2015-06-24 广州远信网络科技发展有限公司 一种皮肤类型和皮肤问题的识别方法
US20170076146A1 (en) * 2015-09-11 2017-03-16 EyeVerify Inc. Fusing ocular-vascular with facial and/or sub-facial information for biometric systems
CN107392866A (zh) * 2017-07-07 2017-11-24 武汉科技大学 一种光照鲁棒的人脸图像局部纹理增强方法
CN108369644A (zh) * 2017-07-17 2018-08-03 深圳和而泰智能控制股份有限公司 一种定量检测人脸抬头纹的方法和智能终端
WO2019014814A1 (zh) * 2017-07-17 2019-01-24 深圳和而泰智能控制股份有限公司 一种定量检测人脸抬头纹的方法和智能终端
CN108932493A (zh) * 2018-06-29 2018-12-04 东北大学 一种面部皮肤质量评价方法
CN109086688A (zh) * 2018-07-13 2018-12-25 北京科莱普云技术有限公司 脸部皱纹检测方法、装置、计算机设备和存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
顼改燕等: "基于Gabor滤波器和BP神经网络的人脸皮肤皱纹区域自动识别", 《计算机应用》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110738678A (zh) * 2019-10-18 2020-01-31 厦门美图之家科技有限公司 脸部细纹检测方法、装置、电子设备及可读存储介质
CN110738678B (zh) * 2019-10-18 2022-05-31 厦门美图宜肤科技有限公司 脸部细纹检测方法、装置、电子设备及可读存储介质
CN110956623A (zh) * 2019-11-29 2020-04-03 深圳和而泰家居在线网络科技有限公司 皱纹检测方法、装置、设备及计算机可读存储介质
CN110956623B (zh) * 2019-11-29 2023-11-07 深圳数联天下智能科技有限公司 皱纹检测方法、装置、设备及计算机可读存储介质
CN111899271A (zh) * 2019-12-16 2020-11-06 西北工业大学 一种图像分割和超声波探伤仪的自动检定的方法及系统
CN111767846A (zh) * 2020-06-29 2020-10-13 北京百度网讯科技有限公司 图像识别方法、装置、设备和计算机存储介质
CN112613459A (zh) * 2020-12-30 2021-04-06 深圳艾摩米智能科技有限公司 一种人脸敏感区域检测方法
CN112712054A (zh) * 2021-01-14 2021-04-27 深圳艾摩米智能科技有限公司 脸部皱纹检测方法
CN112712054B (zh) * 2021-01-14 2024-06-18 深圳艾摩米智能科技有限公司 脸部皱纹检测方法
CN113499034A (zh) * 2021-06-29 2021-10-15 普希斯(广州)科技股份有限公司 一种皮肤检测方法、系统及美容装置
CN113907717A (zh) * 2021-11-01 2022-01-11 南京工程学院 基于妊娠纹严重性客观评价指标的妊娠纹严重性评价方法
CN114612994A (zh) * 2022-03-23 2022-06-10 深圳伯德睿捷健康科技有限公司 训练皱纹检测模型的方法、装置及检测皱纹的方法、装置

Similar Documents

Publication Publication Date Title
CN110110637A (zh) 一种人脸皮肤皱纹自动识别和皱纹严重程度自动分级的方法
Barata et al. A system for the detection of pigment network in dermoscopy images using directional filters
US6151403A (en) Method for automatic detection of human eyes in digital images
US7123783B2 (en) Face classification using curvature-based multi-scale morphology
US7953253B2 (en) Face detection on mobile devices
CN112396573A (zh) 基于图像识别的面部皮肤分析方法及系统
Naji et al. Skin segmentation based on multi pixel color clustering models
CN111524080A (zh) 脸部皮肤特征的识别方法、终端及计算机设备
Abate et al. BIRD: Watershed based iris detection for mobile devices
Martinez et al. Facial component detection in thermal imagery
Monwar et al. Pain recognition using artificial neural network
CN111460950A (zh) 自然驾驶通话行为中基于头-眼证据融合的认知分心方法
Buse et al. A structural and relational approach to handwritten word recognition
Jamil et al. Melanocytic and nevus lesion detection from diseased dermoscopic images using fuzzy and wavelet techniques
Graf et al. Robust recognition of faces and facial features with a multi-modal system
Monwar et al. Eigenimage based pain expression recognition
Takruri et al. Automatic recognition of melanoma using Support Vector Machines: A study based on Wavelet, Curvelet and color features
CN110298815B (zh) 一种皮肤毛孔检测与评价的方法
Karamizadeh et al. Race classification using gaussian-based weight K-nn algorithm for face recognition
Chen et al. Contour detection by simulating the curvature cell in the visual cortex and its application to object classification
KR100596197B1 (ko) 가변 타원 마스크와 형태학적 특징을 이용한 얼굴 인식 방법
Gizatdinova et al. Automatic edge-based localization of facial features from images with complex facial expressions
Chou et al. Toward face detection, pose estimation and human recognition from hyperspectral imagery
Monwar et al. A real-time face recognition approach from video sequence using skin color model and eigenface method
Colak et al. Eyematch: an eye localization method for frontal face images

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190809

WD01 Invention patent application deemed withdrawn after publication