CN108175409B - 一种定量快速锁频磁共振成像方法 - Google Patents

一种定量快速锁频磁共振成像方法 Download PDF

Info

Publication number
CN108175409B
CN108175409B CN201810012513.2A CN201810012513A CN108175409B CN 108175409 B CN108175409 B CN 108175409B CN 201810012513 A CN201810012513 A CN 201810012513A CN 108175409 B CN108175409 B CN 108175409B
Authority
CN
China
Prior art keywords
spin
locking
magnetic resonance
resonance imaging
pulses
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
CN201810012513.2A
Other languages
English (en)
Other versions
CN108175409A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201810012513.2A priority Critical patent/CN108175409B/zh
Publication of CN108175409A publication Critical patent/CN108175409A/zh
Application granted granted Critical
Publication of CN108175409B publication Critical patent/CN108175409B/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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (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

本发明公开了一种定量快速锁频磁共振成像方法,在连续扫描周期内完成不同自旋锁定脉冲对应的磁共振成像信号的数据采集,并将其填充到K空间的同一位置,使得最终得到的不同自旋锁定脉冲对应的Tmapping图像匹配一致度高,可以用于评判组织纤维化指数和诊断组织水肿状态;且使得心律不齐、呼吸不均匀等对Tmapping图像质量影响降到最低,对受试者的身体状态及闭气等能力没有严格要求,使获得Tmapping图像所需扫描时间大大缩短,从而使Tmapping成像技术在临床医学检测领域研究和推广得到进一步发展;此外,相邻RF激励脉冲之间存在的自旋锁定脉冲相位正负交替,从而有效减少因磁场均匀性产生的伪影。

Description

一种定量快速锁频磁共振成像方法
技术领域
本发明属于医学诊断技术领域,涉及临床心脏等需心电、呼吸门控导航的磁共振扫描方法。
背景技术
磁共振成像(MRI,Magnetic Resonance Imaging)近年来被广泛运用于疾病临床检查和实验研究中,它是通过检测不同组织间T1弛豫时间(即自旋-晶格弛豫时间)和T2弛豫时间(自旋-自旋弛豫时间)的差异,并经傅里叶变换得到相关组织的图像,从而以图像形式显示各种组织及其不同病理引起的改变等。
自旋锁定脉冲是指MRI系统下由横向激励磁场B1(B1≥0)作用产生的低幅度射频(RF)脉冲,是一种谐振而且连续的射频脉冲波,其具有特定的持续时间和较低的能量。自旋锁定T成像,是在满足自旋锁定条件下,强制MRI系统用于产生磁共振信号的射频线圈生成的磁化矢量在特定射频激励条件下进行弛豫,得到在自旋锁定时间(Time of SpinLocking,TSL)采集射频脉冲所产生的横向平面磁化矢量,进行傅里叶变换,从而得到相关组织的图像,再进行数学变换,得到不同组织弛豫时间T(即旋转坐标系的自旋-晶格弛豫时间)图像。弛豫时间T主要表现细胞外基质(例如蛋白聚糖)存在条件下的水中氢质子的弛豫特性,可以用于组织中大分子成分及不同分子之间质子交换的分析研究,现已被证实可用于检测组织内胶原和软骨内蛋白多糖含量。与MRI T2成像技术相比,自旋锁定T成像技术通过“锁定”总体磁化矢量,从而防止磁化矢量的能量流失,减少T2(横向弛豫时间)的衰减过程,因此T>T2。如果把图像的每一个点都算出T,这就形成T mapping的图像。
目前,T mapping技术已运用于脑、骨关节、肝脏、肾脏、心脏等部位的检查,用于检测胶原和软骨蛋白多糖含量,这对于判断组织胶原成分、判断组织纤维化程度具有十分重要的价值。然而,T mapping技术由于Tmapping序列(检查时磁场加载顺序)是在一定的周期内完成的(例如心动周期、呼吸周期等),使其应用却受机体心脏运动、呼吸运动等自主运动的影响,使其在心脏、肝脏部位采集时应用受限,尤其是心脏T mapping技术运行极其困难,仅有个别团队进行有关心脏T mapping技术的应用报道(van Oorschot etal.Endogenous assessment of chronic myocardial infarction with T(1ρ)-mappingin patients.Journal of Cardiovascular Magnetic Resonance 2014,16:104;ChunhuaWang et al.Endogenous contrast T1rho cardiac magnetic resonance formyocardial fibrosis in hypertrophic cardiomyopathy patients.J Cardiol.2015Dec;66(6):520-6)。
所报道的心脏T mapping序列为心电呼吸门控触发扫描,在吸气末屏气扫描获得T mapping加权图像,具体实现方式为:在主磁体产生的静态主磁场内,以主磁场所在坐标系为参考系,当横向激励磁场B1=500Hz(12μT)时,按照设定的T mapping序列对被检测对象施加两个射频激励脉冲(90°、180°,将宏观磁化矢量M0偏转90°或180°),连续采集五个TSL时间(TSL1=1msec、TSL2=10msec、TSL3=20msec、TSL4=30msec、TSL5=0msec),获取T加权图像,并将获得的T加权图像代入数学模型拟合,得到T mapping数据(vanOorschot et al.Endogenous assessment of chronic myocardial infarction with T(1ρ)-mapping in patients.Journal of Cardiovascular Magnetic Resonance 2014,16:104)。虽然通过上述方法可以获得T mapping定量数据,然而,上述T mapping成像方法仍存在以下问题:(1)心电、呼吸情况常常因人而异,心电-呼吸触发不良、患者屏气时间不恒定、患者心率不齐等都容易使获取的Tmapping图像存在伪影,影响成像质量;(2)常因设备磁体原因(如不满足B1磁场所需要求),需更改B1值才能继续检测实现图像扫描;(3)为了扫描时的比吸收率(SAR值)符合国际指导准则,成像过程扫描时间较长(当ECG状态良好,获得一幅T mapping图像的扫描总时间大约为40s;获得心脏整体Tmapping图像需要六幅T mapping图像,当ECG紊乱时,则无法获得较好质量图像,且扫描时间将延长,甚至中断扫描,因此总扫描时间远远长于6min),这对患者的闭气能力要求较高,对于医患重疾的患者来讲,一般都难以满足要求,从而使T mapping技术难以有应用价值。而以上Tmapping技术中存在的图像伪影高、患者配合或状态较差等问题也是T mapping技术在研发和推广中难以取得进一步进展的原因。
发明内容
本发明的目的旨在针对上述现有技术中存在的问题,提供一种定量快速锁频磁共振成像方法,不仅可以避免患者配合困难或状态差对时T成像质量的影响,减少图像伪影,而且可以提高不同自旋锁定脉冲条件下的图像匹配度及缩短获得Tmapping图像所需扫描时间。
本发明提供的定量快速锁频磁共振成像方法,采用磁共振成像系统对受试者进行多个循环周期的扫描,所述磁共振成像系统包括射频线圈、脉冲发生器、接收器以及数据处理和图像重建模块;所述射频线圈用于向检查区域的受试者发射射频电磁信号,或者采集来自受试者的射频电磁信号;所述脉冲发生器用于生成经由射频线圈发射的至少两种可形成检查区域内B1磁场的自旋锁定脉冲,以及施加于检查区域受试者的多个RF激励脉冲,所述接收器用于接收经由射频线圈采集的来自受试者的射频电磁信号,并将其作为磁共振成像信号;所述数据处理和图像重建模块用于将采集的磁共振成像信号进行重建得到最终的Tmapping图像;一个循环周期由多个扫描周期构成,一个循环周期包含多个设定的自旋锁定时间TSL和多个设定的自旋锁定脉冲B1
一个扫描周期内的操作步骤包括:
(1)经设定扫描延迟时间后,触发扫描信号;
(2)在设定自旋锁定时间TSL内,按照设定脉冲序列,控制脉冲发生器生成经由射频线圈发射的设定自旋锁定脉冲B1和施加于检查区域受试者的多个RF激励脉冲;
(3)由接收器接收经由射频线圈采集的磁共振成像信号,并将接收的磁共振成像信号读出填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间形成数据点阵;
一个循环周期内的操作步骤包括:
(S1)针对不同自旋锁定脉冲B1,采用同一个设定自旋锁定时间TSL,通过多个扫描周期依次完成所有不同自旋锁定脉冲B1对应的扫描,并将采集的磁共振成像信号填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间的同一位置;
(S2)对设定的不同自旋锁定时间TSL,重复步骤(S1),完成不同自旋锁定时间TSL对应的扫描,即完成一个循环周期;
当完成对一个循环周期的扫描之后依次进入下一循环周期的扫描,循环周期的次数由填满自旋锁定时间和自旋锁定脉冲B1对应的K空间决定,K空间填满之后,结束扫描,并将K空间中的数据进行处理得到不同自旋锁定脉冲B1对应的T mapping图像。
上述定量快速锁频磁共振成像方法,所述步骤(1)中,以心动周期或呼吸周期作为扫描周期,利用采集的心电信号或呼吸信号,根据被检测者心电或呼吸状态设定延迟触发时间,使扫描发生在最小的运动周期内,以极大限度地减少运动伪影。本发明扫描延迟时间由当前心动周期或呼吸周期之前的N个心动周期或呼吸周期的平均值
Figure BDA0001540793400000032
确定,延迟时间
Figure BDA0001540793400000031
Ti为第i个扫描周期,N≥2,T0为一个心动周期内设定的数据扫描时间或一个呼吸周期内设定的数据扫描时间。
上述定量快速锁频磁共振成像方法,所述步骤(1)中,触发扫描信号由磁共振成像系统中的自触发导航模块执行,优选为型号为SG的触发模块,该模块可以测定心脏组织在每次磁共振数据采集时的位置,这样在全部数据采集完后,进行数据后处理,把在不同位置上采集的K空间数据作修正(Pang J,Chen Y,Fan Z,Nguyen C,Yang Q,Xie Y,Li D.Highefficiency coronary MR angiography with nonrigid cardiac motioncorrection.Magn Reson Med.2016;76(5):1345-1353),再按照如下K空间的填充策略,形成一个完整的K空间图,通过傅立叶变换等数据处理,得到最后的Tmapping图像。由于可以利用采集的心脏组织位置数据对K空间数据进行修正,所以可更有效获取扫描图像。本发明采用的K空间的填充策略为:将与自旋锁定时间和自旋锁定脉冲B1对应的K空间数据分为m部分(即扫描次数为m),每部分有n条K-空间线(即一次扫描获得的K-空间线条数),每一扫描周期填充一个设定自旋锁定时间及一个自旋锁定脉冲B1对应的K-空间的一部分。
上述定量快速锁频磁共振成像方法,所述步骤(2)中,在控制脉冲发生器激发脉冲之前,先消除脂肪信号,以避免脂肪信号对待检查组织信号的干扰。本发明采用目前系统带有的常规压脂模块来消除脂肪信号;之后再控制脉冲发生器按照设定脉冲序列生成自旋锁定脉冲和RF激励脉冲。
上述定量快速锁频磁共振成像方法,所述步骤(2)中,设定脉冲序列包括选择性激励核素的多个RF激励脉冲,相邻两个RF激励脉冲之间设置至少一个自旋锁定脉冲B1。多个RF激励脉冲由两个90°激励脉冲和至少两个180°翻转聚焦脉冲组成,多个RF激励脉冲的激发顺序为90°激励脉冲-180°翻转聚焦脉冲-…-180°翻转聚焦脉冲-90°激励脉冲;优选的实施方式中,180°翻转聚焦脉冲的个数为4个。进一步地,当相邻RF激励脉冲之间存在两个及以上的自旋锁定脉冲时,自旋锁定脉冲相位正负交替,以减少由RF激励脉冲产生的磁场到自旋锁定脉冲产生的磁场转变时由于磁场异质性以及非共振效应所产生的伪影。
上述定量快速锁频磁共振成像方法,设定的自旋锁定时间TSL数量和设定的自旋锁定脉冲B1数量均至少为两个。
上述定量快速锁频磁共振成像方法,将不同自旋锁定脉冲B1对应的K空间中的数据进行处理得到不同自旋锁定脉冲B1对应的Tmapping图像的过程包括以下分步骤:
A、将自旋锁定时间和自旋锁定脉冲B1对应的K空间中的数据进行傅里叶变换以及数模转换得到相应的T加权图像,并按照自旋锁定脉冲B1大小对所述T加权图像进行分组;
B、对各组T加权图像分别依据K空间同一位置转换得到的T加权图像中T信号强度S随自旋锁定时间TSL变化关系,拟合得到T值变化曲线
Figure BDA0001540793400000041
从而确定T值;
C、重复步骤B,得到K空间不同位置对应的T值,由K空间不同位置的T值组成Tmapping图像。
通过上述步骤A-步骤C,可以准确获得T的定量值,进而得到成像质量较高的Tmapping图像,为临床医学上不同组织(例如心脏、肝脏等)的纤维化指数(=T,当B1=12μT时;=-T,当B1=0μT时)评判和水肿状态诊断提供准确的分析数据。
与现有技术相比,本发明具有以下有益效果:
1、本发明定量快速锁频磁共振成像方法,在连续扫描周期内完成不同自旋锁定脉冲对应的磁共振成像信号的数据采集,并将其填充到K空间的同一位置,使得最终得到的不同自旋锁定脉冲对应的Tmapping图像匹配一致度高;
2、本发明定量快速锁频磁共振成像方法在连续扫描周期内完成不同自旋锁定脉冲对应的磁共振成像信号的数据采集,将其填充到K空间的同一位置,使得心律不齐、呼吸不均匀等对Tmapping图像质量影响降到最低,并采用自触发导航模块进行修正运动伪影,充分利用每一个采集的数据,对受试者的身体状态及闭气等能力没有要求,从而使获得Tmapping图像所需扫描时间大大缩短,在任何情况下均可以保证6分钟左右完成对心脏的整体扫描获得成像质量较高的Tmapping图像;
3、本发明定量快速锁频磁共振成像方法,由于采用自触发导航模块,控制脉冲发生器进行自触发扫描,从而使扫描过程中可规避心电-呼吸触发不良、患者闭气时间不恒定、患者心律不齐等所带来的图像干扰,减少有效心电信号丢失,可适用于因各种心脏病引起的心律不齐患者;
4、本发明定量快速锁频磁共振成像方法,多个RF激励脉冲中的第二激励脉冲相位正负交替;且当相邻RF激励脉冲之间存在两个及以上的自旋锁定脉冲时,自旋锁定脉冲相位正负交替,从而有效减少因磁场不均匀性产生的伪影。
附图说明
图1为本发明实施例中频磁共振成像过程示意图。
图2为本发明实施例中在自旋锁定时间内设定的脉冲序列示意图。
图3为本发明实施例中一个循环周期流程示意图。
图4为本发明实施例中与自旋锁定时间TSL和自旋锁定脉冲B1对应的K空间填充示意图。
图5为本发明实施例中Tmapping拟合曲线图。
具体实施方式
以下将结合附图通过实施例对本发明的技术方案进行清楚、完整的描述,显然,所描述实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例,都属于本发明所保护的范围。
本实施例基于临床常见3.0T磁共振系统对定量快速锁频磁共振成像过程进行详细说明。所述磁共振成像系统的主要参数包括:宽带=401Hz/pixel、TE=1.6msec、层厚=8mm、FOV=320mm×240mm、分辨率=1.3mm×1.3mm。所述磁共振成像系统包括磁体、磁共振波谱仪、扫描控制器、数据处理和图像重建模块、压脂模块以及自触发导航模块(SG);所述磁体包括主磁体和射频线圈,所述主磁体用于产生静磁场,所述射频线圈用于根据脉冲发生器发出的脉冲信号向人体辐射出指定频率和指定功率的射频电磁信号,或者采集来自受试者的射频电磁信号;所述磁共振波谱仪包括脉冲发生器和接收器,所述脉冲发生器用于生成经由射频线圈发射的至少两种自旋锁定脉冲B1和施加于检查区域的多个RF激励脉冲,所述接收器用于接收经由射频线圈采集的来自受试者的射频电磁信号,并将其作为磁共振成像信号;所述扫描控制器用于控制脉冲发生器按照设定脉冲序列产生脉冲信号,并同时控制接收器采集磁共振成像信号;所述数据处理和图像重建模块用于将采集的磁共振成像信号进行重建得到最终的Tmapping图像;所述压脂模块用于消除脂肪信号;所述自触发导航模块(SG)用于将触发扫描信号发送给扫描控制器。
本实施例以心电监测仪获取受试者的心电(ECG)信号,以获得的一个心动周期为一个扫描周期。
本实施例设定四个自旋锁定时间TSL(M=4),分别为TSL1=0msec,TSL2=30msec,TSL3=50msec,TSL2=70msec。
本实施例设定两种自旋锁定脉冲,对应的磁场B1大小分别为0和12μT。
本实施例提供的定量快速锁频磁共振成像方法目的是得到Tmapping图像,是通过对受试者相应部位(例如心脏、肝脏等)进行多个循环周期的扫描得到的;一个循环周期由多个扫描周期(即这里的心动周期)构成。
如图1所示,每一个心动周期内的操作包括步骤(1)-步骤(3)。
(1)经设定扫描延迟时间后,触发扫描信号。
本实施例根据前五个心动周期(N=5)的平均值,决定下一个扫描周期的延迟时间,延迟时间
Figure BDA0001540793400000061
Ti为第i个扫描周期时长,N=5,
Figure BDA0001540793400000062
为下一个扫描周期之前五个扫描周期的总时长,T0为一个心动周期内设定的数据扫描时间,一般为500~700ms。本实施例根据ECG状态设定扫描延迟时间,从而使扫描发生在心脏的最小运动周期内。
在一个心动周期内,扫描延迟时间结束后,由自触发导航模块(SG)执行触发,将触发扫描信号发送给扫描控制器。
自触发导航模块SG,还可以对所得到的K空间数据进行修正,使T mapping扫描运行过程中更有效获取采集得到的MRI扫描图像,指导后序采集的图像信号在K-空间的填充方式。本实施例中根据ECG情况,采用的K-空间填充方式为:以每幅T加权图K-空间需要生成128根K-空间线为例,每次K-空间填充数为8根,那么就需要16次扫描才能完成K-空间填充。
(2)在设定自旋锁定时间TSL内,按照设定脉冲序列,扫描控制器控制脉冲发生器生成经由射频线圈发射的设定自旋锁定脉冲B1和施加于检查区域受试者的多个RF激励脉冲。
在控制脉冲发生器激发脉冲之前,首先运行磁共振成像系统的压脂模块消除脂肪信号。
在消除脂肪信号之后,进入脉冲发射阶段,本实施例以主磁场X-Y-Z坐标系为参考,脉冲序列包括两个90°激励脉冲(90°表示射频脉冲磁化矢量相对于主磁场磁场方向偏转的角度为90°)和四个180°翻转聚焦脉冲(180°表示射频脉冲磁化矢量相对于主磁场磁场方向偏转的角度为180°)以及多个自旋锁定脉冲B1。在一个自旋锁定时间(TSL)内时间首先沿X轴正方向给入一个90°激励脉冲(作为导航回波),间隔TSL/8后沿Y轴正方向给入绝热180°翻转聚焦脉冲,后间隔TSL/4沿Y轴正方向给入绝热180°翻转聚焦脉冲,随后隔TSL/4时沿Y轴负方向给入绝热180°翻转聚焦脉冲,再间隔TSL/4时沿Y轴负方向给入绝热180°翻转聚焦脉冲,最后在间隔TSL/8时沿X轴负方向给入第二个90°激励脉冲,在第一个90°激励脉冲与第一个180°聚焦脉冲间存在一个沿Y轴正方向的自旋锁定脉冲B1,四个180°聚焦脉冲间均存在两个自旋锁定脉冲B1,其相位正负轮流,第四个180°聚焦脉冲与第二个90°激励脉冲间存在一个沿y轴正方向的自旋锁定脉冲B1。自旋锁定脉冲B1相位方向正负轮流,可以减少从RF射频激励脉冲到自旋锁定脉冲B1异质性所产生的伪影,以及非共振效应所导致的伪影。
(3)由接收器接收经由射频线圈采集的磁共振成像信号,并将接收的磁共振成像信号读出填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间中形成数据点阵。
在发射第二个90°激励脉冲发射之后,脉冲发射阶段结束,在之后的读出间隔通过接收器接收经由射频线圈采集的磁共振信号填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间中。
一个循环周期内的操作包括步骤(S1)-步骤(S2)。
(S1)针对不同自旋锁定脉冲B1,采用同一个设定自旋锁定时间TSL,通过多个扫描周期依次完成所有不同自旋锁定脉冲B1对应的扫描,并将采集的磁共振成像信号填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间的同一位置。
如前所述,本实施例设定了脉冲强度不同的两种自旋锁定脉冲和四个自旋锁定至时间。如图3所示,针对不同自旋锁定脉冲B1,采用同一个自旋锁定时间TSL,通过相邻两个心动周期依次完成两个自旋锁定脉冲B1对应的扫描,并将采集的磁共振成像信号填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间的同一位置。以本实施例自旋锁定时间TSL1为例,在第一个心动周期,B1=0,按照步骤(1)-步骤(3)完成扫描,并将采集的磁共振成像信号填充到B1=0和TSL1对应的K-空间,获得B1=0和TSL1对应的n条K-空间线(如图4所示);在第二个心动周期,B1=12μT,按照步骤(1)-步骤(3)完成扫描,并将采集的磁共振成像信号填充到B1=12μT和TSL1对应的K-空间(且与B1=0μT和TSL1对应的K-空间的填充位置相同),获得B1=12μT和TSL1对应的n条K-空间线(如图4所示)。
(S2)对设定的不同自旋锁定时间TSL,重复步骤(S1),完成不同自旋锁定时间TSL对应的扫描,即完成一个循环周期。
如前所述,已经给出了本实施例扫描过程中针对两个不同自旋锁定脉冲B1,采用同一个自旋锁定时间TSL1,完成第一个心动周期和第二个心动周期的扫描。第三个心动周期和第四个心动周期与第一个心动周期和第二个心动周期类似,不同的是采用的是自旋锁定时间TSL2,具体为:在第三个心动周期,B1=0,按照步骤(1)-步骤(3)完成扫描,并将采集的磁共振成像信号填充到B1=0和TSL2对应的K-空间,获得B1=0和TSL2对应的n条K-空间线(如图4所示);在第四个心动周期,B1=12μT,按照步骤(1)-步骤(3)完成扫描,并将采集的磁共振成像信号填充到B1=12μT和TSL2对应的K-空间,获得B1=12μT和TSL2对应的n条K-空间线(如图4所示);以此类推,直到获得B1=12μT和TSL4对应的n条K-空间线,完成一个循环周期的扫描。之后依次进入下一个循环周期的扫描,循环周期的次数由填满自旋锁定时间和自旋锁定脉冲B1对应的K空间来决定。这样可以使得数据获取最高的空间匹配度,使得心律不齐、呼吸不均匀对获取T mapping图像影响降到最低。本实施例提供的填充方式是沿平行于Kx方向进行填充,本领域技术人员也可以根据实际情况选择沿平行于Ky的方向填充的填充方式或者放射状填充方式等。
通过以上方法可以获得8个K-空间填充数据。本实施例利用数据处理和图像重建模块将8个K-空间填充数据经傅里叶变换和数模转换,获得8个T加权图像。数据处理方法可以参考(Chunhua Wang et al.Endogenous contrast T1rho cardiac magneticresonance for myocardial fibrosis in hypertrophic cardiomyopathy patients.JCardiol.2015 Dec;66(6):520-6)。
继后将8个T加权图像中同一位置的信号强度S提取出来,并将T信号强度S按照自旋锁定脉冲B1=0或12μT分为两组,两组分别按照T信号强度S随自旋锁定时间TSL变化关系,利用合适的指数拟合算法进行拟合,获得如图5所示的曲线拟合图像,拟合公式为
Figure BDA0001540793400000091
从而确定该位置的T值。重复上述过程,由T加权图像不同位置(即K空间不同位置)计算得到K空间不同位置对应的的T值,由K空间不同位置的T值组成相应的Tmapping图像。对于不同自旋锁定脉冲B1,采用相同的方式处理,最终得到B1=0及B1=12μT对应的Tmapping图像。

Claims (8)

1.一种定量快速锁频磁共振成像方法,采用磁共振成像系统对受试者进行多个循环周期的扫描,所述磁共振成像系统包括射频线圈、脉冲发生器、接收器以及数据处理和图像重建模块;所述射频线圈用于向检查区域的受试者发射射频电磁信号,或者采集来自受试者的射频电磁信号;所述脉冲发生器用于生成经由射频线圈发射的至少两种可形成检查区域内B1磁场的自旋锁定脉冲,以及施加于检查区域受试者的多个RF激励脉冲,所述接收器用于接收经由射频线圈采集的来自受试者的射频电磁信号,并将其作为磁共振成像信号;所述数据处理和图像重建模块用于将采集的磁共振成像信号进行重建得到最终的Tmapping图像;其特征在于一个循环周期由多个扫描周期构成,一个循环周期包含多个设定的自旋锁定时间TSL和多个设定的自旋锁定脉冲B1
一个扫描周期内的操作步骤包括:
(1)经设定扫描延迟时间后,触发扫描信号;
(2)在设定自旋锁定时间TSL内,按照设定脉冲序列,控制脉冲发生器生成经由射频线圈发射的设定自旋锁定脉冲B1和施加于检查区域受试者的多个RF激励脉冲;
(3)由接收器接收经由射频线圈采集的磁共振成像信号,并将接收的磁共振成像信号读出填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间形成数据点阵;
一个循环周期内的操作步骤包括:
(S1)针对不同自旋锁定脉冲B1,采用同一个设定自旋锁定时间TSL,通过多个扫描周期依次完成所有不同自旋锁定脉冲B1对应的扫描,并将采集的磁共振成像信号填充到数据处理和图像重建模块中与自旋锁定时间和自旋锁定脉冲B1对应的K空间的同一位置;
(S2)对设定的不同自旋锁定时间TSL,重复步骤(S1),完成不同自旋锁定时间TSL对应的扫描,即完成一个循环周期;
当完成对一个循环周期的扫描之后依次进入下一循环周期的扫描,循环周期的次数由填满自旋锁定时间和自旋锁定脉冲B1对应的K空间决定,K空间填满之后,结束扫描,并将K空间中的数据进行处理得到不同自旋锁定脉冲B1对应的Tmapping图像。
2.根据权利要求1所述定量快速锁频磁共振成像方法,其特征在于所述步骤(1)中,以心动周期或呼吸周期作为扫描周期,扫描延迟时间由当前心动周期或呼吸周期之前的N个心动周期或呼吸周期的平均值
Figure FDA0002824440450000011
确定,延迟时间
Figure FDA0002824440450000012
Ti为第i个扫描周期,N≥2,T0为一个心动周期内设定的数据扫描时间或一个呼吸周期内设定的数据扫描时间。
3.根据权利要求1所述定量快速锁频磁共振成像方法,其特征在于所述步骤(1)中,触发扫描信号由磁共振成像系统中的自触发导航模块执行。
4.根据权利要求1所述定量快速锁频磁共振成像方法,其特征在于所述步骤(2)中,在控制脉冲发生器激发脉冲之前,消除脂肪信号。
5.根据权利要求1至4任一权利要求所述定量快速锁频磁共振成像方法,其特征在于所述步骤(2)中,设定脉冲序列包括选择性激励核素的多个RF激励脉冲,相邻两个RF激励脉冲之间设置至少一个自旋锁定脉冲B1
6.根据权利要求5所述定量快速锁频磁共振成像方法,其特征在于多个RF激励脉冲由两个90°激励脉冲和至少两个180°翻转聚焦脉冲组成,多个RF激励脉冲的激发顺序为90°激励脉冲-180°翻转聚焦脉冲-…-180°翻转聚焦脉冲-90°激励脉冲。
7.根据权利要求6所述定量快速锁频磁共振成像方法,其特征在于相邻RF激励脉冲之间两个及以上的自旋锁定脉冲相位正负交替。
8.根据权利要求5所述定量快速锁频磁共振成像方法,其特征在于将不同自旋锁定脉冲B1对应的K空间中的数据进行处理得到不同自旋锁定脉冲B1对应的Tmapping图像的过程包括以下分步骤:
A、将自旋锁定时间和自旋锁定脉冲B1对应的K空间中的数据进行傅里叶变换以及数模转换得到相应的T加权图像,并按照自旋锁定脉冲B1大小对所述T加权图像进行分组;
B、对各组T加权图像分别依据K空间同一位置转换得到的T加权图像中T信号强度S随自旋锁定时间TSL变化关系,拟合得到T值变化曲线
Figure FDA0002824440450000021
从而确定T值;
C、重复步骤B,得到K空间不同位置对应的T值,由K空间不同位置的T值组成Tmapping图像。
CN201810012513.2A 2018-01-05 2018-01-05 一种定量快速锁频磁共振成像方法 Active CN108175409B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810012513.2A CN108175409B (zh) 2018-01-05 2018-01-05 一种定量快速锁频磁共振成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810012513.2A CN108175409B (zh) 2018-01-05 2018-01-05 一种定量快速锁频磁共振成像方法

Publications (2)

Publication Number Publication Date
CN108175409A CN108175409A (zh) 2018-06-19
CN108175409B true CN108175409B (zh) 2021-03-23

Family

ID=62549995

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810012513.2A Active CN108175409B (zh) 2018-01-05 2018-01-05 一种定量快速锁频磁共振成像方法

Country Status (1)

Country Link
CN (1) CN108175409B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111568420B (zh) * 2019-02-15 2021-05-14 深圳先进技术研究院 一种旋转坐标内弛豫时间确定方法、装置及存储介质
CN109998548B (zh) * 2019-04-17 2020-09-18 清华大学 定量心肌磁共振成像方法、设备及存储介质
US11280867B2 (en) * 2019-11-08 2022-03-22 The Chinese University Of Hong Kong System and method for quantitative magnetization transfer imaging based on spin-lock
CN111090069B (zh) * 2019-11-21 2022-03-29 深圳先进技术研究院 定量磁共振成像参数确定方法、装置、设备及存储介质
CN113050009B (zh) * 2019-12-27 2022-12-13 深圳先进技术研究院 三维磁共振快速参数成像方法和装置
WO2021217509A1 (zh) * 2020-04-29 2021-11-04 深圳先进技术研究院 一种磁共振旋转坐标系下的自旋晶格弛豫成像方法和系统
CN113567901B (zh) * 2020-04-29 2024-08-20 深圳先进技术研究院 一种磁共振旋转坐标系下的自旋晶格弛豫成像方法和系统
CN112132825A (zh) * 2020-10-10 2020-12-25 青岛市市立医院 一种基于T2 mapping及T1ρ关节软骨定量成像分析方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5420510A (en) * 1993-08-13 1995-05-30 Wisconsin Alumni Research Foundation Suppression of oscillations in NMR measurement using off-resonance spin-locking pulses
US6836114B2 (en) * 2002-03-15 2004-12-28 The Trustees Of The University Of Pennsylvania Pulse imaging sequences and methods for T1p-weighted MRI
CN103328999A (zh) * 2011-01-25 2013-09-25 皇家飞利浦有限公司 交错的自旋锁定成像
US20120226136A1 (en) * 2011-03-04 2012-09-06 The Chinese University Of Hong Kong Methods and devices of detection, grading, monitoring, and follow-up of fibrosis
JP5808659B2 (ja) * 2011-12-12 2015-11-10 株式会社日立メディコ 磁気共鳴イメージング装置及びT1ρイメージング法

Also Published As

Publication number Publication date
CN108175409A (zh) 2018-06-19

Similar Documents

Publication Publication Date Title
CN108175409B (zh) 一种定量快速锁频磁共振成像方法
Usman et al. Free breathing whole-heart 3D CINE MRI with self-gated Cartesian trajectory
JP6084573B2 (ja) マルチポイントディクソン技術を用いるmr撮像
US7863895B2 (en) System, program product, and method of acquiring and processing MRI data for simultaneous determination of water, fat, and transverse relaxation time constants
CN102914754B (zh) 用于产生磁共振照片的方法
CN105793722B (zh) 用于稳态mr序列的实时自适应生理同步和门控
CN110031786B (zh) 磁共振图像重建方法、磁共振成像方法、设备及介质
US9081073B2 (en) System for suppression of artifacts in MR imaging
US20150309147A1 (en) Systems and Methods For Designing Magnetic Resonance Imaging Radio Frequency Pulses That Are Robust Against Physiological Motion Errors
CN105283774B (zh) 在不利的心脏条件下改进受试者心脏成像的系统与方法
EP2610632A1 (en) MRI with Dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field
US7501821B2 (en) Method for determination of the transverse relaxation time T2* in MR data
CN114076911B (zh) 磁共振成像的方法及系统
CN108742626B (zh) T1参数图成像方法及磁共振成像系统
JP2006501919A (ja) 磁気共鳴方法及び装置
US9081055B2 (en) Method for reducing local specific absorption rate in magnetic resonance imaging using radio frequency coil array dark modes
US8928317B2 (en) System and method for controlling apparent timing dependencies for T2-weighted MRI imaging
JP2003325479A (ja) 撮像データの補正方法
CN106725508B (zh) 生理运动数据采集方法、磁共振成像方法以及装置
Kubach et al. Cardiac phase-specific shimming (CPSS) for SSFP MR cine imaging at 3 T
US20170115367A1 (en) Magnetic resonance imaging method and apparatus with motion-corrected model-based acceleration of parameter mapping
WO2021012972A1 (zh) 核磁共振成像控制方法、装置及计算机可读存储介质
US20050237059A1 (en) Method and apparatus for suppression of artifacts in MRT imaging
US20240201302A1 (en) Systems and methods for non-selective stimulated echo multislice diffusion imaging
CN112986878B (zh) 冠脉成像方法以及磁共振成像系统

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