CN106133789B - 用于分割感兴趣区域的图像处理装置和方法 - Google Patents

用于分割感兴趣区域的图像处理装置和方法 Download PDF

Info

Publication number
CN106133789B
CN106133789B CN201580015322.6A CN201580015322A CN106133789B CN 106133789 B CN106133789 B CN 106133789B CN 201580015322 A CN201580015322 A CN 201580015322A CN 106133789 B CN106133789 B CN 106133789B
Authority
CN
China
Prior art keywords
image data
deformable model
interest
region
processing apparatus
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201580015322.6A
Other languages
English (en)
Other versions
CN106133789A (zh
Inventor
I·韦希特尔施特勒
J·威斯
C·比格尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of CN106133789A publication Critical patent/CN106133789A/zh
Application granted granted Critical
Publication of CN106133789B publication Critical patent/CN106133789B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4494Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer characterised by the arrangement of the transducer elements
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/483Diagnostic techniques involving the acquisition of a 3D volume of data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/149Segmentation; Edge detection involving deformable models, e.g. active contour models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/100764D tomography; Time-sequential 3D tomography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20116Active contour; Active surface; Snakes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20124Active shape model [ASM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30056Liver; Hepatic

Abstract

公开了一种用于在目标(12)的多维图像数据中分割感兴趣区域(15)的图像处理装置(16)。所述图像处理装置包括:接口,其用于接收包括要被分割的所述感兴趣区域的所述目标的图像数据;选择单元,其选择对应于所述图像数据中的所述目标的解剖结构的可变形模型30;处理单元,其通过基于所述图像数据(xt)和所述目标的额外的解剖信息调整所述可变形模型来分割所述感兴趣区域。

Description

用于分割感兴趣区域的图像处理装置和方法
技术领域
本发明涉及用于在目标的多维图像数据中分割感兴趣区域的图像处理装置。另外,本发明还涉及用于在目标的多维图像数据中分割感兴趣区域的方法。本发明还涉及用于在目标的多维图像数据中分割感兴趣区域的计算机程序。本发明最终还涉及包括图像采集单元的医学成像系统。
背景技术
在医学图像处理的领域中,各种处理任务通常是在如超声图像、MRT图像、计算机断层摄影图像等的医学图像上执行的。一种特定的处理任务(其是在许多图像处理应用中的基本任务)是对例如特定器官的感兴趣区域的分割。所述分割对于识别器官或者对于例如基于体积量化的特别诊断是必要的,以便改善对处置参数的确定。
对于许多器官而言,图像分割能够利用可变形模型而被成功地执行,所述可变形模型是基于具有拓扑的网格结构的,所述拓扑在针对要被分割的图像进行调整期间保持不变。基于模型的分割已经被认为对于各种各样的简单到复杂的器官(例如,骨骼、肝脏、具有嵌套结构的心脏)非常有效。例如从WO 2007/072363获知用于使用可变形网格促进图像的对应的方法。
WO 2011/070464 A2公开了一种用于通过选择体积图像中被成像的感兴趣解剖结构的可变形模型来执行自动分割的系统和方法,其中,所述可变形模型是由包括顶点和边的多个多边形来形成的,其中,检测对应于所述多个多边形中的每个的感兴趣解剖结构的特征点,并且通过朝向对应的特征点移动顶点中的每个直到所述可变形模型变体到所述感兴趣解剖结构的边界来调整所述可变形模型。
WO 2006/029037 A2公开了一种用于定义和跟踪三维图像中的候选解剖结构壁的可变形形状的系统和方法,其中,所述候选解剖结构的形状由多个标记的3D界标点来表示,并且定义在图像帧中的所述可变形形状的至少一个3D界标点。
尽管基于模型的分割总体上提供对器官的可靠且准确的识别和分割,但是如心脏的左心室的心尖或者肝脏的顶端或者脑的脑室的角的一些解剖结构常常难以在医学图像中检测,使得对这些特征的分割是具有挑战性的。
发明内容
本发明的目的是提供一种用于在目标的多维图像数据中分割感兴趣区域的图像处理装置和对应的方法,其具有改进的可靠性和准确性。
根据本发明的一个方面,提供了一种用于在目标的多维图像数据中分割感兴趣区域的图像处理装置,包括:
-接口,其用于接收包括要被分割的所述感兴趣区域的所述目标的图像数据,
-选择单元,其用于选择对应于所述图像数据中的所述目标的解剖结构的可变形模型,
-处理单元,其用于通过基于所述图像数据和所述目标的额外的解剖信息调整所述可变形模型来分割所述感兴趣区域。
根据本发明的另一方面,提供了一种用于在目标的多维图像数据中分割感兴趣区域的方法,包括如下步骤:
-接收包括要被分割的所述感兴趣区域的所述目标的图像数据,
-选择对应于所述图像数据中的所述目标的解剖结构的可变形模型,
-通过将所述可变形模型应用到所述图像数据来分割所述感兴趣区域,并且
-基于所述图像数据和所述目标的额外的解剖信息来调整所述可变形模型。
根据本发明的另一方面,提供了一种用于在目标的多维图像数据中分割感兴趣区域的计算机程序,所述计算机程序包括用于当所述计算机程序在计算机上被执行时令计算机执行如下步骤的程序单元:
-接收包括要被分割的所述感兴趣区域的所述目标的图像数据,
-选择对应于所述图像数据中的所述目标的解剖结构的可变形模型,
-通过将所述可变形模型应用到所述图像数据来分割所述感兴趣区域,并且
-基于所述图像数据和所述目标的额外的解剖信息来调整所述可变形模型。
根据本发明的另一方面,提供了一种医学成像系统,包括:图像采集单元,其用于采集医学图像;以及根据本发明的图像处理装置,其用于在由所述图像采集单元采集的目标的多维图像数据中分割感兴趣区域。
在从属权利要求中定义了本发明的优选实施例。应当理解,所请求保护的方法具有与所请求保护的设备以及如在从属权利要求中所定义的相似和/或相同的优选实施例。
本发明基于这样的思想,即:将解剖结构的可变形模型应用到图像数据并基于图像数据内的目标的已知解剖信息来调节可变形模型。换言之,要被分割的目标的已知解剖特征被用于调节可变形模型以便改进对所选择的可变形模型的分割。借助于目标的已知解剖特征,能够将额外的参数用于调整,使得分割变得更为可靠且更为准确,并且也能够分割复杂的结构以及在图像数据中难以检测的结构。
所述处理单元还被设计为基于所述图像数据和所述目标的解剖特征的已知位置来调整所述可变形模型的元素的位置。这能够将所述可变形模型调节到已知位置,使得降低不确定性并且能够实现更为可靠且更为准确的分割。
在优选实施例中,所述处理单元被设计为基于所述图像数据和所述图像数据的视场内的已知位置来调整所述可变形模型的元素的位置。这能够将已知位置以低的技术努力添加到所述分割中,这是因为所述图像数据的所述视场通常包括边条件,其能够被用作针对所述分割的额外的参数。
在优选实施例中,所述处理单元被设计为基于所述图像数据和所述图像数据的连续时间帧中的对应位置来调整所述可变形模型的元素的位置。这能够利用能从解剖结构的已知运动或已知稳定性导出的目标的解剖特征,使得能够将额外的参数用于增加分割的可靠性和准确性。优选地,对应位置是解剖结构,其通常是随时间稳定的,使得各自的解剖结构能够被相关并且所述可变形模型能够相应地进行调整。
在优选实施例中,所述对应位置是所述目标的已知的对应解剖特征。这能够将不同时间帧的解剖特征相关,使得能够添加时间依赖的特征作为额外的参数来改善分割。
在优选实施例中,所述处理单元被设计为基于所述图像数据和从所述图像数据的视场导出的矢量或平面来调整所述可变形模型的元素的位置。这能够在所述位置在所有三个维度中未知时包括外部信息。
在优选实施例中,所述处理单元被设计为基于所述元素的所述位置与所述已知位置之间的距离或者所述对应位置之间的距离来调整所述元素的所述位置。这能够减少用于调整所述可变形模型的技术努力,这是因为能够在可变形模型中容易地计算距离。
在优选实施例中,所述可变形模型被调整为使包括从所述距离导出的弹性能量(spring energy)的模型能量最小化。
在优选实施例中,所述模型能量包括从所述可变形模型与所述感兴趣区域中的对应图像特征的差异导出的外部能量以及所述弹性能量。这能够将来自不同源的信息进行组合,使得实现对可变形模型形状的总体改善。
在优选实施例中,所述弹性能量通过权重因子进行加权。这能够在总体模型能量内调节弹性能量,使得能够实现在模型能量内的不同能量之间的平衡。
在优选实施例中,所述处理单元被设计为进一步基于所述目标的期望形状来调整所述可变形模型。这能够改善分割,这是因为所述可变形模型针对包括要被分割的感兴趣区域的医学图像进行调整。
在优选实施例中,所述可变形模型是通过包括具有顶点和边的多边形的网格来形成的。这能够实现具有高分辨率的分割,使得能够改善对解剖结构的诊断。
在优选实施例中,所述可变形模型是通过调节所述网格的所述顶点的位置来调整的。这能够以低的技术努力来调整可变形模型。
如上文所提到的,本发明基于这样的思想,即:通过使用关于图像数据的采集或解剖结构或特征的外部知识来改善基于模型的分割的可靠性和准确性,以引导分割。通常,解剖特征是解剖结构、解剖距离或解剖点的位置,其随时间侧向稳定,使得能够添加额外的已知参数,以便降低分割的不确定性并改善模型的可靠性或准确性。外部已知参数被实施为弹性能量,所述弹性能量是从可变形模型的点与图像数据内的对应点的距离导出的,使得能够通过对元素的逐步移动来调整所述可变形模型。
对所述模型调整的总体收敛能够通过使模型能量最小化来实现,所述模型能量包括外部模型能量、内部模型能量和弹性能量。因此,所有的相关调整方法能够被考虑,使得能够实现最优的分割。
附图说明
参考下文描述的(一个或多个)实施例,本发明的这些方面和其他方面将是明显的并且得到阐明。在如下附图中:
图1示出了用于扫描患者的身体的体积的医学成像系统的示意性表示;
图2示出了用于分割感兴趣区域的包括多边形和顶点的可变形模型的示意图;
图3示出了可变形模型的截面视图以便解释对所述模型的调整;
图4示出了基于使用的场中的外部知识来调整可变形模型的实施例;
图5a和图5b示出了基于外部知识的改善的分割的范例;并且
图6示出了用于分割感兴趣区域的方法的流程图。
具体实施方式
图1示出了根据一个实施例的医学成像系统10的示意性图示,具体为医学三维(3D)超声成像系统10。所述医学成像系统10被应用于检查解剖部位(尤其是患者12的解剖部位)的体积。医学成像系统10包括超声探头14,所述超声探头14具有至少一个换能器阵列,所述至少一个换能器阵列具有用于发射和/或接收超声波的多个换能器元件。所述换能器元件优选被布置在二维阵列中,尤其用于提供多维图像数据。
3D超声扫描通常涉及发射超声波,所述超声波照射患者12内的具体体积或目标,其可以被设计为靶体积或感兴趣区域15。这能够通过以多个不同角度发出超声波来实现。然后能够通过接收并处理所反射的波来获得体积数据的集合。所述体积数据的集合是患者12内的感兴趣区域15的表示。
医学成像系统10包括图像处理装置16,所述图像处理装置16用于经由医学成像系统10来提供图像。图像处理装置16控制图像处理,所述图像处理形成通过超声探头14的换能器阵列接收的超声波束的回波的图像输出。图像处理装置16包括控制单元18,所述控制单元18经由超声探头14的换能器阵列来控制图像数据的采集,并且充当用于接收图像数据的接口18。控制单元18被连接到选择单元20,所述选择单元20选择对应于图像数据的感兴趣区域15中的目标的解剖结构的可变形模型,并且将所选择的可变形模型和图像数据提供到处理单元22,所述处理单元22用于通过将从选择单元20接收的可变形模型应用到图像数据来分割感兴趣区域15中的目标。在对感兴趣区域15内的目标进行分割之后,感兴趣区域15的分割图像和超声图像可以被叠加并显示在显示屏24上。另外,输入设备26可以被提供为连接到显示单元24,或者可以被提供为连接到图像处理装置16,以便控制所述超声图像的图像采集、分割和/或显示以及分割结果。
在备选实施例中,医学成像系统10被连接到MR设备、CT设备或X射线设备,以便提供患者的身体12的对应的医学图像,所述医学图像能够通过图像处理装置16而被分割并且被显示在显示单元24上。
图2示出了表示解剖结构的表面以及其空间群集的可变形模型30。可变形模型30是通过选择单元20进行选择的,并且例如从图像处理装置16的存储器加载。从接口18接收的所选择的可变形模型30和图像数据然后被提供到处理单元22以针对图像数据中所识别的解剖结构调整可变形模型。
可变形模型30由包括三角形的网格形成,所述三角形由三个邻近的顶点以及通过边被连接到彼此的模型的表面形成。所述网格形成解剖结构的表面。在备选实施例中,所述解剖结构的表面能够通过替换网格的三角形的多边形来定义。
可变形模型30包括初始空间形状,所述初始空间形状形成平均网格或形状模型,并且如下文所描述的针对图像数据的感兴趣区域内的解剖结构来调整可变形模型30。为了针对图像数据的感兴趣区域15中的解剖结构调整可变形模型30,执行对所述图像数据的模式检测并且相应地将对应位置相关联。
图3示出了可变形模型30的截面视图。可变形模型30包括多个三角形T1到T5。三角形T由顶点v形成,所述顶点v通过边被连接到彼此。每个三角形T包括三角形中心c1到c5
在已经执行对图像数据的模式检测之后,针对所述图像数据内所检测的解剖结构调整可变形模型30,其中,所述三角形T的位置被调节,使得将所谓的总的能量函数最小化:
E=E外部+αE内部
E外部一般是外部能量,其将顶点和三角形T移动到通过模式检测在所述图像数据中所检测到的解剖结构的对应位置,其中,E内部是内部能量,其将顶点和三角形T移动到所述可变形模型的最初的位置,即,移动到平均网格或形状模型。
所述外部能够通过下式来计算:
Figure GDA0001938608730000071
其中,
Figure GDA0001938608730000072
是由所述模式检测所发现的靶点,
Figure GDA0001938608730000073
是从各自的三角形T到靶点
Figure GDA0001938608730000074
的方向,并且ci要被调节到靶点
Figure GDA0001938608730000075
的三角形Ti的中心。在所述三角形T的每个上形成加和,并且每个都通过权重因子wi来进行加权。当通过使外部能量E外部最小化以刚性仿射的方式调整所述模型时,其也被称为对模型30的参数调整。
所述内部能够通过如下公式来计算:
Figure GDA0001938608730000076
其中,v是经调整的可变形模型30的顶点,并且m是如最初提供或选择的形状模型或平均网格的顶点。所述加和在可变形模型30的所有顶点上形成,其中,确定到邻近顶点vj的距离中的每个。当所述模型通过对外部能量和内部能量的最小化进行自由变形时,其也被称为对模型30的可变形的调整。
对总的能量E的最小化引起分割,其是在平均网格与通过模式检测在图像数据中确定的所检测到的靶点之间的介质。对可变形模型30的调整的这一流程也在EP 1251462B1中有详细描述。
根据本发明,可变形模型30也基于解剖结构的已知特征进行调整,以便改善对可变形和参数调整的分割。
所述解剖结构的解剖特征可以是在图像数据内的特定点、位置或距离,其是已知的并且其能够被识别。这可以例如是心脏的心尖,该位置通常是已知的,并且到胸部的皮肤的该距离是公知的并且由于超声探头14的位置而能够从图像数据导出。备选地,所述解剖结构的解剖特征可以是从视场导出的关于超声探头14的位置的平面,或者可以是已知随时间稳定的解剖结构,例如,心脏的心尖,使得在图像数据的不同时间帧中的对应点或位置能够相关联。通常,根据关于所述可变形模型的解剖特征,计算额外的弹性能量,并且将其包括在总的能量E的函数中。
如果在图像数据中的特定解剖结构被识别或确定,则所述额外的弹性能量能够将所述三角形T的特定中心c拉到选定点
Figure GDA0001938608730000081
对已知解剖结构的选择能够借助于图像处理装置16来自动地执行或者可以由用户手动地执行。关于特定选定点形成的弹性能量能够通过如下公式来计算:
Figure GDA0001938608730000082
其中,cj是三角形T的各自的中心,
Figure GDA0001938608730000083
是在图像数据内的选定的对应点,并且wj是权重因子。在图3中,选定点
Figure GDA0001938608730000084
是示意性示出的,其中,弹性能量Es将第一三角形T1的中心c1移动到选定的点
Figure GDA0001938608730000085
因此,弹性能量Es是基于三角形中心c与通过权重因子wj进行加权的对应选定解剖结构
Figure GDA0001938608730000086
之间的差来计算的,以对总的能量的公式内的三种能量进行加权。
在图4中,示出了视场和可变形模型30的示意性图示,以解释针对医学成像系统10的视场内的平面对可变形模型30的元素的调整。在图4中,医学成像系统10的视场以及尤其是超声探头14的视场总体上由32指代。
要被分割的目标的解剖结构是如在图4中所示的患者12的心脏。尤其地,对所述心脏的心尖的分割是困难的,其中,根据本发明,心尖到超声探头14的顶端的距离是大致已知的。超声探头14到心尖的距离是已知的,这是因为超声探头14位于患者12的皮肤处并且心尖到皮肤的解剖距离是大致已知的。在这种情况下,心尖的已知距离能够被用于调整可变形模型30。靶点
Figure GDA0001938608730000091
被实施为在距超声探头14的位置的特定距离处的平面34,并且确定可变形模型30的心尖位置的三角形中心cj到平面34的方向
Figure GDA0001938608730000092
基于这些信息,能够确定弹性能量ES,以调整可变形模型30。关于平面34的弹性能量能够通过如下公式来计算:
Figure GDA0001938608730000093
其中,
Figure GDA0001938608730000094
是平面34的方向,cj是可变形模型30的要被调整的三角形的中心,并且
Figure GDA0001938608730000095
是所述解剖结构实际所位于的靶点。这样确定的弹性能量ES将可变形模型30的心尖cj拉到平面34。在这种情况下,超声探头14的位置被估计为描述在中心切片中的视场的截头椎体的两个边缘射线36、38的交叉点。
在另外的实施例中,使用在图像数据的不同时间帧中的可变形模型30的特定元素的位置。特定解剖结构,尤其是心脏的心尖,是在心脏周期中空间稳定的,使得在图像数据的不同时间帧中的可变形模型30内的心尖的位置能够被用于调整可变形模型30。
在这种情况下,不同的时间帧被一起分割,并且基于可变形模型30的对应元素的距离,尤其是基于图像数据的不同时间帧之间的三角形中心cj的距离,来计算弹性能量。用于调整可变形模型30的时间依赖的弹性能量是通过如下公式来计算的:
Figure GDA0001938608730000096
其中,cj,t1和cj,t2是不同时间帧t1和t2的三角形中心cj,其中,wj是权重因子。如果在不同时间帧t1、t2中的三角形中心cj的位置具有大的差异,则时间依赖的弹性能量Es是高的,并且相应地调整可变形模型30的各自的三角形,以便降低弹性能量Es
应当理解,不仅仅是空间稳定的解剖结构能够被用于调整所述可变形模型以改善分割,在周期内的解剖结构的已知运动也能够被用于调整所述可变形模型以改善分割,使得要被分割的目标的一般解剖行为能够被用作额外的形成。
要在调整期间最小化的可变形模型的总的能量E能够借助于如下公式来计算:
E=E外部+αE内部+ES
所述额外的弹性能量和关于目标的解剖特征的额外的信息能够改善对视场52的分割,尤其是对心脏的心尖的分割。如果心尖在从超声探头14接收的图像数据中不可见,则这是尤其有用的。在图5a中示出了基于额外的解剖信息改善的对心脏的所得到的分割。其已经示出了额外的解剖信息能够改善体积量化,尤其是改善在舒张末期和收缩末期图像中的左心室的体积量化,使得能够改善总体诊断可靠性。
图6示出了图示用于在患者12的多维图像数据中分割感兴趣区域的方法的示意性流程图。所述方法总体上由40指代。
方法40以借助于超声探头14对图像数据的采集以及借助于在步骤42处所示的图像处理装置16的接口18对图像数据的接收开始。在接收到所述图像数据之后,由选择单元20选择对应于如在步骤44中所示的要被分割的感兴趣区域15中的目标的解剖结构的可变形模型30。可变形模型30是从处理装置16的存储器中选择的。在步骤44处的选择之后,通过将可变形模型30应用到在步骤46中所示的图像数据而对感兴趣区域15进行分割。可变形模型30最终基于在步骤48处所示的图像数据中的目标的已知特征而被调整。
在步骤48中,可变形模型30可以基于由模式检测所确定的图像特征并且基于如最初选择的可变形模型的平均网格而被进一步调整,以便使如上所述的总的能量E最小化。
应当理解,图像处理装置16能够被形成为在计算机上执行的计算机程序,并且方法40能够通过所述计算机程序来执行。
计算机程序可以被存储/分布在合适的(非瞬态)介质上,例如与其他硬件一起或作为其他硬件的部分供应的光学存储介质或固态介质,但是也可以被以其他形式分布,例如经由互联网或其他有线或无线的电信系统。此外,不同的实施例能够采取计算机程序产品的形式,所述计算机程序产品可从计算机可用介质或计算机可读介质存取,所述计算机可用介质或计算机可读介质提供用于由计算机或运行指令的任何设备或系统使用或者与计算机或运行指令的任何设备或系统结合使用的程序代码。就本公开内容的目的而言,计算机可用介质或计算机可读介质总体上能够是能够包含、存储、通信、传播或输送用于由指令运行设备使用或与指令运行设备结合使用的程序的任何有形装置。
此外,不同的实施例能够采取计算机程序产品的形式,所述计算机程序产品可从计算机可用介质或计算机可读介质存取,所述计算机可用介质或计算机可读介质提供用于由计算机或运行指令的任何设备或系统使用或者与计算机或运行指令的任何设备或系统结合使用的程序代码。就本公开内容的目的而言,计算机可用介质或计算机可读介质总体上能够是能够包含、存储、通信、传播或输送用于由指令运行设备使用或与指令运行设备结合使用的程序的任何有形装置。
至于公开内容的实施例已经被描述为至少部分通过软件控制的数据处理设备进行实施,应当认识到,非瞬态机器可读介质执行这样的软件,例如,光盘、磁盘、半导体存储器等,其也被认为表示本公开内容的实施例。
计算机可用介质或计算机可读介质能够例如是非限制性的,电学、磁性、光学、电磁、红外或半导体系统或传播介质。计算机可读介质的非限制性范例包括半导体或固态存储器、磁带、可移除计算机软盘、随机存取存储器(RAM)、只读存储器(ROM)、硬磁盘以及光盘。光盘可以包括压缩盘-只读存储器(CD-ROM)、压缩盘-读/写(CD-R/W)以及DVD。
另外,计算机可用介质或计算机可读介质可以包含或存储计算机可读程序代码或计算机可用程序代码,使得当所述计算机可读程序代码或所述计算机可用程序代码在计算机上运行时,这种计算机可读程序代码或计算机可用程序代码的运行令所述计算机通过通信链路发射另一计算机可读程序代码或计算机可用程序代码。这种通信链路可以使用介质,亦即,例如,非限制性地,是物理的或无线的介质。
适合于存储和/或运行计算机可读程序代码或计算机可用程序代码的数据处理系统或设备将包括通过诸如系统总线的通信结构直接或间接地耦合到存储器元件的一个或多个处理器。所述存储器元件可以包括在所述程序代码的实际运行期间采用的本地存储器、大容量存储设备以及缓存存储器,其提供对至少一些计算机可读程序代码或计算机可用程序代码的暂时存储,以降低可能在代码的运行期间从大容量存储设备检索代码的次数。
输入/输出或I/O设备能够被直接地耦合到所述系统或者通过介于中间的I/O控制器被耦合到所述系统。这些设备可以包括,例如,非限制性地,键盘、触摸屏显示器以及指向设备。不同的通信适配器也可以被耦合到所述系统,以使得数据处理系统能够通过介于中间的私有网络或公共网络变得耦合到其他数据处理系统、远程打印机或存储设备。非限制性的范例是调制解调器和网络适配器,并且仅仅是数种当前可用类型的通信适配器。
尽管已经在附图和前面的描述中详细图示和描述了本发明,但是这样的图示和描述应当被认为是图示性或示范性的,而非限制性的;本发明不限于所公开的实施例。本领域技术人员通过研究附图、公开内容以及权利要求,在实践请求保护的发明时能够理解并实现对所公开的实施例的其他变型。
在权利要求中,“包括”一词不排除其他元件或步骤,并且词语“一”或“一个”不排除多个。单个元件或其他单元可以实现在权利要求中记载的若干项的功能。尽管某些措施被记载在互不相同的从属权利要求中,但是这并不指示不能有利地使用这些措施的组合。
权利要求中的任何附图标记都不应被解释为对范围的限制。

Claims (17)

1.一种用于在目标(12)的多维图像数据中分割感兴趣区域(15)的图像处理装置(16),包括:
-接口(18),其用于接收包括要被分割的所述感兴趣区域的所述目标的图像数据,所述图像数据是经由超声探头的换能器阵列采集的,
-选择单元(20),其用于选择对应于所述图像数据中的所述目标的解剖结构的可变形模型(30),
-处理单元(22),其用于通过基于所述图像数据(xt)和所述目标的额外的解剖信息调整所述可变形模型来分割所述感兴趣区域,
其中,所述处理单元被设计为基于所述图像数据和从所述图像数据的视场导出的关于所述超声探头的位置的矢量(n)或平面来调整所述可变形模型的元素(c)的位置。
2.根据权利要求1所述的图像处理装置,其中,所述处理单元被设计为基于所述图像数据和所述目标的解剖特征在所述图像数据中的已知位置(xs)来调整所述可变形模型的所述元素(c)的所述位置。
3.根据权利要求1所述的图像处理装置,其中,所述处理单元被设计为基于所述图像数据和所述图像数据的视场(32)内的已知位置(34)来调整所述可变形模型的元素(c)的位置。
4.根据权利要求1所述的图像处理装置,其中,所述处理单元被设计为基于所述图像数据和所述图像数据的连续时间帧中的对应位置来调整所述可变形模型的元素(c)的位置。
5.根据权利要求4所述的图像处理装置,其中,所述对应位置是所述目标的已知的对应解剖特征。
6.根据权利要求2所述的图像处理装置,其中,所述处理单元被设计为基于所述元素的所述位置与所述已知位置之间的距离或者所述对应位置之间的距离来调整所述元素的所述位置。
7.根据权利要求3所述的图像处理装置,其中,所述处理单元被设计为基于所述元素的所述位置与所述已知位置之间的距离或者所述对应位置之间的距离来调整所述元素的所述位置。
8.根据权利要求6所述的图像处理装置,其中,所述可变形模型被调整为使包括从所述距离导出的弹性能量(ES)的模型能量(E)最小化。
9.根据权利要求8所述图像处理装置,其中,所述模型能量(E)包括从所述可变形模型与所述感兴趣区域中的对应图像特征的差异导出的外部能量(E外部)以及所述弹性能量(ES)。
10.根据权利要求9所述的图像处理装置,其中,所述弹性能量通过权重因子(w)进行加权。
11.根据权利要求1所述的图像处理装置,其中,所述处理单元被设计为基于所述目标的期望形状来调整所述可变形模型。
12.根据权利要求1所述的图像处理装置,其中,所述可变形模型是通过包括具有顶点(v)和边的多边形(T)的网格来形成的。
13.根据权利要求12所述的图像处理装置,其中,所述可变形模型是通过调节所述网格的所述顶点的位置来调整的。
14.一种用于在目标(12)的多维图像数据中分割感兴趣区域(15)的方法(40),包括如下步骤:
-接收(42)包括要被分割的所述感兴趣区域的所述目标的图像数据,所述图像数据是经由超声探头的换能器阵列采集的,
-选择(44)对应于所述图像数据中的所述目标的解剖结构的可变形模型(30),
-通过将所述可变形模型应用到所述图像数据来分割(46)所述感兴趣区域,并且
-基于所述图像数据(xt)和所述目标的额外的解剖信息调整所述可变形模型,
其中,基于所述图像数据和从所述图像数据的视场导出的关于所述超声探头的位置的矢量(n)或平面来调整所述可变形模型的元素(c)的位置。
15.根据权利要求14所述的方法(40),还包括
-基于所述图像数据(xt)和所述目标的额外的解剖信息并且基于所述目标的解剖特征在所述图像数据中的已知位置来调整(48)所述可变形模型的所述元素(c)的所述位置。
16.一种用于在目标(12)的多维图像数据中分割感兴趣区域(15)的计算机程序,所述计算机程序包括用于当所述计算机程序在计算机上被执行时令所述计算机执行如下步骤的程序代码单元:
-接收(42)包括要被分割的所述感兴趣区域的所述目标的图像数据,所述图像数据是经由超声探头的换能器阵列采集的,
-选择(44)对应于所述图像数据中的所述目标的解剖结构的可变形模型(30),
-通过将所述可变形模型应用到所述图像数据来分割(46)所述感兴趣区域,并且
-基于所述图像数据(xt)和所述目标的额外的解剖信息调整所述可变形模型,
其中,基于所述图像数据和从所述图像数据的视场导出的关于所述超声探头的位置的矢量(n)或平面来调整所述可变形模型的元素(c)的位置。
17.根据权利要求16所述的计算机程序,还包括用于当所述计算机程序在所述计算机上被执行时令所述计算机执行如下步骤的程序代码单元
-基于所述图像数据(xt)和所述目标的额外的解剖信息并且基于所述目标的解剖特征在所述图像数据中的已知位置来调整(48)所述可变形模型的所述元素(c)的所述位置。
CN201580015322.6A 2014-03-21 2015-02-27 用于分割感兴趣区域的图像处理装置和方法 Active CN106133789B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP14161050.1 2014-03-21
EP14161050 2014-03-21
PCT/EP2015/054189 WO2015139937A1 (en) 2014-03-21 2015-02-27 Image processing apparatus and method for segmenting a region of interest

Publications (2)

Publication Number Publication Date
CN106133789A CN106133789A (zh) 2016-11-16
CN106133789B true CN106133789B (zh) 2022-11-18

Family

ID=50396869

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580015322.6A Active CN106133789B (zh) 2014-03-21 2015-02-27 用于分割感兴趣区域的图像处理装置和方法

Country Status (5)

Country Link
US (1) US10043270B2 (zh)
EP (1) EP3120323B1 (zh)
JP (1) JP6400725B2 (zh)
CN (1) CN106133789B (zh)
WO (1) WO2015139937A1 (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018108923A1 (en) * 2016-12-15 2018-06-21 Koninklijke Philips N.V. X-ray apparatus having a composite field of view
CN106887000B (zh) * 2017-01-23 2021-01-08 上海联影医疗科技股份有限公司 医学图像的网格化处理方法及其系统
US10537261B2 (en) * 2017-05-10 2020-01-21 Boston Scientific Scimed Inc. Region-of-interest representations for electroanatomical mapping
EP3432262A1 (en) * 2017-07-18 2019-01-23 Koninklijke Philips N.V. Method and system for dynamic multi dimensional images of an object
US11006926B2 (en) 2018-02-27 2021-05-18 Siemens Medical Solutions Usa, Inc. Region of interest placement for quantitative ultrasound imaging
EP3683773A1 (en) * 2019-01-17 2020-07-22 Koninklijke Philips N.V. Method of visualising a dynamic anatomical structure
EP4012650A1 (en) 2020-12-14 2022-06-15 Koninklijke Philips N.V. Segmentation of anatomical structure in image

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1350681A (zh) * 1999-03-18 2002-05-22 纽约州立大学研究基金会 实行三维虚拟检查、导引行进和可视化的系统和方法
CN102202576A (zh) * 2008-10-10 2011-09-28 皇家飞利浦电子股份有限公司 用于为了在微创的x射线引导的介入中减少x射线剂量而产生覆盖经分割的目标结构或病变的减小视场的、具有自动快门适应的血管造影图像采集系统和方法
CN102203630A (zh) * 2008-11-05 2011-09-28 皇家飞利浦电子股份有限公司 Mr扫描的自动顺序规划
CN102711617A (zh) * 2009-11-16 2012-10-03 皇家飞利浦电子股份有限公司 扫描规划视场调整器、确定器和/或质量评估器

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10111661A1 (de) 2001-03-09 2002-09-12 Philips Corp Intellectual Pty Verfahren zum Segmentieren einer in einem Objekt enthaltenen dreidimensionalen Struktur, insbesondere für die medizinische Bildanalyse
WO2002073536A2 (en) * 2001-03-09 2002-09-19 Koninklijke Philips Electronics N.V. Image segmentation
US6961454B2 (en) * 2001-10-04 2005-11-01 Siemens Corporation Research, Inc. System and method for segmenting the left ventricle in a cardiac MR image
US7343365B2 (en) * 2002-02-20 2008-03-11 Microsoft Corporation Computer system architecture for automatic context associations
US20060008143A1 (en) 2002-10-16 2006-01-12 Roel Truyen Hierachical image segmentation
US7555151B2 (en) 2004-09-02 2009-06-30 Siemens Medical Solutions Usa, Inc. System and method for tracking anatomical structures in three dimensional images
US7672492B2 (en) 2005-01-31 2010-03-02 Siemens Medical Solutions Usa, Inc. Method of incorporating prior knowledge in level set segmentation of 3D complex structures
JP4918048B2 (ja) * 2005-02-11 2012-04-18 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 画像処理装置及び方法
US8320652B2 (en) * 2005-04-01 2012-11-27 Koninklijke Philips Electronics N.V. Method, a system and a computer program for segmenting a structure in a Dataset
US7545979B2 (en) * 2005-04-12 2009-06-09 General Electric Company Method and system for automatically segmenting organs from three dimensional computed tomography images
CN101203883B (zh) * 2005-06-21 2010-06-23 皇家飞利浦电子股份有限公司 渐进的基于模型的适配
US7876938B2 (en) * 2005-10-06 2011-01-25 Siemens Medical Solutions Usa, Inc. System and method for whole body landmark detection, segmentation and change quantification in digital images
CN101310305A (zh) * 2005-11-18 2008-11-19 皇家飞利浦电子股份有限公司 用于在3d图像中描绘预定结构的方法
JP5113077B2 (ja) * 2005-12-19 2013-01-09 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 変形可能なメッシュを用いて画像の後処理を容易にする方法
US8023703B2 (en) 2006-07-06 2011-09-20 The United States of America as represented by the Secretary of the Department of Health and Human Services, National Institues of Health Hybrid segmentation of anatomical structure
US10111726B2 (en) * 2007-04-26 2018-10-30 Koninklijke Philips N.V. Risk indication for surgical procedures
US8150116B2 (en) * 2007-07-02 2012-04-03 Siemens Corporation Method and system for detection of deformable structures in medical images
EP2446416B1 (en) * 2009-06-24 2014-04-30 Koninklijke Philips N.V. Establishing a contour of a structure based on image information
JP5736386B2 (ja) * 2009-12-10 2015-06-17 コーニンクレッカ フィリップス エヌ ヴェ 外傷性脳損傷の迅速かつ正確な定量的評価システム
WO2013028937A1 (en) * 2011-08-23 2013-02-28 Parmar Jaywant Philip Em guidance device for a device enabled for endovascular navigation placement including a remote operator capability and em endoluminal imaging technique
US10282846B2 (en) * 2013-12-10 2019-05-07 Koninklijke Philips N.V. Model-based segmentation of an anatomical structure
KR20150098119A (ko) * 2014-02-19 2015-08-27 삼성전자주식회사 의료 영상 내 거짓양성 병변후보 제거 시스템 및 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1350681A (zh) * 1999-03-18 2002-05-22 纽约州立大学研究基金会 实行三维虚拟检查、导引行进和可视化的系统和方法
CN102202576A (zh) * 2008-10-10 2011-09-28 皇家飞利浦电子股份有限公司 用于为了在微创的x射线引导的介入中减少x射线剂量而产生覆盖经分割的目标结构或病变的减小视场的、具有自动快门适应的血管造影图像采集系统和方法
CN102203630A (zh) * 2008-11-05 2011-09-28 皇家飞利浦电子股份有限公司 Mr扫描的自动顺序规划
CN102711617A (zh) * 2009-11-16 2012-10-03 皇家飞利浦电子股份有限公司 扫描规划视场调整器、确定器和/或质量评估器

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Medical Image Segmentation Using Minimal Path Deformable Models With Implicit Shape Priors》;P. Yan等;《IEEE Transactions on Information Technology in Biomedicine》;20061231;第10卷(第4期);第677-684页 *
《分割位置提示的可变形部件模型快速目标检测》;杨扬等;《自动化学报》;20120430;第38卷(第4期);第540-548页 *

Also Published As

Publication number Publication date
JP2017507754A (ja) 2017-03-23
EP3120323A1 (en) 2017-01-25
US10043270B2 (en) 2018-08-07
CN106133789A (zh) 2016-11-16
WO2015139937A1 (en) 2015-09-24
EP3120323B1 (en) 2020-05-20
US20170084023A1 (en) 2017-03-23
JP6400725B2 (ja) 2018-10-03

Similar Documents

Publication Publication Date Title
CN106133789B (zh) 用于分割感兴趣区域的图像处理装置和方法
KR101908520B1 (ko) 메디컬 이미징에서 공간 및 시간 제약들을 이용하는 랜드마크 검출
EP3100236B1 (en) Method and system for constructing personalized avatars using a parameterized deformable mesh
CN109589170B (zh) 医学成像中的左心耳闭合引导
US10102622B2 (en) Processing apparatus, processing method, and non-transitory computer-readable storage medium
US8582854B2 (en) Method and system for automatic coronary artery detection
KR101835873B1 (ko) 의료 이미지들의 분할 불확실성의 계산 및 시각화를 위한 시스템들 및 방법들
US20120253170A1 (en) Method and apparatus for generating medical image of body organ by using 3-d model
US20140241606A1 (en) Apparatus and method for lesion segmentation in medical image
JP2015500083A (ja) 心エコー検査のための自動的な画像化平面選択
US10762648B2 (en) Image processing apparatus, image processing method, image processing system, and program
US9886755B2 (en) Image processing device, imaging system, and image processing program
CN107787203B (zh) 图像配准
US20130346050A1 (en) Method and apparatus for determining focus of high-intensity focused ultrasound
US20190362548A1 (en) Apparatus and method for creating biological model
CN111563496A (zh) 用于图像采集的自动视图规划的连续学习
CN109863534B (zh) 用于分割解剖结构的二维图像的方法和装置
RU2505860C2 (ru) Одновременная основанная на модели сегментация объектов, удовлетворяющих заранее заданным пространственным соотношениям
JP6498316B2 (ja) 超音波撮像装置、および画像処理方法
US20220375099A1 (en) Segmentating a medical image
JP7165541B2 (ja) ボリュームデータ処理装置、方法及びプログラム
CN116152143A (zh) 剂量形变评估方法、系统及设备

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant