CN106560163B - 手术导航系统及手术导航系统的配准方法 - Google Patents
手术导航系统及手术导航系统的配准方法 Download PDFInfo
- Publication number
- CN106560163B CN106560163B CN201510650347.5A CN201510650347A CN106560163B CN 106560163 B CN106560163 B CN 106560163B CN 201510650347 A CN201510650347 A CN 201510650347A CN 106560163 B CN106560163 B CN 106560163B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- space coordinate
- scanning
- point
- 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.)
- Active
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种手术导航系统及一种手术导航系统的配准方法,所述手术导航系统包括动态基准、扫描装置、立体相机及处理装置。所述扫描装置用于在术前扫描所述手术部位以获取三维扫描图像。所述立体相机用于在术中获取所述手术部位实时图像。所述处理装置用于,处理所述三维扫描图像及所述实时图像,以计算所述实时图像上的任一点在所述三维扫描图像的特征匹配点。本发明实施方式的手术导航系统在配准过程中能够减少人工操作的误差,即使在病患部位的面积很小的情况下也能快速实现自动配准,缩短了手术时间,进而降低手术风险,符合微创手术的需求。
Description
技术领域
本发明涉及医疗领域,尤其涉及一种手术导航系统及一种手术导航系统的配准方法。
背景技术
在手术过程中,往往需要对实际手术部位与术前扫描获取手术部位的三维图像进行配准以引导手术进行。目前的手术导航系统的配准方法中,医生要么人工选取并确认配准过程中的特征匹配点,从而实现配准,但是,由于医生的操作误差,导致配准精度低;要么在手术中大面积暴露患者的手术部位以实现自动配准,但是,不符合微创手术的需求。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此,本发明提供一种手术导航系统及一种手术导航系统的配准方法。
本发明实施方式的手术导航系统包括扫描装置、立体相机、处理装置及动态基准。所述扫描装置用于在术前手术部位以获取所述扫描图像。所述扫描图像包括扫描图像坐标系。所述立体相机用于在术中获取所述手术部位的实时图像。所述立体相机包括相机空间坐标系。所述处理装置连接所述扫描装置、所述立体相机及动态基准。所述处理装置用于处理所述扫描图像及所述实时图像,以对所述相机空间坐标系与所述扫描图像坐标系进行坐标系转换,及计算所述实时图像上的任一点在所述扫描图像的特征匹配点。
在本发明实施方式的手术导航系统中,处理装置预存有所述实时图像上的任一点在所述扫描图像的特征匹配点的计算算法,在手术进行时,手术部位的实时图像上的任一点都能够在手术部位的扫描图像中找到特征匹配点,从而实现术中手术部位的实际位置与术前的扫描图像之间的配准,进而准确地引导手术的进行。本发明实施方式的手术导航系统在配准过程中能够减少人工操作的误差,即使在手术部位的面积很小的情况下也能够快速实现自动配准,缩短了手术时间,进而降低手术风险,符合微创手术的需求。在某些实施方式中,所述扫描装置包括扫描空间坐标系,所述动态基准包括动态空间坐标系;
所述处理装置还用于,对所述扫描图像坐标系、所述相机空间坐标系、所述扫描空间坐标系及所述动态空间坐标系进行坐标系转换,进而计算所述相机空间坐标系中的任一点在所述扫描图像坐标系中的特征匹配点,以得到所述实时图像上的任一点在所述扫描图像的特征匹配点。
在某些实施方式中,所述扫描装置上设有至少三个不共线的第一标志,所述第一标志用于建立所述扫描空间坐标系;
所述动态基准上设有至少三个不共线的第二标志,所述第二标志用于建立所述动态空间坐标系;
所述立体相机用于根据双目成像原理,通过追踪所述第一标志及所述第二标志发出的信号确认所述相机空间坐标系与所述扫描空间坐标系及所述动态空间坐标系之间的坐标系转换关系。
在某些实施方式中,所述处理装置还用于,计算所述扫描空间坐标系中任一点在所述扫描图像坐标系中的特征匹配点,计算所述相机空间坐标系中任一点在所述扫描空间坐标系及所述动态空间坐标系中的特征匹配点,计算所述动态空间坐标系中任一点在所述扫描空间坐标系中的特征匹配点,计算所述动态空间坐标系中任一点在所述扫描图像坐标系中的特征匹配点,及计算任意时刻所述动态空间坐标系中任一点在所述相机空间坐标系中的特征匹配点。在某些实施方式中,所述相机空间坐标系中的任一点在所述扫描图像坐标系中的特征匹配点根据以下公式计算:
式中:
为所述扫描图像坐标系中的任一点的坐标;
为所述相机空间坐标系中任一点的坐标;
R为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的旋转矩阵;
R12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的旋转矩阵;
为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的旋转矩阵的逆矩阵;
R13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的旋转矩阵;
T13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的平移向量;
T13为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的平移向量;
T12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的平移向量;
T为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的平移向量。
在某些实施方式中,所述手术导航系统还包括导航探针,所述导航探针包括探头,所述探头用于在所述手术部位上点取比较点,所述导航探针包括探针空间坐标系,所述处理装置用于,对所述探针空间坐标系与所述相机空间坐标系进行坐标系转换,进而计算获取所述比较点在所述扫描图像坐标系中的特征匹配点,以实现所述手术导航系统术中导航。
在某些实施方式中,所述比较点在所述扫描图像坐标系中的特征匹配点根据以下公式计算:
式中:
为所述扫描图像坐标系中的任一点的坐标;
为所述比较点在所述探针空间坐标系中的坐标;
R为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的旋转矩阵;
R12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的旋转矩阵;
为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的旋转矩阵的逆矩阵;
R13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的旋转矩阵;
为任意时刻所述探针空间坐标系与所述相机空间坐标系之间的旋转矩阵的逆矩阵;
T14(t)为任意时刻所述探针空间坐标系与所述相机空间坐标系之间的平移向量;
T13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的平移向量;
T13为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的平移向量;
T12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的平移向量;
T为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的平移向量。
本发明实施方式的手术导航系统的配准方法,包括以下步骤:
在受施手术者的手术部位固定设置动态基准;
在术前获取所述手术部位的扫描图像;
对扫描图像坐标系与相机空间坐标系进行坐标系转换;
在术中获取所述手术部位的实时图像;及
计算所述实时图像上的任一点在所述扫描图像的特征匹配点。
本发明实施方式的手术导航系统的配准方法中,在手术进行时,手术部位的实时图像上的任一点都能够在手术部位的扫描图像中找到特征匹配点,从而实现术中手术部位的实际位置与术前的扫描图像之间的配准,进而准确地引导手术的进行。本发明实施方式的手术导航系统的配准方法在配准过程中能够减少人工操作的误差,即使在手术部位的面积很小的情况下也能够快速实现自动配准,缩短了手术时间,进而降低手术风险,符合微创手术的需求。
在某些实施方式中,所述对扫描图像坐标系与相机空间坐标系进行坐标系转换的步骤包括:
对所述扫描图像坐标系、所述相机空间坐标系、扫描空间坐标系及动态空间坐标系进行坐标系转换;及
计算所述相机空间坐标系中的任一点在所述扫描图像坐标系中的特征匹配点。
在某些实施方式中,所述对所述扫描图像坐标系、所述相机空间坐标系、扫描空间坐标系及动态空间坐标系进行坐标系转换的步骤包括:
计算所述扫描空间坐标系中任一点在所述扫描图像坐标系中的特征匹配点;
计算所述相机空间坐标系中任一点在所述扫描空间坐标系及所述动态空间坐标系中的特征匹配点;
计算所述动态空间坐标系中任一点在所述扫描空间坐标系中的特征匹配点;
计算所述动态空间坐标系中任一点在所述扫描图像坐标系中的特征匹配点;及
计算任意时刻所述动态空间坐标系中任一点在所述相机空间坐标系中的特征匹配点。
在某些实施方式中,所述相机空间坐标系中的任一点在所述扫描图像坐标系中的特征匹配点根据以下公式计算:
式中:
为所述扫描图像坐标系中的任一点的坐标;
为所述相机空间坐标系中任一点的坐标;
R为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的旋转矩阵;
R12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的旋转矩阵;
为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的旋转矩阵的逆矩阵;
R13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的旋转矩阵;
T13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的平移向量;
T13为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的平移向量;
T12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的平移向量;
T为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的平移向量。
在某些实施方式中,所述手术导航系统的配准方法还包括步骤:
在所述手术部位上点取比较点;
对探针空间坐标系与所述相机空间坐标系进行坐标系转换;及
计算获取所述比较点在所述扫描图像坐标系中的特征匹配点。
在某些实施方式中,所述比较点在所述扫描图像坐标系中的特征匹配点根据以下公式计算:
为所述扫描图像坐标系中的任一点的坐标;
为所述比较点在所述探针空间坐标系中的坐标;
R为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的旋转矩阵;
R12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的旋转矩阵;
为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的旋转矩阵的逆矩阵;
R13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的旋转矩阵;
为任意时刻所述探针空间坐标系与所述相机空间坐标系之间的旋转矩阵的逆矩阵;
T14(t)为任意时刻所述探针空间坐标系与所述相机空间坐标系之间的平移向量;
T13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的平移向量;
T13为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的平移向量;
T12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的平移向量;
T为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的平移向量。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施方式的描述中将变得明显和容易理解,其中:
图1是本发明实施方式的手术导航系统的模块示意图。
图2是本发明实施方式的手术导航系统的原理示意图。
图3是本发明实施方式的手术导航系统的配准方法的流程示意图。
图4是本发明实施方式的手术导航系统的配准方法的步骤3的流程示意图。
图5是本发明实施方式的手术导航系统的配准方法的步骤31的流程示意图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
请参阅图1及图2,本发明实施方式的手术导航系统100包括扫描装置10、立体相机20、处理装置30及动态基准40。扫描装置10用于在术前扫描患者的手术部位以获取扫描图像,扫描图像包括扫描图像坐标系O。立体相机20用于在术中获取手术部位的实时图像,立体相机20包括相机空间坐标系O1。处理装置30连接扫描装置10、立体相机20及动态基准40。处理装置30用于处理扫描图像及实时图像,以对相机空间坐标系O1与扫描图像坐标系O进行坐标系转换,及计算实时图像上的任一点在扫描图像的特征匹配点。动态基准40用于固定设置在受施手术者的手术部位上。
在本发明实施方式的手术导航系统100中,处理装置30预存有计算实时图像上的任一点在扫描图像的特征匹配点的算法,在手术进行时,手术部位的实时图像上的任一点都能够在手术部位的扫描图像中找到特征匹配点,从而实现术中手术部位的实际位置与术前的扫描图像之间的配准,进而准确地引导手术的进行。本发明实施方式的手术导航系统100在配准过程中能够减少人工操作的误差,即使在手术部位的面积很小的情况下也能够快速实现自动配准,缩短了手术时间,进而降低手术风险,符合微创手术的需求。
具体地,扫描装置10可以为计算机断层扫描(Computed Tomography,CT)设备或核磁共振成像(Magnetic Resonance Imaging,MRI)设备等人体扫描设备。
在本实施方式中,扫描装置10包括扫描空间坐标系O2,动态基准40包括动态空间坐标系O3。处理装置30还用于,对扫描图像坐标系O、相机空间坐标系O1、扫描空间坐标系O2及动态空间坐标系O3进行坐标系转换,进而计算相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点,以得到实时图像上的任一点在扫描图像的特征匹配点。
本实施方式中,扫描装置10上设有至少三个不共线的第一标志11,第一标志11用于建立扫描空间坐标系O2。
动态基准40上设有至少三个不共线的第二标志41,第二标志41用于建立动态空间坐标系O3。
立体相机20用于根据双目成像原理,通过追踪所述第一标志41及所述第二标志42发出的信号确定相机空间坐标系O1与扫描空间坐标系O2及动态空间坐标系O3之间的坐标系转换关系。
具体地,双目成像原理类似人眼的三维感知机制,通过不同视角来观察场景,而不同视角观察场景存在一定的视觉差,根据这一视差就可以得到场景中物体距离相机的距离,也即场景的三维结构信息。
在本实施方式中,立体相机20提前完成了标定以及校正,得到了相机的焦距、主点以及相机之间基线长度等信息。
第一标志11、第二标志41为发光二极管,例如为近红外发光二极管,如此,第一标志11及第二标志41发出不可见光,从而避免对人眼产生刺激。
第一标志11、第二标志41发出的光经过滤光片滤过杂光在立体相机20传感器上得到清晰的图像,然后通过图像处理以及极线校正得到各个发光二极管在立体相机20上的视差,最后利用前面立体相机20标定得到的参数采用三角测量原理计算出各个发光二极管在相机空间坐标系O1中的坐标,进而确定相机空间坐标系O1与扫描空间坐标系O2及动态空间坐标系O3之间的坐标系转换关系。
在手术前及手术中,将第一标志11及第二标志41设置于立体相机20的视野范围内,并且在扫描装置10扫描手术部位的时候保持第一标志11、第二标志41及立体相机20的位置不变,以便立体相机20追踪第一标志11及第二标志41发出的信息,从而保证扫描图像坐标系O、相机空间坐标系O1、扫描空间坐标系O2及动态空间坐标系O3的转换关系唯一。
如此,通过坐标系的转换计算,可准确地得到手术部位的实时图像上的任一点在扫描图像的特征匹配点。
在本实施方式中,相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点根据以下公式(一)计算:
(一)
式中:
为扫描图像坐标系O中的任一点的坐标;
为相机空间坐标系O1中任一点的坐标;
R为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的旋转矩阵;
R12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的旋转矩阵;
为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的旋转矩阵的逆矩阵;
R13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的旋转矩阵;
T13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的平移向量;
T13为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的平移向量;
T12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的平移向量;
T为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的平移向量。
如此,可准确地计算相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点。
本实施方式中,处理装置30还用于,计算扫描空间坐标系O2中任一点在扫描图像坐标系O中的特征匹配点,计算相机空间坐标系O1中任一点在扫描空间坐标系O2及动态空间坐标系O3中的特征匹配点,计算动态空间坐标系O3中任一点在扫描空间坐标系O2中的特征匹配点,计算动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点,及计算任意时刻动态空间坐标系O3中任一点在相机空间坐标系O1中的特征匹配点。
如此,准确确认扫描图像坐标系O、相机空间坐标系O1、扫描空间坐标系O2及动态空间坐标系O3之间的转换关系。
本实施方式中,具体地,扫描空间坐标系O2中任一点在扫描图像坐标系O中的特征匹配点通过以下公式(二)计算:
(二)
式中:
为扫描空间坐标系O2中的任一点的坐标。
其中,R和T可根据扫描装置10的成像原理计算得到,如此,确定扫描空间坐标系O2与扫描图像坐标系O之间的转换关系。
相机空间坐标系O1中任一点在扫描空间坐标系O2及动态空间坐标系O3中的特征匹配点可通过以下公式(三)计算:
(三)
式中:
为相机空间坐标系O1中任一点的坐标。
其中,R12和T12可根据立体相机20的双目成像原理计算而得。如此,可确定相机空间坐标系O1与扫描空间坐标系O2之间的转换关系。
相机空间坐标系O1中任一点在动态空间坐标系O3中的特征匹配点可通过以下公式(四)计算:
(四)
为动态空间坐标系O3中的任一点的坐标。
其中,R13和T13可根据立体相机20的双目成像原理计算而得。如此,确定相机空间坐标系O1与扫描空间坐标系O2之间的转换关系。
进一步地,根据上述公式(三)及(四),动态空间坐标系O3中任一点在扫描空间坐标系O2中的特征匹配点通过以下公式计算:
(五)
进一步地,根据公式(二)及(五),动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点通过以下公式计算:
(六)
在扫描装置10扫描受施手术者的手术部位以获取扫描图像的过程中,由于第一标志11、第二标志41及立体相机20的位置保持不变,因此,扫描图像坐标系O、相机空间坐标系O1、扫描空间坐标系O2及动态空间坐标系O3之间的转换关系唯一,即旋转矩阵R、R12、R13,以及平移向量T、T12、T13恒为常量,也就是说,手术部位在动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点的对应关系确认唯一。
进一步地,任意时刻动态空间坐标系O3中任一点在相机空间坐标系O1中的特征匹配点通过以下公式计算:
(七)
具体地,在手术过程中,动态基准40与立体相机20的相对位置可能发生移动,但是立体相机20通过追踪动态基准40上的第二标志41发出的信息可确定任意时刻相机空间坐标系O1与动态空间坐标系O3的旋转矩阵R13(t)及平移向量T13(t),如此,在手术过程中,当动态基准40与立体相机20的相对位置发生改变时,处理装置30也能快速地确认任意时刻手术部位在动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点。
如此,根据上述公式(六)及公式(七)的计算得到公式(一),即可确定相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点且关系唯一,也就是说,立体相机20获取手术部位的实时图像上的任一点都能够在扫描图像坐标系O中确认唯一的特征匹配点,即手术部位的实时图像上的任一点在扫描图像中的特征匹配点唯一。
在本实施方式中,手术导航系统100还包括导航探针50,导航探针50包括探头51,探头51用于在手术部位上点取比较点,导航探针50包括探针空间坐标系O4,处理装置30用于,对探针空间坐标系O4与相机空间坐标系O1进行坐标系转换,进而计算获取比较点在扫描图像坐标系O中的特征匹配点,以实现手术导航系统100术中导航。
如此,医生可通过观察比对导航探针50的探头51点取的比较点在手术部位的实时图像中的位置和手术部位的实际位置,进而验证手术导航系统100的配准关系是否正确,进而降低手术导航出现的风险。
本实施方式中,导航探针50上设有至少三个不共线的第三标志52,第三标志52用于建立探针空间坐标系O4。第三标志51为近红外发光二极管,如此,第三标志52发出不可见光,从而避免可见光对人眼产生刺激。立体相机20根据双目成像原理便可计算跟踪第三标志52的位置,从而确定相机空间坐标系O1与探针空间坐标系O4之间的坐标系转换关系。
在本实施方式中,比较点在扫描图像坐标系O中的特征匹配点根据以下公式(八)计算:
(八)
式中:
为扫描图像坐标系O中的任一点的坐标;
为比较点在探针空间坐标系O4中的坐标;
R为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的旋转矩阵;
R12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的旋转矩阵;
为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的旋转矩阵的逆矩阵;
R13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的旋转矩阵;
为任意时刻探针空间坐标系O4与相机空间坐标系O1之间的旋转矩阵的逆矩阵;
T14(t)为任意时刻探针空间坐标系O4与相机空间坐标系O1之间的平移向量;
T13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的平移向量;
T13为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的平移向量;
T12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的平移向量;
T为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的平移向量。
具体地,相机空间坐标系O1中的任一点在探针空间坐标系O4中的特征匹配点通过以下公式(九)计算:
(九)
为探头51点取的任一点在探针空间坐标系O4中的坐标。
其中,及T14(t)通过追踪导航探针50上的第三标志52发出的信息计算确认。
进一步地,由于导航探针50的出厂设置,导航探针50的探头51在探针空间坐标系
O4中的坐标确认唯一,即为常量,因此,根据上述公式(九),导航探针50的
探头51在相机空间坐标系O1中的坐标通过以下公式(十)计算:
(十)
根据上述公式(一)及(十)可计算出公式(八),从而可确定任意时刻导航探针50的探头51点取的比较点在扫描图像坐标系O中的特征匹配点,也就是说,探头51点取的比较点都能够在扫描图像中确定关系唯一的特征匹配点,因此,实现手术导航系统100在手术过程中准确地导航,降低手术的风险。
请参阅图3,本发明实施方式的手术导航系统的配准方法,包括以下步骤:
S1,在受施手术者的手术部位固定设置动态基准;
S2,在术前获取手术部位的扫描图像;
S3,对扫描图像坐标系与相机空间坐标系进行坐标系转换;
S4,在术中获取手术部位的实时图像;及
S5,计算实时图像上的任一点在扫描图像的特征匹配点。
本发明实施方式的手术导航系统的配准方法可通过上述实施方式的手术导航系统100实现。
本发明实施方式的手术导航系统100的配准方法中,处理装置30预存有计算实时图像上的任一点在扫描图像的特征匹配点的算法,在手术进行时,手术部位的实时图像上的任一点都能够在手术部位的扫描图像中找到特征匹配点,从而实现术中手术部位的实际位置与术前的扫描图像之间的配准,进而准确地引导手术的进行。本发明实施方式的手术导航系统100在配准过程中能够减少人工操作的误差,即使在手术部位的面积很小的情况下也能够快速实现自动配准,缩短了手术时间,进而降低手术风险,符合微创手术的需求。
请参阅图4,在本实施方式中,步骤S3包括:
S31,对扫描图像坐标系O、相机空间坐标系O1、扫描空间坐标系O2及动态空间坐标系O3进行坐标系转换;及
S32,计算相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点。
如此,通过坐标系的转换计算,可准确地得到实时图像上的任一点在扫描图像的特征匹配点。
请参阅图5,本发明实施方式中,进一步地,步骤S31包括:
S311,计算扫描空间坐标系O2中任一点在扫描图像坐标系O中的特征匹配点;
S312,计算相机空间坐标系O1中任一点在扫描空间坐标系O2及动态空间坐标系O3中的特征匹配点;
S313,计算动态空间坐标系O3中任一点在扫描空间坐标系O2中的特征匹配点;
S314,计算动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点;及
S315,计算任意时刻动态空间坐标系O3中任一点在相机空间坐标系O1中的特征匹配点。
具体地,在步骤S311中,通过以下公式计算扫描空间坐标系O2中任一点在扫描图像坐标系O中的特征匹配点:
式中:
为扫描图像坐标系O中的任一点的坐标;
为扫描空间坐标系O2中的任一点的坐标;
R为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的旋转矩阵;
T为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的平移向量。
其中,R和T可根据扫描装置10的成像原理计算得到,如此,确定扫描空间坐标系O2与扫描图像坐标系O之间的转换关系。
在步骤S312中,通过以下公式计算相机空间坐标系O1中任一点在扫描空间坐标系O2的特征匹配点;
式中:
为相机空间坐标系O1中任一点的坐标;
R12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的旋转矩阵;
T12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的平移向量。
其中,R12和T12可根据立体相机20的双目成像原理计算而得。
通过以下公式计算相机空间坐标系O1中任一点在动态空间坐标系O3中的特征匹配点;
式中:
为动态空间坐标系O3中的任一点的坐标;
R13为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的旋转矩阵;
T13为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的平移向量。
在步骤S313中,通过以下公式计算动态空间坐标系O3中任一点在扫描空间坐标系O2中的特征匹配点;
式中:
为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的旋转矩阵的逆矩阵;
在步骤S314中,通过以下公式计算动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点:
在步骤S315中,通过以下公式计算任意时刻动态空间坐标系O3中任一点在相机空间坐标系O1中的特征匹配点:
式中:
R13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的旋转矩阵;
T13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的平移向量。
如此,在手术过程中,当动态基准40的位置发生改变的时候,也能快速地确认任意时刻手术部位在动态空间坐标系O3中任一点在扫描图像坐标系O中的特征匹配点。
在本实施方式中,相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点根据以下公式计算:
式中:
为扫描图像坐标系O中的任一点的坐标;
为相机空间坐标系O1中任一点的坐标;
R为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的旋转矩阵;
R12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的旋转矩阵;
为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的旋转矩阵的逆矩阵;
R13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的旋转矩阵;
T13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的平移向量;
T13为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的平移向量;
T12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的平移向量;
T为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的平移向量。
如此,可确定相机空间坐标系O1中的任一点在扫描图像坐标系O中的特征匹配点且关系唯一,也就是说,立体相机20获取手术部位的实时图像上的任一点都能够在扫描图像坐标系O中确认唯一的特征匹配点,即手术部位的实时图像上的任一点在扫描图像中的特征匹配点唯一。
请参阅图3,在本实施方式中,手术导航系统的配准方法还包括步骤:
S6,在手术部位上点取比较点;
S7,对探针空间坐标系O4与相机空间坐标系O1进行坐标系转换;及
S8,计算获取比较点在扫描图像坐标系O中的特征匹配点。
如此,医生可通过观察比对比较点在手术部位的实时图像中的位置和实际中的手术部位的位置,进而验证手术导航系统100的配准关系是否正确,进而降低手术导航出现的风险。
在本实施方式中,比较点在扫描图像坐标系O中的特征匹配点根据以下公式计算:
式中:
为扫描图像坐标系O中的任一点的坐标;
为比较点在探针空间坐标系O4中的坐标;
R为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的旋转矩阵;
R12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的旋转矩阵;
为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的旋转矩阵的逆矩阵;
R13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的旋转矩阵;
为任意时刻探针空间坐标系O4与相机空间坐标系O1之间的旋转矩阵的逆矩阵;
T14(t)为任意时刻探针空间坐标系O4与相机空间坐标系O1之间的平移向量;
T13(t)为任意时刻相机空间坐标系O1与动态空间坐标系O3的平移向量;
T13为在获取扫描图像时相机空间坐标系O1相对于动态空间坐标系O3的平移向量;
T12为在获取扫描图像时相机空间坐标系O1相对于扫描空间坐标系O2的平移向量;
T为在获取扫描图像时扫描图像坐标系O与扫描空间坐标系O2之间的平移向量。
如此,实现手术导航系统100在手术过程中准确地导航,降低手术的风险。
本发明实施方式的手术导航系统的配准方法未展开的部分请参照以上实施方式的手术导航系统100相同或相类似的部分,在此不再赘述。
在本发明的实施方式的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明的实施方式和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的实施方式的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个所述特征。在本发明的实施方式的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
在本发明的实施方式的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接或可以相互通讯;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明的实施方式中的具体含义。
在本发明的实施方式中,除非另有明确的规定和限定,第一特征在第二特征之“上”或之“下”可以包括第一和第二特征直接接触,也可以包括第一和第二特征不是直接接触而是通过它们之间的另外的特征接触。而且,第一特征在第二特征“之上”、“上方”和“上面”包括第一特征在第二特征正上方和斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”包括第一特征在第二特征正上方和斜上方,或仅仅表示第一特征水平高度小于第二特征。
下文的公开提供了许多不同的实施方式或例子用来实现本发明的实施方式的不同结构。为了简化本发明的实施方式的公开,下文中对特定例子的部件和设置进行描述。当然,它们仅仅为示例,并且目的不在于限制本发明。此外,本发明的实施方式可以在不同例子中重复参考数字和/或参考字母,这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施方式和/或设置之间的关系。此外,本发明的实施方式提供了的各种特定的工艺和材料的例子,但是本领域普通技术人员可以意识到其他工艺的应用和/或其他材料的使用。
在本说明书的描述中,参考术语“一个实施方式”、“一些实施方式”、“示意性实施方式”、“示例”、“具体示例”或“一些示例”等的描述意指结合所述实施方式或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施方式或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施方式或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施方式或示例中以合适的方式结合。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。
应当理解,本发明的实施方式的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本发明的各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (6)
1.一种手术导航系统,其特征在于,包括:
动态基准,所述动态基准用于固定设置在受施手术者的手术部位;
扫描装置,所述扫描装置用于在术前扫描所述手术部位以获取扫描图像,所述扫描图像包括扫描图像坐标系,所述扫描装置包括扫描空间坐标系,所述动态基准包括动态空间坐标系;
立体相机,所述立体相机用于在术中获取所述手术部位的实时图像,所述立体相机包括相机空间坐标系;及
与所述动态基准、所述扫描装置及所述立体相机连接的处理装置,所述处理装置用于处理所述扫描图像及所述实时图像,以对所述相机空间坐标系与所述扫描图像坐标系进行坐标系转换,及计算所述实时图像上的任一点在所述扫描图像的特征匹配点;
所述扫描装置上设有至少三个不共线的第一标志,所述第一标志用于建立所述扫描空间坐标系;所述动态基准上设有至少三个不共线的第二标志,所述第二标志用于建立所述动态空间坐标系;所述立体相机用于根据双目成像原理,通过追踪所述第一标志及所述第二标志发出的信号确认所述相机空间坐标系与所述扫描空间坐标系及所述动态空间坐标系之间的坐标系转换关系。
2.如权利要求1所述的手术导航系统,其特征在于,
所述处理装置还用于,对所述扫描图像坐标系、所述相机空间坐标系、所述扫描空间坐标系及所述动态空间坐标系进行坐标系转换,进而计算所述相机空间坐标系中的任一点在所述扫描图像坐标系中的特征匹配点,以得到所述实时图像上的任一点在所述扫描图像的特征匹配点。
3.如权利要求1所述的手术导航系统,其特征在于,所述处理装置还用于:
计算所述扫描空间坐标系中任一点在所述扫描图像坐标系中的特征匹配点;
计算所述相机空间坐标系中任一点在所述扫描空间坐标系及所述动态空间坐标系中的特征匹配点;
计算所述动态空间坐标系中任一点在所述扫描空间坐标系中的特征匹配点;
计算所述动态空间坐标系中任一点在所述扫描图像坐标系中的特征匹配点;及
计算任意时刻所述动态空间坐标系中任一点在所述相机空间坐标系中的特征匹配点。
4.如权利要求1所述的手术导航系统,其特征在于,所述相机空间坐标系中的任一点在所述扫描图像坐标系中的特征匹配点根据以下公式计算:
式中:
为所述扫描图像坐标系中的任一点的坐标;
为所述相机空间坐标系中任一点的坐标;
R为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的旋转矩阵;
R12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的旋转矩阵;
为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的旋转矩阵的逆矩阵;
R13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的旋转矩阵;
T13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的平移向量;
T13为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的平移向量;
T12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的平移向量;
T为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的平移向量。
5.如权利要求1所述的手术导航系统,其特征在于,所述手术导航系统包括导航探针,所述导航探针包括探头,所述探头用于在所述手术部位上点取比较点,所述导航探针包括探针空间坐标系,所述处理装置用于对所述探针空间坐标系与所述相机空间坐标系进行坐标系转换,进而计算获取所述比较点在所述扫描图像坐标系中的特征匹配点,以实现所述手术导航系统在术中导航。
6.如权利要求5所述的手术导航系统,其特征在于,所述比较点在所述扫描图像坐标系中的特征匹配点根据以下公式计算:
式中:
为所述扫描图像坐标系中的任一点的坐标;
为所述比较点在所述探针空间坐标系中的坐标;
R为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的旋转矩阵;
R12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的旋转矩阵;
为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的旋转矩阵的逆矩阵;
R13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的旋转矩阵;
为任意时刻所述探针空间坐标系与所述相机空间坐标系之间的旋转矩阵的逆矩阵;
T14(t)为任意时刻所述探针空间坐标系与所述相机空间坐标系之间的平移向量;
T13(t)为任意时刻所述相机空间坐标系与所述动态空间坐标系的平移向量;
T13为在获取所述扫描图像时所述相机空间坐标系相对于所述动态空间坐标系的平移向量;
T12为在获取所述扫描图像时所述相机空间坐标系相对于所述扫描空间坐标系的平移向量;
T为在获取所述扫描图像时所述扫描图像坐标系与所述扫描空间坐标系之间的平移向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510650347.5A CN106560163B (zh) | 2015-09-30 | 2015-09-30 | 手术导航系统及手术导航系统的配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510650347.5A CN106560163B (zh) | 2015-09-30 | 2015-09-30 | 手术导航系统及手术导航系统的配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106560163A CN106560163A (zh) | 2017-04-12 |
CN106560163B true CN106560163B (zh) | 2019-11-29 |
Family
ID=58485323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510650347.5A Active CN106560163B (zh) | 2015-09-30 | 2015-09-30 | 手术导航系统及手术导航系统的配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106560163B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108186117A (zh) * | 2018-02-28 | 2018-06-22 | 安徽大中润科技有限公司 | 一种分布式光学定位跟踪系统及方法 |
CN108969099B (zh) * | 2018-07-26 | 2020-02-18 | 艾瑞迈迪医疗科技(北京)有限公司 | 一种校正方法、手术导航系统、电子设备及存储介质 |
CN109700550B (zh) * | 2019-01-22 | 2020-06-26 | 雅客智慧(北京)科技有限公司 | 一种用于牙科手术的增强现实方法及装置 |
CN109785374B (zh) * | 2019-01-23 | 2020-12-04 | 北京航空航天大学 | 一种牙科增强现实手术导航的自动实时无标记图像配准方法 |
WO2020181498A1 (zh) * | 2019-03-12 | 2020-09-17 | 上海复拓知达医疗科技有限公司 | 体内导航系统和方法 |
CN110215284B (zh) * | 2019-06-06 | 2021-04-02 | 上海木木聚枞机器人科技有限公司 | 一种可视化系统和方法 |
CN112043382B (zh) * | 2020-09-27 | 2024-02-23 | 安徽埃克索医疗机器人有限公司 | 一种手术导航系统 |
CN112006776A (zh) * | 2020-09-27 | 2020-12-01 | 安徽埃克索医疗机器人有限公司 | 一种手术导航系统及手术导航系统的配准方法 |
CN112155732B (zh) * | 2020-09-29 | 2022-05-17 | 苏州微创畅行机器人有限公司 | 可读存储介质、骨建模配准系统及骨科手术系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130281821A1 (en) * | 2011-01-13 | 2013-10-24 | Koninklijke Philips Electronics N.V. | Intraoperative camera calibration for endoscopic surgery |
CN103110429B (zh) * | 2012-06-11 | 2015-12-16 | 大连理工大学 | 超声波探头的光学标定方法 |
CN104287830A (zh) * | 2013-07-18 | 2015-01-21 | 中国科学院深圳先进技术研究院 | 基于Kinect相机的术中实时注册方法 |
CN104434313B (zh) * | 2013-09-23 | 2019-03-01 | 中国科学院深圳先进技术研究院 | 一种腹部外科手术导航方法及系统 |
CN104008269B (zh) * | 2014-04-03 | 2017-05-03 | 北京航空航天大学 | 一种基于人工标记物的手术导航系统空间自动注册方法 |
-
2015
- 2015-09-30 CN CN201510650347.5A patent/CN106560163B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106560163A (zh) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106560163B (zh) | 手术导航系统及手术导航系统的配准方法 | |
JP7433932B2 (ja) | 定量的3次元イメージング及び手術用インプラントのプリント | |
US11357593B2 (en) | Endoscopic imaging with augmented parallax | |
CN112867427A (zh) | 使用方位和方向(p&d)跟踪辅助的光学可视化的计算机化断层摄影(ct)图像校正 | |
KR101572487B1 (ko) | 환자와 3차원 의료영상의 비침습 정합 시스템 및 방법 | |
Shahidi et al. | Implementation, calibration and accuracy testing of an image-enhanced endoscopy system | |
US20180116732A1 (en) | Real-time Three Dimensional Display of Flexible Needles Using Augmented Reality | |
Chu et al. | Registration and fusion quantification of augmented reality based nasal endoscopic surgery | |
US11026747B2 (en) | Endoscopic view of invasive procedures in narrow passages | |
CN109549689A (zh) | 一种穿刺辅助引导装置、系统及方法 | |
WO2017117517A1 (en) | System and method for medical imaging | |
EP2680755B1 (en) | Visualization for navigation guidance | |
KR20160138502A (ko) | 정량적인 3차원 모델의 3차원 영상과의 정렬 | |
CN109124764A (zh) | 手术引导装置以及手术系统 | |
JP6116754B2 (ja) | 低侵襲手術において画像データを立体視表示するための装置およびその装置の作動方法 | |
US20110105895A1 (en) | Guided surgery | |
US20100208958A1 (en) | Image processing device, image processing system, and computer readable medium | |
CN106456267A (zh) | 器械在视野中的定量三维可视化 | |
CA2808757A1 (en) | System and method for determining camera angles by using virtual planes derived from actual images | |
CN103948432A (zh) | 术中立体内窥视频与超声图像增强现实算法 | |
WO2007115825A1 (en) | Registration-free augmentation device and method | |
US20130293690A1 (en) | Medical device navigation system stereoscopic display | |
US9437003B2 (en) | Method, apparatus, and system for correcting medical image according to patient's pose variation | |
KR20170057141A (ko) | Ct 이미지를 위한 국소 적용 투명성 | |
WO2022150767A1 (en) | Registration degradation correction for surgical navigation procedures |
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 |