专利名称：一种基于张量空间偏最小二乘回归的颅面重构方法
技术领域：
本发明涉及面皮重构方法，特别涉及一种通过在张量空间进行偏最小二乘回归来实现对未知颅骨进行面皮重构的方法。主要用于刑侦、考古、法医人类学、医学手术、美容整容等领域。
背景技术：
颅骨是人类面貌的内在生物特征，根据颅骨形态实现面貌重构是法医人类学领域的典型应用项目。传统的手工重构是在待定颅骨的基础上，借助解剖测量和针刺测量途径获得人脸软组织的厚度规律，由艺术家雕塑而成，属于人类学家、艺术家和医生的领域，较多的依赖重构者个人对颅面特征和艺术的把握，其结果受主观因素影响较大，不具备科学推广的价值，并且重构一个颅面需要几天的时间。目前，计算机辅助的颅面三维重构已引起国际信息学、人类学、法医学等相关领域的高度重视，成为国际上研究的热点。典型技术可分为两类一类是基于软组织厚度测量的方法。通过测量一定数量的颅面样本在颅骨上一些关键特征点处的平均软组织厚度作为经验知识，对于目标颅骨，首先根据这些平均软组织厚度的经验值，获得颅骨关键特征点处的面皮位置，然后通过插值、 变形等技术得到整个颅骨的面皮。这类方法只要求在颅骨上标注少量关键特征点，并获得这些点的软组织厚度，所需的数据较少。然而，这些关键特征点处的软组织厚度和颅面形状之间并没有直接的相关关系；另一类是基于统计变形模型的方法。这类方法将一套颅骨和面皮数据组织为一个整体形成一个向量，通过对大量的训练样本进行主成份分析(PCA)，建立参数化的颅面统计变形模型，通过优化的方法将统计模型匹配到目标颅骨，获得模型参数实现颅面重构。这类方法本质上是对仅具有小部分数据的颅面样本求解数据缺失这一病态问题，优化过程容易陷入局部极小，很难体现颅骨对面皮的决定性作用。
发明内容
本发明的目的在于提供一种基于张量空间偏最小二乘回归的颅面重构方法。为实现上述目的，本发明的技术解决方案是分别将颅骨和面皮变换到不同的张量特征空间，利用偏最小二乘技术分性别进行面皮张量特征对颅骨张量特征和颅面属性的回归建模。对待重构的目标颅骨，利用建立的回归模型得到相应的面皮张量特征，根据张量特征重构出目标面皮。方法简单易行，速度快，充分体现了颅骨决定面貌这一本质联系，重构准确度高。本发明所述颅面重构方法的主要步骤包括1模型训练1. 1数据处理，把训练样本中的三维颅面统一到法兰克福坐标系(Frankfurt Horizontal)并将大小规格化，去除数据中不属于面部的部分，对颅骨和面皮分别进行三维数据配准，将配准后的颅骨和面皮转化为二维深度图。1. 2张量特征提取，将颅骨和面皮分别变换到颅骨张量特征空间和面皮张量特征空间；1. 3回归模型建立，采用偏最小二乘回归分性别建立面皮张量特征对颅骨张量特征和属性(年龄、体重指数BMI)的回归模型。2颅面重构2. 1数据处理把待重构的未知颅骨统一到法兰克福坐标系并将大小规格化，将未知颅骨与训练样本中的颅骨进行数据配准，将配准后的未知颅骨转化为二维深度图；2. 2张量特征提取将未知颅骨变换到颅骨张量特征空间；2. 3特征回归，估计颅骨的性别和属性(年龄、体重指数BMI)，根据其张量特征和属性，利用1.3)中建立的回归模型回归出相应的面皮张量特征；2. 4三维面皮重构，根据回归到的面皮张量特征恢复其二维深度图，将得到的二维深度图转化为三维点云数据。
图1为本发明基于张量空间偏最小二乘回归的颅面重构方法的总体流程图。
具体实施例方式下面结合附图对本发明所述一种基于张量空间偏最小二乘回归的颅面重构法作进一步详细的说明。下面参考图1对本发明基于张量空间偏最小二乘回归的颅面重构方法进行详尽的描述。实现本发明所述方法的具体步骤是1模型训练1. 1数据处理把三维颅面统一到法兰克福坐标系，按颅骨左耳孔Lp和右耳孔Rp这两个特征点间的距离将颅骨和面皮的大小规格化，人工去除数据中不属于面部的部分，对颅骨训练样本和面皮训练样本分别进行三维数据配准，将配准后的颅骨和面皮转化为二维深度图。1. 2张量特征提取按照颅面的生理特征将每个颅骨和每个面皮的深度图分为k个大小相同且部分重叠的小区域，分别将训练样本中全体颅骨和全体面皮的深度图按身份、区域位置、区域横向和纵向信息变化的方向组织为四阶张量Ds e Rnxkxhx^P Df e Rmxkxhxw,其中，η和m分别为训练样本中颅骨和面皮的个数，k为小区域个数，h和w分别是每个小区域的行数和列数， 通过执行张量分解和空间降维得
权利要求
1. 一种基于张量空间偏最小二乘回归的颅面重构方法，其特征在于，包括下列步骤1)模型训练1.1)数据处理，把训练样本中的三维颅面统一到法兰克福坐标系并将大小规格化，去除数据中不属于面部的部分，对颅骨和面皮分别进行三维数据配准，将配准后的颅骨和面皮转化为二维深度图；1. 2)张量特征提取，将颅骨和面皮分别变换到颅骨张量特征空间和面皮张量特征空间；1.3)回归模型建立，采用偏最小二乘回归分性别建立面皮张量特征对颅骨张量特征和年龄、体重指数BMI属性的回归模型；2)颅面重构2.1)数据处理把待重构的未知颅骨统一到法兰克福坐标系并将大小规格化，将未知颅骨与训练样本中的颅骨进行数据配准，将配准后的未知颅骨转化为二维深度图；2. 2)张量特征提取将未知颅骨变换到颅骨张量特征空间；2. 3)特征回归，估计颅骨的性别和年龄、体重指数BMI属性，根据其张量特征和属性， 利用1.3)中建立的回归模型回归出相应的面皮张量特征；2. 4)三维面皮重构，根据回归到的面皮张量特征恢复其二维深度图，将得到的二维深度图转化为三维点云数据。
2.按照权利要求1所述的方法，其特征在于，所述第1. 2)步中的张量特征空间是由训练样本集中全体样本所确定的特征空间，其确定过程如下按照颅面的生理特征将每个颅骨和每个面皮的深度图分为大小相同且部分重叠的小区域，分别将训练样本集中全体颅骨和全体面皮的深度图按身份、区域位置、区域的纵向和横向信息变化的方向组织为一个四阶张量Ds e Rnxkxhx^P Df e Rmxkxhxw,其中，η为训练样本集中颅骨的个数，m为面皮的个数，k为每个深度图所确定的小区域个数，h和w分别是每个小区域的行数和列数，通过执行张量分解和空间降维得下式(1)和O)Ds ^ Cs X iUpeople X 2UregionX 3Urow ^ ^column — ^s X iUpeople ......(丄)其中CS是核张量，控制着身份、区域位置、区域纵向和横向这四个模式之间的相互作用，U-fe eiT ,Mregwn^Rkyk' ,VrW、，是四个行正交的矩阵，由其行向量分别确定了相应的模式特征空间，ns，ks，hs，ws分别是相应空间经过降维以后的空间维数，Ts — CsX2UregionX3UrowX4Ucolumn 是颅骨身份张量； Df ^ Cf X ιVpeople X 2Vregion X 3Vrow X 4VCoiumn = TfX1Ypec5ple (2)其中cf是核张量，W NregwnWf ’V_GRh ，V-_ eiT，是四个行正交的矩阵，由其行向量分别确定了相应的模式特征空间，nf, kf，hf，wf分别是相应空间经过降维以后的空间维数，Tf = CfX2VregionX3VrowX4Vcolm是面皮身份张量；定义颅骨原始特征向量
3.按照权利要求1所述的方法，其特征在于，所述第2.4)步中根据回归到的面皮张量特征恢复其深度图如下J = ^ypeople ……(5)其中为回归到的面皮张量特征，$为恢复出的面皮原始特征向量，根据面皮原始特征向量的定义得到面皮的深度图。
4.按照权利要求2所述的方法，其特征在于，所述训练样本集可以包含一部分仅有颅骨或仅有面皮的样本，但必须包含一定数量的既有颅骨又有面皮的颅面样本。
全文摘要
本发明涉及计算机应用技术领域的一种通过在张量空间进行偏最小二乘回归来实现对未知颅骨进行面皮重构的方法。所述方法分别将训练样本中全体颅骨和全体面皮组织为两个四阶张量，通过张量分解将颅骨和面皮变换到不同的张量特征空间，利用偏最小二乘技术分性别进行面皮张量特征对颅骨张量特征的回归建模。对待重构的目标颅骨，按照与训练样本相同的方式变换到颅骨张量特征空间，利用建立的回归模型得到相应的面皮张量特征，根据张量特征重构出目标面皮。本发明所述方法充分体现了颅骨决定面貌这一本质联系，避免了以往基于统计模型的重构中求解数据缺失这一病态问题。在刑侦、考古、法医人类学、医学手术、美容整容等领域具有重要的应用价值。
文档编号G06T17/00GK102521875SQ20111038209
公开日2012年6月27日 申请日期2011年11月25日 优先权日2011年11月25日
发明者周明全, 武仲科, 段福庆 申请人:北京师范大学