CN116363184B - 一种髋关节位置配准方法、装置、电子设备及存储介质 - Google Patents
一种髋关节位置配准方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN116363184B CN116363184B CN202310615797.5A CN202310615797A CN116363184B CN 116363184 B CN116363184 B CN 116363184B CN 202310615797 A CN202310615797 A CN 202310615797A CN 116363184 B CN116363184 B CN 116363184B
- Authority
- CN
- China
- Prior art keywords
- registration
- dimensional model
- hip joint
- area
- transformation
- 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
Links
Classifications
-
- 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/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/344—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving models
-
- 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/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
- G06T7/75—Determining position or orientation of objects or cameras using feature-based methods involving models
-
- 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)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Image Processing (AREA)
Abstract
本发明实施例提供了一种髋关节位置配准方法、装置、电子设备及存储介质,涉及图像处理技术领域。髋关节位置配准方法包括:确定三维模型上的第一配准区域;实时获得第一候选位置;根据第一变换关系变换第一候选位置,得到第二候选位置;判断第二候选位置是否位于第一配准区域内;当第二候选位置位于第一配准区域内时,响应于位置获得指令,获得工具的当前位置,作为第一配准位置;根据三维模型表面各个点的位置、第一配准位置与第一变换关系,确定目标髋关节上的点与三维模型表面的点间位置的第二变换关系;根据第二变换关系,对目标髋关节上的实际位置进行位置配准。应用本发明实施例提供的方案可以准确地对髋关节的实际位置进行配准。
Description
技术领域
本发明涉及图像处理技术领域,特别是涉及一种髋关节位置配准方法、装置、电子设备及存储介质。
背景技术
髋关节置换手术是一种常见的人工关节置换术,用于重建患者髋关节的正常功能。目前,医生需要根据患者的髋关节CT(Computed Tomography,电子计算机断层扫描)图像规划手术路径,然后手术机器人根据医生规划好的手术路径与对髋关节的实际位置的配准结果,控制机械臂移动,从而协助医生进行髋关节置换手术。
现有技术中,需要基于髋关节CT图像以及实体髋关节获得进行配准的配准点,然后根据所获得的配准点的位置对髋关节进行位置配准。这样得到的配准结果准确度低。
发明内容
本发明实施例的目的在于提供一种髋关节位置配准方法、装置、电子设备及存储介质,以实现准确地对患者髋关节的实际位置进行配准。具体技术方案如下:
根据本发明实施例的第一方面,提供了一种髋关节位置配准方法,所述方法包括:
确定三维模型上的第一配准区域,其中,所述三维模型为目标髋关节对应的三维模型;
实时获得用于在所述目标髋关节上采集配准位置的工具的第一候选位置;
根据第一变换关系,对所述第一候选位置进行位置变换,得到第二候选位置,其中,第一变换关系为预先获得的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系;
判断所述第二候选位置是否位于所述第一配准区域内;
当所述第二候选位置位于所述第一配准区域内时,响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置;
根据所述三维模型表面各个点的位置、第一配准位置与所述第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系;
根据所述第二变换关系,对所述目标髋关节上的实际位置进行变换,完成位置配准。
可选的,在所述第二候选位置位于所述第一配准区域内时,所述方法还包括:
获得所述第二候选位置与所述三维模型表面间的最短距离;
判断所述最短距离是否小于预设误差距离;
若为是,则执行所述响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置的步骤。
可选的,所述方法还包括:
获得采集配准位置的工具在所述目标髋关节上对应于所述第一配准区域的区域内采集的第二区域更新位置;
根据所述第二区域更新位置与所述第一变换关系,确定所述三维模型表面上与所述第二区域更新位置对应的第一区域更新位置;
根据所确定的第一区域更新位置,更新所述第一配准区域。
可选的,所述根据所述三维模型表面各个点的位置、已获得的第一配准位置与第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系,包括:
根据所述第一变换关系对各第一配准位置进行位置变换,得到第二配准位置;
针对每一第二配准位置,根据该第二配准位置与所述三维模型表面各个点的位置,确定所述三维模型表面上与该第二配准位置距离最近的第三配准位置;
根据各第一配准位置以及与其对应的第三配准位置,确定所述目标髋关节与所述三维模型间的第二变换关系。
可选的,所述确定所述三维模型上的第一配准区域,包括:
获得用户在所述三维模型上选择的多个髋关节的表征性位置,并获得各表征性位置的描述信息;根据所获得的表征性位置与描述信息,确定第一配准区域;
和/或
获得用户在所述三维模型上框选的区域;将所述框选的区域确定为第一配准区域。
可选的,所述获得所述工具的当前位置,作为第一配准位置,包括:
获得与双目相机相连的带有反光阵列的手术探针的当前位置,作为第一配准位置。
根据本发明实施例的第二方面,提供了一种髋关节位置配准装置,所述装置包括:
第一配准区域确定模块,用于确定所述三维模型上的第一配准区域;
第一候选位置获得模块,用于实时获得用于在所述目标髋关节上采集配准位置的工具的第一候选位置;
第二候选位置变换模块,用于根据第一变换关系,对所述第一候选位置进行位置变换,得到第二候选位置,其中,第一变换关系为预先获得的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系;
位置判断模块,用于判断所述第二候选位置是否位于所述第一配准区域内;
第一配准位置确定模块,用于当所述第二候选位置位于所述第一配准区域内时,响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置;
第二变换关系确定模块,用于根据所述三维模型表面各个点的位置、第一配准位置与所述第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系;
目标髋关节位置配准模块,用于根据所述第二变换关系,对所述目标髋关节上的实际位置进行变换,完成位置配准。
可选的,所述装置还包括:
最短距离获得模块,用于获得所述第二候选位置与所述三维模型表面间的最短距离;
最短距离判断模块,用于判断所述最短距离是否小于预设误差距离,若为是,则触发第一配准位置确定模块。
可选的,所述装置还包括:
第二区域更新位置获得模块,用于获得采集配准位置的工具在所述目标髋关节上对应于所述第一配准区域的区域内采集的第二区域更新位置;
第一区域更新位置确定模块,用于根据所述第二区域更新位置与所述第一变换关系,确定所述三维模型表面上与所述第二区域更新位置对应的第一区域更新位置;
第一配准区域更新模块,用于根据所确定的第一区域更新位置,更新所述第一配准区域。
可选的,所述第二变换关系确定模块,包括:
第二配准位置获得子模块,用于根据所述第一变换关系对各第一配准位置进行位置变换,得到第二配准位置;
第三配准位置确定子模块,用于针对每一第二配准位置,根据该第二配准位置与所述三维模型表面各个点的位置,确定所述三维模型表面上与该第二配准位置距离最近的第三配准位置;
第二变换关系确定子模块,用于根据各第一配准位置以及与其对应的第三配准位置,确定所述目标髋关节与所述三维模型间的第二变换关系。
可选的,所述第一配准区域确定模块,具体用于获得用户在所述三维模型上选择的多个髋关节的表征性位置,并获得各表征性位置的描述信息;根据所获得的表征性位置与描述信息,确定第一配准区域;和/或 获得用户在所述三维模型上框选的区域;将所述框选的区域确定为第一配准区域。
可选的,所述第一配准位置确定模块,具体用于获得与双目相机相连的带有反光阵列的手术探针的当前位置,作为第一配准位置。
根据本发明实施例的第三方面,提供了一种电子设备,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的计算机程序时,实现上述第一方面所述的髋关节位置配准方法。
根据本发明实施例的第四方面,提供了一种计算机可读存储介质,所述计算机可读存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现上述第一方面所述的髋关节位置配准方法。
本发明实施例有益效果:
本发明实施例提供的髋关节位置配准方案,可以确定三维模型上的第一配准区域,并获得采集配准位置的工具在目标髋关节上采集的第一候选位置。对第一候选位置进行位置变换,得到第二候选位置。判断第二候选位置是否位于第一配准区域内,将位于第一配准区域内的第二候选位置对应的第一候选位置确定为第一配准位置,这样,可以引导使用者采集经过位置变换后的位置在第一配准区域内的位置作为第一配准位置,保证了获得的第一配准位置的准确度。由于获得的第一配准位置的准确度高,根据第一配准位置确定的第二变换关系的准确度高,使用第二变换关系对目标髋关节进行位置配准的准确度高。因此,应用本发明实施例提供的方案可以准确地对髋关节的实际位置进行配准。
当然,实施本发明的任一产品或方法并不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其他的实施例。
图1为本发明实施例提供的第一种髋关节位置配准方法的流程示意图;
图2a为本发明实施例提供的三维模型上对应于目标髋关节上髋臼窝内侧的第一配准位置示意图;
图2b为本发明实施例提供的三维模型上对应于目标髋关节上髋臼窝外侧的第一配准位置示意图;
图3为本发明实施例提供的一种更新第一配准区域的流程示意图;
图4为本发明实施例提供的第二种髋关节位置配准方法的流程示意图;
图5为本发明实施例提供的髋关节位置配准的整体流程示意图;
图6为本发明实施例提供的一种髋关节位置配准装置的结构示意图;
图7为本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员基于本申请所获得的所有其他实施例,都属于本发明保护的范围。
髋关节置换手术是一种常见的人工关节置换术。目前,可以使用手术机器人辅助医生进行髋关节置换手术。手术机器人可以根据医生规划好的路径控制机械臂移动,从而帮助医生完成手术。
由于手术过程中髋关节的位置和形态与采集髋关节CT图像时髋关节的位置和形态往往不同,手术机器人难以精准地按照医生规划的手术路径控制机械臂移动。鉴于此,本发明实施例提供了一种髋关节位置配准方法,实现准确地对髋关节的实际位置进行配准。
下面对本发明实施例的执行主体进行说明。
本发明实施例的执行主体可以是电子设备,具体的,可以为与手术机器人相连接的电子设备,例如台式机等。
下面通过具体实施例对本发明实施例提供的髋关节位置配准方法进行说明。
参见图1,提供了第一种髋关节位置配准方法的流程示意图,该方法包括以下步骤S101-步骤S107。
步骤S101:确定三维模型上的第一配准区域。
其中,三维模型为对目标髋关节对应的三维模型,目标髋关节为患者术侧髋关节。
下面对第一配准区域进行说明。
第一配准区域为三维模型表面的区域,可以包括三维模型上髋臼窝内侧区域与髋臼窝外侧区域两部分区域。
如图2a与图2b中的黑色区域所示,第一配准区域为三维模型上的对应于目标髋关节上髋臼窝内侧与外侧的区域。其中,图2a提供了三维模型上对应于目标髋关节上髋臼窝内侧的第一配准位置示意图,图2b提供了三维模型上对应于目标髋关节上髋臼窝外侧的第一配准位置示意图。
下面对电子设备确定第一配准区域的方式进行说明。
第一种实现方式中,电子设备可以获得用户在三维模型上选择的多个髋关节的表征性位置,并获得各表征性位置的描述信息;根据所获得的表征性位置与描述信息,确定第一配准区域。
其中,髋关节的表征性位置是使用者根据医学知识确定的用于描述髋关节结构的位置,各表征性位置的信息用于表示各表征性位置所描述的髋关节结构。
具体的,电子设备可以在通信连接的显示设备上显示三维模型,这样使用者可以对三维模型进行观察,并使用输入设备向电子设备输入其选择的目标髋关节的表征性位置。
一种情况下,电子设备可以预设各表征性位置的获得顺序以及各表征性位置的描述信息,这样,使用者可以按预设的顺序向电子设备输入其选择的髋关节的表征性位置,电子设备根据获得顺序,确定各表征性位置对应的描述信息。
另一种情况下,使用者可以在输入其选择的目标髋关节的表征性位置时,也输入各表征性位置的描述信息。这样,电子设备可以获得各表征性位置与各表征性位置的描述信息。
第二种实现方式中,电子设备可以获得用户在三维模型上选择的框选区域;将框选区域确定为第一配准区域。
具体的,电子设备可以在通信连接的显示设备上显示三维模型,这样使用者可以对三维模型进行观察,使用输入设备向电子设备输入其选择的区域,电子设备将输入设备输入的区域确定为第一配准区域。具体的,上述输入设备可以是鼠标。触控板等等。使用者可以使用鼠标、触控板等输入设备以框选的方式确定其所选择的区域。
一种情况下,使用者可以选择三维模型上对应目标髋关节的实际暴露区域的区域作为框选区域。
这样,电子设备可以根据髋关节在医学上的结构特点,高效而准确地确定第一配准区域;也可以根据目标髋关节的实际暴露情况,确定能够实际用于配准的第一配准区域。通过上述方式,可以提高第一配准区域确定的效率和准确性。
第三种实现方式中,电子设备内可以预设髋关节上第一配准区域内的位置对应的形状参数,这样,电子设备可以分析计算三维模型表面各个位置的形状参数,并确定形状参数符合预设的第一配准区域内的位置对应的形状参数的位置属于第一配准区域,从而确定第一配准区域。例如,预设的形状参数可以为髋臼窝的内侧区域的曲率向量与外侧区域的曲率向量,电子设备可以通过分析计算各个位置的曲率向量,确定属于髋臼窝内侧区域的位置与属于髋臼窝外侧区域的位置,从而确定包括髋臼窝内侧区域与髋臼窝外侧区域在内的第一配准区域。
步骤S102:实时获得用于在目标髋关节上采集配准位置的工具的第一候选位置。
其中,上述第一候选位置为采集配准位置的工具的位置。一种情况下,可以为电子设备获得的采集配准位置的工具的位置的光学坐标。
步骤S103:根据第一变换关系,对第一候选位置进行位置变换,得到第二候选位置。
其中,第一变换关系为预先获得的目标髋关节上的点与三维模型表面的点间位置的变换关系。第一种实现方式中,第一变换关系可以为根据以往配准经验确定的目标髋关节上的点与三维模型表面的点间位置的变换关系。第二种实现方式中,第一变换关系可以为根据在三维模型上规划的点的位置,以及这些点在目标髋关节上对应点的位置预先确定的变换关系。
其中,第二候选位置为第一候选位置变换后得到的位置,为三维模型上的位置。当第一候选位置为空间坐标系下的坐标时,第二候选位置为模型坐标系下的坐标。第一候选位置与第二候选位置一一对应。
一种情况下,第二候选位置为电子设备获得的采集配准位置的工具的位置的光学坐标经位置变换后得到的模型坐标。
步骤S104:判断第二候选位置是否位于第一配准区域内。
具体的,可以通过以下方式判断第二候选位置是否位于第一配准区域内。
一种实现方式中,可以判断第二候选位置是否为第一配准区域包含的点的位置。若是,则第二候选位置位于第一配准区域内;若不是,则第二候选位置位于第一配准区域内。
另一种实现方式中,由于采集配准位置的工具可以为手术探针,此时可以根据第二候选位置以及采集配准位置的工具的几何形状,确定采集配准位置的工具,即手术探针在三维模型上对应的直线,并判断该直线是否与第一配准区域相交。若是,则第二候选位置位于第一配准区域内;若不是,则第二候选位置位于第一配准区域内。
步骤S105:当第二候选位置位于第一配准区域内时,响应于位置获得指令,获得工具的当前位置,作为第一配准位置。
上述工具的当前位置即为对应于位置获得指令的工具的第一候选位置,为目标髋关节上的位置。一种情况下,可以为电子设备获得的采集配准位置的工具的位置的光学坐标。上述位置获得指令可以为使用者通过工具向电子设备发动的指令。例如,可以在工具上设置采集按钮,当使用者按下采集按钮时,向电子设备发送位置获得指令;又例如,可以设置采集踏板,当使用者踩下踏板时,向电子设备发送位置获得指令。
所获得的工具的当前位置为目标髋关节上对应于第一配准区域的的位置,可以将该区域称为第二配准区域。由于第一配准区域包括三维模型上髋臼窝内侧区域与髋臼窝外侧区域两部分区域,因此第二配准区域也包括目标髋关节上髋臼窝内侧区域与髋臼窝外侧区域两部分区域。
第一配准位置为在第二配准区域内采集的目标髋关节上的点的位置,并且这些点位于第二配准区域内。由于目标髋关节存在于实际空间中,第一配准位置为这些点在实际空间中的位置。具体的,第一配准位置可以为这些点实际空间对应的空间坐标系下的坐标。
由于存在多个第一配准位置,因此,可以给各第一配准位置设置序号。序号可以为从1开始的连续正整数。如1、2、3……等。
由于第二配准区域也包括目标髋关节上髋臼窝内侧区域与髋臼窝外侧区域两部分区域,第一配准位置可以包括在髋臼窝内侧与髋臼窝外侧采集的点的位置。
第一预设数量可以为1、3、5等,本发明实施例对此并不进行限定。
下面对电子设备获得第一配准位置的方式进行说明。
第一种实现方式中,电子设备可以获得与双目相机相连的带有反光阵列的手术探针的当前位置,作为第一配准位置。具体的,可以由使用者将手术探针针尖移动至第二配准区域内的位置,并按下采集按钮,此时双目相机将针尖所在的位置发送至电子设备,电子设备将接收到的位置信息确定为第一配准位置。
此时,上述空间坐标系可以为双目相机采集的图像对应的光学坐标系,第一配准位置可以为在第二配准区域内采集的目标髋关节上的点在光学坐标系下的坐标。
这样,双目相机通过反光阵列的特定结构确定手术探针针尖所在的位置,可以准确地获得第一配准位置。
第二种实现方式中,电子设备可以获得与带有深度传感器的相机相连的手术探针采集的位置,作为第一配准位置。具体的,使用者将手术探针针尖移动至第二配准区域内的位置,并按下采集按钮,相机会将确定的位置发送至电子设备,电子设备将接收到的位置信息依次确定为第一配准位置。
第三种实现方式中,电子设备可以将相机采集得到的目标髋关节的图像与三维模型进行内容匹配,根据内容匹配结果,结合目标髋关节的图像,在目标髋关节上第二配准区域内确定第一配准位置。
步骤S106:根据三维模型表面各个点的位置、第一配准位置与第一变换关系,确定新的目标髋关节上的点与三维模型表面的点间位置的变换关系,作为第二变换关系。
三维模型表面的各个点的位置可以为三维模型表面的各个点在模型坐标系下的模型坐标。例如,各个点的CT坐标。
三维模型与目标髋关节对应,可以认为三维模型上的点与目标髋关节上的点之间只发生了刚体变换。因此,所确定的第二变换关系可以通过以下述形式表示。
一种情况下,所确定的第二变换关系可以以表达式进行表示,由于三维模型上的点与目标髋关节上的点之间发生了刚体变换,因此,表达式参数包括平移矩阵与旋转矩阵。假设目标髋关节上的点的坐标为P,该点对应的三维模型上的点的坐标为Q,R为第二变换关系对应的第一旋转矩阵,T为第二变换关系对应的第一平移矩阵,那么,使用表达式表示第二变换关系时,该表达式为:
;
另一种情况下,所确定的第二变换关系也可以由平移矩阵与旋转矩阵确定的变换矩阵表示。
同样假设目标髋关节上的点的坐标为P,该点对应的三维模型上的点的坐标为Q,R为第二变换关系对应的第一旋转矩阵,T为第二变换关系对应的第一平移矩阵,使用变换矩阵表示第二变换关系时,该变换矩阵为:。
类似的,第一变换关系也可以由表达式或变换矩阵的形式表示。
下面对电子设备确定第二变换关系的实现方式进行说明。
第一种实现方式中,可以根据第一变换关系,从三维模型表面各个点中确定对应于各第一配准位置的位置,从而根据各第一配准位置与三维模型表面对应的位置确定第二变换关系。具体的实现方式将在后文进一步说明,此处暂不详述。
第二种实现方式中,可以每采集第一预设数量个第一配准位置,就根据三维模型表面各个点的位置、已获得的第一配准位置与目标变换关系,确定第二变换关系。具体的实现方式与目标变换关系的获得方式将在后文进一步说明,此处暂不详述。
第三种实现方式中,可以根据第一变换关系,从三维模型表面各个点中确定与第一配准位置相对应的位置的点的位置,并对第一配准位置与第一配准位置对应的点的位置使用NDT(Normal Distributions Transform,正态分布变换)算法,确定第二变换关系。
步骤S107:根据第二变换关系,对目标髋关节上的实际位置进行变换,完成位置配准。
假设目标髋关节上的实际位置为,配准后对应在三维模型上的位置为/>,如前文对步骤S106的说明,当第二变换关系的形式为表达式时,有
;
当第二变换关系的形式为变换矩阵时,记该变换矩阵为,从而有
。
这样,电子设备即可根据第二变换关系,完成目标髋关节的位置配准。
上述目标髋关节的实际位置可以为根据后续具体应用中涉及的目标髋关节上的位置。例如,在完成位置配准后,需要根据配准结果,引导手术机器人的手术路径,此时,上述目标髋关节的实际位置可以为手术机器人的手术路径涉及到的目标髋关节上的位置。
一种情况下,电子设备可以将目标髋关节上的各实际位置进行变换,例如,电子设备可以接收双目相机采集的目标髋关节上的各实际位置的光学坐标,并变换每一光学坐标,得到对应的模型坐标。另一种情况下,电子设备可以变换目标髋关节上的部分实际位置的坐标,例如,电子设备可以获得采集位置的工具所采集的目标髋关节上的位置的光学坐标,并变换得到的光学坐标,得到对应的模型坐标。本发明实施例对此并不进行限定。
应用本实施例提供的方案,可以获得采集配准位置的工具在目标髋关节上采集的第一候选位置,并对第一候选位置进行位置变换,得到第二候选位置。通过判断第二候选位置是否位于第一配准区域内,将位于第一配准区域内的第二候选位置对应的第一候选位置确定为第一配准位置,可以引导使用者采集经过位置变换后的位置在第一配准区域内的位置作为第一配准位置,保证了获得的第一配准位置的准确度,这样,根据第一配准位置确定的第二变换关系的准确度高,使用第二变换关系对对目标髋关节进行位置配准的准确度高。因此,应用本发明实施例提供的方案可以准确地对髋关节的实际位置进行配准。
此外,由于获得的第一配准位置准确度高,因此,使用第一配准位置进行配准时,达到同样的配准精度需要的第一配准位置的数量少,提高了配准效率。
进一步的,本发明的一个实施例中,当第二候选位置位于第一配准区域内的情况下,在步骤S105前,还包括以下步骤A-步骤B。
步骤A:获得第二候选位置与三维模型表面间的最短距离。
上述第二候选位置与三维模型表面间的最短距离可以为第二候选位置到三维模型表面的最近距离,该距离是三维模型表面法线方向的距离。具体的,电子设备可以确定三维模型表面上法线经过第二候选位置的点,从而确定该点与第二候选位置之间的距离,作为第二候选位置与三维模型表面间的最短距离。
步骤B:判断最短距离是否小于预设误差距离。
若最短距离小于预设误差距离,说明第二候选位置的准确度能够满足配准需要,第二候选位置对应的第一候选位置的准确度能够满足配准的需要。此时,执行步骤S105。
另外,电子设备还可以提示使用者第一候选位置的准确度能够满足配准的需要。例如,在与其通信连接的显示设备上显示满足配准需要的信号提示使用者第一候选位置的准确度能够满足配准的需要;又例如,发出提示声音信号提示使用者第一候选位置的准确度能够满足配准的需要。
上述预设误差距离可以为5mm、10mm等,本发明实施例对此并不进行限定。
这样,可以在确保第一配准位置经过位置变换后的位置在第一配准区域内的前提下,进一步保证第一配准位置经过位置变换后的位置与三维模型表面间的最短距离小于预设误差距离,也就保证了获得的第一配准位置的准确度。
由于获得的第一配准位置准确度高,因此,使用第一配准位置进行配准时,达到同样的配准精度需要的第一配准位置的数量少;使用同样数量的第一配准位置进行配准时,达到满足要求的配准精度的成功率高。
将第一候选位置作为第一配准位置,所获得的第一配准位置的准确度能够满足配准的需要。
当第二候选位置与三维模型表面间的最短距离不小于预设误差距离时,说明第二候选位置的准确度不能满足配准需要,第二候选位置对应的第一候选位置的准确度也不能满足配准的需要。此时,不能将第一候选位置作为第一配准位置。
如前文对步骤S105的说明,电子设备可以获得与带有深度传感器的相机相连的手术探针采集的位置,作为第一配准位置。具体的,使用者将将手术探针针尖移动至第二配准区域内的位置,并按下采集按钮,此时双目相机将针尖所在的位置发送至电子设备,电子设备将接收到的位置信息确定为第一配准位置。
这样,可以提高获得的第一配准位置的准确度,根据第一配准位置确定第二变换关系的准确度高。基于这样确定的第二变换关系进行位置配准,可以提高配准的准确度。
本发明的一个实施例中,可以根据目标髋关节的实际情况,更新第一配准区域。具体的,参见图3,提供了一种更新第一配准区域的流程示意图,可以通过以下步骤S301-步骤S303完成。
步骤S301:获得采集配准位置的工具在目标髋关节上对应于第一配准区域的区域内采集的第二区域更新位置。
上述第二区域更新位置为采集配准位置的工具的位置。一种情况下,可以为电子设备获得的采集配准位置的工具的位置的光学坐标。
步骤S302:根据目标髋关节上采集的第二区域更新位置与第一变换关系,确定三维模型表面上与目标髋关节上采集的第二区域更新位置对应的第一区域更新位置。
其中,第一区域更新位置为第二区域更新位置变换后得到的位置,为三维模型上的位置。当第二区域更新位置为空间坐标系下的坐标时,第一区域更新位置为模型坐标系下的坐标。第一区域更新位置与第二区域更新位置一一对应。
步骤S303:根据所确定的第一区域更新位置,更新第一配准区域。
具体的,电子设备可以将所确定的第一区域更新位置作为更新后的第一配准区域的边界位置,从而依据边界位置确定更新后的第一配准区域,完成对第一配准区域的更新。
由于第一配准区域与第二配准区域相对应,而第二配准区域在目标髋关节上,可以认为,步骤S301-步骤S303是根据更新后的第二配准区域更新第一配准区域。步骤S301获得的第二区域更新位置是更新后的第二配准区域的边界位置,因此,第二区域更新位置经位置变换后,得到的第一区域更新位置,是更新后的第一配准区域的边界位置。这样,电子设备可以根据第一区域更新位置更新第一配准区域。
上述步骤S301-步骤S303可以在步骤S105,即获得多个第一配准位置之前进行,也可以在获得第一配准位置中进行,还可以在获得多个第一配准位置后进行,本发明对此并不进行限定。
这样,可以根据目标髋关节的实际情况,更新第一配准区域,从而保证第一配准区域内的点可以被采集获得,也就保证了获得的第一配准位置的准确度。
如前文对步骤S106的说明,本发明的一个实施例中,可以根据各第一配准位置与三维模型表面对应的位置确定第二变换关系。具体的,参见图4,提供了第二种髋关节位置配准方法的流程示意图。本实施例中,上述步骤S106可以通过以下步骤S106A-步骤S106B实现。
步骤S106A:根据第一变换关系对各第一配准位置进行位置变换,得到第二配准位置。
其中,第二配准位置为第一配准位置变换后得到的位置,为三维模型上的位置。当第一配准位置为空间坐标系下的坐标时,第二配准位置为模型坐标系下的坐标。第一配准位置与第二配准位置一一对应。
一种情况下,上述第二配准位置为第一配准位置的光学坐标经位置变换后得到的模型坐标,例如,第一配准位置的光学坐标经位置变换后得到的CT坐标。
一种实现方式中,假设第一变换关系为变换矩阵matrixC,第一配准位置为,变换得到的第二配准位置为/>,那么,电子设备可以通过以下关系式对第一配准位置进行位置变换:
。
步骤S106B:针对每一第二配准位置,根据该第二配准位置与三维模型表面各个点的位置,确定三维模型表面上与该第二配准位置距离最近的第三配准位置。
如对步骤S106A的说明,第二配准位置可以为模型坐标系下的坐标。而三维模型表面各个点的位置也可以为模型坐标系下的坐标。这样,可以通过坐标之间的关系确定三维模型表面上与该第二配准位置距离最近的第三配准位置。第二配准位置与第三配准位置一一对应。
上述第三配准位置为距离第一配准位置的光学坐标经位置变换后得到的模型坐标最近的点的模型坐标,可以认为每一第三配准位置对应的点是该第三配准位置对应的第二配准位置对应的点的对应近点。
下面对电子设备确定第三配准位置的实现方式进行说明。
一种实现方式中,针对每一第二配准位置,可以计算模型表面各个点的位置与该第二配准位置之间的距离,并从计算得到距离中确定与该第二配准位置最近的距离,将该距离对应的位置确定为第三配准位置。
另一种实现方式中,可以使用kd-tree(k-dimensional tree,k维树)方法确定第三配准位置。具体的,可以将三维模型所在的空间分割成多个子空间,每个子空间中只存在一个提取到的三维模型表面的点。这样,电子设备可以根据第二配准位置所对应的子空间,确定该子空间中存在的三维模型表面的点为与该第二配准位置距离最近的点,该三维模型表面的点的位置为第三配准位置。
步骤S106C:根据各第一配准位置以及与其对应的第三配准位置,确定目标髋关节与三维模型间的第二变换关系。
其中,由于第一配准位置与第二配准位置一一对应,第二配准位置又与第三配准位置一一对应,因此第一配准位置与第三配准位置一一对应。
下面对电子设备确定第二变换关系的实现方式进行说明。
一种实现方式中,电子设备可以对第一配准位置与第三配准位置使用SVD(Singular Value Decomposition, 奇异值分解)法确定第二变换关系。
具体的,电子设备可以根据以下关系式对第一配准位置与第三配准位置使用SVD法进行位置配准,求取平移矩阵与旋转矩阵。
;
;
;
;
;
;
;
;
其中,n表示确定的第三配准位置的数量,由于第一配准位置与第三配准位置一一对应,n也表示第一配准位置的数量;p表示第一配准位置对应的第一质心位置,表示第三配准位置对应的第二质心位置;/>表示第i个第一配准位置,/>表示第i个第三配准位置;表示第i个第一配准位置相对于第一质心位置的第一质心位移向量,/>表示第i个第三配准位置相对于第二质心位置的第二质心位移向量;H为构造的中间矩阵,/>表示对中间矩阵进行奇异值分解的过程,U为第一分解矩阵,V为第二分解矩阵;R表示解得的旋转矩阵,T表示解得的平移矩阵。
这样,电子设备可以基于上述解得的旋转矩阵与平移矩阵,确定第二变换关系。
另一种实现方式中,电子设备可以基于预设的刚体变换约束条件与预设变换调整结束条件,对各第一配准位置以及与其对应的第三配准位置进行迭代计算,从而确定第二变换关系。具体的实现方式将在后文实施例中进一步说明,此处暂不详述。
应用上述方式,可以使用第一变换关系对第一配准位置进行变换,得到第二配准位置,并确定三维模型表面上与每一第二配准位置距离最近的第三配准位置,从而根据第一配准位置与对应的第三配准位置确定第二变换关系。这样确定的第二变换关系是在第一变换关系的基础上进一步确定的,因此,确定的第二变换关系比第一变换关系更准确,基于第二变换关系进行位置配准,可以提高配准的准确度。
本发明的一个实施例中,如前文对步骤S106C的说明,可以基于预设的刚体变换约束条件与预设变换调整结束条件,对各第一配准位置以及与其对应的第三配准位置进行迭代计算,从而确定第二变换关系。具体的,本实施例中,上述步骤S106C可以通过以下步骤C-步骤E完成。
步骤C:基于每一配准位置对中两位置间的距离,确定满足预设的刚体变换约束条件的旋转参数和平移参数,并获得基于旋转参数和平移参数确定的目标髋关节与三维模型间的第二变换关系。
其中,每一配准位置对中包括:一个第一配准位置和与该第一配准位置对应的第二配准位置距离最近的第三配准位置。
具体的,上述旋转参数可以以旋转矩阵进行表示,上述平移参数可以以平移矩阵进行表示。这样,与前文对步骤S104的说明相同,可以根据旋转矩阵与平移矩阵确定第二变换关系。
上述预设的刚体变换约束条件可以有多种情况。
一种情况下,预设的刚体变换约束条件可以为:使得对各第一配准位置进行刚体变换后的位置与各第一配准位置对应的第三配准位置间平均距离最小的约束条件。
这种情况下,通过上述预设的刚体变换约束条件确定第二变换关系,可以使各第一配准位置进行刚体变换后的位置与各第一配准位置对应的第三配准位置间平均距离最小,也就使第二变换关系精确,基于第二变换关系进行位置配准,可以提高配准的准确度。
具体的,可以根据以下表达式确定旋转参数和平移参数。
;/>
其中,表示预设的刚体变换约束条件,/>表示第一配准位置的数量,由于第一配准位置与第三配准位置一一对应,/>也表示第三配准位置的数量;/>表示确定的旋转参数,/>表示确定的平移参数;/>表示第j个第三配准位置,/>表示第j个第一配准位置。
另一种情况下,预设的刚体变换约束条件可以为:使得对各第一配准位置进行刚体变换后的位置与各第一配准位置对应的第三配准位置间最大距离最小的约束条件。
步骤D:采用平移参数和旋转参数对各第一配准位置进行刚体变换,得到各第一配准位置对应的新的第二配准位置。
此处得到新的第二配准位置的方式与前文对步骤S106A的说明相类似,因此不再详细说明。
步骤E:针对每一第一配准位置,计算该第一配准位置对应的新的第二配准位置与该第一配准位置对应的第三配准位置的距离。
当距离不满足预设变换调整结束条件时,返回步骤S106B。
对应于对步骤C的说明中预设的刚体变换约束条件的有多种情况,上述预设变换调整结束条件也存在多种情况。
一种情况下,预设变换调整结束条件为:各第一配准位置对应的新的第二配准位置与该第一配准位置对应的第三配准位置的平均距离小于预设的距离阈值。
这种情况下,通过上述预设变换调整结束条件调整第二变换关系,可以使各第一配准位置进行刚体变换后的位置与各第一配准位置对应的第三配准位置间平均距离小于预设的距离阈值,也就使第二变换关系精确,基于第二变换关系进行位置配准,可以提高配准的准确度。
另一种情况下,预设变换调整结束条件为:各第一配准位置对应的新的第二配准位置与该第一配准位置对应的第三配准位置的最大距离小于预设的距离阈值。
通过以上方式,基于预设的刚体变换约束条件与预设变换调整结束条件确定并调整第二变换关系,可以使第二变换关系的准确度不断提高,从而得到的最终确定的第二变换关系具有高准确度。这样,基于第二变换关系进行位置配准,可以提高配准的准确度。
本发明的一个实施例中,如前文对步骤S106的说明,可以根据三维模型表面各个点的位置、已获得的第三配准位置与目标变换关系,确定第二变换关系。具体的,步骤S106可以通过以下步骤S106D实现。
步骤S106D:当获得第一预设数量个第一配准位置时,根据三维模型表面各个点的位置、已获得的第一配准位置与目标变换关系,确定目标髋关节上的点与三维模型表面的点间位置的第二变换关系,若已获得的第一配准位置的数量未达到第二预设数量,将目标变换关系更新为第二变换关系,并返回执行步骤S105。
其中,第一预设数量可以为1、3、5等,本发明实施例对此并不进行限定。
目标变换关系的初始关系为第一变换关系。第二预设数量可以为15、20等,本发明实施例对此并不进行限定。
电子设备可以在通信连接的显示设备上显示已获得的第一配准位置数量,也可以显示还需要获得的第一配准位置的数量,即第二预设数量与已获得的第一配准位置数量的差值。当已获得的第一配准位置的数量达到第二预设数量时,电子设备可以在通信连接的显示设备上显示提醒结束采集信息,提醒使用者结束采集。
此外,一种实现方式中,在每次更新目标变换关系后,电子设备可以计算使用目标变换关系对目标髋关节进行配准时存在的误差,并根据计算所得的误差,调整第一预设数量。具体的,可以设置上述误差与第一预设数量的对应关系。例如,当计算所得误差为5mm时,第一预设数量为6;当计算所得误差为4mm时,第一预设数量为3。
下面对第二变换关系的确定方式进行说明。
一种情况下,可以根据目标变换关系,从三维模型表面各个点中确定对应于各第一配准位置的位置,从而根据各第一配准位置与三维模型表面对应的位置确定第二变换关系。
另一种情况下,可以根据目标变换关系,从三维模型表面各个点中确定与第一配准位置相对应的位置的点的位置,并对第一配准位置与第一配准位置对应的点的位置使用NDT算法,确定第二变换关系。
通过上述方式,可以在获得第一预设数量个第一配准位置后,根据初始关系为第一变换关系的目标变换关系更新第二变换关系,并将第二变换关系确定为目标变换关系。这样,随着已获得的第一配准位置数量增加,第二变换关系会不断迭代更新,从而提高确定的第二变换关系的准确度。基于第二变换关系进行位置配准,可以提高配准的准确度。
下面对髋关节位置配准的整体流程进行说明。
参见图5,图5提供了髋关节位置配准的整体流程示意图。具体的,包括以下步骤S501-步骤S512。
步骤S501:获取目标髋关节的三维模型。
步骤S502:确定三维模型上的第一配准区域。
步骤S503:使用者在目标髋关节上移动采集配准位置的工具。
步骤S504:根据第一变换关系,实时获得采集配准位置的工具在三维模型中的位置。
步骤S505:判断采集配准位置的工具在三维模型中的对应位置是否在第一配准区域内,若是,执行步骤S506;若否,执行步骤S511。
步骤S506:判断采集配准位置的工具在三维模型中的对应位置与三维模型表面的最短距离是否小于预设误差距离。
若是,执行步骤S507;若否,执行步骤S512。
步骤S507:通过声光信息,提醒使用者在目标髋关节上对应于第一配准区域的第二配准区域内采集第一配准位置,获得采集得到的第一配准位置。
步骤S508:判断在第二配准区域内获得的第一配准位置数量是否不小于预设数量。
若是,执行步骤S509;若否,执行步骤S503。
步骤S509;根据第一变换关系与第一配准位置,确定三维模型与目标髋关节之间的第二变换关系。
步骤S510:根据第二变换关系,对目标髋关节上的实际位置进行变换。
步骤S511:移动采集配准位置的工具,使采集配准位置的工具在三维模型中的对应位置位于第一配准区域内,并执行步骤S504。
步骤S512:移动采集配准位置的工具接近目标髋关节,并执行步骤S504。
一种情况下,上述步骤S511可以通过以下步骤S511A-步骤S511D完成。
步骤S511A:判断采集配准位置的工具能否移动到第一配准区域对应的第二配准区域。若是,执行步骤S511D;若否,执行步骤S511B。
步骤S511B:使用采集配准位置的工具获取第二区域更新位置。
步骤S511C:根据第一变换关系与第二区域更新位置,确定第一区域更新位置,并基于第一区域更新位置更新第一配准区域。
步骤S511D:移动采集配准位置的工具,使采集配准位置的工具在三维模型中的对应位置位于第一配准区域内,并执行步骤S504。
与上述髋关节位置配准方法相对应,本发明实施例还提供了一种髋关节位置配准装置。
参见图6,提供了一种髋关节位置配准装置的结构示意图,所述装置包括:
第一配准区域确定模块601,用于确定所述三维模型上的第一配准区域;
第一候选位置获得模块602,用于实时获得用于在所述目标髋关节上采集配准位置的工具的第一候选位置;
第二候选位置变换模块603,用于根据第一变换关系,对所述第一候选位置进行位置变换,得到第二候选位置,其中,第一变换关系为预先获得的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系;
位置判断模块604,用于判断所述第二候选位置是否位于所述第一配准区域内;
第一配准位置确定模块605,用于当所述第二候选位置位于所述第一配准区域内时,响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置;
第二变换关系确定模块606,用于根据所述三维模型表面各个点的位置、第一配准位置与所述第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系;
目标髋关节位置配准模块607,用于根据所述第二变换关系,对所述目标髋关节上的实际位置进行变换,完成位置配准。
应用本实施例提供的方案,可以获得采集配准位置的工具在目标髋关节上采集的第一候选位置,并对第一候选位置进行位置变换,得到第二候选位置。通过判断第二候选位置是否位于第一配准区域内,将位于第一配准区域内的第二候选位置对应的第一候选位置确定为第一配准位置,可以引导使用者采集经过位置变换后的位置在第一配准区域内的位置作为第一配准位置,保证了获得的第一配准位置的准确度,这样,根据第一配准位置对目标髋关节进行位置配准的准确度高。因此,应用本发明实施例提供的方案可以准确地对髋关节的实际位置进行配准。
本发明的一个实施例中,所述装置还包括:
最短距离获得模块,用于获得所述第二候选位置与所述三维模型表面间的最短距离;
最短距离判断模块,用于判断所述最短距离是否小于预设误差距离,若为是,则触发第一配准位置确定模块。
这样,可以提高获得的第一配准位置的准确度,根据第一配准位置确定第二变换关系的准确度高。基于这样确定的第二变换关系进行位置配准,可以提高配准的准确度。
本发明的一个实施例中,所述装置还包括:
第二区域更新位置获得模块,用于获得采集配准位置的工具在所述目标髋关节上对应于所述第一配准区域的区域内采集的第二区域更新位置;
第一区域更新位置确定模块,用于根据所述第二区域更新位置与所述第一变换关系,确定所述三维模型表面上与所述第二区域更新位置对应的第一区域更新位置;
第一配准区域更新模块,用于根据所确定的第一区域更新位置,更新所述第一配准区域。
这样,可以根据目标髋关节的实际情况,更新第一配准区域,从而保证第一配准区域内的点可以被采集获得,也就保证了获得的第一配准位置的准确度。
本发明的一个实施例中,所述第二变换关系确定模块606,包括:
第二配准位置获得子模块,用于根据所述第一变换关系对各第一配准位置进行位置变换,得到第二配准位置;
第三配准位置确定子模块,用于针对每一第二配准位置,根据该第二配准位置与所述三维模型表面各个点的位置,确定所述三维模型表面上与该第二配准位置距离最近的第三配准位置;
第二变换关系确定子模块,用于根据各第一配准位置以及与其对应的第三配准位置,确定所述目标髋关节与所述三维模型间的第二变换关系。
应用上述方式,可以使用第一变换关系对第一配准位置进行变换,得到第二配准位置,并确定三维模型表面上与每一第二配准位置距离最近的第三配准位置,从而根据第一配准位置与对应的第三配准位置确定第二变换关系。这样确定的第二变换关系是在第一变换关系的基础上进一步确定的,因此,确定的第二变换关系比第一变换关系更准确,基于第二变换关系进行位置配准,可以提高配准的准确度。
本发明的一个实施例中,所述第一配准区域确定模块601,具体用于获得用户在所述三维模型上选择的多个髋关节的表征性位置,并获得各表征性位置的描述信息;根据所获得的表征性位置与描述信息,确定第一配准区域;和/或 获得用户在所述三维模型上框选的区域;将所述框选的区域确定为第一配准区域。
这样,电子设备可以根据髋关节在医学上的结构特点,高效而准确地确定第一配准区域;也可以根据目标髋关节的实际暴露情况,确定能够实际用于配准的第一配准区域。通过上述方式,可以提高第一配准区域确定的效率和准确性。
本发明的一个实施例中,所述第一配准位置确定模块605,具体用于获得与双目相机相连的带有反光阵列的手术探针的当前位置,作为第一配准位置。
这样,双目相机通过反光阵列的特定结构确定手术探针针尖所在的位置,可以准确地获得第一配准位置。
本发明实施例还提供了一种电子设备,如图7所示,包括处理器701、通信接口702、存储器703和通信总线704,其中,处理器701,通信接口702,存储器703通过通信总线704完成相互间的通信,
存储器703,用于存放计算机程序;
处理器701,用于执行存储器703上所存放的程序时,实现前述方法实施例所述的髋关节位置配准方法。
上述电子设备提到的通信总线可以是外设部件互连标准(PeripheralComponentInterconnect,PCI)总线或扩展工业标准结构(Extended Industry StandardArchitecture,EISA)总线等。该通信总线可以分为地址总线、数据总线、控制总线等。为便于表示,图中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。
通信接口用于上述电子设备与其他设备之间的通信。
存储器可以包括随机存取存储器(Random Access Memory,RAM),也可以包括非易失性存储器(Non-Volatile Memory,NVM),例如至少一个磁盘存储器。可选的,存储器还可以是至少一个位于远离前述处理器的存储装置。
上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(Network Processor,NP)等;还可以是数字信号处理器(DigitalSignalProcessor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
在本发明提供的又一实施例中,还提供了一种计算机可读存储介质,该计算机可读存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现前述方法实施例所述的髋关节位置配准方法。
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘Solid State Disk (SSD))等。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置、电子设备及存储介质实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的较佳实施例,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (12)
1.一种髋关节位置配准方法,其特征在于,所述方法包括:
确定三维模型上的第一配准区域,其中,所述三维模型为目标髋关节对应的三维模型;
实时获得用于在所述目标髋关节上采集配准位置的工具的第一候选位置;
根据第一变换关系,对所述第一候选位置进行位置变换,得到第二候选位置,其中,第一变换关系为预先获得的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系;
判断所述第二候选位置是否位于所述第一配准区域内;
当所述第二候选位置位于所述第一配准区域内时,响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置;
根据所述三维模型表面各个点的位置、第一配准位置与所述第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系;
根据所述第二变换关系,对所述目标髋关节上的实际位置进行变换,完成位置配准;
所述根据所述三维模型表面各个点的位置、第一配准位置与第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系,包括:
根据所述第一变换关系对各第一配准位置进行位置变换,得到第二配准位置;
针对每一第二配准位置,根据该第二配准位置与所述三维模型表面各个点的位置,确定所述三维模型表面上与该第二配准位置距离最近的第三配准位置;
根据各第一配准位置以及与其对应的第三配准位置,确定所述目标髋关节与所述三维模型间的第二变换关系。
2.根据权利要求1所述的方法,其特征在于,在所述第二候选位置位于所述第一配准区域内时,所述方法还包括:
获得所述第二候选位置与所述三维模型表面间的最短距离;
判断所述最短距离是否小于预设误差距离;
若为是,则执行所述响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置的步骤。
3.根据权利要求1所述的方法,其特征在于,所述方法还包括:
获得采集配准位置的工具在所述目标髋关节上对应于所述第一配准区域的区域内采集的第二区域更新位置;
根据所述第二区域更新位置与所述第一变换关系,确定所述三维模型表面上与所述第二区域更新位置对应的第一区域更新位置;
根据所确定的第一区域更新位置,更新所述第一配准区域。
4.根据权利要求1-3任一项所述的方法,其特征在于,所述确定所述三维模型上的第一配准区域,包括:
获得用户在所述三维模型上选择的多个髋关节的表征性位置,并获得各表征性位置的描述信息;根据所获得的表征性位置与描述信息,确定第一配准区域;
和/或
获得用户在所述三维模型上框选的区域;将所述框选的区域确定为第一配准区域。
5.根据权利要求1-3任一项所述的方法,其特征在于,所述获得所述工具的当前位置,作为第一配准位置,包括:
获得与双目相机相连的带有反光阵列的手术探针的当前位置,作为第一配准位置。
6.一种髋关节位置配准装置,其特征在于,所述装置包括:
第一配准区域确定模块,用于确定三维模型上的第一配准区域;
第一候选位置获得模块,用于实时获得用于在目标髋关节上采集配准位置的工具的第一候选位置;
第二候选位置变换模块,用于根据第一变换关系,对所述第一候选位置进行位置变换,得到第二候选位置,其中,第一变换关系为预先获得的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系;
位置判断模块,用于判断所述第二候选位置是否位于所述第一配准区域内;
第一配准位置确定模块,用于当所述第二候选位置位于所述第一配准区域内时,响应于位置获得指令,获得所述工具的当前位置,作为第一配准位置;
第二变换关系确定模块,用于根据所述三维模型表面各个点的位置、第一配准位置与所述第一变换关系,确定新的所述目标髋关节上的点与所述三维模型表面的点间位置的变换关系,作为第二变换关系;
目标髋关节位置配准模块,用于根据所述第二变换关系,对所述目标髋关节上的实际位置进行变换,完成位置配准;
所述第二变换关系确定模块,包括:
第二配准位置获得子模块,用于根据所述第一变换关系对各第一配准位置进行位置变换,得到第二配准位置;
第三配准位置确定子模块,用于针对每一第二配准位置,根据该第二配准位置与所述三维模型表面各个点的位置,确定所述三维模型表面上与该第二配准位置距离最近的第三配准位置;
第二变换关系确定子模块,用于根据各第一配准位置以及与其对应的第三配准位置,确定所述目标髋关节与所述三维模型间的第二变换关系。
7.根据权利要求6所述的装置,其特征在于,所述装置还包括:
最短距离获得模块,用于获得所述第二候选位置与所述三维模型表面间的最短距离;
最短距离判断模块,用于判断所述最短距离是否小于预设误差距离,若为是,则触发第一配准位置确定模块。
8.根据权利要求6所述的装置,其特征在于,所述装置还包括:
第二区域更新位置获得模块,用于获得采集配准位置的工具在所述目标髋关节上对应于所述第一配准区域的区域内采集的第二区域更新位置;
第一区域更新位置确定模块,用于根据所述第二区域更新位置与所述第一变换关系,确定所述三维模型表面上与所述第二区域更新位置对应的第一区域更新位置;
第一配准区域更新模块,用于根据所确定的第一区域更新位置,更新所述第一配准区域。
9.根据权利要求6-8任一项所述的装置,其特征在于,所述第一配准区域确定模块,具体用于获得用户在所述三维模型上选择的多个髋关节的表征性位置,并获得各表征性位置的描述信息;根据所获得的表征性位置与描述信息,确定第一配准区域;和/或 获得用户在所述三维模型上框选的区域;将所述框选的区域确定为第一配准区域。
10.根据权利要求6-8任一项所述的装置,其特征在于,所述第一配准位置确定模块,具体用于获得与双目相机相连的带有反光阵列的手术探针的当前位置,作为第一配准位置。
11.一种电子设备,其特征在于,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的计算机程序时,实现权利要求1-5任一所述的方法步骤。
12.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质内存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1-5任一所述的方法步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310615797.5A CN116363184B (zh) | 2023-05-26 | 2023-05-26 | 一种髋关节位置配准方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310615797.5A CN116363184B (zh) | 2023-05-26 | 2023-05-26 | 一种髋关节位置配准方法、装置、电子设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116363184A CN116363184A (zh) | 2023-06-30 |
CN116363184B true CN116363184B (zh) | 2023-10-03 |
Family
ID=86922466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310615797.5A Active CN116363184B (zh) | 2023-05-26 | 2023-05-26 | 一种髋关节位置配准方法、装置、电子设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116363184B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109890281A (zh) * | 2016-08-30 | 2019-06-14 | 马科外科公司 | 用于术中骨盆配准的***和方法 |
CN111476832A (zh) * | 2020-03-23 | 2020-07-31 | 杭州柳叶刀机器人有限公司 | 一种髋臼杯配准方法 |
CN113570648A (zh) * | 2021-07-30 | 2021-10-29 | 武汉联影智融医疗科技有限公司 | 多骨骼影像配准方法、电子装置以及医学导航*** |
CN113558766A (zh) * | 2021-07-19 | 2021-10-29 | 北京纳通医学研究院有限公司 | 图像配准方法、装置、手术机器人以及手术机器人*** |
CN114463393A (zh) * | 2021-12-31 | 2022-05-10 | 浙江大华技术股份有限公司 | 图像配准方法、计算机设备及存储装置 |
CN114938995A (zh) * | 2022-05-27 | 2022-08-26 | 骨圣元化机器人(深圳)有限公司 | 应用于髋关节置换手术的骨盆配准***和医疗设备 |
CN115005987A (zh) * | 2022-07-01 | 2022-09-06 | 北京大学第三医院(北京大学第三临床医学院) | 髋关节翻修手术中骨盆配准的方法及*** |
CN115100258A (zh) * | 2022-08-29 | 2022-09-23 | 杭州三坛医疗科技有限公司 | 一种髋关节图像配准方法、装置、设备以及存储介质 |
CN115153835A (zh) * | 2022-07-22 | 2022-10-11 | 华南理工大学 | 基于特征点配准与增强现实的髋臼假体放置引导***及方法 |
CN115607279A (zh) * | 2022-09-26 | 2023-01-17 | 苏州微创畅行机器人有限公司 | 硬组织的配准方法、装置及*** |
CN115624385A (zh) * | 2022-09-19 | 2023-01-20 | 重庆生物智能制造研究院 | 术前空间配准方法及装置、计算机设备、存储介质 |
CN115661214A (zh) * | 2022-09-16 | 2023-01-31 | 杭州三坛医疗科技有限公司 | 一种配准精度验证方法及装置 |
CN115969515A (zh) * | 2022-12-31 | 2023-04-18 | 杭州键嘉医疗科技股份有限公司 | 基于快速配准的机械臂辅助全髋关节置换手术 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8538108B2 (en) * | 2005-12-20 | 2013-09-17 | University Of Maryland, Baltimore | Method and apparatus for accelerated elastic registration of multiple scans of internal properties of a body |
-
2023
- 2023-05-26 CN CN202310615797.5A patent/CN116363184B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109890281A (zh) * | 2016-08-30 | 2019-06-14 | 马科外科公司 | 用于术中骨盆配准的***和方法 |
CN111476832A (zh) * | 2020-03-23 | 2020-07-31 | 杭州柳叶刀机器人有限公司 | 一种髋臼杯配准方法 |
CN113558766A (zh) * | 2021-07-19 | 2021-10-29 | 北京纳通医学研究院有限公司 | 图像配准方法、装置、手术机器人以及手术机器人*** |
CN113570648A (zh) * | 2021-07-30 | 2021-10-29 | 武汉联影智融医疗科技有限公司 | 多骨骼影像配准方法、电子装置以及医学导航*** |
CN114463393A (zh) * | 2021-12-31 | 2022-05-10 | 浙江大华技术股份有限公司 | 图像配准方法、计算机设备及存储装置 |
CN114938995A (zh) * | 2022-05-27 | 2022-08-26 | 骨圣元化机器人(深圳)有限公司 | 应用于髋关节置换手术的骨盆配准***和医疗设备 |
CN115005987A (zh) * | 2022-07-01 | 2022-09-06 | 北京大学第三医院(北京大学第三临床医学院) | 髋关节翻修手术中骨盆配准的方法及*** |
CN115153835A (zh) * | 2022-07-22 | 2022-10-11 | 华南理工大学 | 基于特征点配准与增强现实的髋臼假体放置引导***及方法 |
CN115100258A (zh) * | 2022-08-29 | 2022-09-23 | 杭州三坛医疗科技有限公司 | 一种髋关节图像配准方法、装置、设备以及存储介质 |
CN115661214A (zh) * | 2022-09-16 | 2023-01-31 | 杭州三坛医疗科技有限公司 | 一种配准精度验证方法及装置 |
CN115624385A (zh) * | 2022-09-19 | 2023-01-20 | 重庆生物智能制造研究院 | 术前空间配准方法及装置、计算机设备、存储介质 |
CN115607279A (zh) * | 2022-09-26 | 2023-01-17 | 苏州微创畅行机器人有限公司 | 硬组织的配准方法、装置及*** |
CN115969515A (zh) * | 2022-12-31 | 2023-04-18 | 杭州键嘉医疗科技股份有限公司 | 基于快速配准的机械臂辅助全髋关节置换手术 |
Non-Patent Citations (3)
Title |
---|
A METHOD FOR FRAME-BY-FRAME US TO CT REGISTRATION IN A JOINT CALIBRATION AND REGISTRATION FRAMEWORK;Matthias Peterhans 等;《ISBI 2008》;1131-1134 * |
基于CT图像的全髋关节置换辅助诊断算法;王琪 等;《科学技术与工程》;第19卷(第1期);183-189 * |
骨科手术机器人及其导航技术;宋国立 等;《中国科学》;第58卷;8-19 * |
Also Published As
Publication number | Publication date |
---|---|
CN116363184A (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10507002B2 (en) | X-ray system and method for standing subject | |
US20210012492A1 (en) | Systems and methods for obtaining 3-d images from x-ray information for deformed elongate bones | |
CN110464463B (zh) | 手术器械尖端标定方法、装置、计算机设备和存储介质 | |
WO2023078309A1 (zh) | 目标特征点提取方法、装置、计算机设备和存储介质 | |
CN112509119B (zh) | 针对颞骨的空间数据处理及定位方法、装置及电子设备 | |
CN110909580B (zh) | 数据处理方法及装置、电子设备及存储介质 | |
EP3414736B1 (en) | Method and system for registering a patient with a 3d image using a robot | |
CN112382359B (zh) | 患者注册方法、装置、电子设备及计算机可读介质 | |
US11628012B2 (en) | Patient positioning using a skeleton model | |
CN113397704B (zh) | 机器人定位方法、装置、***及计算机设备 | |
CN114332378B (zh) | 基于二维医学影像的人体骨骼三维模型获取方法及*** | |
CN115100258B (zh) | 一种髋关节图像配准方法、装置、设备以及存储介质 | |
CN114792326A (zh) | 一种基于结构光的手术导航点云分割与配准方法 | |
CN115345928A (zh) | 关键点的获取方法、计算机设备和存储介质 | |
CN112581422A (zh) | 用于确定用于在整形外科手术中使用的输入图像的空间取向的技术 | |
CN116363184B (zh) | 一种髋关节位置配准方法、装置、电子设备及存储介质 | |
JP6643416B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
CN116342672B (zh) | 髋关节实际位置配准方法、装置、电子设备及存储介质 | |
CN114881930A (zh) | 基于降维定位的3d目标检测方法、装置、设备和存储介质 | |
CN116645403A (zh) | 一种髋关节位置配准方法、装置、电子设备及存储介质 | |
CN115702983A (zh) | 放疗设备的摆位控制方法、设备、***及介质 | |
CN117618107A (zh) | 髋关节位置配准方法、***、装置、设备及存储介质 | |
CN116012526B (zh) | 一种基于二维影像的三维ct影像病灶重建方法 | |
CN114176779B (zh) | 一种手术机器人导航定位方法及装置 | |
CN117679060B (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 |