一种三维人脸相似性度量方法
【专利摘要】本发明涉及一种三维人脸相似性度量方法，所述方法包括：在两个待比较的三维人脸模型上，分别以各自的鼻尖点为源点，按相同的初始方向和相等的角度间隔均匀地提取m条测地线；分别计算两个三维人脸模型上方向相同的两条测地线的Fréchet距离；计算m个Fréchet距离值的平均值，得到两个三维人脸模型的Fréchet距离；计算两个三维人脸模型的相似度。应用本发明计算出的同一个人的两个不同的三维人脸模型的相似度值高于不同人的相似度值，而且不同人的相似度值与人的主观判断基本一致。本发明能够有效区分同一个人不同表情的人脸模型，给出不同人的三维模型之间量化的相似度值，能够用于表情变化的三维人脸相似性度量。
【专利说明】
一种三维人脸相似性度量方法
技术领域
[0001]本发明属于人脸识别技术领域，具体涉及一种基于测地线的Frgchet(弗雷歇)距 离的三维人脸相似性度量方法。
【背景技术】
[0002] 三维人脸相似性度量属于计算机图形学、计算机视觉等领域，是三维人脸识别的 基础，在公安刑侦、医学整容、人类学研究等领域都有许多重要的应用。虽然人类的视觉系 统可以在大脑中快速对两个人脸的相似程度做出判断，但要让计算机自动地判断出两个人 脸的相似程度还是比较困难的，因为一方面人脸的整体形状大致上是相同的（如都有鼻子、 眼睛、嘴等），只能通过细节来区分不同的人脸;另一方面人脸相似性度量还要解决同一个 人在不同表情、不同姿态等情况下的辨别，这更加大了三维人脸相似性度量的难度。
[0003] 目前，关于人脸相似性的研究还比较少，早期多采用主观定性的方法，即集合一定 量的被试对一组面貌的相似程度做出主观判断。这种方法虽然符合人的认知理论，但评价 过程费时费力，而且评价结果容易受到人的主观因素的影响。近年来，人们开始采用客观方 法对三维人脸的相似性进行研究，如采用在深度图像和三维人脸纹理图像上提取Gabor特 征进行相似度评价，并同主观评价结果对比的方法;采用相对角直方图、半径相对角直方 图、Procrustes距离及Principal warps对三维颜面的相似性进行研究的方法等。这些方法 主要是针对没有表情变化的三维颅面模型或深度图像进行研究，均没有考虑到表情变化对 三维人脸相似性的影响，人脸相似性度量的精度低。
[0004] 中国发明专利申请号201510008187.4公开了一种基于半刚性区域面部轮廓线的 三维人脸识别方法，其特征在于该方法包括以下步骤:步骤1:从输入的三维人脸点云中提 取出感兴趣的人脸区域，该区域为以鼻尖为球心，90mm为半径的球所包含的人脸区域;步骤 2:对由步骤1所提取的三维人脸点云进行网格化操作，采用基于网格的平滑算法对三维人 脸模型进行平滑去噪处理，然后将经过9~12次迭代处理得到的平滑的三维人脸网格恢复 成三维人脸点云;步骤3:对由步骤2得到的三维人脸点云，采用主成分分析法PCA，得到三个 互相垂直的主轴方向，将最大特征值对应的特征向量作为轴，最小特征值对应的特征向量 作为轴，建立右手坐标系，在该坐标系下人脸具有相同的正面姿态;该坐标系被称为姿态坐 标系PCS;以鼻尖点为PCS的坐标原点，将切割平滑后的人脸点云数据都转换到PCS中，从而 实现人脸姿态的归一化，得到姿态矫正后的人脸点云;步骤4:对由步骤3得到的姿态矫正后 的人脸点云，根据坐标信息，在人脸的半刚性区域提取垂直方向的若干条面部轮廓线并对 面部轮廓线进行重采样，得到测试人脸的19条面部轮廓线；步骤5:对所有库集人脸模型按 照步骤1-步骤4所示方法来提取每一个库集人脸的19条面部轮廓线，计算测试人脸与库集 人脸对应的轮廓线间的相似度向量;步骤6:采用决策级融合中的加权求和规则对不同轮廓 线的相似度向量进行融合，根据最近邻法则得出识别结果。上述三维人脸识别方法任然没 有考虑到人脸的情绪变化对于人脸识别精确度的影响。
【发明内容】
[0005] 为了解决现有技术中存在的上述问题，本发明从全局的角度提供一种基于测地线 的Fr6chet距离的三维人脸相似性度量方法。
[0006] 本发明通过两个三维人脸模型上对应测地线的相似度计算所述两个三维人脸模 型的相似度。每对测地线的相似度采用测地线的Frgchet距离表示:两条测地线的Frgchet 距离越小，两条测地线越相似;反之，越不相似。
[0007] 一种三维人脸相似性度量方法，包括以下步骤：
[0008] 步骤1，在两个待比较的三维人脸模型上，分别以各自的鼻尖点为源点，按相同的 初始方向和相等的角度间隔均匀地提取m条测地线。
[0009] 步骤2,分别计算两个三维人脸模型上方向相同的两条测地线的Frgchet距离，得 至丨Jm个Fr6chet距离值。
[0010] 步骤3,计算m个Frgchet距离值的平均值，得到两个三维人脸模型的Frgchet距离 d〇
[0011] 步骤4,计算两个三维人脸模型的相似度S，公式如下：
[0012] S = l-d/dmax
[0013]式中，dmaX为对多人的三维人脸模型进行计算得到的d中的最大值。
[0014] 进一步地，步骤1所述提取m条测地线的方法如下：
[0015] 步骤1.1，求三维人脸模型中纵坐标最大的点即脸部最高的点，该点即为鼻尖点；
[0016] 步骤1.2，求解最外层等测地线；
[0017] 步骤1.3，在所述最外层等测地线上求m个分点；
[0018] 步骤1.4,以所述鼻尖点为源点，分别以所述m个分点为目标点求解测地线，得到m 条测地线。
[0019] 更进一步地，步骤1.2求解所述最外层等测地线的方法如下：
[0020] 步骤1.2.1，提取三维人脸模型的边界；
[0021 ]步骤1.2.2，计算从鼻尖点到所述边界上所有点的测地距离；
[0022]步骤1.2.3,求所述测地距离的最小值，提取所有到鼻尖点的测地距离等于这该最 小值的点得到最外层等测地线。
[0023] 更进一步地，步骤1.3在所述最外层等测地线上求m个分点的方法如下：
[0024] 步骤1.3.1，以所述最外层等测地线上与所述鼻尖点横坐标相同的点为第一个分 占.
[0025] 步骤1.3.2，将所述最外层等测地线投影到所述鼻尖点的切平面得到曲线1，以所 述鼻尖点为起点，作m条角度间隔为231/m的射线，与曲线1相交得到 m个交点，其中第一个交 点为第一个分点的投影；
[0026] 步骤1.3.3，对所述m个交点进行反投影变换，在所述最外层等测地线上得到m个分 点。
[0027] 优选地，所述步骤1 . 2 . 2采用MMP算法（MMP算法是由Mitchell、Mount和 Papadimitriou三人于 1987年在"The discrete geodesic problem"一文中提出的测地线 算法，简称为MMP算法)求解所述测地距离。
[0028]优选地，所述步骤1.4采用MMP算法求解测地线。
[0029] 优选地，采用离散变体的Fr6chet距离算法求解所述Fr6chet距离。
[0030] 与现有技术相比，本发明具有以下有益效果：
[0031] 本发明通过在三维人脸模型上自鼻尖出发以等角度间隔均匀地向四周提取若干 条测地线，将三维曲面的比较转化为三维曲线的比较，以度量曲线相似性的Frgchet距离度 量两个三维人脸模型上对应曲线的相似程度。在Gavadb三维人脸公开数据库、Texas 3D人 脸识别库及变形的三维人脸数据和三维真实人脸数据库上进行实验的结果表明，应用本发 明所述方法计算出的同一个人的两个不同的三维人脸模型的相似度值，高于不同人的三维 人脸模型的相似度值，而且不同人的三维人脸模型的相似度值与人的主观判断基本一致。 因此，本发明不仅能够有效区分同一个人不同表情的人脸模型，而且能够给出不同人的三 维模型之间量化的相似度值，能够用于表情变化的三维人脸相似性度量。
【附图说明】
[0032]图1为本发明所述方法的流程图；
[0033]图2为本发明所述方法求解原理示意图。
[0034]图中：1一目标点，2-第一目标点，3-二维人脸模型，4一鼻尖点，5-测地线，6- 最外层等测地线。
【具体实施方式】
[0035]下面结合附图和实施例对本发明做进一步说明。
[0036] 一种三维人脸相似性度量方法，其流程图如图1所示，包括以下步骤：
[0037] 步骤S100,提取测地线。
[0038]在两个待比较的三维人脸模型心和5上，分别以各自的鼻尖点为源点，如图2所示， 按相同的初始方向和相等的角度间隔均匀地提取m条测地线，即：
[0042]式中，紀、gf分别为在fjPf2上提取的方向角为α的测地线，Gi、G2分别为由m条测 地线d、m条测地线gf组成的集合，m的值根据精度要求选取，A为角度的集合。
[0043] 步骤S200,计算两个三维人脸模型上对应测地线的Fr6chet距离，得到m个Fr6chet 距离值3dF。
[0044] 步骤S300,计算两个三维人脸模型的Fr6chet距离d(fi，f2)。
[0045] 两个三维人脸模型的Frgchet距离定义为其上所有对应测地线的Frgchet距离的 平均值，公式如下：
[0047]步骤S400,计算两个三维人脸模型的相似度。
[0048] 两个三维人脸模型f#Pf2的相似度通过这两个模型上对应测地线的相似性来计 算，每对测地线的相似度用两个模型fdPf 2上角度均为α的对应测地线的Frgchet距离表示。 两条曲线的Frgchet距离越小，说明这两条曲线越相似。Frgchet距离越大，说明两条曲线越 不相似。因此，两个三维人脸模型心和5的相似度能够用两个三维人脸模型的Frgchet距离d (f^h)来表示:Fr6Chet距离越小，这两个模型越相似;相反，距离越大，说明越不 相似。计算公式如下：
[0049] S(fl，f2) = l-d(fl，f2)/dmax
[0050] 式中，dmax为对多人的三维人脸模型行计算得到的d ( f 1，f 2 )中的最大值。很明显，S (fl，f2) ￡[0,1]。8 = 0表示两个三维人脸模型完全不相似，S = 1表示完全相同。S值越大表示 越相似。上述定义的相似性函数满足自反性、对称性和伪三角不等式性质，即：
[0051] 8(乜，；^) = 1(自反性）
[0052] = s(fj，fi)(对称性）
[0054] 其中，i、j、k的值均为1或2。
[0055] 所述步骤S100提取测地线的方法包括以下步骤：
[0056] S110,求鼻尖点；
[0057] 处于标准姿态的三维人脸模型，鼻尖点是整个人脸中最高的点，因此标准姿态下 整个三维人脸模型中纵坐标最大的点即为鼻尖点，并将其作为测地线的源点；
[0058] S120,求最外层等测地线；
[0059] S130,在最外层等测地线上求m个分点；
[0060] S140，以鼻尖点为源点，分别以所述m个分点为目标点求解测地线，得到m条测地 线。
[0061] 所述步骤S120求最外层等测地线的方法包括以下步骤：
[0062] S121，提取三维人脸模型的边界；
[0063] S122,计算从鼻尖点到所述边界上所有点的测地距离；
[0064] S123,求所述测地距离的最小值，提取所有到鼻尖点的测地距离等于这该最小值 的点得到最外层等测地线。
[0065]所述步骤S130在最外层等测地线上求m个分点的方法包括以下步骤：
[0066] S131，以所述最外层等测地线上与所述鼻尖点横坐标相同的点为第一个分点；
[0067] S132,将所述最外层等测地线投影到所述鼻尖点的切平面得到曲线1，以所述鼻尖 点为起点，作m条角度间隔为2Vm的射线，与曲线1相交得到m个交点，其中第一个交点为第 一个分点的投影；
[0068] S133,对所述m个交点进行反投影变换，在所述最外层等测地线上得到m个分点。
[0069]作为一种最佳实施例，步骤S122采用MMP算法求解所述测地距离。
[0070]作为一种最佳实施例，步骤S140采用MMP算法求解测地线。
[0071 ]作为一种最佳实施例，采用离散变体的Fr6chet距离算法求解所述Fr6chet距离。 [0072]本发明不限于上述实施方式，本领域技术人员所做出的对上述实施方式任何显而 易见的改进或变更，都不会超出本发明的构思和所附权利要求的保护范围。
【主权项】
1. 一种三维人脸相似性度量方法，其特征在于，包括以下步骤： 步骤1，在两个待比较的三维人脸模型上，分别以各自的鼻尖点为源点，按相同的初始 方向和相等的角度间隔均匀地提取m条测地线;m的值根据精度要求选取； 步骤2，分别计算两个三维人脸模型上方向相同的两条测地线的Frgchet距离，得到m个 Fr6chet距离值； 步骤3,计算m个Frgchet距离值的平均值，得到两个三维人脸模型的Frgchet距离d; 步骤4,计算两个三维人脸模型的相似度s，公式如下： S - 1-d/dmax 式中，dmax为对多人的三维人脸模型进行计算得到的d中的最大值。2. 根据权利要求1所述的三维人脸相似性度量方法，其特征在于，步骤1所述提取m条测 地线的方法如下： 步骤1.1，求三维人脸模型中纵坐标最大的点即脸部最高的点，该点即为鼻尖点； 步骤1.2，求解最外层等测地线； 步骤1.3，在所述最外层等测地线上求m个分点； 步骤1.4,以所述鼻尖点为源点，分别以所述m个分点为目标点求解测地线，得到m条测 地线。3. 根据权利要求2所述的三维人脸相似性度量方法，其特征在于，步骤1.2求解所述最 外层等测地线的方法如下： 步骤1.2.1，提取三维人脸模型的边界； 步骤1.2.2，计算从鼻尖点到所述边界上所有点的测地距离； 步骤1.2.3，求所述测地距离的最小值，提取所有到鼻尖点的测地距离等于这该最小值 的点得到最外层等测地线。4. 根据权利要求2所述的三维人脸相似性度量方法，其特征在于，步骤1.3在所述最外 层等测地线上求m个分点的方法如下： 步骤1.3.1，以所述最外层等测地线上与所述鼻尖点横坐标相同的点为第一个分点； 步骤1.3.2,将所述最外层等测地线投影到所述鼻尖点的切平面得到曲线1，以所述鼻 尖点为起点，作m条角度间隔为2 π /m的射线，与曲线1相交得到m个交点，其中第一个交点为 第一个分点的投影； 步骤1.3.3，对所述m个交点进行反投影变换，在所述最外层等测地线上得到m个分点。5. 根据权利要求3所述的三维人脸相似性度量方法，其特征在于，所述步骤1.2.2采用 MMP算法求解所述测地距离。6. 根据权利要求2所述的三维人脸相似性度量方法，其特征在于，所述步骤1.4采用MMP 算法求解测地线。7. 根据权利要求1所述的三维人脸相似性度量方法，其特征在于，采用离散变体的Fr6 chet距离算法求解所述Fr6chet距离。
【文档编号】G06T15/10GK105868727SQ201610219246
【公开日】2016年8月17日
【申请日】2016年4月8日
【发明人】赵俊莉, 武仲科, 段福庆, 王康, 陈郁葱, 刘翠婷
【申请人】青岛大学