CN114451997A - 一种解决光学遮挡的手术导航装置及导航方法 - Google Patents

一种解决光学遮挡的手术导航装置及导航方法 Download PDF

Info

Publication number
CN114451997A
CN114451997A CN202210238950.2A CN202210238950A CN114451997A CN 114451997 A CN114451997 A CN 114451997A CN 202210238950 A CN202210238950 A CN 202210238950A CN 114451997 A CN114451997 A CN 114451997A
Authority
CN
China
Prior art keywords
coordinate system
endoscope
transformation matrix
tracking
camera
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
Application number
CN202210238950.2A
Other languages
English (en)
Other versions
CN114451997B (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.)
Changchun University of Science and Technology
Original Assignee
Changchun University of Science and Technology
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 Changchun University of Science and Technology filed Critical Changchun University of Science and Technology
Priority to CN202210238950.2A priority Critical patent/CN114451997B/zh
Publication of CN114451997A publication Critical patent/CN114451997A/zh
Application granted granted Critical
Publication of CN114451997B publication Critical patent/CN114451997B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • A61B2034/2046Tracking techniques
    • A61B2034/2055Optical tracking systems
    • A61B2034/2057Details of tracking cameras
    • 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
    • A61B2034/2046Tracking techniques
    • A61B2034/2059Mechanical position encoders

Abstract

一种解决光学遮挡的手术导航装置及导航方法属于医疗器械技术领域,目的在于解决现有技术存在的跟踪光线受遮挡,导航系统丢失内窥镜和手术工具的空间位置信息,无法进行手术操作的问题。本发明的一种解决光学遮挡的手术导航装置包括:光学跟踪设备;组合装置,组合装置包括被动式跟踪刚体、内窥镜以及相机,被动式跟踪刚体刚性固定于内窥镜尾端上方,相机刚性固定于内窥镜尾端下方;支架;支架尾端固定于手术部位附近的手术台上;以及棋盘标定板;棋盘标定板刚性固定于支架顶端。本发明在少量地增加了手术导航系统成本与复杂性的基础上,达到了保持跟踪信息的连续性、减少跟踪光线对医生的束缚、解决因跟踪光线遮挡问题带来的手术风险等效果。

Description

