CN107550567A - 一种骨折复位的计算机辅助设计方法 - Google Patents

一种骨折复位的计算机辅助设计方法 Download PDF

Info

Publication number
CN107550567A
CN107550567A CN201710700551.2A CN201710700551A CN107550567A CN 107550567 A CN107550567 A CN 107550567A CN 201710700551 A CN201710700551 A CN 201710700551A CN 107550567 A CN107550567 A CN 107550567A
Authority
CN
China
Prior art keywords
fracture
fracture site
distal end
pieces
disconnected
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.)
Pending
Application number
CN201710700551.2A
Other languages
English (en)
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.)
Peking University
Beijing Friendship Hospital
Original Assignee
Peking University
Beijing Friendship Hospital
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 Peking University, Beijing Friendship Hospital filed Critical Peking University
Priority to CN201710700551.2A priority Critical patent/CN107550567A/zh
Publication of CN107550567A publication Critical patent/CN107550567A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Surgical Instruments (AREA)

Abstract

本发明公开了一种骨折复位的计算机辅助设计方法。本发明确定骨折断端的断截面长短主轴、建立三维直角坐标系、确定中心主轴及其方向上的矢量在三维坐标系中的四元数表示、选定易操作的骨折断端作为远心端骨折断端、分别进行转动、自旋和平移的操作得到骨折断端共轴最优路径点,实现断截面吻合,完成骨折复位;本发明应用四元数表示两块骨折断端中心主轴方向的矢量,进而得到每次转动的转轴与转动角度,四元数在描述旋转变换时求解速度更快,并可有效避免万向锁等问题的出现;本发明通过对骨折断端进行旋转、自旋和平移操作,得到最优的骨折断端共轴轨迹点,可为医生诊断与操作提供有效信息。

Description

