基于左右对称信息的肩关节缺损参数统计方法和系统的制作方法
【技术领域】
[0001] 本发明涉及医学技术领域，具体涉及一种基于左右对称信息的肩关节缺损参数统 计方法和系统。
【背景技术】
[0002] 肩关节是人体活动度最大的关节，与上肢的日常运动密切相关。肩关节不稳定指 在活动时肱骨头在肩关节盂窝内存在有症状的移位。肩关节不稳定是临床常见的肩部疾 病，发生在前下关节盂或者后上部肱骨头的骨缺损都可以导致肩关节不稳定，其会严重影 响患者的生活及工作，降低患者的生活质量。目前治疗骨缺损的有效方法主要有人工关节 置换术和骨移植术。CT图像能够比较清晰的显示骨骼解构，可以提高计算肩关节骨缺损的 体积和深度等参数的精度。因此磨损参数的计算可为外科医生提供可靠依据，是进行骨移 植前最重要的工作之一。
【发明内容】
[0003] 本发明实施例所要解决的技术问题在于，提供一种基于左右对称信息的肩关节缺 损参数统计方法和系统，能够完成肩关节缺损部分信息的精确统计。
[0004] 为解决上述技术问题，本发明实施例提供了一种基于左右对称信息的肩关节缺损 参数统计方法，包括：
[0005] 分割得到CT图像中左右两部分的肩关节表面点云信息；
[0006] 对两个肩关节表面点云进行三角剖分；
[0007] 提取位于两个肩关节表面的四面体；
[0008] 寻找两个肩关节最优变换矩阵；
[0009] 根据最优变换矩阵，配准两个肩关节点云；
[0010] 提取肩关节缺损部分点云信息。
[0011] 在本发明的其他方案中，分割得到的CT图像中的肩关节表面点云信息是肩关节 表面的点云，是离散的空间点的集合。
[0012] 在本发明的其他方案中，对两个肩关节表面点云进行三角剖分以获得构成由肩关 节表面点云构成的四面体集合。
[0013] 在本发明的其他方案中，提取位于两个肩关节表面的四面体的步骤中，该四面体 不包含位于三角剖分后肩关节空间结构内部的点。
[0014] 在本发明的其他方案中，最优变换矩阵中存在7个参数值。
[0015] 在本发明的其他方案中，根据最优变换矩阵，改变一个肩关节点云的空间位置，与 另外一个肩关节点云融合显不。
[0016] 在本发明的其他方案中，根据融合的两个肩关节点云信息，分离两个肩关节点云 互不包含点云，该部分点云为肩关节缺损部分，统计该部分点云信息。
[0017] 本发明实施例还提供了一种基于左右对称信息的肩关节缺损参数统计系统，包 括：
[0018] 肩关节分割模块，其分割得到CT图像中左右两部分的肩关节表面点云信息；
[0019] 三角剖分模块，其用于对两个肩关节表面点云进行三角剖分；
[0020] 表面信息提取模块，其用于提取位于两个肩关节表面的四面体；
[0021] 变换矩阵计算模块，其用于寻找两个肩关节最优变换矩阵；
[0022] 点云配准模块，其根据最优变换矩阵，配准两个肩关节点云；
[0023] 缺损信息统计模块，其提取肩关节缺损部分点云信息。
[0024] 本发明能够快速精准的统计得到肩关节缺损部分的信息，为临床病灶诊疗提供理 论指导。
【附图说明】
[0025] 图1是本发明所提出的基于左右对称信息的肩关节缺损参数统计方法的工作流 程图；
[0026] 图2是本发明所提出的点云三角剖分的示意图；
[0027] 图3是本发明所提出的两个肩关节点云融合的示意图；
[0028] 图4是本发明所提出的基于左右对称信息的肩关节缺损参数统计系统的结构框 架图。
【具体实施方式】
[0029] 下面结合附图和具体实施例对本发明作进一步详细描述，但不作为对本发明的限 定。
[0030] 图1为本发明实施例的基于左右对称信息的肩关节缺损参数统计方法的流程示 意图，具体步骤包括：
[0031] 步骤 S1，
[0032] 分割得到CT图像中左右两部分的肩关节表面点云信息。
[0033] 从CT图像中分割得到左右两部分的肩关节，分割的结果为两个肩关节表面的点 云，是由两个肩关节表面的空间中的点组成的点集。
[0034] 步骤 S2，
[0035] 对两个肩关节表面点云进行三角剖分。
[0036] 三角剖分是指把一个离散点集合剖分成不均匀的三角形网格，如图2所示，示意 性地示出了三角剖分的结果。对于二维实数域内的有限点集N，将点集中的点作为端点构成 封闭线段记为边e，E为边e的集合。该点集N的三角剖分T = (N，E)是一个平面图，该平 面图满足以下条件：
[0037] (1)该平面图中的边彼此不相交；
[0038] (2)平面图的边中，除去端点，不包含点集N中的任何点；
[0039] (3)该平面图中只含有三角形，同时所有三角形的合集是离散点集N的凸包。
[0040] 若集合E中的边e (假设短点为a，b)满足条件：存在一个圆经过a，b两点，且圆内 不含点集N中的任何点，则该三角剖分T称为Delaunay三角剖分[18]。
[0041 ] Delaunay三角剖分的主要特性：
[0042] (I)Delaunay三角剖分是唯一的。
[0043] (2)空圆特性。在点集N的Delaunay三角剖分中，任意一个三角形的外接圆内不 含有N中的任何点。
[0044] (3)最小角最大化性质。在点集N的所有三角剖分中，Delaunay三角剖分将最大 化所有三角形的最小角。
[0045] 步骤 S3，
[0046] 提取位于两个肩关节表面的四面体。
[0047] 对于两个肩关节点云的三角剖分的结果，是由两个肩关节点云中的四个点组成的 四面体的集合，提取肩关节三角剖分结果的表面的四面体，其他四面体位于表面四面体组 成的空间的内部。
[0048] 在提取位于两个肩关节表面的四面体的步骤中，该四面体可以不包含位于三角剖 分后肩关节空间结构内部的点
[0049] 步骤 S4，
[0050] 寻找两个肩关节最优变换矩阵。
[0051] 给定任意初始位置的两组点云，分别为原始点云
与目标点云
:，寻求最优的刚性变换作用于X，使得点云X与Y最为匹配。其中，N r^Ny 为点云X和Y中三维点的数量。最佳变换Tbest的目标是对点云X经变换后的每一个点求取 与点云Y中每个点的最小距离，求所有N1个最小距离的和S i，使S1最小。
[0052] 首先，定义点云X与点云Y的三角剖分结果为D 2，分别提取三角剖分D 2 表面的四面体集合？1与？2。
[0053] 然后，选择四面体集合F1中一个四面体F1，将四面体与四面体集合F 2中所以四 面体进行刚性配准，四面体全等必然存在一下关系：
[0055] 其中I Ia1^1I I表示&1与b i两点之间的欧氏距离，每一组相等判定允许小于误差 限S的误差存在。
[0056] 如果配准的两组点云间存在尺度缩放，我们只需将全等判定改为相似判定，原算 法依然有效。f,与f /相似的判定条件为：
[0058] 其中s为这一对潜在匹配计算出的尺度变换因子，判定同样允许小于δ的误差存 在。
[0059] 步骤 S5，
[0060] 根据最优变换矩阵，配准两个肩关节点云。
[0061 ] 可以根据各种方式配准两个肩关节点云。优选地，可以根据计算得到的最优变换 矩阵，改变一个肩关节点云的空间位置，与另外一个肩关节点云融合。
[0062]步骤 S6，
[0063] 提取肩关节缺损部分点云信息。
[0064] 提取可以由多种方式实现。优选地，可以根据融合的两个肩关节点云信息，分离两 个肩关节点云互不