一种解决光学遮挡的手术导航装置及导航方法
技术领域
本发明属于医疗器械技术领域,具体涉及一种解决光学遮挡的手术导航装置以及导航方法。
背景技术
近年来,计算机软硬件技术、空间定位技术、医学影像技术的飞速发展,使得手术导航系统可以应用在更多类型的外科手术中。手术导航系统根据患者术前采集的影像数据(MRI、CT等)重建手术部位的三维模型,该模型在术前协助医生制定手术方案,在术中向医生反馈手术器械与病灶的相对位置。
在传统的外科手术中,病患受到的创伤较大、恢复时间较长。微创外科手术通过轻微的创伤,将手术工具或药物送入人体中,进行一系列的外科手术操作,从而达到治疗的效果。与传统外科手术相比,微创手术具有创伤小、恢复时间短的特点。因此,微创外科手术逐渐被普及到更多类型的外科手术中。但是,微创外科手术也大大地增加了手术难度。现今,应用于临床的手术导航系统大多是基于光学跟踪方法的,这种基于光学跟踪的手术导航系统具有较高的跟踪精度。但是,一旦跟踪光线受到了遮挡,导航系统就会丢失内窥镜和手术工具的空间位置信息。微创手术的操作空间非常狭小,在操作过程中遮挡是无法避免的。
发明内容
本发明的目的在于提出一种解决光学遮挡的手术导航装置以及导航方法,解决现有技术存在的跟踪光线受遮挡,导航系统丢失内窥镜和手术工具的空间位置信息,无法进行手术操作的问题。实现在光学跟踪光线被遮挡的情况下,实时的跟踪定位内窥镜和手术工具的空间位置信息。
为实现上述目的,本发明的一种解决光学遮挡的手术导航装置包括:
光学跟踪设备;
组合装置,所述组合装置包括被动式跟踪刚体、内窥镜以及相机,所述被动式跟踪刚体刚性固定于内窥镜尾端上方,相机刚性固定于内窥镜尾端下方;
支架;所述支架尾端固定于手术部位附近的手术台上;
以及棋盘标定板;所述棋盘标定板刚性固定于支架顶端。
基于一种解决光学遮挡的手术导航装置的导航方法包括以下步骤:
步骤一:在手术前,分别获取以下变换矩阵:
1)棋盘标定板坐标系至相机坐标系的变换矩阵,记为
Figure BDA0003536028220000021
并记录当前位姿下被动式跟踪刚体坐标系至世界坐标系的变换矩阵,记为
Figure BDA0003536028220000022
相机的内参矩阵,记为M;
2)被动式跟踪刚体坐标系至内窥镜坐标系的变换矩阵,记为
Figure BDA0003536028220000023
3)被动式跟踪刚体坐标系至相机坐标系的变换矩阵,记为
Figure BDA0003536028220000024
4)内窥镜坐标系至相机坐标系的变换矩阵,记为
Figure BDA0003536028220000025
5)棋盘标定板坐标系至世界坐标系的变换矩阵,记为
Figure BDA0003536028220000026
步骤二:在手术中,将光学跟踪设备捕获的内窥镜的位姿信息0与利用组合装置捕获的内窥镜的位姿信息B修正后的位姿信息B’输入扩展卡尔曼滤波器中,由扩展卡尔曼滤波器估计出较优的位姿信息,并以此作为内窥镜的跟踪结果;其中,当光学跟踪设备的跟踪光线没有被遮挡时,扩展卡尔曼滤波器主要根据位姿信息0与修正后的位姿信息B’来估计出较优的位姿信息;当跟踪光线被遮挡时,扩展卡尔曼滤波器根据修正后的组合装置捕获的内窥镜的位姿信息B′来估计出较优的位姿信息。
步骤一中第1)条获取的变换矩阵的方法具体为:
利用相机、被动式跟踪刚体、棋盘标定板、光学跟踪设备根据张正友标定法,获取相机的内参矩阵M与棋盘标定板坐标系至相机坐标系的变换矩阵
Figure BDA0003536028220000027
并通过光学跟踪设备记录当前位姿下被动式跟踪刚体坐标系至世界坐标系的变换矩阵,记为
Figure BDA0003536028220000028
步骤一中第2)条获取的变换矩阵的方法具体为:
利用内窥镜、被动式跟踪刚体、棋盘标定板、光学跟踪设备根据手眼标定法,获取被动式跟踪刚体坐标系至内窥镜坐标系的变换矩阵,记为
Figure BDA0003536028220000029
步骤一中第3)条获取的变换矩阵的方法具体为:
利用相机、被动式跟踪刚体、棋盘标定板、光学跟踪设备根据手眼标定法,获取被动式跟踪刚体坐标系至相机坐标系的变换矩阵
Figure BDA0003536028220000031
步骤一中第4)条获取的变换矩阵的方法具体为:
根据式(一)求得内窥镜坐标系至相机坐标系的变换矩阵
Figure BDA0003536028220000032
Figure BDA0003536028220000033
其中:
Figure BDA0003536028220000034
为术前获得的
Figure BDA0003536028220000035
的逆矩阵。
步骤一中第5)条获取的变换矩阵的方法具体为:
根据公式(二)求解出棋盘标定板坐标系至世界坐标系的变换矩阵
Figure BDA0003536028220000036
Figure BDA0003536028220000037
其中:
Figure BDA0003536028220000038
为术前获得的
Figure BDA0003536028220000039
的逆矩阵。
步骤二中所述的光学跟踪设备捕获的内窥镜的位姿信息0具体包括以下步骤:
步骤I、利用光学跟踪设备、被动式跟踪刚体根据光学跟踪法,获取当前被动式跟踪刚体至世界坐标系的变换矩阵,记为
Figure BDA00035360282200000310
步骤II、根据公式(三)计算获得当前位姿下内窥镜坐标系至世界坐标系的变换矩阵
Figure BDA00035360282200000311
Figure BDA00035360282200000312
其中:
Figure BDA00035360282200000313
为术前获得的
Figure BDA00035360282200000314
的逆矩阵;
变换矩阵
Figure BDA00035360282200000315
即为由光学跟踪设备获取的内窥镜在世界坐标系下的位姿信息,记为0。
步骤二中所述的组合装置捕获的内窥镜的位姿信息B具体包括以下步骤:
步骤I、从相机a3在当前位姿下获取的图像与棋盘标定板板坐标系中提取四对匹配点{Pa,Pb,Pc,Pd}与{PA,PB,PC,PD},PA与Pa、PB与Pb、PC与Pc、PD与Pd,一一对应;其中小写字母表示的点是大写字母表示的点在相机成像平面上的投影,根据这四对匹配点利用P3P方法即可求出点Pa、Pb、Pc在相机坐标系的坐标;其中,{PD,Pd}为验证点对,所以不需要求出Pd在相机坐标系的坐标。
步骤II、根据点Pa、Pb、Pc在相机a3坐标系的坐标利用ICP方法即可求出棋盘标定板坐标系到相机坐标系的变换矩阵
Figure BDA0003536028220000041
步骤III、根据公式(四)求得当前位姿下内窥镜坐标系至世界坐标系的变换矩阵
Figure BDA0003536028220000042
Figure BDA0003536028220000043
其中:
Figure BDA0003536028220000044
Figure BDA0003536028220000045
的逆矩阵;
变换矩阵
Figure BDA0003536028220000046
即为由组合装置获取的内窥镜在世界坐标系下的位姿信息,记为B。
步骤二中所述修正组合装置捕获的位姿信息B具体包括以下步骤:
步骤I、求取偏置误差E;
当光学跟踪设备的跟踪光线没有被遮挡时,利用公式(五)求得偏置误差E:
E=O-B (五)
当跟踪光线被遮挡时,利用公式(六)求得偏置误差E;其中,0w为跟踪光线没有被遮挡时光学跟踪设备最后捕获的内窥镜的位姿信息,Bw为跟踪光线没有被遮挡时组合装置最后捕获的内窥镜的位姿信息。
E=Ow-Bw (六)
步骤II、根据公式(七)修正位姿信息B;
B′=B+E (七)
从而得到修正后的位姿信息B′。
本发明的有益效果为:本发明在内窥镜末端固定一个相机,凭借该相机与一个已知世界坐标位置的棋盘标定板根据张正友标定法与手眼标定法获得内窥镜坐标系与棋盘标定板坐标系间的变换矩阵。因为棋盘标定板坐标系与世界坐标系间的变换矩阵已知,所以内窥镜坐标系与世界坐标系间的变换矩阵就可以被求出,从而能够在跟踪光线被遮挡的情况下,获取到内窥镜在世界坐标系下的位姿信息。以此来解决跟踪光线被遮挡时内窥镜的空间位姿信息丢失的问题。该方法保持了跟踪信息的连续性,减少了跟踪光线对医生的束缚,解决了因跟踪光线遮挡问题带来的手术风险。本发明在少量地增加了手术导航系统成本与复杂性的基础上,达到了保持跟踪信息的连续性、减少跟踪光线对医生的束缚、解决因跟踪光线遮挡问题带来的手术风险等效果。
附图说明
图1为本发明的一种解决光学遮挡的手术导航装置中组合装置以及棋盘标定板结构示意图;
图2为本发明的一种解决光学遮挡的手术导航装置整体结构示意图;
图3为术前求出的各个变换矩阵示意图;
图4为术中光学跟踪设备与组合装置捕获内窥镜位姿信息的示意图,
图5为棋盘标定板上的点投影到相机成像平面的示意图;
其中:a、组合装置,a1、被动式跟踪刚体,a2、内窥镜,a3、相机,b、棋盘标定板,c、支架,d、光学跟踪设备。
具体实施方式
下面结合附图对本发明的实施方式作进一步说明。
参见附图1至附图5,本发明的一种解决光学遮挡的手术导航装置包括:
光学跟踪设备d;光学跟踪设备是领域内比较常见的设备,我们在用的设备有多种型号,如NDI Polaris Vicra,0ptiTrack V120:Duo;
组合装置a,所述组合装置a包括被动式跟踪刚体a1、内窥镜a2以及相机a3,所述被动式跟踪刚体a1刚性固定于内窥镜a2尾端上方,相机a3刚性固定于内窥镜a2尾端下方;
支架c;所述支架c尾端固定于手术部位附近的手术台上;
以及棋盘标定板b;所述棋盘标定板b刚性固定于支架c顶端。棋盘标定板是在领域内一个具有公知性的工具。
基于一种解决光学遮挡的手术导航装置的导航方法包括以下步骤:
步骤一:在手术前,为了得到棋盘标定板坐标系至世界坐标系的变换矩阵,分别需要获取以下变换矩阵:利用张正友标定法、手眼标定发求取各个变换矩阵,在本领域内是公知的方法;
1)棋盘标定板b坐标系至相机a3坐标系的变换矩阵,记为
Figure BDA0003536028220000061
并记录当前位姿下被动式跟踪刚体a1坐标系至世界坐标系的变换矩阵,记为
Figure BDA0003536028220000062
相机a3的内参矩阵,记为M;
2)被动式跟踪刚体a1坐标系至内窥镜a2坐标系的变换矩阵,记为
Figure BDA0003536028220000063
3)被动式跟踪刚体a1坐标系至相机a3坐标系的变换矩阵,记为
Figure BDA0003536028220000064
4)内窥镜a2坐标系至相机a3坐标系的变换矩阵,记为
Figure BDA0003536028220000065
5)棋盘标定板b坐标系至世界坐标系的变换矩阵,记为
Figure BDA0003536028220000066
上述各个变换矩阵如图3所示,该步骤的目的是为了得到棋盘标定板b坐标系至世界坐标系的变换矩阵
Figure BDA0003536028220000067
世界坐标系由光学跟踪设备的制造商定义,其大多建立在光学跟踪设备的左相机处;
步骤二:在手术中,将光学跟踪设备d捕获的内窥镜a2的位姿信息0与利用组合装置a捕获的内窥镜a2的位姿信息B修正后的位姿信息B’输入扩展卡尔曼滤波器中,由扩展卡尔曼滤波器估计出较优的位姿信息,并以此作为内窥镜a2的跟踪结果;其中,当光学跟踪设备的跟踪光线没有被遮挡时,扩展卡尔曼滤波器主要根据位姿信息0与修正后的位姿信息B’来估计出较优的位姿信息;当跟踪光线被遮挡时,扩展卡尔曼滤波器根据修正后的组合装置a捕获的内窥镜a2的位姿信息B′来估计出较优的位姿信息。
图4中,实线为光学跟踪设备d捕获内窥镜a2位姿信息所用的变换矩阵,虚线为组合装置a捕获内窥镜a2位姿信息所用的变换矩阵;需要说明的是,加粗字体的变换矩阵为术前求得的,非加粗字体的变换矩阵为术中求得的;
步骤一中第1)条获取的变换矩阵的方法具体为:
利用相机a3获取支架c顶端的棋盘标定板b图像,并利用光学跟踪设备d获取当前位姿下内窥镜a2上方被动式跟踪刚体a1的坐标系与世界坐标系的变换关系,记为
Figure BDA0003536028220000068
重复上述操作N次,获取N张棋盘标定板图像与N个跟踪刚体坐标系至世界坐标系的变换矩阵,记为
Figure BDA0003536028220000069
根据N张棋盘标定板图像利用张正友标定法,求解出相机a3的内参矩阵与N个外参矩阵,并利用L-M算法优化求得的矩阵;外参矩阵,既棋盘标定板坐标系与相机a3的坐标系之间的变换矩阵,记为
Figure BDA0003536028220000071
其中,世界坐标系由光学跟踪设备的制造商定义,其大多建立在光学跟踪设备的左相机处;
步骤一中第2)条获取的变换矩阵的方法具体为:
利用内窥镜a2、被动式跟踪刚体a1、棋盘标定板b、光学跟踪设备d根据手眼标定法,获取被动式跟踪刚体a1坐标系至内窥镜a2坐标系的变换矩阵,记为
Figure BDA0003536028220000072
步骤一中第3)条获取的变换矩阵的方法具体为:
利用相机a3、被动式跟踪刚体a1、棋盘标定板b、光学跟踪设备d根据手眼标定法,获取被动式跟踪刚体a1坐标系至相机a3坐标系的变换矩阵
Figure BDA0003536028220000073
步骤一中第4)条获取的变换矩阵的方法具体为:
根据式(一)求得内窥镜a2坐标系至相机a3坐标系的变换矩阵
Figure BDA0003536028220000074
Figure BDA0003536028220000075
其中:
Figure BDA0003536028220000076
为术前获得的
Figure BDA0003536028220000077
的逆矩阵。
步骤一中第5)条获取的变换矩阵的方法具体为:
根据公式(二)求解出棋盘标定板b坐标系至世界坐标系的变换矩阵
Figure BDA0003536028220000078
Figure BDA0003536028220000079
其中:i={1,2,...,N},
Figure BDA00035360282200000710
为术前获得的
Figure BDA00035360282200000711
的逆矩阵。
根据公式(三)求解出N个
Figure BDA00035360282200000712
的均值
Figure BDA00035360282200000713
Figure BDA00035360282200000714
作为棋盘标定板坐标系至世界坐标系的变换矩阵。
Figure BDA00035360282200000715
步骤二中所述的光学跟踪设备d捕获的内窥镜a2的位姿信息0具体包括以下步骤:
步骤I、利用光学跟踪设备d、被动式跟踪刚体a1根据光学跟踪法,获取当前被动式跟踪刚体a1至世界坐标系的变换矩阵,记为
Figure BDA00035360282200000716
步骤II、根据公式(四)计算获得当前位姿下内窥镜a2坐标系至世界坐标系的变换矩阵
Figure BDA0003536028220000081
Figure BDA0003536028220000082
其中:
Figure BDA0003536028220000083
为术前获得的
Figure BDA0003536028220000084
的逆矩阵;
变换矩阵
Figure BDA0003536028220000085
即为内窥镜a2在世界坐标系下的位姿信息,记为0。
步骤二中所述的组合装置(a)捕获的内窥镜(a2)的位姿信息B具体包括以下步骤:
步骤I、从相机a3在当前位姿下获取的图像与棋盘标定板板坐标系中提取四对匹配点{Pa,Pb,Pc,Pd}与{PA,PB,PC,PD},PA与Pa、PB与Pb、PC与Pc、PD与Pd,一一对应;PA、PB和PC为棋盘标定板b上的点,Pa,Pb和Pc为相机a3成像平面上的点,其中小写字母表示的点是大写字母表示的点在相机a3成像平面上的投影,前三对匹配点的投影关系如图5所示,根据这四对匹配点利用P3P方法即可求出点Pa、Pb、Pc在相机a3坐标系的坐标;其中,{PD,Pd}为验证点对,所以不需要求出Pd在相机a3坐标系的坐标。
步骤II、根据点Pa、Pb、Pc在相机a3坐标系的坐标利用ICP方法即可求出棋盘标定板坐标系到相机a3坐标系的变换矩阵
Figure BDA0003536028220000086
步骤III、根据公式(五)求得当前位姿下内窥镜a2坐标系至世界坐标系的变换矩阵
Figure BDA0003536028220000087
Figure BDA0003536028220000088
其中:
Figure BDA0003536028220000089
Figure BDA00035360282200000810
的逆矩阵;
变换矩阵
Figure BDA00035360282200000811
即为内窥镜a2在世界坐标系下的位姿信息,记为B;
步骤二中所述修正组合装置a捕获的位姿信息B具体包括以下步骤:
步骤I、求取偏置误差E;
当光学跟踪设备的跟踪光线没有被遮挡时,利用公式(六)求得偏置误差E;
E=O-B (六)
当跟踪光线被遮挡时,利用公式(七)求得偏置误差E;其中,0w为跟踪光线没有被遮挡时光学跟踪设备d最后捕获的内窥镜a2的位姿信息,Bw为跟踪光线没有被遮挡时组合装置a最后捕获的内窥镜a2的位姿信息;
E=Ow-Bw (七)
需要说明的是,光学跟踪设备d拥有判断是否跟踪到被动式跟踪刚体a1的功能,当光学跟踪设备d判断追踪到了被动式跟踪刚体a1时,说明光学跟踪设备的跟踪光线没有被遮挡;当光学跟踪设备d判断没有追踪到被动式跟踪刚体a1时,说明光学跟踪设备的跟踪光线被遮挡;
步骤II、根据公式(八)修正位姿信息B;
B′=B+E (八)
从而得到修正后的位姿信息B′。

