基于CT配准结果的光学导航定位系统及其导航方法技术方案

一种信息处理技术领域的基于CT配准结果的光学导航定位系统及其导航方法，该系统包括：术前CT图像导入模块、图像分割模块、体表初始配准模块、术前CT图像与术中二维超声图像模块和术中导航模块；本发明专利技术采用虚拟现实和术中超声相结合的方式，补偿由于呼吸等因素造成的术中定位误差，从而实现对冠脉搭桥术的目标点的精确定位导航。通过对术前心脏CT图像数据中的心脏和冠脉血管树进行手动分割重建，然后在光学导航仪和基于CT与超声的术中配准误差校正的帮助下搭建一个将内窥镜和虚拟内窥镜相融合的增强虚拟现实环境，从而可实现对冠脉搭桥术的目标点的精确定位导航。

【技术实现步骤摘要】
基于CT配准结果的光学导航定位导航方法
本专利技术涉及的是一种信息处理
的系统及方法，具体是一种用于辅助冠状动脉搭桥手术的基于术中超声与术前CT配准结果的导航定位系统及其导航方法。
技术介绍
近年来我国心血管病的发病率逐年上升，其中冠心病是最常见的心血管疾病。针对冠心病的治疗方法中，冠状动脉搭桥术是目前除了药物治疗和介入治疗之外的最主要、成熟治疗方式。然而传统冠状动脉搭桥手术需要胸骨正中切口，必要时还需要借助体外循环等操作才能完成手术。其包括切口大、恢复慢、并发症多等缺点。作为新型治疗方式的导航辅助微创冠脉搭桥手术，仅需要通过在胸壁上开几个手指粗细的切口、利用特殊的手术器械即可完成手术，达到了切口美观、创伤小、恢复快、并发症少等要求。进行导航辅助微创冠脉搭桥的主要手术难点之一是如何快速而精确定位目标点，不当或者错误的定位可直接影响到手术的成败和远期疗效。目前，微创冠脉搭桥手术的定位主要依赖术前影像，未能利用反映手术区域真实情况的实时信息，无法解决术中由于呼吸、体表标记点位移和体位变化等因素导致的误差问题，从而使得定位效果不理想。经过对现有技术的检索发现，“3D-imageguidanceforminimallyinvasiveroboticcoronaryarterybypass，HeartSurgForum2000-9732(3)”（三维图像指导机器人辅助微创冠状动脉搭桥术，心脏外科论坛，2000年）中TM.Peter最先尝试应用三维影像图像进行心脏外科手术术前规划的研究，他通过对患者术前采集的CT图像进行分割，得到心脏和骨架的表面模型，初步建立了一个微创搭桥手术的手术导航系统。该系统同时提供了一个虚拟内窥镜，利用术前手术配准的技术，使得真实内窥镜与模型试验中模型的相对位置关系和虚拟内窥镜与系统中建模出的心脏骨架模型的相对位置关系一致。但该系统只是一个原型系统，存在许多不足，如真实内窥镜2D图像没有与虚拟的3D场景有效融合等问题。“Flexiblecalibrationofactuatedstereoscopicendoscopeforoverlayinrobotassistedsurgery，MICCAI2002（1）：LNCS2488，T.DohiandR.Kikinis(eds)，25-34.”（机器人辅助手术中立体内窥镜的灵活标定，2002年医学影像计算与计算机辅助介入大会）中Mourguess和Coste-Maniere尝试着将一个动物的心脏在内窥镜下的影像与该心脏的冠脉树三维模型的虚拟内窥镜各自以一定透明度重叠进行融合，对手术中定位目标点起到了一定的效果。但是术中冠脉树的模型是静止的，而且冠脉树的三维模型与内窥镜影像的融合是依靠识别内窥镜中观察到的标示点来使二者较为准确的融合，这种方法的鲁棒性以及准确性并不十分令人满意。中国专利文献号CN1650813，记载了“外科手术导航系统基于光学定位的机器人手术定位方法”，该方法首先在机器人基座上设计三个标志位，通过光学跟踪仪定位指针来选取这三个标志，让机器人探针和光学跟踪仪定位指针对接时同时采集空间同一坐标系，最后通过光学跟踪仪与机器人系统的坐标互换。该技术主要使用与脑外科手术中，由于脑外科手术解剖学形变小，精确度问题上更容易把握；但该技术无法应用于解剖学形变较大的心血管系统领域。中国专利文献号CN101703409A，记载了“一种超声引导机器人辅助治疗的系统和方法”，该系统包括手术机器人、二维超声仪、磁定位仪和工作站，通过将他们相结合，实现了对病患区域的有效检测及自动化治疗，从而降低医生的劳动强度，并提高手术实施的精确度。但是此方法中使用的是磁定位仪，磁定位仪在使用时受到场地、器械、其他仪器等各种干扰，容易产生误差。同时三维重建的超声信息会产生图像变形问题，用其图像信息指导手术可能会造成偏差，而且产生信息量有限，无法用于一些较为精细复杂的手术，如心脏手术的精确操作。可见，虽然机器人辅助治疗系统已经开始应用于临床治疗中，但是其在手术目标定位上的发展尚不完善，很多临床难题仍有待解决。
技术实现思路
本专利技术针对现有技术存在的上述不足，提出了一种基于CT配准结果的光学导航定位系统及其导航方法，采用虚拟现实和术中超声相结合的方式，补偿由于呼吸等因素造成的术中定位误差，从而实现对冠脉搭桥术的目标点的精确定位导航。通过对术前CT图像数据中的心脏和冠脉血管树进行手动分割重建，然后在光学导航仪和基于术前CT图像与二维超声图像的术中配准误差校正的帮助下搭建一个将内窥镜和虚拟内窥镜相融合的增强虚拟现实环境，从而可实现对冠脉搭桥术的目标点的精确定位导航。本专利技术是通过以下技术方案实现的：本专利技术涉及一种基于CT配准结果的光学导航定位系统，包括：术前CT图像导入模块、图像分割模块、体表初始配准模块、术前CT图像与术中二维超声图像配准模块和术中导航模块，其中：术前CT图像导入模块接收术前影像学检查获得的DICOM格式图像文件，生成术前影像包并分别输出至图像分割模块和体表初始配准模块，术前CT图像与术中二维超声图像配准模块分别与图像分割模块和体表初始配准模块相连并接收带有手术目标点的三维动态冠脉树以及体表初始配准结果，术前CT图像与术中二维超声图像配准模块输出术前CT转换矩阵至术中导航模块，通过术中导航模块输出精确手术导航信息。所述的术前影像包包括：目标区域内一个或若干个心动周期的若干幅DIOM格式的术前CT图像以及根据DICOM数据重建得到的三维立体影像，其中：一个心动周期术前CT图像数据由若干个相位的术前图像组成。所述的带有手术目标点的三维动态冠脉树包括：由图像分割模块生成的动态心脏模型和血管树模型。所述的术前影像导入模块包括：DICOM图像读入单元和术前影像三维绘制单元，其中：DICOM图像读入单元将术前影像学检查获得的心脏一个心动周期的一系列DICOM格式图像文件导入解析出DICOM数据并传输至术前影像三维绘制单元，术前影像三维绘制单元根据DICOM数据重建成三维立体影像并合并为术前影像包并分别输出至图像分割模块和体表初始配准模块。所述的图像分割模块包括：心脏分割单元和血管树分割单元，其中：心脏分割单元对术前影像包中一个心动周期内的心脏CT图像基于临床经验逐帧手动勾勒出心脏轮廓，并根据获得的心脏CT图像分割结果重建得到动态心脏模型；血管树分割单元对术前影像包中一个心动周期内的血管树CT图像基于临床经验逐帧手动勾勒出血管树轮廓并根据临床需求手动标出手术目标点，从而根据获得的血管树CT图像分割结果重建获得血管树模型并与动态心脏模型一并作为带有手术目标点的三维动态冠脉树输出至术前CT图像与术中二维超声图像配准模块。所述的体表初始配准模块包括：配准标记点选择单元和转换矩阵计算单元，其中：配准标记点选择单元选择设置于手术对象体表上的配准标记点，并获得其术前CT图像数据在图像坐标系坐标以及真实手术对象的空间坐标，转换矩阵计算单元根据图像坐标系坐标以及空间坐标，利用刚体配准算法实现这两个坐标之间的配准，即实现图像空间与手术对象空间的配准，并利用这两组坐标计算出初始转换矩阵作为体表初始配准结果，输出至术前CT图像与术中二维超声配准模块。所述的术前CT图像与术中二维超声图像配本文档来自技高网...

