专利名称：Rf消融计划器的制作方法
技术领域：
本申请涉及射频消融(RFA)，具体地涉及将RF电极"探针"插入到病人体内，以便 覆盖预期的消融区(ablation zone)，并通过提高其温度来杀死肿瘤。然而，会意识到，所述 的技术还可以应用于其它消融技术、其它外科手术技术及其它放射治疗技术。
背景技术：
探针可以在预定时间段(例如大约15分钟，或一些其它适合的时间段)内连接到 电源，并且以各种形状(但通常是球体或椭圆体)进行消融。当前，外科医生可以在内心里 做出与要消融的病变的位置有关的记号，并可以利用各种图像向导来插入探针。因为探针 非常昂贵，因此阻止了外科医生使用多种探针尺寸和结构来支持尝试着使用最少数量的探 针来消融组织块(mass)。 当肿瘤或病变超过了由单次消融所能够成功覆盖的尺寸时，治疗成功率就急剧下 降。可以认为在肿瘤周围的、经常是lcm的边缘是为了杀死极微小的癌症细胞所必要的。遗 留下肿瘤的任何部分不进行治疗，都会导致复发，并且常常是强烈的复发。此外，对肿瘤的 覆盖区域的确定涉及复杂的3D几何运算和呈现，这即使是对于最好的外科医生来说也是 困难的。再另外，每一次额外的消融都增加了外科手术、麻醉时间和花费，并提高了对病人 的风险。 会发生的其它困难涉及在3维空间中单个目标位置在内心里的复杂呈现，以及在 控制探针以使其准确地到达目标位置时的困难度，这种困难度增加了发生错误的额外可能 性。另外，消融形状常常与肿瘤的形状或大小不匹配，这就导致消融了围绕肿瘤周围的健康 组织。肿瘤还可以由两个或多个更小的肿瘤区组成，这些肿瘤区在位置上不相连的，但彼此 足够靠近，从而应该一起治疗它们。消融还会损害"关键区域"，这会导致对病人严重的伤 害。通常对探针进入角度进行选择，因为它们匹配成像系统，而不是因为它们使消融的数量 最少或者减小了过程的风险。此外，由于每一个外科医生都在内心里创建一幅画面和计划， 并手动地引导探针，因此就不存在再现的机会，而再现性是"基于证据的医学"(evidence based madicine)的一个重要方面。本申请提供了新的且改进的肿瘤消融系统和方法，克服 了以上提及的以及其它的问题。
发明内容
根据一个方面， 一种为用以对在病人体内的组织块进行消融的消融过程制定计划 的系统，包括对目标块的表示。目标体积包括肿瘤以及任何可任选的边缘。这个表示可 以是诸如从CT图像分割出的区域之类的输入的体积。可替换地，其可以通过用户界面来指 定，在用户界面中，用户使用绘图工具从图像数据的不同视图描绘出边界轮廓，并且系统计 算相关的3D体素。优化组件接收与这个块相关的图像数据，并确定用以覆盖这个块的消融 区域。可以对消融区域进行选择，以便它们使消融的数量最少(从而使过程时间最短)，或 者可以对消融区域进行选择以便使对相邻于目标体积的健康组织的损害最小。
5
在另一个方面，确定关键区域以便该计划制定能够避开它们。 在另一个方面，确定目标体积的由于与关键区域的相对位置的原因而难以达到的 一个子集。这个子集可以用于建议外科医生考虑例如借助于盐水注射将目标与关键区域分 离。 在另一个方面，可以确定可能的、避开了禁止区域的进入点和方向的组。其优点是 将方向和进入向量的组縮小到仅包括切实可行的那些组。这允许系统能够为外科医生提供 对于切实可行的入口的选择。可以提供对于每一个入口的消融及附带损害的数量，或者可 以仅建议切实可行的入口中的最佳入口 。 —个优点是可以确定在消融区域内的关键区域(例如骨骼、肠等)，可以采取动作 来避免其消融。 另一个优点在于使外科手术持续时间最短。
另一个优点在于使消融体积与目标准确地匹配。 另一个优点是确定了特定的、定量的目标位置及方向。这些特定值可以用于以诸 如磁性追踪器的追踪设备指导外科医生。它们还可以用于控制需要定量数据的、诸如机器 人之类的其它设备。 本领域普通技术人员在阅读并理解了以下详细描述后，会意识到本发明更多的优 点。
本发明可以采取各种组件和组件排列的形式，并可以采取各种步骤和步骤排列的
形式。附图仅是为了图示说明多个方面的目的，不应解释为限制本发明。 图1示出了根据本文阐述的一个或多个方面的一种射频(RF)消融计划器系统，其
用于产生一个用以执行一个或多个消融协议来治疗病人体内的肿瘤块或病变的计划。 图2示出了一种算法方法，其为每一次消融的消融探针和目标位置产生包括探针
进入角度在内的定量计划。 图3A、3B和3C显示了主动脉附近的示例性PTV，以及根据与主动脉腔之间距离而 造成的温度变化的曲线图。 图4A、4B、4C和4D示出了与消融计划制定相关的各种动作的图形表示。 图5示出了根据多个特征，一种用于产生消融计划制定解决方案的方法。 图6A和6B示出了在执行与消融计划制定相关的多个动作时，可以经由GUI显示
给用户的图像。 图7示出了多个椭球形消融体积的图像。 图8示出了可以用于产生消融椭圆体的消融探针。 图9A和9B示出了被计划目标体积(PTV)以及与诸如骨骼之类的关键区域重叠的 候选椭球形消融体积的计算机建模图像。 图10示出了一种使用多个重力点来进一步减小消融数量的方法。
具体实施例方式
通常将诸如加长的细长探针之类的消融器械插入到肿瘤、病变或要消融的其它组织中，使用高射频来加热探针尖端，以便将周围的组织加热到足以杀死在其中的细胞的温 度，该温度常常认为是50摄氏度。尽管本申请主要描述了可以用于包括肝脏、肾脏、胸部、 肺部等在内的许多部位的射频(RF)消融技术，但会理解也可以类似地为冷冻消融术、微波 及其它消融和治疗过程制定计划。 通常相对于探针尖端来定位消融区域，并且消融区域是球状体或椭圆体形状，注 意球体是具有相等a、b、c轴的椭圆体。当肿瘤大于指定探针尺寸的消融区域时，外科医生 就选择一个以上的探针位置来产生相互交叠以覆盖整个肿瘤块的多个消融区域。通常的消 融过程包括定义目标区域，将探针插入到预期位置，并对探针供电大约15分钟以使得探 针尖端变热。 规定覆盖了整个肿瘤块以及在肿瘤周围的缓冲区域(例如通常是一厘米等)的被 计划目标体积(PTV)。这确保了消融全部肿瘤细胞和在缓冲区域中发现的极微小的肿瘤细 胞，以便减少肿瘤的复发。 根据本文提出的各个特征，通过确定探针在不穿过关键区域或骨骼组织的情况下 到达肿瘤块所沿着的可能轨迹，来选择探针位置。随后由外科医生选择或者自动地选择一 条或多条轨迹，以便优化用以消融整个肿瘤块而执行的消融的数量。 定义一系列椭圆体体积来包围PTV。如以下详细描述的，基于计算机的优化技术根 据椭圆体的大小来覆盖PTV。另外，在仍然会覆盖整个PTV的同时，移动椭圆体以使得延伸 到PTV之外的健康组织中的椭圆体消融区域的量最小。随着对椭圆体的移动，可以发现多 余的椭圆体，并可以删除这些多余的椭圆体以减小消融PTV所需的消融的数量。此外，可以 重复这个技术，以使得用于填充PTV空间的椭圆体的数量最少，从而使得所执行的消融数 量最少。可以为用于探针插入的几条或全部可能的轨迹执行优化，以便确定使得所需消融 数量最少或者减小了健康组织的附带损害的一个(或多个)插入角度。另外，这个技术可 以确保所计划的消融椭圆体体积不会与诸如骨骼或其它敏感组织之类的关键区域相交。
图1示出了射频(RF)消融系统IO，其有助于产生用于执行一个或多个消融协议
以治疗病人体内的肿瘤块或病变的计划。可以通过准确地计划消融探针位置以使得没有剩 余肿瘤的任何未治疗部分，并且准确地执行这个计划，来实现对于大肿瘤的成功治疗。系统 10产生定量的消融计划，包括每一次消融的目标位置和方向。可任选地确定在身体外部上 通向目标的一个或多个进入点。该消融计划确保覆盖了肿瘤的所有区域，并报告使用特定 探针完成完整的消融所需的消融数量。系统10还利用优化技术来使得消融数量最少。由 于该计划是定量的，因此可以用机器人和/或用登记的图像向导(例如通过定量地追踪消 融探针)来执行该计划。 系统10包括消融组件12，其可操作地连接到优化组件14和诸如CT扫描器之类的 成像组件16中的每一个。系统10可以进一步或者可替换地包括超声成像组件、X光透视 成像组件、磁共振成像系统、影像存档和传输系统(picture archiving and communication system) (PACS)或任何其它适合的成像组件或传递系统。在一个实施例中，消融组件12 是RF消融系统，其包括电源、射频发生器、可操作地与之耦合的探针等，以及有助于将探针 插入肿瘤块中并将在相对于探针尖端的区域内的肿瘤块加热到足以杀死肿瘤细胞的温度 (例如大约50摄氏度)的任何其它适合的组件。在可替换的实施例中，消融组件12包括 高强度聚焦的超声组件(HIFU)，其通过使用超声的机械振动和/或加热特性来消融特定区域内的组织。在一些实施例中，即使超声不是严格"聚焦的"，但使用换能器元件阵列也可以 预测超声消融区域。成像组件16产生被重构为3-D图像的数据。可以使用多种算法来自 动地分割诸如病变、器官、关键区域之类的对象，或者可以用绘图工具手动地沿着每一个轴 进行分割。所述分割产生了与特定对象相关联的体积区域的描述。具体的，可任选地，可以 经由图形用户界面18(GUI)将体积呈现给用户。可以以预期的距离对该体积进行"增长"， 以便将肿瘤加上边缘都包括在结果产生的体积中。在本文中无论何时使用词语"肿瘤"，尤 其是与优化有关的时候，都假定为意思是"被计划目标体积"(PTV)，它覆盖了意图结合起来 实现完全覆盖的特定肿瘤加上边缘。诸如Pi皿acle之类的系统提供了从医学图像手动或 半自动地分割出对象的能力。Pinnacle还允许用户能够设定边缘，在所述边缘上定义了新 的体积(PTV)。优化组件14分析与PTV相关的信息，尤其是尺寸，并还以方向为给定消融 探针定义了一组消融位置。在一个实施例中，优化组件14确定可以覆盖PTV的最少数量 的消融。在另一个实施例中，优化组件14以方向确定避开了最健康的组织(即，使附带损 害最小)的消融位置。在另一个示例中，分割出额外的对象体积，其表示不进行消融的组织 的"关键区域"或骨骼，并且优化组件14在避开这些区域的同时，尝试产生最少的消融或使 附带损害最小。然而在一些情况下，优化组件14产生未消融的区域，用户在这些区域上受 到警告，并可以将这些区域显示在用户界面上。与消融技术、消融时间段、探针尺寸、温度、 PTV、椭圆体体积、椭圆体中心、PTV减小和/或调整、用于使PTV最小化的协议、探针进入角 度和/或位置相关的信息，以及与系统性能有关的任何其它适合的信息都存储在存储器20 中，存储器20也可操作地耦合到优化组件14。另外，优化组件14包括处理器21，用于执行 各种分析动作，执行一个或多个例程等。例如，处理器21可以执行诸如结合图2、5和10所 述的那些例程。根据一些方面，优化组件包括计算机，其包括一个或多个处理器和用于存储 与消融技术有关的信息的存储器。 可以由医生选择探针进入角度，并经由用户界面手动地或以图形方式输入。可替 换的，随着探针穿过身体到达肿瘤，系统10可以有助于选择探针进入角度以避开"关键区 域"。理论上，应避免受到消融的关键区域可以与会被探针进入路径穿透的其它关键区域不 同。在RF消融中，当前探针需要在皮肤与目标之间的物理通路。然而在诸如HIFU之类的 其它消融技术中，可以在不干扰中间组织的情况下到达目标。然而实际上，消融和进入路径 的关键区域常常是相同的。优化组件14选择不与关键区域相交的消融体积。在一些情况 下，关键区域可能会非常靠近肿瘤从而不容易到达肿瘤。在此情况下，高亮显示关键区域以 警告用户希望有或建议诸如盐水注射之类的额外准备过程，以便将肿瘤与关键区域分离 或热隔绝。 由优化组件14执行对探针进入点(可任选的)和消融点的选择(例如为了实现 预期的消融形状而以给定的偏移设置探针尖端的位置)。根据其它方面，优化组件14提供 了多个建议的进入点和/或消融点，经由GUI 18向用户呈现这些进入点和/或消融点以便 进行选择。以此方式，RF消融计划器系统10有助于通过精确地计划消融位置以便不遗留 未治疗的肿瘤部分且不损害关键组织，来成功地治疗大肿瘤，这个计划制定产生了定量的 目标，从而使得消融组件12能够使用追踪或控制来提高精度，以及以更多的可重复过程来 改进总体系统。 根据其它实施例，预先产生并存储肿瘤块的PTV，并稍后将其输入到优化组件14，以便进行消融过程计划制定。输入的PTV是PTV的三维表示，其是用诸如CT扫描器、超声 成像组件、X光透视成像组件、磁共振成像系统、影像存档和传输系统(PACS)或任何其它适 合的成像组件或传递系统之类的成像组件16产生的。 一旦接收到PTV，就为完全包围PTV 的一个或多个消融体积确定消融体积中心。可以将目标消融点(例如这些中心)输出到追 踪系统，追踪系统可以包括GUI，用以辅助外科医生沿着预期的插入轨迹来设置探针。附加 地或者可替换地，可以将目标消融点输出到GUI 18 (例如在有或没有追踪系统的情况下)， 用于在执行消融过程时由外科医生认可和/或使用。 图2示出了一种用于产生定量计划的算法方法22，该定量计划包括消融探针的进 入角度和每一次消融的目标位置。该算法确保全部肿瘤区域都被覆盖，并确定使用特定探 针进行完整肿瘤消融的多次消融。另外，该算法包括优化技术，其使得对给定肿瘤体积的消 融数量最少。在一个示例中，该算法确定覆盖被计划目标体积(PTV)的多个消融椭圆形。将 椭圆形的中心拉向置于中心的"重力点"，其位于质心处。通常，其会落在PTV内，然而对于 例如"C"形的一些不常见的形状，这个点实际上会位于PTV之外。随着将椭圆形拉向重力 点，删除了没有覆盖PTV的唯一性部分的那些消融椭圆体，以便确定最少数量的消融椭圆 体及其各自的中心。在群集了足够多的椭圆形以致于一个椭圆形的关键覆盖区域已经被其 它椭圆形覆盖的情况下，就会发生这种情况。 一旦将这些椭圆形拉到一起，就可以使附带损 害最小。 根据该方法，在24处，分割并确定肿瘤体积。肿瘤可以是任何形状，并且可以是不 相连的区域的集合。随后为消融确定PTV，例如其可以是肿瘤体积加上约几厘米厚的周围体 积。以下，图3显示了示例性肿瘤，它在所有维度上都增长了固定的尺寸。尽管为这个2D 示例显示了一个圆形，但在3D示例中可以使用球形。沿着肿瘤的边缘设置该中心，从而构 成了PTV边界。这也可以用其它技术来完成。如果是多个不相连的区域靠近在一起，理想 情况下是通过边缘将这些区域相连接，就可以获得最佳结果。如果它们不相连，那么出于计 划的目的，可以优选地将边缘相连的体积的每一个群集认为是单独的肿瘤。PTV表示具有围 绕肿瘤而添加的边缘的体积。在24处，通常由本文献中所述的以及在产品中可获得的多种 分割技术中的任何一种来确定接近肿瘤的关键区域和吸热体(heat sink)。每一个关键区 域、吸热体、肿瘤、PTV等都是感兴趣区(ROI)。机器分割可以是自动进行的、由手工描绘轮 廓执行的，或者采用这二者的组合执行的。关键区域是器官或生命关键的组织，热会对其造 成不利影响。实际上有两类关键区域。第一类是接近PTV的不应被消融的那些关键区域。 第二类是例如在RF消融探针进入身体时不应被穿透的那些关键区域。肠子是不应被穿透 或消融的结构。然而主要影响是在工作流程上。外科医生可以计算一组消融，随后进行关于 仅分割相关的(接近PTV或者接近插入路径)关键结构的额外工作，以确定他们如何改变 所建议的过程。例如，可以随后增加接近所建议的消融区域的神经，以便查明该计划必须如 何改变，并且确定该过程是否会严重威胁这些神经。分割其它结构，并确定为"吸热体"。例 如，尽管看起来应保护诸如主动脉之类的关键区域免于消融，但实际上它可以保护其自身，
因为血液流动非常迅速，以致于它"自我冷却"了。然而这种冷却为需要消融的附近结构带 来了问题，因为冷却与加热过程是竞争的，结果导致未消融的区域和癌症复发的风险。在当 前RF消融探针中，15分钟后接近探针处的温度高于消融区域边缘处的温度。因此，可取的 是将探针移近主动脉3mm，以使得主动脉表面附近的温度更接近55摄氏度。可以使用基于描述液体和气体流动的Navier-Stokes等式的计算流体动力学来估计由于接近具有特定 尺寸和血流速度的特定吸热体而损失的热量。 在26处，在病人皮肤上定义进入角度和/或一个或多个进入点。在一个实施例中， 利用光线行进协议(ray marching protocol)来确定进入点。例采用二相体积，将CT体素 标记为"自由的"或者标记为"关键区域"。使用诸如由Perlin介绍的光线行进协议算法 (例如，见K. Perl in禾口 E. M. Hoffert， "Hypertexture"， Computer Graphics, vol. 23， issue 3，第253-261页，1989)来确定在皮肤上的允许沿着一条路径将探针插入PTV的多个位置， 其中所述路径不穿过诸如骨骼之类的敏感区域或者关键区域。直观地，这类似于将光设置 在肿瘤的中心处，使得关键区域(例如，诸如骨骼等之类的实体)阻碍该光，并确定该光到 达皮肤的点。光线以线状射线形式从PTV的质量中心(质心)穿过3D图像"行进"，直到出 现以下三种情形之一 1)光线到达图像体积的边缘，于是它从PTV中心处在新方向上重新开始
2)光线到达皮肤或被认可作为进入点的另一个位置，于是指明x、y、z位置和光线 方向。这是可能的进入点，其可以以图形方式显示，或者存储在列表以供选择，或者可以对 其进行评估以确定从这个角度进行覆盖所需的消融数量。 3)光线到达被标记为"关键区域"的体素，于是从PTV中心处以一个新的方向开始 新的光线。 这个过程持续进行，直到评估了所有预期的角度为止。对可能角度的选择可以基 于对方向的强力离散化、随机选择、用传统协议縮小范围，或者仅由外科医生在GUI上进行 简单的选取。该光线行进过程类似于将微小的光置于PTV质心位置处的观念，其中，例如， 关键区域阻挡该光到达皮肤。因此，皮肤上被照亮的区域就是可能的进入点，而有阴影的区 域就不是可能的进入点。接入角度在临床上和计算上都是重要的。进入角度应是安全的， 如该步骤所述的，避开关键区域。在临床上，例如在肝脏中，优选的从一个外部位置仅剌穿 肝脏囊一次，并产生基本上平行的消融来覆盖PTV。这同时减小了该问题的计算复杂度，因 为可以在假定一个固定方向的情况下计算覆盖。由于光线行进比覆盖分析快，因此预先计 算切实可行的进入点的子集限制了覆盖计算的数量。 在28处，定义一组椭圆中心，以使得与各个椭圆中心相对应的多个3-D椭圆体消 融体积共同包围整个肿瘤块。椭圆具有由满足x7a+y7b+z7c = 1的标准a、b、c半轴值定 义的、被平移到所确定的中心并被旋转了一个给定进入角度(例如由26定义的角度)的几 何形状。例如，给定消融探针具有基于其a、b、c值和方向的已知的消融体积X。在一个实 施例中，Y个具有已知体积X的消融完全覆盖PTV。可以由用户选择探针尺寸(例如根据探 针的类型)，且系统可以针对所分配的或所推导出的方向来优化探针的位置设置。用于确定 椭圆中心的这个方法也防止了消融关键区域，并在图5中进一步进行描述。
在30处，将为给定探针尺寸确定的椭圆中心经由如图7所示的图形用户界面输出 给诸如外科医生等的用户，以允许用户评估和/或执行该模型。以此方式，方法22有助于 为消融过程制定计划，以便在不同点处执行多次消融，这进而允许对在空间上交叠的多个 椭圆体体积进行消融，以确保对非均匀形状的肿瘤体积进行消融，同时避开了消融对其而 言是有害的区域。外科医生还可以为没有在模型中表示的因素来调整椭圆。可以将椭圆中 心连同所指定的消融探针几何形状和方向传输到消融组件12。消融组件12随后可以协同追踪系统或机器人系统一起使用。追踪系统向外科医生实时提供关于工具或工具尖端的位 置与方向的可视化反馈。将目标连同实况追踪的位置一起进行显示，并且还常常是连同登 记的病人图像一起进行显示。可替换地，使用来自追踪系统的反馈和外科医生批准的前进， 为机器人提供了目标点，并且机器人可以实施机动。 设置多个消融椭圆的位置以使得它们完全覆盖PTV。这些椭圆通常与PTV重叠，并 且可以延伸超出PTV的边缘。 一旦例程28使得椭圆数量最少，这些椭圆的中心，即探针尖 端位置就是被可视化传输给外科医生的探针位置的其他标识。 图3A显示了靠近主动脉44的示例性PTV 40。在此情况下，主动脉具有随着如42 所示的距离而逐渐减弱的热效应。这个热变化是RF消融探针加热到至少50摄氏度的能力 以及随着消融时间过去主动脉的37摄氏度的冷却效应的最终结果。图3B显示了作为与主 动脉内腔之间的距离的函数的结果产生的温度变化的曲线图。内腔是血管的内部空间。计 划器可以通过确保PTV实际上向吸热区域伸展一直到对应于较高预期温度的距离，来部分 地应对这种吸热结构。在图3C中，调整PTV 40使其移近主动脉，造成了 PTV的扩张46。这 还可以通过手动地增加补偿"可视PTV"来完成，该补偿"可视PTV"根据外科医生的专家意 见，在主动脉方向上将PTV扩展到一个位置。 图4A-D示出了与为具有另外的不宜手术的肿瘤的病人计划消融过程有关的各种 操作的多个图示60。 图示60包括诸如用磁共振成像工具等产生的、具有腹部肿瘤的病人的多个视 图。在图4A中，CT图像62示出了轴切面；图像64显示了病人胸腔和上腹部的弧矢视图 (sagital view);图像66示出了病人胸腔和上腹部的冠状视图(coronal view)。
分割出肿瘤块体积68，并且其由被计划目标体积(PTV) 70包围，PTV70是可视化 的或者计算的肿瘤块68体积加上额外的安全边缘或者缓冲区72，所述安全边缘或者缓冲 区72如在肿瘤块68与PTV 70的周边之间的空间所示的。在手术过程中，这个边缘72通 常是超出肿瘤大约lcm，并可以由病理学家进行验证。 一个良好的外科手术结果是"边缘干 净的"情况的，换句话说，在病理学家检查时，通过手术去除的肿瘤的外边缘不存在癌细胞。 然而当使用RF或其他最小侵入式消融技术时，病理学家就不能够对边缘进行评估。
实质上，边缘试图补偿不同的可能的错误来源，但可以针对特定病人解剖结构和 诸如治疗是否是减轻的之类的问题对其进行调整。 一些错误来源包括
在成像时间与介入时间之间的病人运动，
定义肿瘤的准确边缘的困难性，
不能够成像围绕肿瘤的极微小的癌症， 难以量化覆盖肿瘤并使得附带损害最小的目标消融位置，以及
难以执行针对特定目标的过程。 这个边缘由外科医生定义，但常常是大约lcm，无论如何都在O到大约1. 5cm之间。 以这些方式中的任意一种来减小错误都可以减小边缘。 由图2的例程26来计算消融探针对肿瘤68中心的一组可能的皮肤进入点78和 方位角，或者可替换的，由外科医生手工输入该位置和方向。在图4B中示出了单个进入点 74。在图像76的箭头78处的加黑的(以前是红色的)区域来显示表示成功地(例如安全 地)穿过皮肤到达肿瘤68的皮肤表面点。
图5示出了由优化组件14执行的一种方法，用于产生消融计划制定解决方案，如 本文结合在前附图并根据多个特点所述。在80处，例程28产生用以包围用于病人体内肿 瘤的PTV的最初的一组候选消融椭圆体。用贪婪算法(greedy algorithm)来产生这组最 初的候选消融区域，以便从N个边界侧获得最大的覆盖。例如，在N = 6的情况下，PTV由 具有6个侧面的盒体或立方体来限定。然而，会意识到PTV可以由任何适合的多面体来限 定(包括但不限于四面体、八面体、十二面体、二十面体、k印ler-Poinsot固体、阿基米德固 体、正多面体、非正多面体等)，以便进一步优化消融计划制定和为根除肿瘤和/或PTV而执 行的消融的数量。 在82处，为PTV选择中心点或质心，它可以是由PTV(例如，肿瘤块加上缓冲区) 包围的组织的质量中心。在84处，判断是否已经访问或评估了所有候选消融椭圆体。如果 在84处确定已经评估了所有候选消融椭圆体，则该方法就跳到94，在94处，判断是否所有 候选消融椭圆体都需要重新评估。如果在94处的判断是否定的，那么该方法结束。如果在 94处的判断是肯定的，那么该方法就返回到84，重复进行候选消融椭圆体评估。
如果在84处判断一个或多个候选消融区域需要进行评估，那么就在86处选择后 续的候选，并为其确定唯一性覆盖区(UCA)。 UCA可以视为PTV中的仅仅或"唯一性地"由所 评估的特定候选椭圆体体积所包围的一个范围或区域。在88处，判断UCA是否等于0。如 果UCA等于O，那么所考虑的候选消融椭圆体就不会覆盖PTV中未被另一个候选体积覆盖的 任何部分，并在90处将其从候选列表中删除。该方法随后返回到84，以便判断是否存在任 何剩余的要评估的候选消融体积。 如果在88处的判断表明该候选体积的UCA不是0，那么该候选就是有效的，并将其 保留在候选列表中，并且该方法前进到92。在92处，确定一个点，该点表示在该候选椭圆体 仍包围UCA的同时，该候选椭圆体的中心可以占据的、距在82处所选择的中心点最近的位 置。会意识到，可以进一步约束对于仍允许椭圆体覆盖UCA的、距离中心点最近的椭圆体中 心位置的确定，以确保该椭圆体消融区域不覆盖诸如骨骼或其它如果被消融就会伤害病人 的区域之类的关键区域。 如果椭圆体中心不在所确定的点处，那么就移动该候选椭圆体，这触发了对于仍 在候选列表中的所有候选椭圆体的重新评估。在94处，判断在列表中的所有候选是否都要 进行重新评估。如果是(例如如果候选列表中的一个或多个候选被移动)，则就在96处重 新访问所有候选，该方法返回到84，以便对消融计划制定解决方案进行进一步重复和优化。 当所有候选的中心都尽可能地靠近中心点，同时仍覆盖它们各自的UCA时，该解决方案就 稳定了，且该例程结束。 根据其它特征，在计划并优化消融点时，可以考虑例如与肿瘤密度和/或增长活 动有关的功能数据。例如，如果在一个区域中的肿瘤的密度大于另一个区域，那么这个更 密的区域就会比不太密的区域需要更多的消融。根据另一个示例，可以确定肿瘤中与其余 部分相比呈现出更快增长的部分，并将该部分作为进行比肿瘤的剩余部分更强烈消融的目 标，以确保根除它。因此，在优化消融计划时，可以使用与肿瘤有关的功能数据，以便应对肿 瘤块中的相对强和弱。 再其它的变化例包括在优化期间执行自适应计划制定方案，其包括但不限于，考 虑关于在关键区域与PTV之间的接近程度和/或距离的先验知识，作为与肿瘤块有关的功能数据的函数的、对于给定探针的消融温度和持续时间，在靠近PTV处的吸热体的影响(例 如，会从消融区域吸走热量，从而影响PTV的消融的结构)等。 图6A，图像100示出了具有定义PTV 102的边缘的经3D分割的不规则形状的肿 瘤。该PTV的形状的大小和不规则性造成了传统消融过程的问题，因为其不能被这种大小 的单次消融所覆盖。即使足够大的消融是可能的，但如果没有理想地定形，它就会同时造成 大量的附带损害、杀死健康的组织。由于消融探针数量有限，其每一个都具有固定消融带， 因此该任务就是模拟每一种可能性，以便能够选择最佳的探针。然而，如图6B中图像106 所示，为给定消融探针计算多个椭圆体消融区域108。确定椭圆中心以覆盖PTV 102与上述 的图2的例程28有关。 产生第一组椭圆108，以覆盖与PTV的消融体积相对应的最可能的"叮咬 (bites)"。在一个实施例中，针对PTV和PTV边界框的6个面每一个上的相切点来定义边 界框或立方体。从这6个点中选择覆盖最大数量PTV体素(例如表示3-D空间中的点的体 积像素)的最大的叮咬，并从要覆盖的PTV中删除相应的体素。重复这个动作，直到PTV不 再有未覆盖的像素。至此产生的椭圆是用于消融的候选椭圆。 选择PTV 102中的一点，例如是质量中心，称为中心点104。例程28依次访问候 选椭圆列表中的每一个椭圆。对于每一个候选，确定椭圆对覆盖做出贡献的唯一性体积，称 为唯一性覆盖区(UCA)。如果UCA是零(O)，那么该椭圆就没有覆盖PTV 102的任何唯一性 的部分，这可以是由于其它椭圆的移动造成的。当出现这种情况时，就从候选列表中删除具 有零UCA的椭圆，并评估下一个椭圆。如果椭圆具有非零UCA，就执行二分法搜索(binary search)，以测试在当前椭圆中心与PTV 102的中心点104之间的位置，以找到椭圆可以移 动到的、最接近中心点104的、同时仍覆盖UCA的位置。向中心点104移动椭圆中心确保 了消融整个肿瘤块，同时使得对围绕肿瘤的组织的不希望有的消融最小。就是说，向中心点 104移动椭圆中心收縮了PTV 102以匹配肿瘤，同时使得对接近PTV 102周边的健康组织的 损害最小。如果任何椭圆都能够在仍覆盖其UCA的同时更加移近PTV102的中心点104，那 么就重新检查所有候选椭圆的位置。当在不损害椭圆各自的UCA的覆盖的情况下，没有椭 圆的中心能够更移近PTV 102的中心点104时，该消融计划解决方案就稳定了。可任选的， 可以以初始的覆盖椭圆组重复这个过程，从而得到了更少数量的椭圆。
在图7中是多个椭圆体消融体积108的图像IIO，每一个椭圆体都具有中心114， 如上结合在前附图所述的来确定该中心114，并且每一个椭圆体都具有作为所选择的探针 大小和消融持续时间的函数的已知体积。参考图8，由椭圆形区域124描绘了消融探针120 的消融图案。沿着从皮肤上的插入点到PTV的所选择轨迹将消融探针120插入到病人体内 的肿瘤块或PTV中，以便于消融全部或一部分PTV。根据一个示例，探针尖端的目标位置超 出了椭圆体的中点，因此将探针尖端的位置控制到超出椭圆中心114的距离122，以确保消 融出现在预期的位置处。会意识到距离122通常随着探针与探针的不同而改变，并不是旨 在限于图6中表示的实际和/或按比例縮放的距离。 图9A示出了 PTV的计算机建模图像和覆盖了诸如骨骼之类的关键区域136的多 个候选椭圆体消融体积。例如，在关键区域非常难以接近以致于不做出特别的努力就不能 到达肿瘤的一些较小部分的情形下，就执行动作以确保在将消融数量保持为最小的同时， 不消融关键区域136。在图像130中，区域134包含了多个区域，这些区域过于接近关键区域136，以致于不能进行消融。在图像132中(图9B)，将不能消融的区域高亮呈现给外科 医生(例如经由GUI 18等)，以使得外科医生能够采取适当的动作。例如，外科医生可以选 择另一不同的消融探针。在这个示例中，选择具有更小消融椭圆体体积的更小的探针，其椭 圆体体积不会与高亮的关键体积相交。根据另一个示例，外科医生将盐水注射到肿瘤与关 键区域136之间的空间中，以将肿瘤与关键区域136分离。如果为一个或多个消融选择了 另一不同的探针，那么就有利地重复图2的优化例程28。类似的，如果注射盐水以增大在肿 瘤与关键区域136之间的空间，那么就重新评估PTV，并重复优化例程28，以处理由盐水注 射而引起的变化。 图10示出了根据一些实施例，由优化组件14执行的一种方法，其用于为消融过程 制定计划。该方法允许将消融体积拉向质心点，就如同被重力拉一样，这有助于使得在典型 病变之外的组织的过度消融最小，同时提供了对区域进行交叠以便进行合并的可能性。该 方法可以通过创建额外的"重心"或点来进一步减小消融体积的数量，这进一步减小了椭圆 的数量。 在140处，确定多个"重力点"。在142处，判断是否要分析更多候选消融体积。如 果没有更多候选消融体积要进行评估，那么该方法结束。然而，如果剩余一个或多个其它消 融体积要进行评估，那么在144处，就选择下一个候选，并且定义和/或确定其UCA。在146 处，判断UCA是否等于O(例如，候选不具有唯一性覆盖区域)。如果是，那么就在148处从 列表中删除该候选，并且方法返回到142。如果该候选具有UCA，那么就在150处，就确定与 该候选消融体积的质心最接近的重力点。 在152处，确定最接近于重力点、同时仍覆盖UCA的消融体积。如果移动该消融体 积以确保仍覆盖UCA，那么条件"重新访问所有候选"就为真。在154处，判断条件"重新访 问所有候选"是否为真，如果是，则就在156处对候选消融体积的列表全部重新访问。在此 情况下，该方法重复进行。如果没有移动候选消融体积，那么该方法就结束。
以下论述提供了用于确定重力点的一种方式的示例。"要消融的肿瘤"因子 (TAfactor)可以定义为球形肿瘤半径/球形消融半径。因此，当TAfactor小于或等于1时， 一个消融就会覆盖肿瘤。当1 < TAfactor <= 1. 25时，那么6个消融覆盖肿瘤。当1. 25 < TAfactor <= 1. 66时，那么14个消融覆盖肿瘤，以此类推。由于该方法允许对任意形状 的肿瘤进行消融计划，因此就希望将更普遍且更灵活的椭圆体用于消融形状，而不是球形， 因为消融区域常常是椭圆形，但也可以是球形(具有相等a、b、 c轴的椭圆形)。为了简单 起见，以下示例以用于肿瘤和消融的两个不同大小的球形来描述该过程。
在该方法与图5的方法之间的一个重要差别是"重力点"的确定，这将消融从中心 拉向区域节点的集合，在这里这些消融可以进一步聚集，以使得消融数量最小。例如，将一 组6个球形(1 < TAfactor <= 1. 25)的重力点放置在沿着每一个X、Y和Z轴(同时在+ 和_方向上)距离肿瘤中心0. 76*消融-半径的距离Q处。例如，如果肿瘤的半径为2. 5， 消融的半径为2，那么TAfactor实际上就是1. 25。重力点的位置因此就是沿着3个轴的每 一个，距离肿瘤中心Q二 0.76*2 = 1.52处。如果肿瘤以(O，O，O)为中心，那么重力点就位 于:(0，0， 1. 52) 、 (0，0， -1. 52) 、 (0， 1. 52，0) 、 (0， -1. 52，0) 、 (1. 52，0，0) 、 (-1. 52，0，0)。
—组14个球形(1. 25 < TAfactor < = 1. 66)的重力点是沿x、 y、 z轴放置的6 个球形以及填充这6个球形之间的"角落"的8个球形的组合。前6个重力点在距离肿瘤中心较远的距离W = 0. 90*消融-半径处。例如，如果肿瘤的半径为3. 32并且消融的半径 为2，那么TAfactor实际上就是1. 66。类似地计算，前6个重力点位于距离半径3. 32的肿 瘤的中心W = 0. 9*a = 1. 8处。剩余的8个重力点在对角线上。其由沿着两个不同轴的相 等的绝对距离来构成。绝对距离L = W2sqrt(l/2)。在这个示例中，L = (1. 8)2sqrt(l/2) =3. 24*0. 07071 = 2. 291。对角线位于球形的以下位置:(L， L， L) 、 (L， L， _L) 、 (L， _L， L)、 仏，_L， _U 、 K， L， U 、 K， L， _U 、 K， _L， U 、 K， _L， _U 。在(1. 25 < TAfactor < = 1.66)的情况下，可以执行第一个合并，以便如上所述地将消融中心拉向真正的中心，随后 可以执行图10的步骤，以便将它们拉向最近的重力点。 如果以球体来限定肿瘤的范围，从而使其更普遍一些，且在消融探针产生球形消 融体积时，则使用该技术就是有利的。当TAfactor > 1. 66时也可以使用这个技术。当肿 瘤由球体来限定范围时，这尤其是正确的，对于消融数量的当前实际限制不仅限于每个消 融所需的15分钟，还在于放置的精度。随着消融技术变得更快，可以用成像或追踪来引导 探针放置，以便提高精度。 实际上，肿瘤可以是任意形状，消融形状可以是任意的。在这些情况下，可以通过 将肿瘤内切于球体中并且将消融内切于球体中，来确定重力点。半径的比值可以用于定义 TAfactor及用于计算重力点。这个简易性允许该技术更快地进行计算，同时适应于任意形 状的肿瘤和消融。快速计算随后允许评估多种情形，以便可以确定具有最少消融和/或最 少附带损害的进入点。 可替换的，可以使用更复杂的计算来定义重力点，在这里PTV被球体包围，计算椭
圆形消融以覆盖该球体。覆盖球体的可能的许多椭圆的质心是重力点。 根据另一个实施例，可以执行图5的步骤80，以创建初始消融，并且随后执行图10
的步骤。这个技术在大量情况下都是有效的，然而在其它情况下，优选的执行图5的步骤，
随后是图10的步骤。例如，对于针状(钉状)病变(一些乳腺癌及其它癌症的特性)，将覆
盖压縮到中心并且随后重新分布它，就可以减小消融的总数。
权利要求
一种用于为对病人体内的组织块(68)进行消融的消融过程制定计划的系统，包括图形用户界面(18)，向用户呈现所述块(68)的表示；以及优化组件(14)，接收与所述块(68)相关的图像数据，并确定用以覆盖所述块(68)的消融区域(108)。
2. 如权利要求l所述的系统，其中，所述优化组件(14)产生包含所述块(68)的被计划 目标体积(PTV) (40 、70 、102)。
3. 如权利要求2所述的系统，其中，所述优化组件(14)产生候选消融区域(112)的初 始列表，每一个候选消融区域(112)都具有特征点(114)，在所述特征点(114)处设置与消 融组件(12)相关联的探针(120)，以治疗所述块(68)在每一个消融区域(108)中的部分。
4. 如权利要求3所述的系统，其中，所述优化组件(14)为预定数量的候选消融区域 (108)确定唯一性覆盖区域(UCA)，以使得给定候选消融区域(108)的所述UCA包括所述 PTV(40、70、102)中的仅由所述给定候选消融区域(108)覆盖的部分。
5. 如权利要求4所述的系统，其中，如果一个候选消融区域(108)不具有UCA，则所述 优化组件(104)从候选消融区域(112)的列表中删除该候选消融区域(108)。
6. 如权利要求4所述的系统，其中，所述优化组件(14)执行二分法搜索算法(82)，以 确定所述给定候选消融区域(108)的所述特征点(114)是否可以在不损害所述给定候选消 融区域(108)的所述UCA的覆盖的情况下，移动靠近所述PTV(40、70、102)的中心点(104)。
7. 如权利要求6所述的系统，其中，所述优化组件(14)继续优化肿瘤块消融的模型，直 到没有候选消融区域(108)可以被移动靠近所述中心点(104)且所有UCA都被覆盖。
8. 如权利要求l所述的系统，其中，所述优化组件(14)利用光线行进技术在病人的皮 肤上产生至少一个候选插入点(74)，通过所述候选插入点沿着到所述PTV(40、70、102)的 轨迹插入与所述消融组件(12)相关联的消融探针(120)，同时避开了一个或多个关键区域 (136)。
9. 如权利要求1所述的系统，其中，所述优化组件(14)包括用于选择覆盖了被计划目标体积(PTV) (40、70、102)的初始的一组候选消融区域 (112)的例程(80);用于选择所述PTV(40、70、102)的中心点(104)的例程(82);用于确定在所述初始的一组候选消融区域(112)中是否还剩余一个或多个候选消融 区域(108)要进行评估的例程(84);用于选择候选消融区域(108)，并评估该候选消融区域(108)的唯一性覆盖区(UCA)的 例程(86);用于确定该候选消融区域(108)的所述UCA是否等于O的例程(88);用于在所述UCA等于O的情况下，从所述初始的一组候选消融区域(112)中删除该候 选消融区域(108)的例程(90);用于在所述UCA不等于O的情况下，确定这样的与所述中心点(104)最接近的一个位 置的例程(92) :S卩，所述候选消融区域(108)能够在仍覆盖所述UCA的同时被移动到该位 置；用于在已经将一个或多个候选消融区域(108)移动靠近所述中心点(104)的情况下， 确定所有候选消融区域(108)已经准备好进行重新评估的例程(94);以及用于使得所有候选消融区域(108)都被进行重新评估的例程(96)。
10. 如权利要求l所述的系统，其中，所述优化组件(14)根据与以下各项有关的功能数 据中的至少一个，来对消融计划进行优化所述块(68)、接近所述块(68)的一个或多个吸热体的位置、被选择来执行所述消融的消融探针(12、 120)的形状、被选择来执行所述消融 的消融探针(12、120)的大小、以及在消融期间所使用的不同消融探针(12、120)的数量。
11. 如权利要求l所述的系统，其中，所述优化组件(14)根据以下各项中的至少一个， 来对消融计划进行优化候选消融点、消融探针进入角度、消融温度、消融时间和关键区域。
12. 如权利要求l所述的系统，其中，所述优化组件(14)使用了关于以下至少一项的 先验知识PTV(70、102)与关键区域(136)的接近程度、与所选择的消融探针(12、120)有 关的消融时间、与所选择的消融探针(12、120)有关的消融温度、以及位于所述PTV(40、70、 102)附近的一个或多个吸热体。
13. 如权利要求l所述的系统，其中，消融区域(108)的最佳数量是覆盖所述块(68)的 消融区域(108)的最小数量。
14. 如权利要求l所述的系统，其中，用户通过在所述块(68)的所述表示的一个或多个 不同视图上勾画出边界，来输入包围了所述块(68)的被计划目标体积(PTV) (40、70、102)， 并且所述优化组件(14)产生表示被计划目标体积(40、70、102)的3D体素。
15. —种用于为消除病人体内的组织块(68)的消融过程制定计划的方法，包括 确定所述病人体内的块体积(68); 产生所述病人体内的所述块体积(68)的图像表示； 将与所述图像有关的信息呈现给用户；确定覆盖了所述块体积(68)的消融区域(108)的最佳数量。
16. 如权利要求15所述的方法，还包括产生包含所述块体积(68)的被计划目标体积 (PTV) (40、70、102)。
17. 如权利要求16所述的方法，还包括确定一组候选消融区域(112)，每一个候选消 融区域(112)都具有已知的近似体积，并且交叠以包围所述PTV(40、70、102)。
18. 如权利要求17所述的方法，还包括为预定数量的消融区域(108)评估唯一性覆 盖区(UCA)，其中，特定消融区域(108)的所述UCA是所述PTV(40、70、102)中仅由该椭圆体 覆盖的部分。
19. 如权利要求18所述的方法，还包括如果一个消融区域(108)不具有UCA，就从所 述一组候选消融区域(112)中删除该消融区域(108)。
20. 如权利要求19所述的方法，还包括确定在所述一组候选消融区域(112)中的一 个或多个消融区域(108)的中心(114)是否可以被移动靠近与所述PTV(40、70、102)有关 的中心点(104)。
21. 如权利要求20所述的方法，还包括将所述一个或多个消融区域(108)的所述中 心(114)移动到尽可能靠近所述中心点(104)的位置，同时仍保持与所述一个或多个消融 区域(108)相关的所述UCA的覆盖。
22. 如权利要求21所述的方法，还包括定义一个或多个重力点，并将至少一个消融区 域(108)拉向所述一个或多个重力点，以便在覆盖与所述至少一个消融区域(108)相关的 所述UCA的同时，进一步合并消融区域。
23. 如权利要求15所述的方法，还包括向所述用户提供指示由于所述PTV(40、70、102)的一部分极其接近所确定的关键区 域(136)，因而不能消融该部分。
24. 如权利要求23所述的方法，还包括允许所述用户将盐水注射到接近所述关键区 域(136)的空间中，以便增加在所述关键区域(136)与所述PTV(40、70、102)的所指示部分 之间的距离。
25. 如权利要求23所述的方法，还包括允许所述用户在使用具有第一尺寸的消融探 针(120)的第一消融过程期间，避开所述PTV(40、70、102)的所指示部分，并且用具有第二 尺寸的消融探针(120)消融所述PTV(40、70、102)的所指示部分，其中，所述第二尺寸小于 所述第一尺寸。
26. —种被编程执行如权利要求14所述的方法的处理器(21)或计算机介质(20)。
27. —种有助于计划用以根除病人体内的肿瘤的肿瘤消融过程的系统，包括 用于产生病人体内的肿瘤块(68)的图像的模块(16);用于接收与所述肿瘤块(68)有关的信息和与所述肿瘤块(68)的去除有关的信息的模 块(18);用于产生包围用于消融的被计划目标体积(PTV) (40、70、102)的多个体积的模块 (28)，其中，所述PTV(40、70、102)包含所述肿瘤块(68);以及用于通过确定覆盖整个所述PTV(40、70、102)的体积的最佳数量来优化所述体积的位 置设置的模块(14)。
28. —种用于计划消融过程的系统，包括优化组件(14)，从成像组件(16)输入三维PTV(40、70、102);并且 其中，所述优化组件根据所选择的消融探针(120)的已知消融体积尺寸，确定将覆盖 所述PTV(40、70、102)的消融体积(108)的一个或多个消融体积中心(114)。
29. 如权利要求28所述的系统，其中，所述优化组件(14)向图形用户界面(18)或追踪 系统之中的至少一个输出所述一个或多个消融体积中心(114)。
30. 如权利要求28所述的系统，其中，所述成像组件(18)是影像存档和传输系统 (PACS)。
全文摘要
在为消融过程制定计划时，输入被计划目标体积(PTV)(70、102)，它通常是由医生选择的，但也可以是计算机确定的。计算覆盖PTV(40、70、102)的初始的一组候选消融区域(112)。计算避开关键结构(136)的消融探针轨迹，以便为安全进入点和角度提供外科手术的5个选择。优化组件(14)确定完全覆盖PTV(40、70、102)的消融区域的最小数量，并为每一个消融产生位置和方向。可以将这些消融目标提供给追踪系统，以使得外科医生能够更准确地控制探针达到预期位置。
文档编号G06F19/00GK101795636SQ200880003120
公开日2010年8月4日 申请日期2008年1月10日 优先权日2007年1月24日
发明者D·萨默斯-斯特耶, K·I·特罗瓦托 申请人:皇家飞利浦电子股份有限公司