CN104217418A

CN104217418A - 钙化血管的分割

Info

Publication number: CN104217418A
Application number: CN201410224743.7A
Authority: CN
Inventors: M.凯尔姆; F.卢高尔; 张竞丹; 郑冶枫
Original assignee: Siemens AG
Current assignee: Siemens Medical Ag
Priority date: 2013-05-31
Filing date: 2014-05-26
Publication date: 2014-12-17
Anticipated expiration: 2034-05-26
Also published as: CN104217418B; US20140355850A1; EP2808844A2; KR20140141526A; EP2808844A3; US9430827B2

Abstract

本发明涉及一种用于在图像数据中分割钙化的血管的方法。该方法包括步骤：提供该血管的血管树表示，提供该血管的多个横截面的多个初步边界表示，提供在所述多个横截面中的图像数据中的多个强度一概图，基于强度一概图确定横截面中的钙化，将每个初步边界表示校正为校正的边界表示，其排除了来自于血管内部的钙化。本发明还涉及一种用于该目的的分割系统。

Description

钙化血管的分割

技术领域

本发明涉及一种用于分割图像数据中的钙化血管的方法。本发明还涉及一种用于该目的的分割系统。

背景技术

据美国心脏协会称，冠状动脉疾病(CAD)是导致西方世界死亡的主要原因。对于CAD的当前诊断标准是通常的介入血管造影(ICA)，其涉及相当大的风险和成本。新型心脏计算机断层成像(CT)扫描仪使得能够采集具有前所未有的质量的冠状动脉CT血管造影(CCTA)。在CCTA中可以以高的灵敏度探测到冠状动脉狭窄，这允许该方法被用作为诸如ICA的介入诊断和外科过程的守卫。

用于在CCTA中自动探测冠状动脉狭窄的方法已经被提出用于临床试验。近年来，CCTA也已经被提出用于仿真沿着冠状动脉狭窄的压力分布和用于计算所谓的冠脉血流储备分数(FFR)，其表示了局部缺血导致的损害。冠状动脉狭窄的自动探测以及其血液动力学相关度的仿真(即冠状动脉内的压力下降的基于仿真的探测)取决于在提供的图像数据中的冠状动脉腔的精确分割。这是一项具有挑战性的任务，因为冠状动脉比较小(在其远端部分中仅延伸至在图像数据中的一些体素)而CCTA图像体积的质量又是变化的(特别地与噪声、伪影、对比度均匀性等有关)。精确的分割是更复杂的，因为脉管腔的对比度仅稍微地高于未钙化斑块的对比度和低于钙化斑块的对比度。因此，未钙化斑块在视觉上看起来与脉管的背景非常相似，特别是在对比度增强的图像采集过程中。另一方面，钙化斑块在这样的对比度增强的图像采集过程中看起来是脉管腔的部分，因为它们在图像数据中显示了与造影剂大约相同的外观。因此，在钙化斑块和血液能够流通的血管腔之间的区分是很难作出的。

发明内容

由此本发明的目的是提供一种提高的可能性，用于在图像数据中分割钙化的血管。

该目的通过按照本发明的方法和按照本发明的分割系统实现。

按照本发明，上面提到的方法包括步骤：

-提供血管的血管树表示。

这样的血管树表示例如可以包括血管(树)的多个走向线，诸如中线等。

-提供血管的多个横截面的多个初步边界表示。这样的初步边界表示例如可以包括横截面的边界图。

-提供在多个横截面中在图像数据中的多个强度一概图(intensityprofile)。强度一概图表示了横截面中图像数据的强度并且由此取决于用于图像数据采集的采集方法。例如，在CT中，强度一概图基于(即对应于)图像数据中的豪恩斯菲尔德单位值。

-基于强度一概图确定横截面中的钙化。强度一概图现在被用于探测钙化，例如钙化斑块。

-将每个初步边界表示校正为校正的边界表示，其排除了来自于血管内部的钙化。

在该上下文中，可以注意到，图像数据可以是虚拟的任何医学图像数据，特别是断层成像图像数据。它们可以包括二维和/或三维和/或四维图像数据(在此其一个维度可以是时间的维度，其提示了图像数据包括随时间移动的图像)。

