CN117257459B - 一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 - Google Patents
一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 Download PDFInfo
- Publication number
- CN117257459B CN117257459B CN202311561950.7A CN202311561950A CN117257459B CN 117257459 B CN117257459 B CN 117257459B CN 202311561950 A CN202311561950 A CN 202311561950A CN 117257459 B CN117257459 B CN 117257459B
- Authority
- CN
- China
- Prior art keywords
- airway
- point
- map
- registration
- preparation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000000241 respiratory effect Effects 0.000 title claims description 28
- 238000013276 bronchoscopy Methods 0.000 title claims description 16
- 238000002360 preparation method Methods 0.000 claims abstract description 109
- 230000009466 transformation Effects 0.000 claims abstract description 74
- 239000011159 matrix material Substances 0.000 claims abstract description 67
- 238000010276 construction Methods 0.000 claims abstract description 16
- 238000013507 mapping Methods 0.000 claims description 49
- 230000004927 fusion Effects 0.000 claims description 17
- 238000003860 storage Methods 0.000 claims description 12
- 230000007704 transition Effects 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000009499 grossing Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 abstract description 13
- 230000029058 respiratory gaseous exchange Effects 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 11
- 230000003097 anti-respiratory effect Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 210000000621 bronchi Anatomy 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000001574 biopsy Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013425 morphometry Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 210000003437 trachea Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
- A61B1/00—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
- A61B1/267—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor for the respiratory tract, e.g. laryngoscopes, bronchoscopes
- A61B1/2676—Bronchoscopes
-
- 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
- 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/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2051—Electromagnetic tracking systems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Robotics (AREA)
- Pulmonology (AREA)
- Otolaryngology (AREA)
- Physiology (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Endoscopes (AREA)
Abstract
本发明公开了一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统,该方法将电磁导航支气管镜在工作过程中的状态分为三种,即常规状态、建图预备状态和建图状态,常规状态下进行实时配准并实时更新配准坐标变换矩阵,确定当前配准点在地图中的坐标,在建图预备状态同时形成预备配准点列表,建图状态下利用保留气道最近点、预备配准点列表和实时记录的配准点更新气道点列表;将更新的气道点列表与地图信息进行融合,得到新的地图信息;本发明的方法能够在有呼吸及其他干扰的情况下区分器械是在已有地图信息的气道内还是在无地图信息的气道内,并且在记录新气道信息时排除呼吸干扰对新气道形态的影响。
Description
技术领域
本发明属于手术导航技术领域,涉及一种支气管镜导航方法,尤其涉及一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统。
背景技术
支气管镜是呼吸系统疾病的重要诊疗器械之一,用于对支气管腔进行观察以及进行肺部组织活检等。支气管镜导航技术可以在支气管镜手术过程中指明器械末端在患者支气管中的位置,并引导到待观察目标以及活检目标处。
目前的支气管镜导航技术,按照利用的信息模态,可分为电磁导航支气管镜、图像导航支气管镜,以及混合模态导航支气管镜等。电磁导航支气管镜是相对可靠并被广泛使用的技术。
电磁导航支气管镜在器械中装有磁定位传感器,其实现过程包括:
1)对器械进行手眼标定,得到器械末端坐标系和磁定位传感器坐标系/>之间的坐标变换矩阵(以下简称“手眼标定坐标变换矩阵”)/>。两个坐标系的示意图如图1所示。
和/>分别是手眼标定中对支气管镜末端位姿和磁定位传感器位姿的采样信息,n为收集的位姿数量。对/>的求解如公式(1)所示。
(1)
的计算方式为首先用奇异值分解的方法进行初始估计,然后用迭代最近点的方法进行优化求解。
若已知器械末端的位姿,即可由公式(2)计算出磁定位传感器的位姿/>
(2)
2)对术前CT数据进行处理,得到支气管树的三维形态描述信息(以下简称“地图信息”),其形式为支气管树的表面信息,或支气管树的表面信息和支气管树中各气道中心线信息。根据地图信息,可以在计算机图形界面上显示支气管树的三维形态图(以下简称“地图”)。
3)手术开始时,进行初始配准。用支气管镜携带器械在支气管中移动,将器械末端依次放置在多个位置,按照磁定位传感器记录的多个坐标,以及与器械末端放置位置最接近的地图特征的坐标/>,和手眼标定坐标变换矩阵/>,计算出初始配准坐标变换矩阵/>。其坐标变换关系如公式(3)所示。
(3)
m为初始配准收集的坐标数量。的计算方式与/>的计算方式相同。
初始配准坐标变换矩阵描述磁定位传感器位姿/>和地图坐标系中的磁定位传感器位姿/>之间的变换,如公式(4)所示,/>为器械末端位姿。
(4)
4)在手术过程中,由公式(5-6)进行实时配准。
得到手眼标定坐标变换矩阵和初始配准坐标变换矩阵/>后,即可根据公式(5),由手术过程中磁定位传感器输出的实时坐标/>(以下简称“磁导航轨迹点”),和支气管树中各气道中心线上的点/>进行配准,得到配准坐标变换矩阵/>。k为参与配准的点数量。使用icp算法求解公式(5),以/>作为/>的初始估计。
得到后,由公式(6)计算得到当前磁导航轨迹点对应的配准后的地图坐标系坐标/>(以下简称“配准点”),并在计算机的图形界面显示。此配准点在地图中的坐标,即代表器械末端在患者支气管腔中的位置。
(5)
(6)
电磁导航支气管镜使用的配准技术有多种,例如构建代价函数来计算最佳的配准坐标变换矩阵,或者综合磁定位传感器输出的坐标和角度信息,与支气管树中各气道中心线信息进行ICP配准等。
支气管镜手术由支气管镜镜体和器械耗材配合完成。器械耗材包括可视探头、活检钳等(下面简称“器械”)。器械由母镜的器械通道伸出,可以携带磁定位传感器进入更细的支气管。如果器械能够进入对CT数据进行三维重建时无法分辨的气道,就能够在手术中记录新的地图信息(下面简称“新气道信息”)。
由于呼吸及其他干扰,导航过程中支气管会发生形变,与地图信息有差异。如果没有考虑呼吸干扰,按照配准点在地图中的位置来判断是否进入了新气道以及记录新气道的形状,就出现误判和记录不准确。
本发明提出一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法,能够在有呼吸及其他干扰的情况下区分器械是在已有地图信息的气道内还是在无地图信息的气道内,并且在记录新气道信息时排除呼吸干扰对新气道形态的影响。
发明内容
本发明实施例提供了抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统,可以根据经配准坐标变换矩阵变换后的磁导航轨迹,在术中记录新气道信息,更新地图信息。
第一方面,本发明实施例提供一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法,该方法包括:
针对电磁导航支气管镜,进行初始配准得到配准坐标变换矩阵;
在常规状态下,进行实时配准并实时更新配准坐标变换矩阵;由实时的配准坐标变换矩阵及磁导航轨迹点确定当前配准点在地图中的坐标;当满足建图预备判定条件,则记录满足建图预备判定条件之前最后时刻支气管树表面或支气管树各气道中心线中与配准点最近的点,定义为保留气道最近点,之后转移到建图预备状态;
在建图预备状态下,继续进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算当前配准点在地图中的坐标;同时记录从进入建图预备状态时刻开始所有的配准点,并保存在预备配准点列表中;当满足建图预备结束判定条件,则删除对保留气道最近点和预备配准点列表的记录信息,并转移到常规状态;当满足建图判定条件,则记录满足建图判定条件之前最后时刻的配准坐标变换矩阵,以及由此刻磁传感器记录导出的代表器械朝向的方向向量,并转移到建图状态;
在建图状态下,配准坐标变换矩阵不再更新,并由当前保留的配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标,利用保留气道最近点、预备配准点和实时记录的配准点,维护一个表示气道的三维点列表即气道点列表;每次更新气道点列表,将其与地图信息进行融合,得到新的地图信息;当满足建图结束判定条件,则停止对气道点列表的更新,将气道点列表和地图信息最后一次融合的结果作为新的地图信息;转移到常规状态,并恢复对配准坐标变换矩阵的更新。
第二方面,本发明实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现任一实施例所述的抗呼吸干扰的电磁导航支气管镜术中地图拓展方法。
第三方面,本发明实施例提供了一种电磁导航支气管镜系统,所述系统包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如任一实施例所述的抗呼吸干扰的电磁导航支气管镜术中地图拓展方法。
本发明实施例所提供的技术方案,通过将电磁导航支气管镜在工作过程中的状态分为三种,即常规状态、建图预备状态和建图状态,通过判定条件来确定状态是否转移;由于对建图预备状态和建图状态的判定不仅考虑了配准点和最近气道点的距离,还考虑了配准器械朝向向量和从气道最近点指向配准点的向量V之间的夹角,当配准点和最近气道点的距离较大时,如果夹角接近直角则认为器械末端仍在有地图信息的气道内,如果夹角与直角差异较大则认为器械末端朝向与最近气道朝向夹角较大或器械末端已超出最近气道末端一段距离,进而认为器械进入了无地图信息的气道,排除了当器械末端在有地图信息的气道内而由于呼吸干扰使配准点和最近气道点距离增大造成对建图预备状态和建图状态出现误判的情况;由于在记录新气道信息时将预备配准点列表中的点和实时记录的配准点向气道朝向向量上做投影,得到的投影点组成气道点列表,排除了呼吸干扰引起的预备配准点和实时记录的配准点位置震荡对新气道形态的影响。
附图说明
图1为器械末端坐标系和磁定位传感器坐标系的示意图;
图2为本发明实施例提供的一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法流程图;
图3为不加抗呼吸干扰机制的地图拓展结果(图中圆点为初始配准使用的配准点);
图4为具备抗呼吸干扰机制的地图拓展结果(图中圆点为初始配准使用的配准点)。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图2是本发明实施例提供的一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法,该方法可以由软件和/或硬件的方式来实现。
如图2所示,抗呼吸干扰的电磁导航支气管镜术中地图拓展方法包括:
针对电磁导航支气管镜,在手术开始时进行初始配准得到配准坐标变换矩阵;用一个状态机来管理电磁导航支气管镜系统的工作过程,所述状态机中的状态有常规状态、建图预备状态和建图状态三种。
所述状态机的状态转移条件有建图预备判定条件、建图判定条件、建图结束判定条件、建图预备结束判定条件四种。系统在常规状态下,满足建图预备判定条件,即转移到建图预备状态,系统在建图预备状态下,满足建图判定条件,即转移到建图状态。系统在建图状态下,满足建图结束判定条件,即转移到常规状态。
在常规状态下,进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统在计算机图形界面上实时显示地图和配准点位置。
如果满足建图预备判定条件,则记录满足建图预备判定条件之前的最后时刻支气管树表面或支气管树各气道中心线中与配准点最近的点即保留气道最近点,然后进入建图预备状态。
在建图预备状态下,继续进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统记录从进入建图预备状态时刻开始所有的配准点,并保存在预备配准点列表中。系统并在计算机图形界面上实时显示地图和配准点位置。
如果满足建图预备结束判定条件,则删除对保留气道最近点和预备配准点列表的记录信息,并转移到常规状态;如果满足建图判定条件,则记录满足建图判定条件之前最后时刻的配准坐标变换矩阵,以及由此刻磁传感器记录导出的代表器械朝向的方向向量,并转移到建图状态。
在建图状态下,配准坐标变换矩阵不再更新,并由当前保留的配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统利用保留气道最近点、预备配准点列表和实时记录的配准点,维护一个表示气道的三维点列表即气道点列表。系统每次更新气道点列表,将其与地图信息进行融合,得到新的地图信息。系统在计算机图形界面上实时显示新的地图和配准点。
在建图状态下,如果满足建图结束判定条件,则停止对气道点列表的更新,将气道点列表和地图信息最后一次融合的结果作为新的地图信息,供接下来的系统工作过程使用。系统转移到常规状态,并恢复对配准坐标变换矩阵的更新。
以当前配准点为起点,以代表器械朝向的方向向量为方向,建立一个向量即配准器械朝向向量。建立从气道最近点指向当前配准点的向量V,其中气道最近点是指当前地图信息中支气管树表面或支气管树各气道中心线中与配准点最近的点;
状态机中的状态转移条件有两种定义方式,第一种定义方式是根据和V之间的点积的大小来判断状态是否转移。根据大量实验研究,经验证一种可行的判定条件具体如下:
建图预备判定条件为,在常规状态下,满足:
(7)
建图判定条件为,在建图预备状态下,满足:
(8)
建图预备结束判定条件为,在建图预备条件下,满足:
(9)
建图结束判定条件为,在建图条件下,满足:
(10)
其中,为常数,通常可取1-15,优选为3,/>为常数,通常可取3-20,优选为7,/>的取值要大于/>的取值。
第二种定义方式是根据和V之间的余弦相似度的大小,以及气道最近点与配准点之间的欧几里得距离的大小来共同判断状态是否转移。根据大量实验研究,经验证一种可行的判定条件具体如下:
当气道最近点为支气管树各气道中心线上的点时:
建图预备判定条件为,在常规状态下,满足:
(11)
建图判定条件为,在建图预备状态下,满足:
(12)
建图预备结束判定条件为,在建图预备条件下,满足:
(13)
建图结束判定条件为,在建图条件下,满足:
(14)
其中,,/>为0与1之间的常数,二者取值优选为0.36;/>,/>为常数,通常可取1-10,/>取值优选为1,/>取值优选为4,/>的取值要大于/>的取值;/>为根据当前气道最近点所在的气道分支级数对应的预计气道半径,/>为器械半径。
当气道最近点为支气管树表面上的点时:
建图预备判定条件为,在常规状态下,满足:
(15)
建图判定条件为,在建图预备状态下,满足:
(16)
建图预备结束判定条件为,在建图预备条件下,满足:
(17)
建图结束判定条件为,在建图条件下,满足:
(18)
其中,,/>为0与1之间的常数,二者取值优选为0.36;/>,/>为常数,通常可取1-10,/>取值优选为1,/>取值优选为4,/>的取值要大于/>的取值;/>为根据当前气道最近点所在的气道分支级数对应的预计气道半径,/>为器械半径。
气道点列表的计算和更新方式如下:
在进入建图状态时,以保留气道最近点为起点,以代表器械朝向的方向向量为方向,建立一个向量即气道朝向向量,将预备配准点列表中的点往气道朝向向量上做投影,得到的投影点作为气道点列表。在建图状态中,实时记录的配准点继续往气道朝向向量上做投影,得到的投影点加入到气道点列表中。
实施例1
对器械进行手眼标定,得到手眼标定坐标变换矩阵。其中的坐标变换关系如公式(2)所示。
的计算方式为首先用奇异值分解的方法进行初始估计,然后用迭代最近点的方法进行优化求解。
对术前CT数据进行处理,得到地图信息,其形式为支气管树的表面信息。
手术开始时,进行初始配准,得到配准坐标变换矩阵。已知手眼标定坐标变换矩阵/>和配准坐标变换矩阵/>,即可根据公式(4),由磁定位传感器输出的实时坐标/>,计算得到配准点在地图中的坐标/>,并在计算机的图形界面显示。此配准点在地图中的位置,即代表器械末端在患者支气管腔中的位置。
初始配准后,用一个状态机来管理所述电磁导航支气管镜系统的工作过程,如图1所示。系统在工作过程中的状态分为常规状态、建图预备状态和建图状态。
系统在常规状态下,进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统在计算机图形界面上显示地图和配准点位置。
系统在常规状态下,如果满足建图预备判定条件(公式7),则记录满足建图预备判定条件之前的最后时刻支气管树表面与配准点最近的点,作为保留气道最近点,然后进入建图预备状态。本实施例使用的参数为3,此参数也可以为其他1与15之间的值;系统在建图预备状态下,继续进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统记录从进入建图预备状态时刻开始所有的配准点,并保存在预备配准点列表中。系统并在计算机图形界面上显示地图和配准点位置。
系统在建图预备状态下,如果满足建图预备结束判定条件(公式9),则删除对保留气道最近点和预备配准点列表的记录信息,并转移到常规状态。
系统在建图预备状态下,如果满足建图判定条件(公式8),则记录满足建图判定条件之前最后时刻的配准坐标变换矩阵,以及由此刻磁传感器记录导出的代表器械朝向的方向向量,并转移到建图状态。本实施例使用的参数为7,此参数也可以为其他3与20之间的值。
系统在建图状态下,配准坐标变换矩阵不再更新,并由当前保留的配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统利用保留气道最近点、预备配准点和实时记录的配准点,维护气道点列表。系统每次更新气道点列表,将其与地图信息进行融合,得到新的地图信息。系统在计算机图形界面上显示新的地图和配准点。
系统在建图状态下,如果满足建图结束判定条件(公式10),则停止对气道点列表的更新,将气道点列表和地图信息最后一次融合的结果作为新的地图信息,供接下来的系统工作过程使用。系统转移到常规状态,并恢复对配准坐标变换矩阵的更新。
气道点列表的计算和更新方式如下:
在进入建图状态时,以保留气道最近点为起点,以代表器械朝向的方向向量为方向,建立一个向量即气道朝向向量,将预备配准点列表中的点往气道朝向向量上做投影,得到的投影点作为气道点列表。在建图状态中,实时记录的配准点继续往气道朝向向量上做投影,得到的投影点加入到气道点列表中。
每次更新气道点列表时,将气道点列表与支气管树的表面信息进行融合并在计算机图形界面显示,融合的方式如下:
1. 利用气道点列表,进行平滑和插值,得到一个新的点列表。
2. 依次连接上述新的点列表,得到一个线结构。将此线结构作为轴线,设置半径r,生成一个两端封闭的管道结构。
3. 将此管道结构与支气管树的表面信息做布尔并运算,得到的结果即是气道点列表与地图信息融合的结果。
实施例2
步骤同实施例1,不同之处在于:
建图预备判定条件如公式(15)。本实施例使用的参数为0.36,此参数也可以为其他0与1之间的值;使用的参数/>为1,此参数也可以为其他1与10之间的值。
为根据当前气道最近点所在的气道分支级数对应的预计气道半径,所依据的气道分支模型参见Morphometry of the human lung[M]. Berlin: Springer, 1963.,/>的定义方式如下表所示:
表 1 分支级数与对应预计气道半径的关系
分支级数 | 气道半径(mm) |
1 | 8 |
2 | 6 |
3 | 4 |
4 | 3 |
5 | 2.5 |
>5 | 2 |
建图预备结束判定条件如公式(17)。
建图判定条件如公式(16)。本实施例使用的参数为0.36,此参数也可以为其他0与1之间的值;使用的参数/>为4,此参数也可以为其他1与10之间的值。
建图结束判定条件如公式(18)。
实施例3
对器械进行手眼标定,得到手眼标定坐标变换矩阵。其中的坐标变换关系如公式(2)所示。
的计算方式为首先用奇异值分解的方法进行初始估计,然后用迭代最近点的方法进行优化求解。
对术前CT数据进行处理,得到地图信息,其形式为支气管树中各气道中心线信息和支气管树的表面信息。
手术开始时,进行初始配准,得到配准坐标变换矩阵。已知手眼标定坐标变换矩阵/>和配准坐标变换矩阵/>,即可根据公式(4),由磁定位传感器输出的实时坐标/>,计算得到配准点在地图中的坐标/>,并在计算机的图形界面显示。此配准点在地图中的位置,即代表器械末端在患者支气管腔中的位置。
初始配准后,用一个状态机来管理所述电磁导航支气管镜系统的工作过程,如图1所示。系统在工作过程中的状态分为常规状态、建图预备状态和建图状态。
系统在常规状态下,进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统在计算机图形界面上显示地图和配准点位置。
系统在常规状态下,如果满足建图预备判定条件(公式7),则记录满足建图预备判定条件之前的最后时刻支气管树各气道中心线中与配准点最近的点,作为保留气道最近点,然后进入建图预备状态。本实施例使用的参数为3,此参数也可以为其他1与15之间的值。
系统在建图预备状态下,继续进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统记录从进入建图预备状态时刻开始所有的配准点,并保存在预备配准点列表中。系统并在计算机图形界面上显示地图和配准点位置。
系统在建图预备状态下,如果满足建图预备结束判定条件(公式9),则删除对保留气道最近点和预备配准点列表的记录信息,并转移到常规状态。
系统在建图预备状态下,如果满足建图判定条件(公式8),则记录满足建图判定条件之前最后时刻的配准坐标变换矩阵,以及由此刻磁传感器记录导出的代表器械朝向的方向向量,并转移到建图状态。本实施例使用的参数为7,此参数也可以为其他3与20之间的值。
系统在建图状态下,配准坐标变换矩阵不再更新,并由当前保留的配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标。系统利用保留气道最近点、预备配准点和实时记录的配准点,维护气道点列表。系统每次更新气道点列表,将其与地图信息进行融合,得到新的地图信息。系统在计算机图形界面上显示新的地图和配准点。
系统在建图状态下,如果满足建图结束判定条件(公式10),则停止对气道点列表的更新,将气道点列表和地图信息最后一次融合的结果作为新的地图信息,供接下来的系统工作过程使用。系统转移到常规状态,并恢复对配准坐标变换矩阵的更新。
气道点列表的计算和更新方式如下:
在进入建图状态时,以保留气道最近点为起点,以代表器械朝向的方向向量为方向,建立一个向量即气道朝向向量,将预备配准点列表中的点往气道朝向向量上做投影,得到的投影点作为气道点列表。在建图状态中,实时记录的配准点继续往气道朝向向量上做投影,得到的投影点加入到气道点列表中。
每次更新气道点列表时,将气道点列表与支气管树的表面信息、支气管树中各气道中心线信息进行融合并在计算机图形界面显示,融合的方式如下:
1. 气道点列表与支气管树的表面信息的融合
1)利用气道点列表,进行平滑和插值,得到一个新的点列表。
2)依次连接上述新的点列表,得到一个线结构。将此线结构作为轴线,设置半径r,生成一个两端封闭的管道结构。
3)将此管道结构与支气管树的表面信息做布尔并运算,得到的结果即是气道点列表与地图信息融合的结果。
2. 气道点列表与支气管树中各气道中心线信息的融合
1) 利用气道点列表,进行平滑和插值,得到一个新的点列表。
2) 依次连接上述新的点列表,得到一个线结构。
3) 将此线结构加入到支气管树中各气道中心线信息中,得到的结果即是气道点列表与支气管树中各气道中心线信息的融合结果。
实施例4
步骤同实施例3,不同之处在于:
建图预备判定条件如公式(11)。本实施例使用的参数为0.36,此参数也可以为其他0与1之间的值;使用的参数/>为1,此参数也可以为其他1与10之间的值。
为根据当前保留气道最近点所在的气道分支级数而假设的气道半径,其定义方式如表1所示。
建图预备结束判定条件如公式(13)。
建图判定条件如公式(12)。本实施例使用的参数为0.36,此参数也可以为其他0与1之间的值;使用的参数/>为4,此参数也可以为其他1与10之间的值。
建图结束判定条件如公式(14)。
为了验证本发明方法的有效性,设计了实验对本发明的抗呼吸干扰的术中地图拓展方法重建气道的误差进行评价。其实验方法如下:
首先,从算法用气道树中心线数据中删除需要重建的中心线分支,并把删除的中心线分支数据作为标准值。
然后用支气管镜携带磁定位传感器在假体内移动,进入需要重建的中心线分支,然后退出,此过程中录制磁导航数据。
在磁导航数据上加以主隆嵴为中心的呼吸干扰,每点的运动幅度和此点与主隆嵴距离成正比。假设呼吸过程中气管没有被伸长,因此在主隆嵴的气管一侧不加横断面方向的呼吸干扰。
公式(19-20)为在磁导航数据上加呼吸干扰的公式。公式(19)中,为磁导航坐标,X为冠状面的方向,Y为横断面的方向,Z为矢状面的方向。为主隆嵴位置坐标。n为呼吸干扰幅度与点/>和主隆嵴距离的比例,i为磁导航数据编号。
(19)
M为一个呼吸周期的磁导航数据量,其计算如公式(20),为呼吸频率,/>为磁导航数据采样频率。
(20)
使用加呼吸干扰后的磁导航数据作为导航算法输入,在导航过程中重建出气道,将气道点列表和对应的中心线分支标准值相比较,为气道点列表中每一点找到中心线分支标准值中的最近点,并计算平均距离,作为气道重建的误差值。
在1条中心线分支上进行此实验过程,每个呼吸幅度下进行3次重建取平均值,并把无抗呼吸干扰的算法结果和加抗呼吸干扰的算法结果相比较。这里无抗呼吸干扰的算法是指状态机的状态转移条件只考虑气道最近点与配准点之间的欧几里得距离,并且在建立和更新气道点列表时直接使用预备配准点列表中的点和实时记录的配准点的原始坐标。对应的气道重建误差如表2所示,可以看出加抗呼吸干扰的算法气道重建误差比无抗呼吸干扰的算法气道重建误差有明显的减小。
表 2 有呼吸干扰下的气道重建误差
在呼吸干扰幅度为0.2时,对无抗呼吸干扰的算法结果和加抗呼吸干扰的算法结果各取截图如图3-4所示,可以看出抗呼吸干扰算法的效果是积极的,一是新气道起始位置不再有横向偏移,二是不会因为呼吸运动而使新气道出现折叠的形状。
本申请实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述任一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法。
本申请实施例还提供一种电磁导航支气管镜系统,所述系统包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现上述任一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明并不限于上述实施例,在本领域技术人员所具有的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。比如:将其扩展应用至多机协同探索系统。
Claims (14)
1.一种电磁导航支气管镜系统,其特征在于,所述系统包括:
一个或多个处理器;
存储器,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法,所述方法包括:
针对电磁导航支气管镜,进行初始配准得到配准坐标变换矩阵;
在常规状态下,进行实时配准并实时更新配准坐标变换矩阵;由实时的配准坐标变换矩阵及磁导航轨迹点确定当前配准点在地图中的坐标;当满足建图预备判定条件,则记录满足建图预备判定条件之前最后时刻支气管树表面或支气管树各气道中心线中与配准点最近的点,定义为保留气道最近点,之后转移到建图预备状态;
在建图预备状态下,继续进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算当前配准点在地图中的坐标;同时记录从进入建图预备状态时刻开始所有的配准点,并保存在预备配准点列表中;当满足建图预备结束判定条件,则删除对保留气道最近点和预备配准点列表的记录信息,并转移到常规状态;当满足建图判定条件,则记录满足建图判定条件之前最后时刻的配准坐标变换矩阵,以及由此刻磁传感器记录导出的代表器械朝向的方向向量,并转移到建图状态;
在建图状态下,配准坐标变换矩阵不再更新,并由当前保留的配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标,利用保留气道最近点、预备配准点和实时记录的配准点,维护一个表示气道的三维点列表即气道点列表;每次更新气道点列表,将其与地图信息进行融合,得到新的地图信息;当满足建图结束判定条件,则停止对气道点列表的更新,将气道点列表和地图信息最后一次融合的结果作为新的地图信息;转移到常规状态,并恢复对配准坐标变换矩阵的更新。
2.根据权利要求1所述的电磁导航支气管镜系统,其特征在于,方法中所涉及的判定条件按下述方法得到:
以当前配准点为起点,以此刻磁传感器记录导出的代表器械朝向的方向向量为方向,建立一个向量即配准器械朝向向量,V为从气道最近点指向当前配准点的向量,所述气道最近点是指当前支气管树表面或支气管树各气道中心线中与配准点最近的点;不同状态下的判定条件由配准器械朝向向量/>和从气道最近点指向配准点的向量V之间的点积、余弦相似度和相关性,以及气道最近点和配准点之间的欧几里得距离中的一种或几种的组合与阈值相比较而得到。
3.根据权利要求2所述的电磁导航支气管镜系统,其特征在于,所述的判定条件为:
建图预备判定条件为,在常规状态下,满足:
,
建图判定条件为,在建图预备状态下,满足:
,
建图预备结束判定条件为,在建图预备条件下,满足:
,
建图结束判定条件为,在建图条件下,满足:
,
其中,,/>为常数,/>的取值要大于/>的取值。
4.根据权利要求2所述的电磁导航支气管镜系统,其特征在于,所述的判定条件为:
当气道最近点为支气管树各气道中心线上的点时:
建图预备判定条件为,在常规状态下,满足:
,
建图判定条件为,在建图预备状态下,满足:
,
建图预备结束判定条件为,在建图预备条件下,满足:
,
建图结束判定条件为,在建图条件下,满足:
,
其中,,/>为0与1之间的常数;/>,/>为常数,/>的取值要大于的取值;/>为根据当前气道最近点所在的气道分支级数对应的预计气道半径,/>为器械半径;
当气道最近点为支气管树表面上的点时:
建图预备判定条件为,在常规状态下,满足:
,
建图判定条件为,在建图预备状态下,满足:
,
建图预备结束判定条件为,在建图预备条件下,满足:
,
建图结束判定条件为,在建图条件下,满足:
,
其中,,/>为0与1之间的常数;/>,/>为常数,/>的取值要大于的取值;/>为根据当前气道最近点所在的气道分支级数对应的预计气道半径,/>为器械半径。
5.根据权利要求1所述的电磁导航支气管镜系统,其特征在于,所述的气道点列表及其更新具体为:
在进入建图状态时,以保留气道最近点为起点,以代表器械朝向的方向向量为方向,建立一个向量即气道朝向向量,将预备配准点列表中的点向气道朝向向量上做投影,得到的投影点作为气道点列表,在建图状态中,实时记录的配准点继续向气道朝向向量上做投影,得到的投影点加入到气道点列表中。
6.根据权利要求1所述的电磁导航支气管镜系统,其特征在于,每次更新气道点列表时,将气道点列表与地图信息进行融合得到新的地图信息,当地图信息是支气管树的表面信息,融合的方式如下:
利用气道点列表,进行平滑和插值,得到一个新的点列表;
依次连接上述新的点列表,得到一个线结构,将此线结构作为轴线,设置半径r,生成一个两端封闭的管道结构;
将此管道结构与支气管树的表面信息做布尔并运算,得到的结果即是气道点列表与地图信息融合的结果,即新的地图信息。
7.根据权利要求1所述的电磁导航支气管镜系统,其特征在于,每次更新气道点列表时,将气道点列表与地图信息进行融合得到新的地图信息,如果地图信息是支气管树中各气道中心线信息,融合的方式如下:
利用气道点列表,进行平滑和插值,得到一个新的点列表;
依次连接上述新的点列表,得到一个线结构;
将此线结构加入到支气管树中各气道中心线信息中,得到的结果即是气道点列表与地图信息融合的结果,即新的地图信息。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法,所述方法包括:
针对电磁导航支气管镜,进行初始配准得到配准坐标变换矩阵;
在常规状态下,进行实时配准并实时更新配准坐标变换矩阵;由实时的配准坐标变换矩阵及磁导航轨迹点确定当前配准点在地图中的坐标;当满足建图预备判定条件,则记录满足建图预备判定条件之前最后时刻支气管树表面或支气管树各气道中心线中与配准点最近的点,定义为保留气道最近点,之后转移到建图预备状态;
在建图预备状态下,继续进行实时配准,更新配准坐标变换矩阵,并由配准坐标变换矩阵和磁导航轨迹点计算当前配准点在地图中的坐标;同时记录从进入建图预备状态时刻开始所有的配准点,并保存在预备配准点列表中;当满足建图预备结束判定条件,则删除对保留气道最近点和预备配准点列表的记录信息,并转移到常规状态;当满足建图判定条件,则记录满足建图判定条件之前最后时刻的配准坐标变换矩阵,以及由此刻磁传感器记录导出的代表器械朝向的方向向量,并转移到建图状态;
在建图状态下,配准坐标变换矩阵不再更新,并由当前保留的配准坐标变换矩阵和磁导航轨迹点计算配准点在地图中的坐标,利用保留气道最近点、预备配准点和实时记录的配准点,维护一个表示气道的三维点列表即气道点列表;每次更新气道点列表,将其与地图信息进行融合,得到新的地图信息;当满足建图结束判定条件,则停止对气道点列表的更新,将气道点列表和地图信息最后一次融合的结果作为新的地图信息;转移到常规状态,并恢复对配准坐标变换矩阵的更新。
9.根据权利要求8所述的计算机可读存储介质,其特征在于,方法中所涉及的判定条件按下述方法得到:
以当前配准点为起点,以此刻磁传感器记录导出的代表器械朝向的方向向量为方向,建立一个向量即配准器械朝向向量,V为从气道最近点指向当前配准点的向量,所述气道最近点是指当前支气管树表面或支气管树各气道中心线中与配准点最近的点;不同状态下的判定条件由配准器械朝向向量/>和从气道最近点指向配准点的向量V之间的点积、余弦相似度和相关性,以及气道最近点和配准点之间的欧几里得距离中的一种或几种的组合与阈值相比较而得到。
10.根据权利要求9所述的计算机可读存储介质,其特征在于,所述的判定条件为:
建图预备判定条件为,在常规状态下,满足:
,
建图判定条件为,在建图预备状态下,满足:
,
建图预备结束判定条件为,在建图预备条件下,满足:
,
建图结束判定条件为,在建图条件下,满足:
,
其中,,/>为常数,/>的取值要大于/>的取值。
11.根据权利要求9所述的计算机可读存储介质,其特征在于,所述的判定条件为:
当气道最近点为支气管树各气道中心线上的点时:
建图预备判定条件为,在常规状态下,满足:
,
建图判定条件为,在建图预备状态下,满足:
,
建图预备结束判定条件为,在建图预备条件下,满足:
,
建图结束判定条件为,在建图条件下,满足:
,
其中,,/>为0与1之间的常数;/>,/>为常数,/>的取值要大于的取值;/>为根据当前气道最近点所在的气道分支级数对应的预计气道半径,/>为器械半径;
当气道最近点为支气管树表面上的点时:
建图预备判定条件为,在常规状态下,满足:
,
建图判定条件为,在建图预备状态下,满足:
,
建图预备结束判定条件为,在建图预备条件下,满足:
,
建图结束判定条件为,在建图条件下,满足:
,
其中,,/>为0与1之间的常数;/>,/>为常数,/>的取值要大于的取值;/>为根据当前气道最近点所在的气道分支级数对应的预计气道半径,/>为器械半径。
12.根据权利要求8所述的计算机可读存储介质,其特征在于,所述的气道点列表及其更新具体为:
在进入建图状态时,以保留气道最近点为起点,以代表器械朝向的方向向量为方向,建立一个向量即气道朝向向量,将预备配准点列表中的点向气道朝向向量上做投影,得到的投影点作为气道点列表,在建图状态中,实时记录的配准点继续向气道朝向向量上做投影,得到的投影点加入到气道点列表中。
13.根据权利要求8所述的计算机可读存储介质,其特征在于,每次更新气道点列表时,将气道点列表与地图信息进行融合得到新的地图信息,当地图信息是支气管树的表面信息,融合的方式如下:
利用气道点列表,进行平滑和插值,得到一个新的点列表;
依次连接上述新的点列表,得到一个线结构,将此线结构作为轴线,设置半径r,生成一个两端封闭的管道结构;
将此管道结构与支气管树的表面信息做布尔并运算,得到的结果即是气道点列表与地图信息融合的结果,即新的地图信息。
14.根据权利要求8所述的计算机可读存储介质,其特征在于,每次更新气道点列表时,将气道点列表与地图信息进行融合得到新的地图信息,如果地图信息是支气管树中各气道中心线信息,融合的方式如下:
利用气道点列表,进行平滑和插值,得到一个新的点列表;
依次连接上述新的点列表,得到一个线结构;
将此线结构加入到支气管树中各气道中心线信息中,得到的结果即是气道点列表与地图信息融合的结果,即新的地图信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311561950.7A CN117257459B (zh) | 2023-11-22 | 2023-11-22 | 一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311561950.7A CN117257459B (zh) | 2023-11-22 | 2023-11-22 | 一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117257459A CN117257459A (zh) | 2023-12-22 |
CN117257459B true CN117257459B (zh) | 2024-03-12 |
Family
ID=89212822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311561950.7A Active CN117257459B (zh) | 2023-11-22 | 2023-11-22 | 一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117257459B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101879092A (zh) * | 2010-06-30 | 2010-11-10 | 珠海和佳医疗设备股份有限公司 | 一种用于电磁定位手术导航系统的配准组件及其校准方法 |
CN111356414A (zh) * | 2017-12-05 | 2020-06-30 | 柯惠有限合伙公司 | 电磁导航与计算机断层摄影体积的多刚性配准 |
CN113288429A (zh) * | 2021-05-07 | 2021-08-24 | 山东商业职业技术学院 | 一种乳腺微创介入手术机器人空间配准与实时导航方法 |
JP2022039906A (ja) * | 2020-08-28 | 2022-03-10 | 中国計量大学 | マルチセンサによる複合キャリブレーション装置及び方法 |
WO2022123577A1 (en) * | 2020-12-10 | 2022-06-16 | Magnisity Ltd. | Dynamic deformation tracking for navigational bronchoscopy |
CN115120346A (zh) * | 2022-08-30 | 2022-09-30 | 中国科学院自动化研究所 | 目标点定位方法、装置、电子设备及支气管镜系统 |
CN116026318A (zh) * | 2022-07-20 | 2023-04-28 | 合肥科大智能机器人技术有限公司 | 基于多传感器融合slam和磁导的组合导航方法及存储介质 |
CN116019558A (zh) * | 2023-03-07 | 2023-04-28 | 内蒙古精劢医疗科技有限公司 | 电磁导航穿刺机器人系统及其定位方法 |
CN116999166A (zh) * | 2023-10-07 | 2023-11-07 | 杭州先奥科技有限公司 | 术中地图拓展方法、存储介质及电磁导航支气管镜系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10152786B2 (en) * | 2016-10-11 | 2018-12-11 | Biosense Webster (Israel) Ltd. | Registration of a magnetic tracking system with an imaging device |
CN106539624B (zh) * | 2016-11-23 | 2019-12-03 | 常州朗合医疗器械有限公司 | 医疗路径导航方法、规划方法及系统 |
US12004849B2 (en) * | 2017-12-11 | 2024-06-11 | Covidien Lp | Systems, methods, and computer-readable media for non-rigid registration of electromagnetic navigation space to CT volume |
-
2023
- 2023-11-22 CN CN202311561950.7A patent/CN117257459B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101879092A (zh) * | 2010-06-30 | 2010-11-10 | 珠海和佳医疗设备股份有限公司 | 一种用于电磁定位手术导航系统的配准组件及其校准方法 |
CN111356414A (zh) * | 2017-12-05 | 2020-06-30 | 柯惠有限合伙公司 | 电磁导航与计算机断层摄影体积的多刚性配准 |
JP2022039906A (ja) * | 2020-08-28 | 2022-03-10 | 中国計量大学 | マルチセンサによる複合キャリブレーション装置及び方法 |
WO2022123577A1 (en) * | 2020-12-10 | 2022-06-16 | Magnisity Ltd. | Dynamic deformation tracking for navigational bronchoscopy |
CN113288429A (zh) * | 2021-05-07 | 2021-08-24 | 山东商业职业技术学院 | 一种乳腺微创介入手术机器人空间配准与实时导航方法 |
CN116026318A (zh) * | 2022-07-20 | 2023-04-28 | 合肥科大智能机器人技术有限公司 | 基于多传感器融合slam和磁导的组合导航方法及存储介质 |
CN115120346A (zh) * | 2022-08-30 | 2022-09-30 | 中国科学院自动化研究所 | 目标点定位方法、装置、电子设备及支气管镜系统 |
CN116019558A (zh) * | 2023-03-07 | 2023-04-28 | 内蒙古精劢医疗科技有限公司 | 电磁导航穿刺机器人系统及其定位方法 |
CN116999166A (zh) * | 2023-10-07 | 2023-11-07 | 杭州先奥科技有限公司 | 术中地图拓展方法、存储介质及电磁导航支气管镜系统 |
Also Published As
Publication number | Publication date |
---|---|
CN117257459A (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11631174B2 (en) | Adaptive navigation technique for navigating a catheter through a body channel or cavity | |
CN110831481B (zh) | 管状网络的基于路径的导航 | |
US20210358142A1 (en) | Feature-based registration method | |
US10575907B2 (en) | Registration with trajectory information with shape sensing | |
JP4388958B2 (ja) | 内視鏡的パスプランニングのための方法及びシステム | |
CN116999166B (zh) | 术中地图拓展方法、存储介质及电磁导航支气管镜系统 | |
Soper et al. | In vivo validation of a hybrid tracking system for navigation of an ultrathin bronchoscope within peripheral airways | |
CN110432986B (zh) | 用于从ct数据构建虚拟径向超声图像的系统和方法 | |
Wang et al. | A visual SLAM-based bronchoscope tracking scheme for bronchoscopic navigation | |
CN114948199A (zh) | 一种外科手术辅助系统及手术路径规划方法 | |
CN114332013B (zh) | 基于肺动脉树分级的ct图像靶肺段识别方法 | |
CN108175502A (zh) | 一种支气管镜电磁导航系统 | |
JP2023552577A (ja) | ナビゲーション気管支鏡検査のための動的変形追跡 | |
CN117257459B (zh) | 一种抗呼吸干扰的电磁导航支气管镜术中地图拓展方法及系统 | |
WO2022098912A1 (en) | System and method for annotation of anatomical tree structures in 3d images | |
CN116327362A (zh) | 磁探头辅助支气管术中导航方法、装置、介质及电子设备 | |
CN114073581B (zh) | 一种支气管电磁导航系统 | |
Gil et al. | Intraoperative extraction of airways anatomy in videobronchoscopy | |
CN111281535B (zh) | 双腔支气管导管匹配和模拟插管方法及设备 | |
CN113616336A (zh) | 手术机器人仿真系统、仿真方法及可读存储介质 | |
CN116649875A (zh) | 基于表面重建技术的虚拟气管内窥镜导航系统及方法 | |
EP3454293B1 (en) | Method and apparatus for enhancement of bronchial airways representations using vascular morphology | |
US12008760B2 (en) | Systems and methods for estimating the movement of a target using a universal deformation model for anatomic tissue | |
CN117717413A (zh) | 一种用于支气管镜导航术的实时路径矫正方法 | |
CN116823894A (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 |