一种骨折复位的计算机辅助设计方法
技术领域
本发明涉及计算机手术导航领域,具体涉及一种骨折复位的计算机辅助设计方法。
背景技术
骨折是指骨结构的连续性完全或部分断裂,是和人们生活有密切关系的常见创伤。骨折患者的正常解剖结构被破坏,影响了肢体功能,降低生活质量,给社会和家庭带来了负担。随着现代影像技术的不断发展,X线、计算机体层摄影(CT)、核磁共振成像能将绝大多数骨折明确诊断,清晰的显示骨折类型和移位情况等,指导临床治疗。CT扫描获取的医学图像能清楚反映复杂关节内和脊柱等骨折的详细信息,重建的三维图像能够更便捷地指导外科手术治疗方案,在临床上得到了广泛应用。
多数移位的骨折患者需要行骨折复位手术治疗。手术过程中,骨科医师依赖X线透视明确骨折断端位置、复位情况和内固定物的置入位置等,而过程中的射线暴露对患者和医师造成潜在的风险。目前常用的术中透视设备多为C形臂和G形臂,获取的图像为二维图像,骨科医师需结合临床经验判断骨折复位情况,手术操作难度大,因此需发展新的手术导航系统解决上述问题。
发明内容
为实现在减少射线透视次数的同时保证骨折断端信息的充分性与准确性,本发明提出了一种骨折复位的计算机辅助设计方法,利用计算机模拟给出骨折复位的最优路径,进而辅助医生操作,为医生操作带来便利的同时也能在一定程度上促进数字时代的医疗卫生发展。
本发明的一种骨折复位的计算机辅助设计方法,包括以下步骤:
1)通过CT扫描图像重建出两块骨折断端的三维的刚体点系模型,并确定两块骨折断端的髓腔内壁和骨外壁的空间刚体点系;
2)利用主成分分析(Principal Component Analysis,PCA)方法分别确定两块骨折断端的断截面的长短主轴,以此为基础给出各自对应的三维直角坐标系,并确定骨折断端的中心主轴;
3)判断两块骨折断端的可移动度,远心端骨折断端为可移动度大的骨折断端,将其作为可移动的骨折断端,近心端骨折断端为可移动度小的骨折断端,以近心端骨折断端为参照,向近心端骨折断端转动远心端骨折断端,实现两块骨折断端的共轴;
4)两块骨折断端共轴后,确定两断截面的断截面方程,确定远心端骨折断端的旋转角度,达到两个断截面的长主轴和短主轴均重合;
5)沿中心主轴无旋转的平移远心端骨折断端,实现两块骨折断端的完全复位,完成骨折复位的最优路径设计。
其中,在步骤2)中,利用主成分分析PCA方法分别确定两块骨折断端的断截面的长主轴和短主轴,并以此为基础给出各自对应的三维直角坐标系中的x1O1y1平面和x2O2y2平面;长主轴与短主轴的交点为断截面的中心点,也是三维直角坐标系中的原点O1和O2,将过O1点且垂直于x1O1y1平面的轴作为三维直角坐标系中的z1轴,同理将过O2点且垂直于x2O2y2平面的轴作为三维直角坐标系中的z2轴,并设定由断截面指向骨折断端外部的方向为正方向;确定平行于骨外壁且经过断截面中心点的直线,将其作为骨折断端的中心主轴,骨折断端的中心主轴与其各自的z轴之间的夹角为θ,θ≥0。
在步骤3)中,通过与心脏间的距离大小判断两块骨折断端的可移动度。以近心端骨折断端为参照,向近心端骨折断端移动远心端骨折断端,具体包括以下步骤:分别取在两块骨折断端各自的中心主轴方向上的矢量,并给出矢量在三维直角坐标系中的四元数表示,其中对于远心端骨折断端,转动前其中心主轴方向的矢量表示为r1=(m1,m2,m3),近心端骨折断端的中心主轴方向的矢量表示为恒定的r2=(n1,n2,n3);将r1=(m1,m2,m3)与r2=(n1,n2,n3)作为转动过程中的初始姿态与终止姿态,通过四元数球面线性插值的方法得到转动过程中姿态变换的空间位置表示,实现两块骨折断端的中心主轴的共轴。
在步骤4)中,确定两断截面的断截面方程,确定远心端骨折断端的旋转角度,达到两个断截面的长主轴和短主轴均重合,具体包括以下步骤:远心端骨折断端和近心端骨折断端的三维直角坐标系分别表示为C1(x1,y1,z1)和C2(x2,y2,z2),进一步由断截面的中心点与共同的中心主轴分别确定两断截面的断截面方程,其中远心端骨折断端的断截面方程表示为π1:a1x1+b1y1+c1z1=0,近心端骨折断端的断截面方程表示为π2:a2x2+b2y2+c2z2=0,且确定该截面对应的法向量为h,以向量h作为远心端骨折断端的截面方程的法向量并结合满足方程π2的一个点,给出两断截面法向量平行时的远心端骨折断端的断截面方程π3:a3x1+b3y1+c3z1=0,通过与方程π1的对比确定远心端骨折断端的旋转角度,达到两个断截面的长短主轴均重合。
在步骤5)中,沿中心主轴无旋转的平移远心端骨折断端,具体包括以下步骤:确定两块骨折断端共同的中心主轴与其各自断截面的交点分别为p1和p2,以远心端骨折断端的截面中的p1点为起点得到两点间的矢量L并沿中心主轴无旋转的平移远心端骨折断端,实现两块骨折断端的完全复位,完成骨折复位的最优路径设计。
本发明的优点:
本发明确定骨折断端的断截面长短主轴、建立三维直角坐标系、确定中心主轴及其方向上的矢量在三维坐标系中的四元数表示、选定易操作的骨折断端作为远心端骨折断端、分别进行转动、自旋和平移的操作得到骨折断端共轴最优路径点,完成骨折复位;本发明应用四元数表示两块骨折断端的中心主轴方向的矢量,进而得到每次转动的转轴与转动角度,四元数相对于欧拉角、旋转矩阵等方法,在描述旋转变换时求解速度更快,并可有效避免万向锁等问题的出现;本发明通过对骨折断端进行旋转、自旋和平移操作,得到最优的骨折断端共轴轨迹点,可为医生诊断与操作提供有效信息。
附图说明
图1为本发明的骨折复位的计算机辅助设计方法的流程图;
图2为本发明的骨折复位的计算机辅助设计方法的断截面建立的三维直角坐标系与中心轴线的示意图。
具体实施方式
下面结合附图,通过具体实施例,进一步阐述本发明。
如图1所示,本实施例的骨折复位的计算机辅助设计方法,包括以下步骤:
1)通过CT扫描图像重建出两块骨折断端的三维的刚体点系模型,并确定两块骨折断端的髓腔内壁和骨外壁的空间刚体点系;
2)利用主成分分析(Principal Component Analysis,PCA)方法分别确定两块骨折断端的断截面的长短主轴,并以此为基础给出各自对应的三维直角坐标系中的x1O1y1平面和x2O2y2平面;长主轴与短主轴的交点为断截面的中心点,也是三维直角坐标系中的原点O1和O2,将过O1点且垂直于x1O1y1平面的轴作为三维直角坐标系中的z1轴,同理将过O2点且垂直于x2O2y2平面的轴作为三维直角坐标系中的z2轴,并设定由断截面指向骨折断端外部的方向为正方向;确定平行于骨外壁且经过断截面中心点的直线,将其作为骨折断端的中心主轴,骨折断端的中心主轴与其各自的z轴之间的夹角为θ,θ≥0;
3)通过与心脏间的距离大小判断两块骨折断端的可移动度,可移动度大的骨折断端为远心端骨折断端,作为可移动的骨折断端,可移动度小的骨折断端为近心端骨折断端,以近心端骨折断端为参照,向近心端骨折断端移动远心端骨折断端;分别取在两块骨折断端各自的中心主轴方向上的矢量,并给出矢量在三维直角坐标系中的四元数表示,其中对于远心端骨折断端,转动前其中心主轴方向的矢量表示为r1=(m1,m2,m3),近心端骨折断端的中心主轴方向的矢量表示为恒定的r2=(n1,n2,n3);将r1=(m1,m2,m3)与r2=(n1,n2,n3)作为转动过程中的初始姿态与终止姿态,通过四元数球面线性插值的方法得到转动过程中姿态变换的空间位置表示,实现两块骨折断端的中心主轴的共轴;
4)两块骨折断端共轴后,其三维直角坐标系分别表示为C1(x1,y1,z1)和C2(x2,y2,z2),进一步由断截面的中心点与共同的中心主轴分别确定两断截面的断截面方程,其中远心端骨折断端的断截面方程表示为π1:a1x1+b1y1+c1z1=0,近心端骨折断端的断截面方程表示为π2:a2x2+b2y2+c2z2=0,且确定该截面对应的法向量为h,以向量h作为远心端骨折断端的截面方程的法向量并结合满足方程π2的一个点,给出两断截面法向量平行时的远心端骨折断端的断截面方程π3:a3x1+b3y1+c3z1=0,通过与方程π1的对比确定远心端骨折断端的旋转角度,达到两个断截面的长短主轴均重合;
5)确定两块骨折断端共同的中心主轴与其各自断截面的交点分别为p1和p2,以远心端骨折断端的截面中的p1点为起点得到两点间的矢量L并沿中心主轴无旋转的平移远心端骨折断端,实现两块骨折断端的完全复位,完成骨折复位的最优路径设计。
最后需要注意的是,公布实施例的目的在于帮助进一步理解本发明,但是本领域的技术人员可以理解:在不脱离本发明及所附的权利要求的精神和范围内,各种替换和修改都是可能的。因此,本发明不应局限于实施例所公开的内容,本发明要求保护的范围以权利要求书界定的范围为准。

Claims (6)

1.一种骨折复位的计算机辅助设计方法,其特征在于,所述计算机辅助设计方法包括以下步骤:
1)通过CT扫描图像重建出两块骨折断端的三维的刚体点系模型,并确定两块骨折断端的髓腔内壁和骨外壁的空间刚体点系;
2)利用主成分分析PCA方法分别确定两块骨折断端的断截面的长短主轴,以此为基础给出各自对应的三维直角坐标系,并确定骨折断端的中心主轴;
3)判断两块骨折断端的可移动度,远心端骨折断端为可移动度大的骨折断端,将其作为可移动的骨折断端,近心端骨折断端为可移动度小的骨折断端,以近心端骨折断端为参照,向近心端骨折断端转动远心端骨折断端,实现两块骨折断端的共轴;
4)两块骨折断端共轴后,确定两断截面的断截面方程,确定远心端骨折断端的旋转角度,达到两个断截面的长主轴和短主轴均重合;
5)沿中心主轴无旋转的平移远心端骨折断端,实现两块骨折断端的完全复位,完成骨折复位的最优路径设计。
2.如权利要求1所述的计算机辅助设计方法,其特征在于,在步骤2)中,利用主成分分析PCA方法分别确定两块骨折断端的断截面的长主轴和短主轴,并以此为基础给出各自对应的三维直角坐标系中的x1O1y1平面和x2O2y2平面;长主轴与短主轴的交点为断截面的中心点,也是三维直角坐标系中的原点O1和O2,将过O1点且垂直于x1O1y1平面的轴作为三维直角坐标系中的z1轴,同理将过O2点且垂直于x2O2y2平面的轴作为三维直角坐标系中的z2轴,并设定由断截面指向骨折断端外部的方向为正方向;确定平行于骨外壁且经过断截面中心点的直线,将其作为骨折断端的中心主轴,骨折断端的中心主轴与其各自的z轴之间的夹角为θ,θ≥0。
3.如权利要求1所述的计算机辅助设计方法,其特征在于,在步骤3)中,通过与心脏间的距离大小判断两块骨折断端的可移动度。
4.如权利要求1所述的计算机辅助设计方法,其特征在于,在步骤3)中,以近心端骨折断端为参照,向近心端骨折断端移动远心端骨折断端,具体包括以下步骤:分别取在两块骨折断端各自的中心主轴方向上的矢量,并给出矢量在三维直角坐标系中的四元数表示,其中对于远心端骨折断端,转动前其中心主轴方向的矢量表示为r1=(m1,m2,m3),近心端骨折断端的中心主轴方向的矢量表示为恒定的r2=(n1,n2,n3);将r1=(m1,m2,m3)与r2=(n1,n2,n3)作为转动过程中的初始姿态与终止姿态,通过四元数球面线性插值的方法得到转动过程中姿态变换的空间位置表示,实现两块骨折断端的中心主轴的共轴。
5.如权利要求1所述的计算机辅助设计方法,其特征在于,在步骤4)中,确定两断截面的断截面方程,确定远心端骨折断端的旋转角度,达到两个断截面的长主轴和短主轴均重合,具体包括以下步骤:远心端骨折断端和近心端骨折断端的三维直角坐标系分别表示为C1(x1,y1,z1)和C2(x2,y2,z2),进一步由断截面的中心点与共同的中心主轴分别确定两断截面的断截面方程,其中远心端骨折断端的断截面方程表示为π1:a1x1+b1y1+c1z1=0,近心端骨折断端的断截面方程表示为π2:a2x2+b2y2+c2z2=0,且确定该截面对应的法向量为h,以向量h作为远心端骨折断端的截面方程的法向量并结合满足方程π2的一个点,给出两断截面法向量平行时的远心端骨折断端的断截面方程π3:a3x1+b3y1+c3z1=0,通过与方程π1的对比确定远心端骨折断端的旋转角度,达到两个断截面的长短主轴均重合。
6.如权利要求1所述的计算机辅助设计方法,其特征在于,在步骤5)中,沿中心主轴无旋转的平移远心端骨折断端,具体包括以下步骤:确定两块骨折断端共同的中心主轴与其各自断截面的交点分别为p1和p2,以远心端骨折断端的截面中的p1点为起点得到两点间的矢量L并沿中心主轴无旋转的平移远心端骨折断端,实现两块骨折断端的完全复位,完成骨折复位的最优路径设计。
CN201710700551.2A 2017-08-16 2017-08-16 一种骨折复位的计算机辅助设计方法 Pending CN107550567A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710700551.2A CN107550567A (zh) 2017-08-16 2017-08-16 一种骨折复位的计算机辅助设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710700551.2A CN107550567A (zh) 2017-08-16 2017-08-16 一种骨折复位的计算机辅助设计方法

