[go: up one dir, main page]

CN111888059B - 基于深度学习与x线的全髋关节图像处理方法及装置 - Google Patents

基于深度学习与x线的全髋关节图像处理方法及装置 Download PDF

Info

Publication number
CN111888059B
CN111888059B CN202010707817.8A CN202010707817A CN111888059B CN 111888059 B CN111888059 B CN 111888059B CN 202010707817 A CN202010707817 A CN 202010707817A CN 111888059 B CN111888059 B CN 111888059B
Authority
CN
China
Prior art keywords
image
hip joint
determining
femoral
neural network
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
CN202010707817.8A
Other languages
English (en)
Other versions
CN111888059A (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.)
Longwood Valley Medtech Co Ltd
Original Assignee
Beijing Changmugu Medical 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 Beijing Changmugu Medical Technology Co Ltd filed Critical Beijing Changmugu Medical Technology Co Ltd
Publication of CN111888059A publication Critical patent/CN111888059A/zh
Application granted granted Critical
Publication of CN111888059B publication Critical patent/CN111888059B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/46Special tools for implanting artificial joints
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61FFILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F2/00Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
    • A61F2/02Prostheses implantable into the body
    • A61F2/30Joints
    • A61F2/46Special tools for implanting artificial joints
    • A61F2002/4632Special tools for implanting artificial joints using computer-controlled surgery, e.g. robotic surgery
    • A61F2002/4633Special tools for implanting artificial joints using computer-controlled surgery, e.g. robotic surgery for selection of endoprosthetic joints or for pre-operative planning

