CN110547799B - 磁共振成像方法、计算机设备和计算机可读存储介质 - Google Patents

磁共振成像方法、计算机设备和计算机可读存储介质 Download PDF

Info

Publication number
CN110547799B
CN110547799B CN201910762725.7A CN201910762725A CN110547799B CN 110547799 B CN110547799 B CN 110547799B CN 201910762725 A CN201910762725 A CN 201910762725A CN 110547799 B CN110547799 B CN 110547799B
Authority
CN
China
Prior art keywords
data
acquisition
data set
region
imaging sequence
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
CN201910762725.7A
Other languages
English (en)
Other versions
CN110547799A (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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201910762725.7A priority Critical patent/CN110547799B/zh
Publication of CN110547799A publication Critical patent/CN110547799A/zh
Application granted granted Critical
Publication of CN110547799B publication Critical patent/CN110547799B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本申请涉及一种磁共振成像方法、计算机设备和计算机可读存储介质,所述方法包括:在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,激发第一成像序列,获得目标区域的第一数据集;第一数据集包括欠采样区域、全采样区域;第一成像序列包括第一采样模块;在第二采集周期内激发第二成像序列,获得目标区域的第二数据集;第二数据集包括欠采样区域;第二成像序列包括准备模块和第二采样模块;根据第一数据集中的全采样区域的数据恢复第一数据集或/和第二数据集中的欠采样区域的数据,以获得校正的数据集,并获得目标区域的磁共振图像。上述方法中,在第一成像序列内采集的全采样区域数据(即参考线对应的数据)的信噪比较高。

Description

磁共振成像方法、计算机设备和计算机可读存储介质
技术领域
本申请涉及图像技术领域,特别是涉及一种磁共振成像方法、计算机设备和计算机可读存储介质。
背景技术
磁共振成像(Magnetic Resonace Imaging,MRI)是断层成像的一种,通过对静磁场中的人体施加某种特定频率的射频脉冲,使人体中的氢质子受到激励而发生磁共振现象。停止脉冲后,质子在弛豫过程中,通过对MR信号的接收、空间编码和图像重建等处理过程,产生MR信号。
为了缩短采集数据所需要的时间,通常采用并行加速的方法实现数据的采集;然而现有技术中通过采用并行加速的方式所采集的数据中,采集的并行加速线的信噪比较低。
发明内容
基于此,有必要针对上述技术问题,提供一种磁共振成像方法、计算机设备和计算机可读存储介质。
一种磁共振成像方法,所述方法包括:
在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,利用第一成像序列,获得所述目标区域的第一数据集;所述第一数据集包括欠采样区域和全采样区域;所述第一成像序列包括第一采样模块;
在所述一次数据采集过程中的第二采集周期内,利用第二成像序列,获得所述目标区域的第二数据集;所述第二数据集包括欠采样区域;所述第二成像序列包括准备模块和第二采样模块;
根据所述第一数据集中的全采样区域的数据恢复所述第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集;
根据所述校正的数据集,获得所述目标区域的磁共振图像。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述方法的步骤。
上述磁共振成像方法、计算机设备和计算机可读存储介质,在对检测对象的目标区域的一次采集过程中,第一采集周期内利用第一成像序列采集获得第一数据集,而在第二采集周期内,利用第二成像序列以获得第二数据集;其中,第一成像序列中不包含准备模块,第二成像序列中包含准备模块;第一数据集包括欠采样区域和全采样区域,第二数据集包括欠采样区域;然后通过第一数据集中的全采样区域的数据恢复所有欠采样区域的数据,获得校正的数据集,根据校正的数据集获得磁共振图像。如此,由于第一成像序列中不包含准备模块,在期间采集的全采样区域数据(即参考线对应的数据)的信噪比较高。
附图说明
图1为一个实施例中磁共振成像方法的流程示意图;
图2为一个具体实施例中首过灌注的第一采样周期的序列结构示意图;
图3为一个具体实施例中首过灌注的第二采样周期的序列结构示意图;
图4为另一个实施例中具体实施例中心脏延迟增强序列PSIR中的序列结构示意图;
图5为一个具体实施例中T2Mapping中的序列结构示意图;
图6为一个实施例中磁共振成像装置的结构框图;
图7为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
在一个实施例中,如图1所示,提供了一种磁共振成像方法,磁共振成像是断层成像的一种,它利用磁共振现象从人体中获得电磁信号,并重建出人体信息。在本实施例中磁共振成像方法,包括步骤S110至步骤S140。
步骤S110,在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,利用第一成像序列,获得目标区域的第一数据集;第一数据集包括欠采样区域和全采样区域;第一成像序列包括第一采样模块。
在本实施例中,将有待进行磁共振成像的对象记为检测对象,检测对象可以是健康试者,也可以是患者,或者还可以是动物。在一个实施例中,目标区域可以是任何部位或组织,如心脏、血管或其他存在搏动区域的器官或者组织。在利用序列激发目标区域之前,可向检测对象注入对比剂,该对比剂可影响目标区域的成像对比度。激发序列的类型可以时成像序列或者饱和脉冲序列。可选地,利用成像序列可激发目标区域,以获取磁共振信号,该磁共振信号经过相位编码、频率编码等,可得到K空间数据集。
一次数据采集过程对应一次磁共振成像的过程。在一次数据采集过程中,包括多个采集周期,采集周期与采集周期之间所采用的成像序列可以是相同的,也可以是不相同的,在本实施例中,将激发第一成像序列的采集周期记为第一周期。将激发第一成像序列所获得的数据记为第一数据集,并且在本实施例中,第一数据集包括部分欠采样区域和部分全采样区域,欠采样区域表示这一数据区域中包含未完全采集的数据,而相对应地,全采样区域则表示这一数据区域中包含采集完整的数据。通过全采样区域的数据计算获得欠采样区域中未采集完全的数据,从而可以获得所有数据对应的磁共振图像;这一方法可称为并行加速方法,第一数据集中的全采样区域对应并行加速参考线,欠采样区域对应非参考线。
可选地,全采样区域以奈奎斯特速率采样;而欠采样区域则每采集一条相位编码线或成像数据线后可连续跳过R–1条相位编码步。由于填充K空间不同区域的回波信号所采用的相位编码梯度场强度不同,因此其信号强度也存在差别。越靠近K空间中心的相位编码线所施加的相位编码梯度场越弱,其信号强度越高,因此对图像的对比影响越大,但缺乏空间信息;越靠近K空间周边的相位编码线,所使用的相位编码梯度场强度越强,所提供相位编码方向的空间信息越丰富,但磁共振信号的幅度越小,对图像的对比贡献越小。在一个实施例中,全采样区域为K空间中的部分区域,且在全采样区域的全部K空间填充位点填充K空间数据,而欠采样区域分为位于欠采样区域的两侧。在另一实施例中,全采样区域和欠采样区域可交错/间隔分布。
其中,第一成像序列包括第一采样模块,根据第一数据集中包括部分欠采样区域和部分全采样区域可以确定,第一采样模块包括欠采样序列和全采样序列;欠采样序列可以非奎斯特(Nyquist)频率进行稀疏采样以对第一数据集的欠采样区域进行数据填充;全采样序列可以奎斯特(Nyquist)频率进行顺序采样以对第一数据集的全采样区域进行数据填充。在上述实施例中,第一成像序列中的第一采样模块前未施加其他准备模块,如此,能够提高第一采集周期中全采样区域中的数据的信噪比。
步骤S120,在一次数据采集过程中的第二采集周期内,利用第二成像序列,获得目标区域的第二数据集;第二数据集包括欠采样区域;第二成像序列包括准备模块和第二采样模块。
在一个实施例中,第二成像序列中包含脉冲序列模块,例如包括饱和恢复脉冲、反转恢复脉冲或者T2准备脉冲;在一个具体实施例中,准备模块包括饱和恢复脉冲,在一个第二采集周期内施加一饱和恢复脉冲,然后使用一成像序列。成像序列可以是自由感应衰减(free induction decay,FID)类序列、自旋回波类序列、梯度回波类序列或者自旋回波与梯度回波组成的杂合序列、快速自旋梯度回波序列与平面回波序列组成的杂合序列等。
在与步骤S110相同的一次数据采集过程中,将利用第二成像序列激发目标区域的采集周期记为第二采集周期,将第二采集周期期间采集的数据记为第二数据集,在本实施例中,第二数据集只包含欠采样区域。
在一个实施例中,目标区域中包含心脏,即磁共振成像是对心脏进行检测,在本实施例中,一个采集周期即为一个心动周期。其中,心动周期是指从一次心跳的起始到下一次心跳的起始,心血管系统所经历的过程。在一个心动周期内使用第一成像序列采集数据获得第一数据集,在另一个心动周期内使用第二成像序列采集数据获得第二数据集。
其中,第二成像序列包括第二采样模块,根据第二数据集中仅包括欠采样区域可以确定,第二采样模块仅包括欠采样部分;在并行加速方法中,欠采样部分为非参考线对应的数据。在上述实施例中,第二成像序列中的第二采样模块采用加速采集/稀疏采集的方式获取磁共振信号,如此,实际采集过程中磁共振信号的K空间数据线减少,能够提高数据采集效率,节省采集时间。
在一个实施例中,第二成像序列的激发过程可以采用包括基于图像的SENSE技术、基于K空间的SMASH(simultaneous acquisition of spatial harmonics)或者GRAPPA(generalized autocalibrating partially acquisition)技术,还可以包括基于混合K空间(Hybrid space)的ARC(auto-calibrating reconstruction for cartesian sampling)技术等并行加速采集序列。上述并行加速采集序列利用多通道表面线圈的空间位置差异,来区分信号的位置来源,以此通过每个线圈彼此之间的配合,减少成像过程中相位编码的步数。
进一步地,在一个实施例中,第一采样模块的施加起始时间与第二采样模块的施加起始时间一致。即在第一采集周期与第二采集周期中所施加的采样模块对应不同心动周期内的同一期相,如两采样模块同时对应心脏的舒张期;也即,第一采样模块与第二采样模块所对应的成像序列在检测对象生理运动的同一期相内激发目标区域。如此,可以利用第一采样模块中采集的参考线来恢复第二采样模块中采集的非参考线中未采集的数据,从而获得本次数据采集过程中的所有数据。
步骤S130,根据第一数据集中的全采样区域的数据恢复第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集。
在一个实施例中,根据第一数据集中的全采样区域的数据恢复第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集,包括:通过第一数据集中的全采样区域的数据获取恢复系数;根据恢复系数,恢复第一数据集和/或第二数据集的欠采样区域的数据,获得校正的数据集。
进一步地,在一个实施例中,根据第一数据集中的全采样区域的数据恢复第二数据集中的欠采样区域的数据,包括:根据第一数据集中全采样区域的数据获取恢复系数,根据恢复系数对第二数据集中欠采样区域的数据进行恢复,以获取填充完整的第二数据集。
在另一个实施例中,根据第一数据集中的全采样区域的数据恢复第二数据集中的欠采样区域的数据,包括:根据第一数据集中全采样区域的数据获取部分校准系数,根据该部分校准系数可对第一数据集中欠采样区域的数据进行恢复,获取填充完整的第一数据集。进一步地,根据填充完整的第一数据集,获取完整校准系数,并根据完整校准系数对第二数据集中欠采样区域的数据进行恢复,以获得完整填充的第二数据集。
进一步地,根据填充完整的第一数据集,获取完整校准系数,并根据完整校准系数对第二数据集中欠采样区域的数据进行恢复,计算获得校正的数据集。需要说明的是,本申请实施例中的“恢复系数”也可用“校准系数”代替,表征K空间中相邻编码线之间的关系,通过该系数可计算缺失的数据,该系数与接收线圈分布有关。
在一个实施例中,校准系数可通过如下方式获取:
在第一数据集中全采样区域的任意一个数据点PA,以PA作为校准数据点,建立PA分别与第一数据集中全采样区域相邻数据点之间的关系,并基于该关系生成滤波器。同样的,在第一数据集中全采样区域内,可构建全采样区域的中心处的数据点与其余相邻各数据点之间的关系。仅作为示例,滤波器可以是基于全采样区域的中心处的数据点和其他数据点来生成的卷积核。
在本实施例中,将根据全采样区域的数据进行恢复后获得的数据记为校正的数据集;在一个实施例中,校正的数据集包括填充完整的第二数据集;在另一个实施例中,校正的数据集包括填充完整的第一数据集以及填充完整的第二数据集。
步骤S140,根据校正的数据集,获得目标区域的磁共振图像。
在根据全采样区域中的数据采用并行加速的方法获得一次数据采集过程中的所有数据的数据集之后,根据这一数据集即可获得目标区域的磁共振图像。
进一步地,根据校正的数据集,获得目标区域的磁共振图像可以是通过任意一种方式实现。
上述磁共振成像方法,在对检测对象的目标区域的一次采集过程中,第一采集周期内利用第一成像序列激发目标区域,获得第一数据集,而在第二采集周期内,利用第二成像序列激发目标区域以获得第二数据集;其中,第一成像序列中不包含准备模块,第二成像序列中包含准备模块;第一数据集包括欠采样区域和全采样区域,第二数据集包括欠采样区域;然后通过第一数据集中的全采样区域的数据恢复所有欠采样区域的数据,获得校正的数据集,根据校正的数据集获得磁共振图像。如此,由于第一成像序列中不包含准备模块,在期间采集的全采样区域数据(即参考线对应的数据)的信噪比较高。
在一个实施例中,第一采集周期包括一次数据采集过程中前预设数目个采集周期,第二采集周期包括第一采集周期以外的采集周期;准备模块包括饱和恢复模块。饱和恢复模块对应的为饱和恢复脉冲。
其中,饱和恢复(saturation recovery,SR)序列也称部分饱和(partialsaturation)序列。SR序列利用连续的90°脉冲进行激发,在每个90°脉冲后采集FID信号。可以理解地,第二成像序列中的采样模块为欠采样模块,对应欠采样的脉冲序列。
进一步地,在本实施例中,在一次数据采集过程内,第一成像序列和第二成像序列均包括四个成像片层参数。其中,成像片层是指待扫描组织的分层,四个成像片层参数指四个成像片层同时射频激发或者连续射频激发的参数设置。磁共振成像是多切面的断层显像。要使某一待扫描的人体组织分层面显示,就要进行层面定位,人为地分解组织器官成为许多具有一定厚的断面。在一个实施例中,该四个成像片层参数相关的脉冲可同时激发。
在上述第一成像序列和第二成像序列均包括四个成像片层参数的实施例中,可以是心脏首过灌注的数据采集过程;心脏首过灌注是用来观察对比剂逐渐进入心脏的一种成像技术,序列的结构为心电触发的饱和恢复脉冲加上梯度回波序列,饱和恢复脉冲饱和组织信号,随着时间推移,对比剂逐渐进入心脏和心肌,图像信噪比逐渐增强。心脏首过灌注通常在一个心动周期需要采集4幅图像(其中包括3幅短轴图像和1幅4腔图像;对应4个成像片层),重复采集60个心动周期左右。为了在一个心动周期内完成4幅图像的采集,一般需要使用并行加速技术。而首过灌注扫描时,由于有饱和恢复脉冲的作用,在对比剂未进入心脏时,采集的图像信噪比低,采集的并行加速参考线的信噪比也同样低。
在本实施例中,在一次数据采集过程中,将前面的预设数目个采集周期定义为第一采集周期,将该次数据采集过程中其它的采集周期定义为第二采集周期;即,在本次数据采集过程中,前预设数目个采集周期中激发第一成像序列,采集的数据包括全采样区域和欠采样区域;而在本次数据采集过程中的剩余其它采集周期中,激发第二成像序列,采集的数据仅包括欠采样区域。在另一个实施例中,也可以称为在本次数据采集过程中的前预设数目个采集周期中,采集参考线和非参考线,参考线填充在K空间的全采样区域,非参考线填充在K空间的欠采样区域;而在本次数据采集过程中的后面的采集周期中,仅采集非参考线,即在本次数据采集过程中的后面的采集周期采用欠采样。其中,预设数目可以根据实际情况进行设置。
在上述磁共振成像方法中,在一次数据采集过程中的前预设数目个采集周期中采集获得全采样区域的数据,从而在后面的采集周期中激发第二成像序列采集获得的欠采样区域的数据时,可以在数据采集的同时,利用前面采集获得的全采样区域的数据来恢复欠采样区域中未采集部分的数据,即数据采集完成后即可获得对应的磁共振图像,如此,可以提高磁共振成像的效率。
如图2所示,为一个具体实施例中首过灌注的第一采样周期的序列结构示意图;如图3所示,为本实施例中首过灌注的第二采样周期的序列结构示意图。其中,S#(S1、S2、S3和S4)表示成像片层,SR(Saturation-Recovery)表示饱和恢复脉冲。在本实施例中,首过灌注的前预设数目个心动周期内采集不施加饱和恢复脉冲的数据,直接使用第一采样模块,其中包含并行加速的参考线和非参考线;而后面的心动周期内,先施加饱和恢复脉冲,然后使用第二采样模块,仅采集非参考线。
然后将前预设数目个心动周期中采集的数据中并行加速的参考线(第一采集周期中全采样区域的数据)作为后续施加饱和恢复脉冲的数据的参考线,恢复其中的非参考线对应的数据(第二采集周期中欠采样区域的数据),以获得本次首过灌注中的所有数据,进一步可以得到磁共振图像。如此,由于前预设数目个心动周期中未施加饱和恢复脉冲,其中采集的参考线的信噪比较高,使得并行加速的计算结果更好;且后续心动周期中无需再采集参考线,可以节省后续每帧的采集时间,从而可用于心率过快的病人。
进一步地,在另一个实施例中,准备模块包括饱和恢复模块;在本实施例中,第二采集周期包括一次数据采集过程中前预设数目个采集周期,第一采集周期包括第二采集周期以外的采集周期,在时间上第二采集周期之后/后面的采集周期。饱和恢复模块对应的为饱和恢复脉冲。
在本实施例中,在一次数据采集过程中,将前面的预设数目个采集周期定义为第二采集周期,将该次数据采集过程中其它的采集周期定义为第一采集周期;即,在本次数据采集过程中,前预设数目个采集周期中激发第二成像序列,采集的数据仅包括欠采样区域;而在本次数据采集过程中的剩余其它采集周期中,激发第一成像序列,采集的数据包括全采样区域和欠采样区域。在另一个实施例中,也可以称为在本次数据采集过程中的前预设数目个采集周期中,仅采集非参考线;而在本次数据采集过程中的后面的采集周期中,采集参考线和非参考线。其中,预设数目可以根据实际情况进行设置。
在上述磁共振成像方法中,是在一次数据采集过程中,在前预设数目个采集周期中仅采集非参考线,而在该次数据采集过程中的剩余采集周期中,采集参考线和非参考线;可以理解地,在本实施例中,在后面的第一采集周期内完成采集获得全采样区域的数据之后,可以计算恢复前预设数目个第二采集周期中欠采样区域的数据。由于前预设数目个第二采集周期中仅需要采集欠采样区域的数据,从整体数据采集过程来说,可以提高数据采集的效率,节省采集时间。
在一个实施例中,准备模块包括反转恢复模块;进一步地,在本实施例中,第二采集周期包括一次数据采集过程中的第一个采集周期;第一采集周期包括一次数据采集过程中的第二个采集周期。
进一步地,本实施例中表示心脏延迟增强序列PSIR(相位敏感反转恢复)的数据采集过程。反转恢复模块对应的为反转恢复序列。
反转恢复(inversion recovery,IR)序列的特点是利用180°射频脉冲把组织的宏观纵向磁化矢量偏转180°,即反转到与主磁场相反的方向上,在组织发生纵向弛豫的过程中施加90°脉冲,来记录不同组织间纵向弛豫的差别。90°脉冲后可以采集FID信号,也可以利用180°聚焦脉冲采集自旋回波信号。可以理解地,第二成像序列中的采样模块为欠采样模块,对应欠采样的脉冲序列(非参考线)。
如图4所示,为一个具体实施例中将上述磁共振成像方法应用于心脏延迟增强序列PSIR中的序列结构示意图。其中,IMG#(IMG1、IMG2)表示成像序列,IR表示反转恢复脉冲。
在本实施例中,以2个心动周期为单位为例说明,该2个心动周期包括时间上相邻的第1心动周期和第2心动周期,第1个心动周期内,施加一反转恢复脉冲,后面使用一成像序列(反转恢复脉冲与成像序列共同组成上述第二成像序列);第2个心动周期内,不施加反转恢复脉冲,直接使用一成像序列(上述第一成像序列),成像序列施加时间与第1个心动周期内的成像序列起始时间一样.在本实施例中,可以以第2个心动周期的成像序列的参考线(第一采集周期中全采样区域的数据)来重建第1个心动周期的缺失线(第二采集周期中欠采样区域的数据)。
上述磁共振成像方法中,第1个心动周期中仅需采集非参考线,可以提高采集效率,节省采集时间;第2个心动周期中未施加反转恢复脉冲时采集的参考线,可以提高参考线的信噪比,使得并行加速的计算结果更好。
在一个实施例中,准备模块包括T2准备模块;进一步地,在本实施例中,第一采集周期包括一次数据采集过程中的第一个采集周期;第二采集周期包括一次数据采集过程中的第四个采集周期以及第七个采集周期。
在一个实施例中,T2准备模块对应T2准备脉冲,可以理解地,第二成像序列中的采样模块为欠采样模块,对应欠采样的脉冲序列。
如图5所示,为一个具体实施例中将上述磁共振成像方法应用于T2Mapping中的序列结构示意图。其中,IMG#(IMG1、IMG2和IMG3)表示成像序列,TR_PreP#表示T2准备脉冲。
可选地,第一成像序列的施加时间与第二成像序列的施加时间在时间分布上也可间隔设置。在本实施例中,以7个心动周期为单位,第1个心动周期内,直接使用一成像序列(上述第一成像序列);第4个心动周期内,施加一T2准备脉冲,再使用一成像序列(T2准备脉冲与成像序列共同组成上述第二成像序列),成像序列施加时间与第1个心动周期内的成像序列起始时间一样;第7个心动周期内,施加另一T2准备脉冲,再使用一成像序列,成像序列施加时间与第1个心动周期内的成像序列起始时间一样;即第1个心动周期对应第一采样周期,第4、7个采样周期对应第二采样周期。
则在本实施例中,可以以第1个心动周期的成像序列中采集的参考线(第一采集周期中全采样区域的数据)来重建第4,第7个心动周期的缺失线(第二采集周期中欠采样区域的数据),从而获得本次数据采集过程中所有数据,进一步可以得到磁共振图像。
根据第一数据集中全采样区域的数据可获取部分校准系数,根据校准系数可对第一数据集中欠采样区域的数据进行恢复,获取填充完整的第一数据集。可选地,根据填充完整的第一数据集,获取完整校准系数,并根据完整校准系数对第二数据集中欠采样区域的数据进行恢复,计算获得校正的数据集(全部或局部为填充完整的第二数据集)。
上述磁共振成像方法中,第1个心动周期中未施加T2准备脉冲时采集的参考线信噪比较高,使得并行加速的计算结果更好;且第4、第7个心动周期中仅需采集非参考线,可以提高采集效率,节省采集时间。
上述磁共振成像方法可以应用于任何参考线信噪比不足的实施例中,并不仅限于上述实施例中举例说明的示例。
应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中,如图6所示,提供了一种磁共振成像装置,包括:第一成像序列激发模块、第二成像序列激发模块、数据恢复模块和成像模块,其中:
第一成像序列激发模块610,用于在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,利用第一成像序列,获得目标区域的第一数据集;第一数据集包括欠采样区域和全采样区域;第一成像序列包括第一采样模块。
第二成像序列激发模块620,用于在一次数据采集过程中的第二采集周期内利用第二成像序列,获得目标区域的第二数据集;第二数据集包括欠采样区域;第二成像序列包括准备模块和第二采样模块。
数据恢复模块630,用于根据第一数据集中的全采样区域的数据恢复所述第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集。
成像模块640,用于根据校正的数据集,获得目标区域的磁共振图像。
上述磁共振成像装置,在对检测对象的目标区域的一次采集过程中,第一采集周期内利用第一成像序列激发目标区域采集获得第一数据集,而在第二采集周期内,利用第二成像序列激发目标区域以获得第二数据集;其中,第一成像序列中不包含准备模块,第二成像序列中包含准备模块;第一数据集包括欠采样区域和全采样区域,第二数据集包括欠采样区域;然后通过第一数据集中的全采样区域的数据恢复所有欠采样区域的数据,获得校正的数据集,根据校正的数据集获得磁共振图像。如此,由于第一成像序列中不包含准备模块,在期间采集的全采样区域数据(即参考线对应的数据)的信噪比较高。
关于磁共振成像装置的具体限定可以参见上文中对于磁共振成像方法的限定,在此不再赘述。上述磁共振成像装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图7所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储磁共振图像数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种磁共振成像方法。
本领域技术人员可以理解,图7中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,利用第一成像序列,获得目标区域的第一数据集;第一数据集包括欠采样区域、全采样区域;第一成像序列包括第一采样模块;
在一次数据采集过程中的第二采集周期内利用第二成像序列,获得目标区域的第二数据集;第二数据集包括欠采样区域;第二成像序列包括准备模块和第二采样模块;
根据第一数据集中的全采样区域的数据恢复第一数据集或/和第二数据集中的欠采样区域的数据,以获得校正的数据集;
根据校正的数据集,获得目标区域的磁共振图像。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:目标区域中包含心脏,采集周期为心动周期。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:准备模块包括饱和恢复模块;第一采集周期包括一次数据采集过程中前预设数目个采集周期,第二采集周期包括第一采集周期以外的采集周期;
或者,第二采集周期包括一次数据采集过程中前预设数目个采集周期,第一采集周期包括第二采集周期以外的采集周期。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:在一次数据采集过程内,第一成像序列和第二成像序列均包括四个成像片层参数。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:准备模块包括反转恢复模块;第二采集周期包括一次数据采集过程中的第一个采集周期;第一采集周期与所述第二采集周期在时间上邻接。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:准备模块包括T2准备模块;第一采集周期包括一次数据采集过程中的第一个采集周期;第二采集周期包括一次数据采集过程中的第四个采集周期以及第七个采集周期。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:第一采样模块的施加起始时间与第二采样模块的施加起始时间一致。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:根据第一数据集中的全采样区域的数据恢复第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集,包括:
通过第一数据集中的全采样区域的数据获取恢复系数;根据恢复系数,恢复第一数据集和/或第二数据集的欠采样区域的数据,获得校正的数据集。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,利用第一成像序列,获得目标区域的第一数据集;第一数据集包括欠采样区域、全采样区域;第一成像序列包括第一采样模块;
在一次数据采集过程中的第二采集周期内,利用第二成像序列,获得目标区域的第二数据集;第二数据集包括欠采样区域;第二成像序列包括准备模块和第二采样模块;
根据第一数据集中的全采样区域的数据恢复第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集;
根据校正的数据集,获得目标区域的磁共振图像。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:目标区域中包含心脏,采集周期为心动周期。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:第一采集周期包括一次数据采集过程中前预设数目个采集周期,第二采集周期包括第一采集周期以外的采集周期;
或者,第二采集周期包括一次数据采集过程中前预设数目个采集周期,第一采集周期包括第二采集周期以外的采集周期;
准备模块包括饱和恢复模块。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:在一次数据采集过程内,第一成像序列和第二成像序列均包括四个成像片层参数。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:准备模块包括反转恢复模块;第二采集周期包括一次数据采集过程中的第一个采集周期;第一采集周期与第二采集周期在时间上邻接。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:准备模块包括T2准备模块;第一采集周期包括一次数据采集过程中的第一个采集周期;第二采集周期包括一次数据采集过程中的第四个采集周期以及第七个采集周期。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:第一采样模块的施加起始时间与第二采样模块的施加起始时间一致。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:根据第一数据集中的全采样区域的数据恢复第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集,包括:
通过第一数据集中的全采样区域的数据获取恢复系数;根据恢复系数,恢复第一数据集和/或第二数据集的欠采样区域的数据,获得校正的数据集。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种磁共振成像方法,包括:
在对检测对象的目标区域的一次数据采集过程中的第一采集周期内,利用第一成像序列,获得所述目标区域的第一数据集;所述第一数据集包括欠采样区域和全采样区域;所述第一成像序列包括第一采样模块;
在所述一次数据采集过程中的第二采集周期内,利用第二成像序列,获得所述目标区域的第二数据集;所述第二数据集包括欠采样区域;所述第二成像序列包括准备模块和第二采样模块;
根据所述第一数据集中的全采样区域的数据恢复所述第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集;
根据所述校正的数据集,获得所述目标区域的磁共振图像。
2.根据权利要求1所述的方法,其特征在于,所述目标区域中包含心脏,所述采集周期为心动周期。
3.根据权利要求1所述的方法,其特征在于,所述准备模块包括饱和恢复模块;
所述第一采集周期包括一次数据采集过程中前预设数目个采集周期,所述第二采集周期包括所述第一采集周期以外的采集周期;
或者,所述第二采集周期包括一次数据采集过程中前预设数目个采集周期,所述第一采集周期包括所述第二采集周期以外的采集周期。
4.根据权利要求3所述的方法,其特征在于,在所述一次数据采集过程内,所述第一成像序列和所述第二成像序列均包括四个成像片层参数。
5.根据权利要求1所述的方法,其特征在于,所述准备模块包括反转恢复模块;
所述第二采集周期包括一次数据采集过程中的第一个采集周期;第一采集周期与所述第二采集周期在时间上邻接。
6.根据权利要求1所述的方法,其特征在于,所述准备模块包括T2准备模块;
所述第一采集周期包括一次数据采集过程中的第一个采集周期;第二采集周期包括一次数据采集过程中的第四个采集周期以及第七个采集周期。
7.根据权利要求1至6任意一项所述的方法,其特征在于,所述第一采样模块的施加起始时间与所述第二采样模块的施加起始时间一致。
8.根据权利要求1至6任意一项所述的方法,其特征在于,所述根据所述第一数据集中的全采样区域的数据恢复所述第一数据集和/或第二数据集中的欠采样区域的数据,以获得校正的数据集,包括:
通过所述第一数据集中的全采样区域的数据获取恢复系数;
根据所述恢复系数,恢复所述第一数据集和/或第二数据集的欠采样区域的数据,获得所述校正的数据集。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至8中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至8中任一项所述的方法的步骤。
CN201910762725.7A 2019-08-19 2019-08-19 磁共振成像方法、计算机设备和计算机可读存储介质 Active CN110547799B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910762725.7A CN110547799B (zh) 2019-08-19 2019-08-19 磁共振成像方法、计算机设备和计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910762725.7A CN110547799B (zh) 2019-08-19 2019-08-19 磁共振成像方法、计算机设备和计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN110547799A CN110547799A (zh) 2019-12-10
CN110547799B true CN110547799B (zh) 2022-10-28

Family

ID=68737635

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910762725.7A Active CN110547799B (zh) 2019-08-19 2019-08-19 磁共振成像方法、计算机设备和计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN110547799B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023122919A1 (zh) * 2021-12-27 2023-07-06 中国科学院深圳先进技术研究院 心脏磁共振弥散张量成像方法、装置、设备及存储介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9103898B2 (en) * 2011-04-18 2015-08-11 General Electric Company Motion correction in accelerated T1-weighted magnetic resonance imaging
CN106491131B (zh) * 2016-12-30 2019-05-07 深圳先进技术研究院 一种磁共振的动态成像方法和装置
EP3376246B1 (de) * 2017-03-14 2023-04-26 Siemens Healthcare GmbH Beschleunigtes erzeugen einer serie von magnetresonanzbildern mit simultaner magnetresonanz-mehrschichtbildgebung

Also Published As

Publication number Publication date
CN110547799A (zh) 2019-12-10

Similar Documents

Publication Publication Date Title
US8581583B2 (en) Method and apparatus for magnetic resonance imaging to create T1 maps
CN111656392B (zh) 用于合成磁共振图像的系统和方法
CN110244246B (zh) 磁共振成像方法、装置、计算机设备和存储介质
US8588889B2 (en) Method and apparatus for breath-held MR data acquisition using interleaved acquisition
JP5599893B2 (ja) ナビゲータを使用するmrイメージング
CN107510458B (zh) 一种磁共振成像方法和设备
US7570985B2 (en) Method and apparatus for magnetic resonance imaging using contrast agent
JP2009160378A (ja) 磁気共鳴イメージング装置
JP3556176B2 (ja) Mr投影画像の自動修正
CN109219757A (zh) Dixon型水/脂肪分离MR成像
CN111257810B (zh) 磁共振信号采集系统、方法
US9791532B2 (en) Method and apparatus for magnetic resonance imaging
US8143891B2 (en) System for image acquisition with fast magnetic resonance gradient echo sequences
CN110547799B (zh) 磁共振成像方法、计算机设备和计算机可读存储介质
JP2006021023A (ja) 磁気共鳴診断装置
US20170261582A1 (en) Method and apparatus for generating a magnetic resonance image data set
US7706855B1 (en) System and method for MR data acquisition with uniform fat suppression
CN112649775B (zh) 磁共振成像控制方法、装置
US10712415B2 (en) Method and apparatus for recording magnetic resonance data
CN114114115A (zh) 心脏图像获取方法、装置、磁共振设备和存储介质
US10061006B2 (en) System and method for accelerated angiographic magnetic resonance imaging
US20170231523A1 (en) Method and magnetic resonance apparatus scar quantification in the myocardium
JP4086544B2 (ja) 磁気共鳴イメージング装置
WO2009047690A2 (en) Segmented multi-shot mri involving magnetization preparation
WO2015031408A1 (en) Multi-spectral mri scan with magnetization recovery

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant