基于高斯混合模型和轮廓描述子的断骨模型配准方法技术

本发明专利技术公开了一种基于高斯混合模型和轮廓描述子的断骨模型配准方法，包括以下步骤：S1：利用高斯混合模型对低曲率点进行聚类；S2：对每个聚类进行椭圆拟合，并根据椭圆参数提取断面点集；S3：根据断面轮廓构造轮廓描述子；S4：利用卷积自编码器对轮廓描述子进行降维，得到几何特征向量；S5：根据几何向量提取匹配点，并筛选得到基准点，然后根据基准点进行断骨模型配准。本方法实现了骨骼断面的精确分割，可以将碎裂为多块的骨骼模型进行自动精确配准，可以大幅度提高骨折手术的治疗效果。

【技术实现步骤摘要】
基于高斯混合模型和轮廓描述子的断骨模型配准方法
本专利技术涉及断骨匹配
，尤其涉及一种基于高斯混合模型和轮廓描述子的断骨模型配准方法。
技术介绍
目前骨折手术缺少良好的术前规划，需要医生在手术过程中凭经验手工制作钢板，这种手术方法时间长、出血多、容易造成感染，并且手工制作的钢板不够精确。因此，我们利用计算机对断骨的三维模型进行配准，然后根据配准后的模型进行钢板拟合，从而在手术前获取所需钢板，可以大幅度缩短手术时间，提高治疗效果；现有技术仅能实现断裂为两部分的骨骼的配准，而不能将碎裂为多部分的骨骼进行精确配准。
技术实现思路
根据现有技术存在的问题，本专利技术公开了一种基于高斯混合模型和轮廓描述子的断骨模型配准方法，具体包括如下步骤：S1:提取断骨模型中低曲率的点，利用高斯混合模型对低曲率的点进行聚类；S2:对每个聚类进行椭圆拟合、得到该聚类对应的长轴和短轴，根据每个聚类的长轴和短轴的长度找到互相匹配的聚类，将聚类中包含的点集定义为断面点集；S3:以断面点集中的每个点为起始点、与该断面点集的轮廓点之间进行均匀连线，将所有线段按照长短顺序保存在一个向量中，将该向量定义为轮廓描述子；S4:利用卷积自编码器对断面点集中每个点的轮廓描述子进行降维，其中降维后得到的向量为该点的几何特征向量；S5:根据几何特征向量提取多组匹配点、并利用遗传算法对匹配点进行筛选得到基准点，根据获取的基准点采用RANSAC算法对断骨模型进行配准。S1具体采用如下方式：S11：计算断骨模型中每个点的曲率，将曲率小于1.5的点提取出来作为低曲率的点；S12：计算断骨模型的体积V，然后利用公式nc＝0.0014V+4计算得到聚类的数目nc；S13：建立一个高斯混合模型，并将其聚类数目设置为nc，然后利用该高斯混合模型对低曲率的点进行聚类，共得到nc个聚类点集。S2具体采用如下方式：S21：将每个聚类点集的边缘点提取出来，作为该聚类点集的轮廓；S22：对每个聚类点集的轮廓进行椭圆拟合，将椭圆的长轴和短轴的长度提取出来作为该聚类点集的形状参数；S23：将两个断骨模型的所有聚类点集分别放在两个数组中，然后对两个数组中的聚类点集进行两两比较，分别计算两个聚类点集之间形状参数的差值，差值最小的两个聚类点集就是互相匹配的断面点集；S24：遍历所有的断骨模型，将所有互相匹配的断面点集提取出来。S3具体采用如下方式：S31：计算稳定法向量：对断面点集中的每个点，都采集其周围邻域内的点组成邻域点集，然后利用PCA算法对邻域点集的坐标矩阵进行特征分离，得到3个特征值，其中最小的特征值对应的特征向量即为该点的稳定法向量。计算断面点集中每个点的稳定法向量，同时计算断面点集的整体法向量；S32：提取轮廓控制点：以断面点集的整体法向量为轴，按照逆时针的方向每隔每隔一定角度(该角度值可以是12)选一个轮廓点作为初始种子点，共得N(其中N可以取30)个初始种子点，然后根据这些初始种子点使用K-means算法对轮廓点进行聚类，将聚类后的中心点作为轮廓控制点；S33：根据轮廓控制点筛选轮廓点：轮廓控制点是按照提取的顺序保存在数组中的，第i个控制点的空间坐标是Ci，控制点i处的方向向量vi可以用如下公式计算：vi＝Ci-Ci-1。对每个轮廓控制点i，都用如下公式计算γ：γ＝vi×vi-1，如果γ<0，则将轮廓控制点i剔除；S34：构建轮廓描述子：对断面点集中的每个点，以它的稳定法向量为轴，按照逆时针方向每隔1度向轮廓发射一条射线，在射线周围随机选取三个邻近的轮廓点，然后计算这三个轮廓点与该断面点之间的欧式距离并取平均值，按照这种方法共计算得到360个距离值，将其保存在一个向量中，该向量就是该断面点的轮廓描述子。S4具体采用如下方式：S41：建立卷积自编码器模型：构建12层的网络模型，其中编码器部分和解码器部分各6层，编码器由卷积层和池化层构成，解码器由反卷积层和上采样层构成，该网络模型的损失函数是输出向量和输入向量之间的均方误差；S42：利用卷积自编码器模型进行降维：对断面点集中的每个点，将其轮廓描述子作为卷积自编码器模型的输入，然后将编码器部分的输出提取出来作为该断面点的几何特征向量。S5具体采用如下方式：S51：均匀选取匹配点：在断面点集中进行均匀选点，然后在匹配的断面中找到所选点的匹配点，形成多个匹配点对；S52：利用遗传算法筛选得到基准点：随机抽取两组匹配点对可以构成一个DNA的结构，构造多个含有不同DNA的个体，在进化过程中将DNA中两组匹配点之间相对距离大的个体淘汰掉，最后种群中比例较大的匹配点即为基准点；S53：利用RANSAC算法进行配准：选取多组与基准点相对距离较小的匹配点对，然后利用RANSAC算法根据这些匹配点对进行配准，就可以实现两个断面点集之间的配准。按照这种方法将所有互相匹配的断面点集都进行配准，就可以实现断骨模型的整体配准。采用了上述技术方案，本专利技术提供的一种基于高斯混合模型和轮廓描述子的断骨模型配准方法，本方法利用高斯混合模型对低曲率的点进行聚类，然后对每个聚类进行椭圆拟合，并根据椭圆的参数找到互相匹配的断面点集；根据断面的轮廓构建轮廓描述子，然后利用卷积自编码器进行降维，得到点的几何特征向量；最后利用RANSAC算法根据几何特征向量进行配准；因此该方法可以将碎裂为多块的骨骼进行精确配准，在医疗领域中可以大幅度提升骨折手术的治疗效果。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术方法的流程图；图2为本专利技术中利用高斯混合模型进行聚类的结果图；图3为本专利技术中提取得到的两个断面点集示意图；图4为本专利技术中构造轮廓描述子的示意图；图5为本专利技术中卷积自编码器的结构图；图6为本专利技术中断骨模型整体配准的结果图。具体实施方式为使本专利技术的技术方案和优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述：如图1所示的一种基于高斯混合模型和轮廓描述子的断骨模型配准方法，具体采用如下方式：S1：将断骨模型中低曲率的点提取出来，然后利用高斯混合模型对低曲率的点进行聚类。如图2所示，该步骤的具体过程为：S11：计算断骨模型中每个点的曲率，将曲率小于设定阈值(这里可以取1.5)的点提取出来作为低曲率的点；S12：计算断骨模型的体积V，然后利用公式nc＝0.0014V+4计算得到聚类的数目nc；S13：建立一个高斯混合模型，并将其聚类数目设置为nc，然后利用该高斯混合模型对低曲率的点进行聚类，本文档来自技高网...

【技术保护点】
1.一种基于高斯混合模型和轮廓描述子的断骨模型配准方法，其特征在于包括：/n提取断骨模型中低曲率的点，利用高斯混合模型对低曲率的点进行聚类；/n对每个聚类进行椭圆拟合，得到该聚类对应的长轴和短轴，根据每个聚类的长轴和短轴的长度找到互相匹配的聚类，将聚类中包含的点集定义为断面点集；/n以断面点集中的每个点为起始点，与该断面点集的轮廓点之间进行均匀连线，将所有线段按照长短顺序保存在一个向量中，将该向量定义为轮廓描述子；/n利用卷积自编码器对断面点集中每个点的轮廓描述子进行降维，其中降维后得到的向量为该点的几何特征向量；/n根据几何特征向量提取多组匹配点，并利用遗传算法对匹配点进行筛选得到基准点，根据获取的基准点采用RANSAC算法对断骨模型进行整体配准。/n

【技术特征摘要】
1.一种基于高斯混合模型和轮廓描述子的断骨模型配准方法，其特征在于包括：
提取断骨模型中低曲率的点，利用高斯混合模型对低曲率的点进行聚类；
对每个聚类进行椭圆拟合，得到该聚类对应的长轴和短轴，根据每个聚类的长轴和短轴的长度找到互相匹配的聚类，将聚类中包含的点集定义为断面点集；
以断面点集中的每个点为起始点，与该断面点集的轮廓点之间进行均匀连线，将所有线段按照长短顺序保存在一个向量中，将该向量定义为轮廓描述子；
利用卷积自编码器对断面点集中每个点的轮廓描述子进行降维，其中降维后得到的向量为该点的几何特征向量；
根据几何特征向量提取多组匹配点，并利用遗传算法对匹配点进行筛选得到基准点，根据获取的基准点采用RANSAC算法对断骨模型进行整体配准。


2.根据权利要求1所述的方法，其特征还在于：所述对低曲率的点进行聚类时：
首先计算断骨模型中每个点的曲率，将曲率小于设定阈值的点设定为低曲率的点；计算断骨模型的体积V，采用公式nc＝0.0014V+4获得聚类的数目nc；建立高斯混合模型，并将其聚类数目设置为nc，利用该高斯混合模型对低曲率的点进行聚类，共得到nc个聚类点集。


3.根据权利要求1所述的方法，其特征还在于：所述断面点集采用如下方式获取：
将每个聚类点集的边缘点提取出来作为该聚类点集的轮廓；对每个聚类点集的轮廓进行椭圆拟合，将椭圆的长轴和短轴的长度提取出来作为该聚类点集的形状参数；将两个断骨模型的所有聚类点集分别放在两个数组中，再对两个数组中的聚类点集进行两两比较，分别计算两个聚类点集之间形状参数的差值，其中差值最小的两个聚类点集为互相匹配的断面点集；遍历所有的断骨模型、提取所有互相匹配的断面点集。


4.根据权利要求1所述的方法，其特征还在于：获取所述轮廓描述子时：计算稳定法向量
对断面点集中的每个点采集其周围邻域内的点组成邻域点集，利用PCA算法对邻域点集的坐标矩阵进行特征分离得到3个特征值，其中最小的特征值对应的特征向量为该点的稳定法向量，计算断面点集中每个点的稳定法向量，同...

【专利技术属性】
技术研发人员：刘斌，张松，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21