[发明专利]一种基于张量空间偏最小二乘回归的颅面重构方法

说明书

技术领域

本发明涉及面皮重构方法，特别涉及一种通过在张量空间进行偏最小二乘回归来实现对未知颅骨进行面皮重构的方法。主要用于刑侦、考古、法医人类学、医学手术、美容整容等领域。

背景技术

颅骨是人类面貌的内在生物特征，根据颅骨形态实现面貌重构是法医人类学领域的典型应用项目。传统的手工重构是在待定颅骨的基础上，借助解剖测量和针刺测量途径获得人脸软组织的厚度规律，由艺术家雕塑而成，属于人类学家、艺术家和医生的领域，较多的依赖重构者个人对颅面特征和艺术的把握，其结果受主观因素影响较大，不具备科学推广的价值，并且重构一个颅面需要几天的时间。目前，计算机辅助的颅面三维重构已引起国际信息学、人类学、法医学等相关领域的高度重视，成为国际上研究的热点。

典型技术可分为两类：一类是基于软组织厚度测量的方法。通过测量一定数量的颅面样本在颅骨上一些关键特征点处的平均软组织厚度作为经验知识，对于目标颅骨，首先根据这些平均软组织厚度的经验值，获得颅骨关键特征点处的面皮位置，然后通过插值、变形等技术得到整个颅骨的面皮。这类方法只要求在颅骨上标注少量关键特征点，并获得这些点的软组织厚度，所需的数据较少。然而，这些关键特征点处的软组织厚度和颅面形状之间并没有直接的相关关系；另一类是基于统计变形模型的方法。这类方法将一套颅骨和面皮数据组织为一个整体形成一个向量，通过对大量的训练样本进行主成份分析(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张量特征提取：