CN103430211B - 用于灌注比较和量化的数据动态归一化 - Google Patents

用于灌注比较和量化的数据动态归一化 Download PDF

Info

Publication number
CN103430211B
CN103430211B CN201280014040.0A CN201280014040A CN103430211B CN 103430211 B CN103430211 B CN 103430211B CN 201280014040 A CN201280014040 A CN 201280014040A CN 103430211 B CN103430211 B CN 103430211B
Authority
CN
China
Prior art keywords
mrow
msub
contrast agent
blood vessel
image data
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
CN201280014040.0A
Other languages
English (en)
Other versions
CN103430211A (zh
Inventor
R·弗洛朗
C·迪富尔
V·奥夫雷
O·博纳富
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 CN103430211A publication Critical patent/CN103430211A/zh
Application granted granted Critical
Publication of CN103430211B publication Critical patent/CN103430211B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/48—Diagnostic techniques
    • A61B6/481—Diagnostic techniques involving the use of contrast agents
    • 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/48—Diagnostic techniques
    • A61B6/486—Diagnostic techniques involving generating temporal series of image 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/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/50—Clinical applications
    • A61B6/507—Clinical applications involving determination of haemodynamic parameters, e.g. perfusion 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/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
    • G06T5/94—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • 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/12—Devices for detecting or locating foreign bodies
    • 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/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/464—Displaying means of special interest involving a plurality of displays
    • 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/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • 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
    • 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
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion

Abstract

本发明涉及X射线成像技术以及图像后期处理。具体而言，本发明涉及通过对灌注图像数据进行绝对或相对归一化，对X射线成像装置采集的灌注图像数据进行后期处理，以就不同采集以及不同患者对图像数据进行优选比较。为了允许对灌注图像数据进行归一化，可能希望知道保留在冠状动脉中的已注入造影剂的量。接下来，可以基于冠状动脉之内已知的造影剂的量对图像参数进行调整或归一化，以对灌注图像数据进行归一化。为了获得注入的造影剂的精确量，可以估计流经血管的界定区域或截面的造影剂的注入体积。于是可以从这个位置处总体积流量的估计，推断出造影剂的所述注入体积。因此，提供了一种方法（10），用于对数据进行动态归一化，以进行灌注比较和量化，包括如下步骤：确定（20）血管中总体积流量或造影剂的量，以及基于所确定的总体积流量或造影剂的量对灌注数据进行归一化（34）。

Description

用于灌注比较和量化的数据动态归一化
技术领域
本发明涉及X射线成像技术以及图像后期处理。
具体而言，本发明涉及针对用于灌注比较和量化的数据动态归一化的方法和装置、计算机可读介质、程序单元以及处理设备。
具体而言，本发明涉及通过对灌注图像数据进行绝对或相对归一化，对X射线成像装置采集的灌注图像数据进行后期处理，以就不同采集以及不同患者对图像数据进行优选比较。
背景技术
当今，所谓的经皮冠状动脉介入或PC I是一种常见的手术程序。在这种介入期间，将导管插入患者进入部位的血管系统中。相应地沿大血管之内向需要处置的特定血管结构推进导管。
接下来经由导管注射造影剂，使用X射线成像装置采集血管造影图像序列，所述血管造影图像序列示出填充造影剂时的血管结构。
通常，例如，通过C型臂X射线成像装置获得二维图像数据。在其他实施例中，也可以采用三维医学成像模态获得三维血管结构。
将这样获得的血管造影图像数据提供给医生，接下来用于规划介入程序。不过，医生感兴趣的不仅仅是冠状动脉血管造影照片数据，而且还有关于心肌灌注的信息。
在这里，将灌注定义为流体通过淋巴系统或血管通向并进入器官或组织。于是，灌注扫描是一种可以观察、记录并量化灌注的过程。
例如，患有严重心脏病的患者通常会进行心脏导管插入术。在这种程序期间，可以确定冠状动脉狭窄的程度和/或动脉瘤尺寸。
不过，仅从冠状动脉血管造影照片提供的冠状动脉形状不能获得关于心肌灌注的信息。例如，这可能是因为一旦发生狭窄症，其余正常的冠状血管结构可以开始为局部缺血的肌肉提供血液，使得冠状动脉形状和心肌灌注之间几乎没有关系。
在采集冠状动脉血管造影照片图像数据期间，也可以获得心肌灌注图像数据。
为了改善心肌灌注图像数据的可见度，例如，可以从采集的图像减去，例如消除本身不属于冠状动脉血管造影照片或心肌灌注数据的背景特征或图像部分，从而获得仅可看到被注入的冠状动脉和心肌灌注的图像。
现有技术中这种过程被称为DSA过程。例如，为了检查患者的心脏，可以在注射造影剂之前采集与一个心搏周期对应的掩模图像，接下来将其从冠状动脉血管造影照片期间采集的图像数据中减掉。在与相同心搏相位相关的掩模图像和造影图像之间进行这种减法，由此从减法图像消除由于心搏导致的图像成分。
US 2004/0082846 A1描述了一种用于体积心脏计算机断层摄影扫描的方法和装置，用于估计选定区域中造影剂的密度。
US 2009/274358描述了一种方法和断层摄影扫描机，用于针对由造影剂产生的图像数据中的对比度对图像数据进行归一化。
发明内容
本发明的一个目的可以被视为提供用于分析、解释和量化灌注图像数据的优选模块。
可以由独立权利要求的主题实现这个目的。在从属权利要求中描述了本发明的优选实施例。
可以认为例如展开程序前后的心脏血管造影术中心肌灌注图像数据的比较或量化是复杂的，因为难以控制冠状动脉中注射的造影剂量。具体而言，从冠状动脉泄漏回流到主动脉中的造影剂不能忽略。
泵送到主动脉中的大量血液实际上创建压差，容易将经由注射导管通过对应的冠状动脉心门注射的造影剂的一部分吸回主动脉中。必须要考虑这种自然现象，以能够就不同的图像采集甚至不同患者对灌注图像数据，进行有意义的比较，甚至在心肌层之内进行量化。
尽管在造影剂注射期间的其他措施是已知的，例如具有受控量、稀释率和注射速率的ECG门控注射，但是仍然必须要考虑回流到主动脉中的造影剂影响。在比较灌注情形时，所考虑图像信息的灰度级标定的轻微修改都可能带来状态评估中的巨大差异。
在比较冠状动脉结构已经变化的灌注情形时，例如在展开程序前后，由于冠状血管结构变化造成灌注参数大大不同，更会出现这种情况。
可以进一步从下文描述的优选实施例的详细描述推断出本发明的方面、特征和优点，参考附图对其进行了解释。
可以利用类似的附图标记表示类似元件。
附图不是按比例绘制的，不过可以描绘定性的比例。
附图说明
图1示出了心肌层血管结构的示意图，
图2a-d示出了冠状动脉血管造影照片的示范性实施例，
图3示出了根据本发明用于对灌注数据进行动态归一化以进行比较和量化的方法的示范性实施例，以及
图4示出了用于对灌注数据进行动态归一化以比较和量化的示范性装置和X射线成像装置。
附图标记列表：
2 心肌层
4 主动脉
5 分叉点
6 右冠状动脉
7 主冠状动脉收集部位
8 左冠状动脉
8a 旋动脉
8b 冠状动脉左前降支
9 冠状动脉心门
10 用于动态归一化以进行灌注比较和量化的方法
20 步骤：造影剂体积积分
20a 步骤：体积流量估计
20b 步骤：造影剂测量
20c 步骤：校准
20d 步骤：确定3D血管几何结构/2D采集几何结构
21 步骤：对时间积分
22 步骤：受控采集
24 步骤：受控注射
26 步骤：输入数据
28 步骤：注射尖端跟踪
30 步骤：主要收集部位跟踪
32 步骤：心脏DSA过程
34 步骤：数据归一化
36 步骤：量化或比较
40 X射线系统
41 悬挂式C型臂
42 装置
44 X射线探测器
46 X射线源
48 对象
50 支撑体
51 X辐射
53 注射导管系统
55 可视化元件
具体实施方式
为了允许对灌注图像数据进行归一化，可能希望知道保留在冠状动脉中的已注入的造影剂的量。接下来，可以基于冠状动脉之内已知的造影剂的量来调整图像参数，以对灌注图像数据进行归一化。
为了获得注入的造影剂的精确量，可以估计流经血管的界定区域或截面的造影剂的注入体积。于是可以从这个位置处总体积流量的估计，推断出造影剂的所述注入体积。
例如，在血管造影程序期间，可以在实时采集的X射线图像上跟踪注射导管的尖端。接下来可以将注射导管的尖端引入冠状动脉中，以注射造影剂。接下来可以在相应冠状动脉的部分中布置注射导管的尖端，所述部分基本没有分叉点，例如，体现为具有定义几何结构的单一血管。后面将冠状动脉的所述部分称为主要收集部位，这样例示造影剂被注入的冠状动脉之内的段或区域。
常规上，不将导管尖端引入冠状动脉中远超过冠状动脉心门之外，使得造影剂的注射还发生在布置在相应冠状动脉心门附近的主要收集部位处。
在现在向冠状动脉中注入造影剂时，造影剂的主要部分经由冠状动脉流向心肌组织结构，而由于回流效应的缘故，少量，但是不能忽略的量的造影剂可以被视为泄漏回主动脉中。为了精确获得归一化的灌注图像数据，因而需要精确地知道或估计流经主要收集部位的、保留在相应冠状动脉中的、因此对采集心肌灌注图像数据有贡献的造影剂部分或量。
由于主动脉中大量的血流，泄漏回主动脉中的那部分造影剂对于局域化的灌注图像没有任何贡献，从而可以被完全忽略。
尽管可以在特定精确程度上知道造影剂的注射量，但上述主动脉回流需要造影剂量的更新确定。
为了获得流经主要收集部位的造影剂的精确量，本发明提出估计流经主要收集部位的造影剂积分体积。接下来可以使用这种估计，可能是估计到可以相乘的项，对整个图像轮次上的灌注图像数据进行归一化，允许灌注图像数据的适当量化或比较。
在已知血管结构的几何参数的情况下，可以获得绝对值，允许绝对的归一化。在不精确知道血管结构的情况下，至少可以获得造影剂的相对值，至少允许进行相对归一化。
不过，仅在比较不同患者时，才可能需要绝对值，因为对于相同的患者，任何相乘项在每次个体采集程序中都表现为相同，从而相对于不同的图像采集而消除其自身。
即便是实际流经心肌层的造影剂可变时，尤其是由于回流现象，所提出的方法也能够获得可比较的结果。
不过，在经由注射导管尖端向冠状血管中施予造影剂时，仍然可以考虑ECG触发和控制的注射。
可以从这个位置处总体积流量的估计，推断出流经主要收集部位的造影剂的注入体积估计值。
灌注比较的两个阶段之间的另一种可能偏差可能来自那些阶段期间发生的不同血流速度。例如，在冠状动脉修复(搭支架)前后，由于狭窄位置处修复的血管形状，血流可能不同。现在，由于通过浑浊化观察灌注，并且由于更高的血流量意味着由于更高的稀释导致的更少浑浊化，所以任何流动偏差都可能导致灌注评估偏差。如果在不同速度下注入造影剂，这同样是适用的。但在两种情况下，在主要收集部位中随时间测量血液体积使得能够对两种情形进行归一化。
图像采集程序期间的若干考虑可能对于后继归一化程序有益。
尽管基于流经主要收集部位的造影剂注入体积的估计进行归一化能够消除所述不同造影剂体积或造影剂的量所导致的图像数据中的差异，但可能有益的是还控制X射线成像系统的采集参数以在不同的图像采集程序之间，可能甚至在不同患者之间，进行接近的匹配。换言之，特定的采集设置可能对于灌注量化和/或比较是有用的。例如，可能有益的是，关闭对图像采集系统的连续自动剂量控制，至少为了获得适当的数字相减图像结果。
而且，在采集不同的灌注数据集期间，采用X射线成像系统的相同或可比拟的系统设置或图像采集参数，例如用于扩展前/后的比较，可能是有益的。
此外，为了为后续的造影剂体积估计提供充分的依据，可以精确控制向主要收集部位中注入造影剂的程序，以进行精细的灌注分析。这样的受控注射可以包括利用或根据特定的心搏周期来触发注射，即由ECG触发，并进一步控制注射量、稀释、或总血管体积与造影剂体积之比、以及注射速度。
为了获得灌注图像数据的有意义结果，可能需要前述掩模，使得仅剩余具有随后的灌注图像数据的冠状动脉血管造影照片图像数据。换言之，输入数据(可以是注射血管造影照片)应当包含注射阶段本身和可以用于掩模目的的注射前阶段。
在注射前阶段期间，可能已经采集ECG触发或依赖于ECG的图像，接下来可以采用其在所获得的灌注图像数据中实现不相关组织结构的掩模目的。两个阶段，即注射前阶段和注射阶段本身，应当与后续的心脏DSA过程兼容，这是本领域中公知的。
为了在主要收集部位中进行造影剂的优选施予，可以沿输入序列，例如通过实时X射线图像，跟踪导管尖端或注射尖端，直到正确定位。可以采用已知的图像处理分割和跟踪工具，辅助注射尖端跟踪过程。注射尖端的容易识别的特性，例如已知的一段有造影剂管，可以使这个任选步骤简单，这显著方便了主要收集部位的跟踪。
在这里，跟踪主要血管收集部位，即位于注射尖端输出处，恰在心门之后的血管段或区域。对于主要血管收集部位，可以仅需要短段。这个血管段不应该包含相当大的分叉点。对于冠状动脉而言，对于左冠状动脉树，收集部位可以是所谓的“左主干”。对于右冠状动脉树，主要血管收集部位可以例如是主动脉从其出现的血管桩(vessel stump)。
可以由注射尖端跟踪步骤的结果大大促进这个步骤，因为可以认为主要收集部位与注射导管对准。
神经学介入是可以应用本发明的另一种情形。在注射造影剂之后在大脑中评估灌注。同样，为了进行适当的灌注比较，在注射点处或在本地观测的区域上游任一点处随时间测量血液体积流量可能是必要的。甚至可以应用相同的原理对动脉瘤之内的造影剂活动测量值进行归一化，通常在流动转向器植入前后。尽管动脉瘤腔中的造影剂活动精确地讲不是灌注效果，但是从更宽的角度而言，可以应用相同原理。可以将灌注的概念扩展到远离馈送点的区域中造影剂浑浊化的研究中。
在图像采集程序的这一点，可能已经采集了注射前阶段的图像，具体而言是依赖于ECG的图像。现在，可以开始造影剂的注射。
可以在所谓的心脏DSA过程中处理注射造影剂之后采集的血管造影图像。在这个步骤中，实现心脏DSA。可以将所述心脏DSA过程视为相当于针对每个当前注射帧或图像从注射前阶段图像发现非注射帧，可以采用其进行相减，以便消除背景图像信息，仅留下造影剂贡献。在非心脏情况下，例如在神经学中，DSA常常简单地相当于仅在无造影剂图像(掩模)和当前注射造影剂的图像系列之间进行相减。可能需要对未注射图像进行几何学翘曲，以针对优选的减法和背景图像信号的消除进行可能的运动补偿。于是，除了相应冠状动脉中的造影剂贡献之外，运动诱发的伪影可能保留在DSA已处理图像中。
接下来，作为本发明的中心方面，进行对时间积分的造影剂体积估计或造影剂的量估计，实际上在几个步骤中实现。更一般地，通过随时间估计总体积流量来进行归一化。
首先，进行总体积流量估计。主要收集部位之内的总体积流量涉及总流体(既有血液又有造影剂)的体积流量。在DSA图像中的主要收集部位的界定精细剖面或截面上估计所述总流体流量。
对于总体积流量估计，本领域已知几种方法(例如，J.Waechter，J.Bredno，PMJRongen，J.Weese，D.C.Barratt，D.J.Hawkes：“Model-based blood flow quantificationfrom rotational angiography”；Medical Image Analysis，卷12，第5期，第586-602页(2008年10月)，在此通过引用将其公开并入)。
而且，可以通过在主要收集部位本地布置的流量计估计或确定总体积流量。
总体积流量估计可以是绝对估计，或可以将体积流量至少估计到相乘因数在后一种情况下，可以不必知道或估计主要收集部位的透视缩短和直径，即血管结构的几何特征。
这一总体积流量估计可以提供总绝对体积流量qT(t)或相对总体积流量q′T(t),它们根据方程1彼此相关。
方程1
接下来，作为下一步，在投影中看到的主要收集部位的同一精细截面上并在单位时间样本上测量造影剂流量。接下来将这个值称为时间t处的C(t)。为了获得C(t)，考虑注射的造影剂的体积流量。注射的造影剂仅仅指保留在相应血管中且未泄漏回主动脉中的造影剂部分，即流经主要收集部位的部分。
最优地，校准图像采集系统，以获得常量或相乘因子α，将观察的对比度，即观察的图像对比度与比吸收水平相关。常量α也可以取决于所用造影剂的性质及它的适当稀释。
因此，考虑α，观察到的与特定吸收水平相关的图像对比度允许接下来确定注入到主要收集部位中的造影剂的体积流量。于是，可以通过方程2确定观察到的图像对比度C(t)。
方程2
qI(t)与主要收集部位中注射造影剂的体积流量相关，而比例与造影剂的稀释相关，因此，直接与图像中观察到的吸收水平相关。
从系统校准和造影剂类型确定α。系统校准简单地表示定义的物理吸收可能与之相关的定义的灰度级。假设可以仅通过相乘因数来实现这个目的，但偏移也可能干预。可以利用已知的吸收体模或测试模式精确的执行这样的校准。
C(t)是观察到的对比度，或等价地，已消除背景的DSA图像中观察到的灰度级值，其根据注射物和总流量的比例随时间变化。例如，如果总流量突然相对于注射流量增长10倍，那么观察到的对比度将简单地减小10倍。于是，该比值是血液中造影剂的稀释倍数。
接下来，通过从开始注射时间t0到当前时间t进行造影剂体积积分，随时间确定心肌层中注射造影剂的量。随时间的估计值特别重要，因为在整个序列中需要它来进行归一化。
可以获得体积，从而获得该量，要么到相乘项，因此得到Q′I(t)，要么获得绝对值QI(t)。相应的确定采用的是方程3a或3b。
方程3a
方程3b
在获得造影剂体积积分信息，从而获得注射造影剂的量之后，执行数据归一化过程。在这里，可以采用量QI(t)或Q′I(t)以对应于时间t的DSA图像信息的当前图像数据D(t)进行归一化。
可以通过方程4A和4B实现两个数据集D1(t)和D2(t)之间的相对归一化，例如，用于比较搭支架前的情形与搭支架后的情形，从而获得可比的图像信息。
可比于
方程4a、b
在考虑同一患者时，用于比较的一致相对归一化可以是足够的。
为了考虑和比较不同患者或者为了进行实际量化，可能需要根据方程5进行绝对归一化。
方程5
在实践中，D(t)通常指峰值浑浊化(POP)，通常可以将其用于灌注评估。在每个像素x处，在所考虑的采样时间t上估计对比度。例如，在心脏应用中，通常针对给定的心脏时相定义每次心跳一个样本，而在神经应用中，可以定义采样速率，尤其是对应于实际帧速率。这定义了对比度图像序列C(x,t)，并根据方程6定义峰值浑浊化图像。
方程6
现在，可以考虑归一化方案包括归一化因数的并入，可以认为其独立于像素，因为仅在主要收集部位层次上定义，通过采用方程7进入这一最大计算。
方程7
峰值浑浊化指标可能不是评估灌注的唯一手段，但通常可以将相同种类的推理应用于灌注归一化。
可以通过采用通常的量化或比较手段，例如基于并排颜色的有色表示(blushrepresentation)，来进行后续的可视化、比较或量化。对于绝对归一化而言，可以计算并显示绝对量化数字/值，无论精确的注射条件如何。
通过结合对数据进行动态归一化以进行灌注比较和量化的方法、装置、计算机可读介质、程序单元、处理设备以及用于操作设备的方法描述的特征描述了本发明。不过，应该理解，不应认为特征仅仅涉及特定范畴，而是可互换的，甚至可以用于用于对数据进行动态归一化以进行灌注比较和量化的方法、装置、计算机可读介质、程序单元、处理设备以及用于操作设备的方法的每一相应种类中。
现在参考图1，示出了心肌层的血管结构示意图。
图1示出了心肌层2的示意图。示出了主动脉4，右冠状动脉6以及左冠状动脉8从其延伸出来。相应的冠状动脉6，8在从主动脉4进入相应冠状动脉6，8的进口处具有冠状动脉心门9。左冠状动脉8接下来分成旋动脉8a和冠状动脉左前降支8b。
可以将每个心门9和相应第一分叉点5之间的相应冠状动脉6、8的段落称为主冠状动脉收集部位7。从所述主冠状动脉收集部位7，仅将小的界定截面用于根据本发明的方法，尤其是总体积流量估计。
在血管造影程序期间，通常在任一时刻仅对冠状动脉6、8之一成像，因为在特定时间仅可以在冠状动脉6、8之一中注射造影剂。在DSA过程之后，只有相应冠状动脉6、8的血管结构保留在灌注图像中，因为DSA过程应确保从图像中消除不相关信息，例如心肌组织以及主动脉。
现在参考图2a到d，示出了冠状动脉灌注图像的示范性实施例。这些图对应于给定心脏阶段的峰值浑浊化(POP)可视化。
图2a到d例示了不采用标准程序时的困难。
对于图2a、b，两幅图都涉及基本相同的本情况，甚至可能是相同的图像数据，不过图2a和2b是以不同的比例或不同的参考峰浑浊化值示出的，因此导致相同解剖情形的显著不同描绘。实际上，被成像的左冠状动脉8未被绘示到与图2a、b相同的程度。于是，没有归一化，图2a和2b之间的有意义比较因而是不可行的。
对于图2c和2d，图2c示出了搭支架前的灌注情形，而图2d示出了相应的相关联的搭支架后的情形。
可以关于图2c和2d手动地调节参考值或比例，因此允许对解剖情形搭支架前后进行有意义的比较。
通过采用根据本发明的归一化方法，可能不再需要比例的手动调节。相反，可以获得直接可比较的图像信息或心脏灌注图像数据。示范性地，图2a到2d示出了左冠状动脉的情形。
现在参考图3，示出了根据本发明的用于动态归一化以进行心脏灌注比较和量化的方法的示范性实施例。
图3中示出了用于动态归一化以进行灌注比较和量化的核心方法10，嵌入于进一步的步骤中。
如前所述，造影剂的受控采集22以及受控注射24都允许采集有意义的输入数据或冠状动脉血管造影照片数据。
任选的注射尖端跟踪28以及主要收集部位跟踪30都允许后续的造影剂积分步骤20。具体而言，主要收集部位跟踪30便于确定和监测要估计并关于时间跟踪体积流量的位置。
后续心脏DSA过程32允许最终获得灌注图像信息。
根据本发明的方法10包括造影剂体积积分步骤20以及后续的数据归一化步骤34，得到有意义的图像信息，以进行灌注图像信息的比较或甚至量化。
造影剂体积积分步骤20包括在主要收集部位中确定或估计总体积流量的个体子步骤。此外，执行造影剂流量测量20b，从而提供所有所需的参数，用于在时间上进行后续积分21。
为估计相对归一化Q′I(t)，校准步骤20c是必要的。不过，为了确定绝对的归一化QI(t)，需要来自校准的α和绝对总体积流量qT(t)。为了获得绝对归一化QI(t)，还需要3D血管几何结构和采集几何结构，即关于3D血管以及该血管投影到2D平面上的方式的知识。没有这些值，仅可以估计q′T(t)，因此仅估计Q′I(t)。
可以从3D数据采集、介入前或附近，并从3D中的血管分割，获得血管几何结构。可以从该系统采集到采集几何结构。此外，为了以绝对的方式分析DSA中观察到的造影剂，可能需要将3D数据对准到考虑的2D投影。典型范例是在空间中倾斜并在血管造影照片图像中以一定程度的透视缩小投影的血管段。3D血管几何结构提供了血管的直径及其在空间中的倾斜角。接下来假设这个血管段是精确的馈送点(即，主要收集部位)。如果与血管造影照片正确对准，在血管造影照片中观察到的这一主要收集部位可以与3D血管段相关。知道了2D采集几何结构，可以推断出3D和2D之间的透视缩短比例。
现在，因为透视缩短小成更多吸收，因为更多造影剂被通过倾斜血管的X射线束通过，并因为相比于3D实际血管长度，2D投影中的血管段长度减小，所以应当引入若干校正以将2D中观察到的体积流量转换成3D中的实际体积流量。通常，对于这样的校正而言，仅仅透视缩小比例β就足够了，该比例是3D血管段和图像平面之间角度的余弦倒数。实际上，要将观察到的长度乘以β以获得实际3D长度，这意味着在2D中估计的流量速度也应乘以β以获得3D中对应的速度。对于观察到的血管截面，相同种类的推理通常是成立的。
在观察到的2D和3D血管截面之间还有近似β因子。不过，在有3D时，可以直接从3D中的估计速度直接确定体积流量，即2D中估计的速度的β倍，乘以实际的3D截面积。
换言之，为了获得绝对总血液体积流量或绝对归一化，需要校准步骤20c(对于α)和步骤20d(对于β),它们涉及到确定3D血管几何结构、确定2D采集几何结构以及确定它们之间的关系。于是，步骤20c的输出可以被视为是α,而步骤20d的输出可以被视为是
现在参考图4，示出了一种示范性的X射线成像装置，用于对图像数据进行动态归一化，以进行灌注比较和量化，尤其是心脏灌注成像。
图4示出了X射线成像系统40，包括用于对数据进行动态归一化以进行灌注比较和量化的装置42。
X射线成像系统40示范性地包括悬挂式C型臂41，其具有X射线源46和X射线探测器44。对象48位于支撑体50上，可以布置于X辐射51的路径中，X辐射51从X射线源发出、穿透对象48，接下来到达X射线探测器44，以采集对象48内部结构的X射线图像。
仅示意性地示出了注射导管系统53，其用于在PCI程序期间应用造影剂，用于获得冠状动脉血管造影照片并接下来获得心肌灌注图像数据。
装置42可以采用事先采集的血管造影照片图像数据进行归一化，或者可以采用从X射线系统40实时采集的X射线图像信息。
装置42适于执行如图3所示的方法，尤其是心脏DSA过程32的步骤，具有一些或全部相关联的步骤20a、b、c、21的造影剂体积积分20以及数据归一化34。
为了进行可视化、比较或量化，可以利用可视化元件55向用户显示图像信息。例如，可以利用可视化元件55进行搭支架前后图像数据的并排比较。
装置42也可以向用户提供从归一化的灌注图像信息导出的绝对量化数据。

Claims (18)

1.一种用于对图像数据进行动态归一化以进行灌注比较和量化的装置，包括处理元件，所述处理元件适于执行用于对图像数据进行动态归一化以进行灌注比较和量化的方法，
其中，所述处理元件适于确定流经血管的定义截面的造影剂的积分体积，其中，造影剂的所述积分体积是通过所述血管的所述定义截面中的总体积流量和观察到的所述血管的所述定义截面上的平均对比度量而确定的，并且
其中，所述处理元件适于基于造影剂的所述积分体积对所述图像数据进行归一化。
2.根据权利要求1所述的装置，还包括：
可视化元件；
其中，所述可视化元件适于向用户可视地提供第一图像采集的第一归一化灌注X射线图像数据和第二图像采集的第二归一化灌注X射线图像数据中的至少一个。
3.根据权利要求1或2所述的装置，还包括：
X射线源，以及
X射线探测器，
其中，对象能够布置于所述X射线源和所述X射线探测器之间，并且
其中，操作性耦合所述X射线源和所述X射线探测器，以采集所述对象的感兴趣区域的X射线图像数据，并且其中，要被归一化的所述图像数据是所述X射线图像数据。
4.根据权利要求3所述的装置，其中，所述X射线图像数据是灌注X射线图像数据。
5.一种用于对图像数据进行动态归一化以进行灌注比较和量化的方法，包括如下步骤：
确定流经血管的定义截面的造影剂的积分体积，其中，造影剂的所述积分体积是通过所述血管的所述定义截面中的总体积流量和观察到的所述血管的所述定义截面上的平均对比度量而确定的；以及
基于造影剂的所述积分体积对所述图像数据进行归一化。
6.根据权利要求5所述的方法，
其中，通过以下中的至少一项来确定所述总体积流量：由布置于所述血管的所述定义截面处的流量计做出的总体积流量估计，以及估计X射线图像信息中所述血管的所述定义截面处的总流体流量。
7.根据权利要求5所述的方法，
其中，所述方法还包括以下至少一个步骤：
a)通过如下方程确定造影剂的绝对积分体积QI(t)，
<mrow> <msub> <mi>Q</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mi>t</mi> </msubsup> <msub> <mi>q</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>&amp;alpha;</mi> </mfrac> <mo>&amp;CenterDot;</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mi>t</mi> </msubsup> <msub> <mi>q</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> </mrow>
以及
b)通过如下方程确定造影剂的相对积分体积Q＇I(t)，
<mrow> <msub> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mi>t</mi> </msubsup> <msub> <msup> <mi>q</mi> <mo>&amp;prime;</mo> </msup> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> </mrow>
其中，qI(t)对应于所述血管的所述定义截面处造影剂的绝对体积流量；
其中，qT(t)对应于所述血管的所述定义截面处的绝对总体积流量；
其中，q＇T(t)对应于所述血管的所述定义截面处造影剂的相对体积流量；
其中，C(t)对应于观察到的所述定义截面上的平均对比度；并且
其中，α是常数，t0是开始注射时间，而t是当前时间。
8.根据权利要求7所述的方法，
其中，所述方法还包括以下至少一个步骤：
通过如下方程确定绝对归一化图像数据Dn(t)，
<mrow> <msub> <mi>D</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>Q</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
以及
通过如下方程确定相对归一化图像数据D′n(t)，
<mrow> <msub> <msup> <mi>D</mi> <mo>&amp;prime;</mo> </msup> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
其中，D(t)对应于要被归一化的所述图像数据。
9.根据权利要求8所述的方法，其中，要被归一化的所述图像数据在X射线图像数据中。
10.根据权利要求5所述的方法，
其中，所述血管是在预先采集的X射线图像数据中的血管的表示，并且/或者
其中，所述方法是用于预先采集的X射线图像数据的图像后期处理的方法。
11.一种用于对图像数据进行动态归一化以进行灌注比较和量化的装置，包括：
用于确定流经血管的定义截面的造影剂的积分体积的模块，其中，造影剂的所述积分体积是通过所述血管的所述定义截面中的总体积流量和观察到的所述血管的所述定义截面上的平均对比度量而确定的；以及
用于基于造影剂的所述积分体积对所述图像数据进行归一化的模块。
12.根据权利要求11所述的装置，
其中，通过以下中的至少一项来确定所述总体积流量：由布置于所述血管的所述定义截面处的流量计做出的总体积流量估计，以及估计X射线图像信息中所述血管的所述定义截面处的总流体流量。
13.根据权利要求11所述的装置，
其中，所述装置还包括以下至少之一：
a)用于通过如下方程确定造影剂的绝对积分体积QI(t)的模块，
<mrow> <msub> <mi>Q</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mi>t</mi> </msubsup> <msub> <mi>q</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>&amp;alpha;</mi> </mfrac> <mo>&amp;CenterDot;</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mi>t</mi> </msubsup> <msub> <mi>q</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> </mrow>
以及
b)用于通过如下方程确定造影剂的相对积分体积Q′I(t)的模块，
<mrow> <msub> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mi>t</mi> </msubsup> <msub> <msup> <mi>q</mi> <mo>&amp;prime;</mo> </msup> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>C</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> </mrow>
其中，qI(t)对应于所述血管的所述定义截面处造影剂的绝对体积流量；
其中，qT(t)对应于所述血管的所述定义截面处的绝对总体积流量；
其中，q′T(t)对应于所述血管的所述定义截面处造影剂的相对体积流量；
其中，C(t)对应于观察到的所述定义截面上的平均对比度；并且
其中，α是常数，t0是开始注射时间，而t是当前时间。
14.根据权利要求13所述的装置，
其中，所述装置还包括以下至少之一：
用于通过如下方程确定绝对归一化图像数据Dn(t)的模块，
<mrow> <msub> <mi>D</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>Q</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
以及
用于通过如下方程确定相对归一化图像数据D＇n(t)的模块，
<mrow> <msub> <msup> <mi>D</mi> <mo>&amp;prime;</mo> </msup> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
其中，D(t)对应于要被归一化的所述图像数据。
15.根据权利要求14所述的装置，其中，要被归一化的所述图像数据在X射线图像数据中。
16.根据权利要求11所述的装置，
其中，所述血管是在预先采集的X射线图像数据中的血管的表示，并且/或者
其中，所述装置是用于预先采集的X射线图像数据的图像后期处理的装置。
17.一种在其中执行计算机程序的处理设备，其中，所述处理设备适于控制根据权利要求1至4中任一项所述的装置来执行根据权利要求5至10中任一项所述的方法，用于对灌注X射线图像数据进行归一化。
18.一种在其中存储了用于执行根据权利要求5至10中任一项所述的方法的计算机程序的计算机可读介质。
CN201280014040.0A 2011-03-18 2012-02-28 用于灌注比较和量化的数据动态归一化 Active CN103430211B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP11305304 2011-03-18
EP11305304.5 2011-03-18
PCT/IB2012/050921 WO2012127339A1 (en) 2011-03-18 2012-02-28 Dynamic normalization of data for perfusion comparison and quantification