Claims (10)

1.一种解决光学遮挡的手术导航装置,其特征在于,包括:
光学跟踪设备(d);
组合装置(a),所述组合装置(a)包括被动式跟踪刚体(a1)、内窥镜(a2)以及相机(a3),所述被动式跟踪刚体(a1)刚性固定于内窥镜(a2)尾端上方,相机(a3)刚性固定于内窥镜(a2)尾端下方;
支架(c);所述支架(c)尾端固定于手术部位附近的手术台上;
以及棋盘标定板(b);所述棋盘标定板(b)刚性固定于支架(c)顶端。
2.基于权利要求1所述的一种解决光学遮挡的手术导航装置的导航方法,其特征在于,包括以下步骤:
步骤一:在手术前,分别需要获取以下变换矩阵:
1)棋盘标定板(b)坐标系至相机(a3)坐标系的变换矩阵,记为
Figure FDA0003536028210000011
并记录当前位姿下被动式跟踪刚体(a1)坐标系至世界坐标系的变换矩阵,记为
Figure FDA0003536028210000012
相机(a3)的内参矩阵,记为M;
2)被动式跟踪刚体(a1)坐标系至内窥镜(a2)坐标系的变换矩阵,记为
Figure FDA0003536028210000013
3)被动式跟踪刚体(a1)坐标系至相机(a3)坐标系的变换矩阵,记为
Figure FDA0003536028210000014
4)内窥镜(a2)坐标系至相机(a3)坐标系的变换矩阵,记为
Figure FDA0003536028210000015
5)棋盘标定板(b)坐标系至世界坐标系的变换矩阵,记为
Figure FDA0003536028210000016
步骤二:在手术中,将光学跟踪设备(d)捕获的内窥镜(a2)的位姿信息O与利用组合装置(a)捕获的内窥镜(a2)的位姿信息B修正后的位姿信息B’输入扩展卡尔曼滤波器中,由扩展卡尔曼滤波器估计出较优的位姿信息,并以此作为内窥镜(a2)的跟踪结果;其中,当光学跟踪设备(d)的跟踪光线没有被遮挡时,扩展卡尔曼滤波器主要根据位姿信息O与修正后的位姿信息B’来估计出较优的位姿信息;当跟踪光线被遮挡时,扩展卡尔曼滤波器根据修正后的组合装置(a)捕获的内窥镜(a2)的位姿信息B'来估计出较优的位姿信息。
3.根据权利要求2所述导航方法,其特征在于,步骤一中第1)条获取的变换矩阵的方法具体为:
利用相机(a3)、被动式跟踪刚体(a1)、棋盘标定板(b)、光学跟踪设备(d)根据张正友标定法,获取相机(a3)的内参矩阵M与棋盘标定板(b)坐标系至相机(a3)坐标系的变换矩阵
Figure FDA0003536028210000021
并通过光学跟踪设备(d)记录当前位姿下被动式跟踪刚体(a1)坐标系至世界坐标系的变换矩阵,记为
Figure FDA0003536028210000022
4.根据权利要求2所述导航方法,其特征在于,步骤一中第2)条获取的变换矩阵的方法具体为:
利用内窥镜(a2)、被动式跟踪刚体(a1)、棋盘标定板(b)、光学跟踪设备(d)根据手眼标定法,获取被动式跟踪刚体(a1)坐标系至内窥镜(a2)坐标系的变换矩阵,记为
Figure FDA0003536028210000023
5.根据权利要求2所述导航方法,其特征在于,步骤一中第3)条获取的变换矩阵的方法具体为:
利用相机(a3)、被动式跟踪刚体(a1)、棋盘标定板(b)、光学跟踪设备(d)根据手眼标定法,获取被动式跟踪刚体(a1)坐标系至相机(a3)坐标系的变换矩阵
Figure FDA0003536028210000024
6.根据权利要求2所述导航方法,其特征在于,步骤一中第4)条获取的变换矩阵的方法具体为:
根据式(一)求得内窥镜(a2)坐标系至相机(a3)坐标系的变换矩阵
Figure FDA0003536028210000025
Figure FDA0003536028210000026
其中:
Figure FDA0003536028210000027
为术前获得的
Figure FDA0003536028210000028
的逆矩阵。
7.根据权利要求2所述导航方法,其特征在于,步骤一中第5)条获取的变换矩阵的方法具体为
根据公式(二)求解出棋盘标定板(b)坐标系至世界坐标系的变换矩阵
Figure FDA0003536028210000029
Figure FDA00035360282100000210
其中:
Figure FDA0003536028210000031
为术前获得的
Figure FDA0003536028210000032
的逆矩阵。
8.根据权利要求2所述导航方法,其特征在于,步骤二中所述的光学跟踪设备(d)捕获的内窥镜(a2)的位姿信息O具体包括以下步骤:
步骤Ⅰ、利用光学跟踪设备(d)、被动式跟踪刚体(a1)根据光学跟踪法,获取当前被动式跟踪刚体(a1)至世界坐标系的变换矩阵,记为
Figure FDA0003536028210000033
步骤Ⅱ、根据公式(三)计算获得当前位姿下内窥镜(a2)坐标系至世界坐标系的变换矩阵
Figure FDA0003536028210000034
Figure FDA0003536028210000035
其中:
Figure FDA0003536028210000036
为术前获得的
Figure FDA0003536028210000037
的逆矩阵;
变换矩阵
Figure FDA0003536028210000038
即为由光学跟踪设备(d)获取的内窥镜(a2)在世界坐标系下的位姿信息,记为O。
9.根据权利要求2所述导航方法,其特征在于,步骤二中所述的组合装置(a)捕获的内窥镜(a2)的位姿信息B具体包括以下步骤:
步骤Ⅰ、从相机(a3)在当前位姿下获取的图像与棋盘标定板(b)板坐标系中提取四对匹配点{Pa,Pb,Pc,Pd}与{PA,PB,PC,PD},PA与Pa、PB与Pb、PC与Pc、PD与Pd,一一对应;其中小写字母表示的点是大写字母表示的点在相机a3成像平面上的投影,根据这四对匹配点利用P3P方法即可求出点Pa、Pb、Pc在相机a3坐标系的坐标;其中,{PD,Pd}为验证点对,所以不需要求出Pd在相机a3坐标系的坐标;
步骤Ⅱ、根据点Pa、Pb、Pc在相机a3坐标系的坐标利用ICP方法即可求出棋盘标定板坐标系到相机a3坐标系的变换矩阵
Figure FDA0003536028210000039
步骤Ⅲ、根据公式(四)求得当前位姿下内窥镜(a2)坐标系至世界坐标系的变换矩阵
Figure FDA00035360282100000310
Figure FDA00035360282100000311
其中:
Figure FDA0003536028210000041
Figure FDA0003536028210000042
的逆矩阵;
变换矩阵
Figure FDA0003536028210000043
即为由组合装置(a)获取的内窥镜(a2)在世界坐标系下的位姿信息,记为B。
10.根据权利要求9所述导航方法,其特征在于,所述的修正后的位姿信息B′具体包括以下步骤:
步骤I、求取偏置误差E;
当光学跟踪设备的跟踪光线没有被遮挡时,利用公式(五)求得偏置误差E;
E=O-B (五)
当跟踪光线被遮挡时,利用公式(六)求得偏置误差E;其中,Ow为跟踪光线没有被遮挡时光学跟踪设备(d)最后捕获的内窥镜(a2)的位姿信息,Bw为跟踪光线没有被遮挡时组合装置(a)最后捕获的内窥镜(a2)的位姿信息;
E=Ow-Bw (六)
步骤II、根据公式(七)修正位姿信息B;
B′=B+E (七)
从而得到修正后的位姿信息B′。
CN202210238950.2A 2022-03-08 2022-03-08 一种解决光学遮挡的手术导航装置及导航方法 Active CN114451997B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210238950.2A CN114451997B (zh) 2022-03-08 2022-03-08 一种解决光学遮挡的手术导航装置及导航方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210238950.2A CN114451997B (zh) 2022-03-08 2022-03-08 一种解决光学遮挡的手术导航装置及导航方法

Publications (2)

Publication Number Publication Date
CN114451997A true CN114451997A (zh) 2022-05-10
CN114451997B CN114451997B (zh) 2023-11-28

Family

ID=81416511

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210238950.2A Active CN114451997B (zh) 2022-03-08 2022-03-08 一种解决光学遮挡的手术导航装置及导航方法

Country Status (1)

Country Link
CN (1) CN114451997B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080097156A1 (en) * 2006-10-23 2008-04-24 Pentax Corporation Camera calibration for endoscope navigation system
US20140049629A1 (en) * 2011-04-29 2014-02-20 The Johns Hopkins University Sytem and method for tracking and navigation
US20160000516A1 (en) * 2014-06-09 2016-01-07 The Johns Hopkins University Virtual rigid body optical tracking system and method
US20180049622A1 (en) * 2016-08-16 2018-02-22 Insight Medical Systems, Inc. Systems and methods for sensory augmentation in medical procedures
US20190000564A1 (en) * 2015-12-30 2019-01-03 The Johns Hopkins University System and method for medical imaging
CN109674534A (zh) * 2017-10-18 2019-04-26 深圳市掌网科技股份有限公司 一种基于增强现实的手术导航图像显示方法和系统
CN111590594A (zh) * 2020-06-22 2020-08-28 南京航空航天大学 基于视觉引导的机器人轨迹跟踪控制方法
CN112991457A (zh) * 2021-02-22 2021-06-18 北京理工大学 手术导航中投影仪空间位置和内外参标定方法、装置
CN113940755A (zh) * 2021-09-30 2022-01-18 南开大学 一种术像一体的外科手术规划与导航方法
CN114027980A (zh) * 2021-10-30 2022-02-11 浙江德尚韵兴医疗科技有限公司 一种介入手术机器人系统及其标定与误差补偿方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080097156A1 (en) * 2006-10-23 2008-04-24 Pentax Corporation Camera calibration for endoscope navigation system
US20140049629A1 (en) * 2011-04-29 2014-02-20 The Johns Hopkins University Sytem and method for tracking and navigation
US20160000516A1 (en) * 2014-06-09 2016-01-07 The Johns Hopkins University Virtual rigid body optical tracking system and method
US20190000564A1 (en) * 2015-12-30 2019-01-03 The Johns Hopkins University System and method for medical imaging
US20180049622A1 (en) * 2016-08-16 2018-02-22 Insight Medical Systems, Inc. Systems and methods for sensory augmentation in medical procedures
CN109674534A (zh) * 2017-10-18 2019-04-26 深圳市掌网科技股份有限公司 一种基于增强现实的手术导航图像显示方法和系统
CN111590594A (zh) * 2020-06-22 2020-08-28 南京航空航天大学 基于视觉引导的机器人轨迹跟踪控制方法
CN112991457A (zh) * 2021-02-22 2021-06-18 北京理工大学 手术导航中投影仪空间位置和内外参标定方法、装置
CN113940755A (zh) * 2021-09-30 2022-01-18 南开大学 一种术像一体的外科手术规划与导航方法
CN114027980A (zh) * 2021-10-30 2022-02-11 浙江德尚韵兴医疗科技有限公司 一种介入手术机器人系统及其标定与误差补偿方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
WEI HE, KUMSOK KANG, HUAMIN YANG: "A hand-eye calibration method for computer assisted endoscopy", 11TH INTERNATIONAL CONFERENCE ON NATURAL COMPUTATION *
易新华;钱晋武;张震;章亚男;沈林勇;: "基于视觉的内窥镜端部位姿模型分析", 高技术通讯, no. 12 *
李维维: "基于CBCT 影像与双目视觉技术的颞下颌关节运动可视化系统开发与实现", 中国优秀硕士学位论文全文数据库 *
李超;吴玉敬;徐飞飞;赵创社;张慧;孟立新;: "一种目标遮挡情况下的自动跟踪控制方法", 应用光学, no. 05 *
贺长宇;刘越;王涌天;: "用于增强现实手术导航系统的光学-惯性混合跟踪方法", 计算机辅助设计与图形学学报, no. 03 *

Also Published As

Publication number Publication date
CN114451997B (zh) 2023-11-28

Similar Documents

Publication Publication Date Title
CN109416841B (zh) 影像增强真实度的方法与应用该方法在可穿戴式眼镜的手术导引
EP2967297B1 (en) System for dynamic validation, correction of registration for surgical navigation
US20220054200A1 (en) Calibration method and device for dental implant navigation surgery, and tracking method and device for dental implant navigation surgery
US6782287B2 (en) Method and apparatus for tracking a medical instrument based on image registration
Shahidi et al. Implementation, calibration and accuracy testing of an image-enhanced endoscopy system
US10166078B2 (en) System and method for mapping navigation space to patient space in a medical procedure
US20150366628A1 (en) Augmented surgical reality environment system
JP2002186603A (ja) 対象物の案内のための座標変換法
WO2017117517A1 (en) System and method for medical imaging
CN111281540B (zh) 基于虚实融合的骨科微创术中实时可视化导航系统
CN110464462B (zh) 腹部外科介入手术的图像导航配准系统及相关装置
WO2020129034A1 (en) Robotic surgery systems and surgical guidance methods thereof
WO2001057805A2 (en) Image data processing method and apparatus
CN114451997B (zh) 一种解决光学遮挡的手术导航装置及导航方法
CN112381750A (zh) 一种超声波图像和ct/mri图像的多模态配准融合方法
Harders et al. Multimodal augmented reality in medicine
CN113317874B (zh) 一种医学图像处理装置及介质
CN114283179A (zh) 基于超声图像的骨折远近端空间位姿实时获取与配准系统
WO2016042297A1 (en) Computer and computer-implemented method for supporting laparoscopic surgery
Konen et al. An image-based navigation support system for neuroendoscopic surgery
CN113925611A (zh) 对象三维模型和对象实体的匹配方法、装置、设备及介质
CN114073579A (zh) 手术导航方法、装置、电子设备及存储介质
TWI741889B (zh) 手術空間註冊系統與方法
US11832895B2 (en) Method and system for register operating space
Westermann et al. Online head motion tracking applied to the patient registration problem

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