一种基于主测地线分析的三维下颌骨重建方法技术

本发明专利技术涉及一种基于主测地线分析的三维下颌骨重建方法,包括：步骤1，构建形状空间，将原始的三维头骨去除平移、缩放和旋转的李群作用，投影到形状空间中；步骤2，训练阶段，将头骨训练样本投影到形状空间后，接着在内蕴均值处对所有训练样本进行主测地线分析，最后从训练样本得到上颅骨形状和下颌骨形状之间的关系，用于重建阶段；步骤3，重建阶段，将输入的上颅骨模型投影到形状空间中，根据训练阶段得到的先验知识来构建下颌骨形状空间，最后将其还原到欧氏空间中。本发明专利技术的优点在于使用本发明专利技术所述基于主测地线分析的三维下颌骨重建方法，能有效地重建出与真实值接近的下颌骨。有效地重建出与真实值接近的下颌骨。有效地重建出与真实值接近的下颌骨。

【技术实现步骤摘要】
一种基于主测地线分析的三维下颌骨重建方法


[0001]本专利技术涉及三维下颌骨的重建领域，具体涉及一种基于主测地线分析的三维下颌骨重建方法。

技术介绍

[0002]随着CT和DICOM技术的出现，大量的头骨重建技术涌现出来，传统的手工重建技术逐渐被计算机辅助设计和制造技术替代，比如基于镜像的重建方法，表面插值方法，基于数据驱动的方法等。由于人类颅骨大致是对称的结构，基于镜像的方法在头骨重建的虚拟规划过程中表现出极大的便利。然而，这种方法要求缺陷的部分严格限制在头骨的单侧，不适用于双侧或穿过面部中线的缺陷。表面插值方法对于穿过中线的缺陷是有用的，这种方法确保了缺陷区域与头骨形状的光滑近似。但是对于大的重建区域，该方法缺乏对重建区域的约束，使得重建部分表面过于平坦。另一类基于数据驱动的方法是处理大的缺陷或双侧缺陷的一种通用工具，这种方法利用形变模型从头骨数据库中提取相应的模板，来指导感兴趣区域的重建，比如Marreirors等人提出了一种半自动化的几何形态重建方法，该方法使用薄板样条从数据库中生成模板来指导大面积曲面的重建；Fuessinger等人设计了基于完整颅骨的统计形状模型，将捕获的形状变化作为先验知识来估计缺失部分。数据驱动的方法克服了其它头骨重建方法的前提条件，不受到头骨的对称性或缺陷/孔洞的封闭边界等约束，可以结合模板和系数重建出具有个人特点的头骨形状。
[0003]现有的数据驱动方法绝大部分是直接在欧氏空间中计算头骨模板，忽略了形状本身的特性。欧氏空间中的笛卡尔坐标表示法并没有完全反映出头骨形状本身的信息，因其还包含了与头骨形状无关的其他信息。而Kendall等人提出的形状空间是一种只关注形状本身的内蕴表示法，在形状空间中可以更加精确地刻画头骨形状本身的细节和局部变化特性。
[0004]在现有已经公开的专利文献中，例如,中国专利申请公开号CN111563953A公开了基于机器学习的颌骨缺损重建方法、装置、终端和介质，包括：采集预设人群范围内的多个被采样人员的颌骨CT数据；基于所述颌骨CT数据，确立每个所述被采样人员的上、下颌骨表面特征性的多个颌骨特征点；通过机器学习算法获取各所述颌骨特征点之间的相关性，以基于各所述颌骨特征点之间的相关性来进行颌骨缺损重建。本专利技术利用机器学习进行颌骨特征点的还原，可在复杂颌骨重建过程中提供精确且个性化的方案，解决了跨中线大范围颌骨缺损病例重建仅凭经验、无参照可依的临床难题。此外，本专利技术基于预设人群范围内的颌骨CT数据，有针对性地对特定族群提供颌骨重建策略，在保证骨段血供、保留种植位点的同时，更加贴近该特定族群的面部外形。
[0005]又例如，中国专利公开号CN112785691A公开了一种下颌骨重建方法、装置、电子设备及存储介质，所述下颌骨缺损重建方法通过获取被采样人员的CT数据并进行三维重建，按照预设方案选取下颌骨相关特征点并保存其坐标以建立正常下颌骨数据库；虚拟截骨获得下颌骨缺损模型；根据所述下颌骨相关特征点对所述下颌骨缺损模型标准化处理，并通
过预设的下颌骨数据库检索匹配算法在所述正常下颌骨数据库中检索匹配得到最相似下颌骨；将检索匹配获得的最相似下颌骨与缺损下颌骨模型进行配准。为跨中线下颌骨缺损、二期下颌骨缺损以及双侧面部不对称的下颌骨缺损修复重建的数字化设计提供了除医生经验外的、可重复性高的参考，有效解决了镜像技术在临床运用中的缺陷。
[0006]再例如，中国专利申请公开号CN113397702A公开了一种下颌骨修复体及其制造方法。下颌骨包括颌骨本体，颌骨本体上有缺口，下颌骨修复体包括修复体本体和连接件，修复体本体用于填补颌骨本体上的缺口，修复体本体包括外网格体和内网格体，内网格体套设在外网格体内并连接于外网格体的内壁，修复体本体通过连接件连接于颌骨本体。下颌骨修复体的制造方法包括S1、对颌骨本体进行建模，得到颌骨本体的缺口结构数据；S2、根据缺口的结构数据设置本体模型，本体模型包括外网格体模型和内网格体模型，将内网格体模型套设在外网格体模型内并连接于外网格体模型；S3、将本体模型通过连接件模型连接于颌骨本体。
[0007]上述公开的专利申请均未解决形状空间中可以更加精确地刻画头骨形状本身的细节和局部变化特性，也没有涉及采用主测地线分析的三维下颌骨重建方法。有鉴于此，本专利技术结合数据驱动方法及Kendall等人提出的形状空间，提出一种基于主测地线分析(Principal Geodesic Analysis,PGA)的三维下颌骨重建方法。