根据本发明，上面提到的分割系统包括以下单元：

-用于图像数据的输入接口。经由该接口将图像数据馈入到分割系统，例如从PACS或从成像设备。

-第一供应(provision)单元，实现为用于提供血管的血管树表示。这样的供应单元可以实现为从图像数据中导出血管树表示的计算单元，但是也可以仅仅包括输入接口-其也可以与第一输入接口组合，例如作为一个公用的接口。

-第二供应单元(也是输入接口)，其实现为用于提供血管的多个横截面的多个初步边界表示。该第二供应单元可以与第一供应单元组合或单独存在。同样，其可以实现为从图像数据导出初步边界表示的计算单元或其可以包括输入接口，该输入接口又可以与任意上述其他输入接口以相应方式自由组合。

-第三供应单元，其实现为用于提供在多个横截面中在图像数据中的多个强度一概图。该第三供应单元优选实现为基于图像数据提供强度一概图的计算单元。

-确定单元，其在操作中基于强度一概图确定横截面中的钙化。该确定单元又优选实现为基于强度一概图确定钙化的计算单元。

-校正单元，其在操作中将每个初步边界表示校正为校正的边界表示，其排除了来自于血管的内部的钙化。

此外，本发明还涉及一种具有采集单元和按照本发明的分割系统的成像设备。

按照本发明的分割系统，特别是其第一和/或第二和/或第三供应单元、确定单元、和校正单元(以及以下提到的分割系统的其他组件)可以部分地或全部地通过硬件组件实现，例如使用半导体芯片，诸如ASIC(专用集成电路)、FPGA(现场可编程门控阵列)，或PLA(可编程逻辑阵列)。然而它们也可以由软件组件或硬件和软件组件的组合来组成。因此，本发明也涉及一种计算机程序产品。计算机程序产品可以直接加载到可编程分割系统的处理器中，包括当计算机程序产品在分割系统上被执行时用于执行按照本发明的方法的所有步骤的程序代码资源。

本发明的特别优选的实施方式和特征通过独立权利要求如在以下说明书中解释地那样给出。不同的权利要求类别的特征可以被组合为适合于给出在此没有描述的其他实施方式。

按照一种优选实施方式，所述多个强度一概图在多个横截面中、沿着从在每个横截面中位于血管内部的一个种子点出发向着和可能超出初步边界表示的多个一概图射线被确定。

也优选的是，钙化在其强度高于预定的阈值水平的区域中被确定。由此，优选地，图像数据是基于X射线采集(例如像在CT情况下那样)的并且阈值水平高于576HU，也就是钙化的典型UV水平。

更优选的是，将来自于钙化确定的结果与初步边界表示的延伸进行比较，以验证该结果。这提供分割结果的更高精度并且构建了一种结果的交叉检验。

优选地，如上所述，血管树表示包括走向线，特别是血管的中线。由此更优选地是，中线至少一次，优选在校正初步边界表示之后至少一次被再中心化。这有助于提高在迭代处理中中线的精度。

此外，优选的是，图像数据被变形(warped)，从而血管树表示沿着直线被对齐并且所述多个初步的和/或校正的边界表示互相平行地对齐。这使得处理的其他步骤更容易处理并且降低计算开销。

优选地，初步边界表示的数量通过相应的横截面中的候选点的导出似然性值来产生。这隐含了，初步边界表示并不固定地代表血管的边界而是包括对于按照概率来说血管边界位于何处的似然性指示。因此，似然性值优选通过使用诸如概率推进树或随机森林算法等训练过的分类器算法来得到。

一般地，按照本发明的方法可以是分割步骤的任何较大数学框架的部分。然而优选的是，其被包括在特别优选的较大分割方法中，该方法包括步骤：

a)提供钙化的血管的图像表示，

b)提供初始表面模型，其包括具有通过边、特别是无向边连接的多个顶点的网格，

c)对于多个顶点的每个顶点定义在(相对应的)顶点的位置处与表面模型正交的一条射线。

d)将多于两个(例如五个)标签分配给每个顶点，每个标签表示该顶点在射线上的候选位置。

e)基于在校正的边界表示中，提供对于每个候选位置的似然性的表示，该似然性是指候选位置是否与器官的、即在图像表示中的血管的表面点相对应。

