CN109091145A - 心肌定量磁共振成像方法、设备及存储介质 - Google Patents

心肌定量磁共振成像方法、设备及存储介质 Download PDF

Info

Publication number
CN109091145A
CN109091145A CN201810923440.2A CN201810923440A CN109091145A CN 109091145 A CN109091145 A CN 109091145A CN 201810923440 A CN201810923440 A CN 201810923440A CN 109091145 A CN109091145 A CN 109091145A
Authority
CN
China
Prior art keywords
signal
heart
acquisition
picture signal
bat
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
Application number
CN201810923440.2A
Other languages
English (en)
Other versions
CN109091145B (zh
Inventor
丁海艳
郭瑞
罗建文
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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201810923440.2A priority Critical patent/CN109091145B/zh
Publication of CN109091145A publication Critical patent/CN109091145A/zh
Priority to PCT/CN2019/086400 priority patent/WO2020034676A1/zh
Application granted granted Critical
Publication of CN109091145B publication Critical patent/CN109091145B/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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/023—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the heart

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Signal Processing (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Radiology & Medical Imaging (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提供了一种心肌定量磁共振成像方法、设备和存储介质。该方法包括：每隔恢复时间段，在心电门控信号和呼吸导航信号的控制下，执行图像信号采集操作；根据所采集的图像信号以及其对应的饱和脉冲的延迟时间确定参数T1；根据所述参数T1生成心肌定量磁共振图像。该方案可以在受试者自由呼吸的情况下完成扫描，而无需闭气。同时也允许进一步扩大成像视野，提高空间分辨率。此外，通过k‑空间在各采样点间完全交错分段采集，从而实现原始图像的内在配准，无需在后期进行另外的图像处理。

Description

心肌定量磁共振成像方法、设备及存储介质
技术领域
本发明涉及医疗成像领域，更具体地涉及一种心肌定量磁共振成像方法、设备及存储介质。
背景技术
核磁共振成像技术利用核磁共振现象对人体成像，已经是一种常见的医学影像检查方式。
核磁共振的基本物理参数T1(自旋晶格弛豫时间，或称作纵向弛豫时间)是描述纵向磁化向量恢复的时间常数。T1由生物组织的组成成分和存在的结构形式以及磁场强度决定。在确定的磁场强度下，不同组织具有特定的T1数值。当生物组织发生改变，T1也会随之发生改变。因此T1可以作为特征参数识别心肌组织特征。
近年来发展起来的心肌T1定量磁共振成像技术，对磁共振的基本物理参数T1直接测量，实现定量的心肌组织评价。
在现有心肌T1定量心脏成像的技术中，多数是通过屏气实现扫描过程中对呼吸运动的补偿。这种技术需要较多的T1加权采样点才能够成像，对心率变化敏感。屏气的要求制约了成像分辨率的进一步提高，并且无法用于屏气困难的受试者(这在心脏病患者中比较常见)。采用呼吸导航操作可以突破屏气的限制，进一步采用饱和准备脉冲实现T1加权，可以降低对心率变化的依赖性。但现有的采用呼吸导航操作的技术需要采集的采样点较多。因此，存在扫描时间长、不同采样点间存在错配的问题。并且，由于部分T1加权采样点的信噪比比较低，影响最终的T1拟合质量。如果减少采样点数量，就必须在拟合参数前，对原始加权图像进行滤波处理。甚至还需要对原始加权图像间进行运动校正(例如配准)来消除心肌运动对成像的负面影响。
因此，迫切需要一种新的心肌定量磁共振成像技术，以至少部分地解决上述问题。
发明内容
考虑到上述问题而提出了本发明。
根据本发明一个方面，提供了一种心肌定量磁共振成像方法，包括：
每隔恢复时间段，在心电门控信号和呼吸导航信号的控制下，至少执行如下信号采集操作：
在第一心拍内，在根据呼吸导航信号判断当前时刻符合预定条件的情况下，采集第一图像信号；
在第二心拍内，在使用延迟时间为Tsat2的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第二图像信号；
在第三心拍内，在使用延迟时间为Tsat3的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第三图像信号，其中，Tsat3≠Tsat2；
根据第i图像信号以及第i图像信号对应的饱和脉冲的延迟时间Tsati确定参数T1，其中，i＝1、2、3，且当i＝1时，Tsati为无穷大；
根据所述参数T1生成心肌定量磁共振图像。
根据本发明另一方面，还提供了一种用于心肌定量磁共振成像的设备，包括处理器和存储器，其中，所述存储器中存储有计算机程序指令，所述计算机程序指令被所述处理器运行时用于执行上述心肌定量磁共振成像方法。
根据本发明再一方面，还提供了一种存储介质，在所述存储介质上存储了程序指令，所述程序指令在运行时用于执行上述心肌定量磁共振成像方法。
根据本发明实施例的心肌定量磁共振成像方法、设备及存储介质可以在受试者自由呼吸的情况下完成扫描，而无需闭气。同时也允许进一步扩大成像视野，提高空间分辨率。此外，通过k-空间在各采样点间完全交错分段采集，从而实现原始图像的内在配准，无需在后期进行另外的图像处理。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
通过结合附图对本发明实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显。附图用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与本发明实施例一起用于解释本发明，并不构成对本发明的限制。在附图中，相同的参考标号通常代表相同部件或步骤。
图1示出了根据本发明一个实施例的心肌定量磁共振成像方法的示意性流程图；
图2示出了根据本发明一个实施例的参数T1定量图像；
图3示出了根据本发明一个实施例的成像序列；
图4示出了根据本发明一个实施例的T1估计曲线；
图5示出了根据本发明一个实施例的所采集的受试者的短轴视角下的同一层心脏的第一图像信号至第三图像信号以及拟合得到的参数T1定量图像；
图6示出了根据本发明另一个实施例的成像序列；
图7a示出了根据本发明一个实施例的3D左心室的T1图像；
图7b分别示出了图7a所示的T1图像在长轴视角下的图像；
图7c示出了图7a所示的T1图像中左心室心肌内T1的统计直方图；以及
图8示出了根据本发明一个实施例的在使用T1缩短对比剂后的参数T1定量图像。
具体实施方式
为了使得本发明的目的、技术方案和优点更为明显，下面将参照附图详细描述根据本发明的示例实施例。显然，所描述的实施例仅仅是本发明的一部分实施例，而不是本发明的全部实施例，应理解，本发明不受这里描述的示例实施例的限制。基于本发明中描述的本发明实施例，本领域技术人员在没有付出创造性劳动的情况下所得到的所有其它实施例都应落入本发明的保护范围之内。
根据本发明的实施例，提供了一种心肌定量磁共振成像方法。该方法是一种3D自由呼吸式的定量心肌参数T1成像技术。该技术采用呼吸导航技术实现对呼吸运动的补偿。通过饱和脉冲的使用，结合充分的T1恢复时间获得理想的稳态磁化向量，对心率变化不敏感，可以实现较高的T1拟合准确度。通过k-空间交错分段采集的方式实现原始图像内在的自配准，无需原始图像的配准、滤波等后处理，也不需要参数校正，可用于高场强(3T)3D心肌定量参数T1的准确测量。
图1示出了根据本发明一个实施例的心肌定量磁共振成像方法100的示意性流程图。如图1所示，心肌定量磁共振成像方法100包括以下步骤。
步骤S110，每隔恢复时间段，在心电门控信号和呼吸导航信号的控制下，执行信号采集操作。
该信号采集操作可以以循环的方式采集多个成像序列。每循环一次，采集一个成像序列。每个成像序列中都包括多个图像信号。在一个示例中，每个成像序列包括3个图像信号。k-空间是磁共振采集的数据空间。每完成一个成像序列的图像信号采集即实现了磁共振成像中k-空间的一个分段的填充。需要循环采集成像序列若干次，才能填充完整的k-空间，以用于重建图像。可以理解，可以采用并行采样技术以及其他任何方式的k-空间降采技术。
信号采集操作以心电图(ECG)为基础。可以通过在受试者胸部皮肤表面贴电极并且通过心电监测设备获取心电图。在心电图中，两个R波之间的时间间隔称为心拍(Beat)，即心动周期。可以通过检测R波来确定下一个心拍。成像序列中的每个图像信号都是分别在一个心动周期内采集的。在上述成像序列包括3个图像信号的示例中，需要3个心拍来完成信号采集。
在成像序列中存在一个允许纵向磁化向量充分恢复到稳态的图像信号。为了描述简洁，简称其为稳态图像信号。为了使得纵向磁化矢量能够在上一个成像序列的采集操作之后完全恢复，在稳态图像信号采集前，设置一个恢复时间段，例如n个空闲心拍，或称为恢复心拍。该恢复时间段中不采集图像信号，也不做任何可能扰动磁化向量的恢复过程的操作。
前述空闲心拍个数n可以根据允许磁化向量恢复的时间长度和受试者的心率确定。允许磁化向量恢复的时间长度N(秒)可以由用户在线实时设定，其间接决定了稳态图像信号在T1估计曲线上的时间位置。N越大，磁化向量恢复的越充分。由此，能获得更理想的稳态数据，有利于提高参数T1的准确度。但是N越大，成像所需的扫描时间也越长。因此，需要基于扫描效率和稳态值的准确度二者来设定N。在一个示例中，磁场强度为3T，最小空闲时间可以设定为6秒(即N＝6)。这样能够保证恢复95％以上的磁化向量。如果磁场强度为1.5T，或者是在使用T1缩短对比剂后，N可以相应地减小。通过在线实时设定允许磁化向量恢复的时间长度N而不是空闲心拍个数n，可以在保证信号恢复到稳态的程度的同时，剥离与心率变化的关联。
根据本发明一个实施例，空闲心拍个数n≥N/(60/HR)，其中受试者的心率为HR(心跳/分钟)，允许磁化向量恢复的时间长度为N秒。可选地，n取大于或等于N/(60/HR)的最小整数。利用该公式确定的空闲心拍个数能够保证磁化向量能够充分地恢复。进而保证所生成的心肌定量磁共振图像的准确性。
在每个心拍中，根据心电门控信号，确定采集图像信号的时刻。自R峰开始经过时间段Ttrigger之后，即开始采集图像信号。期望采集图像信号的时刻是心脏相对静止的时刻，例如心脏舒张期末期的一个时刻。通过心电门控技术可以使得所采集的图像信号较少地受到心脏运动的干扰。可以理解，时间段Ttrigger可以由扫描人员根据经验进行设置。
根据本发明的实施例，在每个心拍中，还根据呼吸导航信号来确定是否采集图像信号。通过监测胸膈肌位置随呼吸运动的改变，能够间接估计心脏随呼吸运动的位置改变。在本发明的实施例中，在自心电门控信号的R峰开始经过时间段Ttrigger之前的一小段时间(NAV)内，采集呼吸导航信号。根据该呼吸导航信号判断当前时刻是否符合预定条件，即在当前时刻胸膈肌位置是否在期望的位置。由此判断在自心电门控信号的R峰开始经过时间段Ttrigger后采集的图像信号是否符合呼吸运动补偿的要求，也即确定在本心拍内采集的图像信号是否有效。利用呼吸导航技术，使得在心肌定量磁共振成像过程中，受试者能够自由呼吸。而且也扩大了成像视野，提高了图像的空间分辨率。
步骤S120，根据步骤S110所采集的图像信号分别对应的饱和脉冲的延迟时间确定参数T1。在采集图像信号之前，可以施加饱和脉冲。可以通过施加不同延迟时间的饱和脉冲来改变参数T1的权重。对于不施加饱和脉冲的心拍，可以设该心拍内所采集的图像信号对应延迟时间为无穷大的饱和脉冲。饱和脉冲的使用降低了对图像信号对心率变化的依赖性。
在该步骤中，可以采用不同的信号模型来基于采样点(即图像信号)拟合参数T1。
步骤S130，根据步骤S120所确定的参数T1生成心肌定量磁共振图像。在此步骤中，可以根据通过拟合操作所获得的参数T1生成心肌定量磁共振图像。
图2示出了根据本发明一个实施例的、来自三个受试者的参数T1定量的图像。对于每一个受试者，从心尖、心中、心基分别选一层作为代表图像。
根据本发明的实施例的上述成像方法100，对于受试者的心脏运动和呼吸运动，通过心电门控和呼吸导航来实现对上述运动的补偿，由此保证图像信号在相同的呼吸状态和心脏运动周期上采集。可以在受试者自由呼吸的情况下完成扫描，而无需闭气。同时也允许进一步扩大成像视野，提高空间分辨率。此外，上述成像方法100中，通过k-空间在各采样点间完全交错分段采集，从而实现原始图像的内在配准，无需在后期进行另外的图像处理。
图3示出了根据本发明一个实施例的成像序列。可以理解，在本发明的实施例中以循环的方式获得多个这样的成像序列。在每个成像序列中，共采集了3个图像信号。在各个心拍内，在心电门控信号和呼吸导航信号的控制下，执行不同的信号采集操作，以获得图像序列。图像信号的采集过程即磁共振成像中k-空间的填充过程。
如图3所示，在第一心拍内，在根据呼吸导航信号判断当前时刻符合预定条件的情况下，采集第一图像信号IMG1。第一图像信号IMG1是纵向磁化向量充分恢复到稳态的值。即第一图像信号IMG1为上述稳态图像信号。
该第一心拍内，未使用饱和脉冲。可以设第一图像信号IMG1对应延迟时间为无穷大的饱和脉冲，即Tsat1为无穷大。
可以理解，第一心拍内，根据呼吸导航信号判断当前时刻符合预定条件。在第一心拍之前，可能存在其间根据呼吸导航信号判断当前时刻不符合预定条件的心拍。因此，可选地，在信号采集操作中，在第一心拍之前还包括以下操作：在一个心拍内，在根据呼吸导航信号判断当前时刻不符合预定条件的情况下，等待下一个心拍，以再次根据呼吸导航信号执行判断操作并根据判断结果执行当前心拍的相应图像信号采集操作。为描述方便，称其间根据呼吸导航信号判断当前时刻不符合预定条件的心拍为A心拍。在A心拍中，不进行图像信号采集。在A心拍的下一个心拍中，再次根据呼吸导航信号判断当前时刻是否符合预定条件。如果仍然不符合，那么继续等待。直到在某一心拍中，根据呼吸导航信号判断当前时刻符合预定条件，则该心拍为第一心拍。如上所述，在该第一心拍内，当根据呼吸导航信号判断当前时刻符合预定条件时，采集第一图像信号IMG1。
在上述方案中，在第一心拍前的A心拍中，不采集图像信号。由此，能够保证在第一心拍内采集的磁化向量是在其平衡态的值。
在第二心拍内，首先使用延迟时间为Tsat2的饱和脉冲SAT。饱和脉冲可以将磁化向量置零。如图3所示，饱和脉冲的延迟时间是饱和脉冲到采集图像信号的时刻之间的时间间隔。在使用延迟时间为Tsat2的饱和脉冲SAT之后，在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第二图像信号IMG2。
与第二心拍类似的，在第三心拍内，在使用延迟时间为Tsat3的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第三图像信号IMG3。其中，Tsat3≠Tsat2。
在第二、第三心拍内，利用饱和脉冲实现了T1加权。其中，饱和脉冲的延迟时间不同，T1的权重不同。由此，在这两个心拍内，各获得了一个采样点。饱和脉冲的延迟时间可以是从系统允许的最小时间间隔到系统允许的最大时间间隔之间的任意值。
可选地，第二心拍内饱和脉冲的延迟时间Tsat2为系统允许的最大时间间隔Tmax的35％至70％。可以首先确定时间段Ttrigger期间的信号操作(例如呼吸导航信号NAV)所占用的时间长度与硬件响应延迟时间的和。然后计算时间段Ttrigger与该和的差，该差即系统允许的最大时间间隔Tmax。第三心拍内饱和脉冲的延迟时间Tsat3为系统允许的最大时间间隔Tmax的90％至100％。根据本发明一个实施例，Tsat2为Tmax/2，Tsat3等于Tmax。Tmax越大，磁化向量恢复的时间越长，也即可用于成像的图像信号越强，获得的图像信号的信噪比(SNR)越大，T1的权重越大。图4示出了根据该实施例的T1估计曲线。其中，横轴表示饱和脉冲的延迟时间，纵轴表示可用于数据读取的归一化的纵向磁化向量(Mz)，当Mz＝1时表示纵向磁化向量的稳态值。图4中还示出了分别在第一、第二和第三心拍内获得的采样点IMG1、IMG2和IMG3。Tsat2和Tsat3采用上述取值范围可以使采样点更合理的分布，从而使得在仅获得少量采样点的情况下也能够准确地估计T1值。此外，上述取值范围还使得可用于数据读取的纵向磁化向量比较大，从而提高了信号的信噪比，获得质量比较好的原始加权图像。
与第一心拍类似的，在第二和/或第三心拍之前，可能存在其间根据呼吸导航信号判断当前时刻不符合预定条件的心拍。可选地，在信号采集操作中，在第二和/或第三心拍之前还包括以下操作：在一个心拍内，在根据呼吸导航信号判断当前时刻不符合预定条件的情况下，采集图像信号并将所采集的图像信号设为无效，等待下一个心拍，以再次根据呼吸导航信号执行判断操作并根据判断结果执行当前心拍的相应图像信号采集操作。如果出现误判，所采集的图像信号可以用作原始数据。由此，保证了成像数据的完备性。
以第二心拍为例，假设在第二心拍之前，存在其间根据呼吸导航信号判断当前时刻不符合预定条件的心拍，简称为B心拍。在B心拍的下一个心拍中，再次根据呼吸导航信号判断当前时刻是否符合预定条件。如果仍然不符合，那么采集图像信号并将所采集的图像信号设为无效，继续等待下一个心拍。直到在某一心拍中，根据呼吸导航信号判断当前时刻符合预定条件，则该心拍为第二心拍。如上所述，在该第二心拍内，在使用延迟时间为Tsat2的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第二图像信号IMG2。
上述信号采集操作中，采集了T1权重不同的图像信号(IMG1、IMG2和IMG3)。图5示出了根据本发明一个实施例的所采集的受试者的短轴视角下的心脏的IMG1、IMG2和IMG3(对于IMG1、IMG2和IMG3中的每个，都从心尖、心中、心基分别选一层作为代表图像)以及根据这些图像拟合得到的参数T1定量图像。图像信号IMG1、IMG2和IMG3循环采集。即在完成第三心拍的IMG3采集之后，再返回第一心拍。然后，重复上述过程。换言之，在成像过程中，重复采集上述成像序列，直至完成磁共振成像中k-空间的所有分段的填充。可以理解，上述第一心拍至第三心拍的顺序仅为示例，而非对本发明的限制。可以以任意顺序执行这3个心拍，而不影响本申请技术方案的效果。
如前所述，在第一心拍中采集的是稳态图像信号。为了使得纵向磁化矢量能够从上一次图像信号采集(IMG3)之后完全恢复，在第一心拍前，设置恢复时间段。该恢复时间段中，不采集图像信号。可选地，在该恢复时间段中，只采集呼吸导航信号NAV，以保证呼吸导航信号的连续性。由此，避免干扰呼吸导航信号NAV的参数设定以及它所提供的灵活性。进而保证呼吸导航信号NAV准确地控制信号采集操作，以获得更准确的图像信号。
通过上述优化的成像序列，上述成像方法100扫描效率高、成像分辨率不受限制。由此，只需要较少的采样点，例如仅上述IMG1、IMG2和IMG3这3个采样点，就能够实现3D心肌定量参数T1准确测量，而且无需额外的滤波数据处理。减少采样点的个数，直接的技术效果是缩短成像过程中的扫描时间。
可以理解，上述采集图像信号(IMG1、IMG2和IMG3)可以利用各种合适的数据读取方式。该数据读取方式包括但不限于梯度回波(gradient echo)、回波平面成像(Echo-PlanarImaging，EPI)和自旋回波(Spin Echo)等等。优选地，数据读取方式采用损毁梯度回波(spoiled gradient echo，SPGR)、平衡稳态自由进动(balanced Steady state freeprecession，bSSFP)和梯度自旋回波(Gradient Spin Echo，Grase)技术。这些优选数据读取方式的采用，能够显著降低成像过程对磁场强度均匀性的要求，使得本方案可应用于高场(如3T)磁共振系统。
图6示出了根据本发明另一个实施例的成像序列。图6所示成像序列与图3所示成像序列类似。为了简洁，对两个成像序列中的相同部分不再赘述。如图6所示，在信号采集操作中的采集图像信号(IMG1、IMG2和IMG3)之前，可以分别执行压脂操作(FS)。压脂操作有助于降低呼吸伪影，显著提高成像质量。
根据本发明一个实施例，可以根据上述第i图像信号以及第i图像信号的饱和脉冲延迟时间Tsati确定参数T1，其中，i＝1、2、3。并且当i＝1时，用于数据拟合的Tsat1为无穷大。
在一个示例中，根据如下公式确定参数T1，其中，i＝1、2、3。Si和Tsati分别是信号采集操作所获得的第i图像信号和其对应的饱和脉冲的延迟时间。S0是磁化向量在平衡态时的理论图像信号。S0与T1在此公式中是未知的。根据IMG1、IMG2和IMG3能够确定这二者。根据该公式能够更准确地确定参数T1，从而生成更准确的图像。
图7a和图7b分别示出了根据本发明一个实施例的来自同一个受试者的3D左心室的T1图像和该T1图像在长轴视角下的图像。图7c示出了图7a左心室全部心肌的T1的统计直方图。
图8示出了根据本发明一个实施例的、来自两个受试者的参数T1定量的图像。这些图像采集于这两个受试者使用T1缩短对比剂Gd-DTPA(0.15mmol/kg,Magnevist,BayerPharma AG,Germany)后15分钟。这些图像分别选自心尖、心中、心基。在该实施例中，由于测量心肌的T1值在对比剂的作用下明显缩短，磁化向量的恢复过程也比较快，因此，稳态图像信号(IMG1)采集前的恢复时间段可以设置为3秒。
图2、图5、图7和图8的数据均来自身体健康的受试者。如图7c所示的直方图所示，根据本发明所获得的参数T1的数值分布呈现正态分布。而且参数T1的数值的标准方差较小。因此，根据本发明的实施例所生成的心肌定量磁共振图像较理想地反映了受试者的心肌组织状态。
可选地，上述信号采集操作还包括至少一个以下操作：在第f1心拍内，在使用延迟时间为Tsatf1的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第f1图像信号，其中，Tsatf1不等于其他图像信号对应的饱和脉冲的延迟时间，f1为不等于1、2和3的整数。例如，可以在第四心拍和第五心拍内各采集一个图像信号。该操作与上述采集IMG2和IMG3的操作类似，为了简洁，在此不再赘述。通过该操作，增加了采样点。
可以理解，上述确定参数T1还根据第f1图像信号和Tsatf1。由此可以通过更多采样点参与拟合而获得更准确的参数T1。
可选地，所述信号采集操作还包括至少一个以下操作：在第f2心拍内，在使用延迟时间为Tsat2或Tsat3的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下再次采集与饱和脉冲对应的第二图像信号或第三图像信号，f2为不等于1、2和3的整数。可以理解，如果成像方法中还采集了第f1图像信号，则f2也不等于f1。该操作是重复第二心拍或第三心拍的操作，由此获得了T1权重相同的采样点。在一个示例中，在第六心拍内，在使用延迟时间为Tsat2的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下再次采集与饱和脉冲相应的第二图像信号。在第七心拍内，在使用延迟时间为Tsat3的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下再次采集与饱和脉冲相应的第三图像信号。最后，根据所采集的所有第二图像信号和/或所有第三图像信号确定参数T1。在上述示例中，将第二心拍、第六心拍内采集的第二图像信号和第三心拍、第七心拍内采集的第三图像信号全部输入信号模型进行拟合，以确定参数T1。
上述技术方案的效果相当于平均了采样点(例如第二图像信号)的噪声，从而降低了拟合偏差。总之，上述信号采集操作可以提高参数T1的计算准确率，从而提高图像质量。
根据本发明又一方面，还提供了一种用于心肌定量磁共振成像的设备。该系统包括处理器和存储器。所述存储器存储用于实现根据本发明实施例的心肌定量磁共振成像的方法中的各个步骤的计算机程序指令。所述处理器用于运行所述存储器中存储的计算机程序指令，以执行根据本发明实施例的心肌定量磁共振成像的方法的相应步骤。
根据本发明再一方面，还提供了一种存储介质，在所述存储介质上存储了程序指令，在所述程序指令被计算机或处理器运行时使得所述计算机或处理器执行本发明实施例的心肌定量磁共振成像的方法的相应步骤，并且用于实现根据本发明实施例的用于心肌定量磁共振成像的装置中的相应模块。所述存储介质例如可以包括平板电脑的存储部件、个人计算机的硬盘、只读存储器(ROM)、可擦除可编程只读存储器(EPROM)、便携式紧致盘只读存储器(CD-ROM)、USB存储器、或者上述存储介质的任意组合。所述计算机可读存储介质可以是一个或多个计算机可读存储介质的任意组合。
在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
类似地，应当理解，为了精简本发明并帮助理解各个发明方面中的一个或多个，在对本发明的示例性实施例的描述中，本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该本发明的方法解释成反映如下意图：即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如相应的权利要求书所反映的那样，其发明点在于可以用少于某个公开的单个实施例的所有特征的特征来解决相应的技术问题。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。
本领域的技术人员可以理解，除了特征之间相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
应该注意的是单词“包含”不排除存在未列在权利要求中的元件或步骤。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
以上所述，仅为本发明的具体实施方式或对具体实施方式的说明，本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。本发明的保护范围应以权利要求的保护范围为准。

Claims (12)

1.一种心肌定量磁共振成像方法，包括：
每隔恢复时间段，在心电门控信号和呼吸导航信号的控制下，至少执行如下信号采集操作：
在第一心拍内，在根据呼吸导航信号判断当前时刻符合预定条件的情况下，采集第一图像信号；
在第二心拍内，在使用延迟时间为Tsat2的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第二图像信号；
在第三心拍内，在使用延迟时间为Tsat3的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第三图像信号，其中，Tsat3≠Tsat2；
根据第i图像信号以及第i图像信号对应的饱和脉冲的延迟时间Tsati确定参数T1，其中，i＝1、2、3，且当i＝1时，Tsati为无穷大；
根据所述参数T1生成心肌定量磁共振图像。
2.如权利要求1所述的方法，其中，所述信号采集操作还包括至少一个以下操作：
在第f1心拍内，在使用延迟时间为Tsatf1的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下采集第f1图像信号，其中，Tsatf1不等于其他图像信号对应的饱和脉冲的延迟时间，f1为不等于1、2和3的整数；
其中，所述确定参数T1还根据所述第f1图像信号和所述Tsatf1。
3.如权利要求1所述的方法，其中，所述信号采集操作还包括至少一个以下操作：
在第f2心拍内，在使用延迟时间为Tsat2或Tsat3的饱和脉冲之后，并且在根据呼吸导航信号判断当前时刻符合预定条件的情况下再次采集与饱和脉冲对应的第二图像信号或第三图像信号，f2为不等于1、2和3整数。
4.如权利要求1所述的方法，其中，在所述信号采集操作中，在所述第一心拍之前，还包括以下操作：
在一个心拍内，当根据呼吸导航信号判断当前时刻不符合预定条件时，等待下一个心拍，以再次根据呼吸导航信号执行判断操作并根据判断结果执行当前心拍的相应图像信号采集操作。
5.如权利要求1或4所述的方法，其中，在所述信号采集操作中，在所述第二心拍和/或所述第三心拍之前，还包括以下操作：
在一个心拍内，在根据呼吸导航信号判断当前时刻不符合预定条件的情况下，采集图像信号并将所采集的图像信号设为无效，等待下一个心拍，以再次根据呼吸导航信号执行判断操作并根据判断结果执行当前心拍的相应图像信号采集操作。
6.如权利要求1所述的方法，其中，所述采集第一图像信号、所述采集第二图像信号和所述采集第三图像信号利用损毁梯度回波、平衡稳态自由进动或梯度自旋回波的数据读取方式。
7.如权利要求1所述的方法，其中，在所述信号采集操作中，在所述采集第一图像信号、所述采集第二图像信号和所述采集第三图像信号之前，分别执行压脂操作。
8.如权利要求1所述的方法，其中，
所述Tsat3为系统允许的最大时间间隔Tmax的90％至100％，所述Tsat2为所述Tmax的35％至70％。
9.如权利要求1所述的方法，其中，所述恢复时间段包括n个心拍，其中n是大于或等于N/(60/HR)的最小整数，N为允许磁化向量恢复的时间长度，单位为秒，HR为心率，单位为心跳/分钟。
10.如权利要求1所述的方法，其中，所述确定参数T1包括根据如下公式确定所述参数T1，
其中，S0是纵向磁化向量在平衡态时的信号，Si和Tsati分别为第i图像信号和其对应的饱和脉冲的延迟时间，i＝1、2、3。
11.一种用于心肌定量磁共振成像的设备，包括处理器和存储器，其中，所述存储器中存储有计算机程序指令，所述计算机程序指令被所述处理器运行时用于执行如权利要求1至10任一项所述的心肌定量磁共振成像方法。
12.一种存储介质，在所述存储介质上存储了程序指令，所述程序指令在运行时用于执行如权利要求1至10任一项所述的心肌定量磁共振成像方法。
CN201810923440.2A 2018-08-14 2018-08-14 心肌定量磁共振成像方法、设备及存储介质 Active CN109091145B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201810923440.2A CN109091145B (zh) 2018-08-14 2018-08-14 心肌定量磁共振成像方法、设备及存储介质
PCT/CN2019/086400 WO2020034676A1 (zh) 2018-08-14 2019-05-10 心肌定量磁共振成像方法、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810923440.2A CN109091145B (zh) 2018-08-14 2018-08-14 心肌定量磁共振成像方法、设备及存储介质

Publications (2)

Publication Number Publication Date
CN109091145A true CN109091145A (zh) 2018-12-28
CN109091145B CN109091145B (zh) 2021-08-10

Family

ID=64849622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810923440.2A Active CN109091145B (zh) 2018-08-14 2018-08-14 心肌定量磁共振成像方法、设备及存储介质

Country Status (2)

Country Link
CN (1) CN109091145B (zh)
WO (1) WO2020034676A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020034676A1 (zh) * 2018-08-14 2020-02-20 清华大学 心肌定量磁共振成像方法、设备及存储介质
CN112494072A (zh) * 2019-09-16 2021-03-16 美国西门子医疗系统股份有限公司 定量医学诊断超声的肌肉收缩状态触发

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7689263B1 (en) * 2005-09-15 2010-03-30 General Electric Company Method and apparatus for acquiring free-breathing MR images using navigator echo with saturation RF pulse
CN102143707A (zh) * 2008-09-04 2011-08-03 株式会社东芝 磁共振成像装置
CN103376432A (zh) * 2012-04-20 2013-10-30 西门子公司 快速地位置分辨地确定检查区域内磁共振弛豫参数的方法
US20150123659A1 (en) * 2013-11-07 2015-05-07 Sebastian Weingartner System and method for free-breathing volumetric imaging of cardiac tissue
CN106233154A (zh) * 2014-04-22 2016-12-14 皇家飞利浦有限公司 使用预脉冲和导航器的具有运动校正的磁共振成像
US20170325707A1 (en) * 2016-05-16 2017-11-16 Siemens Healthcare Gmbh Free breathing motion corrected pixel-wise mri myocardial t1 parameter mapping

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006158512A (ja) * 2004-12-03 2006-06-22 Hitachi Medical Corp 磁気共鳴イメージング装置
JP5366370B2 (ja) * 2006-09-06 2013-12-11 株式会社東芝 磁気共鳴イメージング装置
CN102018511A (zh) * 2009-09-18 2011-04-20 株式会社东芝 磁共振成像装置以及磁共振成像方法
CN103479356B (zh) * 2012-06-14 2015-07-29 西门子(深圳)磁共振有限公司 一种扩散张量磁共振成像方法
CN104013404B (zh) * 2014-06-09 2016-04-20 深圳先进技术研究院 心肌t1值测量方法和系统
CN109091145B (zh) * 2018-08-14 2021-08-10 清华大学 心肌定量磁共振成像方法、设备及存储介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7689263B1 (en) * 2005-09-15 2010-03-30 General Electric Company Method and apparatus for acquiring free-breathing MR images using navigator echo with saturation RF pulse
CN102143707A (zh) * 2008-09-04 2011-08-03 株式会社东芝 磁共振成像装置
CN103376432A (zh) * 2012-04-20 2013-10-30 西门子公司 快速地位置分辨地确定检查区域内磁共振弛豫参数的方法
US20150123659A1 (en) * 2013-11-07 2015-05-07 Sebastian Weingartner System and method for free-breathing volumetric imaging of cardiac tissue
CN106233154A (zh) * 2014-04-22 2016-12-14 皇家飞利浦有限公司 使用预脉冲和导航器的具有运动校正的磁共振成像
US20170325707A1 (en) * 2016-05-16 2017-11-16 Siemens Healthcare Gmbh Free breathing motion corrected pixel-wise mri myocardial t1 parameter mapping

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RUI GUO ET AL.: "A three-dimensional free-breathing sequence for simultaneous myocardial T1 and T2 mapping", 《MAGNETIC RESONANCE IN MEDICINE》 *
RUI GUO ET AL.: "Three-dimensional free breathing whole heart cardiovascular magnetic resonance T1 mapping at 3T", 《JOURNAL OF CARDIOVASCULAR MAGNETIC RESONANCE》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020034676A1 (zh) * 2018-08-14 2020-02-20 清华大学 心肌定量磁共振成像方法、设备及存储介质
CN112494072A (zh) * 2019-09-16 2021-03-16 美国西门子医疗系统股份有限公司 定量医学诊断超声的肌肉收缩状态触发

Also Published As

Publication number Publication date
WO2020034676A1 (zh) 2020-02-20
CN109091145B (zh) 2021-08-10

Similar Documents

Publication Publication Date Title
Scott et al. Motion in cardiovascular MR imaging
KR101664138B1 (ko) 멀티 슬라이스 자기 공명 영상을 사용하는 조직 특성 감별 시스템 및 방법
CN103376432B (zh) 快速地位置分辨地确定检查区域内磁共振弛豫参数的方法
CN110031786B (zh) 磁共振图像重建方法、磁共振成像方法、设备及介质
CN106539584B (zh) 磁共振成像方法及系统
US20120108946A1 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
KR101809213B1 (ko) Ｔ2-이완 시간의 다중 슬라이스 영상화를 위한 방법 및 장치
US8154288B2 (en) Method, processor, and magnetic resonance apparatus for selective presentation of lung movement
Becker et al. Fast myocardial T1 mapping using cardiac motion correction
EP2886053A1 (en) Magnetic resonance imaging device
US20180031662A1 (en) System and Method for Mapping and Quantifying In-Vivo Blood Flow Stasis
CN109091145A (zh) 心肌定量磁共振成像方法、设备及存储介质
CN109247938B (zh) 定量心肌磁共振成像方法、设备及存储介质
CN109077728B (zh) 定量心肌磁共振成像方法、设备及存储介质
Stalder et al. Robust cardiac BOLD MRI using an fMRI‐like approach with repeated stress paradigms
CN109998548B (zh) 定量心肌磁共振成像方法、设备及存储介质
JP7337603B2 (ja) 磁気共鳴イメージング装置
US20230270349A1 (en) Magnetic resonance imaging
Duerinckx Cardiac MRI for clinicians: an overview
Guo et al. MyoFold: Joint myocardial tissue composition and wall motion quantification via a highly folded sequence
Abdellatif et al. Basic Principles, Technique and Planning for Cardiac Magnetic Resonance Imaging
CN116350205A (zh) 用于定量心肌磁共振成像的方法、设备及存储介质
Kawanami et al. Flow-weighted MRI of the Lungs with the ECG-gated half-Fourier FSE technique: Evaluation of the Effect of the Cardiac Cycle
Burger A control system approach to subject specific prospective respiratory motion correction in cardiac MRI
Wang Quantitative Evaluation of Left Ventricular Rotational Mechanics using Magnetic Resonance Imaging

Legal Events

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