技术实现思路

[0008]基于现有技术存在的如上述技术问题，本专利技术提出一种基于主测地线分析的三维下颌骨重建方法，包括以下步骤：
[0009]步骤1，构建形状空间，根据Kendall提出的形状定义，将原始的三维头骨去除平移、缩放和旋转的李群作用，投影到形状空间中；
[0010]步骤2，训练阶段，将头骨训练样本投影到形状空间后，使用Fr
é
chet提出的方法找到形状空间中的内蕴均值，接着在内蕴均值处对所有训练样本进行主测地线分析，最后从训练样本得到上颅骨形状和下颌骨形状之间的关系，作为先验知识，用于重建阶段；
[0011]步骤3，重建阶段，将输入的上颅骨模型投影到形状空间中，根据训练阶段得到的先验知识来构建下颌骨形状空间，最后将其还原到欧氏空间中。
[0012]进一步的，步骤1中，构建形状空间，根据Kendall提出的形状定义，将原始的三维头骨去除平移、缩放和旋转的李群作用，投影到形状空间中，包括以下步骤：
[0013]步骤1.1，三维颅骨形状的表示：在三维空间中的一个由m个点组成的中心化三维颅骨点云数据记为这里的x
h
(h＝1,2,
…
,m)是三维颅骨点云中第h点的三维坐标，利用向量化算符将点云数据转化为一个列向量，即构形空间中的一个元素：
[0014][0015]步骤1.2，投影到形状空间，利用Kendall形状空间理论中的赫尔默特矩阵或中心化矩阵为构形元素去除平移作用，再利用质心尺寸去除缩放作用，得到预形状x，接着对预形状进行普氏分析，去除旋转变换群的作用,得到形状空间中的元素[x]，到形状空间中的元素[x]是预形状x在旋转作用下元素的集合。
[0016]进一步的，步骤2中，将头骨训练样本投影到形状空间后，使用Fr
é
chet提出的方法找到形状空间中的内蕴均值，接着在内蕴均值处对所有训练样本进行主测地线分析，最后从训练样本得到上颅骨形状和下颌骨形状之间的关系，作为先验知识，用于重建阶段，包括以下步骤：
[0017]步骤2.1，形状空间的对数映射，通过在形状空间上对形状进行分析，而鉴于形状空间在其上一点的邻域内与欧氏空间微分同胚，将形状空间的局部映射到欧氏空间，在欧氏空间中完成分析过程，将结果映射回形状空间，由于三维颅骨形状在形状空间中的分布非常集中，能够将它们视为在一个点的邻域内，将三维颅骨形状投影到该点的切空间上，就能够方便地分析形状的本质特征，令p和q是预形状球上的两点本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于主测地线分析的三维下颌骨重建方法,其特征在于，包括以下步骤：步骤1.构建形状空间，根据Kendall提出的形状定义，将原始的三维头骨去除平移、缩放和旋转的李群作用，投影到形状空间中；步骤2.训练阶段，将头骨训练样本投影到形状空间后，使用Fr
é
chet提出的方法找到形状空间中的内蕴均值，接着在内蕴均值处对所有训练样本进行主测地线分析，最后从训练样本得到上颅骨形状和下颌骨形状之间的关系，作为先验知识，用于重建阶段；步骤3.重建阶段，将输入的上颅骨模型投影到形状空间中，根据训练阶段得到的先验知识来构建下颌骨形状空间，最后将其还原到欧氏空间中。2.按照权利要求1所述基于主测地线分析的三维下颌骨重建方法,其特征在于，步骤1构建形状空间，根据Kendall提出的形状定义，将原始的三维头骨去除平移、缩放和旋转的李群作用，投影到形状空间中，包括以下步骤：步骤1.1.三维颅骨形状的表示：在三维空间中的一个由m个点组成的中心化三维颅骨点云数据记为x＝(x
1T
,x
2T
,
…
,x
mT
)，这里的x
h
(h＝1,2,
…
,m)是三维颅骨点云中第h点的三维坐标，利用向量化算符将点云数据转化为一个列向量，即下式(1)所示的构形空间中的一个元素：步骤1.2.投影到形状空间，利用Kendall形状空间理论中的赫尔默特矩阵或中心化矩阵为构形元素去除平移作用，再利用质心尺寸去除缩放作用，得到预形状x，接着对预形状进行普氏分析，去除旋转变换群的作用,得到形状空间中的元素[x]，到形状空间中的元素[x]是预形状x在旋转作用下元素的集合。3.按照权利要求1所述基于主测地线分析的三维下颌骨重建方法,其特征在于，步骤2所述将头骨训练样本投影到形状空间后，使用Fr
é
chet提出的方法找到形状空间中的内蕴均值，接着在内蕴均值处对所有训练样本进行主测地线分析，最后从训练样本得到上颅骨形状和下颌骨形状之间的关系，作为先验知识，用于重建阶段，包括以下步骤：步骤2.1形状空间的对数映射，通过在形状空间上对形状进行分析，而鉴于形状空间在其上一点的邻域内与欧氏空间微分同胚，将形状空间的局部映射到欧氏空间，在欧氏空间中完成分析过程，将结果映射回形状空间，由于三维颅骨形状在形状空间中的分布集中，能够将它们视为在一个点的邻域内，将三维颅骨形状投影到该点的切空间上，就能够分析形状的本质特征，令p和q是预形状球上的两点，[p],[q]分别是这两个预形状的形状表示，借助普氏分析方法，计算[p],[q]在预形状空间中的代表元素p,q
ro
，这里的q
ro
代表q向p进行旋转对齐后的元素，则形状空间中的对数映射表示为下式(2)：上式(2)中，Γ
*
是由普氏分析计算出的q向p对齐的最佳旋转矩阵，此时对数映射是等距表示的，p和q对应形状之间的最短测地距离表示为下式(3)：d([p],[q])＝||Log
p
(q
ro
)||
……
(3)；
步骤2.2求三维颅骨模型的内蕴均值，首先，构造颅骨数据的预形状空间然后使用普氏分析消除旋转作用，将所有预形状x
i
向第一个预形状x
i
对齐，得到颅骨形状在预形状空间中的代表元素，分别表示为x
1,ro
,x
2,ro
,
…
,x
N,ro
，所有颅骨形状组成的形状空间表示为根据Fr
é
chet内蕴均值算法和公式(3)中形状空间的测地距离表示，则给出三维颅骨形状空间中的内蕴均值定义，如下式(4)：上式(4)中，μ
*
表示在形状空间中到所有形状的测地距离平方和最小的点，将形状空间中任意一个颅骨形状x
i，ro
作为内蕴均值的初始值，使用梯度下降法求解上式(4)，三...

【专利技术属性】
技术研发人员：武仲科，严伦，赵海川，王醒策，王学松，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：