专利名称：用于分割数据集中的结构的方法、系统和计算机程序的制作方法
技术领域：
本发明涉及用于分割数据集中的结构的方法，所述数据集包括为时间序列中的各个时间相位(time phase)获取的多个多维图像，每个多维图 像包括有关目标组织和有关另外组织的数据。本发明还涉及用于分割数据集中的结构的系统，所述数据集包括为时 间序列中的各个时间相位获取的多个多维图像，每个多维图像包括有关目 标组织和有关另外组织的数据。本发明此外还涉及用于分割数据集中的结构的计算机程序，所述数据 集包括为时间序列中的各个时间相位获取的多个多维图像，每个多维图像 包括有关目标组织和有关另外组织的数据。
背景技术：
在开头段落中阐述的方法的实施例从EP 0 542 414 Al已知。该已知 方法被安排来特别地确定心脏中的心室容积。为此，病人的心脏区域的第 一图像和第二图像以小的时间移位的间隔被拍摄。第 一和第二图像相互从 对方减去以产生差值图像，该差值图像指示组织，特别是目标组织比如心 室组织的位置，该组织在所述间隔期间已经移动。差值图像因此指示心室 边界。这个差值图像加上用户控制的门限被用于感兴趣区(region of interest)计算，接着是将图像体素(voxel )自动分类成心室或者非心室。 对于体素大小来计数和调整连续分割的体素，以创建各个心室容积的指示。已知方法的缺点在于它提供了一种本质上基于低级别像素的分割技 术，其质量大大依赖于原始图像质量。而且，已知方法中的门限需要根据 图像对比度进行调整。发明内容本发明的一个目的是提供一种图像分割方法，该方法甚至对于低质量 的图像也是鲁棒的。为此，根据本发明的方法包括下列步骤 -访问时间上有序的图像的二元分割结果，产生时间上有序的二元编码图像，由此对应于该目标组织和另外组织的图像部分;f皮二元编码；-对于来自时间序列的所有时间相位计算成对的连续二元编码图像之间的连续差，产生包括多个离散值的多维时间特征图；-得到与能选自所述多个离散值的预定离散值相对应的数据集中的空间4立置；-访问所述结构的预定义的可变形形状模型；-通过根据所述空间位置变形所述形状模型来分割该结构。本发明的技术措施是基于以下认识-在目标组织和另外组织之间，特别是在血库和心脏成像的图像的其 余部分之间具有灰度级图像对比。-对于时间序列的不同时间相位，存在目标组织的移动；-当采用所想结构的先验构造的形状模型时，图像分割的质量显著提高。应当看到，这些假设被证明不仅对于心脏成像是有效的，而且对于一 般的胸部成像也是有效的，由此运动是由于呼吸循环引起的。还应当看 到，虽然EP 0 542 414 Al的已知方法公开了在时间上延时的两个图像之 间相减，但是该相减是在原始的灰度级编码图像之间执行并且是为了粗略 分割一组(band)心脏运动的目的。代替这样，根据本发明的方法，借助 于执行二元分割，原始图像首先受到粗略分割。二元分割可以手动进行或 者自动实现。此外，执行根据本发明的方法中的相减步骤是用于产生边缘特征的目 的，以便变形适当的可变形形状模型。可变形形状模型的例子在M.Kaus 等人(在2004年8月)在医学图像分析(Medical Image Analysis )上 发表的"在心脏MRI中自动分割左心室(Automated segmentation of the left ventricle in cardiac MRI )，，的245—254页中纟会出。根据本发明 的方法的 一 个特别好处在于，对于基于形状模型和外观模型的分割算法， 根据本发明的方法消除了对构建外观模型的需要，这样筒化了计算方法。 按定义的可变形形状模型不包括任何灰度级信息并且因此可跨成像模态 (modality)转移。由于对于给定的所寻求结构而言可变形模型必须被只 构造一次的事实，本发明的方法是高度通用的和在计算上是鲁棒的。根据本发明的方法，依据多维时间特征图被计算的事件，先验构造的
可变形形状模型根据特征图的非零值的空间分布被变形。这个技术措施是 基于这样的认识，即特征图的非零体素对应心室边界，而零体素对应另外 的组织(或者心室的非边界区域)，该另外的组织对于不同的时间相位基 本上是静态的。而且，特征图的非零体素的符号也提供信息，如果对应于 时间相位t的图像从时间相位t + l的图像中减去，则心内膜边界在心脏收 缩期间将对应于-1体素而在心脏舒张期间将对应于+ 1体素。使用这个信 息，使得能够实现鲁棒且准确的图像分割方法。在根据本发明的方法的实施例中，该方法还包括确定数据集中的感兴 趣区的准备步骤，由此所述感兴趣区被设想为限制该数据集的 一个数据 域。已发现在对应数据集中的感兴趣区的数据域内实现本发明的方法是 优选的，这样丟弃了超出感兴趣区之外的信息。自动确定感兴趣区可以根据任何适合的本身已知的技术来执行，例如从C. Cocosco等人2004年发 表在Proceedings CARS上的 "Automatic Cardiac Region-of Interest Computation in Cine 3D Structural隨I (在摄影3D结构MRI中自动计 算心脏的感兴趣区)"已知的。在根据本发明的又一个实施例中，分割结构被显示在显示装置上。已发现通过在适当的显示装置上显示所得到的分割来提供反馈给操 作者是特别有利的。优选地，分割的结构覆盖在原始图像数据上，该分割 结构可以以包括时间信息的二、三或四维来呈现。还优选地，所呈现的分 割结构特别地使用适当的颜色码来加以突出。在根据本发明的方法的又一个实施例中，该方法还包括#^亍时间上有 序的图像的自动二元分割的又一个准备步骤。优选地，自动分割通过以下方式来使能，即首先借助例如计算包括 空间和时间维度的整个多维数据域的强度直方图，来区分目标组织(特别 是血液)和另外的组织，然后应用适当的二元取门限方法。还应当看到， 为了自动二元分割的目的，可以使用任何本身已知的技术。适当的自动二 元分割技术的例子在1979年1月由N. Otsu发表在IEEE Transaction ons Systems, Man, And Cybernetics, Vol. SMC-9, No. 1 的 "A threshold selection method from gray-level histograms (根据灰度直方图的门 限选择方法)"上给出。根据本发明的系统包括-输入，用于-访问时间上有序的图像的二元分割结果，产生时间上有序的二元编码图像；-访问所述结构的预定义的可变形形状模型； -计算装置，用于-对于来自时间序列的所有时间相位计算成对的连续二元编码图 像之间的连续差，产生包括多个离散值的多维时间特征图；-得到与能选自所述多个离散值的预定离散值相对应的数据集中 的空间位置；-通过根据所述空间位置变形所述形状模型来分割该结构。 优选地，根据本发明的系统还包括显示装置，用于显示所分割的结构。此外，根据本发明的系统还包括数据获取单元，用于获取多维数据集。适当的数据获取单元的例子包括磁共振成像单元、层析成像单元、超声单元、X射线单元及其任意組合。根据本发明的计算机程序包括用于使处理器执行如下步骤的指令 -访问时间上有序的图像的二元分割结果，产生时间上有序的二元编码图像；-对于来自时间序列的所有时间相位，计算成对的连续二元编码图像 之间的连续差，产生包括多个离散值的多维时间特征图；-得到与能选自所述多个离散值的预定离散值相对应的数据集中的空 间位置；-访问所述结构的预定义的可变形形状模型； -通过根据所述空间位置变形所述形状模型来分割该结构。 优选地，根据本发明的计算机程序还包括用以使处理器执行确定数据 集中的感兴趣区的步骤的指令，该感兴趣区被构想为限制数据集的一个数 据域。又优选地，根据本发明的计算机程序又包括指令，用于使处理器执 行时间上有序的图像的自动二元分割的另外准备步骤。
本发明的这些和其它方面将参照附图更具体地加以解释。 图1示出了根据本发明的方法的实施例的示意图。 图2示出了参照图像的所选方法步骤的示意方式。
图3示出了根据本发明的系统的实施例的示意图。图4示出了根据本发明的系统的又一实施例的示意图。图5示出了根据本发明的计算机程序的实施例的流程图的示意图。
具体实施方式
图1示出了根据本发明的方法的实施例的示意图。根据本发明的方法 将使用心脏成像的例子来解释说明。方法1可以示意性划分成三个主要阶 段。阶段2包括准备步骤，即获取适当数据集的步骤3，该数据集特别地 为包括时间上有序的心脏图像的数据集，所述图像是例如借助平衡式快速 梯度回波(bFFE， balanced Fast Field Echo)技术4吏用磁共振成^(象单元 获取的。所获取的数据集(未示出)然后在步骤4中受到适当的二元分割， 其结果在步骤5中^^访问。该结果包括时间上有序的二元编码图像，由此 对应于血液的图像部分^C标注为1，其余部分一皮设置为0。根据本发明的方法的随后阶段12被致力于执行图像处理以便根据图 像数据来分割结构，该结构特别地为心室。因此，在步骤S中执行计算， 由此对应于来自时间序列的时间相位t的在前二元编码图像8a从对应于 时间相位t+l的随后二元编码图像8b中减去。应当看到，任何随后的图 像都适合于此，而不必是直接紧跟的图像t+l。这个操作对于来自所获取 的时间序列的每个时间相位t，重复进行。结果包括一个多维时间特征图 8c,由此非零体素对应于心室边界。在步骤9中，得到对应于特定体素值 (+1或-1 )的空间位置并且将其用于分割心室。为此，在步骤ll中访问 预定义的可变形形状模型，该模型然后在步骤14中根据在步骤9中得到 的空间位置被变形。优选地，在步骤l6中存储分割结果。最后，在根据 本发明的方法的另外阶段22中，分割结果在步骤18中使用适当的显示装 置来显示。优选地，所分割的表面使用二、三或四维可视化技术被覆盖在 原始数据上。还优选地，为了易于理解本结果，所分割的表面^^以适当的 透明模式呈现为颜色码。图2给出了参照图像的所选方法步骤的示意方式。方法"可以开始 于访问原始多维图4象21a、 21b，其中，优选地，描绘了各个感兴趣区R1、 R2。在步骤23a、 23b中，执行时间上有序的图像的自动二元分割(这里 为了简单起见，仅示出了两个图像)。该结果包括时间上有序的二元编码 图像，由此对应于目标组织和另外组织的图像部分被二元编码。接下来，访问所述二元分割的结果并且计算在前二元编码图像和在后二元编码图 像之间的差，产生包括多个离散值的多维时间特征图。应当看到，对于计 算所述差来说有可能选择来自所述时间序列的直接紧跟的图像。可选地， 可以选择来自所述时间上相间隔的序列的任何随后图像去计算所述差。在计算该差之后所得到的图像由25给出。图像25包括具有离散值即-1、 0、 +1的图像元素。根据本发明的方法使用了这样的认识，即非零图像元素的 符号是含有信息的。即，在心脏收缩阶段期间心内膜边界将是-1，而在心 脏舒张期间它将是+1。根据本发明的方法，计算数据集中与能选自所述多 个离散值(-l，0,+l)的预定离散值(-l，+l)相对应的空间位置。在根据 本发明的方法的步骤26中，访问该结构的预定义的可变形形状模型。应 当看到，预定义的可变形模型可以包括统计形状模型，或者它可以包括简 单的几何/拓朴可变形模型。在可变形形状模型被访问之后，所寻求的结 构在步骤27中通过根据产生所寻求结构S的所述空间坐标变形所述形状 模型来分割。图3给出了根据本发明的系统的实施例的示意图。系统30包括输入 32，用于访问时间上有序的图像的二元分割的结果，这样提供了时间上有 序的二元编码图像，由此对应于目标组织特别是血液的和对应于另外的组 织特别是其余部分的图像部分被编码。应当看到，所述二元分割可以已经 手动或自动得到。优选地，对应于目标组织的像素或体素被设置为1，由 此其它像素或者体素被设置为0。数据集可以从适当的存储单元(未示出) 来访问，该存储单元可以是位于本地的或者远程的。可选地和/或附加地， 输入32可以被安排来接收来自适当数据获取单元(未示出)的数据，该 数据获取单元以适当的计算机可读格式提供数据集。附加地，输入32被 安排来访问所寻求结构的预定义的形状模型。优选地，使用了可变形形状 模型。二元分割步骤的结果然后由输入32变为可用于系统30的计算单元 35，该单元被安排来执行根据本发明的图像分割，产生所寻求的解剖结 构，特别是心脏心室。这些步骤是使用本身已知的、上面所解释的相应计 算算法来实现的。系统30的核心由处理器34形成，该处理器34被安排来操作设备30 的部件，该部件为输入32、计算单元35、工作存储器36、以及后台存储 单元38。适当的处理器34的例子是常规的微处理器或者信号处理器，后 台存储单元38 (通常基于硬盘)和工作存储器36 (通常基于RAM)。后台
存储单元38可以用于当适当的数据集(或其一部分)未被处理时对其进 行存储，并且用于存储二元图像分割步骤、计算各对二元编码图像之差的 步骤等等的结果。工作存储器36通常保持已处理的(部分)数据集和该 结构的分割结果。计算单元35优选地包括适当数量的可执行子例程35a、 35b、 35c、 35d、 35e。子例程35a ^^安排来执行图像特別是心脏图像的二 元图像分割，以区分目标物体特别是血液和其他物体特别是心肌，产生二 元编码图像。子例程35b被安排来计算成对连续二元编码图像之间的连续 差，产生包括多个离散值的多维时间特征图。子例程35c被安排来得到数 据集中与选自所述多个离散值的预定离散值相对应的空间位置。优选地， 预定一个特定值携带关于所述结构的特定信息。例如，如果当相位t的二 元编码图像被从相位t+l的二元编码图像中减去时，预定具有值+1的像素 或体素在心脏舒张期间对应心内膜边界，而具有值-1的像素或体素在心脏 收缩期间对应心内膜边界。子例程35d被安排来访问可变形形状模型，或 者可选地，它可被安排来基于适当的多个形状参数而建立可变形形状模 型。子例程35e被安排来通过根据在步骤35c中确定的空间坐标来变形形 状才莫型35d，而分割结构。根据本发明的系统30还包括覆盖编码器37，被安排来产生用分割步 骤结果适当覆盖原始数据的表现(rendering)。优选地，所计算的覆盖 被存储在文件37a中。优选地，覆盖编码器37、计算单元"和处理器" 可由计算机程序33操作，该程序优选地存储在存储器38中。输出"用 于输出处理结果，如代表覆盖有分割结构的适当表现的、心脏解剖的覆盖 图像数据。图4给出了根据本发明的系统的又一实施例的示意图。系统40一皮安 排来在多维数据集中分割结构，特别是心室，该数据集包括多个时间上间 隔开的心室图像。优选地，系统40包括数据获取单元41，特别是用于获 取数据集的磁共振成像仪、层析成像单元、超声设备、或者X射线单元。 通常该数据被设想为借助适当编码的信号S从数据获取单元41转移到处 理器42。该处理器执行适当的数据分割，如参照图2解释的，由此在它的 输出上可以产生多种可能的数据。例如，有可能数据"a包括左心室的分 割，数据42b提供右心室的分割，且数据42c提供心外膜的分割，数据42c 由数据42a和42b得到。数据42a、 42b、 42c或者其适当组合之一变为可用于适当观察器43 选地，另外的输入45包括适当的另外的处理器，该处 理器被安排来使用适合于控制用户界面48的程序46操作一个适当的接 口，这样使得解剖数据的图像适当地用分割步骤结果，特别是用数据42a、 42b和/或42c覆盖，这样产生图像部分48a、 48b、 48c。优选地，为了用 户方便起见，观察器43配备有高分辨率显示装置47,用户界面借助适当 的交互装置49例如鼠标、键盘或者任何其它适当的用户输入设备可操作。 图5给出了根据本发明的计算机程序的实施例的流程图示意图。计算 机程序50可以示意地被划分成三个主要阶段。阶段52包括准备步骤，即 获取适当数据集的步骤53，所述数据集特别地为包括时间上有序的心脏图 像的数据集。在步骤54中，所获取的数据集(未示出)然后受到适当的 手动或者自动的二元分割，其结果在步骤55中被访问。该结果包括时间 上有序的二元编码图像，由此对应血液的图像部分优选地-故标注为1，其 余被设置为0。根据本发明的计算机程序的操作的随后阶段62致力于执行图像处理 以便根据图像数据分割结构，特别地是心室。因此，在步骤58中执行计 算，由此对应于来自时间序列的相位t的在前二元编码图像5&被从对应 于相位t+l的随后的二元编码图像58b中减去。这个操作为来自所获取的 时间序列的每个相位t'重复进行。结果包括一个多维时间特征图58c,由 此非零体素对应于心室边界。在步骤59中，对应于特定体素值(+l或-1) 的空间位置被得到并且将被用于分割心室。为此，在步骤61中访问预定 义的可变形形状模型，该模型然后在步骤64中根据空间坐标被变形，所 述空间坐标在步骤59中得到。优选地，在步骤66中存储该分割结果。最 后，在根据本发明的方法的另一阶段72中，分割结果在步骤68中使用适 当的显示装置被显示。优选地，所分割的表面使用二或三维可视化技术被 覆盖在原始数据上。还优选地，为了易于理解本结果，所分割的表面被以 适当的透明才莫式呈现为颜色码。
权利要求
1.一种用于分割数据集中的结构的方法，所述数据集包括为时间序列中的各个时间相位获取的多个多维图像，每个多维图像包括有关目标组织和有关另外组织的数据，所述方法包括步骤-访问时间上有序的图像的二元分割结果，产生时间上有序的二元编码图像，由此对应于该目标组织和另外组织的图像部分被二元编码；-对于来自时间序列的所有时间相位计算成对的连续二元编码图像之间的连续差，产生包括多个离散值的多维时间特征图；-得到与能选自所述多个离散值的预定离散值相对应的数据集中的空间位置；-访问所述结构的预定义的可变形形状模型；-通过根据所述空间位置变形所述形状模型来分割该结构。
2. 根据权利要求1的方法，由此所述方法还包括确定数据集中的感 兴趣区的准备步骤，该感兴趣区被构想为限制数据集的一个数据域。
3. 根据以上任一权利要求的方法，由此所述方法还包括在显示装置 上显示所分割的结构的步骤。
4. 根据以上任一权利要求的方法，由此所述方法还包括执行时间上 有序的图像的自动二元分割的另外的准备步骤。
5. —种用于分割数据集(S)中的结构的系统(30)，所述数据集包 括为时间序列中的各个时间相位获取的多个多维图像，每个多维图像包括 有关目标组织和有关另外组织的数据，所述系统包括-输入(32)，用于-访问时间上有序的图像的二元分割结果，产生时间上有序的二 元编码图像；-访问所述结构的预定义的可变形形状模型(35e); -计算装置(35),用于-对于来自时间序列的所有时间相位计算成对的连续二元编码图 像之间的连续差(35a)，产生包括多个离散值的多维时间特征图 (35b);-得到与能选自所述多个离散值的预定离散值相对应的数据集中 的空间位置(35c);-通过根据所述空间位置(35c)变形所述形状模型(35e)来分 割该结构(35d)。
6. 根据权利要求5的系统，由此该系统还包括显示(47)单元用于 显示所分割的结构。
7. 根据权利要求5或权利要求6的系统，由此该系统还包括数据获 取单元(41)用于获取数据集。
8. —种用于分割数据集中的结构的计算机程序，所述数据集包括为 时间序列中的各个时间相位获取的多个多维图像，每个多维图像包括有关 目标组织和有关另外组织的数据，所述计算机程序包括用于使处理器执行 如下步骤的指令-访问时间上有序的图像的二元分割结果，产生时间上有序的二元编 码图像；-对于来自时间序列的所有时间相位，计算成对的连续二元编码图像 之间的连续差，产生包括多个离散值的多维时间特征图；-得到与能选自所述多个离散值的预定离散值相对应的数据集中的空 间位置；-访问所述结构的预定义的可变形形状模型；-通过根据所述空间位置变形所述形状模型来分割该结构。
9. 根据权利要求8的计算机程序，由此该计算机程序还包括指令， 用于使处理器执行确定数据集中的感兴趣区的步骤，该感兴趣区被构想为 限制数据集的一个数据域。
10. 根据权利要求8或权利要求9的计算机程序，由此该计算机程序 还包括指令，用于使处理器执行时间上有序的图像的自动二元分割的步 骤。
全文摘要
根据本发明的方法(1)可以示意性划分成三个主要阶段。阶段(2)包括准备步骤，即获取适当数据集的步骤(3)，该数据集然后在步骤(4)中受到适当的二元分割，其结果在步骤(5)中被访问。该结果包括时间上有序的二元编码图像，由此对应血液的图像部分被标注为1，其余被设置为0。根据本发明的方法的随后阶段(12)致力于执行图像处理以便分割结构。在步骤(8)中执行计算，由此对应于来自时间序列的相位的在前二元编码图像(8a)从对应于一个时间相位的随后二元编码图像(8b)中减去，产生多维时间特征图(8c)。在步骤(9)中，对应于特定体素值的空间位置被得到并且被用于分割该结构。在步骤(11)中访问预定义的可变形形状模型，该模型然后在步骤(14)中根据在步骤(9)中得到的空间坐标被变形。优选地，在步骤(16)中存储分割结果。最后，在根据本发明的方法的另外阶段(22)中，该分割结果在步骤(18)中使用适当的显示装置来显示。优选地，所分割的表面使用二、三或四维可视化技术被覆盖在原始数据上。还优选地，为了易于理解本结果，所分割的表面被以适当的透明模式呈现为颜色码。本发明还涉及用于分割数据集中的结构的设备和计算机程序。
文档编号G06T5/00GK101151638SQ200680010359
公开日2008年3月26日 申请日期2006年3月21日 优先权日2005年4月1日
发明者C·A·科科斯科, T·尼特施 申请人:皇家飞利浦电子股份有限公司