CN116433477B - 骨盆配准方法、装置、存储介质及电子设备 - Google Patents

骨盆配准方法、装置、存储介质及电子设备 Download PDF

Info

Publication number
CN116433477B
CN116433477B CN202310699566.7A CN202310699566A CN116433477B CN 116433477 B CN116433477 B CN 116433477B CN 202310699566 A CN202310699566 A CN 202310699566A CN 116433477 B CN116433477 B CN 116433477B
Authority
CN
China
Prior art keywords
line vector
connecting line
target object
included angle
determining
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
CN202310699566.7A
Other languages
English (en)
Other versions
CN116433477A (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 Yidian Lingdong Technology Co ltd
Original Assignee
Beijing Yidian Lingdong 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 Yidian Lingdong Technology Co ltd filed Critical Beijing Yidian Lingdong Technology Co ltd
Priority to CN202310699566.7A priority Critical patent/CN116433477B/zh
Publication of CN116433477A publication Critical patent/CN116433477A/zh
Application granted granted Critical
Publication of CN116433477B publication Critical patent/CN116433477B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/20Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
    • 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 or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor
    • A61F2/4603Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor for insertion or extraction of endoprosthetic joints or of accessories thereof
    • A61F2/4607Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor for insertion or extraction of endoprosthetic joints or of accessories thereof of hip femoral endoprostheses
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/101Computer-aided simulation of surgical operations
    • A61B2034/105Modelling of the patient, e.g. for ligaments or bones
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B34/00Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
    • A61B34/10Computer-aided planning, simulation or modelling of surgical operations
    • A61B2034/107Visualisation of planned trajectories or target regions
    • 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 or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor
    • A61F2002/4632Special tools or methods for implanting or extracting artificial joints, accessories, bone grafts or substitutes, or particular adaptations therefor using computer-controlled surgery, e.g. robotic surgery
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30008Bone

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Transplantation (AREA)
  • Theoretical Computer Science (AREA)
  • Molecular Biology (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Robotics (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Vascular Medicine (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Prostheses (AREA)

Abstract

本发明公开了一种骨盆配准方法、装置、存储介质及电子设备。其中,该方法包括:基于CT数据确定虚拟骨骼模型,双侧股骨头中心之间的第一连线向量,患侧股骨头中心与患侧髂前上棘之间的第二连线向量;将第一、二连线向量进行坐标转换,得到第三、四连线向量;在髋臼侧骨骼与虚拟骨骼模型之间配准通过的情况下,确定实际骨盆定位坐标系下双侧股骨头中心之间的第五连线向量,患侧股骨头中心与患侧髂前上棘之间的第六连线向量;计算第三、五连线向量,第四、六连线向量之间的夹角;根据夹角确定骨盆配准校验结果。本发明解决了由于髋关节手术中配准不准确、配准结果发生旋转,导致配准结果出现局部最优,假体植入角度与术前规划不一致的技术问题。

Description

骨盆配准方法、装置、存储介质及电子设备
技术领域
本发明涉及医疗技术领域,具体而言,涉及一种骨盆配准方法、装置、存储介质及电子设备。
背景技术
全髋关节置换导航中,骨盆配准注册是影响手术精准性非常重要的一个环节。将术中患者实际骨骼与对应的计算机中骨模型配准,以使实际的物理骨骼的位姿与虚拟骨模型相关联,骨头位姿包括骨头的位置与姿态。通过配准再术中可以根据术前对假体位置和深度的规划辅助髋关节置换过程中的磨削和植入。在传统手术中,患者髋臼暴露有限,不利于术中观察,通过导航系统可以从可视化界面了解当前磨削和植入情况,使得假体植入更精准,且假体植入角度更符合患者实际情况,降低假体磨损与假体脱位风险。而术中可视化精准的前提是真实骨骼与虚拟骨骼的配准精度,只有准确的配准才能使得术中跟踪准确。
传统骨骼点云配准方法中,由于算法本身特性和目标函数的约束,算法常易因噪声数据和不良点对信息影响而陷入局部优化解,丢失优化解。由于骨盆髋臼窝是类似半球形,若术中选点分布不均匀或者分布过于集中的情况下,会导致配准结果发生旋转,进而导致假体植入角度与术前规划不一致。进而容易使得术后假体发生脱位、磨损和活动度受限。
针对上述的问题,目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种骨盆配准方法、装置、存储介质及电子设备,以至少解决由于髋关节手术中配准不准确、配准结果发生旋转,导致配准结果出现局部最优,假体植入角度与术前规划不一致的技术问题。
根据本发明实施例的一个方面,提供了一种骨盆配准方法,包括:基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
根据本发明实施例的另一方面,还提供了一种骨盆配准装置,包括:第一确定模块,用于基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;坐标转换模块,用于将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;第二确定模块,用于在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;计算模块,用于计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;第三确定模块,用于根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
根据本发明实施例的另一方面,还提供了一种非易失性存储介质,上述非易失性存储介质存储有多条指令,上述指令适于由处理器加载并执行任意一项上述的骨盆配准方法。
根据本发明实施例的另一方面,还提供了一种电子设备,包括一个或多个处理器和存储器,上述存储器用于存储一个或多个程序,其中,当上述一个或多个程序被上述一个或多个处理器执行时,使得上述一个或多个处理器实现任意一项上述的骨盆配准方法。
在本发明实施例中,通过基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果,达到了在髋关节手术中对骨盆配准结果进行校验的目的,从而实现了提升患者髋臼侧骨骼和虚拟骨骼模型匹配准确性的技术效果,进而解决了由于髋关节手术中配准不准确、配准结果发生旋转,导致配准结果出现局部最优,假体植入角度与术前规划不一致的技术问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是根据本发明实施例的一种骨盆配准方法的示意图;
图2是根据本发明实施例的一种可选的骨盆配准方法的流程图;
图3是根据本发明实施例的一种可选的髋臼侧骨骼与虚拟骨骼模型之间的匹配结果示意图;
图4是根据本发明实施例的一种可选的骨盆配准方法的示意图;
图5是根据本发明实施例的另一种可选的骨盆配准方法的示意图;
图6是根据本发明实施例的一种骨盆配准装置的示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
根据本发明实施例,提供了一种骨盆配准的方法实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的骨盆配准方法的流程图,如图1所示,该方法包括如下步骤:
步骤S102,基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心。
可选的,上述目标对象可以为将要进行髋关节置换术的患者,上述双侧股骨头中心包括左侧股骨头中心和右侧股骨头中心,在患侧股骨头中心为左侧股骨头中心的情况下,对应的另一侧股骨头中心为右侧股骨头中心;同理若患侧股骨头中心为右侧股骨头中心,对应的另一侧股骨头中心为左侧股骨头中心。上述虚拟骨盆定位坐标系为基于CT数据获取到的目标对象的虚拟骨骼模型对应的骨盆定位坐标系。
在一种可选的实施例中,上述基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,包括:基于上述CT数据,确定上述目标对象的上述虚拟骨骼模型、上述双侧股骨头中心以及上述患侧髂前上棘;基于上述骨骼模型和预设髋关节假体位置,确定虚拟骨盆定位坐标系下上述双侧股骨头中心和上述患侧髂前上棘分别对应的初始定位坐标;根据上述双侧股骨头中心的初始定位坐标,确定上述第一连线向量;根据上述患侧股骨头中心和上述患侧髂前上棘分别对应的初始定位坐标,确定上述第二连线向量。
可选的,上述预设髋关节假体位置为术前规划好的髋关节假体位置,基于确定好目标对象的虚拟骨骼模型以及预设髋关节假体位置,确定双侧股骨头中心和患侧髂前上棘分别对应的初始定位坐标,并进一步计算双侧股骨头中心之间的第一连线向量和第二连线向量。通过虚拟骨骼模型以及预设髋关节假体位置在虚拟骨盆定位坐标系下的双重定位效应,达到提升双侧股骨头中心和患侧髂前上棘分别对应的初始定位坐标获取准确定,进而使得获取到的第一连线向量和第二连线向量更加准确可靠。
可选的,在术前规划阶段,基于患者CT数据获取目标对象的虚拟骨骼模型、双侧股骨头中心分别对应的初始定位坐标L_FHC_V、R_FHC_V,以及患侧髂前上棘对应的初始定位坐标ASP_V;术前策划读取与显示,软件读取术前策划,显示骨骼虚拟模型及策划髋关节假体位置,读取策划的双侧股骨头中心分别对应的初始定位坐标L_FHC_V、R_FHC_V,以及患侧髂前上棘对应的初始定位坐标ASP_V,计算双侧股骨头中心之间的第一连线向量α1,以及患侧股骨头中心和上述患侧髂前上棘之间的第二连线向量β1。
可选的,根据患者CT数据获取的双侧股骨头中心,即左侧股骨头中心和右侧股骨头中心分别在虚拟骨盆定位坐标系下的初始定位坐标L_FHC_V、R_FHC_V。
根据左侧股骨头中心和右侧股骨头中心分别在虚拟骨盆定位坐标系下的初始定位坐标L_FHC_V、R_FHC_V,确定双侧股骨头中心L_FHC_V和R_FHC_V之间的第一连线向量α1,若患侧为左侧,则对应的第一连线向量为:
若患侧为右侧,则对应的第一连线向量如下。
根据患者CT数据获取目标对象的患侧髂前上棘在虚拟骨盆定位坐标系下的初始定位坐标ASP_V,计算虚拟骨盆定位坐标系下患侧股骨头中心与患侧髂前上棘之间的第二连线向量β1,
若患侧为左侧,则对应的第二连线向量为:
若患侧为右侧,则对应的第二连线向量如下。
步骤S104,将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量。
可选的,上述实际骨盆定位坐标系是基于目标对象在手术中的实际骨盆位置确定的。通过以上方式,将虚拟骨盆定位坐标系下的第一连线向量、第二连线向量转换至实际骨盆定位坐标系,用于后续的骨盆配准校验。
步骤S106,在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量。
在一种可选的实施例中,上述在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述目标对象的上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量,包括:在上述目标对象的上述髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,获取上述实际骨盆定位坐标系下的第一点云集,以及上述患侧髂前上棘的当前定位坐标,其中,上述第一点云集是在绕髋关节中心旋转上述目标对象的股骨的情况下,通过安装于预定位置的定位工具采集到,上述患侧髂前上棘的当前定位坐标是通过探针捕获到的;基于上述第一点云集,采用最小二乘拟合法进行拟合处理,得到上述实际骨盆定位坐标系下上述双侧股骨头中心的当前定位坐标;根据上述双侧股骨头中心的当前定位坐标,确定上述第五连线向量;以及根据上述患侧股骨头中心和上述患侧髂前上棘分别对应的当前定位坐标,确定上述第六连线向量。
可选的,分别在目标对象的左右侧股骨远端固定定位工具,在目标对象处于仰卧位下,稳定目标对象骨盆,通过绕髋关节中心旋转目标对象的股骨,采集相应的点云,得到上述第一点云集。可以理解,在实际手术的过程中,双侧股骨头中心无法直接得到,通过以上方式,采集目标对象的股骨绕髋关节中心旋转得到的点云,构成包括有多个点云的第一点云集,基于第一点云集中包括的点云,采用最小二乘拟合法进行拟合处理,可以达到高效获取双侧股骨头中心的当前定位坐标,并且由于最小二乘拟合法是比较成熟的算法,采用该算法获取到的双侧股骨头中心的当前定位坐标具备一定的准确性和可靠性。
可选的,术中获取目标对象的双侧股骨头中心,即左侧股骨头中心和右侧股骨头中心分别在实际骨盆定位坐标系下的当前定位坐标L_FHC和R_FHC;基于左侧股骨头中心和右侧股骨头中心分别在实际骨盆定位坐标系下的当前定位坐标L_FHC和R_FHC,确定术中目标对象的双侧股骨头中心L_FHC和R_FHC之间的第五连线向量α2,其中,
若患侧为左侧,则对应的第五连线向量为:
若患侧为右侧,则对应的第五连线向量如下。
可选的,术中获取目标对象的患侧髂前上棘在实际骨盆定位坐标系下的当前定位坐标ASP,患侧股骨头中心和上述患侧髂前上棘分别对应的当前定位坐标,计算第六连线向量β2,其中,
若患侧为左侧,则对应的第六连线向量为:
若患侧为右侧,则对应的第六连线向量如下。
在一种可选的实施例中,在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述目标对象的上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量之前,上述方法还包括:确定上述目标对象的髋臼侧骨骼中第一数量的特征点;基于上述第一数量的特征点,对上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型进行第一点云匹配,得到初始匹配结果;确定上述目标对象的髋臼侧骨骼中第二数量的特征点,其中,上述第二数量大于上述第一数量;基于上述第二数量的特征点和上述初始匹配结果,对上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型进行第二点云匹配,得到目标匹配结果;在上述目标匹配结果指示上述髋臼侧骨骼与上述虚拟骨骼模型之间对应特征点的位置误差小于预设差值阈值的情况下,确定上述髋臼侧骨骼与上述虚拟骨骼模型之间配准通过。
可选的,上述第一数量的特征点可以但不限于包括在目标对象的髋臼窝前沿、后沿、上沿分别选择的三个特征点;上述第二数量的特征点可以但不限于包括:髋臼窝内及周围预定范围内均匀选择第二数量的特征点。例如,固定骨盆定位工具,在目标对象的髋臼窝前沿、后沿、上沿分别选择三个特征点进行粗配准,得到初始匹配结果。切换到精配准页面在髋臼窝内及周围均匀选择特征点进行精细配准,得到目标匹配结果;根据目标匹配结果计算髋臼侧骨骼与上述虚拟骨骼模型对应特征点之间的位置配准误差,若误差小于等于预设误差阈值(如1毫米),则进入下一步特征点选择,若误差大于1mm,则重新进行粗配准和精配准。通过以上方式,在进行目标对象的髋臼侧骨骼与虚拟骨骼模型之间配准时,首先选取少量关键点进行粗配准,在粗配准的基础上选择更多的点进行髋臼侧骨骼与虚拟骨骼模型的细配准,由此达到提升髋臼侧骨骼与虚拟骨骼模型配准准确性的目的。
可选的,可以但不限于采用迭代最近点算法(Iterative Closest Pointalgorithm,ICP)进行目标对象的髋臼侧骨骼和虚拟骨骼模型之间配准,图2是根据本发明实施例的一种可选的骨盆配准方法的流程图,如图2所示,具体包括:使用探针在目标对象的髋臼侧选择特征点获取输入点云集;对输入点云集进行预处理得到处理后的输入点云集P;根据处理后的输入点云集P和基于CT数据确定的术前策划得到的虚拟骨骼模型点云X进行迭代配准;构造转换矩阵将输入点云集P和虚拟骨骼模型点云X转换至同一坐标系(如实际骨盆定位坐标系)下,并通过目标函数判断髋臼侧骨骼与虚拟骨骼模型的配准结果是否收敛,其中,上述目标函数可以但不限于为阈值判断函数,即通过判断髋臼侧骨骼与虚拟骨骼模型对应特征点之间的位置配准误差是否小于预设差值阈值判断匹配结果是否收敛;如果匹配结果收敛则获取转换矩阵,并在此基础上获取实际骨盆定位坐标系下上述目标对象的上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;如果匹配结果不收敛则重新构造转换矩阵和目标函数判断,直至髋臼侧骨骼与虚拟骨骼模型的配准结果收敛。
步骤S108,计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角。
可选的,对目标对象的髋臼侧骨骼与虚拟骨骼模型之间配准结果进行验证,将术前策划获得的第一连线向量α1转到真实骨盆定位工具坐标系下,得到第三连线向量,利用第三连线向量/>和术中获得的第五连线向量α2,计算第三连线向量/>和第五连线向量α2之间的第一夹角θ,如图3所示,具体的计算方式如下:
将术前策划获得的第二连线向量β1转到实际骨盆定位工具坐标系下,得到第四连线向量,利用第四连线向量/>和术中获得的第六连线向量β2,计算第四连线向量/>和第六连线向量β2之间的第二夹角ω,如图3所示,具体的计算方式如下。
步骤S110,根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
可以理解,在进行髋关节置换术的过程中,需要进行患者髋臼侧骨骼与虚拟骨骼模型的配准,保证患者髋臼侧骨骼与虚拟骨骼模型的一致性。双侧股骨头中心之间的连线向量,以及患侧股骨头中心和患侧髂前上棘之间的连线向量能够很好的反映出目标对象的髋臼侧骨骼的相对位置。在目标对象的髋臼侧骨骼与虚拟骨骼模型之间配准通过的情况下,进一步通过计算双侧股骨头中心在虚拟骨盆定位坐标系和实际骨盆定位坐标系下的向量夹角,以及患侧股骨头中心和患侧髂前上棘在虚拟骨盆定位坐标系和实际骨盆定位坐标系下的向量夹角(即上述第一夹角和第二夹角),进行目标对象的骨盆配准校验,由此达到有效保证患者髋臼侧骨骼与虚拟骨骼模型的一致性的目的。
在一种可选的实施例中,上述根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果,包括:在确定上述第一夹角小于或等于预设第一夹角,并且上述第二夹角小于或等于预设第二夹角的情况下,确定上述目标对象的髋臼侧骨骼中第三数量的特征点;分别确定上述第三数量的特征点处,探针距离上述目标对象的骨面的距离,得到上述第三数量的特征点处分别对应的骨面距;在确定上述第三数量的特征点处分别对应的骨面距满足预设距离条件的情况下,确定上述骨盆配准校验结果为:上述目标对象的骨盆配准通过。
可选的,上述满足预设距离条件可以但不限于为:第三数量的特征点分别对应的骨面距的平均值小于预设距离阈值,或者每一个特征点对应的骨面距均小于预设距离阈值。
可选的,若第一夹角θ角度在可接受误差范围内,则验证第二夹角ω角度是否在可接受范围内,若第一夹角θ和第二夹角ω角度均在可接受范围内(如小于对应的预设夹角阈值)则进入距骨面验证。其中,距骨面验证用于通过确定探针与骨面之间的距离,验证患者髋臼侧骨骼与虚拟骨骼模型的配准情况,若距骨面可接受,则进入髋臼磨削和髋臼杯植入。通过以上方式,在进行匹配结果验证时,通过第一夹角、第二夹角、距骨面综合进行配准结果的判定,并且只有在第一夹角、第二夹角、距骨面均满足对应条件的情况下才确定目标对象的骨盆配准通过。由此达到从多个方面进行配准结果验证的目的,由此提升目标对象骨盆配准验证的准确性和可靠性。
在一种可选的实施例中,上述方法还包括:在确定上述第一夹角大于上述预设第一夹角的情况下,将上述第一夹角投影至多个预定平面,得到上述第一夹角对应的多个第一投影角,其中,上述多个预定平面至少包括:冠状面、矢状面、横断面;确定上述多个第一投影角中最大的第一投影角,以及上述最大的第一投影角对应的预定第一平面;确定上述预定第一平面对应的第一调整方向;按照上述第一调整方向对上述目标对象的上述髋臼侧骨骼进行调整,直至调整后的髋臼侧骨骼对应的新的第一夹角小于或等于上述预设第一夹角。
可选的,若第一夹角θ不在可接受范围内,则将第一夹角θ在实际骨盆坐标系下投影,分别得到第一夹角投影至多个预定平面得到的多个第一投影角,包括冠状面下的θ1投影角、矢状面下的θ2投影角、横断面下的θ3投影角,选择θ1、θ2、θ3中最大的角度在对应视角下进行修正,直到θ在可接受范围内。第一夹角在多个投影平面(即多个预定平面)的投影角度可以体现出在该方向上目标对象的髋臼侧骨骼与虚拟骨骼模型之间的偏差程度,将第一夹角分解至多个投影面进行逐个分析,夹角最大则表明在该方向上目标对象的髋臼侧骨骼与虚拟骨骼模型之间的偏差程度最大。通过以上方式可以快速确定出最大偏差方向,并基于最大偏差方向进行对应的夹角调整,由此达到降低目标对象的髋臼侧骨骼与虚拟骨骼模型之间的位置偏差,提升配准准确度的效果。
在一种可选的实施例中,上述方法还包括:在确定上述第二夹角大于上述预设第二夹角的情况下,将上述第二夹角投影至多个预定平面,得到上述第二夹角对应的多个第二投影角,其中,上述多个预定平面至少包括:冠状面、矢状面、横断面;确定上述多个第二投影角中最大的第二投影角,以及上述最大的第二投影角对应的第二预定平面;确定上述预定第二平面对应的第二调整方向,按照上述第二调整方向对上述目标对象的上述髋臼侧骨骼进行调整,直至调整后的髋臼侧骨骼对应的新的第二夹角小于或等于上述预设第二夹角。
可选的,若第二夹角ω不在可接受范围内,则将第二夹角ω在实际骨盆坐标系下投影,分别得到第二夹角投影至多个预定平面得到的多个第二投影角,包括冠状面下的ω1投影角、矢状面下的ω2投影角、横断面下的ω3投影角,选择ω1、ω2、ω3中最大的角度在对应视角下进行修正,直到ω在可接受范围内。第二夹角在多个投影平面(即多个预定平面)的投影角度可以体现出在该方向上目标对象的髋臼侧骨骼与虚拟骨骼模型之间的偏差程度,将第二夹角分解至多个投影面进行逐个分析,夹角最大则表明在该方向上目标对象的髋臼侧骨骼与虚拟骨骼模型之间的偏差程度最大。通过以上方式可以快速确定出最大偏差方向,并基于最大偏差方向进行对应的夹角调整,由此达到降低目标对象的髋臼侧骨骼与虚拟骨骼模型之间的位置偏差,提升配准准确度的效果。
本发明实施例提供的骨盆配准方法,应用于髋关节置换术中,通过上述步骤S102至步骤S110,可以达到在髋关节手术中对骨盆配准结果进行校验的目的,从而实现提升患者髋臼侧骨骼和虚拟骨骼模型匹配准确性的技术效果,进而解决由于髋关节手术中配准不准确、配准结果发生旋转,导致配准结果出现局部最优,假体植入角度与术前规划不一致的技术问题。
基于上述实施例和可选实施例,本发明提出一种可选实施方式,图4和5是根据本发明实施例的一种可选的骨盆配准方法的流程图,如图4和5所示,该方法包括:
生成术前规划,术前规划包含但不限于基于患者CT数据获取患者的虚拟骨骼模型、双侧股骨头中心分别对应的初始定位坐标L_FHC_V、R_FHC_V,以及患侧髂前上棘对应的初始定位坐标ASP_V。
术前策划读取与显示,软件读取术前策划,显示骨骼虚拟模型及策划髋关节假体位置,读取策划的双侧股骨头中心分别对应的初始定位坐标L_FHC_V、R_FHC_V,以及患侧髂前上棘对应的初始定位坐标ASP_V,计算双侧股骨头中心之间的第一连线向量α1,以及患侧股骨头中心和上述患侧髂前上棘之间的第二连线向量β1。
采用ICP配准方法进行患者的髋臼侧骨骼和虚拟骨骼模型之间配准,即使用探针在患者的髋臼侧选择特征点获取输入点云集P,根据输入点云集P和基于CT数据确定的策划输入骨骼点云X进行迭代配准获取最优配准结果。
在确定点云配准通过的情况下,进行特征点选择,分别在患者左右侧股骨远端固定定位工具,在仰卧位下,稳定患者骨盆,通过绕髋关节中心旋转患者股骨,利用最小二乘法拟合股骨头中心,判断所拟合股骨头中心误差是否可接受,获得双侧股骨头中心分别在实际骨盆定位坐标系下的当前定位坐标L_FHC和R_FHC,并计算双侧股骨头中心在实际骨盆定位坐标系下的第五连线向量α2。
特征点选择,使用探针捕获患侧髂前上棘在实际骨盆定位坐标系下的当前定位坐标ASP,计算患侧股骨头中心和上述患侧髂前上棘之间的第六连线向量β2。
对患者的髋臼侧骨骼与虚拟骨骼模型之间配准结果进行验证,将术前策划获得的第一连线向量α1转到真实骨盆定位工具坐标系(实际骨盆定位坐标系)下,得到第三连线向量,利用第三连线向量/>和术中获得的第五连线向量α2,计算第三连线向量/>和第五连线向量α2之间的第一夹角θ,如图3所示,具体的计算方式如下:
将术前策划获得的第二连线向量β1转到实际骨盆定位工具坐标系下,得到第四连线向量,利用第四连线向量/>和术中获得的第六连线向量β2,计算第四连线向量/>和第六连线向量β2之间的第二夹角ω,如图3所示,具体的计算方式如下:
若第一夹角θ角度在可接受误差范围内,则验证第二夹角角度是否在可接受范围 内,若第一夹角θ和第二夹角ω角度均在可接受范围内(如小于对应的预设夹角阈值)则进 入距骨面验证;若第一夹角θ不在可接受范围内,则将第一夹角θ在实际骨盆坐标系下投影, 分别得到冠状面下θ1、矢状面下θ2、横断面下θ3,选择θ1、θ2、θ3中最大的角度在对应视角下 进行修正,直到θ在可接受范围内;若ω不在可接受范围内,同理修正。
距骨面验证,使用探针在患者髋臼侧骨面选择特征点(扎透软骨),查看探针与骨面的距离,验证患者髋臼侧骨骼与虚拟骨骼模型的配准情况,若距骨面可接受,则进入髋臼磨削和髋臼杯植入,若距骨面不可接受,则重新对骨骼进行配准和特征点选择。
本发明实施例中,基于CT数据进行术前规划;本发明基于点云配准算法解决术中骨盆配准的旋转问题,实现更精准的可视化追踪;针对现有技术中存在的缺陷,本发明旨在提供一种具有更高精准度的软件系统、控制方法以及电子设备,以指导医生更准确地磨削和植入;针对传统手术方法存在的缺陷,本发明旨在实现更精准的配准校准方法及软件系统,实现假体植入的精准性,降低假体磨损和术后脱位率,给患者术后最大的活动空间。至少可以解决髋关节置换术中患者髋臼侧骨骼与虚拟骨骼模型配准不准确、配准结果发生旋转,从而导致配准结果出现局部最优,假体植入角度与术前规划不一致,进而使得术后假体发生脱位、磨损和活动度受限的问题。
基于上述实施例和可选实施例,本发明提出另一种可选实施方式,具体包括:
启动软件读取术前策划计算得到虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量α1,患侧股骨头中心与患侧髂前上棘之间的第二连线向量β1,固定骨盆定位工具,在患者的髋臼窝前沿、后沿、上沿分别选择三个特征点进行粗配准,得到初始匹配结果。切换到精配准页面在髋臼窝内及周围均匀选择特征点进行精细配准,得到目标匹配结果;根据目标匹配结果计算髋臼侧骨骼与上述虚拟骨骼模型对应特征点之间的位置配准误差,若误差小于等于预设误差阈值(如1mm),则进入下一步特征点选择,若误差大于1mm,则重新进行粗配准和精配准。
若精配准通过,则分别在患者左右侧股骨远端固定定位工具,在仰卧位下,稳定患者骨盆,绕髋关节中心旋转患者股骨,获得实际骨盆定位坐标系下的当前定位坐标L_FHC和R_FHC,并计算得到第五连线向量α2;使用探针捕获患侧髂前上棘在实际骨盆定位坐标系下的当前定位坐标ASP,计算第六连线向量β2;将虚拟骨盆定位坐标系下的第一连线向量和第二连线向量分别转换至实际骨盆定位坐标系下,得到对应的第三连线向量和第四连线向量;进一步基于第三连线向量和第五连线向量计算得到第一夹角θ,基于第四连线向量和第六连线向量计算得到第二夹角ω。
软件显示第一夹角θ和第二夹角ω的角度,判断第一夹角θ和第二夹角ω是否在可接受范围内,比如是否在±10°内,若不可接受,则点击修正按钮分别对第一夹角θ和第二夹角ω进行修正,软件界面显示修正后计算的第一夹角θ和第二夹角ω,继续判断第一夹角θ和第二夹角ω的值,直到认为角度可接受则进入距骨面验证。
使用探针在患者骨面采集特征点(扎透软骨),若距骨面显示均小于1mm,且平均误差小于1mm,则配准通过进入磨削和植杯,若距骨面不满足要求则重新进行配准及后续步骤。
根据页面导航提示进行磨削和臼杯植入。
需要说明的是,本发明实施例基于目标对象的髋臼侧骨骼对应的术前CT数据进行术前策划。本发明实施例不指导股骨侧置换,并且所实现的方法不限假体型号及品牌。
在本实施例中还提供了一种骨盆配准装置,该装置用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述。如以下所使用的,术语“模块”“装置”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
根据本发明实施例,还提供了一种用于实施上述骨盆配准方法的装置实施例,图6是根据本发明实施例的一种骨盆配准装置的结构示意图,如图6所示,上述骨盆配准装置,包括:第一确定模块600、坐标转换模块602、第二确定模块604、计算模块606、第三确定模块608,其中:
上述第一确定模块600,用于基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;
上述坐标转换模块602,连接于上述第一确定模块600,用于将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;
上述第二确定模块604,连接于上述坐标转换模块602,用于在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;
上述计算模块606,连接于上述第二确定模块604,用于计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;
上述第三确定模块608,连接于上述计算模块606,用于根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
在本发明实施例中,通过设置上述第一确定模块600,用于基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;上述坐标转换模块602,连接于上述第一确定模块600,用于将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;上述第二确定模块604,连接于上述坐标转换模块602,用于在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;上述计算模块606,连接于上述第二确定模块604,用于计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;上述第三确定模块608,连接于上述计算模块606,用于根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果,达到了在髋关节手术中对骨盆配准结果进行校验的目的,从而实现了提升患者髋臼侧骨骼和虚拟骨骼模型匹配准确性的技术效果,进而解决了由于髋关节手术中配准不准确、配准结果发生旋转,导致配准结果出现局部最优,假体植入角度与术前规划不一致的技术问题。
需要说明的是,上述各个模块是可以通过软件或硬件来实现的,例如,对于后者,可以通过以下方式实现:上述各个模块可以位于同一处理器中;或者,上述各个模块以任意组合的方式位于不同的处理器中。
此处需要说明的是,上述第一确定模块600、坐标转换模块602、第二确定模块604、计算模块606、第三确定模块608对应于实施例中的步骤S102至步骤S110,上述模块与对应的步骤所实现的实例和应用场景相同,但不限于上述实施例所公开的内容。需要说明的是,上述模块作为装置的一部分可以运行在计算机终端中。
需要说明的是,本实施例的可选或优选实施方式可以参见实施例中的相关描述,此处不再赘述。
上述的骨盆配准装置还可以包括处理器和存储器,上述第一确定模块600、坐标转换模块602、第二确定模块604、计算模块606、第三确定模块608等均作为程序模块存储在存储器中,由处理器执行存储在存储器中的上述程序模块来实现相应的功能。
处理器中包含内核,由内核去存储器中调取相应的程序模块,上述内核可以设置一个或以上。存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM),存储器包括至少一个存储芯片。
根据本发明实施例,还提供了一种非易失性存储介质的实施例。可选的,在本实施例中,上述非易失性存储介质包括存储的程序,其中,在上述程序运行时控制上述非易失性存储介质所在设备执行上述任意一种骨盆配准方法。
可选的,在本实施例中,上述非易失性存储介质可以位于计算机网络中计算机终端群中的任意一个计算机终端中,或者位于移动终端群中的任意一个移动终端中,上述非易失性存储介质包括存储的程序。
可选的,在程序运行时控制非易失性存储介质所在设备执行以下功能:基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
根据本发明实施例,还提供了一种处理器的实施例。可选的,在本实施例中,上述处理器用于运行程序,其中,上述程序运行时执行上述任意一种骨盆配准方法。
根据本发明实施例,还提供了一种计算机程序产品的实施例,当在数据处理设备上执行时,适于执行初始化有上述任意一种的骨盆配准方法步骤的程序。
可选的,上述计算机程序产品,当在数据处理设备上执行时,适于执行初始化有如下方法步骤的程序:基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
本发明实施例提供了一种电子设备,该电子设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序,处理器执行程序时实现以下步骤:基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定上述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,上述双侧股骨头中心包括上述患侧股骨头中心,以及除上述患侧股骨头中心之外的另一侧股骨头中心;将上述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将上述第二连线向量旋转至上述实际骨盆定位坐标系,得到第四连线向量;在上述目标对象的髋臼侧骨骼与上述虚拟骨骼模型之间配准通过的情况下,确定上述实际骨盆定位坐标系下上述双侧股骨头中心之间的第五连线向量,以及上述患侧股骨头中心与上述患侧髂前上棘之间的第六连线向量;计算上述第三连线向量和上述第五连线向量之间的第一夹角,以及上述第四连线向量和上述第六连线向量之间的第二夹角;根据上述第一夹角和上述第二夹角,确定上述目标对象的骨盆配准校验结果。
在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。
在本发明所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如上述模块的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,模块或模块的间接耦合或通信连接,可以是电性或其它的形式。
上述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能模块可以集成在一个处理模块中,也可以是各个模块单独物理存在,也可以两个或两个以上模块集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。
上述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取非易失性存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个非易失性存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的非易失性存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (9)

1.一种骨盆配准方法,其特征在于,包括:
基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定所述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,所述双侧股骨头中心包括所述患侧股骨头中心,以及除所述患侧股骨头中心之外的另一侧股骨头中心;
将所述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将所述第二连线向量旋转至所述实际骨盆定位坐标系,得到第四连线向量;
在所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型之间配准通过的情况下,确定所述实际骨盆定位坐标系下所述双侧股骨头中心之间的第五连线向量,以及所述患侧股骨头中心与所述患侧髂前上棘之间的第六连线向量;
计算所述第三连线向量和所述第五连线向量之间的第一夹角,以及所述第四连线向量和所述第六连线向量之间的第二夹角;
根据所述第一夹角和所述第二夹角,确定所述目标对象的骨盆配准校验结果;
其中,在所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型之间配准通过的情况下,确定所述实际骨盆定位坐标系下所述目标对象的所述双侧股骨头中心之间的第五连线向量,以及所述患侧股骨头中心与所述患侧髂前上棘之间的第六连线向量之前,所述方法还包括:确定所述目标对象的髋臼侧骨骼中第一数量的特征点;基于所述第一数量的特征点,对所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型进行第一点云匹配,得到初始匹配结果;确定所述目标对象的髋臼侧骨骼中第二数量的特征点,其中,所述第二数量大于所述第一数量;基于所述第二数量的特征点和所述初始匹配结果,对所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型进行第二点云匹配,得到目标匹配结果;在所述目标匹配结果指示所述髋臼侧骨骼与所述虚拟骨骼模型之间对应特征点的位置误差小于预设差值阈值的情况下,确定所述髋臼侧骨骼与所述虚拟骨骼模型之间配准通过。
2.根据权利要求1所述的方法,其特征在于,所述基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定所述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,包括:
基于所述CT数据,确定所述目标对象的所述骨骼模型、所述双侧股骨头中心以及所述患侧髂前上棘;
基于所述骨骼模型和预设髋关节假体位置,确定虚拟骨盆定位坐标系下所述双侧股骨头中心和所述患侧髂前上棘分别对应的初始定位坐标;
根据所述双侧股骨头中心的初始定位坐标,确定所述第一连线向量;
根据所述患侧股骨头中心和所述患侧髂前上棘分别对应的初始定位坐标,确定所述第二连线向量。
3.根据权利要求1所述的方法,其特征在于,所述在所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型之间配准通过的情况下,确定所述实际骨盆定位坐标系下所述目标对象的所述双侧股骨头中心之间的第五连线向量,以及所述患侧股骨头中心与所述患侧髂前上棘之间的第六连线向量,包括:
在所述目标对象的所述髋臼侧骨骼与所述虚拟骨骼模型之间配准通过的情况下,获取所述实际骨盆定位坐标系下的第一点云集,以及所述患侧髂前上棘的当前定位坐标,其中,所述第一点云集是在绕髋关节中心旋转所述目标对象的股骨的情况下,通过安装于预定位置的定位工具采集到,所述患侧髂前上棘的当前定位坐标是通过探针捕获到的;
基于所述第一点云集,采用最小二乘拟合法进行拟合处理,得到所述实际骨盆定位坐标系下所述双侧股骨头中心的当前定位坐标;
根据所述双侧股骨头中心的当前定位坐标,确定所述第五连线向量;以及根据所述患侧股骨头中心和所述患侧髂前上棘分别对应的当前定位坐标,确定所述第六连线向量。
4.根据权利要求1至3中任意一项所述的方法,其特征在于,所述根据所述第一夹角和所述第二夹角,确定所述目标对象的骨盆配准校验结果,包括:
在确定所述第一夹角小于或等于预设第一夹角,并且所述第二夹角小于或等于预设第二夹角的情况下,确定所述目标对象的髋臼侧骨骼中第三数量的特征点;
分别确定所述第三数量的特征点处,探针距离所述目标对象的骨面的距离,得到所述第三数量的特征点处分别对应的骨面距;
在确定所述第三数量的特征点处分别对应的骨面距满足预设距离条件的情况下,确定所述骨盆配准校验结果为:所述目标对象的骨盆配准通过。
5.根据权利要求4所述的方法,其特征在于,所述方法还包括:
在确定所述第一夹角大于所述预设第一夹角的情况下,将所述第一夹角投影至多个预定平面,得到所述第一夹角对应的多个第一投影角,其中,所述多个预定平面至少包括:冠状面、矢状面、横断面;
确定所述多个第一投影角中最大的第一投影角,以及所述最大的第一投影角对应的预定第一平面;
确定所述预定第一平面对应的第一调整方向;
按照所述第一调整方向对所述目标对象的所述髋臼侧骨骼进行调整,直至调整后的髋臼侧骨骼对应的新的第一夹角小于或等于所述预设第一夹角。
6.根据权利要求4所述的方法,其特征在于,所述方法还包括:
在确定所述第二夹角大于所述预设第二夹角的情况下,将所述第二夹角投影至多个预定平面,得到所述第二夹角对应的多个第二投影角,其中,所述多个预定平面至少包括:冠状面、矢状面、横断面;
确定所述多个第二投影角中最大的第二投影角,以及所述最大的第二投影角对应的预定第二平面;
确定所述预定第二平面对应的第二调整方向,按照所述第二调整方向对所述目标对象的所述髋臼侧骨骼进行调整,直至调整后的髋臼侧骨骼对应的新的第二夹角小于或等于所述预设第二夹角。
7.一种骨盆配准装置,其特征在于,包括:
第一确定模块,用于基于目标对象的髋臼侧骨骼对应的电子计算机断层扫描CT数据,确定所述目标对象的虚拟骨骼模型,虚拟骨盆定位坐标系下双侧股骨头中心之间的第一连线向量,以及患侧股骨头中心与患侧髂前上棘之间的第二连线向量,其中,所述双侧股骨头中心包括所述患侧股骨头中心,以及除所述患侧股骨头中心之外的另一侧股骨头中心;
坐标转换模块,用于将所述第一连线向量旋转至实际骨盆定位坐标系,得到第三连线向量;将所述第二连线向量旋转至所述实际骨盆定位坐标系,得到第四连线向量;
第二确定模块,用于在所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型之间配准通过的情况下,确定所述实际骨盆定位坐标系下所述双侧股骨头中心之间的第五连线向量,以及所述患侧股骨头中心与所述患侧髂前上棘之间的第六连线向量;
计算模块,用于计算所述第三连线向量和所述第五连线向量之间的第一夹角,以及所述第四连线向量和所述第六连线向量之间的第二夹角;
第三确定模块,用于根据所述第一夹角和所述第二夹角,确定所述目标对象的骨盆配准校验结果;
其中,所述装置还用于:确定所述目标对象的髋臼侧骨骼中第一数量的特征点;基于所述第一数量的特征点,对所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型进行第一点云匹配,得到初始匹配结果;确定所述目标对象的髋臼侧骨骼中第二数量的特征点,其中,所述第二数量大于所述第一数量;基于所述第二数量的特征点和所述初始匹配结果,对所述目标对象的髋臼侧骨骼与所述虚拟骨骼模型进行第二点云匹配,得到目标匹配结果;在所述目标匹配结果指示所述髋臼侧骨骼与所述虚拟骨骼模型之间对应特征点的位置误差小于预设差值阈值的情况下,确定所述髋臼侧骨骼与所述虚拟骨骼模型之间配准通过。
8.一种非易失性存储介质,其特征在于,所述非易失性存储介质存储有多条指令,所述指令适于由处理器加载并执行权利要求1至6中任意一项所述的骨盆配准方法。
9.一种电子设备,其特征在于,包括一个或多个处理器和存储器,所述存储器用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现权利要求1至6中任意一项所述的骨盆配准方法。
CN202310699566.7A 2023-06-13 2023-06-13 骨盆配准方法、装置、存储介质及电子设备 Active CN116433477B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310699566.7A CN116433477B (zh) 2023-06-13 2023-06-13 骨盆配准方法、装置、存储介质及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310699566.7A CN116433477B (zh) 2023-06-13 2023-06-13 骨盆配准方法、装置、存储介质及电子设备

Publications (2)

Publication Number Publication Date
CN116433477A CN116433477A (zh) 2023-07-14
CN116433477B true CN116433477B (zh) 2023-09-08

Family

ID=87085899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310699566.7A Active CN116433477B (zh) 2023-06-13 2023-06-13 骨盆配准方法、装置、存储介质及电子设备

Country Status (1)

Country Link
CN (1) CN116433477B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117152256B (zh) * 2023-10-30 2024-02-13 中国人民解放军总医院第一医学中心 一种基于模板的骨盆模型通道定位方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115005987A (zh) * 2022-07-01 2022-09-06 北京大学第三医院(北京大学第三临床医学院) 髋关节翻修手术中骨盆配准的方法及系统
CN116035619A (zh) * 2023-01-30 2023-05-02 复旦大学 一种基于超声骨骼成像的导航配准方法
CN116077177A (zh) * 2022-12-31 2023-05-09 杭州键嘉医疗科技股份有限公司 一种用于全髋置换手术的股骨柄假体自动推荐方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080125630A1 (en) * 2006-09-11 2008-05-29 Caylor Edward J System and method for determining a location of an orthopaedic medical device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115005987A (zh) * 2022-07-01 2022-09-06 北京大学第三医院(北京大学第三临床医学院) 髋关节翻修手术中骨盆配准的方法及系统
CN116077177A (zh) * 2022-12-31 2023-05-09 杭州键嘉医疗科技股份有限公司 一种用于全髋置换手术的股骨柄假体自动推荐方法
CN116035619A (zh) * 2023-01-30 2023-05-02 复旦大学 一种基于超声骨骼成像的导航配准方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陶星光 ; 汪方 ; 王秋根 ; 林艳萍 ; 王会祥.《骨盆三维导航术中A型超声骨面无创定位研究》.《国际骨科学杂志》.2014,第1-4页. *

Also Published As

Publication number Publication date
CN116433477A (zh) 2023-07-14

Similar Documents

Publication Publication Date Title
US10238403B2 (en) Patient-specific total hip arthroplasty
US20200345429A1 (en) Systems And Methods For Determining A Joint Center Of Rotation During A Procedure
US9408617B2 (en) Method for orienting an acetabular cup and instruments for use therewith
EP2667813B1 (en) Computer program for planning the positioning of an implant
US10993817B1 (en) Method for femur resection alignment approximation in hip replacement procedures
US8463004B2 (en) Determining shaft and femur neck axes and three-dimensional reconstruction
WO2012021858A2 (en) Implant alignment
CN111134840B (zh) 膝关节置换手术方案的生成装置和终端
CN116433477B (zh) 骨盆配准方法、装置、存储介质及电子设备
CN114259330B (zh) 髋臼杯假体的角度的测量方法、装置和测量系统
US9042621B2 (en) Intra-operative surgical plan changing
CN113077499B (zh) 骨盆配准方法、骨盆配准装置和骨盆配准系统
CN113729943B (zh) 虚拟骨面的处理方法、装置及翻修手术机器人
CN115005987A (zh) 髋关节翻修手术中骨盆配准的方法及系统
US9889021B2 (en) Method for hip resection alignment approximation in hip replacement procedures
CN116439833B (zh) 骨盆配准处理方法、装置、存储介质及电子设备
CN116348064A (zh) 用于分析髋臼杯位置的系统和方法
McDonald et al. The effect of anatomic landmark selection of the distal humerus on registration accuracy in computer-assisted elbow surgery
EP4014911B1 (en) Artificial-intelligence-based detection of invisible anatomical structures in 2d x-ray images
CN116784978B (zh) 手术导航系统
US20240099775A1 (en) Artificial-intelligence-based determination of implantation curve
EP4014912A1 (en) Artificial-intelligence-based registration of x-ray images
CN114569296A (zh) 一种髋关节手术导板
CN115300102A (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
GR01 Patent grant
GR01 Patent grant