CN101528133B - 心脏相位确定 - Google Patents

心脏相位确定 Download PDF

Info

Publication number
CN101528133B
CN101528133B CN2007800404234A CN200780040423A CN101528133B CN 101528133 B CN101528133 B CN 101528133B CN 2007800404234 A CN2007800404234 A CN 2007800404234A CN 200780040423 A CN200780040423 A CN 200780040423A CN 101528133 B CN101528133 B CN 101528133B
Authority
CN
China
Prior art keywords
image
stage
subtraction image
projected
series
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.)
Expired - Fee Related
Application number
CN2007800404234A
Other languages
English (en)
Other versions
CN101528133A (zh
Inventor
B·莫瓦萨吉
M·格拉斯
D·舍费尔
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 Electronics 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 Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN101528133A publication Critical patent/CN101528133A/zh
Application granted granted Critical
Publication of CN101528133B publication Critical patent/CN101528133B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • 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
    • 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/027—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
    • 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/54—Control of apparatus or devices for radiation diagnosis
    • A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
    • 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

Abstract

为了根据旋转冠状动脉造影数据重建冠状动脉，一个关键点是选择用于数据重建的最佳心脏相位。根据本发明的示范性实施例，通过在延迟的获取协议的基础上完全自动地选择最佳心脏相位提供了一种用于导出最佳重建窗口的自动方法，在所述协议中，必须按照静态投影几何结构获取至少一个心脏相位。

Description

心脏相位确定
本发明涉及医学成像领域。具体而言，本发明涉及一种用于检查感兴趣对象的检查设备、一种检查感兴趣对象的方法、一种图像处理装置、一种计算机可读介质和一种程序单元。
旋转血管造影(RA)是一个越来越受到关注的领域。为了实现对三维心脏血管造影(3D-RCA)的常规临床使用，要求显著减少所需的用户交互或者要求完全自动生成三维冠状动脉造影照片。最近，已经出现了针对用于减少所需的用户交互的新的方法的初步工作。这些方法是以建模[1]或者全自动非运动补偿重建技术[2，3]或运动补偿重建技术[4，5]为基础的。
基于重建的方法似乎格外具有吸引力，因为这些方法可以完全自动地提供对下层解剖结构的准确表示。一个主要问题，尤其是在采用非运动补偿重建技术时遇到的主要问题是选择适当的心脏相位，从而使由冠状动脉的运动所引起的模糊微乎其微。
最近的公开文献[6，7]利用在时间上相继的体积之间的相关性作为帧间运动的指标，并建议根据帧间相关性的最大值选择最佳的心脏相位。尽管，已经证明这一方法提供了有关解剖结构的当前运动的有价值信息，但是其可能无法直接提供对某一帧的图像质量的测度。
在[8]中，提供了一种基于对相应的三维重建图像的直方图分析导出单个时间上的快照的图像质量分析的方法。但是，这一方法要求生成所有心脏相位下的多幅三维重建图像，从而通过图像质量指标的比较确定最佳的心脏相位，所述多幅图像的生成是一项耗时的工作。
人们希望在冠状动脉造影中使对用于图像重建的心脏相位的确定得到改进。
根据本发明的示范性实施例，提供了一种用于检查包括填充了造影剂的冠状动脉的感兴趣对象的检查设备，所述检查设备包括适于在延迟的获取轨迹的基础上获取感兴趣对象的一系列投影图像的获取单元，所述延迟的获取轨迹包括跟随有旋转获取阶段的静态获取阶段，其中，所述静态获取阶段指的是这样的阶段，期间所述获取单元处于静止的角位置下，并且所述旋转获取阶段指的是这样的阶段，期间所述获取单元发生旋转，以及适于在所述一系列投影图像的基础上确定第一减影图像和第二减影图像的确定单元，其中，所述确定单元还适于在所述第一减影图像和所述第二减影图像的基础上确定用于数据重建的心脏相位，其中，在一个心动周期内按照静态投影几何结构在所述静态获取阶段获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
因此，所述检查设备可以适于在对填充了CA(造影剂)的冠状动脉的二维投影图像的分析的基础上执行最佳心脏相位的提取，其中，所述二维图像是在(至少)一个心动周期内在静止的角位置上获取的。这样可以实现对最佳心脏相位的快速精细选择，从而使由心脏运动所引起的模糊微乎其微。
根据本发明的另一示范性实施例，所述一系列投影图像包括第一、第二、第三和第四投影图像，其中，所述第一减影图像对应于从第一投影图像减去第二投影图像所得到的图像，并且其中，所述第二减影图像对应于从第三投影图像减去第四投影图像所得到的图像。
因而，简单地通过使心脏的不同图像彼此相减(所述图像是在静态获取模式下获取的)，可以确定出最佳的心脏相位。
根据本发明的另一示范性实施例，所述第一和第二图像是相继(subsequent)图像，所述第三和第四图像是相继图像。
因此，仅使直接在彼此之后获取的图像彼此相减。
根据本发明的另一示范性实施例，其中，所述确定单元还适于确定对应于所述第一减影图像的第一方差(variance)值和对应于所述第二减影图像的第二方差值，并且其中，所述用于数据重建的心脏相位对应于具有较低的方差值的减影图像。
因而，根据本发明的这一示范性实施例，计算每一减影投影图像的方差，并将计算得到的方差彼此进行比较。之后，通过最低的方差确定用于图像重建的心脏相位。
根据本发明的另一示范性实施例，完全自动地确定用于数据或图像重建的心脏相位。
因此，不需要通过用户交互确定最适合的时间窗口。
根据本发明的另一示范性实施例，所述检查设备还包括心电图仪，其中，数据重建是指心电图门控三维冠状动脉数据重建。
根据本发明的另一示范性实施例，在一个心动周期内，按照静态投影几何结构获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
因而，在单个心动周期内获取确定最佳心脏相位所需的所有图像，随后执行标准的旋转数据获取方案。或者，可以在按照静态投影几何结构获取参考图像之前执行标准的旋转数据获取方案。
因此，可以通过提供用于导出最佳重建相位的自动方案避免通过目测在不同的心脏相位上重建的多个高分辨率数据集而实施耗时的交互式最佳心脏相位选择。
根据本发明的另一示范性实施例，所述检查设备适于作为三维计算机断层摄影设备和三维旋转x射线设备之一。
应当注意，在这一背景下，本发明不限于计算机断层摄影或旋转血管造影，在必须确定适当的重建相位从而对周期性运动的对象成像时，总是可以应用本发明。
根据本发明的另一示范性实施例，将所述检查设备配置成由材料测试设备、医学应用设备和微型CT系统构成的组中的一个。
根据本发明的另一示范性实施例，所述确定单元还适于在一系列减影图像的基础上确定用于数据重建的门控窗口的宽度。
本发明的应用领域可以是医学成像，具体而言，可以是心脏成像。
此外，根据本发明的另一示范性实施例，提供了一种采用检查设备对包括填充了造影剂的冠状动脉的感兴趣对象进行检查的方法，所述方法包括如下步骤：在延迟的获取轨迹的基础上获取感兴趣对象的一系列投影图像，所述延迟的获取轨迹包括跟随有旋转获取阶段的静态获取阶段，其中，所述静态获取阶段指的是这样的阶段，期间所述获取单元处于静止的角位置下，并且所述旋转获取阶段指的是这样的阶段，期间所述获取单元发生旋转；在所述一系列投影图像的基础上确定第一减影图像和第二减影图像；以及在所述第一减影图像和第二减影图像的基础上确定用于数据重建的心脏相位，其中，在一个心动周期内按照静态投影几何结构在所述静态获取阶段获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
此外，根据本发明的另一示范性实施例，所述方法包括如下步骤：确定所述第一减影图像的第一方差和所述第二减影图像的第二方差，并比较所述第一方差和所述第二方差，其中，所述用于数据重建的心脏相位对应于与所述第一方差和所述第二方差中的较低方差相对应的投影图像。
换言之，在静态投影图像获取之后，创建所述静态投影的相继投影的减影图像，并计算每一减影投影图像的方差。之后，确定与属于最低方差值的投影相对应的心脏相位。
根据本发明的另一示范性实施例，提供了一种图像处理装置，其包括存储器，用于存储在延迟的获取轨迹的基础上所获取的感兴趣对象的一系列投影图像，其中，所述静态获取阶段指的是这样的阶段，期间所述获取单元处于静止的角位置下，并且所述旋转获取阶段指的是这样的阶段，期间所述获取单元发生旋转，以及适于在所述一系列投影图像的基础上确定第一减影图像和第二减影图像的确定单元，其中，所述确定单元还适于在所述第一减影图像和第二减影图像的基础上确定用于数据重建的心脏相位，其中，在一个心动周期内按照静态投影几何结构在所述静态获取阶段获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
此外，根据本发明的另一示范性实施例，提供了一种计算机可读介质，其中存储有用于检查感兴趣对象的计算机程序，在由处理器执行时，所述计算机程序使所述处理器实施上述方法步骤。
此外，根据本发明的另一示范性实施例，提供了一种用于检查感兴趣对象的程序单元，在由处理器执行时，所述程序单元使所述处理器实施上述方法步骤。
可以通过计算机程序，即通过软件；或者通过采用一个或多个专用电子优化电路，即通过硬件；或者以混合形式，即借助于软件部件和硬件部件，来实现所述检查感兴趣对象的方法。
优选将根据本发明的实施例的程序单元加载到数据处理器的工作存储器内。因而，可以将所述数据处理器配备为执行本发明的方法的实施例。可以通过任何适当的程序设计语言，例如，C++来编写所述计算机程序，并且可以将所述计算机程序存储到诸如CD-ROM的计算机可读介质上。而且，可以从诸如万维网的网络获得所述计算机程序，通过所述网络，可以将所述计算机程序下载到图像处理单元或处理器或者任何适当的计算机内。
可以看到，本发明的示范性实施例的要旨在于，提供了一种检查设备，所述检查设备适于在延迟的获取协议的基础上通过完全自动的方式选择最佳的心脏相位，在所述协议中，按照静态投影几何结构获取至少一个心脏相位。这样可以快速精细地选择最佳心脏相位，从而使由心脏运动所引起的模糊降至最低。
本发明的这些及其他方面将从下文中描述的各实施例中变得显而易见，并参考其进行阐述。
在下文中，将参照下列附图说明本发明的示范性实施例。
图1示出了根据本发明的示范性实施例的旋转X射线扫描器的示意性表示；
图2示出了延迟的获取轨迹的示意性表示；
图3示出了相继投影图像的三个减影图像；
图4示出了针对三个不同数据集的减影投影图像的方差值的三个示范性绘图的示意性表示；
图5示出了门控三维重建体积图像的最大强度投影图像的结果；
图6示出了根据本发明的示范性方法的流程图；
图7示出了用于执行根据本发明的方法的示范性实施例的根据本发明的图像处理装置的示范性实施例。
附图中的图示是示意性的。在不同的附图中，类似或等同的元件具有相同的附图标记。
图1示出了根据本发明的示范性实施例的旋转X射线扫描器的示意性表示。将X射线源100和具有大的灵敏面积的扁平探测器101安装到C臂102的两端。C臂102由弯轨，即“套筒”103支撑。所述C臂能够在套筒103内滑动，从而执行围绕所述C臂的轴的“滚转(roll)运动”。将套筒103通过可旋转的接头附接到L臂104上，使所述套筒能够执行围绕这一接头的轴的“螺旋桨式运动”。将L臂104通过另一可旋转接头附接到天花板上，使所述L臂能够执行围绕这一接头的轴的旋转。所述的各种旋转运动是通过伺服电动机实现的。所述的三种旋转运动的轴与锥束轴总是在单个固定点，即旋转X射线扫描器的“等中心”105上相交。在所述等中心周围存在一定的体积，由沿着源轨迹的所有锥束投影出所述体积。这一投影体积的形状和尺寸取决于探测器的形状和尺寸以及源轨迹。在图1中，球110表示符合于所述投影体积的最大等中心球。将所要成像的对象(例如，患者或行李物品)放到台111上，从而使所述对象的感兴趣体积符合所述投影体积。因此，所述投影体积限制了感兴趣体积的尺寸。
通过控制单元112控制所述各种旋转运动。每一个由C臂角度、套筒角度和L臂角度构成的三元组都限定了X射线源的位置。通过随时间改变这些角度，能够使所述源沿规定的源轨迹运动。处于所述C臂的另一端的探测器执行对应的运动。所述源轨迹将被限定在等中心球体的表面上。
在下文中，描述了根据本发明的示范性实施例的用于ECG门控三维冠状动脉重建(3D-CR)的最佳心脏相位的全自动提取，该提取是以在(至少)一个心动周期内，在静止角位置下获取的填充了CA的冠状动脉造影照片的二维投影图像的分析为基础的。
图6示出了根据本发明的用于检查感兴趣对象的(当然还用于确定最佳的心脏相位的)示范性方法的流程图。所述方法是以延迟的获取协议为基础的，在该协议中，需要在静态视图中获取至少一个心脏相位。可以通过四个主要步骤来说明用于提取最佳心脏相位的方法：
在步骤1中，在静态设置下获取一系列投影图像。所述图像是以延迟的获取轨迹为基础的填充了CA的冠状动脉的投影图像。
之后，在步骤2中，创建所述静态投影的相继投影的减影图像。在步骤3中，计算每一减影投影图像的方差，在步骤4中，确定与属于最低的方差值的投影相对应的心脏相位。
图2示出了延迟的获取轨迹的示意性表示。大约在开始数据获取之后的1到1.5秒使C臂X射线系统开始旋转(参见附图标记204)。但是，应当注意，如果在更快的获取模式下执行数据获取，则可以更早开始旋转。在静态获取(参见附图标记201、202、203)期间，获取用于确定最佳心脏相位的一系列投影图像。在下述标准获取过程(参见附图标记204、205、206、207)期间，所述C臂系统发生旋转。
下面的一行用于描绘冠状动脉的图像对应于在上面一行图像中描绘的C臂系统的位置。
图3示出了相继投影图像的三个减影图像301、302、303。在旋转运行期间获取的两个相继投影的减影血管造影照片的方差与心脏运动相关。因此，大的方差值对应于在获取之间存在大的运动的相继投影图像(参见图301、302)。另一方面，小的方差值对应于小的运动(参见图像303)。
图4示出了针对体模(图4a)、动物(图4b)和人(图4c)数据集的减影投影图像的方差值的示范性绘图的图示。最佳的心脏相位对应于最小的方差值。
水平轴401、405和409对应于心脏相位，就图4a和图4b而言，其在0到90°的范围内变化，就图4c而言，其在0到大约65°的范围内变化。垂直轴402、406和410对应于方差(任意单位)。
绘图403示出了与体模404相关的方差，绘图407示出了与动物408相关的方差，绘图411示出了与人的心脏412相关的方差。
图5示出了对应于最大方差值(501，502)和最小方差值(503，504)的门控三维重建体积图像的最大强度投影图像的结果。在绘图505中描绘出对应的方差值。
本发明的一种示范性实现可以是利用了旋转获取的血管介入。具体而言，可以在门控三维冠状动脉重建和运动补偿冠状动脉重建领域实现本发明。
图7示出了用于执行根据本发明的方法的示范性实施例的根据本发明的数据处理装置700的示范性实施例。
图7所描绘的数据处理装置700包括连接至存储器702的中央处理单元(CPU)或图像处理器701，所述存储器702用于存储描绘诸如患者或行李物品的感兴趣对象的图像。可以将数据处理器701连接至多个输入/输出网络或诊断装置，例如，X射线C臂设备。还可以将数据处理器701连接至诸如计算机监视器的显示装置703，所述显示装置用于显示在数据处理器701中计算或调整的信息或图像。操作员或用户可以通过键盘704和/或在图7中未示出的其他输入/输出装置与数据处理器701交互。
此外，通过总线系统705，还可能将图像处理和控制处理器701连接至(例如)监视感兴趣对象的运动的运动监视器。例如，在对患者的肺部进行成像的情况下，运动传感器可以是呼气传感器。在对心脏进行成像的情况下，运动传感器可以是心电图仪。
可以将本发明的示范性实施例作为给CT扫描器控制台、成像工作站或PACS工作站的软件选项而出售。
应当注意，术语“包括”不排除其他元件或步骤，“一”或“一个”不排除多个。此外，还可以将与不同的实施例相关联地进行描述的元件组合起来。
还应当注意，权利要求中的附图标记不应理解为对权利要求范围的限定。
文献列表
[1]Movassaghi B、Rasche V、Grass M、Viergever M、Niessen W，“Aquantitative analysis of 3D coronary modeling from two or more projectionimages”，IEEE Trans.Med.Imag.，第12卷第23期，第1517-1531页，2004.
[2]Rasche V、Buecker A、Grass M、Koppe R、Op de Beek J、BertramsR、Suurmond R、Kuehl H、Guenther RW.Ecg-gated 3Drotational coronaryangiography(3DRCA)，in Computer Assisted Radiology and Surgery 2002Proceedings(CAR’02).2002，第827-831页，Springer.
[3]Rasche V、Buecker A、Grass M、Suurmond R、KoPPe R.Kuehl H.ECG-gated 3D Rotational Coronary Angiography，RSNA，83rd ScientificSession，2003
[4]Blondel C、Malandain G、Vaillant R等人，4-D Tomographicrepresentation of coronary arteries from one rotational X-ray sequence，Proc.MICCAI 2003，LNCS 2878，416-423(2003)
[5]Movassaghi B、Rasche V、Viergever M.A、Niessen W.J，“3Dcoronary reconstruction from calibrated motion-compensated 2D projectionsbased on semi-automated feature point detection”，in Proc.SPIE MedicalImaging：Image Processing，第5370卷，San Diego，CA，第1943-1950页，2004.
[6]Movassaghi B、Istel T、Rasche V，“Automatic gating windowpositioning for 3D rotational coronary angiography(3DRCA)”，in Proc.SPIEMedical Imaging：Image Processing，第5370卷，San Diego，CA，第1932-1942页，2004.
[7]Manzke R、Nielsen T.和Th、Hawkes D、Grass M，″Automaticphase determination for retrospectively gated cardiac CT″，Medical Physics，第31(12)卷，第3345-3362页，2004.
[8]Volker Rasche、Babak Movassaghi、Michael Grass、Dirk Schaefer、Arno Buecker，
Automatic Selection of the Optimal Cardiac Phase for GatedThree-Dimensional Coronary X-Ray Angiography，Academic Radiology，第13卷，第630-640页，2006.

Claims (12)

1.一种用于检查包括填充了造影剂的冠状动脉的感兴趣对象的检查设备，所述检查设备包括：
获取单元(101)，其适于在延迟的获取轨迹的基础上获取所述感兴趣对象的一系列投影图像，所述延迟的获取轨迹包括跟随有旋转获取阶段的静态获取阶段，其中，所述静态获取阶段指的是这样的阶段，期间所述获取单元处于静止的角位置下，并且所述旋转获取阶段指的是这样的阶段，期间所述获取单元发生旋转；
确定单元(112)，其适于在所述一系列投影图像的基础上确定第一减影图像和第二减影图像；
其中，所述确定单元(112)还适于在所述第一减影图像和所述第二减影图像的基础上确定用于数据重建的心脏相位，
其中，在一个心动周期内按照静态投影几何结构在所述静态获取阶段获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
2.根据权利要求1所述的检查设备，
其中，所述一系列投影图像包括第一、第二、第三和第四投影图像；
其中，所述第一减影图像对应于从所述第一投影图像减去所述第二投影图像所得到的图像；
其中，所述第二减影图像对应于从所述第三投影图像减去所述第四投影图像所得到的图像。
3.根据权利要求2所述的检查设备，
其中，所述第一和所述第二投影图像是相继投影图像，并且其中，所述第三和所述第四投影图像是相继投影图像。
4.根据权利要求1所述的检查设备，
其中，所述确定单元(112)还适于确定对应于所述第一减影图像的第一方差值和对应于所述第二减影图像的第二方差值；
其中，所述用于数据重建的心脏相位对应于具有较低方差值的所述减影图像。
5.根据权利要求1所述的检查设备，
其中，完全自动地确定所述用于数据重建的心脏相位。
6.根据权利要求1所述的检查设备，
还包括心电图仪；
其中，所述数据重建是ECG门控3D冠状动脉数据重建。
7.根据权利要求1所述的检查设备，其适于作为三维旋转X射线设备。
8.根据权利要求1所述的检查设备，其被配置为由材料测试设备和医学应用设备构成的组中的一个。
9.根据权利要求1所述的检查设备，
其中，所述确定单元(112)还适于在一系列减影图像的基础上确定用于数据重建的门控窗口的宽度。
10.一种采用检查设备对包括填充了造影剂的冠状动脉的感兴趣对象进行检查的方法，所述方法包括如下步骤：
在延迟的获取轨迹的基础上获取所述感兴趣对象的一系列投影图像，所述延迟的获取轨迹包括跟随有旋转获取阶段的静态获取阶段，其中，所述静态获取阶段指的是这样的阶段，期间所述获取单元处于静止的角位置下，并且所述旋转获取阶段指的是这样的阶段，期间所述获取单元发生旋转；
在所述一系列投影图像的基础上确定第一减影图像和第二减影图像；以及
在所述第一减影图像和所述第二减影图像的基础上确定用于数据重建的心脏相位，
其中，在一个心动周期内按照静态投影几何结构在所述静态获取阶段获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
11.根据权利要求10所述的方法，还包括如下步骤：
确定所述第一减影图像的第一方差和所述第二减影图像的第二方差；以及
比较所述第一方差和所述第二方差；
其中，所述用于数据重建的心脏相位对应于与所述第一方差和所述第二方差中的较低方差相对应的所述投影图像。
12.一种图像处理装置，所述图像处理装置包括：
存储器，其用于存储在延迟的获取轨迹的基础上所获取的包括填充了造影剂的冠状动脉的感兴趣对象的一系列投影图像，所述延迟的获取轨迹包括跟随有旋转获取阶段的静态获取阶段，其中，所述静态获取阶段指的是这样的阶段，期间所述获取单元处于静止的角位置下，并且所述旋转获取阶段指的是这样的阶段，期间所述获取单元发生旋转；
确定单元(112)，其适于在所述一系列投影图像的基础上确定第一减影图像和第二减影图像；
其中，所述确定单元(112)还适于在所述第一减影图像和所述第二减影图像的基础上确定用于数据重建的心脏相位，
其中，在一个心动周期内按照静态投影几何结构在所述静态获取阶段获取用于确定所述用于数据重建的心脏相位的所述一系列投影图像。
CN2007800404234A 2006-11-03 2007-10-24 心脏相位确定 Expired - Fee Related CN101528133B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP06123445.6 2006-11-03
EP06123445 2006-11-03
PCT/IB2007/054316 WO2008053401A1 (en) 2006-11-03 2007-10-24 Cardiac phase determination

Publications (2)

Publication Number Publication Date
CN101528133A CN101528133A (zh) 2009-09-09
CN101528133B true CN101528133B (zh) 2012-09-05

Family

ID=39199370

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007800404234A Expired - Fee Related CN101528133B (zh) 2006-11-03 2007-10-24 心脏相位确定

Country Status (4)

Country Link
US (1) US8175356B2 (zh)
EP (1) EP2091435A1 (zh)
CN (1) CN101528133B (zh)
WO (1) WO2008053401A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102008025524A1 (de) * 2008-05-28 2009-07-30 Siemens Aktiengesellschaft Verfahren und Tomographiegerät zur Ermittlung einer Abtastphase in Bezug zu einem EKG-Signal eines Patienten
CN102232293B (zh) 2008-11-24 2016-11-09 皇家飞利浦电子股份有限公司 传输三维视频信息的方法和重放三维视频信息的重放设备
WO2010108146A2 (en) 2009-03-20 2010-09-23 Orthoscan Incorporated Moveable imaging apparatus
US9412044B2 (en) * 2009-06-09 2016-08-09 Siemens Aktiengesellschaft Method of compensation of respiratory motion in cardiac imaging
US9025840B2 (en) 2010-03-12 2015-05-05 Koninklijke Philips N.V. Motion visualisation in angiographic images
US9125611B2 (en) 2010-12-13 2015-09-08 Orthoscan, Inc. Mobile fluoroscopic imaging system
EP2897522B1 (en) * 2012-09-21 2020-03-25 Cardioinsight Technologies, Inc. Physiological mapping for arrhythmia
JP2015130070A (ja) * 2014-01-07 2015-07-16 富士通株式会社 検出プログラム、検出方法および検出装置
CN106102581B (zh) * 2014-01-29 2019-12-31 皇家飞利浦有限公司 图像数据中移动的结构的分割
US9629587B2 (en) 2014-07-10 2017-04-25 General Electric Company Systems and methods for coronary imaging
CN104545973B (zh) * 2014-11-26 2017-08-25 沈阳东软医疗系统有限公司 心脏的ct扫描方法及装置
US9642589B2 (en) 2015-06-30 2017-05-09 General Electric Company Systems and methods for guided selection of acquisition parameters for medical imaging
US10349909B2 (en) 2015-06-30 2019-07-16 General Electric Company Systems and methods for flow rate compensated acquisition parameters for medical imaging
US11024062B2 (en) 2018-06-11 2021-06-01 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for evaluating image quality

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0193712A1 (de) * 1985-02-04 1986-09-10 Siemens Aktiengesellschaft Röntgendiagnostikeinrichtung für Subtraktionsangiographie

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6480560B2 (en) * 2001-03-16 2002-11-12 Ge Medical Systems Global Technology Company, Llc Methods and apparatus for motion gating using CT projection data
DE10129631A1 (de) * 2001-06-20 2003-01-02 Philips Corp Intellectual Pty Verfahren zur Rekonstruktion eines hoch aufgelösten 3D-Bildes
US6754522B2 (en) 2001-09-05 2004-06-22 Medimag C.V.I., Inc. Imaging methods and apparatus particularly useful for two and three-dimensional angiography
CA2533538A1 (en) 2003-07-21 2005-01-27 Paieon Inc. Method and system for identifying an optimal image within a series of images that depict a moving organ
US20050033123A1 (en) * 2003-07-25 2005-02-10 Siemens Medical Solutions Usa, Inc. Region of interest methods and systems for ultrasound imaging
JP4795658B2 (ja) * 2004-03-23 2011-10-19 富士フイルム株式会社 差分画像取得方法、差分画像取得装置、及び、そのプログラム
DE102004048209B3 (de) * 2004-09-30 2005-09-01 Siemens Ag Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie
JP4777007B2 (ja) * 2005-08-03 2011-09-21 東芝メディカルシステムズ株式会社 Ｘ線コンピュータ断層撮影装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0193712A1 (de) * 1985-02-04 1986-09-10 Siemens Aktiengesellschaft Röntgendiagnostikeinrichtung für Subtraktionsangiographie

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MANZKE R ET AL.Automatic phase determination for retrospectively gated cardiac CT.《MEDICAL PHYSICS》.2004,第31卷(第12期),3345-3362. *

Also Published As

Publication number Publication date
WO2008053401A1 (en) 2008-05-08
US8175356B2 (en) 2012-05-08
US20100074485A1 (en) 2010-03-25
CN101528133A (zh) 2009-09-09
EP2091435A1 (en) 2009-08-26

Similar Documents

Publication Publication Date Title
CN101528133B (zh) 心脏相位确定
Weustink et al. Optimal electrocardiographic pulsing windows and heart rate: effect on image quality and radiation exposure at dual-source coronary CT angiography
US10918305B2 (en) Method and system for 4D radiological intervention guidance (4D-cath)
Berger et al. Marker‐free motion correction in weight‐bearing cone‐beam CT of the knee joint
EP1595228B1 (en) Method for the 3d modeling of a tubular structure
Isola et al. Fully automatic nonrigid registration‐based local motion estimation for motion‐corrected iterative cardiac CT reconstruction
US20100014740A1 (en) Optimal rotational trajectory determination for ra based on pre-determined optimal view map
CN101505662B (zh) 用于旋转血管造影术的自动等中心
EP2224851B1 (en) Correction for un-voluntary respiratory motion in cardiac ct
US20210030374A1 (en) Image creation device
US20170196529A1 (en) System and method for motion artifacts reduction
Funama et al. Coronary artery stent evaluation by combining iterative reconstruction and high-resolution kernel at coronary CT angiography
WO2007138526A2 (en) Hierarchical motion estimation
Isola et al. Cardiac motion-corrected iterative cone-beam CT reconstruction using a semi-automatic minimum cost path-based coronary centerline extraction
Müller et al. Evaluation of interpolation methods for surface‐based motion compensated tomographic reconstruction for cardiac angiographic C‐arm data
Li et al. 3D coronary artery reconstruction by 2D motion compensation based on mutual information
Perrenot et al. Motion correction for coronary stent reconstruction from rotational X-ray projection sequences
Wink et al. Intra-procedural coronary intervention planning using hybrid 3-dimensional reconstruction techniques1
Rasche et al. Three-dimensional X-ray coronary angiography in the porcine model: A feasibility study
Nieman et al. Computed tomography
EP1800265A2 (en) Device, method, computer-readable medium and program element for processing image data of a movable object, and rotational x-ray angiography apparatus
Schoonenberg et al. Three-dimensional coronary visualization, Part 2: 3D reconstruction.
Keil et al. Dynamic cone beam reconstruction using a new level set formulation
Liu et al. Improved C-arm cardiac cone beam CT based on alternate reconstruction and segmentation
Yang et al. Simulation environment of X-ray rotational angiography using 3D+ t coronary tree model

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120905

Termination date: 20181024