CN105848577B - 用于心肌灌注病理学表征的***和方法 - Google Patents

用于心肌灌注病理学表征的***和方法 Download PDF

Info

Publication number
CN105848577B
CN105848577B CN201480054511.XA CN201480054511A CN105848577B CN 105848577 B CN105848577 B CN 105848577B CN 201480054511 A CN201480054511 A CN 201480054511A CN 105848577 B CN105848577 B CN 105848577B
Authority
CN
China
Prior art keywords
perfusion
heart
intensity
medical images
computer
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
CN201480054511.XA
Other languages
English (en)
Other versions
CN105848577A (zh
Inventor
M·布雷沃
A·基里比里
E·纳格尔
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.)
Kings College London
Koninklijke Philips NV
Original Assignee
Kings College London
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 Kings College London, Koninklijke Philips NV filed Critical Kings College London
Publication of CN105848577A publication Critical patent/CN105848577A/zh
Application granted granted Critical
Publication of CN105848577B publication Critical patent/CN105848577B/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
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0044—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis 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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56366—Perfusion imaging
    • 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
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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/10072—Tomographic images
    • G06T2207/10076—4D tomography; Time-sequential 3D tomography
    • 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
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Dentistry (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Vascular Medicine (AREA)
  • Signal Processing (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Pulmonology (AREA)
  • Quality & Reliability (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

一种通过分析感兴趣对象(20)的心脏的至少部分的多幅医学图像来表征心肌灌注病理学的方法，所述多幅医学图像是由医学成像模态(10)以连续的方式采集的，所述方法包括以下步骤：‑对来自所述多幅医学图像的选定心肌图像位置的强度进行采样，并且将表示采集顺序的索引分配到所述心肌图像位置的各自的采样强度，以获得强度曲线(60)；并且‑基于所获得的强度曲线(60)，对指示在多个心肌分段中的心肌分段的至少子集之中的时空灌注非均匀性或灌注失相的索引号(64、66)进行计算；一种用于借助于执行这样的方法通过分析多幅医学图像来进行心肌灌注病理学表征的***(52)；一种包括这样的***(52)的医学成像模态(10)；一种用于执行这样的方法的软件模块(48)。

Description

用于心肌灌注病理学表征的***和方法
技术领域
本发明涉及对心脏的医学成像的领域，并且具体涉及对心脏的医学图像进行分析的领域。
背景技术
由心脏磁共振(CMR)成像对心脏进行的首过增强成像以及近来心脏计算机断层摄影(CCT)成像也允许对心肌灌注的量化。该量化包括对时间强度曲线的半量化分析或真实量化分析。半量化分析包括对时间强度曲线的若干特性特征的量化，所述若干特性特征例如为峰值强度、最大上斜、均值传导时间等。在真实量化分析中，根据对在心肌中获得的动脉输入函数(AIF)和时间强度曲线的数学分析来计算实际的心肌血流量。在Jerosch-Herold的“Quantification of myocardial perfusion by cardiovascular magneticresonance”(Journal of Cardiovascular Magnetic Resonance，2010年，第12卷，第57期)中给出了对半量化方法和真实量化方法两者的全面的综述。
国际申请WO 2005/004066 A1描述了一种用于对心脏灌注的量化评估的方法。所述方法包括将被描绘在一系列心脏图像上的心肌划分成包括至少一个图像像素的图像分段，并且确定针对所述图像分段中的每个的心脏灌注参数。此外，选择具有正常灌注参数值的至少一个图像分段，并且随后，剩下的图像分段的心脏灌注参数是基于正常灌注参数值的。在一个实施例中，灌注参数是针对心肌中的造影剂的分布的时间强度概况的最大上斜。正常最大上斜是针对至少一个图像分段来导出的，并且相对最大上斜是针对与正常最大上斜有关的每个分段来计算的。基于这些值，针对每个分段计算在压迫时导出的心肌灌注参数与在休息时导出的心肌灌注参数的比率。作为范例，其被描述为计算针对每个分段的心肌灌注储备索引(MPRI)，所述MPRI被定义为在休息时与在压迫时导出的相对最大上斜的比率。
发明内容
目前，心肌灌注不足使得存在对应用用于表征和分析的有创方法的要求。无创成像方法不能够可靠地区分冠状动脉微脉管功能障碍(MVD)与多支血管冠心病(包括冠状动脉左总干病；CAD)，这是因为这两种状况都可以导致弥漫性心肌缺血，其在视觉上和量化分析上都能够具有相似的表现。因此，患者经历有创血管造影，并且在具有心肌缺血的患者中证实出正常冠状动脉之后才诊断冠状动脉微脉管功能障碍。
对首过增强图像的分析可以实现对心肌血流量的真实量化。尽管能够通过对心肌血流量的真实量化来分析灌注不足，然而，无创方法仅能够评估缺血的分布和严重性，但是无创方法不能够测量心肌的不同分段中的灌注的时空均匀性。例如，MVD和CAD两者都被表征为严重的且在空间上广泛分布的缺血，通常与到心肌的心内膜层的造影剂的延迟到达相关联。
因此，本发明的目的是提供一种通过分析感兴趣对象的心脏的至少部分的多幅医学图像来表征心肌灌注病理学的改进的方法，所述多幅医学图像是由医学成像模态以连续的方式采集的。
该目的通过包括以下步骤的方法来实现：
-在所述多幅医学图像中描绘所述感兴趣对象的所述心脏的选定部分的轮廓，并且将所述选定部分分割成多个分段；
-对来自所述多幅医学图像的选定心肌图像位置的强度进行采样，并且将表示所述医学图像中的每幅的采集顺序的索引分配到所述心肌图像位置的各自的采样强度，以获得针对所述选定心肌图像位置中的每个的强度曲线；
-基于所获得的强度曲线，对指示在多个心肌分段中的心肌分段的至少子集中间的时空灌注非均匀性或灌注失相的索引号进行计算，其中，心肌分段的所述子集包括多个心肌分段。
在本申请中所使用的术语“心脏的选定部分”应当被具体理解为心脏的任何脉管腔，并且也应当被理解为涵盖左心室和升主动脉。
在本申请中所使用的术语“医学成像模态”应当具体涵盖心脏磁共振(CMR)成像设备、心脏计算机断层摄影(CCT)成像设备、冠状动脉血管造影(CA)设备、CCT血管造影(CCTA)设备、脉管内超声(IVUS)设备、单光子发射计算机断层摄影(SPECT)设备、正电子发射断层摄影(PET)设备以及回波心动描记设备。
在本申请中术语“时空灌注非均匀性”也可以指“时空失相”，并且应当被具体理解为在病理学异常的情况下的非均匀心肌血流量的时间和空间分布。
可以人工地、半自动地或全自动地执行描绘心脏的选定部分的轮廓的步骤。适当的分割技术在现有技术中是已知的、商业可获得的，并且因此在本文中不进行更加详细的讨论。
心肌分段的子集可以包括如灌注地带的多个心肌分段的严格子集，或者其也可以包括完整的心肌。
表示对医学图像中的每一幅的采集顺序的索引可以参考时间尺度，并且可以是时间本身。
所述方法的一个优点在于以下事实：即，在对心脏的选定部分的不同区域的灌注的分布上的额外的时间信息被提供，其能够被用于以快速且简单的方式表征特定的心肌病理学。尤其是，所计算的索引号能够被用于在表征上对冠心病与微脉管病进行区别。
所述方法的另一优点在于能够由无创方法来提供额外的信息，使得在许多情况下能够避免有创血管造影评估。
在应用描绘心脏的选定部分的轮廓的步骤之前，多幅医学图像可以已经经历图像配准技术，以对感性对象的呼吸运动进行校正。适当的配准技术在现有技术中是公知的，并且因此应当不在本文中进行详细讨论。
此外，为了更好的结果，所获得的强度曲线可以已经经历滤波。可以采用适合于本领域技术人员的任何滤波技术。
在优选实施例中，所述方法还包括进行对多个分段中的每个分段中的心肌血流量的量化的步骤。进行对心肌血流量的真实量化的步骤可以根据在背景部分中提出的综述文章中描述的现有技术中的一个(例如，去卷积方法)或根据本领域技术人员认为合适的任何其他技术来执行。对在每个分段中的心肌血流量的真实量化能够提供用于表征心肌灌注病理学的补充信息。
在所述方法的另一实施例中，所述多幅医学图像是在向所述感兴趣对象施予造影剂之后由所述医学成像模态采集的。在本申请中所使用的术语“造影剂”应当被具体理解为当由医学成像模态采集时与基线相比，生成比在药剂周围的感兴趣对象的组织更大的信号的任何药剂。如果医学成像模态基于使用伽马射线的工作原理，则术语“造影剂”也应当涵盖放射性示踪剂材料。以这种方式，能够改进信噪比并且能够实现对心肌灌注病理学的更加确切的表征。
在一个实施例中，可以在想感兴趣对象施予造影剂之后或者在造影剂首过通过心脏期间采集多幅医学图像。
在另一实施例中，可以在向感兴趣对象施予造影剂之后并且在当造影剂的浓度在感兴趣对象中已经达到稳定状态浓度时的时间期间采集多幅医学图像。这对于通过压力回波灌注成像来采集医学图像而言是特别重要的。
在一个实施例中，可以使用如动脉自旋标记的内源性对比度并且在造影剂首过通过心脏期间采集多幅医学图像。
在另一优选实施例中，所述方法还包括识别心脏的选定部分中的参考位置的步骤，其中，在计算索引号的步骤中，参考由所识别的参考位置确定的参考时间来评价强度曲线。这样一来，能够提供用于计算索引号的精确参考。优选的参考位置是左心室，这是因为左心室是在首过期间接收造影剂并且先于针对选定的心肌图像位置中的每个的所获得的强度曲线中的上斜的位置。在本申请中所使用的术语“上斜”应当被具体理解为当所获得的强度曲线的强度由于造影剂的施予而与没有造影剂的情况下的强度相比超过预定阈值(例如，10％)时的时刻。
在另一优选实施例中，所述方法还包括针对选定心肌位置中的每个自动确定相对于由所识别的参考位置确定的参考时间直到心肌图像位置中的每个的采样强度的特性特征的出现的个体时间段的步骤。然后，所述个体时间段被用在计算索引号的步骤中。特性特征可以是强度上升的开始的时间点、峰值强度的时间点，或表现为适合于本领域技术人员的任何其他特性特征。以这种方式，能够以自动方式容易地计算指示时空灌注非均匀性的索引号。
优选地，计算所述索引号的步骤包括计算指示相对于在所识别的参考位置处的特性特征的出现的时间直到在个体心肌位置中的每个处的特性特征的出现的时间的方差的统计学量度。这样一来，能够提供以非常显著的方式描述在心肌分段中间的时空灌注非均匀性的索引号。所述统计学量度可以具有方差的形式，方差***方。总体而言，统计学量度可以具有本领域技术人员认为适合于指示直到在个体心肌位置中的每个处的特性特征的出现的时间的方差的任何其他形式。
优选地，将对所述感兴趣对象的所述心脏的至少部分的所述多幅医学图像的所述采集至少部分地与所述感兴趣对象的所述心脏的循环运动同步。例如，能够在如QRS波群的R峰的心电图中的参考事件之前或之后的固定时间量处采集医学图像。本方法的该实施例的优点在于在心脏的相似状态处获取多幅医学图像中的所有医学图像，使得在医学图像中间的心肌心室几乎没有运动，并且相对静止地绘制心肌。
在本方法的又一实施例中，在对心肌图像位置的强度进行采样的步骤中，所述心肌图像位置在沿着心肌的方向上以及在跨心肌的方向上被选择。以这种方式，计算的索引号以特别适当的方式表示心肌分段中间的时空灌注非均匀性。
所述方法还包括生成灌注图并将灌注图显示到用户的步骤。在本申请中所使用的术语“灌注图”应当被具体理解为根据位置和时间的心肌中的图像位置的强度的颜色表示。灌注图的水平轴可以表示时间，并且垂直轴可以表示心肌中的位置(即，分段)。在较早的出版物中，例如在Marcel Breeuwer的“Comprehensive visualization of first-passmyocardial perfusion:The uptake movie and the perfusogram”(InternationalSociety for Magnetic Resonance in Medicine，ISMRM 2002)中已经将“灌注图”描述为将心肌灌注可视的有用工具。例如，灌注图可以被显示在习惯上用于多幅医学图像已经利用其被采集的医学成像模态的监视器单元上。
此外，所述方法优选地可以包括在灌注图中实施至少一个标记物的步骤，所述至少一个标记物指示至少一个特性位置和/或至少一个特性时间点。能够以快速的方式将被包含在灌注图中的信息传达到用户的这样的标记物的范例是冠状动脉灌注地带的位置、左心室的分段的界线(例如，根据AHA17分段模型)以及如开始时间和峰值时间的造影剂通过的特性时刻。标记物可以被形成为直线、曲线，或者为标记心肌的地带的封闭环。
在另一优选实施例中，所述方法还可以包括实施多个计算机链接的步骤，其中，所述多个计算机链接中的每个计算机链接被分配到灌注图中的位置，并且其中，所述多个计算机链接中的每个计算机链接被链接到表示所述多幅医学图像中的医学图像的数据集。以这种方式，可以将关于心肌血流量的详细信息容易地提供到用户，通常是医学工作人员成员。
本发明的另一目标是提供一种用于通过分析感兴趣对象的心脏的至少部分的多幅医学图像来进行心肌灌注病理学表征的***。所述多幅医学图像是由医学成像模态以连续的方式采集的。所述***包括：
-描绘单元，其被提供用于在所述多幅医学图像中描绘所述感兴趣对象的所述心脏的选定部分的轮廓，并且用于将所述选定部分分割成多个分段；
-强度采样器与分析单元。
所述强度采样器与分析单元被配置为：
-用于对来自所述多幅医学图像的心肌图像位置的强度进行采样，并且将表示所述医学图像中的每幅的采集顺序的索引分配到所述心肌图像位置的各自的采样强度；并且
-用于对指示在多个心肌分段中的心肌分段的至少子集中间的时空灌注非均匀性或灌注失相的索引号进行计算。所述***能够提供如以上针对所述方法公开的相同优点。
在本发明的另一方面中，一种用于心肌灌注病理学表征的***是多幅医学图像已经利用其采集的医学成像模态的集成部分。优选地，所述医学图像模态被设计为磁共振成像装置。所述磁共振成像装置可以有利地包括同步器具，所述同步器具用于将对医学图像的采集与感兴趣对象的心脏的循环运动同步。
在本发明的又一方面中，一种软件模块被提供用于执行上面公开的表征心肌灌注病理学的方法中的任何或者它们的组合，其中，要进行的方法步骤被转换成所述软件模块的程序代码，其中，所述程序代码能够实施在所述医学成像模态的控制单元的存储器单元中，并且能够由所述医学成像模态的所述控制单元的处理器单元运行。
所述控制单元可以是习惯上用于控制医学成像模态的功能的控制单元。所述控制单元备选地可以是被特别分配为运行方法步骤的额外的控制单元。
软件模块能够实现对所述方法的鲁棒且可靠的运行，并且能够允许对方法步骤的快速修改和/或对图像配准算法的调整。
附图说明
本发明的这些和其他方面将根据下文描述的实施例而变得显而易见，并且参考下文描述的实施例得到阐述。
在附图中：
图1是被设计为磁共振成像***的根据本发明的医学成像模态的实施例的部分的示意性图示；
图2是针对在造影剂首过通过控制个体的心脏期间(顶部)，通过具有微脉管功能障碍的个体的心脏期间(中间)，以及通过具有三支血管冠心病的个体的心脏期间(底部)的心肌灌注的根据本发明获得的强度曲线的示意性图示；
图3示出了包括根据依照图2的强度曲线计算的，指示时空灌注非均匀性的所计算的索引号的图；并且
图4图示了灌注图的范例。
附图标记列表
10 医学成像模态
12 磁共振扫描器
14 主磁体
16 检查空间
18 中心轴
20 感兴趣对象
22 磁梯度线圈***
24 全身线圈
26 射频屏蔽
28 控制单元
30 射频发射器单元
32 频切换单元
34 图像处理单元
36 监视器单元
38 心电图设备
40 电极
42 同步单元
44 存储器单元
46 处理器单元
48 软件模块
50 触发信号
52 用于心肌灌注病理学表征的***
54 描绘单元
56 血流分析器
58 强度采样器与分析单元
60 随时间的强度曲线
62 灌注图
64 第一索引号
66 第二索引号
AIF 动脉输入函数
LV 左心室
TTPI 到峰值强度的时间
具体实施方式
图1示出了根据本发明的医学成像模态10的实施例的部分的示意性图示，所述医学成像模态被设计为磁共振成像***，用于采集感兴趣人类对象20的心脏的至少部分的医学图像(由磁共振图像表示)。磁共振成像***包括具有主磁体14的磁共振扫描器12，所述主磁体被提供用于生成静态磁场。主磁体14具有中心膛，所述中心膛提供用于要被定位在其内的感兴趣对象20的围绕中心轴18的检查空间16。出于简明的原因，在图1中已经省略了用于支撑感兴趣对象20的常规诊查台。大致静态的磁场定义检查空间16的轴向方向，所述轴向方向平行于中心轴18被对准。此外，磁共振成像***10包括磁梯度线圈***22，所述磁梯度线圈***被提供用于生成被叠加到静态磁场的梯度磁场。如现有技术已知的，磁梯度线圈***22被同轴地布置在主磁体14的膛内。
在原理上，本发明也适用于提供在静态磁场内的检查区域的任何其他类型的磁共振成像***。此外，应当认识到，本发明能够与被配置为采集感兴趣对象的心脏的至少部分的医学图像的任何其他医学成像模态一起使用。本发明能够被应用于的医学成像模态的范例为心脏计算机断层摄影(CCT)成像设备、冠状动脉血管造影(CA)设备、CCT血管造影(CCTA)设备、血管内超声(IVUS)设备、单光子发射计算机断层摄影(SPECT)设备、正电子发射断层摄影(PET)设备以及心电图设备。
此外，磁共振成像***包括被设计为全身线圈24的射频天线，所述全身线圈被提供用于在射频发射阶段期间将射频电磁场施加到检查空间16，以激励感兴趣对象20中或内的原子核。全身线圈24也被提供为在射频接收阶段期间接收来自感兴趣对象20中或内的被激励的原子核的磁共振信号。在磁共振成像***的操作状态中，射频发射阶段和射频接收阶段是以连续的方式发生的。全身线圈24具有中心轴，并且，在操作状态中，全身线圈24被同轴地布置在主磁体14的膛内，使得全身线圈24的中心轴与磁共振成像***的中心轴18重合。如本领域所公知的，圆柱形金属射频屏蔽26被同轴地布置在磁梯度线圈***22与全身线圈24之间。
磁共振成像***还包括控制单元28，所述控制单元被提供用于至少控制磁共振扫描器12和磁梯度线圈***22的功能。控制单元28包括常规的监视器单元36。
此外，磁共振成像***包括射频发射器单元30，所述射频发射器单元被连接到控制单元28并由控制单元28控制。射频发射器单元30被提供为在射频发射阶段期间经由射频切换单元32将磁共振射频的射频功率馈送到全身线圈24。在射频接收阶段期间，射频切换单元32将来自全身线圈24的磁共振信号引导到驻留在控制单元28中的图像处理单元36。图像处理单元34被配置用于处理所采集的磁共振信号，以根据所采集的磁共振信号生成感兴趣对象20的部分的磁共振图像。对于本领域技术人员而言，该技术的许多不同变型是公知的，并且因此不需要在本文中更加详细地描述。
为了采集感兴趣对象20的心脏的磁共振图像，磁共振成像***包括同步器具，所述同步器具用于将对医学图像的采集与感兴趣对象的心脏的循环运动同步。同步器具被形成为心电图设备38和同步单元42。
心电图设备38被提供用于获取感兴趣对象20的心脏的心电图数据的测量结果。对此，心电图设备38的多个电极40可以被布置在感兴趣对象20处。此外，心电图设备38包括用于对心电图数据进行滤波以减少由磁梯度场生成的伪影的器具。合适的滤波器具对于本领域技术人员而言是已知的，并且因此将不在本文中更加详细地描述。
心电图设备38被耦合到同步单元42，所述同步单元被配置用于生成触发信号50，以根据对心脏活动的QRS波群(ORS complex)的R峰的检测来触发采集磁共振信号的采集时期。同步单元42继而被耦合到控制单元28。控制单元28被配置为通过由同步单元42提供的触发信号50来被同步，以用于生成针对生成梯度磁场的磁梯度线圈***22的控制信号。对此，控制单元28被配置为在接收触发信号50后生成多个序列，每个序列包括射频场和磁梯度场。
医学成像模态10包括用于通过分析感兴趣对象20的心脏的至少部分的多幅医学图像来进行心肌灌注病理学表征的***52。用于心肌灌注病理学表征的***52驻留在控制单元28的壳体内，并且包括描绘单元54、血流分析器56以及强度采样器与分析单元58。将在下文中详细解释这些设备的功能和交互。
在向感兴趣对象20施予造影剂之后，在造影剂首过通过感兴趣对象20的心脏期间，由磁共振成像***以连续的方式采集采取感兴趣对象20的心脏的磁共振图像的形式的多幅医学图像。已经分别针对三种不同感兴趣对象采集了多幅医学图像：控制个体、具有三支血管冠心病的个体，以及具有微脉管功能障碍的个体。已经将对感兴趣对象20的心脏的多幅医学图像的采集与感兴趣对象20的各自的心脏的循环运动同步，如上所述。
在对各自的感兴趣对象20的心脏的每种多幅医学图像的采集之后，由根据本方法的表征心肌灌注病理学的方法的实施例来分析医学图像。在以下所述方法被详细描述为被应用于一种多幅医学图像。应当理解，所述方法也以同样的方式适用于其他两种多幅医学图像。
为了能够执行所述方法，用于心肌灌注病理学表征的***52包括软件模块48(图1)。要被进行的方法步骤被转换成软件模块48的程序代码，其中，程序代码可实施在用于心肌灌注病理学表征的***52的存储器单元44中，并且由用于心肌灌注病理学表征的***52的处理器单元46运行。
在预备步骤中，将多幅医学图像中的每幅医学图像服从图像配准算法。
在所述方法的第一步骤中，由用于心肌灌注病理学表征的***52的描绘单元54在多幅医学图像中描绘感兴趣对象20的心脏的左心室的轮廓，并且将左心室分割成多个分段。
在所述方法的下一步骤中，用于心肌灌注病理学表征的***52的血流分析器56在多个分段中的每个分段中进行对心肌血流量的真实量化。
然后，由用于心肌灌注病理学表征的***52的强度采样器与分析单元58选择心肌图像位置。在沿着心肌的方向上以及在跨心肌的方向上选择心肌图像位置。强度采样器与分析单元58然后对来自多幅医学图像的选定的心肌图像位置的强度进行采样，并且将表示由对医学图像中的每一幅的采集时间形成的采集顺序的索引分配到心肌图像位置的各自的采样强度，以获得针对选定的心肌图像位置中的每个的随时间的强度曲线60。
在图2中示意性地示出了针对三种不同感兴趣对象的所获得的随时间的强度曲线60(顶部：控制个体，中间：具有微脉管功能障碍(MVD)的个体，底部：具有三支血管冠心病(CAD)的个体)。
针对图2中的随时间的强度曲线60，为了提供合适的时间尺度，由强度采样器与分析单元58识别左心室中的参考位置。参考位置是在造影剂首过期间并在选定的心肌图像位置的强度的上升的开始之前首次接收造影剂的位置。参考位置处的强度的上升的开始的时间点充当参考时间。参考所识别的参考位置和参考时间来评价随时间的强度曲线60。
在所述方法的下一步骤中，针对选定的心肌位置中的每个，自动确定相对于所识别的参考位置和对应的参考时间直到在心肌图像位置中的每个的采样强度曲线60中的特性特征的出现的个体时间段。在该实施例中，特性特征是随时间的强度曲线的峰值强度。个体到峰值强度时间段TTPI被用在计算索引号64、66的下一步骤中。
在所述方法的后续步骤中，强度采样器与分析单元58基于所获得的随时间的强度曲线60来计算第一索引号64，所述第一索引号指示在多个心肌分段中的心肌分段的子集中间的时空灌注非均匀性或失相。
指示左心室灌注的时间失相的第一索引号64在下一步骤中被计算为在参考时间与直到随时间的强度曲线60的峰值强度出现的时间之间的到峰值强度时间段TTPI的方差。在本文中，方差被理解为时间段TTPI的标准差的平方。
指示左心室灌注的时间失相的第二索引号66被计算为参考时间与直到随时间的强度曲线60的峰值强度出现的时间之间的到峰值强度TTPI的时间段的变异系数。
针对在图3中显示的计算的索引号64、66的结果清楚地示出了控制个体(CTRL，左列)、具有三支血管冠心病的个体(CAD，中间列)以及具有微脉管功能障碍的个体(MVD，右列)之间的显著差异。
根据图2中的随时间的强度曲线60变得明显的是，MVD和三支血管CAD两者都被表征为严重的且在空间上广泛分布的缺血，通常与到心肌的心内膜层的造影剂的延迟到达相关联。然而，添加时空失相分析能够提供关于对左心室的不同区域的灌注的时间分布的额外信息。图3示出了计算的索引号64、66构成CAD与MVD之间的因不同潜在病理生理学的主要差异，并且允许特定表征。
在具有CAD的个体中，心肌中的心肌血流量的时空分布是越来越不均匀的。而在具有MVD的个体中，病理改变涉及微观循环及其在心缩期中与心肌收缩的相互作用。在这种情况下，对心外膜层的灌注未被堵塞并且非常均匀。然而，在首过波的透壁传播中存在延迟。类似于CAD，这引起具有延迟的强度上升开始的广泛分布缺血。然而，该特征在心肌上在时间域中是均匀的，从而允许在CAD与MVD之间进行无创区分。
在所述方法的另一步骤中，灌注图62被生成并由监视器单元36显示到用户。图4图示了具有灰度级或颜色编码的灌注图62的范例。能够在灌注图62(未示出)中实施指示特性位置和/或特性时间点的标记物，如直线。在灌注图62中也能够***其他标记物(未示出)，所述其他标记物被形成为闭合环，并且指示在其中心肌被较不好地灌注的时间和空间方案。
在所述方法的又一步骤中，在灌注图62中实施多个计算机链接。多个计算机链接中的每个计算机链接被分配到灌注图62中的位置，并且多个计算机链接中的每个计算机链接被链接到表示多幅医学图像中的医学图像的数据集。图4中示出的计算机鼠标光标正指向灌注图62的表示特定时间点和心肌中的特定位置的位置。通过敲击鼠标按钮，激活到表示对应于时间点和心肌的位置的医学图像的数据集的计算机链接，并且将所述医学图像显示在磁共振成像***的监视器单元36上。
尽管已经在附图和前面的描述中详细图示和描述了本发明，但是这样的图示和描述应当被认为是说明性或示范性的，而非限制性的；本发明不限于所公开的实施例。本领域技术人员通过研究附图、公开内容以及权利要求，在实践请求保护的发明时能够理解并实现对所公开的实施例的其他变型。在权利要求书中，“包括”一词不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。尽管特定措施被记载在互不相同的从属权利要求中，但是这并不指示不能有利地使用这些措施的组合。权利要求中的任何附图标记都不应被解释为对范围的限制。

Claims (24)

1.一种用于通过分析感兴趣对象(20)的心脏的至少部分的多幅医学图像来进行心肌灌注病理学表征的***(52)，所述多幅医学图像是由医学成像模态(10)以连续的方式采集的，所述***(52)包括：
-描绘单元(54)，其被提供用于在所述多幅医学图像中描绘所述感兴趣对象(20)的所述心脏的选定部分的轮廓，并且用于将所述选定部分分割成多个分段；
-强度采样器与分析单元(58)，其被配置：
-用于对来自所述多幅医学图像的心肌图像位置的强度进行采样，并且将表示所述医学图像中的每幅的采集顺序的索引分配到所述心肌图像位置的各自的采样强度以获得针对所述心肌图像位置中的每个的强度曲线(60)；并且
-用于基于所获得的强度曲线(60)，对指示在多个心肌分段中的心肌分段的至少子集之中的时空灌注非均匀性或灌注失相的一个索引号(64、66)进行计算。
2.如权利要求1所述的***，其中，
实行对在所述多个分段中的每个分段中的心肌血流量的量化。
3.如权利要求1或2所述的***，其中，所述多幅医学图像是在向所述感兴趣对象(20)施予造影剂之后由所述医学成像模态(10)采集的。
4.如权利要求1或2所述的***，其中，
识别在所述心脏的所述选定部分中的参考位置，其中，在对所述索引号(64、66)的计算中，参考由所识别的参考位置确定的参考时间来评价所述强度曲线(60)。
5.如权利要求4所述的***，其中，
针对选定心肌位置中的每个，自动确定相对于由所识别的参考位置确定的参考时间直到所述心肌图像位置中的每个的所述采样强度的特性特征的出现的个体时间段(TTPI)，
其中，所述个体时间段(TTPI)被用在对所述索引号(64、66)的计算中。
6.如权利要求4所述的***，其中，计算所述索引号(64、66)包括计算指示相对于在所识别的参考位置处的特性特征的出现的时间直到在个体心肌位置中的每个处的特性特征的出现的时间(TTPI)的方差的统计学量度。
7.如权利要求1或2所述的***，其中，将对所述感兴趣对象(20)的所述心脏的至少部分的所述多幅医学图像的所述采集至少部分地与所述感兴趣对象(20)的所述心脏的循环运动同步。
8.如权利要求1或2所述的***，其中，在对心肌图像位置的强度的采样中，所述心肌图像位置在沿着心肌的方向上以及在横跨心肌的方向上被选择。
9.如权利要求1或2所述的***，其中，生成灌注图(62)并且将所述灌注图显示到用户。
10.如权利要求1或2所述的***，其中，在灌注图(62)中实施至少一个标记物，所述至少一个标记物指示至少一个特性位置和/或至少一个特性时间点。
11.如权利要求1或2所述的***，其中，实施多个计算机链接，其中，所述多个计算机链接中的每个计算机链接被分配到灌注图(62)中的位置，并且其中，所述多个计算机链接中的每个计算机链接被链接到表示所述多幅医学图像中的医学图像的数据集。
12.一种包括如权利要求1所述的***(52)的医学成像模态(10)。
13.如权利要求12所述的医学成像模态(10)，其被设计为磁共振成像装置。
14.一种计算机可读介质，其上存储有程序代码，所述程序代码能由医学成像模态(10)的控制单元(52)的处理器单元(46)运行以执行一种通过分析感兴趣对象(20)的心脏的至少部分的多幅医学图像来表征心肌灌注病理学的方法，所述多幅医学图像是由医学成像模态(10)以连续的方式采集的，所述方法包括以下步骤：
-在所述多幅医学图像中描绘所述感兴趣对象(20)的所述心脏的选定部分的轮廓，并且将所述选定部分分割成多个分段；
-对来自所述多幅医学图像的选定心肌图像位置的强度进行采样，并且将表示所述医学图像中的每幅的采集顺序的索引分配到所述心肌图像位置的各自的采样强度，以获得针对所述选定心肌图像位置中的每个的强度曲线(60)；
-基于所获得的强度曲线(60)，对指示在多个心肌分段中的心肌分段的至少子集之中的时空灌注非均匀性或灌注失相的一个索引号(64、66)进行计算，其中，心肌分段的所述子集包括多个心肌分段。
15.如权利要求14所述的计算机可读介质，所述方法还包括以下步骤：
-实行对在所述多个分段中的每个分段中的心肌血流量的量化。
16.如权利要求14或15所述的计算机可读介质，其中，所述多幅医学图像是在向所述感兴趣对象(20)施予造影剂之后由所述医学成像模态(10)采集的。
17.如权利要求14或15所述的计算机可读介质，所述方法还包括以下步骤：
-识别在所述心脏的所述选定部分中的参考位置，其中，在计算所述索引号(64、66)的步骤中，参考由所识别的参考位置确定的参考时间来评价所述强度曲线(60)。
18.如权利要求17所述的计算机可读介质，所述方法还包括以下步骤：
-针对选定心肌位置中的每个，自动确定相对于由所识别的参考位置确定的参考时间直到所述心肌图像位置中的每个的所述采样强度的特性特征的出现的个体时间段(TTPI)，
其中，所述个体时间段(TTPI)被用在计算所述索引号(64、66)的步骤中。
19.如权利要求17所述的计算机可读介质，其中，计算所述索引号(64、66)的步骤包括计算指示相对于在所识别的参考位置处的特性特征的出现的时间直到在个体心肌位置中的每个处的特性特征的出现的时间(TTPI)的方差的统计学量度。
20.如权利要求14或15所述的计算机可读介质，其中，将对所述感兴趣对象(20)的所述心脏的至少部分的所述多幅医学图像的所述采集至少部分地与所述感兴趣对象(20)的所述心脏的循环运动同步。
21.如权利要求14或15所述的计算机可读介质，其中，在对心肌图像位置的强度进行采样的步骤中，所述心肌图像位置在沿着心肌的方向上以及在横跨心肌的方向上被选择。
22.如权利要求14或15所述的计算机可读介质，所述方法还包括生成灌注图(62)并且将所述灌注图显示到用户的步骤。
23.如权利要求14或15所述的计算机可读介质，所述方法还包括在灌注图(62)中实施至少一个标记物的步骤，所述至少一个标记物指示至少一个特性位置和/或至少一个特性时间点。
24.如权利要求14或15所述的计算机可读介质，所述方法还包括实施多个计算机链接的步骤，其中，所述多个计算机链接中的每个计算机链接被分配到灌注图(62)中的位置，并且其中，所述多个计算机链接中的每个计算机链接被链接到表示所述多幅医学图像中的医学图像的数据集。
CN201480054511.XA 2013-10-01 2014-10-01 用于心肌灌注病理学表征的***和方法 Active CN105848577B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP13186934.9 2013-10-01
EP13186934 2013-10-01
PCT/EP2014/071109 WO2015049324A1 (en) 2013-10-01 2014-10-01 System and method for myocardial perfusion pathology characterization

Publications (2)

Publication Number Publication Date
CN105848577A CN105848577A (zh) 2016-08-10
CN105848577B true CN105848577B (zh) 2019-11-05

Family

ID=49322196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480054511.XA Active CN105848577B (zh) 2013-10-01 2014-10-01 用于心肌灌注病理学表征的***和方法

Country Status (7)

Country Link
US (1) US10806367B2 (zh)
EP (1) EP3052019B1 (zh)
JP (1) JP6382960B2 (zh)
CN (1) CN105848577B (zh)
BR (1) BR112016006914A2 (zh)
RU (1) RU2677009C2 (zh)
WO (1) WO2015049324A1 (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11415655B2 (en) * 2016-05-27 2022-08-16 University Of Virginia Patent Foundation Reduced field-of-view perfusion imaging with high spatiotemporal resolution
DE102016117889B3 (de) * 2016-09-22 2018-03-15 Tomtec Imaging Systems Gmbh Verfahren und Vorrichtung zur Korrektur von durch Tracking-Verfahren ermittelten dynamischen Modellen
GB201706149D0 (en) * 2017-04-18 2017-05-31 King's College London System and method for medical imaging
EP3437559B1 (de) * 2017-08-03 2020-06-24 Siemens Healthcare GmbH Ermitteln eines funktionsparameters betreffend eine lokale gewebefunktion für mehrere gewebebereiche
JP7002901B2 (ja) * 2017-09-28 2022-02-10 キヤノンメディカルシステムズ株式会社 医用画像診断装置及び医用画像処理装置
CN110428420B (zh) * 2018-09-05 2022-05-17 深圳科亚医疗科技有限公司 基于患者的冠状动脉ct血管造影图像来确定冠状动脉的流动信息的方法、装置和介质
DE102020213938A1 (de) * 2020-11-05 2022-05-05 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Störunterdrückung für MR-Ganzkörperantennen

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7024024B1 (en) * 2000-11-14 2006-04-04 Axle International System for contrast echo analysis
CN102301393A (zh) * 2009-01-29 2011-12-28 皇家飞利浦电子股份有限公司 将生物结构中的时变参数可视化

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7120293B2 (en) * 2001-11-30 2006-10-10 Microsoft Corporation Interactive images
EP1644900B1 (en) * 2003-07-01 2011-10-19 Koninklijke Philips Electronics N.V. Non-invasive quantitative myocardial perfusion assessment
JP4537681B2 (ja) 2003-09-24 2010-09-01 株式会社東芝 血流解析装置
JP4446049B2 (ja) * 2003-11-06 2010-04-07 株式会社三重ティーエルオー 心筋血流の定量化装置
JP4441588B2 (ja) * 2004-03-26 2010-03-31 株式会社三重ティーエルオー 心筋梗塞判定用画像作成装置
JP4599923B2 (ja) * 2004-07-16 2010-12-15 ソニー株式会社 情報処理システム及び情報処理方法、並びにコンピュータプログラム
DE102004055461A1 (de) * 2004-11-17 2006-05-04 Siemens Ag Bildgebendes Verfahren sowie Vorrichtung zur Visualisierung von koronaren Herzkrankheiten
US7587232B2 (en) * 2006-02-28 2009-09-08 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus, magnetic resonance data processing apparatus, magnetic resonance data processing program and magnetic resonance imaging apparatus control method
JP5405046B2 (ja) * 2007-06-28 2014-02-05 株式会社東芝 画像処理装置、画像診断装置および画像処理方法
JP5234905B2 (ja) * 2007-11-20 2013-07-10 東芝メディカルシステムズ株式会社 Ｘ線ｃｔ装置および心筋パーフュージョン像生成システム
JP5361410B2 (ja) 2009-01-22 2013-12-04 株式会社東芝 画像処理装置
JP5658686B2 (ja) 2009-01-29 2015-01-28 コーニンクレッカ フィリップス エヌ ヴェ 貫壁性のかん流の勾配の画像分析
JP5388638B2 (ja) 2009-03-13 2014-01-15 株式会社東芝 画像処理装置および画像処理プログラム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7024024B1 (en) * 2000-11-14 2006-04-04 Axle International System for contrast echo analysis
CN102301393A (zh) * 2009-01-29 2011-12-28 皇家飞利浦电子股份有限公司 将生物结构中的时变参数可视化

Also Published As

Publication number Publication date
WO2015049324A1 (en) 2015-04-09
US10806367B2 (en) 2020-10-20
JP2016531630A (ja) 2016-10-13
EP3052019B1 (en) 2018-03-07
JP6382960B2 (ja) 2018-08-29
CN105848577A (zh) 2016-08-10
BR112016006914A2 (pt) 2017-09-12
RU2677009C2 (ru) 2019-01-14
RU2016116785A (ru) 2017-11-13
RU2016116785A3 (zh) 2018-06-08
EP3052019A1 (en) 2016-08-10
US20160235330A1 (en) 2016-08-18

Similar Documents

Publication Publication Date Title
CN105848577B (zh) 用于心肌灌注病理学表征的***和方法
US20170127972A1 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
Kim et al. Myocardial tissue tracking with two-dimensional cine displacement-encoded MR imaging: development and initial evaluation
US6438403B1 (en) Method and apparatus for cardiac analysis using four-dimensional connectivity
EP1182613B1 (en) Diagnostic imaging
Barkhausen et al. Assessment of ventricular function with single breath-hold real-time steady-state free precession cine MR imaging
US9295406B2 (en) Automatic or semi-automatic whole body MR scanning system
US20190274569A1 (en) Method for determining diastasis timing using an mri septal scout
US7764814B2 (en) Display and analysis of multicontrast-weighted magnetic resonance images
Gensler et al. Myocardial T1: quantification by using an ECG-triggered radial single-shot inversion-recovery MR imaging sequence
US8315450B2 (en) Method and system for display of medical image data
Penzkofer et al. Assessment of myocardial perfusion using multisection first-pass MRI and color-coded parameter maps:: A comparison to 99mTc Sesta MIBI SPECT and systolic myocardial wall thickening analysis
US8036730B1 (en) Temporal magnetic resonance imaging
De Feyter et al. Non-invasive coronary artery imaging with electron beam computed tomography and magnetic resonance imaging
JP2017170146A (ja) 画像処理装置、磁気共鳴イメージング装置及び画像処理方法
US10557905B2 (en) System and method for producing high-resolution magnetic relaxation parameter maps
Wise et al. Magnetic resonance imaging analysis of left ventricular function in normal and spontaneously hypertensive rats
Lin et al. Bi-ventricular assessment with cardiovascular magnetic resonance at 5 Tesla: A pilot study
US11497412B2 (en) Combined oxygen utilization, strain, and anatomic imaging with magnetic resonance imaging
CN108363027B (zh) 采集心脏的诊断测量数据和检查心脏的方法及磁共振设备
Eggen et al. Cardiovascular Magnetic Resonance Imaging and MR-Conditional Cardiac Devices
Robert et al. Uncalibrated Real-Time Stroke Volume Estimation in MRI Using the Magnetohydrodynamic Effect?
Ginami Respiratory Motion Compensation in Coronary Magnetic Resonance Angiography: Analysis and Optimization of Self-Navigation
Klein et al. Magnetic resonance flow measurements in real-time: Comparison with a standard gradient-echo
Fluckiger Absolute quantification of pharmaco-kinetic parameters in dynamic contrast enhanced magnetic resonance imaging

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