Publications (2)

Publication Number Publication Date
CN103430211A CN103430211A (zh) 2013-12-04
CN103430211B true CN103430211B (zh) 2018-05-08

Family

ID=45833486

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280014040.0A Active CN103430211B (zh) 2011-03-18 2012-02-28 用于灌注比较和量化的数据动态归一化

Country Status (7)

Country Link
US (1) US10932739B2 (zh)
EP (1) EP2686826B1 (zh)
JP (1) JP6133794B2 (zh)
CN (1) CN103430211B (zh)
BR (1) BR112013023556A2 (zh)
RU (1) RU2013146539A (zh)
WO (1) WO2012127339A1 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102012205351A1 (de) * 2012-04-02 2013-10-02 Siemens Aktiengesellschaft Darstellung von Blutgefäßen und Gewebe im Herzen
US9247920B2 (en) 2014-02-27 2016-02-02 General Electric Company System and method for performing bi-plane tomographic acquisitions
JP2015217170A (ja) * 2014-05-19 2015-12-07 株式会社東芝 Ｘ線診断装置
GB2533566B (en) * 2014-12-19 2018-02-14 Gen Electric Method for contrast medium injection volume estimation based on image analysis
US9990712B2 (en) 2015-04-08 2018-06-05 Algotec Systems Ltd. Organ detection and segmentation
US10147171B2 (en) * 2016-09-21 2018-12-04 General Electric Company Systems and methods for generating subtracted images
EP3420906A1 (en) * 2017-06-29 2019-01-02 Koninklijke Philips N.V. Image contrast enhancement of an x-ray image
CN110739046A (zh) * 2018-07-18 2020-01-31 安影科技（北京）有限公司 一种多期多参数脑灌注成像平台
CN109669151A (zh) 2019-02-12 2019-04-23 泰山医学院 磁共振动脉自旋标记灌注成像的质量控制体模及评价方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0657208B2 (ja) 1985-12-27 1994-08-03 株式会社東芝 デイジタル・フルオログラフイ−装置
US7289841B2 (en) 2002-10-25 2007-10-30 Koninklijke Philips Electronics N.V. Method and apparatus for volumetric cardiac computed tomography imaging
DK1747477T3 (da) * 2004-05-04 2013-06-17 Stiftelsen Universitetsforskning Bergen Blindbestemmelse af arteriel indgangs- og vævsrestfunktioner i perfusions-mri
JP4795694B2 (ja) * 2004-06-11 2011-10-19 株式会社東芝 Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
CN101448457B (zh) * 2006-05-22 2011-03-30 皇家飞利浦电子股份有限公司 由投影成像得到的经运动补偿的冠状动脉血流
JP5449651B2 (ja) * 2007-03-09 2014-03-19 株式会社東芝 Ｘ線ｃｔ装置および心筋パーフュージョン情報生成システム
US8553832B2 (en) 2007-05-21 2013-10-08 Siemens Aktiengesellschaft Device for obtaining perfusion images
WO2009040742A2 (en) 2007-09-27 2009-04-02 Philips Intellectual Property & Standards Gmbh Methods for imaging the blood perfusion
RU2492812C2 (ru) * 2008-01-23 2013-09-20 Михалакис АВЕРКИУ Оценка лечения с использованием ультразвуковых контрастных веществ
JP5148315B2 (ja) * 2008-02-25 2013-02-20 株式会社東芝 医用画像処理装置、及び医用画像処理プログラム
DE102008021835A1 (de) * 2008-04-30 2009-11-05 Siemens Aktiengesellschaft Verfahren und Tomographiegerät zur Normierung von Bilddaten hinsichtlich eines durch ein Kontrastmittel in den Bilddaten hervorgerufenen Kontrastes
CN102123665B (zh) 2008-08-13 2013-09-11 皇家飞利浦电子股份有限公司 冠状血管和心肌灌注信息的动态可视化
JP5491929B2 (ja) * 2009-04-10 2014-05-14 株式会社東芝 Ｘ線診断装置及びその方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
放射性肺损伤64层螺旋CT灌注成像的临床初步研究;胡晓云;《万方学位论文数据库》;20081106;正文第2,13-19页 *