【技术保护点】
一种基于CT配准结果的光学导航定位系统，其特征在于，包括：术前CT图像导入模块、图像分割模块、体表初始配准模块、术前CT图像与术中二维超声图像配准模块和术中导航模块，其中：术前CT图像导入模块接收术前影像学检查获得的DICOM格式图像文件，生成术前影像包并分别输出至图像分割模块和体表初始配准模块，术前CT图像与术中二维超声图像配准模块分别与图像分割模块和体表初始配准模块相连并接收带有手术目标点的三维动态冠脉树以及体表初始配准结果，术前CT图像与术中二维超声配准模块输出术前CT转换矩阵至术中导航模块，通过术中导航模块输出精确手术导航信息；所述的术前影像包包括：目标区域内一个或若干个心动周期的若干幅DIOM格式的术前CT图像以及根据DICOM数据重建得到的三维立体影像，其中：一个心动周期的术前CT图像数据由若干个相位的术前图像组成；所述的带有手术目标点的三维动态冠脉树包括：由图像分割模块生成的动态心脏模型和血管树模型。

【技术特征摘要】
1.一种基于CT配准结果的光学导航定位导航方法，其特征在于，包括以下步骤：第一步、图像分割模块中的心脏分割单元通过术前CT图像导入模块获得目标区域一个心动周期的一系列DICOM格式的术前CT图像数据后，对该心动周期内的心脏术前CT图像基于临床经验逐帧手动勾勒出心脏轮廓，获得分割结果，并根据分割结果重建得到动态心脏模型；第二步、图像分割模块中的血管树分割单元对术前影像包中一个心动周期内的血管树术前CT图像基于临床经验逐帧手动勾勒出血管树轮廓并根据临床需求手动标出手术目标点，获得分割结果，并根据分割结果重建获得血管树模型，构建出一个包含动态血管树模型的三维虚拟场景，该三维虚拟场景构成了一个虚拟内窥镜图像；第三步、体表初始配准模块利用体表配准来获得术前CT图像数据在图像坐标系坐标与真实手术对象的空间坐标的转换矩阵，即两个不同空间坐标系下的坐标点，通过特征互相映射并实现一一对应关系，达到最终术前CT图像与手术对象实时二维超声图像的对应，即具体通过在术前CT图像中选取若干个配准标记点，在真实的空间找到与图像中配准标记点对应的点并利用光学导航仪得到这些点在真实空间的坐标，利用这两组不同坐标系，但是一一对应的点集的位置坐标，求得两个空间之间的转换矩阵；第四步、通过输入心电图信号将带有手术目标点的三维动态冠脉树中对应的若干个相位与当前心脏跳动相位相匹配，即以转换矩阵T作为初始转化矩阵，设Ti为在转换矩阵T的基础上针对每一个相位的术前CT图像进行校正后得到的转换矩阵，即针对一个相位的术前CT图像i，Ti×T能进一步的校正该相位CT图像与手术对象实时二维超声图像之间的配准误差，具体步骤包括：4.1)通过超声探头对心脏采集一系列手术对象的手术对象实时图像，其中每一幅采集到的手术对象实时图像都对应某一个相应的术前CT图像、每一个相位的术前CT图像都对应了一系列的二维超声图像，对于相位i的术前CT图像数据来说，提取出心脏内壁的表面轮廓，即通过一个标定好的超声探头来获取术中的二维超声图像，通过T将二维超声图像坐标系转换到术前CT图像数据坐标系，从而与术前CT图像进行融合；设Ti，i＝1，2，…N，其中N为一个心动周期的术前CT图像数据中术前图像的相位个数，为在转换矩阵T的基础上针对术前影像包中的每一个相位的术前CT图像进行校正后的转换矩阵，即针对一个相位的、与一组转换矩阵中的Ti相对应的术前CT图像i，Ti×T能够进一步的校正该相位CT图像与手术对象实时二维超声图像之间的配准误差；所述的标定好的超声探头是指：为了能将实时的二维超声图像融合进导航系统，需要求得从二维超声图像坐标系到导航仪坐标系的转换矩阵；设TMtd←ui是从二维超声图像坐标系到固定在超声探头上的追踪设备的坐标系的转换矩阵，TMui←td是从超声探头上的追踪设备坐标系到世界坐标系的转换矩阵，二维超声图像中的一个点的坐标可以通过下面的公式转换到世界坐标系下的坐标；其中，(uk,uv)是这个点在二维超声图像坐标系中的坐标，(sx,sy)是x轴和y轴的比例系数，(xw,yw,zw)是其在世界坐标系中的坐标；超声探头标定就是要求得从二维超声图像坐标系到固定在超声探头上的追踪设备的坐标系的转换矩阵TMtd←ui；4.2)对相位i所对应的每一幅二维超声图像提取出心脏内壁的轮廓，从这一系列的二维超声图像中提取的内壁组成了一组点集，从该点集中提取出特征点点集，从术前CT图像中提取的心脏内壁表面轮廓是另一组点集；4.3)通过迭代最近点算法将二维超声图像上的点集配准到术前CT图像上的点集，得到这两个点集之间的转换矩阵Ti，用于后续的精配准处理的起始矩阵，该迭代最近点算法将二维超声图像上的点集配准到术前CT图像上的点集，得到这两个点集之间的转换矩阵Ti，用于后续的精配准处理的起始矩阵，具体步骤包括：假设两个待配准的点集P和Q，pi和qi分别是两个点集中的点，i＝1…n，配准问题的关键就是求解最优解使得最小时的R和T；在完成了ICP初配准后，每一个相位的术前CT图像i都各自得到了一个转换矩阵Ti×T；之后在导航阶段，将会实时的连续进行精配准，每一次精配准的初始转换矩阵就是Ti×T；对于某一个相位的术前CT图像，目标在于寻找一个使相似性测度最大的最优的转换矩阵，记为T′i，使其能最终满足对由于呼吸等因素造成的误差，从而满足对目标血管的定位要求，在此过程中采用的相似性测度为归一化的互信息；归一化互信息为：其中：M为二维超声图像所在区域在术前CT上采样所得的图像灰度点集，R为实时的二维超声图像灰度点集，H(M)为M的香农熵，iM代表M图像像素点的灰度值，代表M图像中像素点灰度值为iM的概率；H(R)为R的香农熵，iR代表R图像像素点的灰度值，代表R图像中像素点灰度值为iR的概率；H(M,R)为M和R的联合熵，代表M图像中像素点灰度值为iM、R图像中像素点灰度值为iR的概率；对于待配准的M图像和R图像，设IM(XM)、IR(XR)分别为M、R的灰度函数...

【专利技术属性】
技术研发人员：蔡俊锋，罗哲，王肃，赵强，顾立栩，
申请(专利权)人：上海交通大学医学院附属瑞金医院，
类型：发明
国别省市：