CN113729945A - 体表定位装置的配准方法、穿刺引导方法及设备 - Google Patents
体表定位装置的配准方法、穿刺引导方法及设备 Download PDFInfo
- Publication number
- CN113729945A CN113729945A CN202110977438.5A CN202110977438A CN113729945A CN 113729945 A CN113729945 A CN 113729945A CN 202110977438 A CN202110977438 A CN 202110977438A CN 113729945 A CN113729945 A CN 113729945A
- Authority
- CN
- China
- Prior art keywords
- matrix
- tracking
- body surface
- image space
- calculating
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 73
- 230000003287 optical effect Effects 0.000 claims abstract description 37
- 238000013519 translation Methods 0.000 claims abstract description 15
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 8
- 230000001174 ascending effect Effects 0.000 claims description 6
- 238000002591 computed tomography Methods 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 4
- 230000004807 localization Effects 0.000 claims description 2
- 230000000241 respiratory effect Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 230000033001 locomotion Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000000429 assembly Methods 0.000 description 4
- 230000000712 assembly Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000741 silica gel Substances 0.000 description 2
- 229910002027 silica gel Inorganic materials 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012943 hotmelt Substances 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 239000007779 soft material Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B17/34—Trocars; Puncturing needles
- A61B17/3403—Needle locating or guiding means
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/10—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges for stereotaxic surgery, e.g. frame-based stereotaxis
- A61B90/11—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges for stereotaxic surgery, e.g. frame-based stereotaxis with guides for needles or instruments, e.g. arcuate slides or ball joints
- A61B90/13—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges for stereotaxic surgery, e.g. frame-based stereotaxis with guides for needles or instruments, e.g. arcuate slides or ball joints guided by light, e.g. laser pointers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/361—Image-producing devices, e.g. surgical cameras
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/107—Visualisation of planned trajectories or target regions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/108—Computer aided selection or customisation of medical implants or cutting guides
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2055—Optical tracking systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2065—Tracking using image or pattern recognition
Landscapes
- Health & Medical Sciences (AREA)
- Surgery (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Robotics (AREA)
- Optics & Photonics (AREA)
- Physics & Mathematics (AREA)
- Gynecology & Obstetrics (AREA)
- Radiology & Medical Imaging (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供一种体表定位装置的配准方法、穿刺引导方法及设备，所述配准方法包括：获取布置于体表的光学跟踪装置中各个跟踪点的图像空间坐标；采集定位装置根据所述光学跟踪装置中各个跟踪点的当前位置所生成的手术空间坐标；根据所述图像空间坐标和所述手术空间坐标，确定所述各个跟踪点在图像空间和手术空间中的对应关系，并计算旋转矩阵和/或平移矩阵；根据所述各个跟踪点的图像空间坐标、所述旋转矩阵和所述平移矩阵计算光学跟踪装置的配准误差数据。
Description
技术领域
本发明涉及医学图像数据处理领域，具体涉及一种体表定位装置的配准方法、穿刺引导方法及设备。
背景技术
手术导航定位系统可通过光学跟踪的方式来实现定位，系统中设有配套的光学跟踪装置，该装置中设有反光小球，光学系统通过追踪反光小球来实现对目标的实时跟踪定位。
在目前临床使用的手术导航定位系统中，比如在骨科和神经外科领域，会把跟踪装置固定在骨骼或者头颅框架上，与手术部位保持贴合，导航定位系统通过追踪对应的跟踪装置来实现对手术部位的定位。
对于胸腹部位置，人的呼吸运动会导致体表有起伏运动，这种运动会导致定位产生偏差。对于穿刺术来说，需要先在术前根据CT图像确定穿刺路径，然后在术中按照预定的穿刺路径实施穿刺，但是受到人体呼吸导致体表起伏的影响，实际的穿刺路径可能会与预定穿刺路径不同，本申请将其称之为配准误差。配准误差过大将导致穿刺针无法准确抵达靶点，因此医生实施穿刺术的难度较大。
发明内容
有鉴于此，本发明提供一种体表定位装置的配准方法，所述体表定位装置包括多个跟踪点，且任意对跟踪点的间距不同，所述方法包括：
获取布置于体表的光学跟踪装置中各个跟踪点的图像空间坐标；
采集定位装置根据所述光学跟踪装置中各个跟踪点的当前位置所生成的手术空间坐标；
根据所述图像空间坐标和所述手术空间坐标，确定所述各个跟踪点在图像空间和手术空间中的对应关系，并计算旋转矩阵和/或平移矩阵；
根据所述各个跟踪点的图像空间坐标、所述旋转矩阵和所述平移矩阵计算光学跟踪装置的配准误差数据。
可选地，所述图像空间坐标是所述各个跟踪点在计算机断层扫描图像数据中的三维坐标。
可选地，所述手术空间坐标是由双目摄像机对所述各个跟踪点进行拍摄所提供的三维坐标。
可选地，确定所述各个跟踪点在图像空间和手术空间中的对应关系，具体包括：
根据所述图像空间坐标计算图像空间中各个跟踪点的距离矩阵DQ，根据所述手术空间坐标计算手术空间中各个跟踪点的距离矩阵DP；
根据距离矩阵DQ和距离矩阵DP计算其中每一列的距离得到误差矩阵M；
根据误差矩阵M中各列中最小元素确定所述各个跟踪点在图像空间和手术空间中的对应关系。
可选地，距离矩阵DQ和距离矩阵DP中分别包括任意两两跟踪点之间的距离。
可选地，根据间距矩阵DQ和间距矩阵DP计算其中每一列的距离得到误差矩阵M，具体包括：
分别对间距矩阵DQ和间距矩阵DP进行升序排列；
分别计算升序排列后的两个矩阵中每一列的距离，得到误差矩阵M。
可选地，所述根据所述各个跟踪点的图像空间坐标、所述旋转矩阵和所述平移矩阵计算光学跟踪装置的配准误差数据的步骤中，具体利用如下方式计算所述配准误差数据FRE：
其中n表示所述跟踪点的数量，qi表示第i个跟踪点的图像空间坐标，R表示所述平移矩阵，T表示所述旋转矩阵。
可选地，所述多个跟踪点之间通过为柔性连接件串联，使得各个跟踪点的手术空间坐标能够独立地发生变化。
本发明还提供一种穿刺引导方法，包括：根据上述方法计算穿刺对象的至少一个呼吸周期内的配准误差数据，其中所述图像空间坐标为固定值，所述手术空间坐标随穿刺对象呼吸时的体表起伏发生变化，使得所述配准误差数据为动态值；确定配准误差数据的最小值对应的时刻，并以此引导穿刺动作。
本发明还提供一种电子设备，包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行上述体表定位装置的配准方法。
本发明还提供一种电子设备，包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行上述穿刺引导方法。
根据体表定位装置的配准方法、穿刺引导方法及设备，所计算出的配准误差数据可以量化表达人体呼吸体态与术前进行CT扫描时的人体呼吸体态的匹配度，由此可以让医生在合适的时机实施穿刺术，确保实际的穿刺路径与预定的穿刺路径一致，降低穿刺术的难度，提高手术效率。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中的光学跟踪装置的结构示意图；
图2为本发明实施例中的体表定位装置的配准方法的流程图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
光学跟踪装置布置在人体体表，因此也可称之为体表定位装置。手术导航系统中的双目摄像机可以采集光学跟踪装置的位置，从而精准确定人体位置以及要实施穿刺的位置。
如图1所示，本实施例提供一种光学跟踪装置，用于配合光导航定位系统使用，光导航定位系统可以是红外光学导航定位系统，也可以是可见光导航定位系统，本实施例的光学跟踪装置可以通用。具体的，本实施例的光学跟踪装置包括两个以上反光小球组件和若干个体表定位带3。
反光小球组件用于与手术导航定位系统配合定位。反光小球组件用于手术导航定位系统配合定位的工作原理属于本领域的公知技术，在此不做赘述。
若干个体表定位带3为软性结构，体表定位带3数量与反光小球组件数量相同，两个以上反光小球组件依次通过体表定位带3连接形成闭环。
作为本申请的光学跟踪装置的优选实施例，在本实施例中，反光小球组件设置为5组，相应地，体表定位带3设置为5条。应当理解的是，根据实际需求的不同，本领域技术人员还可以采用3组、4组、6组、7组反光小球组件，相应地，体表定位带3的数量为3条、4条、6条、7条。
本实施例中的反光小球组件设置为5组。经过实验发现，5组球能兼顾精度和效率的要求，简化结构和算法。
如图1所示，5组反光小球组件，每2组反光小球组件通过1条体表定位连接，5组反光小球组件依次通过体表定位带连接形成闭环。本申请对于5组反光小球组件的相互位置以及每2组反光小球组件之间连接的体表定位带3的长度不作具体限定，需要根据实际需要进行设定。
本实施例通过采用软性体表定位带，一方面在保证由若干个定位结构构成的整体定位结构的基本形状固定，另一方面由于是软性材质，各个定位结构之间又可以独立运动。这样的好处在于，当体表定位带贴合与皮肤体表时，各个定位结构之间可以动态的随体表运动，这样便可以动态追踪体表皮肤的运动幅度，从而实现动态的导航定位，解决了当前刚体定位装置只能静态追踪的问题。
为了进一步提高动态定位效果，在本实施例中，体表定位带3为具有弹性的软性结构。优选地，在本实施例中，体表定位带3采用硅胶带。但应当理解的是，凡是具有弹性的软性材质制成的体表定位带3均应属于本申请的保护范围之内。体表定位带3为软性结构，采用硅胶制作，具有良好的生物相容性，具有较好的吸附力，在人体上可以灵活调整，所以体表定位带3两端连接的定位结构可以适应于体表，可以贴合肌肤。
为了保证本实施例的体表动态定位装置放置到相应的人体表面时能够更好地贴合人体表面，提高体表动态定位的精度，在本实施例中，将体表定位带3的底面与定位结构的底面平齐，使得体表定位带3和定位结构的底面和体表贴合而不留间隙。此外，将体表定位带3的底面与定位结构的底面平齐还便于后续采用医用胶带对体表定位带3进行固定，然后便于开展实时的追踪定位。
反光小球组件与体表定位带3的连接方式可以是卡接、螺纹连接或者热熔粘接或者是其他固定连接方式，作为本申请的优选实施例，在本实施例中，反光小球组件包括：
体表基座5，用于连接体表定位带3，体表基座5的底面与其连接的体表定位带3的底面平齐，体表基座5成型有限位槽6；
体表基座5为圆柱结构，用于连接体表定位带3，体表基座5和体表定位带3数量一致，依次相互连接，构成一个闭环；体表定位带3连接在体表基座5的外侧面上，同时，为了贴合皮肤，体表基座5的底面与其连接的体表定位带3的底面平齐；
同时，在体表基座5的顶面成型有限位槽6，用来安装小球固定座2；
小球固定座2，用于安装反光小球1，小球固定座2安装到限位槽6中，其结构和限位槽6相匹配；
小球固定柱4，插入到小球固定座2中以固定小球固定座2和体表基座5。小球固定柱4可以为螺栓等组件，将限位槽6槽底打通，小球固定柱4从下方插入和配合在限位槽6内的小球固定座2连接即可。连接方式可以是同时在小球固定座2的底面设置螺纹孔，小球固定柱4从通孔进入而旋入小球固定座2底面设置的螺纹孔中。
利用手术导航系统和上述光学跟踪装置可以追踪人体呼吸动作。本发明实施例提供一种体表定位装置的配准方法，该方法可以由计算机或服务器等电子设备执行，如图2所示该方法包括如下步骤：
S1，获取布置于体表的光学跟踪装置中各个跟踪点的图像空间坐标。本实施例所述跟踪点可以但不限于是上述实施例中的反光小球。具体地，在术前可以通过如计算机断层扫描方式对佩戴着光学跟踪装置的人体部位进行扫描，从而得到CT图像数据(三维)，然后在其中识别各个跟踪点的位置，得到其在CT图像数据中的空间坐标，称之为图像空间坐标，此数据为静态值。在具体实施例中，光学跟踪装置有5个跟踪点，图像空间坐标记为Q＝[q1,q2,q3,q4,q5]，其中q1…q5表示5个跟踪点的三维坐标。
S2，采集定位装置根据光学跟踪装置中各个跟踪点的当前位置所生成的手术空间坐标。具体地，在术中可通过定位装置(导航双目摄像机)追踪各个跟踪点的当前位置，得到其在手术空间中的坐标，称之为手术空间坐标，由于人体的呼吸起伏会导致光学跟踪装置运动，因此手术空间坐标是随人的呼吸发生变化的动态值。在具体实施例中，光学跟踪装置有5个跟踪点，手术空间坐标记为P＝[p1,p2,p3,p4,p5]，其中p1…p5表示5个跟踪点的三维坐标。
S3，根据图像空间坐标和手术空间坐标，确定各个跟踪点在图像空间和手术空间中的对应关系，并计算旋转矩阵和/或平移矩阵。光学跟踪装置的各个跟踪点的图像坐标和手术空间坐标的排列是无序的，为了进行实时跟踪配准，则需要识别两个空间中各个跟踪点的对应关系。光学跟踪装置的各个跟踪点中任意2个点的两两距离都不同，这样通过点间距离的不一致进行匹配，即可确定对应匹配的点。由于CT设备和双目摄像机对同一目标的扫描和拍摄方位不同，在确定上述对应关系后，还需要计算旋转矩阵和/或平移矩阵。针对两个坐标系计算旋转矩阵和平移矩阵的算法属于公知常识，本申请不再赘述。
S4，根据各个跟踪点的图像空间坐标、旋转矩阵和/或平移矩阵计算光学跟踪装置的配准误差数据。这一配准误差数据，用于表征当前(术中)人体呼吸体态与术前进行CT扫描时的人体呼吸体态的匹配度，例如误差数据约小则越匹配，具体计算方法有多种，比如可以分别针对每个跟踪点计算配准误差数据，然后取其中最大值、最小值、平均值等等。
根据本发明实施例提供的体表定位装置的配准方法，所计算出的配准误差数据可以量化表达人体呼吸体态与术前进行CT扫描时的人体呼吸体态的匹配度，由此可以让医生在合适的时机实施穿刺术，确保实际的穿刺路径与预定的穿刺路径一致，降低穿刺术的难度，提高手术效率。
在优选的实施例中，针对柔性材料制成的体表定位带，上述步骤S4中具体利用如下方式计算配准误差数据FRE：
其中n表示跟踪点的数量(比如n＝5)，qi表示第i个跟踪点的图像空间坐标，R表示平移矩阵，T表示旋转矩阵。利用本方案计算出的配准误差数据准确性较高，能够更准确地反映出术中体表的呼吸起伏状态与术前扫描图像时的呼吸起伏状态差异。
关于上述步骤S3中手术空间中的跟踪点与图像空间中的跟踪点的匹配方式，本申请提供一种优选匹配方法：
S31，分别计算手术空间中两两跟踪点的距离和图像空间中两两跟踪点的距离。以5个跟踪点为例，手术空间中的第i个跟踪点与其它跟踪点的距离(欧式距离)记为dpi＝[‖pi-p1‖…‖pi-p5‖]，其中i＝1…5；图像空间中的第i个跟踪点与其它跟踪点的距离记为dqi＝[‖qi-q1‖…‖qi-q5‖]，其中i＝1…5。
由此可以得到图像空间中各个跟踪点的距离矩阵DQ和手术空间中各个跟踪点的距离矩阵DP，DP＝[dp1 dp2 dp3 dp4 dp5]、DQ＝[dq1 dq2 dq3 dq4 dq5]。
S32，根据距离矩阵DQ和距离矩阵DP计算其中每一列的距离得到误差矩阵M。作为优选的实施例，可以先分别对间距矩阵DQ和间距矩阵DP进行升序排列,得到新的矩阵
和
然后再分别计算升序排列后的两个矩阵中每一列的距离mij，计算公式如下：
mij＝‖dpj-dqi‖
其中dpj表示距离矩阵DQ中的第j列，dqi示距离矩阵DP中的第i列。计算出两个矩阵中各列的距离组成误差矩阵M。
S33，根据误差矩阵M中各列中最小元素确定各个跟踪点在图像空间和手术空间中的对应关系。对于5个跟踪点，得到的误差矩阵M为5*5维矩阵，寻找矩阵M中每一列中的最小元素，比如第j列中第i个元素是最小的，则代表点集P中第j个跟踪点pj与点集Q中第i个跟踪点qi是匹配的，这样即可得到两个点集中全部跟踪点的匹配对应关系。
基于上述实施例，本发明另一实施例提供了一种穿刺引导系统，包括定位装置(双目摄像机和处理器等)、穿刺装置(机械臂和穿刺针握持结构件等)和光学跟踪装置，该穿刺引导系统被配置为执行包括以下内容的操作：
利用上述实施例的配准方法，实时计算穿刺对象的至少一个呼吸周期内的配准误差数据，其中图像空间坐标为固定值，手术空间坐标随穿刺对象呼吸时的体表起伏发生变化，使得配准误差数据为动态值；
由此可以得到配准误差数据FRE随着时间周期性变化的误差曲线，该误差曲线变化的时间周期即对应人的呼吸周期。在一个呼吸周期中，当误差值最小时，则代表人当前的呼吸状态与术前扫描图像的呼吸状态最匹配。
确定配准误差数据的最小值对应的时刻，并以此引导穿刺动作。在穿刺装置对准穿刺部位时，系统便可以根据配准误差数据引导医生进行穿刺，从而有效地减少由于呼吸变形引起的定位误差。具体的引导动作可以是通过声音提醒，比如用声音来提示当前时刻的配准误差值的大小，或者是当配准误差数据小于预设阈值时发出声音提醒，以提醒医生在合适的时机将穿刺针插入人体。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (10)
1.一种体表定位装置的配准方法，其特征在于，所述体表定位装置包括多个跟踪点，且任意对跟踪点的间距不同，所述方法包括：
获取布置于体表的光学跟踪装置中各个跟踪点的图像空间坐标；
采集定位装置根据所述光学跟踪装置中各个跟踪点的当前位置所生成的手术空间坐标；
根据所述图像空间坐标和所述手术空间坐标，确定所述各个跟踪点在图像空间和手术空间中的对应关系，并计算旋转矩阵和/或平移矩阵；
根据所述各个跟踪点的图像空间坐标、所述旋转矩阵和所述平移矩阵计算光学跟踪装置的配准误差数据。
2.根据权利要求1所述的方法，其特征在于，所述图像空间坐标是所述各个跟踪点在计算机断层扫描图像数据中的三维坐标。
3.根据权利要求1所述的方法，其特征在于，所述手术空间坐标是由双目摄像机对所述各个跟踪点进行拍摄所提供的三维坐标。
4.根据权利要求1所述的方法，其特征在于，确定所述各个跟踪点在图像空间和手术空间中的对应关系，具体包括：
根据所述图像空间坐标计算图像空间中各个跟踪点的距离矩阵DQ，根据所述手术空间坐标计算手术空间中各个跟踪点的距离矩阵DP；
根据距离矩阵DQ和距离矩阵DP计算其中每一列的距离得到误差矩阵M；
根据误差矩阵M中各列中最小元素确定所述各个跟踪点在图像空间和手术空间中的对应关系。
5.根据权利要求4所述的方法，其特征在于，距离矩阵DQ和距离矩阵DP中分别包括任意两两跟踪点之间的距离。
6.根据权利要求4所述的方法，其特征在于，根据间距矩阵DQ和间距矩阵DP计算其中每一列的距离得到误差矩阵M，具体包括：
分别对间距矩阵DQ和间距矩阵DP进行升序排列；
分别计算升序排列后的两个矩阵中每一列的距离，得到误差矩阵M。
7.根据权利要求1所述的方法，其特征在于，所述根据所述各个跟踪点的图像空间坐标、所述旋转矩阵和所述平移矩阵计算光学跟踪装置的配准误差数据的步骤中，具体利用如下方式计算所述配准误差数据FRE：
其中n表示所述跟踪点的数量，qi表示第i个跟踪点的图像空间坐标，R表示所述平移矩阵，T表示所述旋转矩阵。
8.根据权利要求1-7中任一项所述的方法，其特征在于，所述多个跟踪点之间通过为柔性连接件串联，使得各个跟踪点的手术空间坐标能够独立地发生变化。
9.一种穿刺引导方法，其特征在于，包括：
根据权利要求1-8中任一项所述方法计算穿刺对象的至少一个呼吸周期内的配准误差数据，其中所述图像空间坐标为固定值，所述手术空间坐标随穿刺对象呼吸时的体表起伏发生变化，使得所述配准误差数据为动态值；
确定配准误差数据的最小值对应的时刻，并以此引导穿刺动作。
10.一种电子设备，其特征在于，包括：至少一个处理器；以及与所述至少一个处理器通信连接的存储器；其中，所述存储器存储有可被所述一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器执行如权利要求1-8中任意一项所述的体表定位装置的配准方法和/或权利要求9所述的穿刺引导方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110977438.5A CN113729945B (zh) | 2021-08-24 | 2021-08-24 | 体表定位装置的配准方法、穿刺引导方法及设备 |
| PCT/CN2022/112009 WO2023024932A1 (zh) | 2021-08-24 | 2022-08-12 | 体表定位装置的配准方法、穿刺引导方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110977438.5A CN113729945B (zh) | 2021-08-24 | 2021-08-24 | 体表定位装置的配准方法、穿刺引导方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113729945A true CN113729945A (zh) | 2021-12-03 |
| CN113729945B CN113729945B (zh) | 2022-04-15 |
Family
ID=78732656
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110977438.5A Active CN113729945B (zh) | 2021-08-24 | 2021-08-24 | 体表定位装置的配准方法、穿刺引导方法及设备 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN113729945B (zh) |
| WO (1) | WO2023024932A1 (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114404040A (zh) * | 2021-12-22 | 2022-04-29 | 苏州微创畅行机器人有限公司 | 安装构件、标记构件以及光学跟踪组件和装置 |
| WO2023024932A1 (zh) * | 2021-08-24 | 2023-03-02 | 真健康（北京）医疗科技有限公司 | 体表定位装置的配准方法、穿刺引导方法及设备 |
| CN117017496A (zh) * | 2023-09-28 | 2023-11-10 | 真健康（北京）医疗科技有限公司 | 柔性体表定位装置及穿刺手术导航定位系统 |
| CN117838311A (zh) * | 2024-03-07 | 2024-04-09 | 杭州海沛仪器有限公司 | 基于光学定位的靶点消融呼吸门控方法及系统 |
| CN117838311B (zh) * | 2024-03-07 | 2024-05-31 | 杭州海沛仪器有限公司 | 基于光学定位的靶点消融呼吸门控系统 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1775183A (zh) * | 2005-12-07 | 2006-05-24 | 嘉兴市第一医院 | 有框架立体定向神经外科手术系统注册方法 |
| CN1806771A (zh) * | 2006-01-26 | 2006-07-26 | 清华大学深圳研究生院 | 计算机辅助经皮肾穿刺取石术中的穿刺导航系统及方法 |
| CN104055520A (zh) * | 2014-06-11 | 2014-09-24 | 清华大学 | 人体器官运动监测方法和人体导航系统 |
| CN106361366A (zh) * | 2016-11-02 | 2017-02-01 | 上海联影医疗科技有限公司 | 多模态图像配准方法及系统 |
| CN107874832A (zh) * | 2017-11-22 | 2018-04-06 | 合肥美亚光电技术股份有限公司 | 骨科手术器械导航系统及方法 |
| US20180263707A1 (en) * | 2015-07-21 | 2018-09-20 | Gal Sela | System and method for mapping navigation space to patient space in a medical procedure |
| CN111388092A (zh) * | 2020-03-17 | 2020-07-10 | 京东方科技集团股份有限公司 | 一种定位跟踪件、配准方法、存储介质及电子设备 |
| CN113143459A (zh) * | 2020-01-23 | 2021-07-23 | 海信视像科技股份有限公司 | 腹腔镜增强现实手术导航方法、装置及电子设备 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3443925B1 (en) * | 2014-05-14 | 2021-02-24 | Stryker European Holdings I, LLC | Processor arrangement for tracking the position of a work target |
| US10413377B2 (en) * | 2015-03-19 | 2019-09-17 | Medtronic Navigation, Inc. | Flexible skin based patient tracker for optical navigation |
| CN107468351A (zh) * | 2016-06-08 | 2017-12-15 | 北京天智航医疗科技股份有限公司 | 一种手术定位装置、定位系统及定位方法 |
| CN109091229A (zh) * | 2018-09-13 | 2018-12-28 | 上海逸动医学科技有限公司 | 适用于x光下机器人手术导航的柔性定位装置及导航方法 |
| CN113796959A (zh) * | 2021-08-24 | 2021-12-17 | 真健康（北京）医疗科技有限公司 | 体表动态定位装置和光学跟踪装置 |
| CN113729945B (zh) * | 2021-08-24 | 2022-04-15 | 真健康（北京）医疗科技有限公司 | 体表定位装置的配准方法、穿刺引导方法及设备 |
-
2021
- 2021-08-24 CN CN202110977438.5A patent/CN113729945B/zh active Active
-
2022
- 2022-08-12 WO PCT/CN2022/112009 patent/WO2023024932A1/zh unknown
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1775183A (zh) * | 2005-12-07 | 2006-05-24 | 嘉兴市第一医院 | 有框架立体定向神经外科手术系统注册方法 |
| CN1806771A (zh) * | 2006-01-26 | 2006-07-26 | 清华大学深圳研究生院 | 计算机辅助经皮肾穿刺取石术中的穿刺导航系统及方法 |
| CN104055520A (zh) * | 2014-06-11 | 2014-09-24 | 清华大学 | 人体器官运动监测方法和人体导航系统 |
| US20180263707A1 (en) * | 2015-07-21 | 2018-09-20 | Gal Sela | System and method for mapping navigation space to patient space in a medical procedure |
| CN106361366A (zh) * | 2016-11-02 | 2017-02-01 | 上海联影医疗科技有限公司 | 多模态图像配准方法及系统 |
| CN107874832A (zh) * | 2017-11-22 | 2018-04-06 | 合肥美亚光电技术股份有限公司 | 骨科手术器械导航系统及方法 |
| CN113143459A (zh) * | 2020-01-23 | 2021-07-23 | 海信视像科技股份有限公司 | 腹腔镜增强现实手术导航方法、装置及电子设备 |
| CN111388092A (zh) * | 2020-03-17 | 2020-07-10 | 京东方科技集团股份有限公司 | 一种定位跟踪件、配准方法、存储介质及电子设备 |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023024932A1 (zh) * | 2021-08-24 | 2023-03-02 | 真健康（北京）医疗科技有限公司 | 体表定位装置的配准方法、穿刺引导方法及设备 |
| CN114404040A (zh) * | 2021-12-22 | 2022-04-29 | 苏州微创畅行机器人有限公司 | 安装构件、标记构件以及光学跟踪组件和装置 |
| WO2023116821A1 (zh) * | 2021-12-22 | 2023-06-29 | 苏州微创畅行机器人有限公司 | 安装构件、标记构件以及光学跟踪组件和装置 |
| CN117017496A (zh) * | 2023-09-28 | 2023-11-10 | 真健康（北京）医疗科技有限公司 | 柔性体表定位装置及穿刺手术导航定位系统 |
| CN117017496B (zh) * | 2023-09-28 | 2023-12-26 | 真健康（北京）医疗科技有限公司 | 柔性体表定位装置及穿刺手术导航定位系统 |
| CN117838311A (zh) * | 2024-03-07 | 2024-04-09 | 杭州海沛仪器有限公司 | 基于光学定位的靶点消融呼吸门控方法及系统 |
| CN117838311B (zh) * | 2024-03-07 | 2024-05-31 | 杭州海沛仪器有限公司 | 基于光学定位的靶点消融呼吸门控系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113729945B (zh) | 2022-04-15 |
| WO2023024932A1 (zh) | 2023-03-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113729945B (zh) | 体表定位装置的配准方法、穿刺引导方法及设备 | |
| EP2192855B1 (en) | Patent breathing modeling | |
| US20230044983A1 (en) | Sequential monoscopic tracking | |
| US8792613B2 (en) | Test object for the validation of tracking in the presence of motion | |
| Sayeh et al. | Respiratory motion tracking for robotic radiosurgery | |
| EP3217884B1 (en) | Method and apparatus for determining or predicting the position of a target | |
| US20140093160A1 (en) | 3D Object Tracking in Multiple 2D Sequences | |
| CN101983035B (zh) | 呼吸确定装置 | |
| CN105919612A (zh) | 混合多行检测器和平板成像系统 | |
| US20120310075A1 (en) | Medical measurement system and method | |
| CN105167788A (zh) | 双影像c臂系统 | |
| CN113796959A (zh) | 体表动态定位装置和光学跟踪装置 | |
| JP6692817B2 (ja) | 対象物体の変位を計算する方法及びシステム | |
| Patete et al. | Motion compensation in hand-held laser scanning for surface modeling in plastic and reconstructive surgery | |
| CN107847216A (zh) | 对肺癌辐射的指导 | |
| Kim et al. | Development of an optical‐based image guidance system: Technique detecting external markers behind a full facemask | |
| Li et al. | Geometric performance and efficiency of an optical tracking system for daily pre-treatment positioning in pelvic radiotherapy patients | |
| Hoisak et al. | A history of surface guidance methods in radiation therapy | |
| Talbot et al. | A method for patient set-up guidance in radiotherapy using augmented reality | |
| CN201067408Y (zh) | 可实时监测位置的n型线立体定位框架 | |
| Choi et al. | Real-time tumor motion tracking in 3D using planning 4D CT images during image-guided radiation therapy | |
| Li et al. | C-arm based image-guided percutaneous puncture of minimally invasive spine surgery | |
| JP6799292B2 (ja) | 放射線撮影装置および放射線画像検出方法 | |
| Jang et al. | Smart Glasses Cannula Guide System for Interventional Cardiology Procedures | |
| Chai et al. | Real-time monitoring system with accelerator controlling: An improvement of radiotherapy monitoring based on binocular location and classification |
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: Chen Xiaogang Inventor after: An Xueliang Inventor after: Shi Jipeng Inventor after: He Wanliang Inventor after: Zhang Qiangqiang Inventor after: Chen Xiangqian Inventor before: Chen Xiangqian Inventor before: An Xueliang Inventor before: Shi Jipeng Inventor before: He Wanliang Inventor before: Zhang Qiangqiang Inventor before: Chen Xiaogang |
|
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Xiaogang Inventor after: Chen Xiangqian Inventor after: An Xueliang Inventor after: Shi Jipeng Inventor after: He Wanliang Inventor after: Zhang Qiangqiang Inventor before: Chen Xiaogang Inventor before: An Xueliang Inventor before: Shi Jipeng Inventor before: He Wanliang Inventor before: Zhang Qiangqiang Inventor before: Chen Xiangqian |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230330 Address after: Room 105-75274, No. 6, Baohua Road, Hengqin New District, Zhuhai City, Guangdong Province, 519000 (centralized office area) Patentee after: True Health (Zhuhai) Medical Technology Co.,Ltd. Address before: 100192 Room 308, floor 3, building 2, Tiandi Linfeng, No. 1, yongtaizhuang North Road, Haidian District, Beijing Patentee before: ZHENJIANKANG (BEIJING) MEDICAL TECHNOLOGY Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20231027 Address after: Room 806, Research Headquarters Building, Guangdong Macao Cooperative Traditional Chinese Medicine Technology Industrial Park, No. 1 Doukou Road, Hengqin New District, Zhuhai City, Guangdong Province, 519000 Patentee after: True Health (Guangdong Hengqin) Medical Technology Co.,Ltd. Address before: Room 105-75274, No. 6, Baohua Road, Hengqin New District, Zhuhai City, Guangdong Province, 519000 (centralized office area) Patentee before: True Health (Zhuhai) Medical Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right |