Landscapes

  • Health & Medical Sciences (AREA)
  • Transplantation (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Vascular Medicine (AREA)
  • Cardiology (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本申请公开了一种基于深度学习与X线的全髋关节图像处理方法及装置,本申请的方法包括获取髋关节的X线图像,所述髋关节的X线图像中包含参照物的图像,所述参照物为已知尺寸的参照物;根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原;基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号;根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。本申请是为了提出一种更便捷更准确的术前规划的方式为全髋关节置换手术提供更好的术前支持。

Description

基于深度学习与X线的全髋关节图像处理方法及装置
技术领域
本申请涉及医学技术领域,具体而言,涉及一种基于深度学习与X线的全髋关节图像处理方法及装置。
背景技术
在医学领域中全髋关节置换手术的术前规划主要包括计算所需假体型号及截骨线位置,全髋关节置换手术的术前规划对于手术的成功率起着非常重要的作用,因此如何提供准确的术前规划是非常重要的。目前主要的术前规划方式为人工通过各种工具进行测量,效率低而且准确性无法保证,因此亟需提供一种更便捷更准确的术前规划的方法为全髋关节置换手术提供更好的术前支持。
发明内容
本申请的主要目的在于提出一种基于深度学习与X线的全髋关节图像处理方法及装置,以提供一种更便捷更准确的术前规划的方式为全髋关节置换手术提供更好的术前支持。
为了实现上述目的,根据本申请的第一方面,提供了一种基于深度学习与X线的全髋关节图像处理方法。
根据本申请的基于深度学习与X线的全髋关节图像处理方法包括:
获取髋关节的X线图像,所述髋关节的X线图像中包含参照物的图像,所述参照物为已知尺寸的参照物;
根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原;
基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号;
根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。
可选的,所述基于深度学习模型对还原后的髋关节的X线图像进行识别确定腿长差,包括:
将髋关节的X线图像转化为灰度图;
基于第一神经网络模型对灰度图的每个像素值进行预测,确定泪滴关键点以及股骨小转子关键点位置;
根据泪滴关键点以及股骨小转子关键点位置,确定腿长差。
可选的,所述基于深度学习模型对还原后的髋关节的X线图像进行识别确定髋臼杯位置包括:
将髋关节的X线图像转化为灰度图;
基于第二神经网络模型对灰度图的每个像素值进行预测,确定股骨头位置;
根据平面图像的质心公式计算股骨头旋转中心;
根据股骨头的直径推算髋臼杯直径;
根据骨头旋转中心和髋臼杯直径确定髋臼杯位置。
可选的,所述基于深度学习模型对还原后的髋关节的X线图像进行识别确定股骨柄假体的规格型号包括:
将髋关节的X线图像转化为灰度图;
基于第三神经网络模型对灰度图进行识别,确定髓腔解剖轴线;
基于第四神经网络模型对灰度图进行识别,确定股骨颈中心轴线;
根据髓腔解剖轴线和股骨颈中心轴线确定股骨颈干角;
根据股骨颈干角、在确定髓腔解剖轴线过程中确定的髓腔区域以及股骨头旋转中心确定股骨柄假体的规格型号。
可选的,所述根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置包括:
将股骨柄假体的旋转中心与髋臼杯的旋转中心位置重合,确定股骨柄假体实际位置;
沿股骨柄假体的涂层位置确定截骨线位置。
可选的,所述基于第三神经网络模型对灰度图进行识别,确定髓腔解剖轴线包括:
基于第三神经网络模型对灰度图的每个像素值进行预测,确定股骨头区域和骨皮质区域;
根据股骨头区域、骨皮质区域确定髓腔区域;
对髓腔区域多个中心点坐标进行直线拟合确定髓腔解剖轴线。
可选的,所述基于第四神经网络模型对灰度图进行识别,确定股骨颈中心轴线包括:
基于第四神经网络模型对灰度图的每个像素值进行预测,确定股骨头区域和股骨颈基底区域;
根据平面图像的质心公式计算股骨头区域和股骨颈基底区域对应的股骨头中心坐标和股骨颈基底中心坐标;
根据股骨头中心坐标和股骨颈基底中心坐标确定股骨颈中心轴线。
可选的,根据截骨线位置,计算术后的腿长差,以及偏距。
为了实现上述目的,根据本申请的第二方面,提供了一种基于深度学习与X线的全髋关节图像处理装置。
根据本申请的基于深度学习与X线的全髋关节图像处理装置包括:
比例校准单元,用于根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的真实还原;
腿长差确定单元,用于基于第一神经网络模型对还原后的髋关节的X线图像进行识别,确定腿长差;
髋臼杯位置确定单元,用于基于第二神经网络模型对还原后的髋关节的X线图像进行识别,确定髋臼杯位置;
股骨柄假体规格确定单元,用于基于第三神经网络模型、第四神经网络模型对还原后的髋关节的X线图像进行识别,确定股骨柄假体规格;
截骨线确定单元,用于根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。
为了实现上述目的,根据本申请的第三方面,提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于使所述计算机执行上述第一方面中任意一项所述的基于深度学习与X线的全髋关节图像处理方法。
为了实现上述目的,根据本申请的第四方面,提供了一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的计算机程序,所述计算机程序被所述至少一个处理器执行,以使所述至少一个处理器执行上述第一方面中任意一项所述的基于深度学习与X线的全髋关节图像处理方法。
在本申请实施例中,基于深度学习与X线的全髋关节图像处理方法及装置中,获取髋关节的X线图像,所述髋关节的X线图像中包含参照物的图像,所述参照物为已知尺寸的参照物;根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原;基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号;根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。可以看出,本实施例的全髋关节置换术前规划方式中,将髋关节的X线图像进行了真实尺寸的还原,以实际的尺寸进行后续的位置识别更准确;另外,在对X线图像识别的过程中都是基于深度学习模型进行识别的,进一步的保证了根据识别结果确定的腿长差、髋臼杯位置、股骨柄假体的规格型号、截骨线位置的准确性和快速性,从而为全髋关节置换手术提供了更好的术前支持。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解,使得本申请的其它特征、目的和优点变得更明显。本申请的示意性实施例附图及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是根据本申请实施例提供的一种基于深度学习与X线的全髋关节图像处理方法流程图;
图2是根据本申请实施例提供的一种为髋关节的X线图像的示意图;
图3-4是根据本申请实施例提供的一种确定临床中的实际截骨线位置的示意图;
图5是根据本申请实施例提供的一种确定腿长差的方法流程图;
图6是根据本申请实施例提供的一种自动识别出泪滴关键点以及股骨小转子关键点位置的示意图;
图7是根据本申请实施例提供的一种腿长差确定的示意图;
图8是根据本申请实施例提供的一种确定髋臼杯位置的方法流程图;
图9为根据本申请实施例提供的一种识别股骨头的示意图;
图10为根据本申请实施例提供的股骨头旋转中心的示意图;
图11是根据本申请实施例提供的髋臼杯位置的示意图;
图12是根据本申请实施例提供的一种确定股骨柄假体的规格型号的方法流程图;
图13是根据本申请实施例提供的一种识别股骨头区域、骨皮质区域的示意图;
图14是根据本申请实施例提供的一种髓腔区域的示意图;
图15是根据本申请实施例提供的一种确定髓腔解剖轴线的示意图;
图16是根据本申请实施例提供的一种识别股骨头区域、股骨颈基底区域的示意图;
图17是根据本申请实施例提供的一种股骨颈中心轴线的示意图;
图18是根据本申请实施例提供的一种基于深度学习与X线的全髋关节图像处理装置的组成框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分的实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本申请的实施例。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
根据本申请实施例,提供了一种基于深度学习与X线的全髋关节图像处理方法,如图1所示,该方法包括如下的步骤:
S101.获取髋关节的X线图像,髋关节的X线图像中包含参照物的图像。
髋关节的X线图像是通过对髋关节进行X光片拍摄时获取的,同时在同张照片里拍摄一个已知尺寸的物体,即参照物。因此得到的髋关节的X线图像中包含参照物的图像。如图2所示,为髋关节的X线图像,其中图像的底部中心部位的标示标准尺寸的图像为参照物的图像。是实际应用中,参照物的选取和拍摄时的排放位置可以根据适应性的调整,本实施例不作限制。
S102.根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原。
参照物的尺寸是已知的,参照物的图像尺寸也可以通过测量得到,根据参照物的图像尺寸及其实际尺寸的比例,可以确定出髋关节的X线图像相对于实际的髋关节尺寸的比例(两者比例相同),然后根据比例将髋关节的X线图像进行真实尺寸的还原。将髋关节的X线图像进行真实尺寸的还原是为和后续的图像识别做基础,使后续根据识别结果确定的腿长差、髋臼杯位置、股骨柄假体的规格型号、截骨线位置与实际的对应位置差距更小,保证识别的准确性。
具体的还原操作可以为选取已知尺寸物体的关键部位尺寸。通过计算图像中像素间两点距离,和物体实际尺寸进行比例换算,确定比例,然后根据比例对髋关节的X线图像的比例进行修正。
S103.基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号。
深度学习模型是神经网络模型,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号可能会用到的模型的输入和输出可能是不同的,但是模型训练的原理是相同的。具体的,神经网络模型训练的原理为:将髋关节的X线图像转化为0-255灰度图,然后将图像进行人工选定标注,将图片的每个像素标注划分为几种属性值(根据实际的需求,属性值的种类数不同,比如可以为两种、三种等等)并分别命名,然后将其输入到神经网络模型中进行卷积池化采样一直迭代学习训练得到神经网络模型。
本步骤中的神经网络模型为分类神经网络,是将图像中的不同的区域进行分类,比如在确定腿长差时,应用神经网络模型主要是为了识别出泪滴和股骨小转子的关键点;再比如在确定髋臼杯位置时,应用神经网络模型主要是为了识别出股骨头区域;再比如在确定股骨柄假体的规格型号时,应用神经网络模型主要是为了识别出股骨头、骨皮质区域以及股骨头、股骨颈基底区域。
本实施例中的神经网络可以为卷积神经网络LeNet、卷积神经网络AlexNet、可视化卷积神经网络ZF-Net、卷积神经网络GoogleNet、卷积神经网络VGG、卷积神经网络Inception、卷积神经网络ResNet、卷积神经网络DensNet、卷积神经网络Inception ResNet等。
确定腿长差、髋臼杯位置以及股骨柄假体的规格型号是根据图像的识别结果再进行一些坐标、拟合等计算后确定的。
S104.根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。
具体的,“根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置”为通过移动股骨柄假体,将股骨柄假体的旋转中心与之前计算的髋臼杯旋转中心位置重合,得到股骨柄假体实际位置。沿股骨柄假体的涂层位置可确定临床中的实际截骨线位置,如图3-4所示。图3为移动股骨柄假体到预定位置,使股骨柄假体的旋转中心与之前计算的髋臼杯旋转中心位置重合,图4为根据股骨柄假体的外形确定截骨线位置。
从以上的描述中,可以看出,本申请实施例的基于深度学习与X线的全髋关节图像处理方法中,获取髋关节的X线图像,所述髋关节的X线图像中包含参照物的图像,所述参照物为已知尺寸的参照物;根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原;基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号;根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。可以看出,本实施例的全髋关节置换术前规划方式中,将髋关节的X线图像进行了真实尺寸的还原,以实际的尺寸进行后续的位置识别更准确;另外,在对X线图像识别的过程中都是基于深度学习模型进行识别的,进一步的保证了根据识别结果确定的腿长差、髋臼杯位置、股骨柄假体的规格型号、截骨线位置的准确性和快速性,从而为全髋关节置换手术提供了更好的术前支持。
进一步的,作为上述实施例的进一步细化,步骤S103对于确定腿长差、髋臼杯位置以及股骨柄假体的规格型号的详细步骤进行分别说明。
如图5所示,为确定腿长差的流程图,具体包括如下步骤:
S201.将髋关节的X线图像转化为灰度图。
将髋关节的X线图像转化为0-255灰度图。
S202.基于第一神经网络模型对灰度图的每个像素值进行预测,确定泪滴关键点以及股骨小转子关键点位置。
在进行预测之前,首先要根据样本训练得到第一神经网络模型。具体的,将未标记的原始图像(髋关节的X线样本图像对应的灰度图)以及人工识别标记的泪滴关键点及股骨小转子关键点位置的标记传入到卷积神经网络中,将输入的原始图像与特征点的高斯分布函数进行拟合,进行卷积池化采样一直迭代学习训练得到第一神经网络模型。需要说明的是,本步骤中的卷积神经网络可以为卷积神经网络LeNet、卷积神经网络AlexNet、可视化卷积神经网络ZF-Net、卷积神经网络GoogleNet、卷积神经网络VGG、卷积神经网络Inception、卷积神经网络ResNet、卷积神经网络DensNet、卷积神经网络Inception ResNet等。
得到第一神经网络模型后,将髋关节的X线图像对应的灰度图输入到第一神经网络模型中,可以自动识别出泪滴关键点以及股骨小转子关键点位置。如图6所示,图6为自动识别出泪滴关键点以及股骨小转子关键点位置的示意图。
S203.根据泪滴关键点以及股骨小转子关键点位置,确定腿长差。
具体的,如图7所示,其中水平的直线是由两个泪滴关键点确定的,是两个泪滴关键点的连线,其中两条垂直的线段是由股骨小转子关键点和水平直线确定的,两条垂直直线分别记作A和B,A和B的差值为腿长差。
如图8所示,为确定髋臼杯位置的流程图,具体包括如下步骤:
S301.将髋关节的X线图像转化为灰度图。
将髋关节的X线图像转化为0-255灰度图。
S302.基于第二神经网络模型对灰度图的每个像素值进行预测,确定股骨头位置。
在进行预测之前,首先要根据样本训练得到第二神经网络模型。具体的,将未标记的原始图像(髋关节的X线样本图像对应的灰度图)以及人工识别标记像素属性值的标记图像传入到卷积神经网络中,包括两种属性值,分别命名0、1。数值0代表背景像素,1代表股骨头像素;传入到卷积神经网络中,进行卷积池化采样一直迭代学习训练得到第二神经网络模型。需要说明的是,本步骤中的卷积神经网络可以为卷积神经网络LeNet、卷积神经网络AlexNet、可视化卷积神经网络ZF-Net、卷积神经网络GoogleNet、卷积神经网络VGG、卷积神经网络Inception、卷积神经网络ResNet、卷积神经网络DensNet、卷积神经网络InceptionResNet等。
得到第二神经网络模型后,将髋关节的X线图像对应的灰度图输入到第二神经网络模型中,可以对每个像素值进行预测。自动将X线图像的每个像素值归为一个属性中:0-背景,1-股骨头,完成股骨头区域(即股骨头位置)的自动识别,如图9所示。图9为识别股骨头的示意图。
S303.根据平面图像的质心公式计算股骨头旋转中心。
因为得到的股骨头区域的图像是二值图像,其质量分布是均匀的,所以质心和形心重合,根据平面图像的质心公式可以计算得到股骨头的中心点坐标,即股骨头旋转中心。假设二值图像为B[i,j],则可根据下列公式求得股骨头的中心点坐标:
Figure GDA0002952170040000101
其中:
Figure GDA0002952170040000102
此处得到的是股骨头的中心点的像素坐标,需要将像素坐标转换为图像坐标。图像平面坐标中心坐标为:
Figure GDA0002952170040000111
Figure GDA0002952170040000112
则像素坐标
Figure GDA0002952170040000113
到图像坐标(x’,y’)的变换公式为:
Figure GDA0002952170040000114
Figure GDA0002952170040000115
其中Sx,Sy分别为图像阵列的行列间距。最后通过输出显示模块,得到股骨头旋转中心的位置,如图9所示。图10中圆圈的中心点为股骨头旋转中心。
S304.根据股骨头的直径推算髋臼杯直径。
根据股骨头区域和股骨头旋转中心确定股骨头的直径,根据股骨头的直径推算髋臼杯直径。根据股骨头的直径推算髋臼杯直径的可以参考现有的任意一种推算方式确定髋臼杯直径。
S305.根据骨头旋转中心和髋臼杯直径确定髋臼杯位置。
根据股骨头的直径以及股骨头旋转中心位置自动确定髋臼杯位置,如图11所示。图11中线条勾画的区域为髋臼杯位置。
如图12所示,为确定股骨柄假体的规格型号的流程图,具体包括如下步骤:
S401.将髋关节的X线图像转化为灰度图。
将髋关节的X线图像转化为0-255灰度图。
S402.基于第三神经网络模型对灰度图进行识别,确定髓腔解剖轴线。
具体的,确定髓腔解剖轴线包括如下步骤:
首先,基于第三神经网络模型对灰度图的每个像素值进行预测,确定股骨头区域和骨皮质区域;
在进行预测之前,首先要根据样本训练得到第三神经网络模型。具体的,将未标记的原始图像(髋关节的X线样本图像对应的灰度图)以及人工识别标记像素属性值的标记图像传入到卷积神经网络中,包括三种属性值,分别命名0、1、2。数值0代表背景像素,1代表股骨头像素,2代表骨皮质;传入到卷积神经网络中,进行卷积池化采样一直迭代学习训练得到第三神经网络模型。需要说明的是,本步骤中的卷积神经网络可以为卷积神经网络LeNet、卷积神经网络AlexNet、可视化卷积神经网络ZF-Net、卷积神经网络GoogleNet、卷积神经网络VGG、卷积神经网络Inception、卷积神经网络ResNet、卷积神经网络DensNet、卷积神经网络Inception ResNet等。
得到第三神经网络模型后,将髋关节的X线图像对应的灰度图输入到第三神经网络模型中,可以对每个像素值进行预测。自动将X线图像的每个像素值归为一个属性中:0-背景,1-股骨头,2-骨皮质,完成股骨头区域、骨皮质区域的自动识别,如图13所示。图13为识别股骨头区域、骨皮质区域的示意图。
其次,根据股骨头区域、骨皮质区域确定髓腔区域;
具体的,截取小转子结束处直到股骨末端部位,使用图像中股骨区域减去骨皮质区域得到的是髓腔区域,如图14所示。
最后,对髓腔区域多个中心点坐标进行直线拟合确定髓腔解剖轴线。
具体的,如图15所示,从小转子结束位置以下,每横行与髓腔交点为四个坐标,从左至右分别命名为A1,A2,B1,B2;依据两点可以求出中点,A1(X1,Y1),A2(X2,Y2)的中点坐标:
Figure GDA0002952170040000121
B1,B2同理可算得。每行依次算得髓腔的中点坐标,将这些点拟合成一条直线即为髓腔解剖轴线(也是股骨解剖轴线)。
S403.基于第四神经网络模型对灰度图进行识别,确定股骨颈中心轴线。
具体的,确定股骨颈中心轴线包括如下步骤:
首先,基于第四神经网络模型对灰度图的每个像素值进行预测,确定股骨头区域和股骨颈基底区域;
在进行预测之前,首先要根据样本训练得到第四神经网络模型。具体的,将未标记的原始图像(髋关节的X线样本图像对应的灰度图)以及人工识别标记像素属性值的标记图像传入到卷积神经网络中,包括三种属性值,分别命名0、1、2。数值0代表背景像素,1代表股骨头像素,2代表股骨颈基底像素;传入到卷积神经网络中,进行卷积池化采样一直迭代学习训练得到第四神经网络模型。需要说明的是,本步骤中的卷积神经网络可以为卷积神经网络LeNet、卷积神经网络AlexNet、可视化卷积神经网络ZF-Net、卷积神经网络GoogleNet、卷积神经网络VGG、卷积神经网络Inception、卷积神经网络ResNet、卷积神经网络DensNet、卷积神经网络Inception ResNet等。
在得到第四神经网络模型后,将髋关节的X线图像对应的灰度图输入到第四神经网络模型中,可以对每个像素值进行预测。自动将X线图像的每个像素值归为一个属性中:0-背景,1-股骨头,2-股骨颈基底像素,完成股骨头区域、股骨颈基底区域的自动识别,如图16所示。图16为识别股骨头区域、股骨颈基底区域的示意图。
其次,根据平面图像的质心公式计算股骨头区域和股骨颈基底区域对应的股骨头中心坐标和股骨颈基底中心坐标;
股骨头中心坐标和股骨颈基底中心坐标的计算方式类似,都可以参见步骤S303中计算股骨头中心点坐标的实现方式,此处不在赘述。
最后,根据股骨头中心坐标和股骨颈基底中心坐标确定股骨颈中心轴线。
具体的,股骨头中心坐标和股骨颈基底中心坐标连线即为股骨颈中心轴线,如图17所示。图17中两条斜向下的线段为股骨颈中心轴线。
S404.根据髓腔解剖轴线和股骨颈中心轴线确定股骨颈干角。
具体的,髓腔解剖轴线和股骨颈中心轴线形成的夹角为股骨颈干角。
S405.根据股骨颈干角、在确定髓腔解剖轴线过程中确定的髓腔区域以及股骨头旋转中心确定股骨柄假体的规格型号。
具体的,根据股骨颈干角角度值,再结合髓腔形态,股骨头旋转中心位置可对股骨柄假体型号的选择给出推荐。股骨柄假体型号按照股骨柄假体的形状和尺寸等特征进行区分。
进一步的,作为图1实施例的补充说明,在确定截骨线位置之后,还包括根据截骨线位置,计算术后的腿长差,以及偏距。具体的,偏距包含股骨偏心距:指股骨头旋转中心至股骨干长轴间的垂直距离。还包括联合偏距,具体为股骨和髋臼偏距的累积和。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
根据本申请实施例,还提供了一种用于实施上述图1-17所述方法的基于深度学习与X线的全髋关节图像处理装置,如图18所示,该装置包括:
比例校准单元51,用于根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的真实还原;
腿长差确定单元52,用于基于第一神经网络模型对还原后的髋关节的X线图像进行识别,确定腿长差;
髋臼杯位置确定单元53,用于基于第二神经网络模型对还原后的髋关节的X线图像进行识别,确定髋臼杯位置;
股骨柄假体规格确定单元54,用于基于第三神经网络模型、第四神经网络模型对还原后的髋关节的X线图像进行识别,确定股骨柄假体规格;
截骨线确定单元55,用于根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。
具体的,本申请实施例的装置中各单元、模块实现其功能的具体过程可参见方法实施例中的相关描述,此处不再赘述。
从以上的描述中,可以看出,本申请实施例的基于深度学习与X线的全髋关节图像处理装置中,获取髋关节的X线图像,所述髋关节的X线图像中包含参照物的图像,所述参照物为已知尺寸的参照物;根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原;基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号;根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置。可以看出,本实施例的全髋关节置换术前规划方式中,将髋关节的X线图像进行了真实尺寸的还原,以实际的尺寸进行后续的位置识别更准确;另外,在对X线图像识别的过程中都是基于深度学习模型进行识别的,进一步的保证了根据识别结果确定的腿长差、髋臼杯位置、股骨柄假体的规格型号、截骨线位置的准确性和快速性,从而为全髋关节置换手术提供了更好的术前支持。
根据本申请实施例,还提供了一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于使所述计算机执行上述方法实施例中的基于深度学习与X线的全髋关节图像处理方法。
根据本申请实施例,还提供了一种电子设备,包括:至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的计算机程序,所述计算机程序被所述至少一个处理器执行,以使所述至少一个处理器执行上述方法实施例中的基于深度学习与X线的全髋关节图像处理方法。
显然,本领域的技术人员应该明白,上述的本申请的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本申请不限制于任何特定的硬件和软件结合。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (9)

1.一种基于深度学习与X线的全髋关节图像处理方法,所述方法在计算机中执行,其特征在于,所述方法包括:
获取髋关节的X线图像,所述髋关节的X线图像中包含参照物的图像,所述参照物为已知尺寸的参照物;
比例校准单元根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的还原;
基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号;
所述基于深度学习模型对还原后的髋关节的X线图像进行识别,确定腿长差、髋臼杯位置以及股骨柄假体的规格型号包括:腿长差确定单元基于第一神经网络模型对还原后的髋关节的X线图像进行识别,确定腿长差;髋臼杯位置确定单元基于第二神经网络模型对还原后的髋关节的X线图像进行识别,确定髋臼杯位置;股骨柄假体规格确定单元基于第三神经网络模型、第四神经网络模型对还原后的髋关节的X线图像进行识别,确定股骨柄假体规格;
截骨线确定单元根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置;
所述根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置包括:
将股骨柄假体的旋转中心与髋臼杯的旋转中心位置重合,确定股骨柄假体实际位置;沿股骨柄假体的涂层位置确定截骨线位置。
2.根据权利要求1所述的基于深度学习与X线的全髋关节图像处理方法,其特征在于,所述基于深度学习模型对还原后的髋关节的X线图像进行识别确定腿长差,包括:
将髋关节的X线图像转化为灰度图;
基于第一神经网络模型对灰度图的每个像素值进行预测,确定泪滴关键点以及股骨小转子关键点位置;
根据泪滴关键点以及股骨小转子关键点位置,确定腿长差。
3.根据权利要求1所述的基于深度学习与X线的全髋关节图像处理方法,其特征在于,所述基于深度学习模型对还原后的髋关节的X线图像进行识别确定髋臼杯位置包括:
将髋关节的X线图像转化为灰度图;
基于第二神经网络模型对灰度图的每个像素值进行预测,确定股骨头位置;
根据平面图像的质心公式计算股骨头旋转中心;
根据股骨头的直径推算髋臼杯直径;
根据骨头旋转中心和髋臼杯直径确定髋臼杯位置。
4.根据权利要求1所述的基于深度学习与X线的全髋关节图像处理方法,其特征在于,所述基于深度学习模型对还原后的髋关节的X线图像进行识别确定股骨柄假体的规格型号包括:
将髋关节的X线图像转化为灰度图;
基于第三神经网络模型对灰度图进行识别,确定髓腔解剖轴线;
基于第四神经网络模型对灰度图进行识别,确定股骨颈中心轴线;
根据髓腔解剖轴线和股骨颈中心轴线确定股骨颈干角;
根据股骨颈干角、在确定髓腔解剖轴线过程中确定的髓腔区域以及股骨头旋转中心确定股骨柄假体的规格型号。
5.根据权利要求4所述的基于深度学习与X线的全髋关节图像处理方法,其特征在于,所述基于第三神经网络模型对灰度图进行识别,确定髓腔解剖轴线包括:
基于第三神经网络模型对灰度图的每个像素值进行预测,确定股骨头区域和骨皮质区域;
根据股骨头区域、骨皮质区域确定髓腔区域;
对髓腔区域多个中心点坐标进行直线拟合确定髓腔解剖轴线。
6.根据权利要求4所述的基于深度学习与X线的全髋关节图像处理方法,其特征在于,所述基于第四神经网络模型对灰度图进行识别,确定股骨颈中心轴线包括:
基于第四神经网络模型对灰度图的每个像素值进行预测,确定股骨头区域和股骨颈基底区域;
根据平面图像的质心公式计算股骨头区域和股骨颈基底区域对应的股骨头中心坐标和股骨颈基底中心坐标;
根据股骨头中心坐标和股骨颈基底中心坐标确定股骨颈中心轴线。
7.根据权利要求1所述的基于深度学习与X线的全髋关节图像处理方法,其特征在于,根据截骨线位置,计算术后的腿长差,以及偏距。
8.一种基于深度学习与X线的全髋关节图像处理装置,其特征在于,所述装置包括:
比例校准单元,用于根据参照物的图像尺寸及其实际尺寸的比例,将髋关节的X线图像进行尺寸的真实还原;
腿长差确定单元,用于基于第一神经网络模型对还原后的髋关节的X线图像进行识别,确定腿长差;
髋臼杯位置确定单元,用于基于第二神经网络模型对还原后的髋关节的X线图像进行识别,确定髋臼杯位置;
股骨柄假体规格确定单元,用于基于第三神经网络模型、第四神经网络模型对还原后的髋关节的X线图像进行识别,确定股骨柄假体规格;
截骨线确定单元,用于根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置;
所述根据股骨柄假体的旋转中心和在髋关节的X线图像识别过程中确定的髋臼杯的旋转中心确定截骨线位置包括:
将股骨柄假体的旋转中心与髋臼杯的旋转中心位置重合,确定股骨柄假体实际位置;沿股骨柄假体的涂层位置确定截骨线位置。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于使所述计算机执行权利要求1-7任意一项所述的基于深度学习与X线的全髋关节图像处理方法。
CN202010707817.8A 2020-07-06 2020-07-21 基于深度学习与x线的全髋关节图像处理方法及装置 Active CN111888059B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2020106437135 2020-07-06
CN202010643713 2020-07-06

Publications (2)

Publication Number Publication Date
CN111888059A CN111888059A (zh) 2020-11-06
CN111888059B true CN111888059B (zh) 2021-07-27

Family

ID=73190359

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010707817.8A Active CN111888059B (zh) 2020-07-06 2020-07-21 基于深度学习与x线的全髋关节图像处理方法及装置

Country Status (2)

Country Link
CN (1) CN111888059B (zh)
WO (1) WO2022007972A1 (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112842529B (zh) * 2020-12-31 2022-02-08 北京长木谷医疗科技有限公司 全膝关节图像处理方法及装置
CN113133802B (zh) * 2021-04-20 2022-12-23 四川大学 一种基于机器学习的骨手术线自动定点方法
CN113592821A (zh) * 2021-07-30 2021-11-02 瓴域影诺(北京)科技有限公司 一种股骨髓腔类型检测方法和系统
CN113744214B (zh) * 2021-08-24 2022-05-13 北京长木谷医疗科技有限公司 基于深度强化学习的股骨柄放置装置及电子设备
CN113974920B (zh) * 2021-10-08 2022-10-11 北京长木谷医疗科技有限公司 膝关节股骨力线确定方法和装置、电子设备、存储介质
CN113907775B (zh) * 2021-10-13 2024-12-06 瓴域影诺(北京)科技有限公司 一种髋关节影像质量判断方法及系统
CN114419618B (zh) * 2022-01-27 2024-02-02 北京长木谷医疗科技股份有限公司 基于深度学习的全髋关节置换术前规划系统
CN114742747B (zh) * 2022-02-24 2023-04-18 北京长木谷医疗科技有限公司 基于深度学习的髋关节置换术后影像的评估方法及系统
CN114431957B (zh) * 2022-04-12 2022-07-29 北京长木谷医疗科技有限公司 基于深度学习的全膝关节置换术后翻修术前规划系统
CN115830247B (zh) * 2023-02-14 2023-07-14 北京壹点灵动科技有限公司 髋关节旋转中心的拟合方法和装置、处理器及电子设备
CN116423054A (zh) * 2023-03-09 2023-07-14 中铁九桥工程有限公司 一种u肋板焊接方法及焊接系统
CN116597002B (zh) * 2023-05-12 2024-01-30 北京长木谷医疗科技股份有限公司 基于深度强化学习的股骨柄自动放置方法、装置及设备
CN116650110B (zh) * 2023-06-12 2024-05-07 北京长木谷医疗科技股份有限公司 基于深度强化学习的膝关节假体自动放置方法及装置
CN116993824B (zh) * 2023-07-19 2024-11-15 北京长木谷医疗科技股份有限公司 髋臼旋转中心计算方法、装置、设备及可读存储介质
CN117437459B (zh) * 2023-10-08 2024-03-22 昆山市第一人民医院 基于决策网络实现用户膝关节髌骨软化状态分析方法
CN118121299B (zh) * 2024-03-26 2024-11-15 北京和华瑞博医疗科技有限公司 数据处理方法、装置、设备、介质和程序产品
CN119385684B (zh) * 2024-12-27 2025-03-14 中国人民解放军总医院第四医学中心 一种髋关节置换术中髋关节旋转中心位置的预测方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101815477A (zh) * 2007-09-28 2010-08-25 株式会社力克赛 人工膝关节置换手术的术前计划装置及手术辅助用工具
CN103209652A (zh) * 2010-08-13 2013-07-17 史密夫和内修有限公司 外科手术引导件
CN106456196A (zh) * 2014-02-11 2017-02-22 史密夫和内修有限公司 前参照和后参照的尺寸测量导向器以及切削块和方法
CN107106307A (zh) * 2015-01-06 2017-08-29 沃尔德玛链接有限公司 确定适合患者的膝关节假体的股骨植入物尺寸的量具
CN107252338A (zh) * 2009-05-29 2017-10-17 史密夫和内修有限公司 用于执行膝关节成形术的方法和器械
CN110648337A (zh) * 2019-09-23 2020-01-03 武汉联影医疗科技有限公司 髋关节分割方法、装置、电子设备和存储介质
CN111179350A (zh) * 2020-02-13 2020-05-19 张逸凌 基于深度学习的髋关节图像处理方法及计算设备

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6917827B2 (en) * 2000-11-17 2005-07-12 Ge Medical Systems Global Technology Company, Llc Enhanced graphic features for computer assisted surgery system
EP4252697A3 (en) * 2017-03-14 2024-03-27 Stephen B. Murphy Systems and methods for determining leg length change during hip surgery

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101815477A (zh) * 2007-09-28 2010-08-25 株式会社力克赛 人工膝关节置换手术的术前计划装置及手术辅助用工具
CN107252338A (zh) * 2009-05-29 2017-10-17 史密夫和内修有限公司 用于执行膝关节成形术的方法和器械
CN103209652A (zh) * 2010-08-13 2013-07-17 史密夫和内修有限公司 外科手术引导件
CN106456196A (zh) * 2014-02-11 2017-02-22 史密夫和内修有限公司 前参照和后参照的尺寸测量导向器以及切削块和方法
CN107106307A (zh) * 2015-01-06 2017-08-29 沃尔德玛链接有限公司 确定适合患者的膝关节假体的股骨植入物尺寸的量具
CN110648337A (zh) * 2019-09-23 2020-01-03 武汉联影医疗科技有限公司 髋关节分割方法、装置、电子设备和存储介质
CN111179350A (zh) * 2020-02-13 2020-05-19 张逸凌 基于深度学习的髋关节图像处理方法及计算设备

Also Published As

Publication number Publication date
CN111888059A (zh) 2020-11-06
WO2022007972A1 (zh) 2022-01-13

Similar Documents

Publication Publication Date Title
CN111888059B (zh) 基于深度学习与x线的全髋关节图像处理方法及装置
US20240096508A1 (en) Systems and methods for using generic anatomy models in surgical planning
CN112971981B (zh) 基于深度学习的全髋关节图像处理方法和设备
US10991070B2 (en) Method of providing surgical guidance
CN101189638B (zh) 膝关节形态的特征描述方法和系统
EP1525560B1 (en) Automated measurement of objects using deformable models
CN114742747B (zh) 基于深度学习的髋关节置换术后影像的评估方法及系统
CN102132320A (zh) 图像处理,特别是医学图像处理的方法和设备
CN115252233B (zh) 基于深度学习的全髋关节置换术中髋臼杯的自动化规划方法
CN115456990A (zh) 一种基于ct图像的肋骨计数方法、装置、设备及存储介质
CN113077499B (zh) 骨盆配准方法、骨盆配准装置和骨盆配准系统
CN114708409A (zh) 手术中的定位导航方法、装置、存储介质及设备
CN114299177B (zh) 图像处理方法、装置、电子设备及存储介质
CN117152407A (zh) 一种头影测量标志点自动定位方法
US20070230782A1 (en) Method, a Computer Program, and Apparatus, an Image Analysis System and an Imaging System for an Object Mapping in a Multi-Dimensional Dataset
CN114372970A (zh) 一种手术参考信息生成方法及装置
Kotcheff et al. Shape model analysis of THR radiographs
Redhead et al. An automated method for assessing routine radiographs of patients with total hip replacements
CN117495693B (zh) 用于内窥镜的图像融合方法、系统、介质及电子设备
US20230240762A1 (en) Method and device for segmentation and registration of an anatomical structure
CN118717376B (zh) 基于镜像处理的髋臼杯假体信息的确定方法及装置
CN117745989B (zh) 基于椎管构造的神经根阻滞靶点注入路径规划方法及系统
CN113112560B (zh) 生理点区域标记方法与装置
CN114642444B (zh) 口腔种植精度评价方法、系统和终端设备
JP2024523863A (ja) X線画像計測を術中に使用する患者位置のための調整システム及び方法

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
CB03 Change of inventor or designer information

Inventor after: Zhang Yiling

Inventor after: Liu Xingyu

Inventor before: Zhang Yiling

Inventor before: Liu Xingyu

Inventor before: An Yicheng

Inventor before: Chen Peng

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 101102 room 402, 4th floor, building 28, yard 18, Kechuang 13th Street, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Patentee after: Beijing Changmugu Medical Technology Co.,Ltd.

Address before: 101102 room 402, 4th floor, building 28, yard 18, Kechuang 13th Street, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Patentee before: BEIJING CHANGMUGU MEDICAL TECHNOLOGY Co.,Ltd.

CP01 Change in the name or title of a patent holder