Also Published As

Publication number Publication date
JP2014511713A (ja) 2014-05-19
WO2012127339A1 (en) 2012-09-27
EP2686826A1 (en) 2014-01-22
JP6133794B2 (ja) 2017-05-24
RU2013146539A (ru) 2015-04-27
EP2686826B1 (en) 2018-08-22
US10932739B2 (en) 2021-03-02
US20140005538A1 (en) 2014-01-02
CN103430211A (zh) 2013-12-04
BR112013023556A2 (pt) 2016-12-06

Similar Documents

Publication Publication Date Title
CN103430211B (zh) 用于灌注比较和量化的数据动态归一化
US10803994B2 (en) Vascular flow assessment
US11694339B2 (en) Method and apparatus for determining blood velocity in X-ray angiography images
US9949650B2 (en) Fractional flow reserve estimation
CN105120738B (zh) 狭窄治疗规划
WO2006061814A1 (en) Method and apparatus for finding the tubular organs blood velocity and flow and related parameters
US20210259559A1 (en) Method and device for obtaining vascular pressure difference
JP2021528137A (ja) 定量的血行動態フロー分析のための方法および装置
US20210228171A1 (en) Dynamic angiographic imaging
JP5820549B2 (ja) 灌流イメージング
CN109688908A (zh) 用于确定血流储备分数的装置和方法
Malkasian et al. Quantification of vessel-specific coronary perfusion territories using minimum-cost path assignment and computed tomography angiography: validation in a swine model
US10201321B2 (en) Low-dose CT perfusion technique
US8488861B2 (en) System and method of automatic estimation of arterial input function for evaluation of blood flow
CN104825177B (zh) 修正脑血容积的方法及装置、血管照影设备
US20230102646A1 (en) Provision of a comparison dataset
WO2023023248A1 (en) Systems and methods of identifying vessel attributes using extravascular images
KR101944854B1 (ko) 선택적 컴퓨터 단층촬영을 이용한 혈류 모델링 방법 및 그 장치
Close et al. Vascular parameters from angiographic images

Legal Events

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