f)定义具有离散的多元随机变量的一阶马尔科夫随机场，所述随机变量包括候选位置的标签和似然性的表示。

g)通过使用在该马尔科夫随机场中的最大后验估计来找到，即，搜索和/或识别血管的最优分割(诸如分割模型和/或网格)。因此，最优的分割的特征在于如下事实，即，分割显示了分割模型表示血管腔的真实表面的最优似然性。

附图说明

本发明的其他目的和特征从以下结合附图的详细描述中将变得明显。但是应当理解，附图仅为了解释的目的而设计并且不是对本发明的限制。它们也没有必要按比例画出。

图1示出了按照本发明的一种实施方式的分割过程的示意性框图，

图2示出了在按照图1的处理步骤对走向线的校正过程中，血管结构的剖面图及其该走向线，

图3示意性示出了图1的处理的变形步骤的不同步骤，

图4示出了可以在图1的处理的上下文中被使用的管状坐标系，

图5示出了图4的坐标系内部的层的示意性表示，

图6示出了在图1的处理的上下文中边界图的结构的示意图，

图7示出了在图1的处理的上下文中特征提取处理的示意表示，

图8示出了能够在图1的处理的上下文中被使用的、具有注解(annotated)的边界图的血管的横截面的两个图像，

图9示出了具有相对应的第一曲线的血管的横截面图像，该第一曲线在图1的处理的上下文中用于探测钙化斑块，

图10示出了具有相对应的第二曲线的血管的相同的横截面图像，该第二曲线在图1的处理的上下文中用于探测钙化斑块，

图11示出了用于在图1的处理的上下文中使用的血管的MRF表示，

图12示出了凸函数的三个结果曲线，所述函数能够在图1的处理的上下文中被替换地使用，

图13示出了用于在图1的处理的上下文中的最大流分析的构建的图。

具体实施方式

在该整个框架中钙排除步骤也可以基于与在该例子中描述的不同的技术。特别地，最后的步骤-也就是最终分割步骤-可以按照完全不同的逻辑来实现。

在此示出的腔分割框架在多个阶段中执行，其以类似于管线的方式产生腔表面的网格表示。图1给出了框架的概览并且显示主要控制和数据流。

分割Z在开始点Y开始并且具有第一步骤X，其中体积数据ID，即，图像数据ID和(例如前面跟踪的)血管树VT作为输入，而分割的冠状动脉-按照网格表示，即，包括了网格分割数据MSD的网格分割-在结束点S处形成该框架的输出。

首先，算法的输入数据ID、VT、CVT在第一-可选的-血管树(特别是中线)校正步骤X(其产生校正的血管树CVT)中和在从中得到变形的体积数据WID的体积变形步骤W中被预处理。然后，潜在的腔壁边界在边界探测步骤V中沿着圆柱形坐标被探测并且存储在边界图BM中。为了从分割中排除钙化的区域，在步骤U中对于钙的存在分析该图和如果需要的话修改其-这产生修改的边界图MBM。

最后，在步骤T中，通过将边界图作为优化问题嵌入到具有凸对势的多元马尔科夫随机场中并且然后通过使用标准最大流/最小割算法来求解其而找到最终分割。由此，术语体积(图像-图像数据)与体积数据可互换地被使用。

步骤X：血管树校正

血管树的校正一般地涉及其走向线、更具体地是其中线的校正，其参考图2示出：

腔分割质量高度地取决于血管1的提取的血管树走向线、在此是中线3(左边)的精度。这样的中线3可以通过使用任何可用的走向线(中线)产生算法来计算，有大量这样的算法可用。

该例子的全面分割算法缺省地包括在最终腔分割结果中的中线点。这归因于对中线点的径向距离处的潜在的边界位置进行采样的光线投射方法(比较以下)。在理想情况下，所提取的血管树将总是沿着在腔的中心处的点行进。可是在实践中冠状动脉的提取的血管树可能缺乏精确性并且错误分支，特别是在具有斑5和严重闭塞的区域。此外，在分叉7处可能的是，中线3遵循捷径而不是实际的腔位置(比较图2左边)。再定位组成血管树的点(比较图2右边)可以确保，校正的中线3'总是与腔壁间隔开地延伸，即使出现严重损伤。中线校正X可以利用由Zheng等:Model-Given Centerline Extraction for SeverelyOccluded Major Coronary Arteries.In:Machine Learning in Medical Imaging.2012.pp.10to18.提供的方法来执行。在该参考文献中的算法也提供对腔强度分布(或腔似然性)的估计。这允许对于图像数据的任何图像体素从其强度映射到处于血管腔内部的似然性。

作为血管树校正步骤X的结果，校正的中线3'(至少粗略地-即近似地)通过血管1的腔中心延伸并且避免碰到钙化斑块3。

步骤W：体积变形

步骤W参考图3来解释。在此，产生体积输入图像数据ID的变形的和再采样的版本。

由于血管1的腔边界需要精确地在图像数据ID的体积的、提取的中线3、3'在那里延伸的部分被确定，特别关注那些区域。这些区域通过中线片段来确定，所述中线片段通过在分叉点处分裂(可选地校正的)中线3、3'而产生。

为了获得均匀的层距离和为了由此避免图像失真，中线3首先被再采样为一定的分辨率(例如0.1mm)的再采样中线3''，例如通过使用双三次样条插值。

然后对于再采样中线3''的每个点，体积数据的正交横截面9(层9)通过利用在那个点上居中的平面来切割它们而被提取并且与再采样中线3'的中线方向垂直地延伸。这意味着，对于再采样中线的每个中线点，与再采样中线3''正交的图像层9利用在体素之间的位置处的双线性插值被插值。

这些层9通过一个在另一个顶上平行堆叠而被存储在体积中，其产生变形的体积数据WID或变形的图像数据WID。每个层的尺寸优选被调整为提取的结构的最大预计的尺寸，即，其应当至少覆盖足够大以显示穿过血管1的、具有血管1的最大直径的横截面11的面积。由于冠状动脉的直径不超过大约10mm的尺寸，所以感兴趣平面在该情况下可以安全地被限制到15mm乘15mm的尺寸。如果分割其他血管，则可以相应地应用更大的尺寸。

作为变形步骤W的结果，(再采样)中线3''是在变形的图像空间中的直线并且穿过变形的体积数据WID的每个层的中心延伸。该变换(变形)步骤W的一个优点是，腔分割现在可以在特别地适合于诸如血管等的管状结构的圆柱形坐标系中被执行。

步骤V：边界探测

边界探测步骤的目的是确定在变形的体积中在每个层处将血管腔内部与其壁分离的(初步)边界。为了实现这一点，首先需要找到在潜在的边界位置处的候选并且评估其合适性。换言之，边界点候选被产生并且被分配一个似然性值。

有利的是，在极坐标、分别是圆柱形坐标系中而不是使用笛卡尔坐标系搜索腔壁，因为腔轮廓的探测被降低到沿着一维射线搜索的数量。这样的圆柱形坐标系在图4中画出。在此，层9的高度在变形的体积中通过坐标z表示，而角度k和径向距离r确定了在极坐标中横截面(即层9)中的点。

以下过程也公知为光线投射方法：将层的中心-是腔中心-考虑为一个极，围绕腔中心的密集的径向采样通过使用小的参数空间而变得可行。图6示出了对于变形的体积的层的等距选择z∈[1,Z]而产生的边界候选点，例如对于每个层或对于每第五个层(取决于期望的精度)。在每个这样的层中产生沿着K个射线的R个点(在此5个点)。每个点由此通过层、角度k和沿着其射线的位置1至5被定义。

每个这样产生的候选边界点然后对于其位于血管的腔边界上的似然性而被评估。对于该目的，可以应用很好建立的各种模型和方法，诸如加权的强度差方法、梯度大小方法或概率确定方法，其基于前面之前训练的分类器算法，诸如随机森林或概率推进树算法。预计它们在接近腔边界的位置处产生更高的似然性分值。可以预计，边界值似然性作为在0和1之间的标量来表达。为了方便，候选边界点的获得的似然性值可以被存储在维度为K x Z x R的体积l(k,z,r)中，即边界图B，如图6所示。

边界图中的每个元素应当是非负的值，其在接近真实的腔壁的位置处产生仅高的值。在此可以利用如下事实，即，在腔内部的组织具有比在壁外部的更高的强度(例如HU)值。这意味着，例如，在沿着射线的连续位置处被计算的大的正的前向差可以被用来表示边界。现在存在大量可能性性来将其合并到算法中。最显著的方法中的两个方法是，或者精确地进行求导以获得腔边界分值，或者隐含地使用它们，更精确地，使得不同的导数特征通过在机器学习方法中的分类器来评估以产生所寻求的概率。

相对于机器学习背景的边界探测通常表达为二进制分类问题。由此所需的是一些基本事实注解，例如来自于在训练的数据库中提供的先前的手动分割：训练的分类器被用于在给定特征样本的情况下预测腔壁在每个位置处出现的概率。为了训练分类器，必须使得一组正确分类的观察对于分类器可用，使得其能够从观察中学习区别特性(特征)。这也称为监督的学习。这意味着，需要提供两组特征数据，一组在真实的腔壁处被评估(正的)并且另一组在与边界远离处被计算(负的)。在训练步骤之后，分类器可以被用于对于任何未知的特征样本预测边界概率并且由此能够被用于在边界似然性评估的在此描述的上下文中的评估。

对于在位置(k,z,r)处的每个候选边界点，分类器预测，其是否(即以何种似然性或概率)是血管1的腔边界的部分。其预测是基于从考虑的边界候选点的局部邻居提取的图像数据ID的低水平图像特征F(k,z,r)的，其更详细地参考图7被描述：

对于每个候选边界点CBP，图像特征采样样式SP基于点的CBP局部方向被定义，所述CBP局部方向通过将在相对应的层9中的血管的中点CP和潜在的(候选的)边界点CBP相连的线L来确定。在每个采样位置，然后计算低水平图像特征，诸如强度和梯度。二进制分类器然后使用这些特征、基于手动分割的训练数据的代表性的组而被训练，以确定候选边界点CBP位于血管的边界上的似然性/概率。可以使用任何二进制分类器，诸如上面提到的概率推进树或随机森林分类器，其在测试中都提供相似的高精度的结果。

在训练中，对于每个方向，在射线和基本事实注解11之间的交叉点15被考虑为正的并且在射线上的余下点13被考虑为负的，如图8(左边)所示。图8(右边)示出了基于分类器算法、使用来自于左边以概率图输出的形式的基本事实的边界探测处理的结果。

步骤U：钙排除

对于可靠的腔分割，正确处理钙化斑块是强制的，特别是当边界探测在最宽的意义上基于图像梯度时。由于如下事实，即，在CT图像中的钙的特征在于高的强度值，因此，通常类似于在腔内部捕捉的强度，其通常错误地被分类为腔组织并且边界在斑块和血管背景而不是腔和斑块之间被探测。然而，在冠状动脉中的钙化区域提示威胁生命的狭窄并且从定义来说并不是血管流过的腔的部分并且因此，必须从分割结果中被排除。为此目的，边界要在(钙化)斑块和血管腔之间而不是在(钙化)斑块和血管背景之间被探测。

为了理解，为什么钙化斑块被错误地包括在分割中，发明人针对异常对强度的径向一概图、腔概率和边界概率值进行了分析。它们的一概图展示了当相对应的射线通过钙化区域而不是健康区域时的一定的样式。

图9和10示出血管的相同横截面视图的顶部。在图9中射线17穿过示出的血管的健康部分，而在图10中射线17'在不同的方向在相同的血管内部穿过钙化斑块，所述钙化斑块可以通过在斑块区域中的图像的亮度来区分。在两个图的底部的图上，对于图像强度JJ的一概图，腔似然性LL(比较步骤X)和边界似然性(比较步骤V)已经对于一层和在给定的角度下的特定的射线、诸如在两个图中的两个不同的射线17、17'被提取。水平轴示出了以mm为单位的距离，而垂直轴指按照在1和0之间的数字表示的似然性L并且指没有给出的数字但是具有与576HU的CT图像强度相对应的一组零阈值水平的强度J。

至于图9，图像强度JJ通常低于该组阈值水平。而腔似然性从左到右显示了下降到零的相当稳定的曲线，边界似然性主要显示了在大约位置1.3mm处的峰。

至于图10，图像强度在大约1.5mm的距离处超过该组阈值水平，在大约2.2mm处达到峰值并且在大约2.7mm处下降到低于该组阈值水平。腔似然性LL显示两个峰并且边界似然性BL甚至显示三个峰，最左边的峰与排除了钙化斑块的真实腔边界相对应。边界似然性BL的最右边的峰指的是错误地包括了钙化斑块的边界。校正的边界似然性l_c与排除了钙化斑块的真实边界似然性相对应。这一点在如下的启发式方法中实现：

首先，钙化斑块通过确定高于一定的预定阈值的强度一概图的范围来识别，在此是576HU的该组阈值。为了提高鲁棒性和防止错误响应，阈值可以基于腔强度分布(比较步骤X)、通过向其加上(或从其减去)图像数据依赖方差(优选两次这样的方差)，根据当前图像数据来提高。该组阈值t_CAL由此可以是选择的常数，或者是使用固定的阈值t_f和平均腔强度μ_l加上两次其方差σ_l的、适应性的阈值，从而

t_CAL＝max(t_f,μ_l+2σ_l) (1)

然后，沿着射线17的范围被认为是相关的钙化区域，该射线与中线最接近并且优选具有例如0.3mm的一定的最小长度。除了在中线点处的开始半径r₀的索引之外，还提取了半径r_m，对于所述半径r_m，强度在钙化区域内部是最大的(即在图11中大约2.2mm处图像强度JJ的峰)。然后，沿着射线17'的校正的边界似然性被获得为

l_{c} (k, z, r) = \{\begin{matrix} 1 (k, z, r) & r < r_{0} \\ l (k, z, r) \cdot d_{CAL}^{(r - r_{0}) / Δr} & r_{0} \leq r < r_{m} \\ 0 & r &GreaterEqual; r_{m} \end{matrix} - - - (2)

其中Δ_r是径向采样分辨率(例如Δ_r＝0.1mm)并且d_CAL是衰减系数，其确定了在r₀和r_m之间的边界似然性一概图的衰减速度。例如，d_CAL可以被置为0.6。

如在图10可见，该策略导致仅具有一个峰的校正的边界似然性一概图，所述峰描述了正确排除了钙化斑块的边界。

可以注意到，在血管的边界似然性的校正之后，可以如在步骤X中提出的那样再校正走向线(中线)，因为建模的血管的腔可能具有明显的改变。

步骤T：分割

在最后的步骤中，从校正的边界图找到血管的最优边界。除了(校正的)边界似然性之外，可以考虑三维光滑约束以避免不稳定的轮廓和以实现光滑的表面。该问题可以被表达为具有离散的多元随机变量的一阶马尔科夫随机场(MRF)。找到最优分割的问题然后可以被看作是在该MRF中的最大后验概率(MAP)估计，对于其，已经由Ishikawa,Hiroshi:Exact Optimization for MarkovRandom Fields with Convex Priors.IEEE Transactions on Pattern Analysis andMachine Intelligence.Vol25,No.10,pp.1333-1336提出了基于最大流算法的有效解。该参考文献被明确考虑为本说明书的教导的明确考虑的部分。

更公式化地(比较图11)，利用一组顶点V和一组无向边定义了图G＝(V,E)。MRF由图G，有限组的标签L＝{l₁,...l_R}和在潜在的标签分配X＝L^V(其是所有标签L的数量)的空间上的概率分布P组成。每个标签由此指沿着在顶点的位置处血管的初始表面模型的射线法线的(候选)位置。在此，法线可以考虑为顶点的位置的相邻网格表面的法线的平均法线。

元素X∈X(其也可以被考虑为配置)是图，该图为每个顶点v分配一个在L中的标签X_v。如从Ishikawa的上述参考文献中公知的，可以有效地将以下形式的多标签能量函数最小化

E (X) = \underset{(u, v) &Element; E}{Σ} γ_{uv} g (X_{u} - X_{v}) + \underset{v &Element; V}{Σ} h (v, X_{v}) - - - (3)

其中g(·)是将顶点u和顶点v的标签差加权的凸函数。γ_uv是以下将解释的预定义的光滑系数。

上面给出的和的第一项能够被称为预先(prior)，因为其是取决于数据的，而这些和的第二项可以被称为数据项，因为其收编观察(即在该情况中的边界似然性)。该形式的能量函数可以被转换为等效的图割(graph-cut)问题，其可以精确地在多项式时间内部通过最小割/最大流算法求解，所述算法的几个变型是可用的。

分割问题被如下表达为MRF：MRF的每个顶点v对应于管状网格(比较图11)的网格顶点，其通过角(或射线号)k和高度z表示。

相关联的随机变量X_k,z是多变量，取范围[1...R](对应于标签组L＝{l₁,...l_R})中的整数值并且代表了如下事件，即，表面穿过在高度z、和角度k和半径r处的边界候选。等式(3)的数据项然后被定义为h(v,X_v)＝h(X_k,z)＝-log(l(r,k,z))，在先前的步骤U计算的边界候选的负的对数似然性(负的以便能够随着增加的似然性减小等式(3)的能量函数)。可以注意到，在该上下文中，log(l(r,k,z))也可以指log(l_c(r,k,z))，即，指在钙化排除步骤之后的校正的边界似然性。虽然对于等式(3)中的预先项(光滑项)可以选择任意的凸函数g(·)，但是三种类型的函数是特别合适的：

g(d)＝β|d| (L1范数) (4)

g (d) = β \{\begin{matrix} 0 & | d | \leq α \\ (| d | - α) & | d | > α \end{matrix}

(ε-不敏感)(5)

g (d) = β \{\begin{matrix} (d / α) & | d | \leq α \\ 2 | d / α | - 1 & | d | > α \end{matrix}

(Huber函数)(6)

其中d＝X_u-X_v是在两个整数标签之间的差，所述整数标签在变量X_u和X_v之间，所述变量在MRF图中是相邻的u,v∈E_G。函数通过阈值参数α和斜率参数β被参数化。对于等式(4)的对应曲线C₁、对于等式(6)的C₂和对于等式(5)的C₃在图12中示出。虽然等式(6)可以被考虑为提供最精确的结果，但是来自于等式(4)的结果可以最快地被计算并且等式(5)关于速度和精度位于这两个函数之间的中间。

给定通过无向图G和定义的潜在的势表示的MRF，图割问题对于最小化等式(3)被表达。为此目的，利用KZR+2个顶点和与正容量c(u,v)>0相关联的有向边u,v∈E_H构造有向图H。图以如下方式被构造，即，在H中每个σ–τ切割(即将源σ和汇点τ–cf分离的切割，比较图13)对应于MRF G的配置(即对应于变量分配)X∈X并且切割的成本(即在切割中的H的边的容量c(u,v)的和)是按照等式(3)的该配置的成本。(比较图13)。用于产生图H的伪代码在算法(1)、(2)和(3)中提供，其中H.addCapacity(u,v,a,b)将容量a和b加到在u和v之间的两个有向边，即c(u,v)＝c(u,v)+a和c(v,u)＝c(v,u)+b。

除了两个不同的顶点σ和τ，在H中的每个顶点对应于变量分配X_k,z＝r并且由此通过V_k,z,r表示。代表了来自于等式(3)的数据项的容量被加到在V_k,z,r和V_k,z,r+1之间的边作为H.addCapacity(V_k,z,r,V_k,z,r+1,-log(l(r.k.z)),∞)并且加到如算法(1)中所示的边(σ,V_k,z,1)和(V_k,z,R,τ)。

代表了来自于等式(3)的MRF对势γ_uvg(X_u-X_v)的容量从来自于MRF图G的每个边来构造(比较算法(2))。在相邻的层的变量X_k,z和X_k,z+1之间使用光滑系数γ_uv＝γ_k。每个对势可以要求将几个容量加到图H(比较算法(3))，其作为来自于选择的对势函数(比较等式(4)至(6))的二阶差而被计算为

cap (d) = \frac{g (d + 1) - 2 g (d) + g (d - 1)}{2} - - - (7)

在此，d＝r₁-r₀是在连接的变量X_k1,z1＝r₁和X_k0,z0＝r₀，即，在H中的顶点V_k1,z1,r1和V_k0,z0,r0之间的标签差。在此，对于标签距离d而言cap(d)＝0，其中对势函数g(d)是线性的并且对于这些顶点在H中不创建边。因此，结果的图对于具有线性部分的对势函数，诸如等式(4)至(6)，是更稀疏的，其最终导致更快计算最小割/最大流解。

最后需要指出的是，可以使用任意的最小割/最大流算法来获得等式(3)的最小解，其是对应于最小割的配置X。

尽管以优选实施方式和其变型的形式公开了本发明，等式应当理解的是，在不脱离本发明的范围的情况下可以对其进行大量附加的修改和变型。

为了清楚性起见，应当理解，贯穿该说明书的“一”或“一个”的使用不排除复数形式，并且“包括”不排除其他步骤或元素。

Claims

1.一种用于在图像数据中分割钙化的血管的方法，包括步骤：

-提供所述血管的血管树表示，

-提供所述血管的多个横截面的多个初步边界表示，

-提供在所述多个横截面中在所述图像数据中的多个强度一概图，

-基于所述强度一概图确定所述横截面中的钙化，

-将每个初步边界表示校正为校正的边界表示，其排除了来自于所述血管的内部的钙化。

2.根据权利要求1所述的方法，其中，所述多个强度一概图在所述多个横截面中、沿着从在每个横截面中位于所述血管内部的一个种子点出发向着初步边界表示的多个一概图射线被确定。

3.根据权利要求1所述的方法，其中，所述钙化在如下区域中被确定，所述区域的强度高于预定的阈值水平。

4.根据权利要求3所述的方法，其中，所述图像数据是基于X射线采集的并且所述阈值水平高于576HU。

5.根据权利要求1所述的方法，其中，将来自于对所述钙化进行确定的结果与所述初步边界表示的延伸进行比较，以验证该结果。

6.根据权利要求1所述的方法，其中，所述血管树表示包括所述血管的中线。

7.根据权利要求6所述的方法，其中，所述中线至少一次，优选在校正所述初步边界表示之后至少一次被再中心化。

8.根据权利要求1所述的方法，其中，所述图像数据被变形，从而血管树表示沿着直线被对齐并且所述多个初步边界表示和/或校正的边界表示互相平行地对齐。

9.根据权利要求1所述的方法，其中，所述初步边界表示的数量是通过导出相应的横截面中的候选点的似然性值来产生的。

10.根据权利要求9所述的方法，其中，所述似然性值是通过使用训练过的分类器算法来导出的。

11.根据权利要求1所述的方法，包括步骤：

a)提供钙化的血管的图像表示，

b)提供初始表面模型，其包括具有通过边连接的多个顶点的网格，

c)对于每个顶点定义在顶点的位置处与表面模型正交的一条射线，

