人体器官运动监测方法和人体导航系统的制作方法

人体器官运动监测方法和人体导航系统的制作方法
【专利摘要】人体器官运动监测方法和人体导航系统。提供了手术过程中实时监测人体器官的运动的人体器官运动监测方法和手术导航系统。该人体器官运动监测方法包括：获得从术前三维立体医学图像中识别的各个运动监测工具在影像坐标系中的第一位置和姿态；在手术中实时确定各个运动监测工具在定位坐标系中的第二位置和姿态；基于各个运动监测工具在影像坐标系中的第一位置和姿态以及在定位坐标系中的第二位置和姿态，实时计算定位坐标系到影像坐标系之间的最优坐标转换关系，以及基于该最优坐标转换关系计算各个运动监测工具在从定位坐标系到影像坐标系的坐标转换的总体误差；基于实时确定的各个时刻的各个运动监测工具的坐标转换总体误差，评估人体器官在该各个时刻相对于手术前扫描时刻的运动程度。
【专利说明】
【技术领域】
[0001] 本发明总体地涉及手术导航系统和方法，更具体地涉及医学过程中实时监测人体 器官的运动的人体器官运动监测方法、人体导航系统、计算机可读介质。 人体器官运动监测方法和人体导航系统

【背景技术】
[0002] 介入手术是现代外科手术的发展方法，介入手术与传统外科手术的区别是不需要 开刀，只需要一个很小的创口即可将特制的手术器械例如导管、冷冻针、射频消融振、导丝 等刺入到人体内的病灶部位或者手术靶点位置上，然后通过各种物理/化学作用达到治疗 的目的，从而解决过去需要开放式手术才能解决的肿瘤切除、组织活检、人工器材放置等问 题。
[0003] 在介入手术过程中，医生无法直接观察到病人体内的治疗部位或者病灶，需要借 助于医学影像引导下的计算机辅助手术导航技术。计算机辅助手术导航技术是集计算机科 学、人工智能、自动控制、影像处理、三维图形学、虚拟现实和临床治疗等多方面技术为一体 的交叉研究课题。手术导航技术使用多种模态的医学影像协助医生将手术器械直接穿刺到 病灶进行局部治疗，进而提高手术质量、减少手术创伤、降低患者痛苦。
[0004] 手术导航应用患者的医学影像及由其重构生成的三维模型来实时引导临床手术 的实施。患者的医学影像数据集例如通过CT (computerized tomography,计算机断层摄影） 或MRI (Magnetic Resonance Imaging,核磁共振成像）扫描手段来获得。手术导航系统把 患者术前医学影像数据和术中手术部位通过定位设备联系起来，并能够在软件界面里准确 地显示患者的解剖结构及病灶附近三维空间位置的细节。当手术器械指向患者身体内部的 任意部位时，它的坐标信息都会被导航系统实时获取，并显示在患者的三维模型上。这样即 使不用给患者开刀，医生也能够实时了解手术器械和肿瘤病灶之间的相对位置关系。
[0005] 相对于传统手术手段而言，手术导航具有下述优点：
[0006] 1.可通过三维可视化技术实时重构病灶，并显示手术视野周围的结构特征；
[0007] 2.可通过术前手术方案设计选择最合适的手术路径；
[0008] 3.可显示手术路径上可能遇到的组织结构；
[0009] 4.可显示重要的应回避的组织结构，比如血管、神经和骨骼等；
[0010] 5.可显示病灶的需要治疗的范围；
[0011] 6.可实时精确计算出手术器械的位置姿态并加以显示；
[0012] 7.可显示手术器械与病灶的空间位置关系，并指示手术器械前进的方向；
[0013] 8.可术中实时调整手术入路，从而更加精确地达到病灶部位。
[0014] 为了描述方便，下文以CT图像作为医学影像的例子，不过显然医学影像也可以例 如MRI影像的其它影像。
[0015] 实时手术导航涉及采用定位设备或者跟踪系统，例如，采用电磁 (EM, electromagnetic)跟踪器，来跟踪手术工具例如穿刺针的末端，从而将手术工具的位 置与手术前医学影像例如CT图像相关，并将融合图像显示给临床医生。为了实现CT空间 和跟踪器空间之间的数据整合，通常在导航之前执行基于参考标记的配准过程。从CT图像 中识别这些参考标记（外部参考标记或者内部解剖参考标记），并由经校准的跟踪探头触 碰这些参考标记，从而得到这些参考标记在跟踪器空间，下文称之为定位坐标系中的坐标。 之后，执行基于点的配准，以找到CT空间和跟踪器空间之间的坐标变换矩阵。由坐标变换 矩阵得到配准矩阵，该配准矩阵使得手术工具与手术前CT图像对准，从而能够基于手术工 具在定位系统中的位置信息，实现手术工具的图像与CT图像在CT坐标系中的准确融合，并 可视地显示给临床医生。


【发明内容】

[0016]目前手术导航技术已经广泛地应用于神经外科手术导航和骨科手术导航中。然而 这些手术都假设人体是一个完全静止的刚体，并且这种静止状态会从CT扫描初期保持到 整个手术过程的结束。虽然这种假设在大多数神经外科手术和骨科手术里可以近似假定成 立，但是它在胸腹腔介入手术领域却难以获得成立。由例如呼吸而导致的胸腹腔组织运动 (包括肿瘤组织和正常组织）往往可以达到2-4厘米的范围，而临床手术的精度要求往往是 在毫米级别的，这样的误差会导致目前主流的手术导航系统难以安全稳定地应用于胸腹腔 介入手术的临床实践中。
[0017] 鉴于上述情况，做出了本发明。
[0018] 根据本发明的一个方面，提供了一种手术过程中实时监测人体器官的运动的人体 器官运动监测方法，在手术前通过扫描已经获得了固定有两个或两个以上的运动监测工具 的患者的医疗部位的三维立体医学图像，该三维立体医学图像具有相关联的影像坐标系， 该人体器官运动监测方法包括如下步骤：获得从术前三维立体医学图像中识别的各个运动 监测工具在影像坐标系中的第一位置和姿态；在同样的运动监测工具以与手术前扫描时同 样的位置和姿态固定于患者身体上的状态下：实时确定各个运动监测工具在定位坐标系中 的第二位置和姿态，该定位坐标系是定位手术工具的位置和姿态过程中所参考的坐标系； 基于各个运动监测工具在影像坐标系中的第一位置和姿态以及在定位坐标系中的第二位 置和姿态，实时计算定位坐标系到影像坐标系之间的最优坐标转换关系，以及基于该最优 坐标转换关系计算各个运动监测工具在从定位坐标系到影像坐标系的坐标转换的总体误 差；以及基于实时确定的各个时刻的各个运动监测工具的坐标转换的总体误差，来评估人 体器官在该各个时刻相对于手术前扫描时刻的运动程度。根据本发明的另一方面，提供了 一种手术导航系统，包括：定位设备，用于跟踪手术工具和运动监测工具在定位坐标系中的 位置和姿态；两个或两个以上的运动监测工具，被固定于患者身体上，其在定位坐标系中的 位置和姿态能够被定位设备跟踪到；三维立体医学图像获得单元，用于获得手术前的固定 有运动监测工具的患者的医疗部位的三维立体医学图像，该立体医学图像具有相关联的影 像坐标系；手术导航工作站，进行术前三维立体医学图像与手术中的手术工具图像的配准 和组合，并在所连接的显示设备上进行可视化显示，以指导手术执行者的手术操作；其中， 所述手术导航工作站还通过下述操作监视人体器官的运动状态：从术前三维立体医学图像 中识别各个运动监测工具在影像坐标系中的第一位置和姿态；在手术中实时确定各个运动 监测工具在定位坐标系中的第二位置和姿态；基于各个运动监测工具在影像坐标系中的第 一位置和姿态以及在定位坐标系中的第二位置和姿态，实时计算定位坐标系到影像坐标系 之间的最优坐标转换关系，以及计算各个运动监测工具在从定位坐标系到影像坐标系的坐 标转换的总体误差；以及基于实时确定的各个时刻的坐标转换的总体误差，来评估人体器 官在该各个时刻相对于手术前扫描时刻的运动程度。
[0019] 所述手术导航工作站可以基于所述评估的人体器官在该各个时刻相对于手术前 扫描时刻的运动程度，确定手术工具介入人体病灶的时机，并向执行手术的操作者提示该 时机。
[0020] 所述手术导航工作站还可以配置为：确定坐标转换的总体误差小于预定阈值的时 亥IJ ;确定该时刻的定位坐标系到影像坐标系的坐标转换矩阵；以及基于该坐标转换矩阵， 实时确定手术工具在影像坐标系中的位置，并将手术工具的图像在该位置与三维立体医学 图像组合地显示在显示设备上。
[0021 ] 所述手术工具可以是穿刺针。
[0022] 所述每个运动监测工具可以具有至少四个能够被定位设备跟踪的标志点，其中任 意三个标志点不共线，其中通过在该三维立体医学图像中识别每个运动监测工具的各个标 志点来识别各个运动监测工具在影像坐标系中的第一位置和姿态；以及在手术中定位设备 通过跟踪每个运动监测工具的各个标志点来确定各个运动监测工具在定位坐标系中的第 二位置和姿态。
[0023] 所述运动监测工具可以被固定于受人体器官运动影响大的患者身体部位上。
[0024] 根据本发明的再一方面，提供了一种计算机可读介质，其上记载有计算机程序，该 计算机程序与手术导航系统结合使用，且当被处理装置执行时，执行下述操作：获得从术前 三维立体医学图像中识别的各个运动监测工具在影像坐标系中的第一位置和姿态，该术前 三维立体医学图像基于手术前通过扫描固定有两个或两个以上的运动监测工具的患者的 医疗部位得到的，该三维立体医学图像具有相关联的影像坐标系；在同样的运动监测工具 以与手术前扫描时同样的位置和姿态固定于患者身体上的状态下：实时确定各个运动监测 工具在定位坐标系中的第二位置和姿态，该定位坐标系是定位手术工具的位置和姿态过程 中所参考的坐标系；基于各个运动监测工具在影像坐标系中的第一位置和姿态以及在定 位坐标系中的第二位置和姿态，实时计算定位坐标系到影像坐标系之间的最优坐标转换关 系，以及基于该最优坐标转换关系计算各个运动监测工具在从定位坐标系到影像坐标系的 坐标转换的总体误差；以及基于实时确定的各个时刻的各个运动监测工具的坐标转换的总 体误差，来评估人体器官在该各个时刻相对于手术前扫描时刻的运动程度。
[0025] 该计算机程序当被计算装置执行时还可以执行下述操作：确定坐标转换的总体误 差小于预定阈值的时刻；确定该时刻的定位坐标系到影像坐标系的坐标转换矩阵；基于该 坐标转换矩阵，确定手术工具在影像坐标系中的位置和姿态，并将手术工具的图像基于该 位置和姿态与立体医学图像组合地显示在显示装置上。
[0026] 利用本发明实施例的手术导航系统和人体器官运动监测方法，能够解决或缓解胸 胸腹腔手术导航中人体由于呼吸运动而导致肿瘤位移、进而降低手术导航精度的问题；也 能够简化手术导航的流程，方便手术导航在多个领域的使用和推广。

【专利附图】

【附图说明】
[0027] 从下面结合附图对本发明实施例的详细描述中，本发明的这些和/或其它方面和 优点将变得更加清楚并更容易理解，其中：
[0028] 图1示出了根据本发明实施例的示例性手术导航系统100的配置框图。
[0029] 图2示出了根据本发明实施例的呼吸监测工具120的示例性结构示意图。
[0030] 图3示出了根据本发明实施例的手术导航方法200的示例性过程。
[0031] 图4示出了呼吸运动过程中的坐标转换的总体误差随时间变化曲线的一部分的 示意图。

【具体实施方式】
[0032] 为了使本领域技术人员更好地理解本发明，下面结合附图和【具体实施方式】对本发 明作进一步详细说明。
[0033] 在详细描述本发明实施例之前，为便于本领域技术人员理解和把握本发明，下面 描述以下本发明的主要思想。
[0034] 手术前的CT图像是在患者的某个呼吸状态时刻或者呼吸状态阶段采集到的。在 进行CT扫描的过程中通常不便于跟踪参考标记在定位坐标系中的位置和姿态。扫描CT的 时刻与实际手术的时刻可能间隔数天。
[0035] 在后续的手术过程中患者持续处于呼吸的周期运动过程中，理论上只有当患者的 某时刻的呼吸状态与手术前CT扫描时刻或阶段相一致或相近时，此时的参考标记在定位 坐标系中的位置和姿态才对应于先前从CT图像中识别出来的参考标记在图像坐标系中的 位置和姿态，因而此时得到的配准矩阵才是正确的，从而在此时刻才能精确地将手术工具 的位置变换到CT空间中，提供给临床医生的融合后的三维体模型的可视显示才是精确的。
[0036] 为此，发明人想到，可以监测各个监测工具（上有参考标记）在定位坐标系中的位 置和姿态（取向或方位角），优化一个统一的从定位坐标系到三维影像坐标系的坐标转换 矩阵，确定各个监测工具从定位坐标系到三维影像坐标系的转换的误差的总计（后文中有 时简称之为总体误差），认为误差较小的时刻表明此时的人体器官的位置和姿态接近于在 扫描CT时刻人体器官的位置和姿态，并以可视或者可听等形式将该信息提供给临床医生， 辅助医生进行手术。
[0037] 需要说明的是，下文描述中，以呼吸运动引起肺部及相关组织运动作为人体器官 运动的示例，不过本发明并不局限于此，而是适用于监测人体存在动态内部运动的其他过 程，例如心跳运动。本文中描述的实施例可以优选用于协助对肺或者接近肺的组织进行定 位，不过所述组织也可以位于其他位置，例如心脏、消化器官、血管、肾脏等。
[0038] 在下文描述中，位置和姿态（取向）可以用6个数值来表示，即表示三维位置的位 置坐标和表示取向的三个角度值，这在某些传统技术中也被称为6D位姿。不过根据需要， 位置和姿态可以用更少的数值来描述。
[0039] 在本文中，"配准"表示在两个坐标系之间建立映射关系或者建立坐标转换关系， 例如"定位坐标系和影像坐标系之间的配准"表示在定位坐标系和影像坐标系之间建立映 射关系，或者求得两者之间的坐标转换矩阵，这里的定位坐标系是定位设备定位手术工具 和运动监测工具中所参考的坐标系。
[0040] 图1示出了根据本发明实施例的示例性手术导航系统100的配置框图。
[0041] 如图1所示，手术导航系统100包括定位设备110、两个或两个以上的运动监测工 具120、三维立体医学图像获得单元130、手术导航工作站140。为了便于描述手术导航系统 100工作过程，图1还示出了与手术导航系统100协同工作的控制台210和医疗工具220。
[0042] 定位设备110配置为跟踪运动监测工具在定位坐标系中的位置和姿态。如图所 示，运动监测工具上能够放置定位传感器112,例如磁感应定位传感器。定位设备110借助 于该磁感应定位装置，能够跟踪各个运动监测工具在定位坐标系中的位置和姿态。
[0043] 定位设备110还配置为跟踪手术工具在第二定位坐标系中的位置和姿态。手术工 具例如为穿刺针，穿刺针上安装有定位传感器112,由此定位设备110能够借助于该定位传 感器跟踪穿刺针的位置和姿态。在一个示例中，一个定位设备能够同时跟踪8个定位传感 器的信号。
[0044] 定位设备110可以采用电磁（EM)跟踪技术、光学跟踪技术、光纤光栅跟踪技 术之一。作为商业可用电磁定位系统的例子，具有加拿大NDI公司的Aurora系统，美 国Asension公司的DriveBay,荷兰Philips公司的PercuNav系统、美国GE公司的 InstaTrak3500Plus 系统、美国 Medtronic 公司的 StealthStation AxiEM 系统和美国 Compass 公司的 Cygnus FPS 系统。
[0045] 两个或两个以上的运动监测工具120被配置为分别在手术前（具体地，CT扫描 前）和手术中但在手术工具介入之前以相同的位置和姿态被固定于患者身体上。下文以监 测的运动为呼吸运动为例进行描述。为描述便利，下文可能将运动监测工具120称为呼吸 监测工具。呼吸监测工具在患者体表佩戴的位置最好选择受呼吸影响较大的地方，从而随 着呼吸运动几个呼吸监测工具之间的相对位置可能会发生一定程度的相对变化。
[0046] 在一个示例中，每个运动监测工具具有至少四个能够被定位设备跟踪的标志点， 其中任意三个标志点不共线，其中通过在该三维立体医学图像中识别每个运动监测工具的 各个标志点来识别各个运动监测工具在影像坐标系中的第一位置和姿态；以及在手术中定 位设备通过跟踪每个运动监测工具的各个标志点来确定各个运动监测工具在定位坐标系 中的第二位置和姿态。这里，在一个运动监测工具上设置至少四个标志点并且任意三个标 志点不共线是为了基于该至少四个标志点在影像坐标系和第一定位坐标系的位置和姿态 即能够求解出从定位坐标系到影像坐标系的坐标转换矩阵。在一个示例中，标志点例如为 金属小球，该金属小球在CT中的影像参数明显不同于人体材质的影像参数，很容易人工或 者通过图像处理将其自动识别出来。对于这样的每个呼吸监测工具，可以根据其在影像坐 标系的第一位置和姿态以及在定位坐标系中的第二位置和姿态，来确定一个唯一的坐标转 换矩阵，该坐标转换矩阵可以使得转换误差为零。但是因为有多个呼吸监测工具，要利用单 个坐标转换矩阵来进行该多个呼吸监测工具的每个在影像坐标系的第一位置和姿态和在 定位坐标系中的第二位置和姿态之间的转换，因此会产生转换误差。如前所述，本发明的主 要思想是通过计算各个时刻的最优统一转换矩阵和对应的总体转换误差来监测各个时刻 呼吸监测工具的位姿相对于手术前CT扫描时刻的位姿的差异，从而实现了监测各个时刻 人体器官的状态相对于手术前CT扫描时刻的状态的差异，达到了监测人体器官运动的目 的。
[0047] 需要说明的是，每个呼吸监测工具并不必需具有四个或更多个的标志点。而是只 要基于来自所有呼吸监测工具的全部标志点在定位坐标系和影像坐标系中的位置能够优 化求解两个坐标系之间的转换矩阵和总体转换误差即可。例如可以具有两个呼吸监测工 具，每个呼吸监测工具具有例如三个标志点，这样同样可以由总计六个标志点在定位坐标 系中的位置和在影像坐标系中的位置而求解最优坐标转换矩阵和对应的转换误差。再例 如，可以具有四个呼吸监测工具，每个呼吸监测工具具有2个标志点，这样同样可以由总计 八个标志点在定位坐标系中的位置和在影像坐标系中的位置而求解最优坐标转换矩阵和 对应的转换误差。一个更极端的例子是，具有5个呼吸监测工具，每个呼吸监测工具具有一 个标志点，这样同样可以由总计五个标志点在定位坐标系中的位置和在影像坐标系中的位 置而求解最优坐标转换矩阵和对应的转换误差。
[0048] 图2示出了根据本发明实施例的呼吸监测工具120的示例性结构示意图。该呼吸 监测工具120可以由塑料框架和被识别物组成。在塑料框架中，在标号121所指示的位置 处放置4个影像可识别标记物（比如金属小球，例如铅材质），这些标记物在CT或者MRI影 像里和周边材质的影像参数（主要是CT值和MRI值）有明显的区别特征，从而例如可以在 CT机内得到比较好的显影，可以通过已知方法方便的提取出来。作为标记物的另一个示例， 可以使用心电图制作中采用的电极贴。后续可以通过对CT影像或者MRI影像的后处理自 动从影像里提取可识别标记物，并进一步计算工具的位置和姿态（下文中，为描述方便，可 以简称之为位姿）信息（包括工具中心的坐标信息和方向信息）。
[0049] 图2中例示的呼吸监测工具120的正面可以放置磁感应定位装置，中间的圆孔122 以及左边的凹结构主要是方便磁感应定位装置的拆卸。
[0050] 三个延长的直杆主要是方便在人体表面附着，同时在杆上可以放置能够被CT或 者MR识别的标记物，识别后方能够计算位姿信息
[0051] 直杆末端的小横杆上的圆孔123可以通过记号笔在人体表面记录位置信息，这样 将呼吸监测工具120从患者表面拆卸下来以后，能够依据所记录的位姿信息来准确地将该 呼吸监测工具120再次粘贴到人体表面。
[0052] 图2所示的呼吸监测工具的结构和形状仅为示例，可以根据需要改变其配置。
[0053] 作为示例，呼吸监测工具的利用过程可以分为CT扫描阶段和手术中阶段。具体 地，在进行CT扫描前，将两个或两个以上的呼吸监测工具固定于患者胸部位置，并在人体 表面例如利用记号笔记录下佩戴的位姿信息，然后将患者推入CT室进行CT扫描，得到CT 扫描数据集，从这样的CT扫描数据集中可以识别出呼吸监测工具的位置和姿态。顺便说一 下，由这样的CT扫描数据集经过图像处理可以重构得到患者相关医疗部位的三维立体模 型，并且可以将这样的三维立体模型可视化于显示设备上，并与被跟踪的医疗工具的图像 融合显示，以辅助医生进行手术。在手术中并且在手术工具介入人体之前，按照前述记录在 患者身体表面的位姿信息来佩戴呼吸监测工具，实时监测各个呼吸监测工具在定位坐标系 中的位姿，并且优化求得针对所有呼吸监测工具的坐标转换矩阵和对应的总计误差，基于 总计误差的随时间变化来监测人体器官的运动。
[0054] 三维立体医学图像获得单元130可以是手术导航工作站140的一部分，或者是单 独的设备。三维立体医学图像获得单元130获得手术前的固定有运动监测工具的患者的医 疗部位的三维立体医学图像，该立体医学图像具有相关联的影像坐标系。例如手术前CT扫 描得到的三维数据集被存储于数据库中，三维立体医学图像获得单元130能够操作地耦合 到该数据库，并从该数据库中读取特定的三维数据集，并通过图像处理诸如渲染等可视地 显示到显示设备上。三维立体医学图像获得单元130可以是简单地读取已经重构好的三维 立体模型并显示，也可以是读取重构前的CT图像数据集并进行重构处理来得到三维立体 模型并显示。有关从CT或MRI扫描数据集重构三维立体模型的方法可以参考有关坐标转 换和误差计算的介绍，例如可以参考作者为翟伟明的题为"影像引导下计算机辅助介入手 术导航关键技术的研究"的博士论文中第三章和第四章的介绍，这里将该部分内容通过引 用合并于此。
[0055] 控制台210可以连接至手术导航工作站140或者是手术导航工作站140的一部 分。控制台210可以接收来自定位设备110的呼吸监测工具的位姿信息和手术工具的位姿 信息。控制台210还可以依据临床医生的指示来发出操作穿刺针的指令。
[0056] 手术导航工作站140配置为进行术前三维立体医学图像与手术中的手术工具图 像的配准和组合，并在所连接的显示设备上进行可视化显示，以指导手术执行者的手术操 作。
[0057] 手术导航工作站140可以包括处理器141和存储器142,存储器142中可以存储有 应用程序1421和操作系统1422。应用程序1421可以包括程序代码，该程序代码当被处理 器141执行时可以执行下面将详细描述的人体器官运动监测方法。这里的处理器是广义的 概念，可以通过单个专用处理器、单个共享处理器或者多个个体处理器来实施，其中多个处 理器中的一个处理器可以共享，此外不应将词语"处理器"或"控制器"的使用解释为排他 性地指代能够执行软件的硬件，而是其能够隐含地包括但不限于数字信号处理器（"DSP"） 硬件、用于存储软件的只读存储器（"RAM"）和非易失性存储器。手术导航工作站140还可 以包括用于与显示器150实时交互的接口。
[0058] 显示器150可以在手术导航工作站140的控制下显示医疗目标组织的取向、形状 和/或位置等，也可以显示患者有关部位的三维立体模型和手术工具的图像的融合。
[0059] 手术导航工作站140还可以包括各个人机交互用接口，例如连接到键盘、鼠标、触 摸屏系统等外围设备。手术导航工作站140还可以包括用于与三维立体医学图像获得单元 130交互的接口、用于与控制台210交互的接口等。
[0060] 为便于理解，下面概要描述一下临床手术流程。一般地，临床手术流程可以包括术 前准备、术前影像扫描、术前手术规划、术中手术导航、术后手术评估等处理。术前准备例如 包括导航手术室设备准备、使用真空垫固定患者、往患者体表固定运动监测工具（或者贴 标记点等）。术前影像扫描包括例如对患者进行术前CT扫描或者MRI扫描。术前手术规 划可以包括医学影像的三维可视化、器官及病灶的三维重构、交互式手术路径计划以及手 术效果模拟等，其中手术路径计划例如包括穿刺针的入针点、入针角度的选择；手术效果模 拟可以包括热场模拟、损伤场模型重构等。术中手术导航包括按照手术计划中的路径完成 穿刺等。术后手术评估可以包括术后规定事件重新扫描CT影像或者MRI影像，并在CT数 据上测量实际损伤区域的尺寸并和手术规划计算的数据进行比较，以衡量手术的精度和效 果。
[0061] 根据本发明实施例的人体器官运动监测方法例如可以用在术中手术导航过程中。
[0062] 下面结合图3描述结合人体导航系统进行的示例性手术导航过程。
[0063] 图3示出了根据本发明实施例的手术导航方法200的示例性过程。
[0064] 在步骤S210,在手术前将两个或两个以上的运动监测工具固定于患者身体上。
[0065] 在一个示例中，在患者身体表面放置2-4个呼吸监测工具，佩戴的位置最好选择 受呼吸影响较大的地方，使得随着呼吸运动，几个工具之间的相对位置会发生一定程度的 相对变化，由此对呼吸监测工具的位姿的跟踪更能代表呼吸运动的特性，也就更能描述随 着呼吸运动而运动的目标医疗部分（病灶）的特性。
[0066] 不过上述呼吸监测工具的个数仅为示例，可以根据需要选择更多数目个呼吸监测 工具，例如5个、6个乃至更多个等等。
[0067] 在步骤S220中，在手术前扫描以获得固定有运动监测工具的患者的医疗部位的 三维立体医学图像，该三维立体医学图像具有相关联的影像坐标系。手术前扫描例如为CT 扫描或者MRI扫描。三维立体医学图像例如由图1所示的三维立体医学图像获得单元130 来基于CT扫描图像数据集而重构得到。
[0068] 在进行术前扫描并获得固定有运动监测工具的患者的医疗部位的三维立体医学 图像后，这样的立体医学图像可以可视地显示在显示设备上，临床医生可以参考该立体医 学图像例如借助于导航软件系统进行手术路径规划。
[0069] 在步骤S230中，例如由手术导航工作站140从术前三维立体医学图像中识别各个 运动监测工具在影像坐标系中的第一位置和姿态。
[0070] 如前所述，各个运动监测工具上可以布置有标记点，该标记点在CT影像中具有明 显区别于人体组织的影像参数，从而能够人工或者自动地从CT影像中识别出各个标记点， 从而能够识别出各个运动监测工具在影像坐标系中的位置和姿态。
[0071] 需要说明的是，步骤S220与步骤S230之间以虚线箭头连接，表示两者之间可以间 隔较长时间，例如几个小时或数天。
[0072] 在步骤S240中，在手术中，在同样的运动监测工具以与手术前扫描时同样的位置 和姿态固定于患者身体上的状态下，实时确定各个运动监测工具在定位坐标系中的第二位 置和姿态，该定位坐标系是与定位设备跟踪手术工具的位置和姿态相关联的坐标系。例如， 利用图1中所示定位设备和各个运动监测工具上的定位传感器来得到各个运动监测工具 在定位坐标系中的第二位置和姿态，并传递给手术导航工作站140。
[0073] 在步骤S250中，手术导航工作站基于各个运动监测工具在影像坐标系中的第一 位置和姿态以及在定位坐标系中的第二位置和姿态，实时计算定位坐标系到影像坐标系之 间的最优坐标转换关系，以及基于该最优坐标转换关系计算各个运动监测工具在从定位坐 标系到影像坐标系的坐标转换的总体误差。
[0074] 有关坐标转换和误差计算的介绍，例如可以参考作者为翟伟明的题为"影像引导 下计算机辅助介入手术导航关键技术的研究"的博士论文中第6. 4节的介绍，这里将该部分 内容通过引用合并于此。
[0075] 在步骤S260中，基于实时确定的各个时刻的各个运动监测工具的坐标转换的总 体误差，来评估人体器官在该各个时刻相对于手术前扫描时刻的运动程度。
[0076] 如前所述，理论上，在肺部的呼吸运动周期的各个阶段中，如果某时刻的运动状态 与CT扫描时刻的运动状态最接近，则该时刻的呼吸监测工具在影像坐标系中的映射应该 是最接近于手术前CT扫描获得的呼吸监测工具在影像坐标系中的位姿的，因此与手术中 该时刻优化得到的坐标转换矩阵相对应的坐标转换误差应该是最小的。
[0077] 图4示出了呼吸运动过程中的坐标转换的总体误差随时间变化曲线的一部分的 示意图。其中横轴表示时间，纵轴表示坐标转换总体误差。如图4所示，在时刻tl、t2、t3,坐 标转换总体误差较小，此时的呼吸状态与手术前扫描时刻的呼吸状态较接近，而其它时刻 的呼吸状态不同程度地背离手术前扫描时刻的呼吸状态，坐标转换总体误差也随之变化。
[0078] 通过考察手术过程中坐标转化总体误差随时间的变化情况，能够评估人体器官的 状态在该各个时刻相对于手术前扫描时刻的状态的差异程度。
[0079] 可以以多种形式利用本发明实施例的人体器官运动监测结果。
[0080] 在一个示例中，可以将实时获得的坐标转换总体误差的随时间变化曲线以例如图 4的形式可视地显示在显示设备上。临床医生通过观察该总体误差随时间变化曲线，可以确 定例如穿刺针对于患者入针的时机。
[0081] 在一个示例中，例如图1的手术导航工作站基于所述评估人体器官在该各个时刻 相对于手术前三维扫描时刻的运动程度，确定手术工具介入人体病灶的时机，并向执行手 术的操作者提示该时机。在一个示例中，例如当某时刻的坐标转换总体误差为某一时段的 平均坐标转换总体误差的50%时并且变化趋势为降低时，确定为适于手术工具介入人体病 灶的时机，接着例如以短消息形式显示在显示设备上，或者以语音形式从扬声器输出，来提 醒临床医生。
[0082] 在一个示例中，因为描述呼吸监测工具的位姿所参考的定位坐标系即是手术中定 位手术工具的位置和姿态的过程中参考的坐标系，所以可以基于所评估的人体器官的位姿 在各个时刻相对于术前CT扫描时刻的位姿的运动（或差异）程度，而确定人体器官的位姿 与CT扫描时刻的位姿基本一致的时刻，并且利用所计算的该时刻的从定位坐标系到影像 坐标系的坐标转换矩阵，来将实时跟踪到的手术工具的位置转换到影像坐标系中并与立体 医学图像融合，进而可视地显示在显示设备上。具体地，该人体器官运动监测方法还可以 包括：确定坐标转换的总体误差小于预定阈值的时刻；确定该时刻的定位坐标系到影像坐 标系的坐标转换矩阵；基于该坐标转换矩阵，确定手术工具在影像坐标系中的位置，并将手 术工具的图像在该位置与立体医学图像组合地显示在显示装置上。即，在坐标转换的总体 误差小于预定阈值的时刻，进行手术工具所在的定位坐标系到影像坐标系的配准，并且基 于例如定位设备110跟踪到的手术工具的位置，将手术工具的图像与人体三维立体模型融 合，并可视地显示在显示设备上，以辅助临床医生手术。
[0083] 需要说明的是，该呼吸监测过程可以在手术中持续地进行。
[0084] 另外，由呼吸监测获得的规律，例如前一总体误差最小时刻和总体误差最小的时 刻之间的时间间隔（如图4所示）信息也可以提供给临床医生，让临床医生在下次入针之 前做好心理准备和操作准备。
[0085] 利用上述的手术导航系统和人体器官运动监测方法，能够不借助于专门的人体器 官监测仪器来监测人体器官的状态相对于CT扫描时刻的状态的运动程度（或者说差异程 度），例如对于呼吸运动而言，不需要专门的呼吸门控装置。不过，本发明并不排斥与专门的 人体器官监测仪器的协同使用。
[0086] 利用本发明实施例的手术导航系统和人体器官运动监测方法，不是简单地将人体 视为刚体来进行术前定位坐标系和影像坐标系的配准，而是考虑了人体的内部运动来监测 人体器官的状态相对于CT扫描时刻的状态的差异程度，从而能够提示医生恰当的操作时 机，使得能够更准确地执行术前的手术规划，减少对临床医生个人经验的依赖。
[0087] 此外，本发明的实施例可以采取计算机程序产品的形式。可经由计算机可读介质 访问所述计算机程序产品，以供计算机或任何指令执行系统使用或者与之结合使用。计算 机可读介质可以包括、存储、交换、传播或发送程序的设备。计算机可读介质可以是电子、 磁、光、电磁、红外或半导体系统（或器件或装置）或传播介质。计算机可读介质的例子包 括半导体或者固态存储器、磁带、可拆卸计算机软盘、随机存取存储器（RAM)、只读存储器 (ROM)、磁盘或光盘。
[〇〇88] 以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也 不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技 术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨 在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本【技术领域】的 其它普通技术人员能理解本文披露的各实施例。
【权利要求】
1. 一种手术过程中实时监测人体器官的运动的人体器官运动监测方法，在手术前通过 扫描已经获得了固定有两个或两个以上的运动监测工具的患者的医疗部位的三维立体医 学图像，该三维立体医学图像具有相关联的影像坐标系，该人体器官运动监测方法包括如 下步骤： 获得从术前三维立体医学图像中识别的各个运动监测工具在影像坐标系中的第一位 置和姿态； 在同样的运动监测工具以与手术前扫描时同样的位置和姿态固定于患者身体上的状 态下： 实时确定各个运动监测工具在定位坐标系中的第二位置和姿态，该定位坐标系是定位 手术工具的位置和姿态过程中所参考的坐标系； 基于各个运动监测工具在影像坐标系中的第一位置和姿态以及在定位坐标系中的第 二位置和姿态，实时计算定位坐标系到影像坐标系之间的最优坐标转换关系，以及基于该 最优坐标转换关系计算各个运动监测工具在从定位坐标系到影像坐标系的坐标转换的总 体误差；以及 基于实时确定的各个时刻的各个运动监测工具的坐标转换的总体误差，来评估人体器 官在该各个时刻相对于手术前扫描时刻的运动程度。
2. 根据权利要求1的人体器官运动监测方法，其中所述人体器官的运动包括肺部的呼 吸周期的各个阶段。
3. 根据权利要求1的人体器官运动监测方法，还包括基于所述评估人体器官在该各个 时刻相对于手术前三维扫描时刻的运动程度，确定手术工具介入人体病灶的时机，并向执 行手术的操作者提示该时机。
4. 根据权利要求1的人体器官运动监测方法，还包括： 确定坐标转换的总体误差小于预定阈值的时刻； 确定该时刻的定位坐标系到影像坐标系的坐标转换矩阵； 基于该坐标转换矩阵，确定手术工具在影像坐标系中的位置，并将手术工具的图像在 该位置与立体医学图像组合地显示在显示装置上。
5. 根据权利要求1的人体器官运动监测方法，所述手术工具是穿刺针。
6. 根据权利要求1的人体器官运动监测方法，所述每个运动监测工具具有至少四个能 够被定位设备跟踪的标志点，其中任意三个标志点不共线，其中通过在该三维立体医学图 像中识别每个运动监测工具的各个标志点来识别各个运动监测工具在影像坐标系中的第 一位置和姿态；以及在手术中定位设备通过跟踪每个运动监测工具的各个标志点来确定各 个运动监测工具在定位坐标系中的第二位置和姿态。
7. 根据权利要求1的人体器官运动监测方法，所述运动监测工具被固定于受人体器官 运动影响大的患者身体部位上。
8. 根据权利要求1的人体器官运动监测方法，所述定位设备为电磁跟踪器。
9. 一种手术导航系统，包括： 定位设备，用于跟踪手术工具和运动监测工具在定位坐标系中的位置和姿态； 两个或两个以上的运动监测工具，被固定于患者身体上，其在定位坐标系中的位置和 姿态能够被定位设备跟踪到； 三维立体医学图像获得单元，用于获得手术前的固定有运动监测工具的患者的医疗部 位的三维立体医学图像，该立体医学图像具有相关联的影像坐标系； 手术导航工作站，进行术前三维立体医学图像与手术中的手术工具图像的配准和组 合，并在所连接的显示设备上进行可视化显示，以指导手术执行者的手术操作； 其中，所述手术导航工作站还通过下述操作监视人体器官的运动状态： 从术前三维立体医学图像中识别各个运动监测工具在影像坐标系中的第一位置和姿 态； 在手术中实时确定各个运动监测工具在定位坐标系中的第二位置和姿态； 基于各个运动监测工具在影像坐标系中的第一位置和姿态以及在定位坐标系中的第 二位置和姿态，实时计算定位坐标系到影像坐标系之间的最优坐标转换关系，以及计算各 个运动监测工具在从定位坐标系到影像坐标系的坐标转换的总体误差；以及 基于实时确定的各个时刻的坐标转换的总体误差，来评估人体器官在该各个时刻相对 于手术前扫描时刻的运动程度。
10.根据权利要求9的手术导航系统，其中所述人体器官的运动包括肺部的呼吸周期 的各个阶段。
【文档编号】A61B5/11GK104055520SQ201410259145
【公开日】2014年9月24日 申请日期:2014年6月11日 优先权日:2014年6月11日
【发明者】翟伟明, 宋亦旭 申请人:清华大学