用于确定管状结构的拓扑支撑的评估的方法和系统及其在虚拟内窥镜检查中的使用

摘要

一种用于确定包括内壁及多个不同区域的管状结构的拓扑支撑的评估的方法，该方法包括：（a）获取表示所述管状结构的图像数据；（b）将初始种子放置在从多个不同区域之一中所选择的初始区域中；（c）执行初始区域生长，直到初始生成区域包括内壁的至少一部分和对应于多个不同区域之一的邻近区域的至少一部分为止；（d）开始包括对应于初始区域的初始树节点的树；（e）对于每个邻近区域而言：将随后的种子放置在邻近区域中；执行相应的随后区域生长，直到随后生成区域包括内壁的至少一部分和附加邻近区域的至少一部分为止；并且将对应于邻近区域的树节点添加在该树中；（f）对每个附加邻近区域执行处理步骤（e）；并且（g）根据预定拓扑参数过滤树，从而确定管状结构的拓扑支撑的评估。本发明还公开了用于评估虚拟结肠镜检查的结肠拓扑的方法的应用。

说明书

相关申请的交叉参考

本发明涉及题目为“用于过滤图像数据的方法和系统及其在虚拟内窥 镜检查中的使用”的PCT申请，其全部内容结合于此作为参考。

技术领域

本发明总体涉及图像处理，更具体地，涉及用于确定管状结构的拓扑 支撑的评估的方法和系统。本发明还涉及用于在虚拟结肠镜检查中评估结 肠拓扑的方法的应用。

背景技术

传统内视镜手术通常依靠使用插入病人的身体内的柔性光纤管来可 视地检查内部解剖结构。随后，操作员可以操作解剖结构内的管，以寻找 任何解剖结构异常。

使用该过程进行结肠镜检查虽然可靠，但都很昂贵并耗费时间。此外， 它对于患者来说是一个侵入式、不舒服且有时很痛苦的过程。

非侵入式过程（也被称作虚拟结肠镜检查）已经被用于减少侵入式结 肠镜过程的上述缺陷的至少其中一个。

这些非侵入式过程使用诸如计算机断层扫描（CT）的成像技术获取 表示解剖结构的图像数据来进行分析。

它们还包含三种类型的患者预处理程序：完全导泻准备，旨在通过使 用泻药溶液彻底清洁结肠；适度促泻准备，用于液化结肠材料并标记液体 材料的任何残留固体；最终的无促泻准备，其中，通过患者喝下的溶液来 标记结肠内的材料，诸如基于钡的制剂。

已经提出了不同的自动化技术，用于定位诸如结肠内壁的分析下的解 剖结构。然而，这些技术通常很难正确地定位结肠内壁的结构，尤其是在 其中延伸的空气区域与标记区域之间的界面附近。

事实上，如果空气区域-标记区域界面没有被正确底识别，则会导致 结肠内壁的识别及位置的泄漏，其受到了很大的关注。例如，位于结肠附 近的小肠部分会被分割，并被识别作为结肠的一部分。

此外，很差的结肠内壁分割会导致潜在结肠病变的过评估或低评估， 这个问题也非常受关注。

为了减小上述弊端，已经使用了利用用于识别结肠内壁的相应部分的 局部参数的动态算法。

例如，公布号2008/0008367所公布的美国专利申请描述了在随后设 计分割前执行能够泄漏的初始试验分割的两步分割方法。然而，这种方法 需要在界面太厚或不均匀的情况下，适当地检测出空气区域-标记区域界 面，该方法可以适当地提供结肠内壁的正确的识别和/或位置。

此外，在存在图像获取和/或阻塞性肿瘤期间，在由于患者痉挛而引 起患者结肠衰退的情况下，所述方法不能提供整个结肠的正确识别。

因此，期望提供将减小上述弊端的至少其中一个的用于确定管状结构 的拓扑支撑的评估的改进方法。

发明内容

因此，公开了一种用于确定包括内壁或多个不同区域的管状结构的拓 扑支撑的评估的方法，所述方法包括：（a）获取表示管状结构的图像数据； （b）将初始种子放置在从多个不同区域之一中选择的初始区域中；（c） 执行初始区域生长，直到在初始生成区域包括至少一部分内壁或对应于多 个不同区域之一的至少一部分邻近区域为止；（d）开始包括对应于初始区 域的初始树节点的树；（e）对于每个邻近区域而言，在邻近区域中放置随 后的种子，执行相应的随后区域生长，直到随后生成区域包括至少一部分 内壁和一部分附加邻近区域；并在树中添加对应于邻近区域的树节点； （f）对于每个附加邻近区域执行处理步骤（e）；以及（g）根据预定拓扑 阐述过滤树，从而确定管状结构的拓扑支撑的评估。

该方法提供了不依赖于分割参数的管状结构的拓扑支撑的评估，这是 很大的优势。

所获取的评估允许其需要的更好的最后处理，这也是很大的优势。这 个处理可以为用于非限定性实例的随后的分割。

此外，所获取的评估可以被用于根据体积渲染处理提供管状结构的精 确的3D表示，这是很大的优势。当然，由于不需要分割或图像数据的任 意替换，所以3D表示精确地示出了图像数据的2D信息。

此外，该方法可以被用于提供管状结构的精确的3D表示，而不依赖 于图像数据的预定精确值，这是很大的优势。因此，该方法可以使用多种 图像数据类型和多种扫描设备。

此外，在一个实施方式中，该方法提供管状结构的拓扑支撑的评估， 而不用必须使用不同区域之间的界面，这也是很大的优势。

在一个实施方式中，图像数据的获取包括从CT扫描设备接收图像数 据。

在进一步的实施方式中，图像数据的获取包括从由磁共振成像（MRI） 设备、正电子发射断层扫描（PET）设备、X射线设备、超声设备、及其 任意组合构成的组中选择的设备接收图像数据。

在一个实施方式中，从由体积医学图像数据、体积断层图像数据、及 一组平行连续图像平面构成的组中选择图像数据。

在一个实施方式中，图像数据表示解剖结构。

在一个实施方式中，图像数据包括从由像素和体像素所构成的组中所 选择的多个单位图像元素。

在一个实施方式中，多个不同区域包括多个第一物质区域和多个第二 物质区域。

在进一步的实施方式中，初始种子的放置包括从多个第一物质区域之 一中选择初始区域。初始区域生长的执行进一步包括从多个第二物质区域 之一中选择邻近区域。

在又一实施方式中，在相应随后区域生长的执行过程中，附加邻近区 域被选择，使得邻近区域和附加邻近区域的每一个都分别属于多个第一物 质区域和多个第二物质区域的相应的其中一个。

在一个实施方式中，通过定义比初始生成区域的区域生长具有更大兴 趣的区域特性的处理通过扫描图像数据的邻近部分来执行其他邻近区域 的识别。

在一个实施方式中，执行初始区域生长的执行，直到初始生成区域进 一步包括管状结构内壁外周围的至少一部分为止。

在一个实施方式中，执行相应随后区域生长的执行，直到随后生成区 域进一步包括管状结构内壁的外周的至少一部分为止。

在一个实施方式中，执行初始区域生长的执行，直到初始生成区域包 括初始区域为止。

在进一步的实施方式中，执行相应随后区域生长的执行，直到随后生 成区域包括邻近区域为止。

在一个实施方式中，区域生长的执行定义了能够进行潜在属于所属区 域单位图像元素的处理的给定直径的球形特性。

在进一步实施方式中，通过定义具有比在能够进行潜在属于所属区域 单位图像元素处理的区域生长中所包括的球形直径更大的给定直径的球 形特性的区域生长来进行潜在随后区域的识别。

在又一实施方式中，潜在随后区域的扫描识别用于区域的随后区域生 长的补充种子。

在进一步的实施方式中，补充种子的识别是以密度基标准为基础，元 素的数目定义了相同密度基标准的特性或其组合。

在一个实施方式中，随后补充种子的识别得到属于已经处理区域的种 子元素的识别，其中，仅有拓扑信息被保持并被添加至树，从而防止这样 的种子的随后区域生长。

在进一步的实施方式中，所述方法进一步包括，在初始区域生长执行 前，确定用于第一物质区域的第一物质阈值和用于第二物质区域的第二物 质阈值。

在一个实施方式中，所述方法包括：（i）在从相应区域的其中一个所 选择的补充初始区域中放置补充初始种子；（ii）执行补充初始区域生长， 直到补充初始生成区域包括内壁的至少一部分及对应于其中一个不同区 域的补充邻近区域的至少一部分；（iii）开始包括对应于补充初始区域的 初始树节点的补充树；（iv）对于每个补充邻近区域而言：在补充邻近区域 中放置补充随后种子；执行相应的补充随后区域生长，直到补充随后生成 区域包括内壁的至少一部分和补充附加邻近区域的至少一部分为止；并将 对应于补充邻近区域的树节点添加在补充树中；（v）对补充附加邻近区域 的每一个执行处理步骤（iv）；以及（vi）将补充树归组至树。

在一个实施方式中，所述方法进一步包括，在初始区域中放置初始种 子前，选择初始区域。

在进一步的实施方式中，初始区域的选择包括选择邻近管状结构端部 的初始区域。

在另一个进一步的实施方式中，所述方法进一步包括，在补充初始区 域中放置补充初始种子前，选择补充初始区域。

在一个实施方式中，补充初始区域的选择包括选择邻近管状结构其余 一个端部的补充初始区域。

在进一步的实施方式中，所述方法进一步包括使用相应数目的辅助种 子，用于将相应的树节点添加至所述树，直到所述树包括对应于初始种子 与补充初始种子的每一个的树节点之间的至少一个连续路径为止。

在另一个进一步的实施方式中，所述方法进一步包括使用相应数目的 辅助种子，用于将相应的树节点添加至所述树，直到所述树包括用于每个 不同区域的相应的一个树节点。

在一个实施方式中，由操作者手动执行初始区域的选择。

在另一个实施方式中，自动执行初始区域的选择。

在一个实施方式中，由操作者手动执行补充初始区域的选择。

在另一个实施方式中，自动执行补充初始区域的选择。

在一个实施方式中，图像数据包括多个单位图像元素，并且所述方法 进一步包括，对于对应于相应树节点的每个区域而言，确定对于区域的每 个单位元素的相应分类。

在进一步的实施方式中，多个不同区域包括多个第一物质区域和多个 第二物质区域。分类的确定包括将第一物质级赋予对应于相应的树节点的 每个第一物质区域的每个单位图像元素以及将第二物质级赋予对应于相 应的树节点的每个第二物质区域的每个单位图像元素。

在进一步的实施方式中，不属于任何区域但在随后区域的识别处理期 间被处理的每个剩余单位图像元素被彼此分组作为潜在界面型元素。

在另一个实施方式中，潜在界面型元素根据树的拓扑信息被分成两 组，这两组为非界面元素和界面型元素，界面型元素在树的两个连续节点 之间。

在又一实施方式中，所述方法进一步包括，对于每个区域生长的相应 生成区域中不包括的每个剩余单位图像元素而言：根据相应的单位图像元 素与对应于树节点的至少一个邻近区域之间的距离确定至少一个邻近参 数；并根据相应的至少一个邻近参数确定定义相应的单位图像元素与相应 级的溯源的至少一个溯源参数。所述方法进一步包括：根据邻近于对应于 两个连续节点的相应区域的相应单位图像元素的相应溯源确定所述树的 两个连续节点之间的界面类型；并根据所述树的两个连续节点之间的确定 的界面类型来确定管状结构的拓扑支撑的精确的评估。

在又一实施方式中，根据从由相应单位图像元素的密度基分布、相应 单位图像元素的分布及均匀性、感兴趣结构的拓扑信息、及界面类型的形 态参数构成的组中所选择的至少一个附加参数，进一步执行所述树的两个 连续节点之间的界面类型的确定。

在进一步的实施方式中，根据精确评估来确定管状结构的评估中线。

在一个实施方式中，根据管状结构的拓扑支撑的评估来确定管状结构 的评估中线。

在一个实施方式中，所述树的过滤包括将每个树节点彼此顺序链接。

在一个实施方式中，所述树包括至少一条主路径和闭合回路和附加支 路中的至少一条，所述树的过滤包括从所述树中去除闭合回路的至少一部 分及至少一条附加支路。

在一个实施方式中，根据与相应节点相关的每个不同区域的区域体积 来执行去除。

在一个实施方式中，管状结构包括结肠的至少一部分。

在另一个实施方式中，管状结构包括结肠的至少一部分，并且多个不 同区域包括多个空气型区域和多个标记物质型区域。

在一个实施方式中，所述方法还包括向操作者显示管状结构的拓扑支 撑的评估。

在进一步的实施方式中，显示包括在图像数据中遮蔽管状结构的周 围。

根据另一方面，还公开了对用于确定评估结肠拓扑的管状结构的拓扑 支撑的评估的方法的使用。

根据另一方面，还提供了用于根据感兴趣结构的特性或特征来处理界 面区域及所处理界面区域的随后重构的方法和系统。

在一个实施方式中，提供了用于根据结肠粘膜的特性处理界面区域和 所处理区域的随后重构的装置。

在一个实施方式中，提供了用于处理所述树的两个连续节点之间的界 面型元素的装置。

在另一个实施方式中，提供了用于在虚拟结肠镜检查中处理空气型区 域与标记型区域之间的界面型元素的装置。

在又一实施方式中，通过将值密度归于与典型的结肠粘膜元素的界面 元素不同的每个界面型元素来执行界面型元素的处理。

在一个实施方式中，通过将空气密度分配给每个界面型元素来执行界 面型元素的处理。

在另一实施方式中，在界面型元素处理之前，扩展界面型区域。

在一个实施方式中，执行界面型区域的扩展，同时维持感兴趣结构的 实际拓扑。

在另一个实施方式中，在界面型元素处理之后，根据感兴趣区域结构 的特性或特征重构所处理的界面型区域。

在进一步的实施方式中，根据在扩展处理期间所收集的每个新元素的 隐含信息（implicit information）来执行所处理的界面型区域的重构。表述 “隐含”意味着不通过直接表达的任意的暗示或理解。例如，由于粘膜为 包括未分割单位图像元素的3D区域（即，空气与组织之间的界面区域）， 所以其从未被直接表达，但是可以通过强度值的梯度场提取出诸如这个粘 膜的正向场的隐含特征。

在一个实施方式中，考虑属于每个其他邻近区域的每个元素的隐含信 息，通过将密度值归于每个界面型元素来执行处理区域的重构。

在一个实施方式中，考虑对应于围绕界面型区域的区域的所述区域的 每个元素的隐含信息，通过将密度值归于每个界面型元素来执行处理区域 的重构。

在又一实施方式中，通过考虑从邻近区域元素的密度所确定的法向矢 量所获取的相对空间强度的隐含信息来执行将密度值归于每个界面型元 素。

在一个实施方式中，隐含信息为梯度矢量场、矢量场、及密度基矢量 场的其中一个或组合。

在一个实施方式中，隐含信息为标量场（scalar field）、矢量场、及张 量场的其中一个或组合。

根据另一方面，也提供，用于确定管状结构的拓扑支撑的评估的系统。 该系统包括：数据接收单元，用于接收表示管状结构的图像数据；放置单 元，可操作地耦接至数据接收单元，用于在每个相应的区域中放置每个种 子；处理单元，可操作地耦接至放置单元，用于执行每个区域生长；树建 立单元，可操作地耦接至处理单元，用于建立树；以及过滤单元，可操作 地耦接至树建立单元，用于根据预定的拓扑参数过滤树，从而确定管状结 构的拓扑支撑的评估。

在一个实施方式中，该系统进一步包括显示单元，可操作地耦接至过 滤单元，用于显示管状结构的拓扑支撑的评估。

根据另一方面，还提供其上记录指令的机器可读介质，指令用于执行 确定管状结构的拓扑支撑的评估的方法。

根据另一方面，还提供了对根据先前所述方法确定管状结构的拓扑支 撑的评估进行交易的方法，其中，确定管状结构的拓扑支撑的评估是收费 的。

根据另一方面，还提供了对确定管状结构的拓扑支撑的评估进行交易 的方法，该方法包括：接收图像数据；执行先前所述的方法；通过收费提 供管状结构的拓扑支撑的评估。

根据另一方面，还提供了对确定管状结构的拓扑支撑的评估进行交易 的方法，该方法包括：由供应商将如先前所述的用于确定管状结构的拓扑 支撑的评估的系统提供给第三方；操作该系统，其中，由第三方通过收费 进行操作；并由第三方将收取的费用的至少一部分返还给供应商。

可以通过几种类型的图像数据使用管状结构的拓扑支撑的评估的确 定方法，这具有很大的优势。

此外，在虚拟结肠镜检查的情况下，例如，即使在患者的结肠由于图 像获取期间的痉挛和/或存在阻塞性肿瘤而衰退的情况下，也可以提供适当 的结肠拓扑支撑的评估，这具有很大的优势。

此外，在虚拟结肠镜检查的情况下，由于2D穿越可视化可以被快速 提供给操作者或医生，所以该方法可以快速提供结肠的评估中线，这具有 很大的优势。

表达“区域”指的是在同一小块区域中彼此全部连续的一组邻近单位 图像元素。依赖于所使用的图像数据，区域可以是2D或3D的。

表达“管状结构”应该被理解为具有至少两端的任何中空延长结构。

附图说明

为了容易理解本发明，在附图中利用实例图解说明了本发明的实施方 式。

图1是示出了根据本发明一个实施方式的用于确定管状结构的拓扑 支撑的评估的方法的流程图。

图2是示出了根据本发明的用于确定管状结构的拓扑支撑的评估的 方法的另一实施方式的流程图。

图3是示出了根据本发明的用于确定管状结构的拓扑支撑的评估的 方法的另一实施方式的流程图。

图4是根据本发明的用于确定管状结构的拓扑支撑的评估的系统的 一个实施方式的框图。

图5A示出了表示结肠的图像数据的一部分。

图5B示出了表示图5A的图像数据的一部分的3D图

图5C示出了患者结肠的3D表示。

图6是示出了可以实现用于确定拓扑支撑的评估的方法的处理设备 的实施方式的框图。

图7A是示出了管状结构的实施方式的示意图。

图7B是示出了对应于图7A所示的管状结构的树的示意图。

图8A是示出了管状结构的另一个实施方式的示意图。

图8B是示出了对应于图8A所示的管状结构的树的示意图。

图9A是示出了管状结构的另一个实施方式的示意图。

图9B是示出了对应于图9A所示的管状结构的树的示意图。

图10示出了根据一个实施方式如何执行其他邻近区域的识别。

图11至图15示出了根据实施方式的结肠粘膜的重构。

通过下面所包括的详细描述，本发明的另外的细节及其优势将变得显 而易见。

具体实施方式

在实施方式的下面描述中，通过发明可实行的实例的图解说明来进行 对附图的参考。可以理解，在不背离所公开的发明的范围的情况下，可以 制造或使用各种其他实施方式。

本发明涉及用于确定管状结构的拓扑支撑的评估的方法和系统，特别 用于医学图像处理领域。纵观本描述，将就虚拟结肠镜检查中评估结肠拓 扑的具体应用方面描述所述方法，但相关技术人员将理解，该方法不限于 该具体应用，并且可以考虑许多其他应用，当阅读本描述时，这点将变得 很显而易见。就非限定性实例而言，本方法可以用在CT肠动扫描应用、 用于检测大动脉腹部动脉瘤的应用以及用于肺和脑动脉瘤虚拟内窥镜检 查应用的虚拟内窥镜应用。

相关技术人员将理解，例如，本发明用于确定管状结构的拓扑支撑的 评估的方法通常用于有利于诸如结肠癌筛查的解剖结构的随后检查。相关 技术人员也将理解，该方法也适用于包括至少两个相的解剖结构，例如， 包括内壁及在其中延伸的多个空气区域和标记排泄物区域的结肠结构。

由于该方法不限于特定类型的图像数据，所以其特别具有优势。更适 当地，该方法可以被用在不同类型的图像数据组中，这点将在下面变得更 加显而易见。此外，可以既不依赖于图像数据的严格预定值也不依赖于特 定对比度来实现该方法，这点对于相关技术人员也将变得很显而易见。

当然，相关技术人员将会理解，上面所述的系统和方法由于可以通过 作为非限定性实例的无准备（prepless）CT结肠镜检查、无泻药的CT结 肠镜检查、通过标记药剂的温和制剂CT结肠镜检查、以及通过标记用于 CT结肠镜检查的残余液体/大便的泻药制剂来使用，而特别具有优势。

对比通过低纤维饮食和大便标记所增强的例行程序和非准备 （unprepped）CT结肠镜检查来描述无准备CT结肠镜检查：一项试验性 研究，Abraham H.Dachman and al，Abdom Imaging（2007）32:96-104； 无泻剂的CT结肠成像：可行性研究，Matthew R.Callstrom，Radiology 2001；219:693-698，并且还在无泻剂肠道制剂的CT结肠成像情况的结肠 肿瘤的CAD中，Marius George Linguraru等，30th annual international IEEE  EMBS Conference Vancouver，British Columbia，Canada,Auguest 20-24,2008。这些参考的全部内容结合于此作为参考。

在使用虚拟结肠镜检查的无泻药大肠癌筛查试验研发中描述了无泻 药CT结肠镜检查：可行性研究，Kristina T.Johnson,AJR：188，2007年1 月，p2936；CT结肠成像前作为清洁方法的膳食粪便标记：初步结果-Polyp  Detection and Patient Acceptance 1，Philippe A.Lefere，Radiology 2002；224： 393-403；并且在通过粪便标记的无泻药CT结肠成像中：有无电子大便减 少的性能，C.Daniel Johnson,AJR：190，2008年2月，P361-366。这些参 考的内容都结合于此作为参考。

在Sebastiaan Jensch and al，AJR：191，2008年7月，P158-167具有 不同量的用于CT结肠镜检查的温和泻药的四种治疗方案（four regimen） 的图像质量和患者接受中描述了通过标记制剂的温和制剂CT结肠成像， 其全部内容结合于此作为参考。

在用于CT结肠镜检查的钡基粪便标记的功效中描述了通过用于CT 结肠镜检查的残余液体/大便的标记的泻药制剂：在韩国人口中使用高低密 度钡悬浮液之间的比较-初步的研究，Min Ju Kim等，Korean J Radiol 10(1)， 2009年二月，P25-33；在替代中：粪便标记，Philippe Lefere核 Stefaan Gryspeerdt，虚拟结肠镜检查，斯普林格柏林海德堡，2006年， P35-49；以及在基于标记中，电子清洁CT结肠成像：患者的舒适度和图 像的可读性，Michael E.Zalis等，Radiology：volume 239：1number，2006 年4月，P149-159。这些参考的每一个都结合于此作为参考。

相关技术人员将理解，无泻药制剂可以涉及使用可具有潜在泻药副作 用并且可以提供比仅标记制剂更好的残留标记的碘。

此外，相关技术人员将理解，所公开的方法能够依赖于所使用的处理 源来提供相对较快的周转时间的管状结构的评估。

典型地，如在CT结肠镜检查中精确快速3D结肠分割，Dongqing Chen, Rachid Fahmi,Aly A.Farag,Robert L.Falk,and Gerald W.Dryden,ISBI 2009 p490-493中所提及的一样，相关方法在5分钟至18分钟之间采用每个数 据组，并且对于在奔腾IV 2.6GHz计算机上的512×512×440的单个CT 扫描而言，本方法将在3分钟至5分钟内执行感兴趣区域的拓扑定义，用 于完成结肠癌筛查研究，其为可对比维度的两个数据组并准备用于通过体 积渲染的可视3D检查。因此，本方法至少是当前相关技术方法两倍快。

图5A至图5C示出了表示管状结构的图像数据组的图像500、图解 说明情况下的结肠、及评估的结肠拓扑502的实例。

图7A示出了包括内壁702和多个不同区域的管状结构700的实例。 在管状结构700包括结肠或结肠的至少一部分情况下，内壁702可以包括 结肠粘膜以及软脂肪组织。在图解说明的情况下，多个不同区域包括多个 第一物质区域704（也被称作空气型区域）和对应于标记排泄物区域的多 个第二物质区域706（也被称作标记物质型区域）。

纵观本描述，表达“区域”指的是在同一小块中全部彼此连续的一组 邻近单位图像元素。依赖于所使用的图像数据，区域可以为2D或3D。

参照图1，示出了根据一个实施方式的用于确定管状结构的拓扑支撑 的评估的方法的流程图。

当阅读本描述时将变得很显而易见的是，用于确定管状结构的拓扑支 撑的评估的方法依赖于：用于建立表示管状结构的树的给定步骤；所述树 包括表示在管状结构中延伸的空气区域和标记物质区域的连续的连续节 点。

现在，将参照图1、图7A及图7B描述该方法的实施方式。

根据处理步骤100，提供了表示管状结构的图像数据。作为非限定性 实例，图像数据可以包括相关技术中众所周知的体积医学图像、体积断层 图像、和/或多个平行连续图像平面。

在一个实施方式中，处理步骤100包括从CT扫描设备接收图像数据。 在另一个实施方式中，可以从磁共振成像（MRI）设备、正电子发射断层 扫描（PET）设备、X射线设备、超声设备、或这些设备的任意组合中接 收图像数据。在另一个实施方式中，可以从数据库中检索图像数据，甚至 可以从诸如压缩光盘的可读介质或例如图片归档及通信系统（PACS）中 检索图像数据。

在一个实施方式中，图像数据包括多个单位图像元素，例如，像素或 体像素。然而，相关技术人员将理解，表达“单位图像元素”不限于像素 或体像素，而应该被理解为包括任何同质元素、图像或显示面上的点或圆 点、几何元素、网络网格、网格面、或者强度、颜色、或其他参数与其他 强度、颜色或其他参数分别相关联的顶点。

再次参照图1，根据处理步骤110，初始种子被放置在从多个不同区 域之一中选择的初始区域中。

在图7A所示的示例性实施方式中，初始种子Xc被放置在初始区域 A1中。

在优选实施方式中，在初始种子被放置在初始区域之前，首先选择初 始区域。在一个实施方式中，由操作者手动执行初始区域的选择。可替换 地，在另一个实施方式中，自动执行初始区域的选择。相关技术人员将理 解，可以根据各种参数执行初始区域的自动选择，例如，在Iordanescu G. Pickhardt PJ,Choi JR,Summers RM，Automated Seed placement for colon  segmentation in computed tomography，Acad Radiol.，2005 Feb；12（2）： 182-90中所描述的那样。

在优选实施方式中，初始区域被选择邻居管状结构的端部，其将在下 面变得显而易见。在图7A所示的示例性实施方式中，被选择的初始区域 靠近结肠的盲肠延伸。

再次参照图1，根据处理步骤120，执行初始区域生长，直到初始生 成区域包括至少一部分内壁和对应于多个不同区域之一的邻近区域的至 少一部分。

在优选实施方式中，从第一物质区域之一中选择初始区域，同时从第 二物质区域之一中选择邻近区域。换句话说，如后文将被更详细描述的那 样，从两种类型区域之一中替换地选择所考虑的区域。

在优选实施方式中，执行初始区域生长，直到初始生成区域进一步包 括管状结构的内壁的外周的至少一部分。换句话说，选择初始区域，随后， 还选择紧邻区域，直到生成区域也包括管状结构的内壁的一部分外周。

相关技术人员将理解，在管状结构包括结肠或其一部分的情况下，管 状结构的内壁的外周可以包括软脂肪组织、肌肉、骨骼、或其他邻近结构 的一部分，例如，小肠的一部分。

在进一步实施方式中，执行初始区域生长，直到初始生成区域包括整 个初始区域。

相关技术人员将理解，在区域生长的现有应用中，以均匀性为基础根 据迭代处理对同一区域的单位图像元素分组，从而将图像数据的选择部分 分割成感兴趣的不同区域。在这些先前技术应用中，区域生长被用于从其 他区域中提取出一个特殊区域。

例如，在号为US 2002/0193687并且题目为虚拟内窥镜检查的自动分 析的美国专利申请公布中，区域生长被解释如下：使用三维区域生长技术 和初始静态阈值来分割感兴趣的区域。所选择的阈值应接近通过包括周围 结构作为感兴趣区域的部分来选择并且无分割程序失败的最大阈值。

相关技术人员当阅读本说明书时将理解，在本应用中，区域生长不被 用于提取特殊区域而将图像分割成不同的感兴趣区域的目的。相反，区域 生长被用于选择相应种子周围的图像部分。如先前所述，通过区域生长所 获取的生成区域应包括几种类型的区域部分。

相关技术人员也将理解，可以对体积图像数据执行紧邻的选择。因此， 应该理解，生成区域可以为针对多个连续的二维图像所获取的三维体积。

再次参照图1，根据处理步骤130，开始包括对应于初始区域的初始 树节点的树708.

在图7A和图7B所示的示例性实施方式中，以初始树节点A1开始 树。

再次参照图1，并如下面更详细描述一样，对在处理步骤120中所发 现的每个邻近区域执行处理步骤140、150、及160。

当然，根据处理步骤140，随后的种子被放置在先前所发现的邻近区 域中。

根据处理步骤150，执行相应的随后区域生长，直到随后生成区域包 括内壁的至少一部分和附加邻近区域的至少一部分。

在优选实施方式中，在处理步骤150中，附加邻近区域被选择为使得 邻近区域和附加邻近区域的每一个都分别属于多个第一物质区域和多个 第二物质区域中的相应的一个。换句话说，如前所述，所考虑的区域通过 其相应类型被替换选择。

在优选实施方式中，执行相应的随后区域生长，直到随后生成区域进 一步包括管状结构的内壁的外周的至少一部分。

在进一步的优选实施方式中，执行相应的随后区域生长，直到随后生 成区域包括整个邻近区域。

在一个实施方式中，如图10所示，通过定义比初始生成区域的区域 生长更受关注的区域的特性的处理，通过扫描图像数据的邻近部分来执行 附加邻近区域的识别。相关技术人员将理解，可以使用区域生长或光线投 射（raycast）处理。

根据处理步骤160，将对应于邻近区域的树节点添加至树。

再次参照图1，对在处理步骤150中所发现的每个附加邻近区域执行 处理步骤140、150及160.

在一个实施方式中，如图10所示，区域生长的执行展示了能够处理 潜在属于该区域的单位图像元素的给定参数的球。

在进一步的实施方式中，通过展示具有比在区域生长（能够处理潜在 属于给定区域的单位图像元素）中所包括的球直径更大的给定参数的球的 区域生长，来进行潜在的随后区域的识别。

在又一实施方式中，潜在的随后区域的扫描识别用于区域的随后区域 生长的补充种子。

在进一步的实施方式中，补充种子的识别以基于密度的标准、以相同 密度基标准为特性的元素数、或其组合为基础，由于其能够阻止考虑人造 元素，所以具有很大的优势。

在一个实施方式中，随后补充种子的识别引起属于已经处理区域的种 子元素的识别，在已经处理的区域中，仅拓扑信息被保留并被添加至树， 从而防止这些种子的随后区域生长。

在图7A和图7B所图解说明的示例性实施方式中，一旦初始树节点 A₁已经被添加在树中，则相应种子被放置在区域A₁的邻近区域中，即， 区域T_alpha、T_beta、T₁、T₂、及T_gamma并且相应的树节点T_alpha、T_beta、T₁、 T₂、及T_gamma被添加至树。

随后，如果需要，则为了发现还没有被考虑的其他邻近区域，反过来 考虑这些上述区域的每一个。例如，一旦树节点T₁已经被添加在树中， 则发现对应于相邻小肠部分的区域SB₁，并且相应的节点SB₁被添加在树 中。

一旦每个发现的邻近区域都已经经受了处理步骤140、150、及160， 则根据处理步骤170，根据预定拓扑参数来过滤树，从而确定管状结构的 拓扑支撑的评估。

在一个实施方式中，树的过滤包括将每个树节点一个接一个地顺次链 接。

在进一步的实施方式中，如下面将详细说明的一样，可以从树中消除 与具有预定值以下的面积和体积的区域相对应的死支路和节点。死支路被 定义为不能被用于在管状结构的两端之间提供连续路径的树的一部分，这 点将在下面显而易见。

仍在进一步的实施方式中，也可以从树中去除肯定属于一结构的树的 部分，该结构不属于感兴趣的管状结构（诸如骨骼结构）。当然，对应于 骨骼部分的节点可以已经包括在树中。然而，相关技术人员将理解，可以 执行骨骼去除，如在Automatic vessel extraction by patient motion correction  and bone removal in brain CT angiography,Helen Hong and al，international  Congress Series Volume 1281,May 2005,Pages 369-374所描述的一样。相关 技术人员也将理解，可以考虑其他各种方法。

在一个实施方式中，预定拓扑参数可以以管状结构的众所周知的拓扑 为基础。例如，管状结构包括其两端之间的连续路径，并且管状结构不包 括回路部分。

在这点上，可以获取管状结构的拓扑支撑的粗略评估。

在图7A和图7B所示的实施方式中，一旦节点已经被一个接一个地 链接，则单个的连续路径在直肠与盲肠之间延伸。因此，相关技术人员将 理解，这个路径可以表示管状结构的拓扑支撑的粗略评估。此外，由于树 的其他支路可以不被用于发现连续路径，所以它们可以从树中被随意去 除。

在一个实施方式中，如后文将要详细描述的，一旦已经确定了管状结 构的拓扑支撑的评估，则随后可以确定管状结构的评估中线。

在一个实施方式中，连续路径的每个区域的中心可以被用于粗略地评 估中线。然而，相关技术人员将理解，各种其他手段可以被用于提供中线 的粗略评估。例如，如在骨架化及其应用，Kálmán Palágyi,Dept.Image  Processing & Computer Graphics University of Szeged,Hungary,Summer  School on Image Processing SSIP 2009中所描述的那样，可以使用区域的拓 扑细化或通过沃罗努瓦图从区域的粗略分割中提取的中线。如相关技术人 员众所周知的那样，也可以使用包含水平组处理或基于距离的骨架化的其 他方法。

在进一步的实施方式中，如随后也将被详细描述的那样，管状结构的 拓扑支撑的评估随后可以被显示给操作者。

现在，参照图2，在优选实施方式中，可以使用第二初始种子。

当然，如图7A所示，为了确保已经考虑了全部管状结构，在优选实 施方式中，有利地使用两个初始种子，靠近管状结构的第一端的第一种子， 并且靠近管状结构的第二端的第二种子。在管状结构为结肠的情况下，两 个初始种子被放置在盲肠和直肠附近。由此，发现对应于两个初始种子的 两个节点之间的连续路径可以确保已经完整考虑了管状结构，这具有很大 的优势。

在图7A所示的示例性实施方式中，第二初始种子Xr被放置在直肠 中。由于在对应于两个初始种子的两个节点之间存在连续路径，所以这个 连续路径可以代表整个管状结构的拓扑支撑，如将在下面详细描述的那 样。

因此，仍参照图2，并根据处理步骤200，补充初始种子被放置在从 相应区域之一中所选择的补充初始区域中。

相关技术人员将理解，在一个实施方式中，有利的是，与初始种子相 关联的处理步骤平行地执行与补充初始种子相关联的处理步骤。可替换 地，这两个处理可以被顺序地执行。

在优选实施方式中，在补充初始种子被放置在补充初始区域之前，类 似于处理步骤110的初始区域，首选选择补充初始区域。在一个实施方式 中，由操作者手动执行补充初始区域的选择。可替换地，在另一个实施方 式中，自动执行补充初始区域的选择。

在虚拟结肠镜检查的情况下，在一个实施方式中，可以通过根据图像 数据的全局直方图来评估空气型元素的体积并考虑人体形态学（例如，当 面向患者时，降结肠在患者的右侧腹部，并对应于拉长的空气包）来执行 任意补充种子的自动选择。同时，可以根据臀部形态学并查找邻近于直肠 的显著空气包来识别乙状结肠（在多数情况下，位于图像数据组的“底部” 位置）。这两种方法不是限定性的，并且任何领域的相关技术人员将理解， 在不背离本发明的的范围的情况下，这种形态学方法可以有很多种。

相关技术人员将理解，在进一步的实施方式中，可以同时使用多个补 充种子，同时由此建立每个相应的树。相关技术人员也将理解，在一个实 施方式中，为了不选择已包括在其中一个树中的区域，可以使用补充种子 的选择规则。

在优选实施方式中，如先前所述，被选择的初始区域靠近管状结构的 一端延伸。仍在优选实施方式中，所选择的补充初始区域靠近管状结构的 另一端延伸。

如前所述，在图7A所示的实例中，所选择的初始区域靠近结肠的盲 肠延伸，同时所选择的补充初始区域靠近结肠的直肠延伸。

根据处理步骤210，执行补充初始区域生长，直到补充初始生成区域 包括内壁的至少一部分和与多个不同区域之一对应的补充邻近区域的至 少一部分。

根据处理步骤220，开始包括对应于补充初始区域的初始树节点的补 充树。

再次参照图2并如下面更详细描述的那样，对在处理步骤220中所发 现的每个补充邻近区域执行处理步骤230、240、及250。

当然，根据处理步骤230，补充随后种子被放置在相应的补充邻近区 域中。

根据处理步骤240，执行相应的补充随后区域生长，直到补充随后生 成区域包括内壁的至少一部分和补充附加邻近区域的至少一部分。

根据处理步骤250，对应于补充邻近区域的树节点被添加至补充树 中。

再次参考图2，针对在处理步骤240中所发现的每个补充附加邻近区 域执行处理步骤230、240、及250。

一旦所发现的每个补充邻近区域都经受了处理步骤230、240、及250， 则根据处理步骤260，补充树与树被分组。

在一个实施方式中，补充树和树的分组可以包括在对应于相同区域的 每个树中发现至少一个节点并根据这个共用节点合并两个树。然而，相关 技术人员将理解，如下面将进行详细描述的，可以考虑将树一起分组的各 种其他处理程序。

根据处理步骤270，随后可以根据预定拓扑参数来过滤树，从而确定 管状结构的拓扑支撑的评估。

相关技术人员将理解，对于结肠拓扑支撑确定的具体应用而言，可以 使用至少一个辅助种子。当然，在某些情况下，例如，由于在图像获取期 间的神经性痉挛和/或阻塞性肿瘤的存在，所以患者检查下的结肠可能会衰 退。当结肠衰退时，其内壁可能阻塞盲肠与直肠之间的通道。随后，结肠 内的体积被分割成多个基于管状的部分，并且上文参照图1所描述的方法 不会提供用于结肠整体长度的所评估的拓扑支撑，下文中将被详细描述。 在这种情况下，相应数目的辅助种子可以被用于将相应树节点添加至树， 直到树包括对应于每个初始种子和补充初始种子的树节点之间的至少一 个连续通路为止。

在一个实施方式中，对于每个辅助种子建立相应的辅助树，随后，辅 助树被分组至树。

相关技术人员也将理解，在管状结构包括多余两个的端部的情况下， 可以使用这种辅助种子。在这种情况下，为了确保已经考虑到全部的结构， 可以有利地使用用于管状结构的每端的相应初始种子。

在虚拟结肠镜检查的情况下，在一个实施方式中，结肠的形态学参数 可以被用于放置辅助种子。例如，由于结肠为连续伸长结构，辅助种子可 以被放置在树或辅助树的空间延长部分。

如先前所述，在这点上，可以获取结肠的拓扑支撑的粗略评估。相关 技术人员将理解，不用分割任何区域，就已经获取了这个粗略估计，其具 有很大的优势，将在下面变得更加显而易见。

相关技术人员将理解，这个粗略的估计可以被用于提供管状结构的粗 略中线，其具有很大的优势。当然，根据这个粗略中线，可以提供2D穿 越可视化。换句话说，可以根据所获取的粗略中线屏蔽不感兴趣的全部图 像部分。由于操作者可以以更方便的方式检查图像数据，不会受不感兴趣 的图像的其他区域干扰，所以具有很大的优势。由于其可以大大加快操作 者检查图像的速度同时减少检查图像所通常需要的培训，所以其具有特别 大的优势。对于会被相同图像的不同部分中的残余标记液体或由于仅命名 了一些图像不同部分的主要病变存在而被干扰的读者而言，这种方法具有 很大的价值。

相关技术人员也将理解，也可以提供使用体积渲染处理的3D穿越可 视化。

如先前所说明的一样，已经获取的管状结构的拓扑支撑的粗略评估以 结肠内存在空气型区域和标记型区域的交替并在盲肠与直肠之间延伸以 及这样的交替帮助了确定拓扑支撑的事实为基础。

因此，在上述方法中，空气型区域与标记型区域之间的界面和空气型 区域或标记型区域与结肠组织之间的界面已经不被使用或不被考虑用于 提供粗略评估。

相关技术人员将理解，典型地，在这些界面处，主要由于不同类型区 域可以包括未标识或不均质标识的排泄物，所以它们的识别会很困难。

在所述方法的优选实施方式中，可以进一步获取结肠拓扑支撑的精确 评估。通过定义每个感兴趣界面（即，对应于树的两个连续节点之间延伸 的支路的界面）的属性或分类获取这种精确评估。实际上，如此后详细描 述的那样，一旦感兴趣界面已经被定义为空气型区域-标记型区域界面或 具有结肠组织的界面，则可以改进树的过滤。

现在，参照图3，现在将描述方法的进一步的实施方式。在这个实施 方式中，可以确定管状结构的拓扑支撑的精确评估。如下文中将被相关技 术人员理解一样，在这个实施方式中，各种类型界面可以被识别和分类， 进一步精确拓扑支撑的评估。

因此，根据处理步骤300，一旦已经获取了树，则用于对应于相应树 节点的每个区域的每个单位元素的相应分类被确定。换句话说，一旦在处 理步骤120、150、210、或240期间一个区域已经被确定为第一物质型区 域或第二物质型区域，则这个区域的每个单位图像元素被视为属于这种类 型的区域。

在一个实施方式中，类别的确定包括为对应于相应树节点的每个第一 物质区域的每个单位图像元素赋予第一物质类别并为对应于相应树节点 的每个第二物质区域的每个单位元素赋予第二物质类别。在结肠拓扑评估 的具体应用中，相应单位图像元素已经被“分类”为属于空气型区域或标 记型区域。

在这点上，还没有被考虑的剩余单位元素属于空气区域-标记区域界 面或与结肠组织连接的界面。

为了提供在树的两个连续节点之间延伸的界面类型的识别，可以执行 进一步的处理步骤。一旦感兴趣的界面被适当识别，则可以提供拓扑支撑 的精确评估，当阅读随后的描述时，对于相关技术人员将变得更加显而易 见。

再次参照图3，为了识别每个感兴趣界面的类型，如下面详细描述一 样，对于在每个区域生长的相应生成区域中所包括的每个剩余单位图像元 素可以执行处理步骤310和320。

根据处理步骤310，根据相应的单位图像元素与对应于树节点的至少 一个邻近区域之间的距离，确定每个感兴趣的单位图像元素的至少一个接 近参数。相关技术人员将理解，由于不依赖于分割处理或定量参数，而是 依赖于定性参数，所以使用接近参数具有优势。

在一个实施方式中，为了考虑感兴趣的单位图像元素与在其周围延伸 的每个区域之间的每个距离，可以使用多个接近参数。

根据处理步骤320，将相应单位图像元素的联系定义为相应分类的至 少一个联系参数（affiliation parameter）根据相应的至少一个接近参数被确 定。相关技术人员将理解，由于其不会替换或修改图像数据，所以联系参 数的确定具有优势。

在一个实施方式中，可以根据各种附加参数确定单位图像元素的联系 关系。例如，可以考虑围绕所选择的单位图像元素的邻近区域的形态和形 状以及管状结构的整体拓扑。在进一步的实施方式中，这种形态参数的其 中一个将充分利用身体语言的事实，残余液体将如任意其他液体一样具有 使其表面平坦的趋势。认识到大的平坦区域不可能描述诸如内腔的中空器 官，可以为大的微分器（differentiator）来将标记/空气界面（水平面为残 余标记液体的上部）特征化。

对每个感兴趣的单位图像元素执行处理步骤310和320，及，这些便 于感兴趣界面的识别。

一旦每个感兴趣的单位图像元素已经经受了处理步骤310和320，则 根据处理步骤330，根据邻近对应于两个连续节点的相应区域的相应单位 图像元素的相应联系来确定树的两个连续节点之间的界面类型。

在一个实施方式中，根据从由相应单位图像元素的密度基分布、相应 单位图像元素的分布基均匀性、感兴趣结构的拓扑信息、及界面类型的形 态参数构成的组中所选择的至少一个附加参数进一步执行两个连续节点 之间的界面类型的确定。相关技术人员将理解，表达“基于密度的分布” 应该被理解为包括密度分布、其衍生形式、以及其任意的组合。类似地， 表达“基于分布的均匀性”应该被理解为包括分布均匀性、其衍生形式、 以及其任意组合。

当然，感兴趣界面的单位图像元素的特定密度分布和分布均匀性可以 被用于确定界面类型，如发布号为WO/2007/048091、题目为用于计算断 层结肠镜检查图像的数字清洁的结构分析系统、方法、软件排列、及计算 机可访问介质的PCT申请中所描述的那样。值得注意的是，在PCT申请 中没有讨论管状结构的接近参数和拓扑知识。

另外，也应该考虑感兴趣界面的形态。例如，可以考虑感兴趣界面的 厚度和体积。此外，诸如平面形状或展开形状的界面形状也特别有用。

在一个实施方式中，不属于任意区域但在随后区域的识别处理期间被 处理的每个剩余单位图像元素与其他单位图像元素一起被分组为潜在的 界面型元素。

在另一个实施方式中，根据树的拓扑信息将潜在界面类型元素分成两 组，这两组为非界面元素和界面型元素，界面型元素在树的两个连续节点 之间。

再次参照图3，根据处理步骤340，根据树的两个连续节点之间确定 的界面类型来确定管状结构的拓扑支撑的精确评估。

相关技术人员将理解，在先前所描述的实施方式中，由于树节点之间 的界面类型而增强了树的过滤。

在一个实施方式中，如前所述，树的过滤可以包括一个接一个地顺次 链接每个树节点。

随后，树通常包括至少一条主路径及最终的闭合回路和附加支路的至 少一个。

在一个实施方式中，通过去除闭合回路的一部分和至少一个附加支路 的至少之一（通常为死支路）来过滤树。

在优选实施方式中，从树中去除不能被用于提供管状结构的两个末端 之间的连续通路的死支路（若有的话）。在进一步的实施方式中，通过去 除其第一节点与具有低于定义值的区域体积的区域相关的闭合回路（若有 的话）来进一步过滤树。

换句话说，可以从树中去除更小的区域，通常为与结肠内壁反向延伸 而不会占据结肠整个截面的标记区域。当然，如果小标记区域为闭合回路 的一部分的第一节点，则其可以为对应于结肠外部区域的树部分的一部 分，例如，小肠部或骨骼。如图8A~图9B中所图解说明及如下所详细说 明的一样，其也可能为在结肠中延伸的两个区域之间不正确的捷径。

当然，在图7B中，如先前所述，由于包括节点SB1的支路对应于死 支路，所以其可以被去除。

在图8B中，如图8A所图解说明的一样，由于相应的界面T1-A1应 该被识别为在软组织中延伸，所以可以去除包括节点T1的闭合部分。

在图9A和图9B中，相关技术人员将理解，界面的分类可以帮助选 择区域A1和T1之间的便捷路径。

在优选实施方式中，根据拓扑支撑的精确评估确定管状结构的精确评 估中线。

相关技术人员将理解，由于管状结构的中线可以被充分仿射从而能够 便于结构观察，所以这个实施方式特别受关注。

在一个实施方式中，精确评估中线可以与图像数据一起被提供给用于 3D可视化的体积渲染引擎。在进一步的实施方式中，图像数据可以经受 先前的电子清洁处理程序，用于去除标记物质区域。在通过相同申请人联 合申请、名为“在虚拟内窥镜检查中用于过滤图像数据的方法和系统及其 使用”的PCT申请中公开了这种电子清洁处理程序，该说明书结合于此 作为参考。

相关技术人员将理解，在一个实施方式中，根据所提供的图像数据的 类型来执行初始区域生长，可以确定用于第一物质区域的第一物质阈值和 用于第二物质区域的第二物质阈值。

根据另一方面，一旦已经执行了结构的拓扑，则如图11至图15所图 解说明的一样，可以处理界面区域，并且可以执行根据感兴趣结构特性或 特征的处理界面区域的随后重构。

在一个实施方式中，提供用于处理界面区域和和根据结肠粘膜的特性 的处理区域的随后重构的装置。

在一个实施方式中，提供了用于处理树的两个连续节点之间的界面类 型元素的装置。

在另一个实施方式中，提供了虚拟结肠镜检查中用于处理空气型区域 与标记型区域之间的界面型元素的装置。

在又一实施方式中，通过将值密度归于与典型结肠粘膜元素不同的每 个界面型元素，来执行界面型元素的处理。例如，在一个实施方式中，通 过将空气密度归于每个界面型元素来执行界面型元素的处理。相关技术人 员将理解，可以选择其他值，诸如低于空气型区域值的值。

在进一步的实施方式中，在界面型元素的处理前，界面型区域被扩展。

在一个实施方式中，执行界面型区域的扩展，同时维持感兴趣结构的 实际拓扑。

在另一个实施方式中，界面型元素处理后，进行根据感兴趣结构的特 性或特征的处理的界面型区域的重构。

在进一步的实施方式中，根据在扩展处理期间所收集的每个新元素的 隐含信息来执行所处理界面型区域的重构。值得提出，表达“隐含”意味 着通过非直接表达的任意暗示或理解。例如，由于粘膜为包括未被分割的 单位图像元素（即，空气与组织之间的界面区域）的3D区域，所以其从 未被直接表达，但是可以通过强度值的梯度场提取出隐含特征，诸如这个 粘膜的正向场。

在一个实施方式中，考虑属于每个其他邻近区域的每个元素的隐含信 息，通过将密度值归于每个界面型元素来执行所处理区域的重构。

在一个实施方式中，考虑对应于界面型区域周围的区域的多个区域的 每个元素的隐含信息，通过将密度值归于每个界面型元素来执行所处理区 域的重构。

在另一实施方式中，通过考虑通过邻近区域元素的密度所确定的正向 矢量所获取的相对空间强度的隐含信息来执行密度值对于每个界面型元 素的归因。

在一个实施方式中，隐含信息为梯度矢量场、矢量场、密度基矢量场 的其中一个或组合。在另一个实施方式中，隐含信息为标量场、矢量场、 及张量场的其中一个或组合。

再次参照图11至图15，将详细描述根据隐含信息的处理。在一个实 施方式中，可以执行给定界面的邻近区域与通过界面所连接的两个不同区 域的邻近区域的交叉隔离。界面邻近区域的这个隔离部被称作重构支撑。

隐含几何特征可以被提取，提供关于界面区域的非常限定支撑行为的 信息用于重构。通用特征为根据图像元素强度梯度的隐含正向量场。相关 技术人员将理解，这种隐含信息允许界面型区域中每个单位图像元素的重 构，但是特别对于虚拟结肠镜检查的特殊情况下结肠粘膜周围的区域感兴 趣。

每个界面元素属于目标区域（空气）的可能性可以被计算。为了将给 出的隐含特征的连续性最大化，计算这种可能性或联系。在正向量场的情 况下，每个体像素的联系以通过用于重构的给定体像素的中心及支撑邻近 体像素的中心所形成的矢量；以及在支撑邻近体像素的隐含正向矢量的投 影强度为基础。

界面区域的强度值根据它们的联系在目标区域与支撑区域的分布之 间被重构。

图15示出了为了防止由于感兴趣的精确拓扑的保存而在结肠中产生 切割潜在皱褶可以如何处理邻近界面，其在经常切割皱褶的先前相关技术 方面具有很大的优势。

根据另一方面，还提供了用于提供管状结构的拓扑支撑的评估的系 统。

参照图4，示出了这个系统的实施方式。系统400包括数据接收单元 402，用于接收管状结构典型的图像数据404.

系统400也包括放置单元406，可操作地耦接至数据接收单元402， 用于将每个种子放置在每个相应区域。放置单元406从数据接收单元402 接收图像数据404，并提供所选择的区域408.

在一个实施方式中，如先前详细描述的一样，放置单元406包括一个 模块（没有示出），用于为了放置两个初始种子而识别邻近于管状结构的 端部的初始区域。

在另一个实施方式中，系统400可以包括可选用户界面410，可操作 地连接至放置单元406，用于提供种子放置参数412，从而为操作者提供 用于辅助放置初始种子的装置。在优选实施方式中，用户界面410被省略， 并且种子被自动放置，而没有操作者的任何介入。

系统400包括处理单元，可操作地耦接至放置单元406，用于执行每 个区域生长。处理单元414从放置单元406接收选择区域408，并提供节 点数据416。相关技术人员将理解，处理单元414也将下一个区域数据418 提供给放置单元406，用于将随后的种子放置在相应区域，直到已经考虑 了每个感兴趣的区域。

系统400包括树建立单元420，可操作地耦接至处理单元414，用于 建立树。树建立单元420从处理单元414接收节点数据416，并响应于此 提供树数据422。

系统400包括过滤单元424，可操作地耦接至树建立单元420，用于 根据预定拓扑参数过滤树，从而确定管状结构的拓扑支撑的评估。过滤单 元424从树建立单元420接收树数据422，并响应于此提供评估数据426.

在一个实施方式中，可选用户界面410被可操作地连接至过滤单元 424，用于提供过滤数据428，从而为操作者提供用于辅助过滤树的装置。 在优选实施方式中，用户界面410被省略，并且树被自动地过滤，而没有 操作者的任何干扰。

再次参照图4，在一个实施方式中，系统可以还包括可选显示单元 430，可操作地耦接至过滤单元424，用于接收评估数据426，并显示管状 结构的拓扑支撑的评估。在进一步的实施方式中，用户界面410可以可操 作地连接至显示单元430，用于提供显示参数432.

在没有被图解说明的另一个进一步的实施方式中，系统400可以进一 步包括传输单元，耦接至过滤单元424，用于将评估数据426传输至用于 异常性检测的计算机辅助检测单元（没有示出），或用于异常性诊断的计 算机辅助诊断单元（没有示出）。传输单元可以包括无线模块（没有示出）， 用于提供评估数据426的无线传输。相关技术人员将理解，在不背离本申 请的范围的前提下，可以使用根据各种协议的无线模块来传输评估数据。 相关技术人员也将理解，可以使用有线传输。在一个实施方式中，使用互 联网执行传输。

如先前所述的用于确定管状结构的拓扑支撑的评估的系统由于其可 以远程处理数据而具有很大的优势。当然，可以在配备了成像设备的诊所 或医院的前提下获取图像数据，通过公共或私人数据网络将图像数据发送 至远程处理中心，并在处理中心处被处理。随后，所评估的数据可以被发 送至诊所或医院，用于给定医生的可视化分析。

可选地，相关技术人员将理解，所述系统可以集成至成像设备，或被 可操作地连接至其上。

相关技术人员也将理解，在一个实施方式中，用于确定管状结构的拓 扑支撑的评估的方法可被嵌入在处理设备上运行的计算机程序中。计算机 程序可以包括被记录在机器可读介质上的指令，用于执行确定拓扑支撑的 评估的上述方法。

根据另一方面，还提供了根据先前所述方法确定管状结构的拓扑支撑 的评估的交易方法。

在一个实施方式中，通过收费确定管状结构的拓扑支撑的评估。

在另一个实施方式中，交易方法包括：接收图像数据；执行先前所述 的用于确定管状结构的拓扑支撑的评估的方法；并通过收费提供管状结构 的拓扑支撑的评估。

在另一个实施方式中，确定管状结构的拓扑支撑的评估的交易方法包 括：由供应商将如先前所述用于确定管状结构的拓扑支撑的评估的系统提 供给第三方；对系统操作，其中，由第三方通过收费进行操作；并将收取 的费用中的至少一部分返还给供应商。

将可以理解，可以由系统自身操作此处所描述的用于确定管状结构的 拓扑支撑的评估的系统。可替换地，可以由第三方收费操作所述系统。在 一个实施方式中，收费可以为收入的份额，同时在替换实施方式中，收费 可以包括固定费用。

现在，参照图6，示出了处理设备600的实施方式，其中，可以优势 地使用用于确定管状结构的拓扑支撑的评估的方法。

处理设备600包括中央处理单元602、I/O设备604、网络接口电路 608、数据总线606、及存储器610。中央处理单元602、I/O设备604、网 络接口电路608、及存储器610通过使用数据总线606可操作地耦接。

具体而言，中央处理器602用于处理数据指令。网络接口电路608 用于通过数据网络（没有示出）将处理设备600可操作地连接至另一个处 理设备（没有示出）。相关技术人员将理解，可以提供网络接口电路608 的各种实施方式。此外，相关技术人员也将理解，可以根据诸如例如TCP/IP 的各种通信协议操作网络接口电路608.

I/O设备604被用于使用户能与处理设备600互动。相关技术人员将 理解，可以使用I/O设备604的各种实施方式。例如，I/O设备604可以 包括键盘、屏幕、及鼠标的至少其中一个。

相关技术人员将理解，可以提供数据总线606的各种实施方式。

也将理解，可以提供存储器610的各种实施方式。此外，将理解，在 一个实施方式中，存储器610可被用于存储操作系统612、用于确定管状 结构614的拓扑支撑的评估的模块、以及用于操作用于确定管状结构614 的拓扑支撑的评估的模块的数据库616。

相关技术人员将理解，操作系统612被用于管理中央处理单元602、 I/O设备604、网络接口电路608、数据总线606、及存储器610之间的相 互作用。

尽管上面的描述涉及发明者当前预期的特定优选实施方式，但是应理 解的是，在广义方面，本发明包括本文所描述的元素的机械和功能等价物。 例如，所述方法可以应用于不同人的解剖结构及动物结构的检查。