d)将多于两个标签分配给每个顶点，每个标签表示该顶点在该射线上的候选位置，

e)基于校正的边界表示，提供对于每个候选位置的似然性的表示，该似然性是指候选位置是否与在所述图像表示中的血管的表面点相对应，

f)定义具有离散的多元随机变量的一阶马尔科夫随机场，所述随机变量包括候选位置的标签和似然性的表示，

g)通过使用在该马尔科夫随机场中的最大后验估计来找到钙化的血管的最优分割。

12.一种用于分割图像数据中钙化的血管的分割系统，包括以下单元：

-用于图像数据的输入接口，

-第一供应单元，实现为用于提供所述血管的血管树表示，

-第二供应单元，实现为用于提供所述血管的多个横截面的多个初步边界表示，

-第三供应单元，实现为用于提供在所述多个横截面中在所述图像数据中的多个强度一概图，

-确定单元，其在操作中基于强度一概图确定横截面中的钙化，

-校正单元，其在操作中将每个初步边界表示校正为校正的边界表示，后者排除了来自于所述血管的内部的钙化。

13.一种具有采集单元和按照权利要求12所述的分割系统的断层成像设备。

14.一种计算机程序产品，能够直接加载到可编程的分割系统的处理器中，包括当计算机程序产品在分割系统上被执行时用于执行按照权利要求1的方法的所有步骤的程序代码资源。