基于最短路径方式的骨骼模型轮廓曲线的平铺生成方法技术

提供一种基于最短路径方式的骨骼模型轮廓曲线的平铺生成方法与系统，用于通过一组平行的骨骼切片生成三维曲面模型，包括：扫描图像分割步骤，其用于对骨骼的扫描图像进行分割操作，获得多组平行的骨骼切片图像；切片轮廓点提取步骤，其基于每层骨骼切片图像提取切片轮廓及各切片轮廓上的点；归一化处理步骤，其能够对任一层骨骼切片图像的切片轮廓基于该切片轮廓的中心坐标位置与半径进行归一化处理；和基于最小面积法的平铺处理步骤，其能够在归一化处理步骤之后，搜索每一对相邻的切片轮廓之间的两个最近点作为三角剖分的起点和图形的第一个节点P（0,0），并确定基于一对相邻的切片轮廓上的对应点的所有节点P（i,j）之间的最小路径。的最小路径。的最小路径。

【技术实现步骤摘要】
基于最短路径方式的骨骼模型轮廓曲线的平铺生成方法


[0001]本专利技术涉及医疗图像处理技术，尤其涉及一种基于最短路径方式的骨骼模型轮廓曲线的平铺生成方法与相应的手术导航系统。

技术介绍

[0002]为了提高手术的准确性和精确度，计算机辅助手术导航系统被广泛应用于各种手术中。这些手术通常需要在实施过程中对患者的解剖结构进行配准，例如全膝关节置换术（TKA）和其他骨科手术，必须准确生成骨骼模型。
[0003]为了生成此模型，须对CT扫描图像进行分割操作，以分割出来一组平行的骨骼切片，用于后续的骨骼重建。通常，这些切片的密度为1
‑
2 mm。
[0004]通过这样的一组平行切片生成三维模型或曲面（表面，surface）（以下简称三维曲面模型）的过程是困难的，多数科学文献中均有讨论。
[0005]其中，通过一组轮廓（平行切片的边缘）生成一个曲面的问题可以分解成若干子问题：例如：
‑
相关问题可以通过确定一组切片的轮廓之间相邻的拓扑关系来解决。从而得到最终曲面粗略的拓扑结果。
[0006]‑
平铺（tiling）问题可以通过相邻切片的边缘上对应点构成三角形网格，从而找到最优的相邻拓扑关系来解决。通常确定最优的标准是其能否将表面积最小化。
[0007]‑
当一个对象在相邻区域用不同数量的轮廓表示时会发生分支问题。在该情况下，不能直接使用标准方法来解决平铺问题。然而需要解决平铺和分支问题以确定曲面的拓扑结构和粗略几何外形。
[0008]‑
解决曲面拟合问题，可以通过将一个最佳曲面拟合在解决上述问题后得到的网格上来完成。曲面拟合的结果可用来展现重建的曲面几何体的细节。
[0009]平铺问题的主要目标是找到连续轮廓上点之间的最佳对应关系。这个问题可以简化为在环形图上的搜索问题，如E， Keppel的论文“通过轮廓线的三角剖分逼近复杂曲面”（E. keppel， Approximating Complex Surfaces by Triangulation of Contour Lines）。
[0010]轮廓由有序的数据点列表表示。连接同一轮廓中相邻点的边被称为轮廓线段。连接不同轮廓中两点的边被称为跨距。跨距由图中的节点表示。平铺块由图中节点之间的圆弧完成。通过在图中找到最小路径方式来完成平铺。
[0011]然而，平铺问题会受到严重约束，因为许多曲面都可以作为结果。在以往，研究者使用了若干公认的标准函数来区分曲面的好坏，最广为人知的有：最大体积法：仅对于凸型曲面(对象)而言是一个明显的标准，但使用起来较为困难。其需要计算单个平铺块对基于对象整体体积的目标函数的贡献值。由于曲面构建误差范围大并带来了额外的计算开销，因而效率低并限制了应用范围。
[0012]最小面积法：这本来是一个很好的标准，但是，仍必须特别注意如何处理好某些具
有复杂几何形状或间距较远的切片的情况。
[0013]尤其是，对于至少局部具有复杂几何形状的骨骼模型来说，即便在适用基于最小面积的标准函数时，仍存在克服上述局限性的问题。

技术实现思路

[0014]为解决上述平铺问题，本专利技术的目的在于提供一种适用于从一组平行的骨骼切片中生成三维曲面的算法及相应的系统如下所述。
[0015]根据本专利技术的一方面，提供一种基于最短路径方式的骨骼模型轮廓曲线的平铺生成方法，用于通过一组平行的骨骼切片生成三维曲面模型，其包括如下步骤：扫描图像分割步骤，构成为用于对骨骼的扫描图像进行分割操作，获得多组平行的骨骼切片图像；切片轮廓点提取步骤，构成为基于每层骨骼切片图像提取切片轮廓及各切片轮廓上的点；归一化处理步骤，构成为能够对任一层骨骼切片图像的切片轮廓基于该切片轮廓的中心坐标位置与半径进行归一化处理；和基于最小面积法的平铺处理步骤，构成为能够在归一化处理步骤之后，搜索每一对相邻的切片轮廓之间的两个最近点作为三角剖分的起点和图形的第一个节点，并确定基于一对相邻的切片轮廓上的对应点的所有节点之间的最小路径。
[0016]优选地，通过对骨骼的CT扫描图像进行分割操作，获得平行的骨骼切片图像。
[0017]优选地，在归一化处理步骤中，针对所有层的骨骼切片图像的切片轮廓进行归一化处理；或者通过基于神经网络的分类器或预定阈值来识别需要进行归一化处理的变形程度大的骨骼切片图像，并仅对变形程度大的骨骼切片图像进行归一化处理。
[0018]优选地，在归一化处理步骤中，还包括如下子步骤：S01：计算切片轮廓的中心点坐标；S02：通过将该中心点坐标移动到坐标原点（0，0）作为中心点而移动切片轮廓；S03：定义包围切片轮廓的矩形窗口；S04：基于该矩形窗口的长度和宽度，对切片轮廓的坐标基于预设的比例进行缩放。
[0019]本专利技术还提供一种计算机辅助手术系统，包括用于生成骨骼模型的计算机，该计算机能够执行程序以实现权利要求１至5中任一项的方法的步骤。
[0020]本专利技术还提供一种骨骼手术系统，包括：扫描图像分割单元，构成为用于对骨骼的扫描图像进行分割操作，获得多组平行的骨骼切片图像；切片轮廓点提取单元，构成为基于每层骨骼切片图像提取切片轮廓及各切片轮廓上的点；归一化处理单元，构成为能够对任一层骨骼切片图像的切片轮廓基于该切片轮廓的中心坐标位置与半径进行归一化处理；和基于最小面积法的平铺处理单元，构成为能够搜索每对相邻的切片轮廓之间的两个最近点作为三角剖分的起点和图形的第一个节点，并确定该对相邻的切片轮廓上的所有对应点之间的最小路径。
[0021]在此该半径是指该轮廓的包络半径，尤其是能够包络归一化之后的轮廓的半径，也可以采用与包围该切片轮廓的矩形窗口的长度和宽度相对应的该轮廓的长和高。在此所谓的包围是指该切片轮廓与该矩形窗口的各侧边之间存在至少一个相切点。
[0022]根据本专利技术还提供一种计算机可读存储介质，存储有计算机程序，计算机程序用于上述方法的步骤。
[0023]本专利技术还提供一种电子设备，包括：处理器；用于存储处理器可执行指令的存储器；处理器，用于从存储器中读取可执行指令，并执行指令以实现上述方法的步骤。
[0024]本专利技术还提供一种计算机程序产品，包括计算机程序，其特征在于，该程序被计算机执行时上述方法的步骤。
附图说明
[0025]图1示意性地示出根据最小路径方式完成平铺的原理图。
[0026]图2示意性地示出根据本专利技术的归一化的模式图。
[0027]图3示意性地示出基于向量的三角形面积算法。
[0028]图4示意性地示出根据本专利技术实施例的最短路径算法原理图。
具体实施方式
[0029]下面结合附图详细描述本专利技术的示例性实施例。下文描述的和附图示出的示例性实施例旨在教导本专利技术的原理，使本领域技术人员能够在若干不同环境中和对于若干不同应用实施和使用本专利技术。因此，本专利技术的保护范围由所附的权利要求来限定，示例性实施例并不意在、并且不应该被认为是对本本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于最短路径方式的骨骼模型轮廓曲线的平铺生成方法，用于通过一组平行的骨骼切片生成三维曲面模型，其特征在于，包括如下步骤：扫描图像分割步骤，构成为用于对骨骼的扫描图像进行分割操作，获得多组平行的骨骼切片图像；切片轮廓点提取步骤，构成为基于每层骨骼切片图像提取切片轮廓及各切片轮廓上的点；归一化处理步骤，构成为能够对任一层骨骼切片图像的切片轮廓基于该切片轮廓的中心坐标位置与半径进行归一化处理；和基于最小面积法的平铺处理步骤，构成为能够在所述归一化处理步骤之后，搜索每一对相邻的切片轮廓之间的两个最近点作为三角剖分的起点和平铺图形的第一个节点P（0，0），并确定基于所述一对相邻的切片轮廓上的对应点的所有节点之间的最小路径。2.根据权利要求1所述的方法，其特征在于，通过对骨骼的CT扫描图像进行分割操作，获得所述平行的骨骼切片图像。3.根据权利要求1所述的方法，其特征在于，在所述归一化处理步骤中，针对所有层的骨骼切片图像的切片轮廓进行归一化处理。4.根据权利要求1所述的方法，其特征在于，通过基于神经网络的分类器或预定阈值来识别需要进行归一化处理的变形程度大的骨骼切片图像，并仅对所述变形程度大的骨骼切片图像进行归一化处理。5.根据权利要求1所述的方法，其特征在于，在所述归一化处理步骤中，还包括如下子步骤：S01：计算所述切片轮廓的中心点坐标；S02：通过将该中心点坐标移动到坐标原点（0，0）作为中心点而移动所述切片轮廓；S03：定义包围所述切片轮廓的矩形窗口；S04：基于该矩形窗口的长度和宽度作为所述半径，对所述切片轮廓的坐标按如下比例缩放：，式（1），其中，X...

【专利技术属性】
技术研发人员：M，
申请(专利权)人：杭州素问九州医疗科技有限公司，
类型：发明
国别省市：