基于主向量分析的断骨截面分割及断骨模型配准方法技术

本发明专利技术公开了一种基于主向量分析的断骨截面分割及断骨模型配准方法，包括以下步骤：S1：采用主向量分析算法进行断骨模型的轴线提取；S2：根据断骨三维网格模型上的三角面片法向量突变以及三角面片与断骨轴线夹角提取断骨截面的网格顶点点集；S3：结合断骨轴线对准以及断骨截面网格顶点点集主向量对准的方法对断骨三维网格模型进行空间粗配准；S4：利用迭代最近点算法对两个断骨截面的网格顶点点集进行多次迭代计算，实现断骨三维网格模型的空间精配准；S5：根据断骨三维网格模型配准得到的完整骨骼模型进行骨折钢板模型的拟合。本方法对断骨截面网格顶点点集的精确分割，在粗配准中实现了断面的大致吻合，提升了精配准的精度和成功率。

【技术实现步骤摘要】

本专利技术涉及一种基于主向量分析的断骨截面分割及断骨模型配准方法。
技术介绍
目前骨折手术普遍采用人工复位和伤肢内固定相结合的方法，这种方法存在的问题是创伤大、出血多、容易引发神经血管损伤等并发症。因此，我们可以利用计算机对断骨模型进行虚拟拼接，从而在术前得到钢板的各种几何参数。该技术有望在未来结合3D打印等先进制造方法进行骨折钢板的个性化术前定制，可大大推进骨折手术的微创化进程。
技术实现思路
根据现有技术存在的问题，本专利技术公开了一种基于主向量分析的断骨截面分割及断骨模型配准方法，包括以下步骤：S1：采用主向量分析算法进行断骨模型的轴线提取：对骨折部位的CT序列图像进行重构得到断骨的三维网格模型，断骨的三维网格模型由三维坐标系下一系列的网格顶点点云构成，利用主向量分析算法求得网格顶点点云的主方向向量，将该向量作为断骨的轴线；S2：根据断骨三维网格模型上的三角面片法向量突变和三角面片与断骨轴线的夹角提取断骨截面的网格顶点点集；S3：结合断骨轴线对准以及断骨截面网格顶点点集主向量对准的方法对断骨三维网格模型进行空间粗配准：S4：利用迭代最近点算法对两个断骨截面的网格顶点点集进行多次迭代计算，对断骨三维网格模型的相对位置进行多次微调，直至两个断骨截面网格顶点点集之间的空间距离最小，实现断骨三维网格模型的空间精配准；S5：根据断骨三维网格模型配准得到的完整骨骼模型进行骨折钢板模型的拟合。S1具体采用如下方式：S11：计算断骨三维网格模型中网格顶点点云数据集的中心点坐标：分别计算网格顶点点云数据集中所有顶点的X、Y、Z坐标的平均值，得到中心点坐标；S12：特征中心化：分别计算断骨三维网格模型中每个顶点的坐标与中心点坐标的差值，将计算结果存储在3×n的矩阵A中，其中n为断骨三维网格模型中顶点的总数量；S13：计算协方差矩阵：将矩阵A与它的转置矩阵A′相乘，得到协方差矩阵M；S14：求协方差矩阵M的特征值和特征向量，其中最大的特征值所对应的向量就是该断骨三维网格模型中网格顶点点云数据集的主成分，即断骨三维网格模型轴线的方向向量。S2具体采用如下方式：S21：选择断骨截面上的一个三角面片作为断骨截面顶点点集的种子三角面片；S22：进行三角面片法向量的蔓延分割：从种子三角面片开始，对其相邻的三角面片进行遍历，分别计算种子三角面片的法向量与相邻被遍历到三角面片的法向量之间的夹角余弦cosα，以及被遍历到三角面片的法向量与断骨轴线之间的夹角余弦cosβ；S24：利用公式λ＝2cosα+cosβ来计算被遍历到三角面片的阈值λ，当λ>2.09时，保留该三角面片，否则丢弃该三角面片；S25：选择被保留三角面片作为新的种子三角面片，采用上述方式继续进行蔓延分割，得到断骨截面的网格顶点点集，该点集将作为断骨三维网格模型空间精配准的数据来源。S3具体采用如下方式：S31：断骨轴线对准：通过空间变换矩阵将两个断骨的轴线对齐，使之处于同一直线上；S32：基于断骨截面网格顶点点集主方向对准的断骨截面网格模型对准：通过主向量分析算法分别计算得到两个断骨截面网格顶点点集的主方向，固定一个断骨三维网格模型，将另一个断骨三维网格模型以其轴线为轴进行旋转，直至两个断骨截面网格顶点点集的主方向夹角最小时停止。本专利技术提供了一种基于主向量分析的断骨截面分割及断骨模型配准方法，本方法在断骨截面三维网格模型分割过程中，通过判断三角面片法向量突变及三角面片法向量与断骨三维网格模型轴线的夹角，实现了断骨截面网格顶点点集的精确分割；在断骨三维网格模型预配准的过程中，通过分析断骨截面网格
顶点点集的主方向，实现了断骨截面网格模型的大致吻合，大大提升了精配准的精度和成功率。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术中方法的流程图；图2为本专利技术中两段断骨的轴线提取的效果图；图3(a)为现有技术提取断骨截面的网格顶点点集的效果图；图3(b)为本专利技术中提取断骨截面的网格顶点点集的效果图；图4为已有方法在预配准过程中断面无法实现大致吻合的示意图；图5(a)为断骨截面点集的主方向和断骨模型轴线的示意图；图5(b)为本专利技术在预配准过程中实现断面大致吻合的效果图；图6为在拼接完成的断骨三维网格模型表面选取控制点的示意图；图7为虚拟钢板模型示意图；图8(a)为现有技术缺陷导致断面网格顶点点集配准失败的示意图；图8(b)为现有技术最终的配准效果图；图9(a)为采用本方法实现断面网格顶点点集精确配准的效果图；图9(b)为本方法最终的配准效果图。具体实施方式为使本专利技术的技术方案和优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述：如图1所示的一种基于主向量分析的断骨截面分割及断骨模型配准方法，具体采用如下方式：S1：采用主向量分析算法进行断骨模型的轴线提取：对骨折部位的CT序列图像进行重构得到断骨的三维网格模型，断骨的三维网格模型由三维坐标系下一系列的网格顶点点云构成，利用主向量分析算法求得网格顶点点云的主方向向量，将该向量作为断骨的轴线。如图2所示：该步骤的具体过程为：S11：计算断骨三维网格模型中网格顶点点云数据集的中心点坐标：分别计
算网格顶点点云数据集中所有顶点的X、Y、Z坐标的平均值，得到中心点坐标；S12：特征中心化：分别计算断骨三维网格模型中每个顶点的坐标与中心点坐标的差值，将计算结果存储在3×n的矩阵A中，n为断骨三维网格模型中顶点的总数量；S13：计算协方差矩阵：将矩阵A与它的转置矩阵A′相乘，得到协方差矩阵M；S14：求协方差矩阵M的特征值和特征向量，其中最大的特征值所对应的向量就是该断骨三维网格模型中网格顶点点云数据集的主成分，即断骨三维网格模型轴线的方向向量。现有技术在断骨截面网格顶点点集分割过程中只考虑三角面片法向量突变的因素，而没有考虑三角面片法向量与断骨轴线的夹角，因此当断骨截面较复杂时无法实现断骨截面网格顶点点集的完整提取，如图3(a)所示。本专利技术在提取过程中除三角面片法向量突变外，还引入了三角面片法向量与断骨轴线的夹角这一判断因素，断骨截面网格顶点点集的分割效果大大提升，如图3(b)所示。具体步骤是：S2：根据断骨三维网格模型上的三角面片法向量突变和三角面片与断骨轴线的夹角提取断骨截面的网格顶点点集。S21：选择断骨截面上的一个三角面片作为断骨截面顶点点集的种子三角面片；S22：进行三角面片法向量的蔓延分割：从种子三角面片开始，对其相邻的三角面片进行遍历，分别计算种子三角面片的法向量与相邻被遍历到三角面片的法向量之间的夹角余弦cosα，以及被遍历到三角面片的法向量与断骨轴线之间的夹角余弦cosβ；S24：利用公式λ＝2cosα+cosβ来计算被遍历到三角面片的阈值λ，当λ>2.09时，保留该三角面片，否则丢弃该三角面片；S25：选择被保留三角面片作为新的种子三角面片，采用上述方式继续进行蔓延分割，得到断骨截面的网格顶点点集，该点集将作为断骨三维本文档来自技高网...

【技术保护点】
一种基于主向量分析的断骨截面分割及断骨模型配准方法，其特征在于：包括以下步骤：S1：采用主向量分析算法进行断骨模型的轴线提取：对骨折部位的CT序列图像进行重构得到断骨的三维网格模型，断骨的三维网格模型由三维坐标系下一系列的网格顶点点云构成，利用主向量分析算法求得网格顶点点云的主方向向量，将该向量作为断骨的轴线；S2：根据断骨三维网格模型上的三角面片法向量突变和三角面片与断骨轴线的夹角提取断骨截面的网格顶点点集；S3：结合断骨轴线对准以及断骨截面网格顶点点集主向量对准的方法对断骨三维网格模型进行空间粗配准：S4：利用迭代最近点算法对两个断骨截面的网格顶点点集进行多次迭代计算，对断骨三维网格模型的相对位置进行多次微调，直至两个断骨截面网格顶点点集之间的空间距离最小，实现断骨三维网格模型的空间精配准；S5：根据断骨三维网格模型配准得到的完整骨骼模型进行骨折钢板模型的拟合。

【技术特征摘要】
1.一种基于主向量分析的断骨截面分割及断骨模型配准方法，其特征在于：包括以下步骤：S1：采用主向量分析算法进行断骨模型的轴线提取：对骨折部位的CT序列图像进行重构得到断骨的三维网格模型，断骨的三维网格模型由三维坐标系下一系列的网格顶点点云构成，利用主向量分析算法求得网格顶点点云的主方向向量，将该向量作为断骨的轴线；S2：根据断骨三维网格模型上的三角面片法向量突变和三角面片与断骨轴线的夹角提取断骨截面的网格顶点点集；S3：结合断骨轴线对准以及断骨截面网格顶点点集主向量对准的方法对断骨三维网格模型进行空间粗配准：S4：利用迭代最近点算法对两个断骨截面的网格顶点点集进行多次迭代计算，对断骨三维网格模型的相对位置进行多次微调，直至两个断骨截面网格顶点点集之间的空间距离最小，实现断骨三维网格模型的空间精配准；S5：根据断骨三维网格模型配准得到的完整骨骼模型进行骨折钢板模型的拟合。2.根据权利要求1所述的基于主向量分析的断骨截面分割及断骨模型配准方法，其特征还在于：S1具体采用如下方式：S11：计算断骨三维网格模型中网格顶点点云数据集的中心点坐标：分别计算网格顶点点云数据集中所有顶点的X、Y、Z坐标的平均值，得到中心点坐标；S12：特征中心化：分别计算断骨三维网格模型中每个顶点的坐标与中心点坐标的差值，将计算结果存储在3×n的矩阵A中，其中n为断骨三维网格模型中顶点的总数量；S13：计算协方差矩阵：将矩阵A与它的转置矩阵A′相乘，得到协方差矩阵M；S14：求协方差矩阵...

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