CN102812494A - 血管造影图像中的运动可视化 - Google Patents
血管造影图像中的运动可视化 Download PDFInfo
- Publication number
- CN102812494A CN102812494A CN2011800134316A CN201180013431A CN102812494A CN 102812494 A CN102812494 A CN 102812494A CN 2011800134316 A CN2011800134316 A CN 2011800134316A CN 201180013431 A CN201180013431 A CN 201180013431A CN 102812494 A CN102812494 A CN 102812494A
- Authority
- CN
- China
- Prior art keywords
- image
- angiogram
- motion
- subtraction image
- subtraction
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/254—Analysis of motion involving subtraction of images
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/503—Clinical applications involving diagnosis of heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4464—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit or the detector unit being mounted to ceiling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10121—Fluoroscopy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Abstract
描述了一种用于对血管造影图像中移动对象的运动可视化的方法和设备。在所述方法的优选实施例中,首先采集感兴趣对象的掩模图像,并且采集处在对象的运动的不同时相的所述对象的血管造影图像序列。然后,通过从掩模图像减去血管造影图像来生成第一血管造影相减图像和至少第二血管造影相减图像。随后,通过从第二血管造影相减图像减去第一血管造影相减图像来生成两次相减的图像。通过这种方式,执行双相减,即两次相减的血管造影,来方便对运动的估计。
Description
技术领域
本发明涉及血管造影图像中运动的可视化。具体而言,本发明涉及用于血管造影图像中的移动对象的运动可视化的方法和设备。此外,本发明还涉及当其在计算机上运行时能够控制这种方法的计算机程序单元以及涉及在其上存储有这样的计算机程序单元的计算机可读介质。
背景技术
在众多科学领域中,为观察者提供感兴趣对象的图像或图像序列非常重要。在许多应用中,感兴趣对象可能是移动的。尤其是在一些医学领域中,提供关于对象的运动的准确数据可能非常重要。
为了提供良好的数据,并且尤其为了提供对移动对象的良好可视化,可以对图像进行处理。例如,可以执行在JP 10031745A中所描述的数字剪影血管造影术。
然而,可能难以根据一序列图像对运动进行评估,在该运动快速、不规则和/或与图像的尺寸或周边对象的运动相比较小的情况下尤其如此。
发明内容
实现对感兴趣对象的运动的可视化的改善并由此实现对所述运动的更为可靠和准确的评估将是有利的。
根据本发明的第一方面,提供了一种用于血管造影图像中移动的感兴趣对象的运动可视化的方法。该方法包括如下步骤:采集感兴趣对象的背景的掩膜图像并在所述对象的不同运动阶段采集所述对象的血管造影图像序列;通过从掩膜图像减去血管造影图像来生成第一血管造影相减图像和至少第二血管造影相减图像;通过从第二血管造影相减图像减去第一血管造影相减图像来生成两次相减的图像。
换言之,本发明的第一方面可以被视为基于这样的构想:执行双相减,即两次相减的血管造影,以便于对运动的评估。首先,数字剪影血管造影(DSA)是通过从所采集的血管造影图像中的每幅中减去掩膜图像或相反的操作来执行的。从DSA得到的图像序列包括至少两幅图像,即第一血管造影相减图像和至少第二血管造影相减图像。DSA提供具有感兴趣对象(例如血管结构)的增强的表示的血管造影相减图像。在DSA之后,执行第二次相减,在此过程中,从在DSA中生成的血管造影相减图像中的一幅中减去在DSA中生成的至少一幅另一幅血管造影相减图像。即,从一幅或多幅第二血管造影相减图像减去例如第一血管造影相减图像。所述第二相减提供了两次相减的图像,其仅示出了两幅血管造影相减图像之间的差异。因此,仅有图像中对象的位移可见。在执行双相减之后,可以向观察者提供两次相减的图像。
血管造影血管图像序列的采集可以表示从数据库或存储器检索,或者可以表示通过特定的技术获得图像,所述技术例如是X射线成像,特别是旋转X射线成像、计算机断层摄影(CT)和磁共振(MR)。所述图像可以是数据、帧、帧序列或者其视觉表示。
移动对象还可以被表示为感兴趣对象。在血管造影图像中表示的移动对象可以是血管、心脏、冠状动脉、冠状静脉故障,特别是动脉瘤。
血管造影图像的序列包括以某种方式采集的至少两幅以及优选若干幅图像,使得所述感兴趣对象以增强的方式,例如在注射造影剂之后,被可视化。就周期性运动的对象而言,血管造影图像的序列可以覆盖若干运动周期。所述序列例如是一组或一系列图像,其在时间上彼此相随,并且可以作为电影来观察。
血管造影图像的序列表示处在运动的不同时相(phase)的对象。运动的时相可以与心跳相关或者与呼吸相关。在运动的时相取决于心跳的情况下,血管造影图像可以表示例如处在心脏循环的不同时相,例如处在收缩期和舒张期以及那两个极端时相之间的时相的血管或者动脉瘤。移动对象可以周期性、即脉动式地移动。
掩模图像可以是在诸如血管的感兴趣对象中没有出现造影剂的时间段期间采集的图像。在血管造影图像的采集期间,造影剂可以出现在血管中。所述掩模图像和血管造影图像显示出了相同的对象和感兴趣对象,即所述掩模图像是在造影剂出现在血管中之前采集的。所述血管造影图像可以在造影剂出现在血管中的时间段期间采集。所述掩模图像可以是血管造影图像序列的部分并且可以是在血管造影图像的采集之前采集的。所述掩模图像可以从血管造影图像中减去或者备选地所述血管造影图像可以从掩模图像中减去。
在周期性运动的情况下,掩模图像可以是针对运动的每个阶段采集的,并且从分别在相同阶段采集的对应血管造影图像中减去。
从血管造影图像序列减去掩模图像提供了第一血管造影相减图像和至少一个第二血管相减图像,其中,例如,血管被加强并且诸如骨骼和致密组织的非感兴趣结构被去除。所述掩模图像可以从所述序列中的每幅血管造影图像减去。在第二相减中,通过将在第一相减中接收的血管造影相减图像彼此相减来生成至少一幅两次相减的图像。例如,预确定的图像,例如固定图像,即所述第一血管造影相减图像可以表示对应于心脏收缩的运动状态的感兴趣对象。这种第一血管造影相减图像可以从表示诸如心脏舒张的运动的不同时相的第二血管造影相减图像中减去。所述第一血管相减图像可以从所有其他血管造影相减图像,即第二血管造影相减图像中减去,所述其他血管造影相减图像包括在已经从中减去掩模图像的血管造影图像的序列中。所得到的一幅或多幅两次相减的图像可以仅示出对象的位移。当对象的周边与所述对象相比不移动并且在第一血管造影相减图像和第二血管造影相减图像中是相同的时,尤其是这种情况,。或者,对象的周边结构的运动与对象的运动相比小到可以忽略不计。第一相减和第二相减可以同时执行或者优选相继执行。
通过生成两次相减的图像,可以使得即使小对象的小的运动的检测和评估也成为可能,并且这种运动的评估可以为观察者提供方便。
根据本发明的实施例,所述方法包括还包括对表示处在运动的相同时相的对象的血管造影图像求平均。
就周期性移动对象而言,血管造影图像的序列可以包括表示处于阶段的相同时相的对象的若干幅图像,如果所述序列包括运动的若干次循环的话。例如,就心跳而言,血管造影图像的序列可以包括若干次心跳并且由此包括表示例如心脏收缩的若干幅血管造影图像以及表示心脏舒张的若干血管造影图像。血管造影图像的序列的图像可以首先例如根据其心脏时相进行归类并且然后相同时相的图像可以被求平均。因此,在求平均步骤之后例如表示若干心脏周期的血管造影图像序列仅表示一个心脏周期。
所述求平均通过降低噪声以及此外通过例如减小出现在血管中的造影剂浓度的不规则可以增强所生成的图像的质量和准确度。
根据本发明的第二方面,提供了一种用于血管造影图像中移动对象的运动可视化的设备。所述设备包括第一减法器,其用于通过从掩模图像减去血管造影图像来生成第一血管造影相减图像和至少第二血管造影相减图像。所述设备还包括第二减法器,其用于通过从第二血管造影相减图像减去第一血管造影相减图像来生成两次相减的图像。
所述设备还可以包括图像检测器或图像检索器,其用于采集处在运动的不同时相的对象的图像。此外,所述设备可以包括例如用于提供第二相减的结果,即两次相减的图像的显示器或屏幕,使得由于所述对象的运动造成的差异能够被可视化。所述第一减法器和所述第二减法器可以是所述设备的单元,所述单元被彼此连接和/或连接到存储器和处理器。
根据本发明的第三方面,提供了一种计算机程序单元,当所述计算机程序单元在计算机上执行时,能够控制上文所描述的方法。
根据本发明的第四方面,提供了一种具有上文所述的计算机程序单元的计算机可读介质。
应当注意到,参考不同的主题描述了本发明的实施例。具体而言,一些实施例是参考方法类型的权利要求描述的,而其他实施例是参考设备类型的权利要求加以描述的。然而,本领域技术人员将从上文和下文的描述中认识到,除非另有说明,除了属于一种类型的主题的特征的组合之外,属于不同主题的特征之间,具体而言,设备类型权利要求的特征与方法类型权利要求的特征之间的任意组合,也被认为在本发明中公开了。
从上下文描述的实施例的范例能够导出本发明的上文限定的各方面和其他方面、特征和优点。将参考实施例的范例在下文更为详细的描述本发明,但其并非限制本发明。
附图说明
图1示出了示意性表示用于根据本发明的实施例进行运动可视化的方法的流程图;
图2示出了示意性表示用于根据本发明的另一实施例进行运动可视化的方法的流程图;
图3A示意性示出了要在本发明的实施例中使用的表示处在心脏循环的不同时相的动脉瘤的经平均的血管相减图像;
图3B示意性示出了根据本发明的实施例生成的在图3A中的动脉瘤的两次相减的图像;
图4A和4B示出了类似于图3A和3B的图像序列的另一范例;
图5A和5B示出了类似于图3A和3B的图像序列的另一范例;
图6示意性示出了根据本发明的实施例在血管造影图像中对运动对象进行运动可视化的设备。
附图中的图示仅仅是示意性的。应当认识到,在不同的附图中,相似或相同的要素被赋予相同的附图标记。
附图标记列表
1 对象
3 表示处在心脏循环的不同时相的动脉瘤的经平均的血管造影图
像的序列
5 动脉瘤的两次相减图像的序列
7 显示器
9 X射线探测器
11 用于运动可视化的设备
13 第一减法器
15 第二减法器
17 处理器
a-e 运动的不同时相
S1 采集对象的掩膜图像
S3 采集处在对象的运动的不同时相的所述对象的血管造影图像的
序列
S5 通过从掩膜图像减去血管造影图像来生成第一血管造影相减图
像和至少第二血管造影相减图像
S7a 通过从第二血管造影相减图像减去第一血管造影相减图像来生
成两次相减的图像
S7b 从若干第二血管造影相减图像减去第一血管造影相减图像来生
成若干两次相减的图像
S9 选择其中对象被填充造影剂的血管造影图像
S11 确定在至少两幅血管造影图像中的运动的时相
S13 选择第一血管造影相减图像,使得在第一血管造影相减图像中表
示的运动的阶段是运动的极端时相
S15 配准血管造影图像
S17 对表示处在运动的相同时相的对象的血管造影图像求平均
S19 显示两次相减的图像
具体实施方式
在下文中,经受运动的对象的实施例的范例可以是动脉瘤。在动脉瘤的处置可能的情况下,应当评估在处置期间破裂的风险是否较高或者是否应当对动脉瘤不进行处置。动脉瘤壁运动可能是指示破裂的风险的一个因素。动脉瘤壁运动与动脉瘤的局部生物机械属性相关,其继而与组织经受的损伤的量相关。如果动脉瘤壁的一部分比另一部分移动更多,那里可能有壁的薄弱点,其可能增加破裂的风险。在动脉瘤壁不移动的情况下,这可能指示破裂的低风险。因此,其对于使动脉瘤壁的运动可视化以对状况进行评估可能是很重要的。
动脉瘤可以借助血管造影图像来可视化。根据血管造影图像的序列对动脉瘤壁运动的评估可能是困难的。例如,所述血管造影图像的序列可以以每秒5帧或图像的速率进行采集。优选地,所述血管造影序列是以大约每秒30帧或图像的高速率采集的。如果这样的序列以正常的速度显示,所述运动被显示为闪动但难以评估。如果所述序列是以非常低的速度显示的,所述运动不再可见,因为所述运动过小。此外,需要从所述序列选择帧以增强所述表示的质量的方法要求进行用户交互,并且因此是不精确的并且非常耗时。
图1和图2中所示的方法的实施例可以帮助提供更好的运动可视化并且由此便于对例如动脉瘤壁运动进行估计。
图1示出了示意性表示用于根据本发明的实施例的运动可视化的流程图。在步骤S 1中,采集例如动脉瘤的对象的掩模图像。在步骤S3中,采集处在对象的运动的不同时相的所述对象的血管造影图像的序列。掩模图像的采集发生在造影剂出现在对象中之前。在造影剂出现在对象体内期间采集血管造影图像的序列。掩模图像的采集还可以是步骤S3的部分。或者,步骤S1和步骤S3可以同时执行或相继执行。在步骤S1和S3的同时执行过程中,所述掩模图像可以在血管造影图像序列的采集期间从数据库检索。
在步骤S5中,执行第一相减。其中,通过从掩模图像减去血管造影图像或者通过从血管造影图像减去掩模图像来生成第一血管造影相减图像和至少第二血管造影相减。所述第一血管造影相减图像和第二血管造影相减图像表示对象的运动的不同时相。此外,运动的时相可以表示对象的周期性运动的时相。
运动的时相例如可以是诸如心脏收缩或心脏舒张的心脏时相。此外,运动的时相例如可以是诸如吸气或呼气的呼吸时相。此外,运动的时相可以是诸如心脏跳动和呼吸运动的组合的若干种运动的叠加。
在步骤S7a和S7b中,执行第二相减。在步骤S7a中,通过从第二血管造影相减图像减去第一血管造影相减图像来生成两次相减的图像。或者,在步骤S7b中,通过从若干、优选所有第二血管造影相减图像减去第一血管造影相减图像来生成若干两次相减的图像。例如,根据步骤S7a,两次相减的图像可以表示心脏收缩和心脏舒张运动时相的差别。根据步骤S7b,所述第二相减可以得到两次相减的图像的序列:从其他血管造影相减图像,即第二血管造影相减图像,减去一幅能预定义的图像,例如固定的图像,即表示一个心脏时相的第一血管造影相减图像。
所述方法的步骤可以是自动执行的。自动可以表示没有用户交互的必要。所述方法的自动化可以提供质量更好的并且更为准确的数据,其可以比利用需要用户交互的方法生成的数据有具更好的重复性。
所述移动对象可以是动脉瘤、血管结构、心脏或血管畸形体。此外,所述运动可以是由心脏跳动引起的并且由心脏循环的时相表示,和/或其可以是由呼吸引起的并且由呼吸循环的时相表示。
在图2中,呈现了示意性表示根据本发明的另一实施例进行运动可视化的方法的流程图。
步骤S1、S3、S5、S7a和S7b是类似于图1中所示的范例而执行的。在第一相减之后,在步骤S5中,一些其他的步骤S9至S17可以包括在所述方法中。此外,步骤S19可以跟随步骤S7a和/或S7b。
在步骤S5中生成DSA的血管造影相减图像之后,选择在其中对象(即动脉瘤)被填充造影剂的血管造影相减图像(步骤S9)。或者,在步骤S5中执行第一相减之前,选择血管造影图像的序列中在血管中存在造影剂的图像。即,步骤S9可以在步骤S3之后执行。步骤S9可以增强所得的两次相减的图像的质量,因为例如噪声可以通过在无需造影剂的情况下挑选图像来减少。
在步骤S11中,确定在DSA之后接收的至少两幅血管造影相减图像中的运动的时相。优选地,确定在所有血管造影相减图像中的运动的时相。例如,在步骤S11中,血管造影相减图像可以根据其心脏阶段进行归类。所述归类可以是基于心电图(基于ECG)或者是基于图像的。为了确定运动的时相,采集表示运动的运动信号,如例如ECG信号或呼吸信号。所述ECG信号可以允许对运动阶段的精确确定,并且由此,能够实现对血管造影相减图像的精确归类。ECG信号可以并行于步骤S1和S3来采集,或者可以备选地从数据库中检索。并行于步骤S1的ECG信号的采集可以实现对若干掩模图像的归类,使得例如从在心脏收缩的时相中采集的血管造影相减图像减去在心脏收缩的时相中采集的对应的掩模图像。
或者,基于对血管造影相减图像自身的分析确定血管造影相减图像中运动的时相。基于图像的归类可能需要较少的复杂设备。
步骤S11备选地可以在步骤S3之后并且在步骤S5之前执行。即所述归类可以基于在步骤S5中在DSA之前的血管造影图像的序列。此外,可以基于对在步骤S5中在DSA之前的血管造影图像的分析确定运动的时相。
在步骤S11中确定运动的时相之后,在步骤S13中选择第一血管造影相减图像。可以选择所述第一血管造影相减图像,使得在第一血管造影相减图像中表示的运动的时相是运动的极端时相。运动的极端时相可以是得到对象的最多或最少变形的运动的时相。就心脏周期而言,运动的极端时相可以是心脏收缩和心脏舒张。选择表示运动的极端时相的第一血管造影相减图像可以帮助增强运动可视化的质量,因为在运动的极端时相中,对象的变形与所述序列的其他血管造影图像相比可能是最大的。
此外,在步骤S15中,血管造影图像可以在步骤S03之后被刚性配准。即确定至少两幅血管造影图像之间的空间对应关系。换言之,通过旋转和平移图像之一而对其进行改变以对准所述图像。或者,配准步骤S15可以利用步骤S5之后的血管造影相减图像来实施。即,在步骤S5中的第一相减之后,可以配准所述第一血管造影相减图像和所述第二血管造影相减图像。
在移动感兴趣对象为动脉瘤的情况下,已知两种类型的脉动式运动。动脉瘤以及还有父动脉的脉动式扩张和脉动式摆动是可能的。这两种类型的运动都能够根据本发明在可视化中被增强。在仅脉动式扩张要被可视化的情况下,表示运动的不同时相的血管造影图像在步骤S7a和S7b中的相减之前利用刚性配准与彼此配准。其中,在刚性配准中,通过相对于彼此补偿图像的平移和旋转来对准两幅图像。刚性配准可能是有利的,因为其从图像中去除要被可视化的运动,由于避免了在配准期间图像的变形。
在步骤S17中,执行对表示处在运动的相同时相的对象的血管造影相减图像的求平均。在步骤S17中的求平均可以基于在步骤S11中执行的对运动的时相的确定。例如,相同心脏时相的图像被平均,以便降低造影剂浓度变化的影响以及降低噪声。因此,在步骤S17中,生成血管造影相减图像的经平均的序列,其表示一个心脏周期。处在心脏周期的不同时相的动脉瘤的血管造影相减图像的经平均的序列在图3A、4A和5A中被示出,并且以附图标记3来表示。
或者,求平均步骤S17可以在步骤S11中的第一相减之前执行。即。所述血管造影图像可以被平均。
在求平均步骤S17之后,执行步骤S7a或步骤S7b。其中,从可以作为运动的不同时相的血管造影相减图像的求平均结果的第二血管造影相减图像中减去可以作为表示运动的相同时相的血管造影相减图像的求平均结果的第一血管造影相减图像(步骤S7a)。或者,在步骤S7b中,从每幅第二血管造影相减图像中减去作为求平均结果的第一血管造影相减图像,其中,所述第二血管造影相减图像表示来自运动的其他时相的求平均结果。
换言之,选择第一血管造影相减图像,例如固定图像,例如表示结束心脏收缩时相。从血管造影相减图像的经平均的序列减去这一固定图像。这一相减的结果是两次相减的血管造影图像序列。
或者,在步骤S13中对表示运动的极端时相的第一血管造影相减图像的选择步骤可以在求平均步骤S17之后以及在步骤S7a和S7b之前执行。动脉瘤的两次相减的图像序列的范例在图3B、4B和5B中被示出,并且利用附图标记5来表示。
在步骤S19中,显示两次相减的图像。例如能够在执行相应的步骤之后或者备选地在执行所有步骤之后立即提供,即显示第一相减的结果、求平均结果和第二相减的结果。所提到的结果可以在不同的显示器上被可视化或者备选地在同一显示器上被可视化。此外,它们可以相继或者同时地被可视化。
在图3、4和5中,示出了血管造影序列的可视化的范例。其中,在图3A、4A和5A中,在执行第一相减(步骤S5)和求平均(步骤S17)之后呈现了经平均的血管造影图像。
在图3B、4B和5B中,呈现了在执行步骤S7b之后的两次相减图像序列的可视化。在此处为动脉瘤的对象1被示出为处在运动的不同时相(a、b、c、d、e)。在图3A、4A和5A中示出的经平均的DSA序列中,没有动脉瘤的壁运动可见。在第二相减之后,动脉瘤壁运动在图3B、4B、5B中变得可见。在图5中示出的动脉瘤没有显示出运动,如可以在图5B中的描绘中可见的。对象的运动的缺失是在图5B中的第二相减之后的“空白”表示的原因。
在图3A中的范例中,从表示运动阶段b、c、d和e的第二血管造影相减图像减去表示结束心脏舒张阶段(a)的血管造影相减图像。如图3B中所示,第一图像示出了“空白”表示,因为第一血管造影相减图像被从这一表示自身中被减去。当第一血管造影相减图像被从处在心脏收缩阶段(e)中的第二血管造影相减图像减去时,最大的差异被可视化。
在图6中,根据本发明的实施例,示意性示出了用于在血管造影图像中移动对象的运动可视化的设备11。
设备11包括具有可以同时用于采集所述对象的血管造影序列的X射线探测器9和数据存贮装置的C形臂。其中,所述数据存储装置可以集成在处理器17中。设备11还包括第一减法器13,其用于通过从掩膜图像减去血管造影图像来生成第一血管造影相减图像和至少第二血管造影相减图像。
此外,设备11包括第二减法器15,其用于通过从第二血管造影相减图像减去第一血管造影相减图像来生成两次相减的图像。第一减法器13和第二减法器15可以是处理器17的一部分或者可以与处理器17连接。此外,设备11可以包括一个或若干个显示器7,其用于提供第二相减的结果,使得所述对象的运动的差异被可视化。处理器17、显示器7和X摄像探测器9可以通线缆或无线地连接。
应当认识到,术语“包括”、“包含”等并不排除其他元件或步骤,并且“一”、“一个”并不排除复数。还可以组合结合不同的实施例描述的单元。应当认识到,在权利要求中的参考标记不应当被解释为显示权利要求的范围。此外,计算机程序可以存储和/或分布在合适的介质上,诸如光学存储介质或固态介质,其与其他硬件一起提供或者作为其他硬件的部分,但所述计算机程序也可以以其他形式分布,诸如经由因特网或其他有线或无线电信系统。
Claims (15)
1.一种用于血管造影图像中的移动对象(1)的运动可视化的方法,所述方法包括:
采集所述对象的掩膜图像(S1)以及处在所述对象的运动的不同时相的所述对象的血管造影图像的序列(S3);
通过从所述掩膜图像减去所述血管造影图像来生成第一血管造影相减图像和至少第二血管造影相减图像(S5);
通过从所述第二血管造影相减图像减去所述第一血管造影相减图像来生成两次相减的图像(S7a)。
2.根据权利要求1所述的方法,还包括:
对表示处在运动的相同时相的所述对象的血管造影相减图像求平均(S17)。
3.根据权利要求1所述的方法,还包括如下步骤:
选择所述第一血管造影相减图像,使得在所述第一血管造影相减图像中所表示的运动的所述时相是运动的极端时相(S13)。
4.根据权利要求1所述的方法,
其中,所述第一血管造影相减图像和所述第二血管造影相减图像表示所述对象的运动的不同时相。
5.根据权利要求1所述的方法,
其中,运动的所述时相表示所述对象的周期性运动的时相。
6.根据权利要求1所述的方法,还包括:
在至少两幅血管造影相减图像中确定运动的所述时相(S11)。
7.根据权利要求6所述的方法,还包括:
采集表示所述对象的运动的所述时相的运动信号;
其中,基于所述运动信号确定在所述血管造影相加图像中的运动的所述时相。
8.根据权利要求6所述的方法,
其中,基于对所述血管造影相减图像的分析确定所述血管造影相减图像中的运动的所述时相。
9.根据权利要求1所述的方法,还包括:
选择其中在所述对象体内存在造影剂的血管造影相减图像(S9)。
10.根据权利要求1所述的方法,
其中,所述方法的步骤是自动执行的。
11.根据权利要求1所述的方法,还包括:
确定所述血管造影图像中的至少两幅之间的空间对应关系(S15)。
12.根据权利要求1所述的方法,
其中,所述移动对象是动脉瘤、血管结构、血管畸形体和心脏之一;
其中,运动的所述时相表示心脏周期的时相和/或呼吸周期的时相。
13.一种用于血管造影图像中的移动对象(1)的运动可视化的设备(11),所述设备包括:
第一减法器(13),其用于通过从掩模图像减去血管造影图像来生成第一血管造影相减图像和至少第二血管造影相减图像;
第二减法器(15),其用于通过从所述第二血管造影相减图像减去所述第一血管造影相减图像来生成两次相减的图像。
14.一种计算机程序单元,当所述计算机程序单元在计算机上执行时,能够实现对根据权利要求1所述的方法的控制。
15.一种计算机可读介质,其具有根据权利要求14所述的计算机程序单元。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP10156297.3 | 2010-03-12 | ||
EP10156297 | 2010-03-12 | ||
PCT/IB2011/050926 WO2011110987A1 (en) | 2010-03-12 | 2011-03-04 | Motion visualisation in angiographic images |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102812494A true CN102812494A (zh) | 2012-12-05 |
CN102812494B CN102812494B (zh) | 2016-03-09 |
Family
ID=44247942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180013431.6A Active CN102812494B (zh) | 2010-03-12 | 2011-03-04 | 血管造影图像中的运动可视化 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9025840B2 (zh) |
EP (1) | EP2545528B1 (zh) |
JP (1) | JP6043187B2 (zh) |
CN (1) | CN102812494B (zh) |
WO (1) | WO2011110987A1 (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105264569A (zh) * | 2013-06-05 | 2016-01-20 | 爱克发医疗保健公司 | Cta容积中的自动化主动脉检测 |
CN106033603A (zh) * | 2015-03-10 | 2016-10-19 | 东芝医疗系统株式会社 | 医用图像处理装置和医用图像处理方法 |
CN107004270A (zh) * | 2014-12-17 | 2017-08-01 | 皇家飞利浦有限公司 | 用于计算感兴趣对象的位移的方法和系统 |
CN106999129A (zh) * | 2014-11-25 | 2017-08-01 | 皇家飞利浦有限公司 | 数字减影血管造影术 |
CN108364339A (zh) * | 2017-01-27 | 2018-08-03 | 西门子保健有限责任公司 | 对血管造影数据组进行图像处理的方法以及图像处理装置 |
CN109166097A (zh) * | 2018-07-18 | 2019-01-08 | 北京理工大学 | 从造影图像序列中提取血管的方法及装置 |
CN112006704A (zh) * | 2019-05-29 | 2020-12-01 | 西门子医疗有限公司 | 用于创建3d dsa图像的方法和c形臂x射线设备 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140092437A (ko) * | 2012-12-27 | 2014-07-24 | 삼성전자주식회사 | 엑스선 영상 장치 및 엑스선 영상 장치의 제어 방법 |
KR20150027881A (ko) * | 2013-08-29 | 2015-03-13 | 삼성전자주식회사 | 엑스선 영상 장치 및 그 제어방법 |
KR20150039242A (ko) * | 2013-10-01 | 2015-04-10 | 삼성전자주식회사 | 엑스선 영상 장치 및 그 제어 방법 |
CN107106102B (zh) * | 2015-01-05 | 2020-11-17 | 皇家飞利浦有限公司 | 数字减影血管造影 |
KR20160089762A (ko) | 2015-01-20 | 2016-07-28 | 삼성전자주식회사 | 방사선 촬영 장치 및 방사선 촬영 장치의 제어 방법 |
US10116863B2 (en) * | 2016-06-14 | 2018-10-30 | Goodrich Corporation | Scanning with frame averaging |
JP7233850B2 (ja) * | 2018-04-25 | 2023-03-07 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置、x線診断装置及び医用画像処理プログラム |
US11399791B2 (en) * | 2020-04-01 | 2022-08-02 | Wisconsin Alumni Research Foundation | System and method for flow-resolved three-dimensional imaging |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1305354A (zh) * | 1998-06-08 | 2001-07-25 | 阿米兰·格林瓦尔德 | 用于对血管中各个红血球的运动成像和分析的系统和方法 |
JP2006271484A (ja) * | 2005-03-28 | 2006-10-12 | Toshiba Corp | 画像形成装置及び画像形成方法 |
US7545967B1 (en) * | 2002-09-18 | 2009-06-09 | Cornell Research Foundation Inc. | System and method for generating composite subtraction images for magnetic resonance imaging |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59123385A (ja) * | 1982-12-29 | 1984-07-17 | Shimadzu Corp | X線画像処理装置 |
JPS59155234A (ja) * | 1983-02-23 | 1984-09-04 | 株式会社東芝 | 画像入力装置 |
JPS60142688A (ja) * | 1983-12-29 | 1985-07-27 | Shimadzu Corp | デイジタルサブトラクシヨンシステム |
JPS60150730A (ja) * | 1984-01-19 | 1985-08-08 | 株式会社東芝 | 心機能診断装置 |
JPH02202183A (ja) * | 1989-01-31 | 1990-08-10 | Shimadzu Corp | デジタルサブトラクションアンギオグラフィ装置 |
JPH04307049A (ja) * | 1991-04-03 | 1992-10-29 | Koonan:Kk | 心臓その他脈動体画像の解析装置 |
DE19613342A1 (de) | 1996-04-03 | 1997-10-09 | Philips Patentverwaltung | Automatisches Bildauswertungsverfahren |
JP2006110208A (ja) * | 2004-10-15 | 2006-04-27 | Shimadzu Corp | 医用画像表示装置 |
EP1817743B1 (en) * | 2004-11-24 | 2008-05-14 | Koninklijke Philips Electronics N.V. | Multi-feature time filtering for enhancing structures in noisy images |
DE102005032974B4 (de) * | 2005-07-14 | 2013-11-07 | Siemens Aktiengesellschaft | Verfahren zur 3D-Visualisierung vaskulärer Einsatzstücke im menschlichen Körper mit dem C-Bogen |
DE102005036564A1 (de) * | 2005-08-03 | 2007-02-22 | Siemens Ag | Betriebsverfahren für eine bildgebende medizintechnische Anlage und hiermit korrespondierende Gegenstände |
JP5085031B2 (ja) * | 2005-11-10 | 2012-11-28 | 株式会社東芝 | X線アンギオ撮影装置 |
US7986822B2 (en) * | 2006-04-24 | 2011-07-26 | Siemens Medical Solutions Usa, Inc. | System and method for X-ray based assessment of aneurysm pulsation |
JP5100145B2 (ja) * | 2006-06-20 | 2012-12-19 | 株式会社東芝 | 磁気共鳴イメージング装置 |
WO2008053401A1 (en) | 2006-11-03 | 2008-05-08 | Philips Intellectual Property & Standards Gmbh | Cardiac phase determination |
JP5536974B2 (ja) * | 2006-11-08 | 2014-07-02 | 株式会社東芝 | X線診断装置及び画像処理装置 |
JP5215050B2 (ja) | 2008-06-10 | 2013-06-19 | 株式会社東芝 | 画像処理装置 |
RU2530665C2 (ru) * | 2008-08-13 | 2014-10-10 | Кониклейке Филипс Электроникс Н.В. | Динамическая визуализация информации о состоянии коронарных сосудов и перфузии миокарда |
RU2519353C2 (ru) * | 2008-08-13 | 2014-06-10 | Конинклейке Филипс Электроникс Н.В. | Создание маски для кардиосубстракции |
-
2011
- 2011-03-04 US US13/582,124 patent/US9025840B2/en active Active
- 2011-03-04 CN CN201180013431.6A patent/CN102812494B/zh active Active
- 2011-03-04 WO PCT/IB2011/050926 patent/WO2011110987A1/en active Application Filing
- 2011-03-04 JP JP2012556627A patent/JP6043187B2/ja active Active
- 2011-03-04 EP EP11714426.1A patent/EP2545528B1/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1305354A (zh) * | 1998-06-08 | 2001-07-25 | 阿米兰·格林瓦尔德 | 用于对血管中各个红血球的运动成像和分析的系统和方法 |
US7545967B1 (en) * | 2002-09-18 | 2009-06-09 | Cornell Research Foundation Inc. | System and method for generating composite subtraction images for magnetic resonance imaging |
JP2006271484A (ja) * | 2005-03-28 | 2006-10-12 | Toshiba Corp | 画像形成装置及び画像形成方法 |
Non-Patent Citations (1)
Title |
---|
崔彦: "序列图像中运动目标检测与跟踪技术研究", 《中国优秀硕士学位论文全文数据库》, 15 July 2009 (2009-07-15), pages 2 - 1 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105264569A (zh) * | 2013-06-05 | 2016-01-20 | 爱克发医疗保健公司 | Cta容积中的自动化主动脉检测 |
CN105264569B (zh) * | 2013-06-05 | 2019-06-04 | 爱克发医疗保健公司 | Cta容积中的自动化主动脉检测 |
CN106999129A (zh) * | 2014-11-25 | 2017-08-01 | 皇家飞利浦有限公司 | 数字减影血管造影术 |
CN106999129B (zh) * | 2014-11-25 | 2020-11-03 | 皇家飞利浦有限公司 | 数字减影血管造影术 |
CN107004270A (zh) * | 2014-12-17 | 2017-08-01 | 皇家飞利浦有限公司 | 用于计算感兴趣对象的位移的方法和系统 |
CN107004270B (zh) * | 2014-12-17 | 2021-02-19 | 皇家飞利浦有限公司 | 用于计算感兴趣对象的位移的方法和系统 |
CN106033603A (zh) * | 2015-03-10 | 2016-10-19 | 东芝医疗系统株式会社 | 医用图像处理装置和医用图像处理方法 |
CN106033603B (zh) * | 2015-03-10 | 2019-06-11 | 东芝医疗系统株式会社 | 医用图像处理装置和医用图像处理方法 |
CN108364339A (zh) * | 2017-01-27 | 2018-08-03 | 西门子保健有限责任公司 | 对血管造影数据组进行图像处理的方法以及图像处理装置 |
CN108364339B (zh) * | 2017-01-27 | 2021-12-03 | 西门子保健有限责任公司 | 对血管造影数据组进行图像处理的方法以及图像处理装置 |
CN109166097A (zh) * | 2018-07-18 | 2019-01-08 | 北京理工大学 | 从造影图像序列中提取血管的方法及装置 |
CN109166097B (zh) * | 2018-07-18 | 2021-04-20 | 北京理工大学 | 从造影图像序列中提取血管的方法及装置 |
CN112006704A (zh) * | 2019-05-29 | 2020-12-01 | 西门子医疗有限公司 | 用于创建3d dsa图像的方法和c形臂x射线设备 |
CN112006704B (zh) * | 2019-05-29 | 2024-04-23 | 西门子医疗有限公司 | 用于创建3d dsa图像的方法和c形臂x射线设备 |
Also Published As
Publication number | Publication date |
---|---|
US9025840B2 (en) | 2015-05-05 |
EP2545528A1 (en) | 2013-01-16 |
WO2011110987A1 (en) | 2011-09-15 |
JP2013521849A (ja) | 2013-06-13 |
CN102812494B (zh) | 2016-03-09 |
US20120321156A1 (en) | 2012-12-20 |
JP6043187B2 (ja) | 2016-12-14 |
EP2545528B1 (en) | 2021-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102812494B (zh) | 血管造影图像中的运动可视化 | |
JP6956817B2 (ja) | 画像処理装置及びx線診断装置 | |
JP4468353B2 (ja) | 管状組織の三次元モデリングのための方法 | |
JP6013012B2 (ja) | 画像処理装置及びx線診断装置 | |
US8781257B2 (en) | Synchronization of two image sequences of a periodically moving object | |
JP6021420B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
US20080240536A1 (en) | Method of detection and compensation for respiratory motion in radiography cardiac images synchronized with an electrocardiogram signal | |
CN101242787A (zh) | 用盲去卷积空间增强具有噪声的图像中的结构的系统和方法 | |
US10410347B2 (en) | Digital subtraction angiography | |
US20090149741A1 (en) | Motion correction for tomographic medical image data of a patient | |
CN101663691B (zh) | 不同预采集医学图像的时空扭曲 | |
JP2011530365A (ja) | 心臓サブトラクションのためのマスク構造 | |
EP3359041B1 (en) | Mobile ffr simulation | |
JP2012000135A (ja) | マルチモダリティ動画像診断装置 | |
JP5215050B2 (ja) | 画像処理装置 | |
US20160345926A1 (en) | Acquisition of projection data for motion-corrected computed tomography images | |
JP2001148005A (ja) | 動いている対象物の三次元画像を再構成する方法 | |
JP2006087631A (ja) | 画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体 | |
JP2008228829A (ja) | 周期的運動の同期算出装置および周期的運動の同期算出方法 | |
JP2015097590A (ja) | 画像処理装置及びその作動方法、撮影装置及びその作動方法、撮影システム、並びにコンピュータプログラム | |
JP2012000323A (ja) | 3次元画像形成装置、3次元画像形成プログラム、及び3次元画像形成方法 | |
Schmitt et al. | Spatially resolved automatic cardiac rest phase determination in coronary computed tomography angiography (CTA) | |
JP2010088660A (ja) | 医療画像処理装置、医療画像処理プログラムおよび医療撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |