CN112164443B - 一种交互式模拟骨折复位的方法 - Google Patents

一种交互式模拟骨折复位的方法 Download PDF

Info

Publication number
CN112164443B
CN112164443B CN202011031888.7A CN202011031888A CN112164443B CN 112164443 B CN112164443 B CN 112164443B CN 202011031888 A CN202011031888 A CN 202011031888A CN 112164443 B CN112164443 B CN 112164443B
Authority
CN
China
Prior art keywords
reset
fracture
positioning surface
reduction
index
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011031888.7A
Other languages
English (en)
Other versions
CN112164443A (zh
Inventor
王硕
刘建
赵晓萌
张宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Qingtai Simulation Technology Co ltd
Original Assignee
Shanghai Qingtai Simulation Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Qingtai Simulation Technology Co ltd filed Critical Shanghai Qingtai Simulation Technology Co ltd
Priority to CN202011031888.7A priority Critical patent/CN112164443B/zh
Publication of CN112164443A publication Critical patent/CN112164443A/zh
Application granted granted Critical
Publication of CN112164443B publication Critical patent/CN112164443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/30—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30008—Bone

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Biophysics (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Epidemiology (AREA)
  • Software Systems (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Materials For Medical Uses (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明涉及计算机辅助医疗技术领域，更具体地说，涉及一种交互式模拟骨折复位的方法。本发明的方法，包括以下步骤：S1、扫描并构建骨折体三维模型；S2、选择骨折体三维模型中的一部分骨折体为基准骨折体，另一部分骨折体为待复位骨折体；S3、在基准骨折体选取不少于三个基准点生成基准定位面；S4、根据基准定位面的基准点，通过配准算法生成待复位骨折体的对应基准点，确定对应的待复位定位面；S5、根据基准定位面与待复位定位面的相对位置关系，形成复位评价参数，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，实现骨折体三维模型的复位。本发明提供直观的骨科空间方位表示方法，通过简便交互操作协助医生模拟骨折体复位。

Description

一种交互式模拟骨折复位的方法
技术领域
本发明涉及计算机辅助医疗技术领域，更具体的说，涉及一种交互式模拟骨折复位的方法。
背景技术
骨骼的伤病是人体常见疾病，骨科的伤病治疗包括创伤治疗和人造骨体置换。
骨科疾病治疗中经常需要进行骨折复位，每一例骨科治疗案例中，患者的具体的状况差异较大，这些状况包括受伤情况，疾病情况，基础体质，经济条件，器械选择以及年龄和性别。这些具体状况对于每个患者差异很大，因此，骨科治疗通常不存在一般意义上适用所有患者的唯一最佳复位方案，只存在针对不同患者具体状况的个体最优复位方案。
现代医疗在骨科方向有了很多进展，骨科的伤病治疗几乎百分之百需要医学影像作为治疗的先决条件。
一般而言，多数治疗大夫的主要时间用于治疗方案确定和具体诊疗过程，在医学图形图像处理技术上并非专家。治疗手段需要针对患者当前的身体图像进行虚拟骨折复位，虚拟复位能够协助医生预先准确测量骨骼空间位置关系，演练复位情况。
目前，对骨折复位并没有现成的仿真模拟方法，医学影像难以给医生直观有效的协助。
发明内容
本发明的目的是提供一种交互式模拟骨折复位的方法，解决现有技术中医学影像难以直观有效的仿真模拟骨折复位的问题。
为了实现上述目的，本发明提供了一种交互式模拟骨折复位的方法，包括以下步骤：
S1、扫描并构建骨折体三维模型；
S2、选择骨折体三维模型中的一部分骨折体为基准骨折体，另一部分骨折体为待复位骨折体；
S3、在基准骨折体选取不少于三个基准点生成基准定位面；
S4、根据基准定位面的基准点，通过配准算法生成待复位骨折体的对应基准点，确定对应的待复位定位面；
S5、根据基准定位面与待复位定位面的相对位置关系，形成复位评价参数，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，实现骨折体三维模型的复位。
在一实施例中，所述步骤S3，进一步包括以下步骤：
S31、在基准骨折体的骨折面轮廓线上，选取距离较远且不共线的三个基准点；
S32、将基准点连接形成基准定位面，用于绑定表征基准骨折体。
在一实施例中，所述步骤S4的配准算法，为迭代最近点算法，进一步包括以下步骤：
将基准骨折体的骨折面的外轮廓点云集合P，与待复位骨折体的骨折面的外轮廓点云集合G进行最近点距离的最小二乘法迭代计算，得到配准位移矩阵M；
将基准骨折体的三个基准点，通过叉乘配准位移矩阵M，得到待复位骨折体的对应的三个待复位基准点。
在一实施例中，所述步骤S5的复位评价参数，包括第一复位指标，为直线待复位距离：
所述直线待复位距离，为基准定位面和待复位定位面的中心点的连线距离。
在一实施例中，所述步骤S5中，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，进一步包括以下步骤：
通过连接基准定位面和待复位定位面的中心点的向量，使待复位骨折体沿此向量进行移动，调节第一复位指标的大小，直至第一复位指标满足预设阈值。
在一实施例中，所述步骤S5的复位评价参数，包括第一综合复位评价指标MK1，公式如下：
当ID1≥TolA1时，
Figure BDA0002703958680000031
当TolA1＞ID1≥TolB1时，
Figure BDA0002703958680000032
当TolB1＞ID1时，
Figure BDA0002703958680000033
其中，ID1为第一复位指标，TolA1为第一复位指标的合格容差限值，TolB1为第一复位指标的良好容差限值。
在一实施例中，所述步骤S5的复位评价参数，包括第二复位指标，为轴向待复位角度：
所述轴向待复位角度，为基准定位面的法线与待复位定位面的法线形成的夹角。
在一实施例中，所述步骤S5中，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，进一步包括以下步骤：
通过待复位定位面的法线的向量和基准定位面的法线的向量，用叉乘法计算共同的垂向量；
过待复位定位面的中心点做平行于垂向量的平行线；
待复位骨折体以此平行线为旋转轴，调节第二复位指标的大小，直至第二复位指标满足预设阈值。
在一实施例中，所述步骤S5的复位评价参数，包括第二综合复位评价指标MK2，公式如下：
当ID2≥TolA2时，
Figure BDA0002703958680000034
当TolA2＞ID2≥TolB2时，
Figure BDA0002703958680000041
当TolB2＞ID2时，
Figure BDA0002703958680000042
其中，ID2为第二复位指标，TolA2为第二复位指标的合格容差限值，TolB2为第二复位指标的良好容差限值。
在一实施例中，所述步骤S5的复位评价参数，包括第三复位指标，为周向待复位角度：
所述周向待复位角度，为基准定位面和待复位定位面的中心点到对应顶点的连线所形成的夹角。
在一实施例中，步骤S5中，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，进一步包括以下步骤：
通过基准定位面和待复位定位面的中心点到对应顶点的连线的向量，用叉乘法计算共同的垂向量；
过待复位定位面的中心点做平行于垂向量的平行线；
将待复位骨折体以此平行线为旋转轴，调节第三复位指标的大小，直至第三复位指标满足预设阈值。
在一实施例中，所述步骤S5的复位评价参数，包括第三综合复位评价指标MK3，公式如下：
当ID3≥TolA3时：
Figure BDA0002703958680000043
当TolA3＞ID3≥TolB3时：
Figure BDA0002703958680000044
当TolB3＞ID3时：
Figure BDA0002703958680000045
其中，ID3为第三复位指标，TolA3为第三复位指标的合格容差限值，TolB3为第三复位指标的良好容差限值。
本发明提供的交互式模拟骨折复位方法，使用尽量符合医生的物理直觉和医学传统的方式，提供直观的骨科空间方位(距离和方向)表示方法，通过简便的交互操作以及合理的复位评价参数标准，协助医生更快地模拟骨折体复位。
附图说明
本发明上述的以及其他的特征、性质和优势将通过下面结合附图和实施例的描述而变的更加明显，在附图中相同的附图标记始终表示相同的特征，其中：
图1揭示了根据本发明一实施例的交互式模拟骨折复位方法的流程图；
图2揭示了根据本发明一实施例的第一复位指标的示意图；
图3揭示了根据本发明一实施例的第二复位指标的示意图；
图4揭示了根据本发明一实施例的第三复位指标的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释发明，并不用于限定发明。
为了使医生可以快速便捷的进行术前准备，本发明提出了一种交互式计算机模拟骨折复位的方法，可以用于虚拟骨折复位，也可以用于植入物放置和导板配置等其他相关骨科治疗。
图1揭示了根据本发明一实施例的交互式模拟骨折复位方法的流程图，如图1所示，本发明提出的一种交互式模拟骨折复位方法，包括以下步骤：
S1、扫描并构建骨折体三维模型；
S2、选择骨折体三维模型中的一部分骨折体为基准骨折体，另一部分骨折体为待复位骨折体；
S3、在基准骨折体选取不少于三个基准点生成基准定位面；
S4、根据基准定位面的基准点，通过配准算法生成待复位骨折体的对应基准点，确定对应的待复位定位面；
S5、根据基准定位面与待复位定位面的相对位置关系，形成复位评价参数，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，实现骨折体三维模型的复位。
下面对每一步进行详细的说明。
S1、扫描并构建骨折体三维模型。
将CT断层扫描得到的图像进行滤波和区域分割，根据灰度阈值等特征参数提取出待修复骨折体的骨骼轮廓，并进行去噪和光滑处理，构建骨折体三维模型。
S2、选择骨折体三维模型中的一部分骨折体为基准骨折体，另一部分骨折体为待复位骨折体。
较佳的，选取体积较大的骨折体为基准骨折体，选择另一较小的骨折体作为待复位骨折体。
S3、在基准骨折体选取不少于三个基准点生成基准定位面。
在基准骨折体的骨折面轮廓线上，选取距离较远且不共线的三个基准点A，B，C。
将三个基准点形成三角形ABC，三角形ABC作为基准定位面，用于绑定基准骨折体，并用来表征基准骨折体。
S4、根据基准定位面的基准点，通过配准算法生成待复位骨折体的对应基准点，确定对应的待复位定位面。
使用配准算法将两个骨折断面进行配准，得到待复位骨折体的骨折断面上与A，B，C对应的三点a，b，c。
将A，B，C三个基准点的坐标，乘以配准位移矩阵，分别得到待复位骨折体的对应的三个基准点a，b，c。
a，b，c三个基准点形成三角形abc，三角形abc作为待复位定位面，用于绑定待复位骨折体，并且用来表征待复位骨折体。
本实施例中，配准算法为ICP(Iterative Closest Points，迭代最近点)算法。
ICP算法，是一种广泛使用的点云匹配方法，是求解两个点云集合转换关系的常用方法。ICP配准的目的是将具有相似几何关系的两个点集通过对点集的整体缩放旋转平移，移动到尽可能相近或相同的位置。
ICP算法本质上是基于最小二乘法的最优配准方法。该算法重复进行选择对应关系最近点对，计算最优刚体变换，直到满足正确配准的收敛精度要求。
在本实施例中，配准位移矩阵通过以下方法计算得到：
基准骨折体的骨折面的外轮廓线和待复位骨折体上骨折面的外轮廓线，进行ICP算法配准，得到配准位移矩阵。
假设基准骨折体的骨折面的外轮廓线的点云数据集合P和待复位骨折体上骨折面的外轮廓线的点云数据集合G，同时假设两组点云数据集合存在局部几何特征相似的部分，集合P可以通过叉乘一个配准位移矩阵，进行旋转和平移变换到集合G。
在本实施例中，配准位移矩阵大小为4×4。配准位移矩阵为M。
对于每次旋转平移变换后，计算集合P的所有采样点到集合G对应最近点距离，用最小二乘法求方差得出最小二乘误差值；
判断最小二乘误差值是否满足设定条件；
所述设定条件为最小二乘误差值小于设定阈值，或迭代次数达到上限，或每次重新迭代后最小二乘误差总在一个很小的范围内不再发生变化，如果最小二乘误差值满足设定条件，则迭代计算结束，否则继续进行迭代。
最后，得到的4×4矩阵为配准位移矩阵M。
S5、根据基准定位面与待复位定位面的相对位置关系，形成复位评价参数，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，实现骨折体三维模型的复位。
本实施例中，复位评价参数，根据基准定位面与待复位定位面的相对位置(距离和方向)关系形成，进一步包括第一复位指标ID1、第二复位指标ID2和第三复位指标ID3。三个指标没有主次之分，复位顺序也没有先后之分。
计算特定坐标系下的复位评价参数，使待复位骨折体在特定坐标系下移动；
在特定坐标系下，使三个复位指标逐步减少达到预设阈值，实现所需的骨折位置复位。
在本实施例中，第一复位指标ID1、第二复位指标ID2或第三复位指标ID3的预设阈值，进一步包括合格容差限制和良好容差限值。
更进一步的，所述步骤S5的骨折复位调节方法，进一步包括以下步骤：
定义特定的键盘或者按钮命令，使得用户的骨折复位操作分别单次独立沿着第一复位指标ID1、第二复位指标ID2或者第三复位指标ID3增大或者缩小的方向进行移动和旋转；
在移动和旋转时，实时显示第一复位指标ID1、第二复位指标ID2和第三复位指标ID3的数值；
调节位置使得第一复位指标ID1、第二复位指标ID2或第三复位指标ID3逐渐减小；
当第一复位指标ID1、第二复位指标ID2或第三复位指标ID3小于其对应的合格容差限制或者良好容差限值时，使用高亮数值或其他方式进行提醒；
骨折模型复位完成。
本实施例中，所述合格容差限值，为骨折复位状况为合格状态的容差限值。
所述良好容差限值，为骨折复位状况为良好状态的容差限值，其数值小于合格容差限值。
当复位指标小于对应的合格容差限值时，并且大于对应的良好容差限值时，骨折复位状况为合格状态。
当复位指标小于对应的良好容差限值时，骨折复位情况为良好状态。
下面通过图2-4，分别说明各个复位评价参数以及对应的调节方法。
图2-4分别揭示了根据本发明一实施例的第一复位指标、第二复位指标和第三复位指标的示意图，如图2-4所示，长圆柱代表基准骨折体，短圆柱为待复位骨折体。
A、B、C点为基准骨折体的骨折面上选取得三点，并形成了面ABC代表基准骨折面，通过A，B，C三点坐标信息，计算得到基准定位面三角形ABC的的第一中心点X的坐标。
更进一步的，在本实施例中，三角形ABC的中心点X的坐标公式为：
uX＝(uA+uB+uC)/3；
vX＝(vA+vB+vC)/3；
wX＝(wA+wB+wC)/3；
其中，(uX，vX，wX)为中心点X的坐标，(uA，vA，wA)为A点坐标，(uB，vB，wB)为B点坐标，(uC，vC，wC)为C点坐标。
a、b、c点为待复位骨折体的骨折面上，与面ABC配准得到的对应三点，面abc代表待复位骨折面，通过a，b，c三点坐标信息，计算得到待复位定位面三角形abc的第二中心点x的坐标。
更进一步的，在本实施例中，三角形abc的中心点x的坐标公式为：
ux＝(ua+ub+uc)/3；
vx＝(va+vb+vc)/3；
wx＝(wa+wb+wc)/3；
其中，(ux，vx，wx)为中心点X的坐标，(ua，va，wa)为a点坐标，(ub，vb，wb)为b点坐标，(uc，vc，wc)为c点坐标。
第一复位指标ID1为直线待复位距离，其定义方法如图2所示。
图2揭示了根据本发明一实施例的第一复位指标的示意图，如图2所示的第一复位指标ID1，为连接第一中心点X和第二中心点x得到的线段
Figure BDA0002703958680000091
。
Figure BDA0002703958680000092
的长度为直线待复位距离，定义为第一复位指标ID1。
根据第一复位指标ID1进行调节待复位骨折体的相对位置，进一步包括以下步骤：
通过计算基准定位面和待复位定位面的中心点的向量
Figure BDA0002703958680000093
的长度，使待复位骨折体沿此向量进行移动，调节第一复位指标ID1的大小，直至第一复位指标满足合格容差限值或良好容差限值。
第二复位指标ID2为轴向待复位角度，其定义方法如图3所示。
图3揭示了根据本发明一实施例的第二复位指标的示意图，如图3所示的第二复位指标ID2，为基准定位面的法线F与待复位定位面法线f形成的夹角为轴向待复位角度。
根据基准定位面ABC的任意两边的向量利用向量叉乘法求出基准定位面的法向量，过基准定位面的中心点X做与法向量方向平行的直线，作为基准定位面ABC的法线F。线p为法线F的平行线。
根据待复位定位面abc的任意两边的向量利用向量叉乘法求出待复位定位面的法向量，过待复位定位面的中心点x做与待复位定位面法向量方向平行的直线，作为待复位定位面abc的法线f。线q为法线f的平行线。
在图3中，线p和线q的夹角角度等于法线F和法线f的夹角角度，也就是第二复位指标ID2。
根据第二复位指标ID2进行调节待复位骨折体的相对位置，进一步包括以下步骤：
通过待复位定位面的法线f的向量
Figure BDA0002703958680000101
和基准定位面的法线F的向量/>
Figure BDA0002703958680000102
用叉乘法计算共同的垂向量/>
Figure BDA0002703958680000103
过待复位定位面的中心点x做平行于垂向量
Figure BDA0002703958680000104
的平行线；
待复位骨折体以此平行线为旋转轴，调节第二复位指标ID2的大小，直至第二复位指标满足合格容差限值或良好容差限值。
第三复位指标ID3为周向待复位角度，其定义方法如图4所示。
图4揭示了根据本发明一实施例的第三复位指标的示意图，如图4所示的第三复位指标ID3，为周向待复位角度，即连线
Figure BDA0002703958680000105
和连线/>
Figure BDA0002703958680000106
的夹角。
在基准定位面ABC上连接第一中心点X与对应顶点A，得到连线
Figure BDA0002703958680000107
在待复位定位面abc上连接第二中心点x与对应顶点a，得到连线
Figure BDA0002703958680000108
线k为线
Figure BDA0002703958680000109
的平行线，线j为线/>
Figure BDA00027039586800001010
的平行线。
在图4中，线k和线j的夹角角度等于线
Figure BDA00027039586800001011
与线/>
Figure BDA00027039586800001012
的夹角角度，也就是第三复位指标ID3。
根据第三复位指标ID3进行调节待复位骨折体的相对位置，进一步包括以下步骤：
通过基准定位面和待复位定位面的中心点到对应顶点的连线的向量
Figure BDA00027039586800001013
和/>
Figure BDA00027039586800001014
用叉乘法计算共同的垂向量/>
Figure BDA00027039586800001015
过待复位定位面的中心点x做平行于垂向量
Figure BDA00027039586800001016
的平行线；
待复位骨折体以此平行线为旋转轴，调节第三复位指标ID3的大小，直至第三复位指标满足合格容差限值或良好容差限值。
对三个复位指标分别设置对应的合格容差限值和良好容差限值：
第一复位指标ID1的合格容差限值TolA1，良好容差限值TolB1；
第二复位指标ID2的合格容差限值TolA2，良好容差限值TolB2；
第三复位指标ID3的合格容差限值TolA3，良好容差限值TolB3。
由于人体骨骼的大小尺寸差别极大，一般来说对于同一类型的骨骼，应该具有相近的容差限值，但是对于不同的骨骼，合格容差限值和良好容差限值可以有所调整。
可选的，合格容差限值和良好容差限值，由根据使用者经验或者统计数据标定来制定。
作为较佳实施例，对第一复位指标ID1、第二复位指标ID2和第三复位指标ID3：
对应的合格容差限值，预设值分别为：TolA1＝1mm，TolA2＝1度，TolA3＝1度；
对应的良好容差限值，预设值分别为：TolB1＝0.2mm，TolB2＝0.2度，TolB3＝0.2度。
上述三个复位指标作为复位评价参数，对不同位置的骨折的统一评价精度有差别。实际骨折复位时，绝大多数骨折无法达到解剖复位所要求的精度。按照手术方式分为切开复位和闭合牵引复位的复位方法和精度，需要根据患者伤病、体质和治疗情况决定。比如上肢骨折即使达不到解剖复位精度，也通常使用闭合复位。绝大部分情况下三个复位指标不可能同时为0。
更进一步的，本实施例中，复位评价参数进一步包括综合复位评价指标，综合复位评价指标根据复位指标计算得到，可以提醒使用者骨折复位的评价，从而可以对所有骨科治疗进行统一的复位精度评价。
综合复位评价指标是分别将三个复位指标进行综合计算得到的综合评价，其对应的预设阈值同样可以由根据使用者经验或者统计数据标定来制定。
第一复位指标ID1对应的第一综合复位评价指标MK1，公式如下：
当ID1≥TolA1时，
Figure BDA0002703958680000121
当TolA1＞ID1≥TolB1时，
Figure BDA0002703958680000122
当TolB1＞ID1时，
Figure BDA0002703958680000123
其中，ID1为第一复位指标，TolA1为第一复位指标的合格容差限值，TolB1为第一复位指标的良好容差限值。
第二复位指标ID2对应的第二综合复位评价指标MK2，公式如下：
当ID2≥TolA2时，
Figure BDA0002703958680000124
当TolA2＞ID2≥TolB2时，
Figure BDA0002703958680000125
当TolB2＞ID2时，
Figure BDA0002703958680000126
其中，ID2为第二复位指标，TolA2为第二复位指标的合格容差限值，TolB2为第二复位指标的良好容差限值。
第三复位指标ID3对应的第三综合复位评价指标MK3，公式如下：
当ID3≥TolA3时，
Figure BDA0002703958680000127
当TolA3＞ID3≥TolB3时，
Figure BDA0002703958680000128
当TolB3＞ID3时，
Figure BDA0002703958680000129
其中，ID3为第三复位指标，TolA3为第三复位指标的合格容差限值，TolB3为第三复位指标的良好容差限值。
本发明提供的交互式模拟骨折复位方法，使用尽量符合医生的物理直觉和医学传统的方式，提供直观的骨科空间方位(距离和方向)表示方法，通过简便的交互操作以及合理的复位评价参数标准，协助医生更快地模拟骨折体复位。
尽管为使解释简单化将上述方法图示并描述为一系列动作，但是应理解并领会，这些方法不受动作的次序所限，因为根据一个或多个实施例，一些动作可按不同次序发生和/或与来自本文中图示和描述或本文中未图示和描述但本领域技术人员可以理解的其他动作并发地发生。
如本申请和权利要求书中所示，除非上下文明确提示例外情形，“一”、“一个”、“一种”和/或“该”等词并非特指单数，也可包括复数。一般说来，术语“包括”与“包含”仅提示包括已明确标识的步骤和元素，而这些步骤和元素不构成一个排它性的罗列，方法或者设备也可能包含其他的步骤或元素。
上述实施例是提供给熟悉本领域内的人员来实现或使用本发明的，熟悉本领域的人员可在不脱离本发明的发明思想的情况下，对上述实施例做出种种修改或变化，因而本发明的保护范围并不被上述实施例所限，而应该是符合权利要求书提到的创新性特征的最大范围。

Claims (8)

1.一种交互式模拟骨折复位的方法，其特征在于，包括以下步骤：
S1、扫描并构建骨折体三维模型；
S2、选择骨折体三维模型中的一部分骨折体为基准骨折体，另一部分骨折体为待复位骨折体；
S3、在基准骨折体选取不少于三个基准点生成基准定位面；
S4、根据基准定位面的基准点，通过配准算法生成待复位骨折体的对应基准点，确定对应的待复位定位面；
S5、根据基准定位面与待复位定位面的相对位置关系，形成复位评价参数，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，实现骨折体三维模型的复位；
所述步骤S3，进一步包括以下步骤：
S31、在基准骨折体的骨折面轮廓线上，选取距离较远且不共线的三个基准点；
S32、将基准点连接形成基准定位面，用于绑定表征基准骨折体；
所述步骤S5的复位评价参数，包括第一复位指标，为直线待复位距离：
所述直线待复位距离，为基准定位面和待复位定位面的中心点的连线距离；
所述步骤S5的复位评价参数，包括第二复位指标，为轴向待复位角度：
所述轴向待复位角度，为基准定位面的法线与待复位定位面的法线形成的夹角；
所述步骤S5的复位评价参数，包括第三复位指标，为周向待复位角度：
所述周向待复位角度，为基准定位面和待复位定位面的中心点到对应顶点的连线所形成的夹角。
2.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，所述步骤S4的配准算法，为迭代最近点算法，进一步包括以下步骤：
将基准骨折体的骨折面的外轮廓点云集合P，与待复位骨折体的骨折面的外轮廓点云集合G进行最近点距离的最小二乘法迭代计算，得到配准位移矩阵M；
将基准骨折体的三个基准点，通过叉乘配准位移矩阵M，得到待复位骨折体的对应的三个待复位基准点。
3.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，所述步骤S5中，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，进一步包括以下步骤：
通过连接基准定位面和待复位定位面的中心点的向量，使待复位骨折体沿此向量进行移动，调节第一复位指标的大小，直至第一复位指标满足预设阈值。
4.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，所述步骤S5的复位评价参数，包括第一综合复位评价指标MK1，公式如下：
当ID1≥TolA1时，
Figure FDA0003946137590000021
当TolA1＞ID1≥TolB1时，
Figure FDA0003946137590000022
当TolB1＞ID1时，
Figure FDA0003946137590000023
其中，ID1为第一复位指标，TolA1为第一复位指标的合格容差限值，TolB1为第一复位指标的良好容差限值。
5.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，所述步骤S5中，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，进一步包括以下步骤：
通过待复位定位面的法线的向量和基准定位面的法线的向量，用叉乘法计算共同的垂向量；
过待复位定位面的中心点做平行于垂向量的平行线；
待复位骨折体以此平行线为旋转轴，调节第二复位指标的大小，直至第二复位指标满足预设阈值。
6.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，所述步骤S5的复位评价参数，包括第二综合复位评价指标MK2，公式如下：
当ID2≥TolA2时，
Figure FDA0003946137590000031
当TolA2＞ID2≥TolB2时，
Figure FDA0003946137590000032
当TolB2＞ID2时，
Figure FDA0003946137590000033
其中，ID2为第二复位指标，TolA2为第二复位指标的合格容差限值，TolB2为第二复位指标的良好容差限值。
7.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，步骤S5中，调节待复位骨折体的相对位置，直至复位评价参数满足预设阈值，进一步包括以下步骤：
通过基准定位面和待复位定位面的中心点到对应顶点的连线的向量，用叉乘法计算共同的垂向量；
过待复位定位面的中心点做平行于垂向量的平行线；
将待复位骨折体以此平行线为旋转轴，调节第三复位指标的大小，直至第三复位指标满足预设阈值。
8.根据权利要求1所述的交互式模拟骨折复位的方法，其特征在于，所述步骤S5的复位评价参数，包括第三综合复位评价指标MK3，公式如下：
当ID3≥TolA3时，
Figure FDA0003946137590000034
当TolA3＞ID3≥TolB3时，
Figure FDA0003946137590000041
当TolB3＞ID3时，
Figure FDA0003946137590000042
其中，ID3为第三复位指标，TolA3为第三复位指标的合格容差限值，TolB3为第三复位指标的良好容差限值。
CN202011031888.7A 2020-09-27 2020-09-27 一种交互式模拟骨折复位的方法 Active CN112164443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011031888.7A CN112164443B (zh) 2020-09-27 2020-09-27 一种交互式模拟骨折复位的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011031888.7A CN112164443B (zh) 2020-09-27 2020-09-27 一种交互式模拟骨折复位的方法

Publications (2)

Publication Number Publication Date
CN112164443A CN112164443A (zh) 2021-01-01
CN112164443B true CN112164443B (zh) 2023-03-31

Family

ID=73864207

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011031888.7A Active CN112164443B (zh) 2020-09-27 2020-09-27 一种交互式模拟骨折复位的方法

Country Status (1)

Country Link
CN (1) CN112164443B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113345112B (zh) * 2021-05-25 2023-06-13 上海大学 一种长骨骨折断面点云预处理及配准方法
CN113781495B (zh) * 2021-07-21 2023-03-24 天津大学 面向并联骨折手术机器人的复位轨迹自动式规划方法
CN114974509B (zh) * 2022-05-26 2024-05-10 哈尔滨工业大学 骨折复位路径规划方法、骨折复位方法及电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107174342A (zh) * 2017-03-21 2017-09-19 哈尔滨工程大学 一种计算机辅助骨折整复程度度量方法
CN107967932A (zh) * 2017-12-21 2018-04-27 成都真实维度科技有限公司 基于三维c型臂机的股骨颈骨折闭合复位评价方法及系统
CN109119140A (zh) * 2018-08-27 2019-01-01 北京大学口腔医学院 一种用于准确治疗陈旧性颧骨骨折的计算机辅助导航方法
CN109528305A (zh) * 2018-10-12 2019-03-29 河海大学常州校区 计算机辅助的股骨粗隆间骨折对位对线半自动复位方法
CN109984843A (zh) * 2019-03-20 2019-07-09 中国科学院苏州生物医学工程技术研究所 骨折闭合复位导航系统及方法
US10687899B1 (en) * 2016-07-05 2020-06-23 Smith & Nephew, Inc. Bone model correction angle determination

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10687899B1 (en) * 2016-07-05 2020-06-23 Smith & Nephew, Inc. Bone model correction angle determination
CN107174342A (zh) * 2017-03-21 2017-09-19 哈尔滨工程大学 一种计算机辅助骨折整复程度度量方法
CN107967932A (zh) * 2017-12-21 2018-04-27 成都真实维度科技有限公司 基于三维c型臂机的股骨颈骨折闭合复位评价方法及系统
CN109119140A (zh) * 2018-08-27 2019-01-01 北京大学口腔医学院 一种用于准确治疗陈旧性颧骨骨折的计算机辅助导航方法
CN109528305A (zh) * 2018-10-12 2019-03-29 河海大学常州校区 计算机辅助的股骨粗隆间骨折对位对线半自动复位方法
CN109984843A (zh) * 2019-03-20 2019-07-09 中国科学院苏州生物医学工程技术研究所 骨折闭合复位导航系统及方法

Also Published As

Publication number Publication date
CN112164443A (zh) 2021-01-01

Similar Documents

Publication Publication Date Title
CN112164443B (zh) 一种交互式模拟骨折复位的方法
US10217217B2 (en) Systems and methods for obtaining 3-D images from X-ray information
JP4624173B2 (ja) 幾何学的オブジェクトをデジタルの医学的画像内に自動写像する方法
US20210012492A1 (en) Systems and methods for obtaining 3-d images from x-ray information for deformed elongate bones
US9173716B2 (en) Computer-aided planning with dual alpha angles in femoral acetabular impingement surgery
US20210007806A1 (en) A method for obtaining 3-d deformity correction for bones
CN110214341A (zh) 重建颅骨的方法
US9020223B2 (en) Method for determining bone resection on a deformed bone surface from few parameters
WO2012027150A2 (en) Personalized orthopedic implant cad model generation
TWI711428B (zh) 用於醫療用具的光學追蹤系統及訓練系統
US20220071708A1 (en) Patient positioning using a skeleton model
CN109528305A (zh) 计算机辅助的股骨粗隆间骨折对位对线半自动复位方法
WO2019180746A1 (en) A method for obtaining 3-d deformity correction for bones
CN114332378A (zh) 基于二维医学影像的人体骨骼三维模型获取方法及系统
CN104933672B (zh) 基于快速凸优化算法配准三维ct与超声肝脏图像的方法
Wang et al. Deformable torso phantoms of Chinese adults for personalized anatomy modelling
US7792360B2 (en) Method, a computer program, and apparatus, an image analysis system and an imaging system for an object mapping in a multi-dimensional dataset
King et al. Tracking liver motion using 3-D ultrasound and a surface based statistical shape model
WO2019180747A1 (en) Systems and methods for obtaining patient specific instrument designs
Dobson et al. Feature based models for anatomical data fitting
CN111968236A (zh) 一种基于K-D tree方法的对位对线骨折复位系统
Gottschling et al. A system for performing automated measurements on large bone databases
Vitković et al. Geometrical models of human bones and implants, and their usage in application for preoperative planning in orthopedics
Sinh et al. Integration of computer-assisted fracture reduction system and a hybrid 3-DOF-RPS mechanism for assisting the orthopedic surgery
Bazsó et al. Personalisable vertebral body model development

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant