技术领域
[0001] 以下总体涉及具体应用于冠状动脉分割的医学成像。
背景技术
[0002] 使用经分割的冠状动脉血管的数字表示来准确计算虚拟分数血流储备(FFR)取决于经分割的血管的准确性。与使用基于侵入性
导管的实际FFR测量的更常规的方法相比，该虚拟方法是非侵入性的。FFR测量跨冠状动脉狭窄的压
力差，以确定狭窄阻碍
氧气输送至心肌的可能性。
[0003] 冠状动脉分割的当前方法使用心脏区域的三维(3D)体积图像，例如根据计算机
断层扫描(CT)图像或磁共振(MR)图像所生成的。该模态提供足以识别冠状动脉树的较大直径部分的
分辨率，所述较大直径部分包括冠状动脉和从冠状动脉分支并将氧供应到心脏的组织的血管。由于分辨率和伪影(如运动伪影)上的限制，血管并不总是清晰地描绘的，例如模糊，不明显等。这些成像方面与来自对象的冠状动脉树的解剖结构的变化相结合，使得准确表示的冠状血管树很困难。
[0004] 当前的方法使用各种
算法来定义冠状血管树，例如，具有难以控制的参数的
种子生长算法。例如，算法通常使用邻近的或相邻的
体素与当前生长区域的比较来确定是否将所比较的体素添加到当前生长区域。该算法依赖于图像属性(例如强度的变化)以确定是否添加比较的体素，例如，如果差值小于由参数确定的
阈值量，则添加，否则不添加。
[0005] 在传统的冠状动脉树分割方法中遇到的问题之一是分割
泄漏到附近的结构中。例如，在区域生长算法中，这通常发生在边界模糊或较不明显的地方，并且生长算法将附近结构(例如，心室、心肌或
肺动脉)的体素添加到当前生长的区域中。一些算法尝试基于参数的调整来控制这种泄漏。这些参数控制体素组合中的强度阈值或测量的强度，这在计算上是昂贵的，例如长的运行时间，并且可能错过血管树的较小分支，例如，接近图像分辨率的那些分支，例如直径只有几毫米。
[0006] 心脏模型在本领域是已知的，其通常用于心脏分割或心脏运动校正。例如，心脏的数字表示可以包括被表示为三
角形网格的心脏组织。该模型被适配到被成像的心脏，并且通过将不同阶段的模型匹配到图像心脏来确定运动校正。但是，心脏模型通常仅包括(如果有的话)大直径的血管，并且不包括整个血管树。
发明内容
[0007] 本文中描述的各方面解决以上提及的问题和其他问题。
[0008] 以下描述了用于心脏模型引导的冠状动脉分割的方法和系统。心脏模型被适配到图像数据中的心脏的内部表面和外部表面，并且经适配的心脏模型被用于为冠状动脉分割提供空间限制。经适配的冠状动脉概率图可以被用于引导分割。
[0009] 在一个方面中，一种用于在三维(3D)心脏图像中分割患者心脏的冠状动脉血管树的系统包括冠状动脉体积定义单元和冠状动脉分割单元。所述冠状动脉体积定义单元基于经适配的心脏模型，根据3D心脏图像中的心脏组织的内部表面和外部表面来设置空间边界。所述冠状动脉分割单元使用分割算法，利用由根据心脏组织的内部表面和外部表面设置的空间边界限定的搜索空间来分割3D心脏图像中的冠状动脉血管树。
[0010] 在另一方面中，一种在三维(3D)心脏图像中分割患者心脏的冠状动脉血管树的方法包括基于经适配的心脏模型来根据3D心脏图像中的心脏组织的内表面和外表面来设置空间边界。使用分割算法，利用由根据所述心脏组织的所述内部表面所述和外部表面设置的所述空间边界限定的搜索空间，对所述3D心脏图像中的所述冠状动脉血管树进行分割。
[0011] 在另一方面中，一种用于在医学图像中分割患者心脏的冠状动脉血管树的系统包括一个或多个处理器，所述一个或多个处理器被配置为使用分割算法来分割图像中的冠状动脉血管树，所述分割算法具有由使用经适配的心脏模型根据图像中的心脏组织的内部表面和外部表面设置的空间边界限定的搜索空间。
附图说明
[0012] 本发明可以采取各种部件和各部件的布置以及各种步骤和各步骤的安排的形式。附图仅出于图示优选的
实施例的目的并且不应被解释为对本发明的限制。
[0013] 图1示意性地图示了心脏模型引导的冠状动脉分割系统的实施例。
[0014] 图2图示了适配到患者图像中的内部和外部心脏表面的心脏模型的示例性横截面以及空间边界确定。
[0015] 图3图示了示例性冠状动脉概率图。
[0016] 图4以
流程图示出基于经适配的心脏模型来分割冠状动脉血管树的方法的实施例。
具体实施方式
[0017] 首先参考图1，示意性地图示了心脏模型引导的冠状动脉分割系统100。诸如CT扫描器、MR扫描器、其组合等的医学成像设备110生成三维(3D)体积心脏图像120。3D心脏图像120具有足够的空间分辨率以显示冠状动脉，例如毫米级或更好。生成的图像可以包括使用
造影剂，造影剂对冠状动脉管腔进行造影。生成的3D心脏图像120可以前瞻性地或追溯性地被
门控以减少运动伪影。3D心脏图像120可以包括心脏和周围组织或心脏和周围组织的部分。
[0018] 心脏模型适配单元130将本领域已知的解剖学心脏模型140适配到3D心脏图像中的心脏。心脏模型140是心脏的组织的数字空间表示。心脏模型可以包括标记的解剖点或区域。心脏模型适配单元130将心脏模型弹性地适配到被成像的心脏，例如将心脏模型针对特定患者个体化。例如，组织的表面可以被表示为三角形网格。该适配可以基于心脏模型140中的解剖点或区域与3D心脏图像120中的心脏中的解剖点的匹配，例如瓣膜、
顶点、隔膜、
腱索等的
位置。弹性网格的点或区域可以被锚定在匹配的解剖点或区域处，并且可以调整弹性网格以匹
配对应的组织的表面。
[0019] 在患者心脏的表面被配备以心脏模型的情况下，冠状动脉体积定义单元150根据心脏组织的表面来设置空间边界。搜索空间包括心脏组织的表面之间的心肌组织以及设置在该空间区域中的动脉供给。在一些情况下，这确保了深入心脏组织的动脉血流被包括在搜索空间中。心脏外部的第一边界距离d1针对在外部的心脏外表面被设置，其排除诸如肺脉管系统的组织。针对心脏内部的表面设置第二边界距离d2，其排除诸如心室的体积。例如，第一距离可以垂直于外部心脏表面在外部延伸1.5厘米(cm)，并且第二距离垂直内部心脏表面在内部延伸1.0cm。另一个示例使用3.0cm和2.0cm作为外部和内部距离。在另一个示例中，1.5cm或3.0cm的距离在外部和内部是相同的。预期其他距离。器官内部和外部的边界距离可以不同或相同。在一个实施例中，边界距离可根据沿着心脏表面距离而渐变，朝向心室距离减小。在一个实施例中，渐变距离包括从冠状动脉沿心脏表面开始的空间距离的函数。例如，随着动脉直径减小，例如通过分支，搜索空间可以被移动接近心脏表面。
[0020] 冠状动脉口查找器单元160基于经适配的心脏模型从左冠状动脉口和右冠状动脉口寻找初始种子点。例如，冠状动脉口查找器单元160使用来自网状模型的标记的三角形开始搜索升主动脉模型壁，从而识别对应于左冠状动脉口和右冠状动脉口的从中等强度到高强度的对比转变。
[0021] 冠状动脉辅助单元170生成并保持冠状动脉概率图172，其包括每个体素是冠状动脉血管树的部分的概率。部分基于经适配的心脏模型的较大血管，并且基于血管树跨样本群体的分布可以针对冠状动脉和分支生成概率图172，并且将其与心脏模型匹配或融合。例如，来自健康患者群体的经分割的冠状动脉样本被可
变形地配准到心脏模型中，并且在每个体素中将概率计算为经分割的冠状动脉发生的
频率的函数。根据心脏模型的根据每个体素或空间位置的概率形成概率图172。在一个实施例中，可以将概率图存储为心脏模型的一部分，其中，每个表面元素存储针对围绕表面元素的一组体素的概率。由于心脏模型是弹性地适配的，因而概率图也是弹性地适配的。在一些实施例中，经适配概率图172可以被用于排除分割中的环路，例如在血管节段
水平。当血管表现为只给自己供给时发生环路。在一些实施例中，适配概率图172可以与阈值一起使用，例如在分支的分叉中。例如，动脉在0-70度的角度内分支，其可以被用作阈值以限制相邻体素的局部搜索空间。
[0022] 冠状动脉分割单元180利用由冠状动脉体积定义单元150确定的空间边界界定的分割算法来分割冠状动脉血管树182。合适的算法可以包括树分割
框架或距离图方法，例如T.Buelow等人的“A general framework for tree segmentation and reconstruction from medical volume data”，利用包括血管体素、血管节段、以及整个血管树的度量使用快速前进作为区域扩展。另一个范例是T.Klinder等人的“Robust Peripheral Airway Segmentation”中描述的方法，其可适应于动脉血管，并使用基于
滤波器的血管候选来检测和分割连接功能。距离图方法的范例可以包括S.Svensson等人的“Digital Distance Transforms in 3D images using information from neighborhoods up to 5x5x5”以及R.Calvin等人的“A linear time algorithm for computing exact Euclidean distance transforms of binary images in arbitrary dimensions”。也预期其他算法。在某些情况下，空间边界通过减少搜索空间来缩短分割运行时间。在一些情况下，通过消除通常造成泄漏的结构(例如心室和肺脉管系统)来增加分割鲁棒性。用于包含在经分割的血管树182中的接受准则可以包括基于体素的接受，例如距离测量、概率和/或定向流动、血管分割接受(诸如下一分支的可能性)和/或
全血管树水平接受。
[0023] 心脏模型适配单元130、冠状动脉体积定义单元150、冠状动脉口查找器单元160、冠状动脉辅助单元170和冠状动脉分割单元180包括一个或多个配置的处理器190，例如
微处理器、
中央处理器处理单元、数字处理器等。所述一个或多个配置的处理器190被配置为执行存储在计算机可读存储介质中的至少一条计算机可读指令，所述计算机可读存储介质排除瞬态介质并且包括物理
存储器和/或其它非瞬态介质以执行本文中所描述的技术。所述一个或多个处理器190也可以运行由载波、
信号或其他暂态介质承载的一条或多条计算机可读指令。所述一个或多个处理器190可以包括本地存储器和/或分布式存储器。所述一个或多个处理器190可以包括用于通过网络192的有线和/或无线通信的
硬件/
软件。例如，线指示各种部件之间的可以是有线或无线的通信路径。所述一个或多个处理器190可以包括计算设备194，诸如台式、膝上型、身体穿戴式设备、智能电话、
平板电脑和/或包括一个或多个配置的
服务器的协作/分布式计算设备(未示出)。计算设备194可以包括显示设备196，显示设备196可以显示经分割的血管树182或根据经分割的血管树182进行的计算(诸如FFR)。计算设备194可以包括一个或多个输入设备198，所述一个或多个输入设备198接收命令，诸如识别3D心脏图像120，通过冠状动脉口查找器单元160确认种子点，显示冠状动脉概率图172，和/或分割冠状动脉血管树182。
[0024] 3D心脏图像120，心脏模型140，分割的冠状动脉树182被表示为存储在
电子存储介质或计算机存储器上的数字数据集。3D心脏图像120可以包括医学数字影像和通信(DICOM)格式或其他合适的图像格式。
[0025] 参考图2，图示了根据3D心脏图像120和空间边界确定的适配到心脏的内部表面和外部表面的经适配的心脏模型200的示例性横截面。在经适配的心脏模型200中表示的组织(例如心肌组织205)被图示为
叠加在3D心脏图像120中表示的组织上方。根据经适配的心脏模型200的外部表面来确定第一边界210，d1。第一边界被定义为通过从经适配的心脏模型200的网格表面的
正交投影确定的距离d1的表面。将第一边界210外部的组织从分割算法的搜索空间中排除。
[0026] 根据经适配的心脏模型200的内部表面来确定第二边界220，d2。由第二边界220限定的表面包围的组织也被从分割算法的搜索空间中排除，例如心室体积。两个边界确定是用距离进行的，使得搜索空间包括冠状动脉血管树包括外表面和内表面之间的心肌组织以及它们的相关联的动脉供给。在一些实施例中，边界确定排除根据经适配的心脏模型200识别的基于适配概率图具有包括冠状动脉血管树的低概率的组织。
[0027] 根据经适配的模型200，可以
定位冠状动脉的内腔230，并且确定针对分割算法的初始种子点，例如左冠状动脉和右冠状动脉的口。所述口可以被自动识别并由健康护理从业者视觉地确认。
[0028] 参考图3，图示了示例性冠状动脉概率图172。概率图172在视觉上被图示为对3D图的二维(2D)表示，其具有根据相应的概率变暗的2D
像素，例如具有较高概率的较暗像素属于冠状动脉血管树。冠状动脉概率图172的概率被表示为对应于空间定位的3D体素。
[0029] 冠状动脉概率图172可以包括针对冲突的结构的概率，例如冠状静脉系统。冠状动脉概率图172可以包括距离值，例如最近的动脉与心内膜表面之间的典型距离值，例如，利用经适配的
网格模型的三角形表示。冠状动脉概率图172可以包括基于动脉流的有向加权概率，例如沿着管腔的中心线的相邻体素的概率较高，随着角度远离中心线增加概率逐渐降低。
[0030] 参考图4，以流程图示出了基于经适配的心脏模型分割冠状动脉血管树的方法的实施例。在400处，接收3D心脏图像120。3D心脏图像120可以从可以被存储在诸如图片存档及通信系统(PACS)、放射信息系统(RIS)、电子医学记录(EMR)等的存储子系统150中。3D心脏图像120可以由医学成像设备110生成并且直接从医学成像设备110接收。
[0031] 在410处，心脏模型140被适配到3D心脏图像120中的心脏。经适配的心脏模型200包括在3D心脏图像120中的被适配到心脏的组织表面的限定的组织表面。根据经适配的心脏模型200来确定第一边界210和第二边界220。边界的表面限制了搜索空间。基于经适配心脏模型200来适配冠状动脉概率图172并将其包括在搜索空间中。冠状动脉概率图172根据可变形地配准到心脏模型140的经分割的冠状动脉血管树的
采样来构建。概率可以根据采样分布来计算。冠状动脉概率图172可以包括到最近动脉的距离测量值，动脉血流量的方向指示器和/或分支阈值。
[0032] 在420处，冠状动脉血管树182被分割。分割可以包括基于经适配的心脏模型200通过搜索种子点来识别种子点。分割可以包括基于体素、分割或整个血管树处的接受程度来
修改或更新概率图172。分割包括在空间上将搜索空间限制到由第一边界210和第二边界220限定的表面之间的体积。经分割的冠状动脉血管树182包括数字表示，其针对每个空间定位的体素来定义血管腔中是否包括体素。在一个实施例中，数字表示可以包括狭窄，诸如
钙、斑
块等。
[0033] 在430处，冠状动脉血管树182可以被存储在计算机存储器中和/或被显示在显示设备196上。冠状动脉血管树182或其部分可以被显示为2D投影。投影可以包括内部导航、不同的视角、
颜色对比等。这些投影可以包括在显示中对比的测量值，例如直径。投影可以包括狭窄和/或材料的视觉对比，例如斑块，钙等。
[0034] 在440处，可基于经分割的冠状动脉血管树182来计算分数血流储备值(FFR)。FFR可以用于识别和评估血管树182中的狭窄。
[0035] 己经参考优选的实施例已描述了本发明。通过阅读和理解前述的详细描述，本领域技术人员可以想到各种修改和变型。目的是，本发明被理解为包括所有这样的修改和变动，只要它们落入
权利要求书或其等价方案的范围之内。
