[go: up one dir, main page]

CN110613469B - 一种腿骨下肢力线自动检测方法及装置 - Google Patents

一种腿骨下肢力线自动检测方法及装置 Download PDF

Info

Publication number
CN110613469B
CN110613469B CN201910884469.9A CN201910884469A CN110613469B CN 110613469 B CN110613469 B CN 110613469B CN 201910884469 A CN201910884469 A CN 201910884469A CN 110613469 B CN110613469 B CN 110613469B
Authority
CN
China
Prior art keywords
point
center point
center
cross
coordinates
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
CN201910884469.9A
Other languages
English (en)
Other versions
CN110613469A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201910884469.9A priority Critical patent/CN110613469B/zh
Publication of CN110613469A publication Critical patent/CN110613469A/zh
Application granted granted Critical
Publication of CN110613469B publication Critical patent/CN110613469B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/505Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

一种腿骨下肢力线自动检测方法,包括:(1)获取CT影像数据;(2)对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;(3)针对股骨三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;(4)针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;(5)通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;(6)针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线。还有装置。

Description

一种腿骨下肢力线自动检测方法及装置
技术领域
本发明涉及医疗图像处理的技术领域,尤其涉及一种腿骨下肢力线自动检测方法,以及腿骨下肢力线自动检测装置。
背景技术
人工髋膝关节置换术是目前治疗骨性关节炎、重度类风湿性关节炎和软骨破坏等末期骨性疾病的最有效手段。随着我国已逐步进入老龄化社会,人工髋、膝关节置换量增长迅速,关节置换手术市场份额迅速增长。
髋膝关节置换手术的成功离不开力线的有效检测,精准找到下肢力线是保障整个手术成功的基础。下肢力线包括股骨解剖轴线、股骨机械轴线、胫骨解剖轴线、胫骨机械轴线。
其中,股骨机械轴是指股骨头球心与髁间凹中心点连线;股骨全长中点及股骨远端关节线上10cm处取两个横截面,求取这一不规则髓腔形体的中心,这两中心点的连线为股骨解剖轴。
其中,胫骨机械轴是指胫骨棘中点稍外侧及距骨中点连线为胫骨机械轴。胫骨解剖轴:胫骨全长中点及胫骨远端关节线上10cm处取两个横截面,求取这一不规则髓腔形体的中心,这两中心点的连线为胫骨解剖轴。
在临床中,常规髋膝关节置换术通过术前X光片进行解剖轴线和机械轴线的确定。医生确定解剖轴线和机械轴线的常规做法是,在病人的X光片上进行选点连线,这种确定解剖轴线和机械轴线的方式严重依赖手术者经验,可靠性差,成为制约手术效果的重要因素。
因此,解决该问题的有效方式应该是根据病人CT或X光图像自动提取下肢力线,在保证精度的前提下,减少由于医生手动取点画线造成的人为误差。
发明内容
为克服现有技术的缺陷,本发明要解决的技术问题是提供了一种腿骨下肢力线自动检测方法,其能够保证机械轴线和解剖轴线的精确自动检测,减少了由人为因素引起的主观误差,在减少医生工作量的同时,提高手术的可靠性。
本发明的技术方案是:这种腿骨下肢力线自动检测方法,其包括以下步骤:
(1)获取CT影像数据:利用CT设备对病人的预定部位进行术前扫描,以获取术前CT影像;
(2)对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;
(3)针对股骨三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;
(4)针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;
(5)通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;
(6)针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线。
本发明通过双剖切面获取质心的方法来确定股骨头球心点位置,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标,从而获得股骨机械轴线和胫骨机械轴线,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线,因此能够保证机械轴线和解剖轴线的精确自动检测,减少了由人为因素引起的主观误差,在减少医生工作量的同时,提高手术的可靠性。
还提供了一种腿骨下肢力线自动检测装置,该装置包括:
CT影像数据获取模块,其配置来对病人的预定部位进行术前扫描,以获取术前CT影像;
骨骼三维模型获取模块,其配置来对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;
股骨头球心点获取模块,其配置来针对股骨三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;
配准模块,其配置来针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;
机械轴线获取模块,其配置来通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;
解剖轴线获取模块,其配置来针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线。
附图说明
图1示出了膝关节中心与胫骨平台中心。
图2示出了踝关节中心。
图3是根据本发明的腿骨下肢力线自动检测方法的流程图。
图4是根据本发明的腿骨下肢力线自动检测方法的一个具体实施例的流程图。
具体实施方式
如图3所示,这种腿骨下肢力线自动检测方法,其包括以下步骤:
(1)获取CT影像数据:利用CT设备对病人的预定部位进行术前扫描,以获取术前CT影像;
(2)对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;
(3)针对股骨三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;
(4)针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;
(5)通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;
(6)针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线。
本发明通过双剖切面获取质心的方法来确定股骨头球心点位置,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标,从而获得股骨机械轴线和胫骨机械轴线,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线,因此能够保证机械轴线和解剖轴线的精确自动检测,减少了由人为因素引起的主观误差,在减少医生工作量的同时,提高手术的可靠性。
优选地,所述步骤(3)中,在股骨头球三维模型表面任选一点作为初始点位置,第一步过此点沿着横断面剖切,得到横断剖切面,提取剖切面闭合边缘,然后利用迭代循环搜索方法计算剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为横断剖切面的质心点坐标;第二步过横断剖切面质心点沿着矢状面剖切,得到矢状剖切面,提取剖切面闭合边缘,然后利用迭代循环搜索方法计算剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为矢状剖切面的质心点坐标;第三步,矢状剖切面质心点坐标作为股骨头球心点坐标。
优选地,所述步骤(4)中,选择一例腿骨CT影像作为标准模板,在标准模板上标注3个中心点的位置,3个中心点为膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点,将标准模板CT配准到待测CT上,获得配准矩阵,标准模板上的3个点乘以配准矩阵,来获取待测CT上的3个点位置坐标,重建出待测CT的腿骨三维模型,以便直观观察到3个中心点(即生理特征点)在腿骨模型的位置。
优选地,所述步骤(6)中,提取股骨三维模型的boundingbox的高度值h,选取距离boundingbox最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将此3个中心点利用最小二乘法拟合为一条直线,此直线为股骨解剖轴线。
优选地,所述步骤(6)中,提取胫骨三维模型的boundingbox的高度值h,选取距离boundingbox最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将此3个中心点利用最小二乘法拟合为一条直线,此直线为胫骨解剖轴线。
本领域普通技术人员可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,包括上述实施例方法的各步骤,而所述的存储介质可以是:ROM/RAM、磁碟、光盘、存储卡等。因此,与本发明的方法相对应的,本发明还同时包括一种腿骨下肢力线自动检测装置,该装置通常以与方法各步骤相对应的功能模块的形式表示。该装置包括:
CT影像数据获取模块,其配置来对病人的预定部位进行术前扫描,以获取术前CT影像;
骨骼三维模型获取模块,其配置来对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;
股骨头球心点获取模块,其配置来针对股骨三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;
配准模块,其配置来针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;
机械轴线获取模块,其配置来通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;
解剖轴线获取模块,其配置来针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线。
更进一步地,所述股骨头球心点获取模块执行以下步骤:在股骨头球三维模型表面任选一点作为初始点位置,第一步过此点沿着横断面剖切,得到横断剖切面,提取剖切面闭合边缘,然后利用迭代循环搜索方法计算剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为横断剖切面的质心点坐标;第二步过横断剖切面质心点沿着矢状面剖切,得到矢状剖切面,提取剖切面闭合边缘,然后利用迭代循环搜索方法计算剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为矢状剖切面的质心点坐标;第三步,矢状剖切面质心点坐标作为股骨头球心点坐标。
更进一步地,所述配准模块执行以下步骤:选择一例腿骨CT影像作为标准模板,在标准模板上标注3个中心点的位置,3个中心点为膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点,将标准模板CT配准到待测CT上,获得配准矩阵,标准模板上的3个点乘以配准矩阵,来获取待测CT上的3个点位置坐标,重建出待测CT的腿骨三维模型,以便直观观察到3个中心点(即,生理特征点)在腿骨模型的位置。
更进一步地,所述解剖轴线获取模块执行以下步骤:提取股骨三维模型的boundingbox的高度值h,选取距离boundingbox最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中
x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将此3个中心点利用最小二乘法拟合为一条直线,此直线为股骨解剖轴线。
或者,所述解剖轴线获取模块执行以下步骤:提取胫骨三维模型的boundingbox的高度值h,选取距离boundingbox最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中
x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将此3个中心点利用最小二乘法拟合为一条直线,此直线为胫骨解剖轴线。
以下更详细地说明本发明。
本发明提供了一种腿骨下肢力线自动检测的方法,该方法可用于但不限于下肢腿骨的力线确定。
具体地,如图4所示,该方法包括以下步骤:
S101,获取CT影像数据;
利用CT设备对病人的预定部位进行术前扫描,以获取术前CT影像。预定部位例如人体下肢。
S102,对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型。
S103,针对股骨三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置。股骨头球三维模型不是规则的球形,检测非规则球心的方法实施如下:在股骨头球三维模型表面任选一点作为初始点位置,第一步过此点沿着横断面剖切,得到横断剖切面,提取剖切面闭合边缘,然后利用迭代循环搜索方法计算剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标即为横断剖切面的质心点坐标;第二步过横断剖切面质心点沿着矢状面剖切,得到矢状剖切面,然后同上方法获取矢状剖切面的质心点坐标;第三步,矢状剖切面质心点即为股骨头球心点坐标。
S104,针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标。详述如下:选择一例腿骨CT作为标准模板,在模板上标注3个中心点的位置,将模板CT配准到待测CT上,获得配准矩阵,模板上的3个点乘以配准矩阵,即可获取待测CT上的3个点位置坐标,重建出待测CT的腿骨三维模型,即可直观观察到3个生理特征点在腿骨模型的位置。
S105,针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定解剖轴线。详述如下:按照医学上对股骨解剖轴线的定义自动选择靠近股骨远端面的3个横断切面,方法如下:提取股骨三维模型的boundingbox的高度值h,选取距离boundingbox最顶端x1,x2,x3(x1/h=10%,x2/h=20%,x3/h=30%)高度的横断切面作为医学定义的3个横断切面,获取3个横断切面后,利用S103中的迭代循环搜索方法确定横断切面的中心点坐标,再将此3个中心点利用最小二乘法拟合为一条直线,此直线即为股骨解剖轴线;胫骨解剖轴线确定方法同上。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属本发明技术方案的保护范围。

Claims (8)

1.一种腿骨下肢力线自动检测方法,其特征在于:其包括以下步骤:
(1)获取CT影像数据:利用CT设备对病人的预定部位进行术前扫描,以获取术前CT影像;
(2)对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;
(3)针对股骨三维模型的股骨头球三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;
(4)针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;
(5)通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;
(6)针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线;
所述步骤(3)中,在股骨头球三维模型表面任选一点作为初始点位置,第一步过此初始点沿着横断面剖切,得到横断剖切面,提取横断剖切面闭合边缘,然后利用迭代循环搜索方法计算横断剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为横断剖切面的质心点坐标;第二步过横断剖切面质心点沿着矢状面剖切,得到矢状剖切面,提取矢状剖切面闭合边缘,然后利用迭代循环搜索方法计算矢状剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为矢状剖切面的质心点坐标;第三步,矢状剖切面质心点坐标作为股骨头球心点坐标。
2.根据权利要求1所述的腿骨下肢力线自动检测方法,其特征在于:所述步骤(4)中,选择一例腿骨CT影像作为标准模板,在标准模板上标注3个中心点的位置,3个中心点为膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点,将标准模板配准到待测CT影像上,获得配准矩阵,标准模板上的3个中心点乘以配准矩阵,来获取待测CT影像上的3个中心点位置坐标,重建出待测CT影像的腿骨三维模型,以便直观观察到3个中心点在腿骨模型的位置。
3.根据权利要求2所述的腿骨下肢力线自动检测方法,其特征在于:所述步骤(6)中,提取股骨三维模型的边界框(boundingbox)的高度值h,选取距离边界框(boundingbox)最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中
x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将确定的横断切面的3个中心点利用最小二乘法拟合为一条直线,此直线为股骨解剖轴线。
4.根据权利要求2所述的腿骨下肢力线自动检测方法,其特征在于:所述步骤(6)中,提取胫骨三维模型的边界框(boundingbox)的高度值h,选取距离边界框(boundingbox)最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中
x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将确定的横断切面的3个中心点利用最小二乘法拟合为一条直线,此直线为胫骨解剖轴线。
5.一种腿骨下肢力线自动检测装置,其特征在于:该装置包括:CT影像数据获取模块,其配置来对病人的预定部位进行术前扫描,以获取术前CT影像;
骨骼三维模型获取模块,其配置来对CT影像中的病人骨骼进行分割和三维重建,在模型空间中建立病人的骨骼三维模型;
股骨头球心点获取模块,其配置来针对股骨三维模型的股骨头球三维模型,通过双剖切面获取质心的方法来确定股骨头球心点位置;
配准模块,其配置来针对股骨和胫骨三维模型,采用刚性+弹性的配准方法获取膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点的坐标;
机械轴线获取模块,其配置来通过股骨头球心与膝关节髁间凹中心点的连线获得股骨机械轴线,通过胫骨平台中心点和踝关节中心点的连线获得胫骨机械轴线;
解剖轴线获取模块,其配置来针对股骨或胫骨三维模型,采用迭代搜索最大距离点对的方法确定股骨或胫骨的解剖轴线;
所述股骨头球心点获取模块执行以下步骤:在股骨头球三维模型表面任选一点作为初始点位置,第一步过此初始点沿着横断面剖切,得到横断剖切面,提取横断剖切面闭合边缘,然后利用迭代循环搜索方法计算横断剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为横断剖切面的质心点坐标;第二步过横断剖切面质心点沿着矢状面剖切,得到矢状剖切面,提取矢状剖切面闭合边缘,然后利用迭代循环搜索方法计算矢状剖切面闭合边缘上任意两点的距离值,取所有点对中距离值最大的点对,然后获取这个点对的中点坐标,此中点坐标作为矢状剖切面的质心点坐标;第三步,矢状剖切面质心点坐标作为股骨头球心点坐标。
6.根据权利要求5所述的腿骨下肢力线自动检测装置,其特征在于:所述配准模块执行以下步骤:选择一例腿骨CT影像作为标准模板,在标准模板上标注3个中心点的位置,3个中心点为膝关节髁间凹中心点、胫骨平台中心点、踝关节中心点,将标准模板配准到待测CT影像上,获得配准矩阵,标准模板上的3个中心点乘以配准矩阵,来获取待测CT影像上的3个中心点位置坐标,重建出待测CT影像的腿骨三维模型,以便直观观察到3个中心点在腿骨模型的位置。
7.根据权利要求6所述的腿骨下肢力线自动检测装置,其特征在于:所述解剖轴线获取模块执行以下步骤:提取股骨三维模型的边界框(boundingbox)的高度值h,选取距离边界框(boundingbox)最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将确定的横断切面的3个中心点利用最小二乘法拟合为一条直线,此直线为股骨解剖轴线。
8.根据权利要求6所述的腿骨下肢力线自动检测装置,其特征在于:所述解剖轴线获取模块执行以下步骤:提取胫骨三维模型的边界框(boundingbox)的高度值h,选取距离边界框(boundingbox)最顶端x1,x2,x3高度的横断切面作为医学定义的3个横断切面,其中x1/h=10%,x2/h=20%,x3/h=30%,获取3个横断切面后,利用迭代循环搜索方法确定横断切面的中心点坐标,再将确定的横断切面的3个中心点利用最小二乘法拟合为一条直线,此直线为胫骨解剖轴线。
CN201910884469.9A 2019-09-18 2019-09-18 一种腿骨下肢力线自动检测方法及装置 Active CN110613469B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910884469.9A CN110613469B (zh) 2019-09-18 2019-09-18 一种腿骨下肢力线自动检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910884469.9A CN110613469B (zh) 2019-09-18 2019-09-18 一种腿骨下肢力线自动检测方法及装置

Publications (2)

Publication Number Publication Date
CN110613469A CN110613469A (zh) 2019-12-27
CN110613469B true CN110613469B (zh) 2020-09-15

Family

ID=68923333

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910884469.9A Active CN110613469B (zh) 2019-09-18 2019-09-18 一种腿骨下肢力线自动检测方法及装置

Country Status (1)

Country Link
CN (1) CN110613469B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111311484B (zh) * 2020-02-25 2023-04-07 禹宝庆 一种基于人工智能技术的二维骨骼x光图像的三维转换方法
CN113545956B (zh) * 2021-06-01 2023-10-20 深圳大学 一种个性化定制的高顺应性膝关节外骨骼设计方法
CN113689406B (zh) * 2021-08-24 2022-04-08 北京长木谷医疗科技有限公司 基于运动模拟算法的膝关节股骨后髁点识别方法及系统
CN113974920B (zh) * 2021-10-08 2022-10-11 北京长木谷医疗科技有限公司 膝关节股骨力线确定方法和装置、电子设备、存储介质
CN113907774A (zh) * 2021-10-13 2022-01-11 瓴域影诺(北京)科技有限公司 一种测量下肢力线的方法及装置
CN114098707B (zh) * 2021-10-29 2024-08-16 西安交通大学 一种下肢力线角度的测量系统及测量方法
CN113870261B (zh) * 2021-12-01 2022-05-13 杭州柳叶刀机器人有限公司 用神经网络识别力线的方法与系统、存储介质及电子设备
CN114305697B (zh) * 2021-12-25 2024-08-09 山东航维骨科医疗器械股份有限公司 一种关节置换手术机器人的控制方法
CN115115694B (zh) * 2022-03-03 2023-06-13 中国人民解放军总医院第四医学中心 人体长骨中轴的建立方法、装置、设备及存储介质
CN115005977A (zh) * 2022-05-20 2022-09-06 长春理工大学 一种膝关节置换手术术前规划方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5108451A (en) * 1991-01-31 1992-04-28 Forte Mark R Femoral component of a hip joint prosthesis
FR2841768B1 (fr) * 2002-07-05 2005-05-06 Tornier Sa Prothese d'epaule ou de hanche facilitant l'abduction
CN1179707C (zh) * 2003-08-21 2004-12-15 上海交通大学 机器人全膝关节置换术的股骨定位方法
CN103476363B (zh) * 2011-02-15 2017-06-30 康复米斯公司 改进的适合患者型关节植入物以及处理、评估、校正、修改和/或适应解剖变化和/或非对称性的手术和工具
CN104970904A (zh) * 2014-04-14 2015-10-14 陆声 基于mri的全膝关节置换个体化定位模板设计
CN105139442A (zh) * 2015-07-23 2015-12-09 昆明医科大学第一附属医院 一种结合ct和mri二维图像建立人体膝关节三维仿真模型的方法
CN106963489B (zh) * 2017-05-12 2019-05-10 常州工程职业技术学院 一种个体化股骨骨折复位模型构建方法
CN108042217A (zh) * 2017-12-21 2018-05-18 成都真实维度科技有限公司 一种三维空间下肢生物力线的确定方法

Also Published As

Publication number Publication date
CN110613469A (zh) 2019-12-27

Similar Documents

Publication Publication Date Title
CN110613469B (zh) 一种腿骨下肢力线自动检测方法及装置
JP6862506B2 (ja) 動態位置合わせ並びに新規な大腿骨人工装具及び脛骨人口装具
US9681956B2 (en) Acquiring and utilizing kinematic information for patient-adapted implants, tools and surgical procedures
CN104970904A (zh) 基于mri的全膝关节置换个体化定位模板设计
CN105361883A (zh) 膝关节置换三维空间下肢生物力线的确定方法
CN108478250A (zh) 全膝关节置换术的股骨定位方法、个体化截骨导板及假体
CN111329584B (zh) 外科规划和方法
CN107296651A (zh) 一种数字化辅助确定股骨远端截骨平面的方法
US20240261031A1 (en) Systems and methods for pre-operative visualization of a joint
JP2016187619A (ja) カスタマイズされた整形外科用インプラント及び関連方法並びに変形可能な関節テンプレート
US20150250552A1 (en) Advanced methods of modeling knee joint kinematics and designing surgical repair systems
US9925069B2 (en) Method for femur resection alignment approximation in hip replacement procedures
CN113332008B (zh) 一种髋关节置换术中股骨长度的确定方法及系统
CN114711796A (zh) 双下肢负重位力线测量及诊断系统及智能化测量方法
CN104905869A (zh) 用于椎弓根定位及实时监测钉道的导航模板及其制作方法
US20240277419A1 (en) Patella Tracking
CN105877807A (zh) 一种能够准确定位胫骨高位截骨的导航装置及其制备方法
CN104997562A (zh) 一种能够准确定位半肩置换的导航装置及其制作方法
CN110680353B (zh) 一种股骨头完整率的分析装置
CN115607286B (zh) 基于双目标定的膝关节置换手术导航方法、系统及设备
CN108175378A (zh) 一种膝关节中心点的确定方法
US20240398476A1 (en) Measurement and diagnostic system for weight-bearing lines of lower limbs and intelligent measurement method thereof
CN109106447A (zh) 一种用于x光下的获取股骨标准侧位和标准髁间窝图的方法
CN107913100A (zh) 一种股骨头中心点的确定方法
RU2689032C1 (ru) Способ ацетабулопластики у детей с церебральным параличом

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