Skull identification based on least square canonical correlation analysis

基于最小二乘正则相关性分析的颅骨身份识别 Skull identification based on least square canonical correlation analysis 1 2 corresp first-author E-mail： nwuzmq@163.com E-mail： nwuzmq@163.com 周明全（1954-），男，陕西西安人，教授，博士生导师，陕西省千人计划特聘教授，1976、1988年于西北大学分别获得学士、硕士学位，主要从事虚拟现实、可视化技术、智能信息处理信息处理等领域的基础理论及应用工程创新的研究。E-mail： nwuzmq@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444914&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444918&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444916&type= 22.01330566 32.00404358 周明全（1954-），男，陕西西安人，教授，博士生导师，陕西省千人计划特聘教授，1976、1988年于西北大学分别获得学士、硕士学位，主要从事虚拟现实、可视化技术、智能信息处理信息处理等领域的基础理论及应用工程创新的研究。E-mail： nwuzmq@163.com nwuzmq@163.com 1 杨　稳（1993-），男，陕西安康人，博士研究生，主要从事机器学习与模式识别方面的研究。E-mail： yw@stumail.nwu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444920&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444924&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444922&type= 22.01333618 32.00402832 杨　稳（1993-），男，陕西安康人，博士研究生，主要从事机器学习与模式识别方面的研究。E-mail： yw@stumail.nwu.edu.cn yw@stumail.nwu.edu.cn 1 1 1 1 西北大学 信息科学与技术学院，陕西 西安 710127 西北大学 信息科学与技术学院，陕西 西安 710127 College of Information Science and Technology， Northwest University， Xi'an 710127， China College of Information Science and Technology， Northwest University， Xi'an 710127， China 北京师范大学 信息科学与技术学院，北京 100875 北京师范大学 信息科学与技术学院，北京 100875 College of Information Science and Technology， Beijing Normal University， Beijing 100875， China College of Information Science and Technology， Beijing Normal University， Beijing 100875， China 为了可靠地测量颅骨与面皮之间的相关性，提高颅骨识别能力，本文提出了一种基于最小二乘正则相关性分析的颅骨识别方法。构建颅骨和面皮的统计形状模型，将高维颅骨和面皮映射到低维形状参数空间；基于最小二乘正则相关性分析提取颅骨和面皮的主相关信息，构建整体相关性分析模型，用于度量颅骨与面皮的整体相关性；然后，考虑到颅骨不同区域的相关度不同，将颅骨分为前额、眼睛、鼻子、嘴巴和轮廓5个区域，基于最小二乘正则相关性分析提取颅骨和面皮的主相关信息，分区域构建区域相关性分析模型，用于度量颅骨与面皮之间的局部细节相关性；最后，分别采用整体相关性分析模型和区域相关性分析模型度量颅骨与面皮的匹配关系，计算颅骨与面皮库中每个面皮之间的匹配分数，具有最高匹配分数的面皮即为正确的识别结果，从而实现颅骨身份识别。实验结果表明，整体相关性分析模型的识别准确率为85.2%，单一区域相关性分析模型中轮廓区域的识别准确率最高，鼻子区域的识别准确率最低；通过融合5个区域建立的相关性分析模型，识别准确率高达95.2%，基于区域融合的方法优于整体相关性分析方法。 To measure the correlation between the skull and the face reliably and improve the skull recognition ability， a skull identification method based on least squares canonical dependency analysis （LSCDA） was proposed. First， a statistical shape model of the skull and facial skin was constructed， and the high-dimensional skull and facial skin were mapped to the low-dimensional shape parameter space. Second， the main relevant information about the skull and facial skin was extracted based on LSCDA， and an overall correlation analysis model was constructed to measure the relationship between the skull and the skin. By considering the difference in the correlations between different regions of the skull， the skull was divided into five regions： the forehead， eyes， nose， mouth， and contour. Based on LSCDA， the main relevant information about the skull and facial skin was extracted and constructed by region. A regional correlation analysis model was used to measure the local detail correlation between the skull and the face. Finally， a global correlation analysis model and a regional correlation analysis model were used to measure the matching relationship between the skull and facial skin， and the matching score between each face in the skull and facial skin database was calculated. The face with the highest matching score yielded the correct recognition result for achieving skull identification. The experimental results reveal that the recognition accuracy of the overall correlation analysis model is 85.2%. The recognition accuracy of the contour region in the single region correlation analysis model is the highest， whereas that of the nose region is the lowest. The correlation analysis established by fusing the five regional models indicates that the recognition accuracy rate is as high as 95.2% and that the method based on regional fusion is better than the overall correlation analysis method. 图像重建 颅骨身份识别 统计形状模型 主相关信息 最小二乘正则相关性分析 image reconstruction skull recognition statistical shape model main related information least square canonical dependency analysis 1 引　言 1 颅骨身份识别是法医学中一个重要的研究课题。在很多法医学案例中，技术人员只能获取到受害人的颅骨，并且没有其他证据线索，这种情况下无法使用传统的身份识别技术，颅骨识别 ［ 1 1 ］ 便成为受害者身份认证的重要技术手段之一。随着CT、3D扫描等数字化技术的发展，三维颅面采集变得更加容易，计算机辅助颅骨识别成为研究热点，内容涉及法医学、信息学、人类学和计算机图形学等多个领域。 目前，颅骨识别研究主要集中在颅像重合 ［ 2 2 ］ 和颅面复原 ［ 3 3 ］ 两方面。颅像重合是将未知身源颅骨的2D投影与失踪人2D照片相叠加，根据其相关特征找到一个合适的匹配，从而获取到颅骨人脸叠加像，即可判断颅骨和人脸照片是否属于同一人。颅面复原 ［ 4 4 ］ 是通过分析颅骨和面皮之间的内在关系，复原出未知颅骨面貌的技术。当前主流的颅面复原方法采用统计学习技术，如统计形状模型 ［ 5 5 ］ 和统计回归模型 ［ 6 6 ］ 等。 颅像重合 ［ 7 7 - 11 11 7-11 ］ 和颅面复原 ［ 12 12 - 13 13 12-13 ］ 都需要利用形态学知识准确提取颅骨和面皮之间的内在关系，存在不确定性较高和识别能力较低等问题。随着基于图像的3D人脸建模 ［ 14 14 ］ 和3D数据采集技术 ［ 15 15 ］ 的发展，使用3D颅骨和3D面皮之间的相关性识别颅骨逐渐成为研究热点问题。Duan ［ 16 16 ］ 等人首次提出通过分析颅骨和面皮之间的形态相关度来识别未知颅骨，并使用典型关联分析 ［ 17 17 ］ 提取映射后两组相关性最大的主成分向量，用于匹配颅骨和面皮。但由于颅骨与面皮之间关系复杂，该方法的识别率有待提高。 本文提出了一种基于最小二乘正则相关性分析（Least-Squares Canonical Dependency Analysis，LSCDA）的颅骨识别方法。该方法通过构建颅骨和面皮的相关性分析模型，计算未知颅骨和面皮库中每张面皮的匹配程度，得到识别结果。与颅像重合和颅面复原不同，该方法无需准确提取颅骨和面皮之间内在的复杂关系，仅测量它们之间的相关性即可确定颅骨身份，降低了颅骨识别问题的复杂性，提高了识别率。 2 数据获取与预处理 1 本研究使用了255套颅面数据，这些数据来源于西北大学可视化技术研究所，女性年龄在19~75岁，男性年龄在21~67岁。颅面CT数据由陕西中医药大学附属医院使用西门子多排探测器螺旋CT机采集，每个图像以DICOM格式存储，大小约为512 <math id="M1"><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444928&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444926&type= 3.72533321 2.79399991 512 <math id="M2"><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444934&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444931&type= 3.72533321 2.79399991 250。原始CT切片图像在去噪后，使用Sobel算子提取颅骨和面皮轮廓，采用Marching Cube算法 ［ 17 17 ］ 进行重建，得到颅骨和面皮的三维模型。颅骨和面皮分别表示为大约包含220 000和150 000个顶点的三角形网格。另外，每套颅面还记录了年龄、性别和BMI等属性。 为了消除因位置、姿态和尺度等因素造成的数据不一致，所有样本都被转换在统一的法兰克福坐标系下。法兰克福坐标系是由左右双侧耳门上点、左眼眶下边缘点和眉间点4个颅骨特征点确定的，分别用 L p ， R p ， M p 和 V p 表示。法兰克福平面由 L p ， R p 和 M p 3点确定 ［ 18 18 ］ 。如 图1 图1 所示：以 <math id="M3"><msub><mrow><mi mathvariant="bold-italic">L</mi></mrow><mrow><mi>p</mi></mrow></msub><msub><mrow><mi mathvariant="bold-italic">R</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444938&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444936&type= 7.78933334 3.72533321 为法向量且过点 <math id="M4"><msub><mrow><mi>V</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444942&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444940&type= 4.14866638 3.81000018 的平面与直线 <math id="M5"><msub><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msub><msub><mrow><mi>R</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444950&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444946&type= 7.78933334 3.72533321 的交点记为法兰克福平面的原点 O ′； <math id="M6"><msub><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444955&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444953&type= 3.55599999 3.72533321 ， <math id="M7"><msub><mrow><mi>R</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444959&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444957&type= 3.81000018 3.72533321 以及 <math id="M8"><msub><mrow><mi>M</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444966&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444963&type= 4.57200003 3.72533321 3点所构成的平面为坐标系中的 XO ′ Y 平面，左耳门上点 <math id="M9"><msub><mrow><mi>L</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444971&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444969&type= 3.55599999 3.72533321 到右耳门上点 <math id="M10"><msub><mrow><mi>R</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444979&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444972&type= 3.81000018 3.72533321 的方向为 X 轴方向；过原点 O ′且与 XO ′ Y 平面垂直向上的方向记为 Z 轴；同时垂直于 X 轴、 Z 轴的直线记为坐标系的 Y 轴， Y 轴的正方向由右手法则确定。 统一坐标系后，再进行尺度归一化。设定所有颅骨模型的 L p 到 R p 的距离为单位1，则对颅骨每一个顶点（ x，y，z ）进行尺度变换为（ x/|L p R p | ， y/|L p R p | ， z/|L p R p | ）。 由于面部特征信息主要集中于前半部分 ［ 19 19 ］ ，为了消除不相关数据的影响和减少算法耗时，对颅面样本进行分割，只保留与识别有直接关系的前半部分，分割后的颅骨和面皮分别为包含约160 000和80 000个顶点的三角形网格。为了建立相同的点对应关系，选取一组分割后的颅骨和面皮作为参考模板，如 图2 图2 所示。其余所有颅骨模型和面皮模型通过非刚性配准方法 ［ 20 20 ］ 自动配准到参考模型中。 3 构建统计形状模型 1 在进行LSCDA分析前，需要分别为颅骨和面皮建立统计形状模型。统计形状模型 ［ 21 21 ］ 是一个参数化的可变形模板，每个颅面都有一个参数化的表示，即每个颅面对应着参数空间中的一个点。主成分分析（Principal Components Analysis，PCA）是建立颅面统计形状模型的有力工具。PCA能够反映训练样本数据集真实分布的一组正交轴，这组正交轴原点设在数据集的中心，表示训练数据中具有最大方差的一组方向，数据能够以最小误差重构，从而达到降维的效果。 通过连接所有顶点的坐标，颅骨和面皮可以分别表示为一个高维度的向量。因此，得到一个颅骨的训练数据集 <math id="M11"><mo stretchy="false">{</mo><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>=</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>x</mi></mrow><mrow><mi>i</mi><mn mathvariant="normal">1</mn></mrow><mrow><mi>S</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi mathvariant="normal">y</mi></mrow><mrow><mi>i</mi><mn mathvariant="normal">1</mn></mrow><mrow><mi>S</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>z</mi></mrow><mrow><mi>i</mi><mn mathvariant="normal">1</mn></mrow><mrow><mi>S</mi></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mrow><mi>x</mi></mrow><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>S</mi></mrow></msubsup><mo>,</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444995&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444994&type= 39.36999893 4.91066647 <math id="M12"><msubsup><mrow><mi>y</mi></mrow><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>S</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>z</mi></mrow><mrow><mi>i</mi><mi>m</mi></mrow><mrow><mi>S</mi></mrow></msubsup><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mo stretchy="false">|</mo><mi>i</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo stretchy="false">}</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444997&type= 36.99933624 4.99533367 和一个面皮的训练数据集 <math id="M13"><mo stretchy="false">{</mo><msub><mrow><mi mathvariant="bold-italic">F</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>=</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>x</mi></mrow><mrow><mi>i</mi><mn mathvariant="normal">1</mn></mrow><mrow><mi>F</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>y</mi></mrow><mrow><mi>i</mi><mn mathvariant="normal">1</mn></mrow><mrow><mi>F</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>z</mi></mrow><mrow><mi>i</mi><mn mathvariant="normal">1</mn></mrow><mrow><mi>F</mi></mrow></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mrow><mi>x</mi></mrow><mrow><mi>i</mi><mi>n</mi></mrow><mrow><mi>F</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>y</mi></mrow><mrow><mi>i</mi><mi>n</mi></mrow><mrow><mi>F</mi></mrow></msubsup><mo>,</mo><msubsup><mrow><mi>z</mi></mrow><mrow><mi>i</mi><mi>n</mi></mrow><mrow><mi>F</mi></mrow></msubsup><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup><mo stretchy="false">|</mo><mi>i</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo stretchy="false">}</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445004&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445001&type= 75.01467133 4.99533367 ，其中 N 表示训练样本个数， m 和 n 分别为颅骨和面皮的顶点数，每个坐标索引标记训练集上的对应点。假定 <math id="M14"><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445009&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445005&type= 19.72733307 4.14866638 是均值归一化颅骨协方差矩阵的单位正交特征向量， <math id="M15"><mover><mrow><mi mathvariant="bold-italic">S</mi></mrow><mo>¯</mo></mover><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445013&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445011&type= 16.08666611 8.46666718 是样本颅骨向量的均值。PCA变换可以将每一个颅骨向量 <math id="M16"><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445017&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445015&type= 3.04799986 3.72533321 通过变换投影到以均值 <math id="M17"><mover><mrow><mi mathvariant="bold-italic">S</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445020&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445019&type= 2.37066674 3.38666677 为原点，以这些主方向为坐标轴的坐标系中。那么，颅骨统计形状模型可以表示为： <math display="block" id="M18"><mi mathvariant="bold-italic">S</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">a</mi><mo stretchy="false">)</mo><mo>=</mo><mover><mrow><mi mathvariant="bold-italic">S</mi></mrow><mo>¯</mo></mover><mo>+</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>p</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>a</mi></mrow><mrow><mi>i</mi></mrow></msub><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445031&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445027&type= 32.00400162 8.88999939 ， 其中： <math id="M19"><mi mathvariant="bold-italic">a</mi><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mrow><mi>a</mi></mrow><mrow><mi>p</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445035&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445033&type= 29.04066658 4.65666676 是组合参数， <math id="M20"><msub><mrow><mi>a</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445039&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445037&type= 2.62466669 3.72533321 满足高斯分布，在模型中只要给定组合参数 <math id="M21"><mi mathvariant="bold-italic">a</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445049&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445042&type= 1.94733334 2.79399991 即可产生颅骨参数模型 <math id="M22"><mi>S</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">a</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445056&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445053&type= 7.70466709 3.55599999 。 所有颅骨数据和面皮数据均进行了非刚性配准操作，具有相同的点对应关系，可以通过PCA变换来确定形状模型参数。对于每一个颅骨向量 <math id="M23"><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445064&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445060&type= 3.04799986 3.72533321 ，变换为： <math display="block" id="M24"><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>-</mo><mover><mrow><mi mathvariant="bold-italic">S</mi></mrow><mo>¯</mo></mover><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445071&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445068&type= 26.50066566 4.23333359 ， 其中： <math id="M25"><mi mathvariant="bold-italic">P</mi><mo>=</mo><mo stretchy="false">[</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">P</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445078&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445074&type= 30.14133263 4.48733330 ， <math id="M26"><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>∈</mo><msup><mrow><mi mathvariant="bold">R</mi></mrow><mrow><mi>p</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445084&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445081&type= 10.49866676 4.23333359 称为颅骨向量 <math id="M27"><msub><mrow><mi mathvariant="bold-italic">S</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445091&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445087&type= 3.04799986 3.72533321 的主成分向量。类似地，可通过相同的处理过程来构建面皮统计形状模型： <math display="block" id="M28"><mi mathvariant="bold-italic">F</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">b</mi><mo stretchy="false">)</mo><mo>=</mo><mover><mrow><mi mathvariant="bold-italic">F</mi></mrow><mo>¯</mo></mover><mo>+</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>q</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>b</mi></mrow><mrow><mi>j</mi></mrow></msub><msub><mrow><mi mathvariant="bold-italic">Q</mi></mrow><mrow><mi mathvariant="bold-italic">j</mi></mrow></msub></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445097&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445093&type= 32.85066605 9.31333351 ， 其中： <math id="M29"><mover><mrow><mi mathvariant="bold-italic">F</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445102&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445099&type= 2.62466669 3.38666677 是样本面皮向量的均值， <math id="M30"><msub><mrow><mi mathvariant="bold-italic">Q</mi></mrow><mrow><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><mi>j</mi><mo>=</mo><mn mathvariant="normal">1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>q</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445109&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445106&type= 28.27866745 4.65666676 是均值归一化面皮协方差矩阵的单位正交特征向量。 因此，通过构建统计形状模型，将颅骨和面皮数据分别投影到其低维形状参数空间 <math id="M31"><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445116&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445112&type= 3.13266659 3.72533321 和 <math id="M32"><msub><mrow><mi mathvariant="bold-italic">b</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445122&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445119&type= 3.13266659 3.72533321 中。 4 颅骨和面皮的主相关信息提取 1 通过构建统计形状模型，每个颅骨可以表示为 <math id="M33"><mi>p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445129&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445125&type= 1.86266661 3.64066648 维特征向量，每个面皮表示为 <math id="M34"><mi>q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445134&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445132&type= 1.86266661 3.64066648 维特征向量。LSCDA旨在通过平方损失互信息 ［ 22 22 ］ （Squared-loss Mutual Information，SMI）找到 <math id="M35"><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445141&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445137&type= 3.13266659 3.72533321 和 <math id="M36"><msub><mrow><mi mathvariant="bold-italic">b</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445144&type= 3.13266659 3.72533321 之间的统计相关关系，它可以在线性投影下捕获两个变量的主相关信息。 因此，用LSCDA找到 <math id="M37"><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445154&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445151&type= 3.13266659 3.72533321 和 <math id="M38"><msub><mrow><mi mathvariant="bold-italic">b</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445160&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445158&type= 3.13266659 3.72533321 的低维子空间，投影的相关性通过线性降维达到最大值，即有： <math display="block" id="M39"><mtable columnalign="left"><mtr><mtd><mi mathvariant="bold-italic">u</mi><mo>=</mo><mi mathvariant="bold-italic">U</mi><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></mtd></mtr><mtr><mtd><mi mathvariant="bold-italic">v</mi><mo>=</mo><mi mathvariant="bold-italic">V</mi><msub><mrow><mi mathvariant="bold-italic">b</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445169&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445167&type= 12.86933327 8.46666718 ， 其中 <math id="M40"><mi mathvariant="bold-italic">U</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445175&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445172&type= 2.96333337 2.79399991 和 <math id="M41"><mi mathvariant="bold-italic">V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445184&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445179&type= 2.96333337 2.96333337 是转换矩阵。 <math display="block" id="M42"><mtable columnalign="left"><mtr><mtd><mi mathvariant="bold-italic">U</mi><msup><mrow><mi mathvariant="bold-italic">U</mi></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></mtd></mtr><mtr><mtd><mi mathvariant="bold-italic">V</mi><msup><mrow><mi mathvariant="bold-italic">V</mi></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445190&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445186&type= 15.49400043 8.72066689 ， 其中 <math id="M43"><mi>p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445195&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445191&type= 1.86266661 3.64066648 和 <math id="M44"><mi>q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445199&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445197&type= 1.86266661 3.64066648 分别是 <math id="M45"><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445202&type= 2.11666679 2.79399991 和 <math id="M46"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445210&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445206&type= 2.03200006 2.79399991 的投影维数， <math id="M47"><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445218&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445214&type= 2.62466669 3.72533321 是 <math id="M48"><mi>p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445224&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445221&type= 1.86266661 3.64066648 维单位矩阵， <math id="M49"><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445235&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445228&type= 2.70933342 3.72533321 是 <math id="M50"><mi>q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445245&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445241&type= 1.86266661 3.64066648 维单位矩阵。 由此可知，寻找转换矩阵 <math id="M51"><mi mathvariant="bold-italic">U</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445254&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445251&type= 2.96333337 2.79399991 和 <math id="M52"><mi mathvariant="bold-italic">V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445262&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445259&type= 2.96333337 2.96333337 使得 <math id="M53"><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445272&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445268&type= 2.11666679 2.79399991 和 <math id="M54"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445275&type= 2.03200006 2.79399991 之间的相关关系最大化。这里通过 <math id="M55"><mi mathvariant="normal">S</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445285&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445281&type= 6.94266701 2.87866688 来测量其相关性： <math display="block" id="M56"><mi>S</mi><mi>M</mi><msub><mrow><mi>I</mi></mrow><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>∬</mo><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo><mo>×</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></mrow></mfrac><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445291&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445289&type= 64.17733002 8.72066689 ， 其中： <math id="M57"><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445293&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445292&type= 14.64733315 4.65666676 是联合概率分布密度， <math id="M58"><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445296&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445295&type= 8.80533314 4.65666676 和 <math id="M59"><msub><mrow><mi>p</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445298&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445297&type= 8.80533314 4.65666676 分别是 <math id="M60"><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445299&type= 2.11666679 2.79399991 和 <math id="M61"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445302&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445301&type= 2.03200006 2.79399991 的边缘概率分布密度。 然而， <math id="M62"><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445306&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445304&type= 14.64733315 4.65666676 ， <math id="M63"><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445311&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445308&type= 8.80533314 4.65666676 和 <math id="M64"><msub><mrow><mi>p</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445316&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445314&type= 8.80533314 4.65666676 通常是未知的，不能直接最大化 <math id="M65"><mi mathvariant="normal">S</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445322&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445320&type= 6.94266701 2.87866688 。因此，本文采用最小二乘互信息（Least-Squares Mutual Information， LSMI）方法直接估计密度比，而不用估计每个密度，即有： <math display="block" id="M66"><msup><mrow><mi>g</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo><mo>×</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445326&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445324&type= 40.72466660 8.72066689 . 密度比函数 <math id="M67"><msup><mrow><mi>g</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445333&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445329&type= 12.95399952 3.97933340 由线性组合建模表示为 <math id="M68"><mi>g</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445338&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445336&type= 11.59933376 3.72533321 ，具体如下： <math display="block" id="M69"><mi>g</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>=</mo><msup><mrow><mi mathvariant="bold-italic">α</mi></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup><mi>ϕ</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445344&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445341&type= 31.66533279 3.97933340 ， 其中 <math id="M70"><mi mathvariant="bold-italic">α</mi><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>b</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445348&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445346&type= 29.12533188 4.65666676 是从样本中学习的参数， <math id="M71"><mi>ϕ</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>,</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>b</mi></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445353&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445351&type= 68.91867065 4.91066647 是一个基函数，并且对于所有的 <math id="M72"><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445360&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445356&type= 2.11666679 2.79399991 和 <math id="M73"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445364&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445362&type= 2.03200006 2.79399991 ， <math id="M74"><mi>ϕ</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>≥</mo><msub><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445369&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445366&type= 19.30399895 4.74133348 ， <math id="M75"><msub><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445374&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445372&type= 3.04799986 3.72533321 是一个 b 维零向量。 参数 <math id="M76"><mi mathvariant="bold-italic">α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445379&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445377&type= 2.11666679 2.79399991 是通过最小化平方误差 <math id="M77"><mi>δ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445384&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445381&type= 1.94733334 2.79399991 来估计的，即： <math display="block" id="M78"><mi>δ</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">α</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>∬</mo><mo stretchy="false">(</mo><mi mathvariant="normal">g</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>-</mo><msup><mrow><mi mathvariant="normal">g</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><msup><mrow><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445388&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445386&type= 79.58666992 7.45066643 . 这里，令： <math display="block" id="M79"><mi mathvariant="bold-italic">H</mi><mo>=</mo><mo>∬</mo><mi>ϕ</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mi>ϕ</mi><msup><mrow><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445394&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445391&type= 65.53199768 6.68866634 ， <math display="block" id="M80"><mi mathvariant="bold-italic">h</mi><mo>=</mo><mo>∬</mo><mi>ϕ</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">u</mi><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445397&type= 45.97399902 6.68866634 . 于是有： <math display="block" id="M81"><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">α</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">α</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">H</mi><mi mathvariant="bold-italic">α</mi><mo>-</mo><msup><mrow><mi mathvariant="bold-italic">h</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445408&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445405&type= 37.42266846 7.45066643 . 近似估计 <math id="M82"><mi mathvariant="bold-italic">H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445416&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445411&type= 3.04799986 2.79399991 和 <math id="M83"><mi mathvariant="bold-italic">h</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445422&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445419&type= 2.20133328 2.79399991 的值如下： <math display="block" id="M84"><mover accent="true"><mrow><mi mathvariant="bold-italic">H</mi></mrow><mo>̂</mo></mover><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mi mathvariant="normal">𝓁</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>,</mo><mi>i</mi><mi mathvariant="normal">'</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">𝓁</mi></mrow></munderover></mstyle><mrow><mi>ϕ</mi><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">'</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mi>ϕ</mi><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi><mi mathvariant="bold-italic">'</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">Τ</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445430&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445427&type= 48.17533112 8.88999939 ， <math display="block" id="M85"><mover accent="true"><mrow><mi mathvariant="bold-italic">h</mi></mrow><mo>̂</mo></mover><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">𝓁</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi mathvariant="normal">𝓁</mi></mrow></munderover></mstyle><mrow><mi>ϕ</mi></mrow><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo>,</mo><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="bold-italic">i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445438&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445435&type= 29.12533188 8.46666718 . 将问题转换为： <math display="block" id="M86"><mover accent="true"><mrow><mi mathvariant="bold-italic">α</mi></mrow><mo>̂</mo></mover><mo>=</mo><munder><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">g</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi mathvariant="bold">α</mi><mo>∈</mo><msup><mrow><mi>R</mi></mrow><mrow><mi>b</mi></mrow></msup></mrow></munder><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="bold-italic">α</mi><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">α</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mi mathvariant="bold-italic">α</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445444&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445441&type= 45.21199799 8.88999939 . 为了进一步简化，增加了 <math id="M87"><mfrac><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msup><mrow><mi mathvariant="bold-italic">α</mi></mrow><mrow><mi mathvariant="bold">T</mi></mrow></msup><mi mathvariant="bold-italic">α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445452&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445449&type= 9.22866726 7.53533316 （ <math id="M88"><mi>λ</mi><mo>></mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445458&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445455&type= 8.21266651 2.79399991 ）。 通过计算上述目标函数的导数并将它设为零，可以得到： <math display="block" id="M89"><mover accent="true"><mrow><mi mathvariant="bold-italic">α</mi></mrow><mo>̂</mo></mover><mo>=</mo><mo stretchy="false">(</mo><mover accent="true"><mrow><mi mathvariant="bold-italic">H</mi></mrow><mo>̂</mo></mover><mo>+</mo><mi>λ</mi><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mover accent="true"><mrow><mi mathvariant="bold-italic">h</mi></mrow><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445466&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445461&type= 27.77066803 5.41866684 ， 其中 <math id="M90"><msub><mrow><mi mathvariant="bold-italic">I</mi></mrow><mrow><mi mathvariant="bold-italic">b</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445474&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445470&type= 2.62466669 3.72533321 是 b 阶单位矩阵。 那么 <math id="M91"><mi mathvariant="normal">S</mi><mi mathvariant="normal">M</mi><mi mathvariant="normal">I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445481&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445478&type= 6.94266701 2.87866688 的近似值表示如下： <math display="block" id="M92"><mi>S</mi><mover accent="true"><mrow><mi>M</mi></mrow><mo>̂</mo></mover><mi>I</mi><mo>=</mo><msup><mrow><mover accent="true"><mrow><mi mathvariant="bold-italic">h</mi></mrow><mo>̂</mo></mover></mrow><mrow><mi mathvariant="bold">T</mi></mrow></msup><mover accent="true"><mrow><mi mathvariant="bold-italic">α</mi></mrow><mo>̂</mo></mover><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mover accent="true"><mrow><mi mathvariant="bold-italic">α</mi></mrow><mo>̂</mo></mover><mover accent="true"><mrow><mi mathvariant="bold-italic">H</mi></mrow><mo>̂</mo></mover><mover accent="true"><mrow><mi mathvariant="bold-italic">α</mi></mrow><mo>̂</mo></mover><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445490&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445485&type= 40.89400101 7.45066643 . 综上所述，颅骨和面皮之间的主相关关系提取算法过程如下： 首先，最合适的变换矩阵 <math id="M93"><mi mathvariant="bold-italic">U</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445498&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445494&type= 2.96333337 2.79399991 和 <math id="M94"><mi mathvariant="bold-italic">V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445509&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445502&type= 2.96333337 2.96333337 是以迭代的方式确定的，具体步骤如下： Step1：根据 式（4） 式（4） 初始化 <math id="M95"><mi mathvariant="bold-italic">U</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445517&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445513&type= 2.96333337 2.79399991 和 <math id="M96"><mi mathvariant="bold-italic">V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445527&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445523&type= 2.96333337 2.96333337 ； Step2：根据 式（16） 式（16） 对参数 <math id="M97"><mi mathvariant="bold-italic">α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445536&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445532&type= 2.11666679 2.79399991 进行优化； Step3：采用序列二次规划技术来更新 <math id="M98"><mi mathvariant="bold-italic">U</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445547&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445541&type= 2.96333337 2.79399991 和 <math id="M99"><mi mathvariant="bold-italic">V</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445555&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445551&type= 2.96333337 2.96333337 ，重复此过程直至收敛。 然后，通过线性投影将颅骨形状参数转换成LSCDA空间参数，即： <math display="block" id="M100"><mi mathvariant="bold-italic">u</mi><mo>=</mo><mi mathvariant="bold-italic">U</mi><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="bold-italic">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445565&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445559&type= 12.86933327 3.72533321 . 同理，通过线性投影将面皮形状参数转换为LSCDA空间参数，即： <math display="block" id="M101"><mi mathvariant="bold-italic">v</mi><mo>=</mo><mi mathvariant="bold-italic">V</mi><msub><mrow><mi mathvariant="bold-italic">b</mi></mrow><mrow><mi mathvariant="bold-italic">q</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445575&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445568&type= 12.86933327 3.89466667 . 最后，计算库中的未知颅骨与每张面皮的匹配分数，匹配得分最高的面皮就是识别结果。定义颅骨和面皮之间的匹配程度如下： <math display="block" id="M102"><mi>r</mi><mo stretchy="false">(</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo><</mo><mi mathvariant="bold-italic">u</mi><mo>,</mo><mi mathvariant="bold-italic">v</mi><mo>></mo></mrow><mrow><mfenced open="‖" close="‖" separators="|"><mrow><mi mathvariant="bold-italic">u</mi></mrow></mfenced><mfenced open="‖" close="‖" separators="|"><mrow><mi mathvariant="bold-italic">v</mi></mrow></mfenced></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445588&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445582&type= 31.58066750 8.80533314 ， 其中 <math id="M103"><mo><</mo><mo>⋅</mo><mo>,</mo><mo>⋅</mo><mo>></mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445597&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445593&type= 11.76866722 3.21733332 表示内积。 5 基于LSCDA的颅骨识别方法 1 未知颅骨与面皮之间的识别将分别从整体和区域两个方面进行。首先，构建整体相关性分析模型，从全局来分析颅骨和面皮之间的关系；其次，由于颅面形态非常复杂，不同颅骨和面皮不同区域的相关度不同。例如，前额区域的相关性较强，因为该区域的软组织厚度较薄；而鼻子、眼睛和嘴巴等区域的相关性则较弱，因为这些区域的结构复杂。强关联区域具有较高的识别能力。然而，整体相关性分析模型只反映了全局意义上颅骨和面皮之间的相关性，并没有充分利用它们的特性信息。为了更加可靠地测量颅骨与面皮之间的相关性，提高识别能力，提出了一种基于区域的局部化识别方法。识别流程如 图3 图3 所示。 5.1　基于LSCDA的整体颅骨识别 2 基于LSCDA整体颅骨识别的目的是从三维面皮库中找到一个与未知颅骨最可能匹配的面皮。因此，需要计算未知颅骨与每个3D面皮之间的匹配分数。整体颅骨识别的框架如 图4 图4 所示，识别过程包括两个阶段：训练阶段和识别阶段。其中， <math id="M104"><mi mathvariant="bold-italic">S</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445620&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445616&type= 2.37066674 2.79399991 表示颅骨训练集， <math id="M105"><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445631&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445624&type= 2.62466669 2.79399991 表示面皮训练集， <math id="M106"><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445639&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445635&type= 2.11666679 2.79399991 和 <math id="M107"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445649&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445644&type= 2.03200006 2.79399991 分别表示LSCDA提取后的子空间向量。 训练阶段，建立相关性分析模型的步骤如下： Step1：使用所有训练数据构建颅骨和面皮的统计形状模型，如 式（1） 式（1） 和 式（3） 式（3） 所述； Step2：根据 式（2） 式（2） 将所有训练数据投影到形状参数空间中； Step3：针对训练颅骨和面皮的形状参数特征执行LSCDA，并获得两个映射向量 <math id="M108"><mi mathvariant="bold-italic">u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445669&type= 2.11666679 2.79399991 和 <math id="M109"><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445682&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445678&type= 2.03200006 2.79399991 。 识别阶段，使用相关性分析模型识别未知颅骨，步骤如下： Step1：通过上述的统计形状模型将未知颅骨和面皮库中的面皮分别投影到形状参数空间中； Step2：根据 式（4） 式（4） ，将它们的形状参数特征投影到LSCDA子空间中。 Step3：根据 式（20） 式（20） 计算未知颅骨与面皮库中每个面皮之间的匹配分数，具有最高匹配分数的面皮即为识别结果。 5.2　基于LSCDA的区域颅骨识别 2 如 图5 图5 所示，将颅骨和面皮分成几个相应的生理特征区域，即前额、眼睛、鼻子、嘴巴和轮廓。训练集之间的对应关系已经建立，只需要手动分割一个样本，其他样本可以通过点对应关系实现自动分割。基于LSCDA的区域颅骨识别框架与上述整体颅骨识别框架大致相同，采用该框架分别为5个区域构建相关性分析模型。在识别阶段，可以单区域相关性分析进行未知颅骨的识别；还可以将5个相关性分析模型进行融合，进一步提高未知颅骨的识别率。 对于一个颅骨和面皮匹配对，可以使用5个相关性分析模型获得5个区域匹配分数，最终匹配分数通过加权和计算，如 式（21） 式（21） 所示： <math display="block" id="M110"><mi>R</mi><mo>=</mo><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub><mo>+</mo><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><mo>+</mo><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub><mo>+</mo><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub><mo>+</mo><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445701&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445698&type= 62.06066895 3.81000018 ， 其中： <math id="M111"><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445709&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445705&type= 2.62466669 3.72533321 ， <math id="M112"><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445715&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445712&type= 2.28600001 3.72533321 ， <math id="M113"><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445722&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445718&type= 2.53999996 3.72533321 ， <math id="M114"><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445738&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445730&type= 3.21733332 3.72533321 和 <math id="M115"><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445744&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445741&type= 2.62466669 3.72533321 分别代表轮廓、前额、眼睛、嘴巴和鼻子的区域匹配分数； <math id="M116"><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445752&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445748&type= 4.06400013 3.72533321 ， <math id="M117"><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445760&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445756&type= 3.72533321 3.72533321 ， <math id="M118"><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">e</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445770&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445766&type= 3.97933340 3.72533321 ， <math id="M119"><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445781&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445775&type= 4.65666676 3.72533321 和 <math id="M120"><msub><mrow><mi>w</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445793&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445786&type= 4.06400013 3.72533321 分别表示各个区域相对应的权值。 如何选择 式（21） 式（21） 中每个区域的权值对颅骨识别起到决定性的作用。由于颅骨识别依赖的是颅骨与面皮之间的相关性，因此应该将高权重分配给高相关的区域，如轮廓和前额。相比之下，颅骨和面皮之间的关系在鼻子，眼睛和嘴巴等相关性较弱的区域是不确定的，但是这些区域可以为人脸识别提供更多的判别信息。 6 实验结果与分析 1 本文使用第2节介绍的三维颅面数据进行实验，并采用五折交叉验证来评估本文方法的有效性。将数据集分为5组，每组包括51个三维颅骨和面皮对。对于每一次折叠，选择其中一组51个三维颅骨面皮对作为测试集，剩余的204个颅骨和面皮构成训练集，也就是说，在5次折叠中共有255个三维颅骨作为测试集。对于每一个测试颅骨，所有的255个面皮构成查询集。整个实验的正确识别率是5次验证的平均值。 6.1　整体相关性分析 2 在PCA中大特征值对应的主成分向量表示训练数据集的变化模式，而小特征值对应的主成分向量则表示噪声方向。在这里，本文首先讨论了统计形状模型的方差贡献率对颅骨识别率的影响。由 式（1） 式（1） 所述，模式数量由从累积特征值计算的方差贡献率确定，当方差贡献率从95%变化到99.5%时，整体模型的正确识别率如 表1 表1 所示。从 表1 表1 中可以看出，随着方差贡献率的增加，识别率逐渐提高；在方差贡献率达到99%时，识别率达到最高；当方差贡献率再继续增加时，识别率反而会下降。 VCR/% 95 96 97 98 99 99.5 整体识别率/% 58 69 77 82 85 81 对于每次折叠，构建一个整体统计形状模型，当方差贡献率为99%时，颅骨和面皮数据投影到低维形状参数空间中，分别表示为548和356维特征向量。基于LSCDA的全局颅骨识别通过提取基向量，正确识别率达85.2%，由此表明基于颅骨和面皮相关性分析的颅骨识别方法是有效的。进一步分析识别错误发现，误匹配面皮与真实面皮之间的全局形态相似性较高，从而导致了误匹配。 图6 图6 为两个面皮之间的几何距离比较，从误差分析图可以看出，相似人脸与真实人脸大部分区域的几何距离都很小。由于两个面皮的全局形态相似性程度高，因此具有相似的形状模型参数，使得整体模型产生错误匹配。 6.2　区域相关性分析 2 为了验证区域相关性分析模型，分别为每个区域构建区域统计形状模型，然后进行区域LSCDA映射，完成匹配。当方差贡献率为99%时，与基于整体相关性分析的LSCDA颅骨识别相对比，可以得到区域相关性分析的颅骨识别率，如 表2 表2 所示。 区域 轮廓 前额 眼睛 鼻子 嘴巴 正确识别率 86.5 80 56 26 38 % 通过对比发现，轮廓、前额和整体模型的识别率相当，这是因为轮廓和前额在整体模型相关性测量中起着决定性作用。鼻子主要由软骨组成，嘴唇与口腔区域的牙齿分离，因此与嘴巴、鼻子区域相比，眼睛的形态对眼眶骨形状的依赖性更强，从而识别率也较高。 表3 表3 为区域相关性分析的正确识别率， 图7 图7 显示了不同方差贡献率的局部模型识别率的变化趋势。 VCR/% 95 96 97 98 99 99.5 轮廓 57 67 74 80 86.5 85 前额 34 39.5 61 63 80 83 眼睛 43 52 56.5 61 56 49 鼻子 23.5 30 37 28.5 26 24 嘴巴 29 32 36.5 42 38 23 % 从 表3 表3 和 图7 图7 可以看出，当方差贡献率增加时，前额和轮廓的识别率变化很快，尤其是前额区域，这是因为较多地舍弃了变化模式中包含的一些特性信息。由 图7 图7 可知，当轮廓、前额、眼睛、鼻子和嘴巴区域的方差贡献率分别为99%，99.5%，98%，97%和98%时，这些区域的识别率最高。在使用基于区域相关性分析进行颅骨识别时，对不同区域选择不同的方差贡献率来构建局部统计形状模型，可以提高颅骨识别的准确率。 此外，通过实验发现，各个区域权重经过融合之后将达到最优的识别准确率。区域权重如下： <math display="block" id="M121"><msub><mrow><mi>w</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>r</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mo>∑</mo><msub><mrow><mi>r</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445877&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445874&type= 17.01799965 9.56733322 ， 其中 <math id="M122"><mi>i</mi><mo>=</mo><mi>p</mi><mo>,</mo><mi>f</mi><mo>,</mo><mi>e</mi><mo>,</mo><mi>n</mi><mo>,</mo><mi>m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445887&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445883&type= 22.35199928 3.64066648 ，对应轮廓、前额、眼睛、鼻子和嘴， <math id="M123"><msub><mrow><mi>r</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445895&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9445891&type= 2.28600001 3.72533321 表示对应区域模型的正确识别率。因此，融合权重可以根据 式（22） 式（22） 在实际情况中具体设定。 表4 表4 显示了本实验中选择不同方差贡献率进行颅骨识别时各区域的融合权重。基于得到的融合权重，通过区域融合识别颅骨。在测试阶段，第一次验证时，51个测试颅骨中出现错误识别个数为2，识别准确率为96.0%；第二次验证时，51个测试颅骨中出现错误识别个数为1，识别准确率为98.0%；第三次验证时，51个测试颅骨中出现错误识别个数为3，识别准确率为94.1%；第四次验证时，51个测试颅骨中出现错误识别个数为2，识别准确率为96.0%；第五次验证时，51个测试颅骨中出现错误识别个数为4，识别准确率为92.1%。因此，通过融合5个区域建立的相关性分析模型得到，识别准确率为（96.0%+98.0%+94.1%+96.0%+92.1%）/5=95.2%。 区域 轮廓 前额 眼睛 鼻子 嘴巴 权重 0.34 0.28 0.19 0.07 0.09 为了验证本文方法的优越性，我们分别与文献［ 16 16 ］和文献［ 23 23 ］中的方法进行对比。3种方法的对比结果如 表5 表5 所示。 方法 文献[ 23 23 ] 文献[ 16 16 ] 本文 识别准确率/% 93.3 94.0 95.2 从 表5 表5 可以看出，本文方法的识别准确率最高，文献［ 16 16 ］方法次之，文献［ 23 23 ］方法的识别准确率最低。文献［ 23 23 ］中使用嵌入式隐马尔可夫模型，虽然充分利用了颅骨的形状和整体结构信息，但是使用颅骨二维图像，丢失了颅骨的空间信息；而文献［ 16 16 ］和本文均使用颅骨三维模型，并且两者都是利用颅骨和面皮之间的形态关系进行识别。文献［ 16 16 ］首次提出通过分析颅骨和面皮之间的形态相关度来识别未知颅骨，并使用典型关联分析提取映射后两组相关性最大的主成分向量，用于匹配颅骨和面皮。与文献［ 16 16 ］相比，本文采用LSCDA方法分析颅骨和面皮的相关性，通过两个对变量投影的LSMI测量相关性，并通过线性变换使相关性达到最大。该方法不仅能通过线性投影提取两个变量的最相关信息，还可以消除无关信息，因此识别准确率更高。 7 结　论 1 本文提出了一种基于三维颅骨与三维人脸形态相关性度量的未知颅骨识别技术。使用三维颅骨数据作为探针，三维人脸数据作为匹配库，通过探针和匹配库之间的相关性度量，将未知颅骨与已配准的三维人脸进行匹配，从而实现颅骨身份识别。首先，构建颅骨和面皮的统计形状模型，将高维颅骨和面皮映射到低维形状参数空间；然后，基于LSCDA提取颅骨和面皮的主相关信息，构建整体相关性分析模型；最后，计算颅骨与面皮库中每个面皮之间的匹配分数，具有最高匹配分数的面皮即为正确的识别结果。为了提高匹配的准确性，将颅骨和面皮划分为5个生理特征区域，分别为前额、眼睛、鼻子、嘴巴和轮廓，建立5个区域的相关性分析模型，并通过区域融合进行颅骨识别。与颅像重合和颅面复原相比，本文方法不需要颅骨和面皮之间的精确关系，就可以实现未知颅骨身份的确定。本研究证明了颅骨和面皮之间确实存在着较为紧密的联系，为颅像重合和颅面复原的研究提供了理论支持。此外，在提取颅骨和面皮之间的关系方面，基于区域的策略要优于整体策略，这也为基于区域的颅面复原研究提供了参考。 随着基于图像的三维人脸建模技术和三维数据采集技术的发展，三维颅面模型获取变得越来越容易，扩展三维颅面数据集从而进一步提升本文方法的识别率，使它能够更为广泛地用于实际案例中，是下一步的研究重点。 参考文献 BRINKMANN B . Forensic anthropology ［J］. International Journal of Legal Medicine ， 2007 ， 121 （ 6 ）： 431 - 432 . Forensic anthropology DAMAS S ， CORDÓNOSCAR ， IBÁÑEZÓSCAR ， et al . . Forensic identification by computer-aided craniofacial superimposition： A survey ［J］. ACM Computing Surveys ， 2011 ， 43 （ 4 ）： 1 - 27 . . Forensic identification by computer-aided craniofacial superimposition： A survey HU Y ， DUAN F ， YIN B ， et al . . A hierarchical dense deformable model for 3D face reconstruction from skull ［J］. Multimedia Tools & Applications ， 2013 ， 64 （ 2 ）： 345 - 364 . . A hierarchical dense deformable model for 3D face reconstruction from skull BAI F ， CHENG Z ， QIAO X ， et al .. Face reconstruction from skull based on least squares canonical dependency analysis ［C］. IEEE International Conference on Systems . IEEE ， 2016 ： 3612 - 3615 . Face reconstruction from skull based on least squares canonical dependency analysis LUO L ， WANG MY ， TIAN Y ， et al . . Automatic sex determination of skulls based on a statistical shape model ［J］. Computational & Mathematical Methods in Medicine ， 2013 ： 251628 . . Automatic sex determination of skulls based on a statistical shape model LI Y ， CHANG L ， QIAO X ， et al .. Craniofacial reconstruction based on least square support vector regression ［C］. IEEE International Conference on Systems ， 2014 ： 1147 - 1151 . Craniofacial reconstruction based on least square support vector regression FENTON T W ， HEARD A N ， SAUER N J . Skull-photo superimposition and border deaths： identification through exclusion and the failure to exclude ［J］. Journal of Forensic Sciences ， 2008 ， 53 （ 1 ）： 34 - 40 . Skull-photo superimposition and border deaths： identification through exclusion and the failure to exclude IBÁNEZ O ， VALSECCHI A ， CAVALLI F ， et al . . Study on the criteria for assessing skull-face correspondence in craniofacial superimposition ［J］. Legal Medicine ， 2016 ， 23 ： 59 - 70 . . Study on the criteria for assessing skull-face correspondence in craniofacial superimposition RICCI A ， MARELLA G L ， APOSTOL M A . A new experimental approach to computer-aided face/skull identification in forensic anthropology ［J］. Am J Forensic Med Pathol ， 2006 ， 27 （ 1 ）： 46 - 49 . A new experimental approach to computer-aided face/skull identification in forensic anthropology IBÁNEZ O ， CORDÓN O ， DAMAS S . A cooperative coevolutionary approach dealing with the skull-face overlay uncertainty in forensic identification by craniofacial superimposition ［J］. Soft Computing ， 2012 ， 16 （ 5 ）： 797 - 808 . A cooperative coevolutionary approach dealing with the skull-face overlay uncertainty in forensic identification by craniofacial superimposition BALLERINI L ， CORDON O ， DAMAS S ， et al .. Craniofacial superimposition in forensic identification using genetic algorithms ［C］. International Symposium on Information Assurance & Security . IEEE ， 2007 ： 429 - 434 . Craniofacial superimposition in forensic identification using genetic algorithms 田中金 . 基于分区变形颅面复原算法的研究与实现 ［D］. 西安 ： 西北大学 ， 2011 . 基于分区变形颅面复原算法的研究与实现 TIAN ZH J . Research and Implementation of Craniofacial Reconstruction Algorithm Based on Partition Deformation ［D］. Xi'an ： Northwestern University ， 2011 . （in Chinese） Research and Implementation of Craniofacial Reconstruction Algorithm Based on Partition Deformation DUAN F Q ， YANG S ， HUANG D H ， et al . . Craniofacial reconstruction based on multi-linear subspace analysis ［J］. Multimedia Tools & Applications ， 2014 ， 73 （ 2 ）： 809 - 823 . . Craniofacial reconstruction based on multi-linear subspace analysis SONG M ， TAO D ， HUANG X ， et al . . Three-dimensional face reconstruction from a single image by a coupled RBF network ［J］. IEEE Trans Image Process ， 2012 ， 21 （ 5 ）： 2887 - 2897 . . Three-dimensional face reconstruction from a single image by a coupled RBF network 王俊娟 ， 许磊 ， 黎智辉 ， 等 . 三维人脸图像的数据采集与预处理 ［J］. 刑事技术 ， 2015 ， 40 （ 2 ）： 98 - 101 . 三维人脸图像的数据采集与预处理 WANG J J ， XU L ， LI ZH H ， et al . . Acquisition and pre-processing of 3D facial image ［J］. Forensic Science and Technology ， 2015 ， 40 （ 2 ）： 98 - 101 . （in Chinese） . Acquisition and pre-processing of 3D facial image DUAN F ， YANG Y ， LI Y ， et al . . Skull identification via correlation measure between skull and face shape ［J］. Information Forensics & Security IEEE Transactions on ， 2014 ， 9 （ 8 ）： 1322 - 1332 . . Skull identification via correlation measure between skull and face shape ROSHNARA N N P P ， JABBAR S . Efficient 3D visual hull reconstruction based on marching cube algorithm ［C］. International Conference on Innovations in Information . IEEE ， 2015 ： 1 - 6 . Efficient 3D visual hull reconstruction based on marching cube algorithm 朱艳芳 ， 刘晓宁 ， 耿国华 . 基于体素模型与多尺度约束的三维面皮点对应 ［J］. 西北大学学报：自然科学版 ， 2017 ， 47 （ 1 ）： 43 - 49 . 基于体素模型与多尺度约束的三维面皮点对应 ZHU Y F ， LIU X N ， GENG G H . A correspondence method of 3D facial points based on voxel models and multi-scale constraints ［J］. Journal of Northwest University： Natural Science Edition ， 2017 ， 47 （ 1 ）： 43 - 49 . （in Chinese） A correspondence method of 3D facial points based on voxel models and multi-scale constraints DENG Q Q ， ZHOU M Q ， WU Z K ， et al . . A regional method for craniofacial reconstruction based on coordinate adjustments and a new fusion strategy ［J］. Forensic Science International ， 2016 ， 259 ： 19 - 31 . . A regional method for craniofacial reconstruction based on coordinate adjustments and a new fusion strategy DENG Q Q ， ZHOU M Q ， WU Z K . An automatic non-rigid registration method for dense surface models ［C］. IEEE International Conference on Intelligent Computing and Intelligent Systems . IEEE ， 2010 ： 888 - 892 . An automatic non-rigid registration method for dense surface models 贺毅岳 ， 马自萍 ， 高妮 ， 等 . 基于偏最小二乘回归局部形状关系建模的颅面复原方法 ［J］. 计算机应用 ， 2016 ， 36 （ 3 ）： 820 - 826 . 基于偏最小二乘回归局部形状关系建模的颅面复原方法 HE Y Y ， MA Z P ， GAO N ， et al . . Craniofacial reconstruction method based on partial least squares regression model of local craniofacial morphological correlation ［J］. Journal of Computer Applications ， 2016 ， 36 （ 3 ）： 820 - 826 . （in Chinese） . Craniofacial reconstruction method based on partial least squares regression model of local craniofacial morphological correlation KARASUYAM M ， SUGIYAMA M . Canonical dependency analysis based on squared-loss mutual information ［J］. Neural Networks ， 2012 ， 34 （ 10 ）： 46 - 55 . Canonical dependency analysis based on squared-loss mutual information ZHANG Y ， XU S ， SUN N . Skull recognition based on embedded hidden markov models ［C］. World Congress on Intelligent Control & Automation . IEEE ， 2008 ： 1153 - 1158 . Skull recognition based on embedded hidden markov models 10.37188/OPE.20212901.0201 TP394.1 TH691.9 A 1004-924X（2021）01-0201-10 周明全,杨稳,林芃樾等.基于最小二乘正则相关性分析的颅骨身份识别[J].光学精密工程,2021,29(01):201-210. 周明全,杨稳,林芃樾等.基于最小二乘正则相关性分析的颅骨身份识别[J].光学精密工程,2021,29(01):201-210. ZHOU Ming-quan,YANG Wen,LIN Peng-yue,et al.Skull identification based on least square canonical correlation analysis[J].Optics and Precision Engineering,2021,29(01):201-210. ZHOU Ming-quan,YANG Wen,LIN Peng-yue,et al.Skull identification based on least square canonical correlation analysis[J].Optics and Precision Engineering,2021,29(01):201-210. 国家自然科学基金资助项目(No.61731015);国家重点研发项目(No.2019YFC1521103);陕西省重点研发计划资助项目(No.2019ZDLGY10-01,No.2019ZDLSF07-02);青海省重点研发计划资助项目(No.2020-SF-142);陕西省教育厅自然科学专项(No.18JK0795,No.19JK0842) 2020-08-03 2020-09-30 2021-01-15 2021-03-25 光学精密工程 1 29