Publications (1)

Publication Number Publication Date
CN107550567A true CN107550567A (zh) 2018-01-09

Family

ID=60974576

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710700551.2A Pending CN107550567A (zh) 2017-08-16 2017-08-16 一种骨折复位的计算机辅助设计方法

Country Status (1)

Country Link
CN (1) CN107550567A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111329571A (zh) * 2020-03-10 2020-06-26 河北医科大学第三医院 一种下肢骨折复位方法
CN111590584A (zh) * 2020-05-27 2020-08-28 京东方科技集团股份有限公司 安全限位区的确定方法和装置、复位方法和医用机器人
CN112168347A (zh) * 2020-11-10 2021-01-05 哈尔滨理工大学 一种骨折复位的计算机辅助设计方法
CN112386327A (zh) * 2020-11-13 2021-02-23 毕建平 负压环境骨水泥单孔注射骨折固定系统及工作方法
CN112998854A (zh) * 2021-03-17 2021-06-22 杭州素问九州医疗科技有限公司 胫骨假体与胫骨的拟合方法及机器人手术系统
WO2021142213A1 (en) * 2020-01-09 2021-07-15 Smith & Nephew, Inc. Methods and arrangements to describe deformity of a bone
WO2021164614A1 (zh) * 2020-02-19 2021-08-26 京东方科技集团股份有限公司 骨折复位的计算机辅助方法、装置及电子设备
CN115869072A (zh) * 2023-02-27 2023-03-31 武汉市第四医院 一种骨折复位机器人系统及其控制方法
CN117132747A (zh) * 2023-10-25 2023-11-28 北京爱康宜诚医疗器材有限公司 基于骨骼模型的骨骼复位方法和装置

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021142213A1 (en) * 2020-01-09 2021-07-15 Smith & Nephew, Inc. Methods and arrangements to describe deformity of a bone
US11986249B2 (en) 2020-02-19 2024-05-21 Boe Technology Group Co., Ltd. Computer-aided method for fracture reduction, and electronic device
WO2021164614A1 (zh) * 2020-02-19 2021-08-26 京东方科技集团股份有限公司 骨折复位的计算机辅助方法、装置及电子设备
CN111329571A (zh) * 2020-03-10 2020-06-26 河北医科大学第三医院 一种下肢骨折复位方法
CN111329571B (zh) * 2020-03-10 2021-08-31 河北医科大学第三医院 一种下肢骨折复位系统
CN111590584B (zh) * 2020-05-27 2021-12-10 京东方科技集团股份有限公司 安全限位区的确定方法和装置、复位方法和医用机器人
CN111590584A (zh) * 2020-05-27 2020-08-28 京东方科技集团股份有限公司 安全限位区的确定方法和装置、复位方法和医用机器人
CN112168347A (zh) * 2020-11-10 2021-01-05 哈尔滨理工大学 一种骨折复位的计算机辅助设计方法
CN112386327A (zh) * 2020-11-13 2021-02-23 毕建平 负压环境骨水泥单孔注射骨折固定系统及工作方法
CN112386327B (zh) * 2020-11-13 2022-06-03 毕建平 负压环境骨水泥单孔注射骨折固定系统及工作方法
CN112998854A (zh) * 2021-03-17 2021-06-22 杭州素问九州医疗科技有限公司 胫骨假体与胫骨的拟合方法及机器人手术系统
CN115869072A (zh) * 2023-02-27 2023-03-31 武汉市第四医院 一种骨折复位机器人系统及其控制方法
CN117132747A (zh) * 2023-10-25 2023-11-28 北京爱康宜诚医疗器材有限公司 基于骨骼模型的骨骼复位方法和装置
CN117132747B (zh) * 2023-10-25 2024-03-19 北京爱康宜诚医疗器材有限公司 基于骨骼模型的骨骼复位方法和装置

Similar Documents

Publication Publication Date Title
CN107550567A (zh) 一种骨折复位的计算机辅助设计方法
EP4265214A1 (en) Navigation system and method for joint replacement surgery
KR102014355B1 (ko) 수술도구의 위치 정보 산출 방법 및 장치
CN109124764B (zh) 手术引导装置以及手术系统
US7853058B2 (en) Determining a viewpoint for navigating a virtual camera through a biological object with a lumen
JP2021112588A (ja) 副鼻腔手術用のガイドワイヤ操縦法
CN110946654A (zh) 一种基于多模影像融合的骨科手术导航系统
JP2019528123A (ja) 整形外科の固定制御及び操作
US20150234957A1 (en) Fenestration template for endovascular repair of aortic aneurysms
US6793496B2 (en) Mathematical model and a method and apparatus for utilizing the model
Murphy et al. Design and kinematic characterization of a surgical manipulator with a focus on treating osteolysis
US20030182091A1 (en) Modeling a flexible tube
US20020168618A1 (en) Simulation system for image-guided medical procedures
US20160228075A1 (en) Image processing device, method and recording medium
JPH08131403A (ja) 医用画像処理装置
JP2009545358A (ja) 回転x線走査計画システム
US20230210599A1 (en) Systems and methods for pre-operative visualization of a joint
CN112150543A (zh) 医疗成像设备的成像定位方法、装置、设备及存储介质
BR112013010677A2 (pt) aparelho de geração de imagens para a geração de imagem de um objeto, método de geração de imagens para a geração de imagem de um objeto e programa de computador de geraçã ode imagens para a geração de imagens de um objeto.
CN109350059B (zh) 用于肘部自动对准的组合的转向引擎和界标引擎
Chui et al. Training and pretreatment planning of interventional neuroradiology procedures–initial clinical validation
US20210298981A1 (en) Surgical bed, endoscopic surgical device, endoscopic surgical method, and system
CN115105203A (zh) 一种虚拟颅脑3d建模系统及三轴坐标定位方法
CN114283179A (zh) 基于超声图像的骨折远近端空间位姿实时获取与配准系统
CN111743629A (zh) 医疗导管及其定位方法和三维磁定位系统

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20180109