CN102859386B - 核自旋系统的纵向和横向弛豫时间的同时和动态确定 - Google Patents

核自旋系统的纵向和横向弛豫时间的同时和动态确定 Download PDF

Info

Publication number
CN102859386B
CN102859386B CN201180014334.9A CN201180014334A CN102859386B CN 102859386 B CN102859386 B CN 102859386B CN 201180014334 A CN201180014334 A CN 201180014334A CN 102859386 B CN102859386 B CN 102859386B
Authority
CN
China
Prior art keywords
dynamic
nuclear spin
spin system
longitudinal
baseline
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.)
Expired - Fee Related
Application number
CN201180014334.9A
Other languages
English (en)
Other versions
CN102859386A (zh
Inventor
S·雷梅尔
刘巍
T·R·福格特
C·施特宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN102859386A publication Critical patent/CN102859386A/zh
Application granted granted Critical
Publication of CN102859386B publication Critical patent/CN102859386B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G01R33/50NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5601Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
    • 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
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5613Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • 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
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明涉及一种用于在DCE或DSE MRI的语境下同时和动态确定对象的核自旋系统的纵向弛豫时间T1和横向弛豫时间T2的磁共振成像方法。对此,本发明利用了包括EPI读出模块的稳态梯度回波脉冲序列。

Description

核自旋系统的纵向和横向弛豫时间的同时和动态确定
技术领域
本发明涉及一种用于同时确定对象的核自旋系统的纵向弛豫特性(behavior)和横向弛豫特性的磁共振成像方法、一种计算机程序产品以及一种磁共振成像设备。
背景技术
利用磁场与核自旋之间的交互作用以便形成二维或三维图像的MR成像方法在今天被广泛使用,特别是在医学诊断领域中,因为对于软组织的成像而言,它们在许多方面优于其他成像方法,即它们不需要致电离辐射,并且它们通常是无创的。
根据一般的MR方法,待检查的患者身体或一般对象被布置在强的均匀磁场B0中,磁场B0的方向同时限定了测量所依据的坐标系的轴,通常为z轴。
磁场根据所施加的磁场强度针对个体核自旋产生不同的能级,能够通过施加具有定义频率(所谓的拉莫尔频率或MR频率)的交变电磁场(RF场)来激励个体核自旋(自旋共振)。从宏观角度讲,个体核自旋的分布产生总体磁化,通过施加适当频率的电磁脉冲(RF脉冲)能够使总体磁化偏离平衡状态,同时磁场垂直于z轴延伸,使得磁化绕z轴进行进动运动。
借助接收RF天线能够检测磁化的任何变化,所述接收RF线圈被布置和定向在MR装置的检查体积之内,从而在垂直于z轴的方向测量磁化的变化。
为了在身体中实现空间分辨率,在均匀磁场上叠加沿三个主轴延伸的线性磁场梯度,导致自旋共振频率的线性空间相关性。然后在接收天线中拾取的信号包含可能与身体中的不同位置相关联的不同频率分量。
经由接收天线获得的信号数据对应于空间频率域,并且被称为k空间数据。k空间数据通常包括利用不同相位编码采集的多条线。通过收集若干样本对每条线进行数字化。利用傅里叶变换将一组k空间数据样本变换成MR图像。
磁共振成像(MRI)的一种特定应用是分别评估脉管渗透性和组织灌注、血液体积和流量。获得这种信息的一种方式是执行动态对比增强的MRI(DCE-MRI)和动态灵敏度增强的MRI(DSE-MRI)。DCE-MRI传递关于脉管渗透性、血管外细胞外空间和灌注的信息。DSE-MRI提供关于血液体积和流量的反馈。因此,通过对注射之后造影剂在组织和血液中随时间的浓度的药动学参数建模来导出生理标记。通过由注射药动学造影剂(通常为钆(Gd))诱发的纵向弛豫率R1(DCE)或R2*(DSE)的动态变化来确定浓度的动态变化。DCE最常使用稳态3D扰动梯度回波MR序列和短重复和回波时间(TR、TE),从信号强度的变化来量化R1。然而,忽略了灵敏度改变(R2*效应)对信号强度的影响,以简化该方法,这导致了不准确。另一方面,DSE通常采用单发(shot)EPI自旋回波或梯度回波序列以及长TR、TE来量化R2*,并忽略造影剂的药动学R1效应。
这里必须提到的是,R2*表示自旋系统的横向弛豫率,其包括由于磁场不均匀度的贡献,并且R1表示自旋系统的纵向弛豫。纵观整个说明书,“自旋系统的弛豫特性”被理解为相应的弛豫率或者是该弛豫率的倒数的弛豫时间。
动态氧气或二氧化碳增强的MRI(D(C)O2E-MRI)目前正日益获得关注,用以评估组织氧化和血管反应性。存在对癌症处置进行选择的重要参数。例如,辐射治疗的疗效取决于肿瘤的氧化水平。
该技术通常施加(多)梯度自旋序列,以用于在氧气或二氧化碳呼吸期间的R2*量化。R2*的变化反映血氧或血液流量或体积的变化。由溶解的氧和流量诱发的R1的同时变化是生理上感兴趣的,但幅度过低并且难以测量(耗时)以及因此在大多数情况下不予考虑。
在那些动态方法中R1和R2*的同时测量承载着改善生理性输出的准确性(DCE)或特异性(同时生理性测量)的潜质。
发明内容
从上文容易意识到,需要一种经改进的MR成像方法。因此,本发明的目的是使得MR图像能够同时精确确定对象的核自旋系统的纵向弛豫特性和横向弛豫特性。
根据本发明,提供了一种用于同时确定对象的核自旋系统的纵向弛豫特性和横向弛豫特性的磁共振成像方法。所述方法包括通过施加稳态梯度回波脉冲序列来激励核自旋,其中,所述脉冲序列包括一系列通常(但并非约束)小翻转角的RF激励脉冲,用于将核自旋磁化驱动至稳态水平。所述RF激励脉冲中的每个跟随有生成多个梯度回波的多梯度回波读出串。所述方法包括在施加稳态梯度回波脉冲序列的后续RF激励脉冲期间采集所述梯度回波的MR信号数据,并使用这些数据重建多幅图像。这一图像集—在k空间中或在真实图像空间中—表示在稳态多梯度回波序列的RF激励之后的信号衰减,被称为一个动态多梯度回波数据集。
此外,所述方法包括采集这些动态多梯度回波数据集中的至少两个以在至少两个时间点执行测量。此外,从利用稳态多梯度回波脉冲序列获得的第一动态多梯度回波数据集,针对第一时间点确定自旋系统的基线横向弛豫特性和自旋系统的基线平衡磁化。从利用稳态多梯度回波序列获得的第二动态多梯度回波数据集,获得针对第二时间点的平衡磁化。此外,从第二动态多梯度回波数据集,确定自旋系统的动态横向弛豫特性。此外,从平衡磁化和基线平衡磁化,确定自旋系统在第二时间点的动态纵向弛豫特性。如针对第二动态多梯度回波数据集的确定过程被用于所有可能的随后的动态多梯度回波数据集。
本发明的实施例具有这样的优点,两种弛豫特性,即,R1和R2*或者相应的弛豫时间T1和T2*能够同时或彼此独立地获得。常规动态R1测量结果的准确度得到改善并且加速了用于获得纵向和横向弛豫特性的总数据采集过程。
根据本发明的另一实施例,动态纵向弛豫特性的确定还包括将自旋系统的平衡磁化与基线平衡磁化相关联并从所述关联关系确定动态纵向弛豫特性。这样的优点在于,各种系统依赖关系和例如由于所使用的线圈系统的原因造成的相应伪影被消除,由此进一步改善了所得到的所确定的动态纵向弛豫特性的质量。
根据本发明的另一实施例,稳态多梯度回波脉冲序列是扰动稳态梯度回波脉冲序列,其中,所述序列包括在RF激励脉冲之前的扰动器梯度(梯度扰动),其中,所述扰动器梯度在每个多梯度回波读出串之后被施加。这具有进一步的优点,即任何剩余的横向磁化被去相,使得其不在后续的信号读出中产生任何信号。结果,再次改善了测量的质量。
根据本发明的另一实施例,稳态梯度回波脉冲序列是T1加权的快速场回波脉冲序列(T1FFE),其使用非恒定RF激励脉冲(RF-扰动)。然而,可以施加任何稳态多梯度回波脉冲序列,例如还可以施加涡(turbo)场回波序列(TFE)。所述T1FFE序列具有这样的优点,即自旋系统的纵向弛豫特性的确定被进一步增强和简化。
根据本发明的另一实施例,所述方法还包括通过所述系列的RF激励脉冲对动态横向和/或纵向弛豫特性的动态变化进行药动学建模。结果,可以从例如在钆(Gd)注射之后成相同的测量结果获得由此承载关于组织灌注的联合信息的动态R2*和R1曲线。或者,在氧呼吸期间对R2*和R1的同时测量传递关于血液和组织氧化的联合信息。
根据本发明的另一实施例,所述磁共振成像方法包括动态对比增强的或动态灵敏度增强的MRI方法。换言之,为了提高准确度,同时弛豫时间测量(relaxometry)方面支持同时以及因此有效的DCE-DSE生理测量。
根据本发明的另一实施例,所述方法还包括确定自旋系统的基线纵向弛豫特性,其描述在不存在稳态多梯度回波脉冲序列的情况下的自旋系统的纵向弛豫特性,其中,所述方法还包括从所述动态磁化数据确定所述动态纵向弛豫特性。
这一准确度能够在确定所述动态纵向弛豫特性的过程中通过额外地确定RF激励脉冲的翻转角以及通过采用所述确定的翻转角而被进一步改善。结果,还能够避免对关于RF激励脉冲的翻转角的理论假设。
本发明的方法能够有利地在大多数诊所中现有的MR装置中执行。为此,仅需要利用控制MR装置的计算机程序,使其执行本发明的上述方法步骤。所述计算机程序可以存在于数据载体上或者存在于数据网络上,以供下载从而在MR装置的控制单元中初始化。因此,本发明还涉及一种包括计算机可执行指令以执行上述方法的计算机程序产品。
此外,本发明涉及一种用于同时确定对象的核自旋系统的纵向弛豫特性和横向弛豫特性的磁共振成像设备。
所述设备包括磁共振成像扫描器,其用于采集磁共振图像数据,其中,所述设备还包括控制器,其适于控制扫描器执行以下操作:
-通过施加稳态梯度回波脉冲序列来激励核自旋,其中,所述脉冲序列包括一系列小翻转角的RF激励脉冲,用于将核自旋磁化驱动至稳态水平,其中,所述RF激励脉冲中的每个都跟随有多梯度回波读出串,其生成多梯度回波,
-在稳态梯度回波脉冲序列的后续RF激励脉冲的施加期间采集所述梯度回波的MR信号数据,其中,所述设备还包括数据重建或分析系统,其适于:
-从在第一时间点的所述MR信号数据重建第一动态多梯度回波数据集以及从在第二时间点的所述MR信号数据重建第二动态多梯度回波数据集,
-从所述第一动态多梯度回波数据集确定自旋系统的基线横向弛豫特性以及自旋系统的基线平衡磁化,
-从所述第二动态多梯度回波数据集确定自旋系统的平衡磁化,
-从所述第二动态多梯度回波数据集确定自旋系统的动态横向弛豫特性,
-从自旋系统的基线平衡磁化以及自旋系统的平衡磁化来确定自旋系统的动态纵向弛豫特性。
附图说明
附图公开了本发明的优选实施例。然而,应当理解,所述附图仅仅是设计用于图示说明的目的,而并非作为对本发明的限制。在附图中:
图1示出了执行本发明的方法的MR装置;
图2图示了同时确定自旋系统的纵向和横向弛豫特性的原理;
图3是用于同时确定自旋系统的横向和纵向弛豫特性的动态处理步骤的流程图;
图4描绘了健康志愿者对卡波金(Carbogen)(95%)的R2*和R1反应图。
具体实施方式
参考图1,示出了MR成像系统1。该系统包括超导或电阻式主磁体线圈2,使得沿穿过检查体积的z轴生成大致均匀的时间上恒定的主磁场B0
磁共振生成操控系统施加一系列RF脉冲和切变的磁场梯度以逆转或激励核磁自旋、感生磁共振、重聚焦磁共振、操纵磁共振、对磁共振进行空间或其他方式的编码、使自旋饱和等等,以执行MR程序。
更具体而言,梯度脉冲放大器3施加电流脉冲至沿检查体积的x、y和z轴的全身梯度线圈4、5和6中的选定的一个。RF发射器7经由发送/接收开关8发射RF脉冲或脉冲包至RF天线9,以将RF脉冲发射至检查体积中。典型的MR成像序列包括彼此同时发生的短持续时间的RF脉冲序列包,并且任意施加的磁场梯度实现对核磁共振的选定的操控。RF脉冲被用于使共振饱和、激励共振、逆转磁化、重聚焦共振或者操纵共振以及选择定位在检查体积中的身体10的部分。MR信号还可以由RF天线9来拾取。
为了例如借助并行成像生成身体或者通常为对象10的有限区域的MR图像,邻近所述区域放置一组局部阵列RF线圈11、12和13以进行成像。所述阵列线圈11、12和13能够被用于接收经由RF天线影响的RF发射感生的MR信号。然而,其还可能使用阵列线圈11、12和13来发射RF信号到检查体积。
所得的MR信号由RF天线9和/或由RF线圈11、12和13的阵列来拾取,并由优选包括前置放大器(未示出)的接收器14来解调。接收器14经由发送/接收开关8与RF线圈9、11、12和13连接。
主计算机15控制梯度脉冲放大器3和发射器7以生成多种成像序列中的任意一种,诸如回波平面成像(EPI)、回波体积成像、梯度和自旋回波成像、快速自旋回波成像等。这包括生成如上文所述的成像序列,以确定对象10的核自旋系统的纵向弛豫特性和横向弛豫特性。
针对选定的序列,接收器14以跟随每个RF激励脉冲的快速连续方式接收单条或多条MR数据线。数据采集系统16对所接收的信号执行模数转换并将每条MR数据线转换成适于进行进一步处理的数字格式。在现代的MR装置中,数据采集系统16是分离的计算机,其专用于采集原始图像数据。
最后,所述数字原始图像数据通过重建处理器17被重建为图像表示,其应用傅里叶或其他适当的重建算法,如关于对象10的核自旋系统的纵向弛豫特性和横向弛豫特性的确定所描述的。MR图像可以表示穿过患者的平面切片、平行平面切片的阵列、三维体积等。所述图像然后被存储在图像存储器中,其中,可以访问所述存储器以将切片或图像表示的其他部分转换成适当的格式,以供例如经由视频监视器18可视化,视频监视器18提供所得的MR图像的人可读的显示。
图2a图示了现有技术的横向弛豫特性的状态的原理,尤其是R2*测量。R2*量化通常基于在RF激励之后的不同回波时间处对多梯度回波图像的采集。如图2a中所示,施加具有翻转角α的RF激励脉冲200,随后是一组读出梯度202以及多个图像采集点214。如图2所示,磁化M从初始磁化值M0随时间下降。R2*被量化为所有回波的信号强度S的指数衰减的弛豫率常量:
S=M0·exp(-TE·R2*)=M0·E2    等式1
图2中的灰色阴影块指示多梯度回波序列的各个扫描,其中,由时间TR将所述多RF激励脉冲间隔开。
图2b图示了现有技术的借助稳态扰动梯度回波序列的动态R1量化的原理。所述序列包括一组小翻转角(扰动)(α)的RF脉冲204,其被以短的重复时间TR施加到自旋系统。在每个个体RF脉冲204中,施加读出梯度206和扰动器梯度,以使任何剩余的横向磁化去相。在存在例如切片选择梯度的情况下施加所述小翻转角的RF激励脉冲204,并且在在存在读出梯度206的数据采集窗口期间采集相应的MR信号数据。
为了成像的目的,所述脉冲序列可以以所述间隔TR来重复,并且所述相位编码梯度可以通过一组值被分阶,以采集k空间的完整采样,从k空间的完整采样可以重建图像。因此扫描时间是同相位编码的数量确定的。RF和梯度扰动的目的是,通常TR小于横向弛豫时间,使得横向磁化保持在每次个体扫描210的末端(每次扫描210都在相应的读出梯度206中包括RF脉冲204)。借助扰动器梯度,剩余的横向磁化被去相,使得其不影响读出MR信号。
非常快速的R1量化的一种方法是使用稳态(TR<1/R1)扰动梯度回波序列的取决于R1的平衡磁化。在下文中,所采集的信号由S表示:
其中,E1=exp(-TR·R1),E2=exp(-TE·R2*)
等式2a
通常,使用短回波时间TE<1\R2*并且取决于R2*的指数E2为负,并且因此——如上所述——R2*的影响也改变:
对于TE<1/R2* S &ap; M 0 = &kappa; &CenterDot; sin &alpha; &CenterDot; 1 - E 1 1 - cos &alpha; &CenterDot; E 1 , 等式2b
然后基于在时间点t的信号St相对于基线信号St=0(例如动态序列的首次扫描)的变化来量化R1的动态变化(即R1t):
S t S t = 0 = 1 - E 1 t 1 - cos &alpha; &CenterDot; E 1 t 1 - cos &alpha; &CenterDot; E 1 0 1 - E 1 0 等式3
E10和E1t分别是固有(exp(-TR·R10))指数项和动态(exp(-TR·R1t))指数项。如上所述,扰动(通过在每个RF脉冲之间利用RF扰动和压碎(crusher)进行相位不规则化)是重要的,以消除任何先前的横向磁化对跟随RF脉冲的信号的贡献。
从等式3以及对翻转角α和R10的进一步(例如理论)假设,能够估计动态纵向弛豫率R1t
图2c图示了同时确定核自旋系统的纵向和横向弛豫特性的原理。本发明提出组合图2a和图2b中所示的两种原理,而所述组合能够被理解为动态R1量化原理扩展了多梯度回波读出,或者被理解为动态稳态R2*量化原理扩展了R10量化。
在图2c中所示的这两种原理的组合中。施加一系列RF激励脉冲204,以将核自旋磁化驱动至稳态水平。RF激励脉冲204中的每个跟随有多梯度回波读出串202,其中,这种读出串生成多个梯度回波,其利用多个图像采集步骤被采集为MR信号数据,由箭头214表示。此外,借助扰动器梯度脉冲216,在施加每个个体多梯度回波读出串之后,扰动剩余的横向磁化。在图2c中,个体扫描,每次扫描都包括一系列RF激励脉冲204和读出梯度202,由块214表示。通过连续重复块214(即扫描),能够获得R1和R2*的量化以及动态R1和R2*变化。
任选地,如由参考标记212所表示的,可以确定自旋系统的基线纵向弛豫特性,其描述在不存在稳态梯度回波脉冲序列的情况下的自旋系统的纵向弛豫特性,得到R10。例如,可以在任何造影剂被注射或者氧气、(D(C)O2E-MRI)被吸入之前确定自旋系统的所述基线纵向弛豫特性。所述稳态序列自身对R1没有影响。相应的扫描可以在执行扫描214之前执行或者在执行扫描214之后执行,在任一种情况下,当自旋系统被完全返回到平衡值。
同样任选地,扫描212可以包括确定RF激励脉冲的精确翻转角,其中,之后,所述翻转角在成扫描214确定纵向弛豫特性的过程中被采用。
图3是图示了用于同时确定自旋系统的横向和纵向弛豫特性的动态处理步骤的流程图。每个个体扫描214的多梯度回波信号传递所有回波的信号强度的指数衰减,如相对于等式1所论述的。结果,对于每个动态步骤(对于每次扫描),能够在每个动态时间点例如通过使用指数拟合或数值方法解方程1来获得R2* t(动态横向弛豫率)和M0t(动态磁化)。此外,从在首次扫描中采集的MR信号数据,可以在步骤300中获得自旋系统(R2* 0)的静态横向弛豫特性和自旋系统(M00)的静态平衡磁化。
在下文中,讨论步骤320:根据针对扰动梯度回波短TR序列的稳态信号方程,M0是R1的函数,但独立于R2*。M0t和M00的相除还消除了系统相关的常量κ和正弦α项。
在步骤304中,然后量化R1t的动态变化,其中,通过针对E1t在步骤302中表示的方程进行求解来获得在步骤304中表示的方程。在步骤304中剩余的变量是E10、翻转角α和TR。这些值可以或者已经从额外的在步骤308中(在图2c中的扫描)的先前的测量获得或者是序列特异性常量(TR、α)。
最后,在步骤306中,从E1t的算法中获得R11
应当认识到,在扫描212(步骤308)中的额外的(任选的)B1测量能够通过包含真实的翻转角确定来进一步改善该方法的准确度。由于高的扫描效率,可变翻转角方法常常被用于R10量化,但Look-Locker、逆转恢复TSE或同时B1和R10映射技术当然也是适用的。
本发明的特征是组合R1和R2*量化方法,使得对于例如在DCE-MRI和DSE-MRI中所使用的标准动态R1和R2*测量不同,同时并且相对于彼此独立地传递两种弛豫率值。常规动态R1测量的准确度得到改善,因为在方程2中,E2不再需要为负,而是从多梯度回波信号中导出。在提高准确度之后,同时弛豫时间测量方面支持同时,以及因此有效的DCE-DSE生理测量。从在Gd注射之后以及由此承载关于组织灌注、流量和脉管渗透性的联合信息的相同测量结果获得动态R2*和R1曲线。
本发明可能在D(C)O2E-MRI实验中更为重要。富含氧气或CO2的气体的吸入分别改变了血液和组织氧含量以及血流量和体积。同时的R1和R2*序列被认为提高了特异性并避免了“假阴性”(参见背景技术部分)。R2*对在组织和体液中溶解的分子氧敏感,取决与所述序列,还对流量的变化敏感。
图4中展示了基本原理。该图描绘了健康志愿者对卡波金(95%)的R2*(左)和R1(右)响应图。血氧的增加导致R2*随着时间的下降,以及在血浆和组织中溶解的药动学分子氧的量的升高,以及血流速度的升高导致R1随时间的升高。尽管幅度非常低(比R2*的效果地100倍),但能够在数秒内对R1的动态变化进行量化。在这一2D范例中,时间分辨率是2.1s/帧并且平面内空间分辨率小于2mm。
在组织和血浆中流量和分子氧的增加导致急剧的R1变化,其能够利用在3T时2.1s的时间分辨率进行量化。血氧水平的升高导致更高的R2*变化(≈100倍),其借助上文论述的方法被同时量化。
这种方法能够进一步支持对造影剂的递送和分布的定量检测,其利用R1相关的和灵敏度机制这两者。

Claims (17)

1.一种磁共振成像方法,其用于同时确定对象的核自旋系统的纵向弛豫特性和横向弛豫特性,所述方法包括通过施加稳态梯度回波脉冲序列(210)来激励所述核自旋,其中,所述脉冲序列包括一系列RF激励脉冲(204),用于将核自旋磁化驱动至稳态水平,其中,所述RF激励脉冲(204)中的每个跟随有生成多个梯度回波的多梯度回波读出串(202),其中,所述方法包括在施加所述稳态梯度回波脉冲序列的后续RF激励脉冲(204)期间采集所述梯度回波的MR信号数据,并从在第一时间点的所述MR信号数据重建第一动态多梯度回波数据集以及从在第二时间点的所述MR信号数据重建第二动态多梯度回波数据集,其中,所述方法还包括:
-从所述第一动态多梯度回波数据集确定(300)所述核自旋系统的基线横向弛豫特性以及所述核自旋系统的基线平衡磁化,
-从所述第二动态多梯度回波数据集确定所述核自旋系统的平衡磁化,
-从所述第二动态多梯度回波数据集确定所述核自旋系统的动态横向弛豫特性,
-从所述核自旋系统的所述基线平衡磁化和所述核自旋系统的所述平衡磁化确定所述核自旋系统的动态纵向弛豫特性。
2.根据权利要求1所述的方法,其中,所述动态纵向弛豫特性的所述确定还包括将所述核自旋系统的所述平衡磁化与所述核自旋系统的所述基线平衡磁化关联(302),并从所述关联确定所述动态纵向弛豫特性。
3.根据权利要求1所述的方法,其中,所述稳态梯度回波脉冲序列是经扰动的稳态梯度回波脉冲序列,所述序列包括在所述RF激励脉冲(204)之前的扰动器梯度(216),其中,所述扰动器梯度(216)在每个多梯度回波读出串之后被施加。
4.根据权利要求1所述的方法,其中,所述稳态梯度回波脉冲序列是T1加权的快速场回波脉冲序列。
5.根据权利要求1所述的方法,还包括通过所述系列RF激励脉冲(204)对所述动态横向弛豫特性和/或所述动态纵向弛豫特性的动态变化进行药动学建模。
6.根据权利要求5所述的方法,其中,所述磁共振成像方法包括动态对比增强的MRI方法或动态灵敏度增强的MRI方法。
7.根据权利要求1所述的方法,还包括确定所述核自旋系统的基线纵向弛豫特性,其描述在不存在所述稳态梯度回波脉冲序列的情况下的所述核自旋系统的所述纵向弛豫特性,其中,所述方法还包括从所述基线纵向弛豫特性确定所述动态纵向弛豫特性。
8.根据权利要求1所述的方法,还包括确定(308)所述RF激励脉冲(204)的翻转角,其中,所述翻转角在所述动态纵向弛豫特性的所述确定中被采用。
9.一种磁共振成像装置,其用于同时确定对象的核自旋系统的纵向弛豫特性和横向弛豫特性,所述装置包括用于通过施加稳态梯度回波脉冲序列(210)来激励所述核自旋的模块,其中,所述脉冲序列包括一系列RF激励脉冲(204),用于将核自旋磁化驱动至稳态水平,其中,所述RF激励脉冲(204)中的每个跟随有生成多个梯度回波的多梯度回波读出串(202),其中,所述装置包括用于在施加所述稳态梯度回波脉冲序列的后续RF激励脉冲(204)期间采集所述梯度回波的MR信号数据,并从在第一时间点的所述MR信号数据重建第一动态多梯度回波数据集以及从在第二时间点的所述MR信号数据重建第二动态多梯度回波数据集的模块,其中,所述装置还包括:
-用于从所述第一动态多梯度回波数据集确定(300)所述核自旋系统的基线横向弛豫特性以及所述核自旋系统的基线平衡磁化的模块,
-用于从所述第二动态多梯度回波数据集确定所述核自旋系统的平衡磁化的模块,
-用于从所述第二动态多梯度回波数据集确定所述核自旋系统的动态横向弛豫特性的模块,
-用于从所述核自旋系统的所述基线平衡磁化和所述核自旋系统的所述平衡磁化确定所述核自旋系统的动态纵向弛豫特性的模块。
10.根据权利要求9所述的装置,其中,所述动态纵向弛豫特性的所述确定还包括将所述核自旋系统的所述平衡磁化与所述核自旋系统的所述基线平衡磁化关联(302),并从所述关联确定所述动态纵向弛豫特性。
11.根据权利要求9所述的装置,其中,所述稳态梯度回波脉冲序列是经扰动的稳态梯度回波脉冲序列,所述序列包括在所述RF激励脉冲(204)之前的扰动器梯度(216),其中,所述扰动器梯度(216)在每个多梯度回波读出串之后被施加。
12.根据权利要求9所述的装置,其中,所述稳态梯度回波脉冲序列是T1加权的快速场回波脉冲序列。
13.根据权利要求9所述的装置,还包括用于通过所述系列RF激励脉冲(204)对所述动态横向弛豫特性和/或所述动态纵向弛豫特性的动态变化进行药动学建模的模块。
14.根据权利要求13所述的装置,其中,所述磁共振成像装置包括动态对比增强的MRI装置或动态灵敏度增强的MRI装置。
15.根据权利要求9所述的装置,还包括用于确定所述核自旋系统的基线纵向弛豫特性的模块,所述基线纵向弛豫特性描述在不存在所述稳态梯度回波脉冲序列的情况下的所述核自旋系统的所述纵向弛豫特性,其中,所述装置还包括用于从所述基线纵向弛豫特性确定所述动态纵向弛豫特性的模块。
16.根据权利要求9所述的装置,还包括用于确定(308)所述RF激励脉冲(204)的翻转角的模块,其中,所述翻转角在所述动态纵向弛豫特性的所述确定中被采用。
17.一种磁共振成像设备(1),其用于同时确定对象的核自旋系统的纵向弛豫特性和横向弛豫特性,其中,所述设备包括磁共振成像扫描器,其用于采集磁共振图像数据,其中,所述设备还包括适于控制扫描器执行如下操作的控制器(15):
-通过施加稳态梯度回波脉冲序列来激励所述核自旋,其中,所述脉冲序列包括一系列小翻转角的RF激励脉冲(204),用于将核自旋磁化驱动至稳态水平,其中,所述RF激励脉冲(204)中的每个跟随有生成多个梯度回波的多梯度回波读出串(202),
-在施加所述稳态梯度回波脉冲序列的后续RF激励脉冲(204)期间采集所述梯度回波的MR信号数据,
其中,所述设备还包括数据重建系统(17),其适于执行以下操作:
-从在第一时间点的所述MR信号数据重建第一动态多梯度回波数据集,以及从在第二时间点的所述MR信号数据重建第二动态多梯度回波数据集,
-从所述第一动态多梯度回波数据集确定所述核自旋系统的基线横向弛豫特性以及所述核自旋系统的基线平衡磁化,
-从所述第二动态多梯度回波数据集确定所述核自旋系统的动态横向弛豫特性,
-从所述第二动态多梯度回波数据集确定所述核自旋系统的平衡磁化,
-从所述核自旋系统的所述基线平衡磁化和所述核自旋系统的所述平衡磁化确定所述核自旋系统的动态纵向弛豫特性。
CN201180014334.9A 2010-03-18 2011-03-10 核自旋系统的纵向和横向弛豫时间的同时和动态确定 Expired - Fee Related CN102859386B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US31506210P 2010-03-18 2010-03-18
US61/315,062 2010-03-18
PCT/IB2011/051001 WO2011114264A1 (en) 2010-03-18 2011-03-10 Simultaneous and dynamic determination of longitudinal and transversal relaxation times of a nuclear spin system

Publications (2)

Publication Number Publication Date
CN102859386A CN102859386A (zh) 2013-01-02
CN102859386B true CN102859386B (zh) 2015-08-19

Family

ID=44059096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180014334.9A Expired - Fee Related CN102859386B (zh) 2010-03-18 2011-03-10 核自旋系统的纵向和横向弛豫时间的同时和动态确定

Country Status (5)

Country Link
US (1) US9588205B2 (zh)
EP (1) EP2548039A1 (zh)
JP (1) JP6050126B2 (zh)
CN (1) CN102859386B (zh)
WO (1) WO2011114264A1 (zh)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102859386B (zh) * 2010-03-18 2015-08-19 皇家飞利浦电子股份有限公司 核自旋系统的纵向和横向弛豫时间的同时和动态确定
KR20130025190A (ko) * 2011-09-01 2013-03-11 롬엔드하스전자재료코리아유한회사 신규한 유기 전자재료용 화합물 및 이를 채용하고 있는 유기 전계 발광 소자
US9618596B2 (en) * 2012-07-02 2017-04-11 Syntheticmr Ab Methods and systems for improved magnetic resonance acquisition using a single acquisition type
US10073156B2 (en) 2012-07-02 2018-09-11 Syntheticmr Ab Methods and systems for improved magnetic resonance acquisition
US10203387B2 (en) * 2013-06-06 2019-02-12 Koninklijke Philips N.V. MR imaging with enhanced susceptibility contrast
JP6266325B2 (ja) * 2013-12-06 2018-01-24 東芝メディカルシステムズ株式会社 画像解析装置、画像解析プログラム、及び、磁気共鳴イメージング装置
CN104808161B (zh) * 2014-01-27 2018-11-02 西门子(深圳)磁共振有限公司 磁共振成像方法和装置
JP6783864B2 (ja) 2016-01-22 2020-11-11 シナプティヴ メディカル (バルバドス) インコーポレイテッドSynaptive Medical (Barbados) Inc. 磁気共鳴画像法を使用した磁場に依存するリラクソメトリーのためのシステムおよび方法
CN107219484B (zh) * 2016-03-22 2023-09-22 西门子(深圳)磁共振有限公司 多层并行激发弥散成像的激发翻转角度确定方法和装置
US10290103B2 (en) 2016-05-26 2019-05-14 Synthetic Mr Ab Method, device and non-transitory digital storage medium for non-aqueous tissue volume estimation
EP3336570A1 (en) * 2016-12-15 2018-06-20 Universität Heidelberg Magnetic resonance fingerprinting (mrf) using echo-planar imaging with spoiling
CN109342981B (zh) 2018-09-20 2021-10-12 上海东软医疗科技有限公司 一种基于epi序列的动态成像方法和装置
CN111568420B (zh) * 2019-02-15 2021-05-14 深圳先进技术研究院 一种旋转坐标内弛豫时间确定方法、装置及存储介质
CN110133553B (zh) * 2019-05-10 2020-06-05 浙江大学 一种超短回波时间磁共振指纹弛豫时间测量方法
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
CN111025210B (zh) * 2019-12-20 2022-02-01 东软医疗系统股份有限公司 磁共振成像方法、装置、电子设备、存储介质
US11982728B2 (en) 2022-02-16 2024-05-14 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for magnetic resonance imaging
CN114690101B (zh) * 2022-03-30 2024-08-02 西门子数字医疗科技(上海)有限公司 磁共振成像中的参数定量检测方法、装置及磁共振扫描仪
CN116106806B (zh) * 2023-04-07 2023-08-29 深圳市联影高端医疗装备创新研究院 多核成像参数确定方法、装置和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1039484A (zh) * 1988-07-08 1990-02-07 菲利浦光灯制造公司 磁共振成象方法和设备
US6291996B1 (en) * 1997-08-13 2001-09-18 Btg International Limited Apparatus for and method of determining values of relaxation parameters
EP1136836A2 (en) * 2000-03-21 2001-09-26 Dmitriy A. Yablonskiy T2 contrast in magnetic resonance imaging with gradient echoes
CN1636508A (zh) * 2003-12-22 2005-07-13 Ge医疗系统环球技术有限公司 磁共振成像设备和磁共振成像方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04288141A (ja) 1991-03-18 1992-10-13 Hitachi Ltd 核磁気共鳴イメージング法
US5245282A (en) * 1991-06-28 1993-09-14 University Of Virginia Alumni Patents Foundation Three-dimensional magnetic resonance imaging
JPH06237914A (ja) 1993-02-16 1994-08-30 Yokogawa Medical Syst Ltd Mri装置
JPH06335471A (ja) 1993-05-28 1994-12-06 Hitachi Medical Corp Mri装置
US6205349B1 (en) * 1998-09-29 2001-03-20 Siemens Medical Systems, Inc. Differentiating normal living myocardial tissue, injured living myocardial tissue, and infarcted myocardial tissue in vivo using magnetic resonance imaging
US7081751B2 (en) * 2002-12-31 2006-07-25 The Uab Research Foundation Systems and methods for estimating properties of a sample
RU2434645C2 (ru) * 2006-03-31 2011-11-27 Конинклейке Филипс Электроникс, Н.В. Системы и методы измерения клеток, использующие ультракороткую т2*-релаксометрию
US8143889B2 (en) * 2009-02-24 2012-03-27 University Of Utah Research Foundation Simultaneous acquisitions of spin- and stimulated-echo planar imaging
CN102859386B (zh) * 2010-03-18 2015-08-19 皇家飞利浦电子股份有限公司 核自旋系统的纵向和横向弛豫时间的同时和动态确定
EP3133994A4 (en) * 2014-04-24 2018-06-20 Oregon Health & Science University Contrast reagent leakage correction in dynamic susceptibility contrast magnetic resonance imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1039484A (zh) * 1988-07-08 1990-02-07 菲利浦光灯制造公司 磁共振成象方法和设备
US6291996B1 (en) * 1997-08-13 2001-09-18 Btg International Limited Apparatus for and method of determining values of relaxation parameters
EP1136836A2 (en) * 2000-03-21 2001-09-26 Dmitriy A. Yablonskiy T2 contrast in magnetic resonance imaging with gradient echoes
CN1636508A (zh) * 2003-12-22 2005-07-13 Ge医疗系统环球技术有限公司 磁共振成像设备和磁共振成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Dynamic T1 Estimation of Brain Tumors Using Double-Echo Dynamic MR Imaging;Yoshiyuki Ishimori et al.;《JOURNAL OF MAGNETIC RESONANCE IMAGING》;20031231;第18卷;第113-120页 *
Novel Method for Rapid, Simultaneous T1, T2*, and Proton Density Quantification;J.B.M. Warntjes et al.;《Magnetic Resonance in Medicine》;20071231;第57卷;第528-537页 *
Separation of T1 and T2* shortening effects using simultaneous T1/T2* gradient dual echo pulse sequence;E.Kim et al.;《Proc. Intl. Soc. Mag. Reson. Med.》;20030710;第1246页 *
Simultaneous dynamic T1 and T2* measurement for AIF assessment combined with DCE MRI in a mouse tumor model;Melanie Heilmann et al.;《Magn Reson Mater Phy》;20071011;第20卷(第4期);第193-203页 *

Also Published As

Publication number Publication date
EP2548039A1 (en) 2013-01-23
JP2013521922A (ja) 2013-06-13
CN102859386A (zh) 2013-01-02
US9588205B2 (en) 2017-03-07
US20120326721A1 (en) 2012-12-27
WO2011114264A1 (en) 2011-09-22
JP6050126B2 (ja) 2016-12-21

Similar Documents

Publication Publication Date Title
CN102859386B (zh) 核自旋系统的纵向和横向弛豫时间的同时和动态确定
CN103238082B (zh) 使用多点Dixon技术和低分辨率校准的MR成像
CN102959388B (zh) 利用压缩感测重建的动态对比度增强mr成像
CN111505553B (zh) 磁共振成像系统和方法
US8558547B2 (en) System and method for magnetic resonance radio-frequency field mapping
CN102859387B (zh) 采用谱模型的化学物类的磁共振成像
CN105814449B (zh) 具有水/脂肪分离的零回波时间mr成像
US6064203A (en) Method and apparatus for determining or imaging longitudinal spin relaxation time or producing images which substantially reflect longitudinal spin relaxation time contrast
CN103649765A (zh) 具有b1绘制的mr成像
CN103229069A (zh) 使用多点狄克逊技术的mr 成像
CN103376432A (zh) 快速地位置分辨地确定检查区域内磁共振弛豫参数的方法
CN103857331A (zh) 用于量化人体中肝脏脂肪的方法和系统
CN102713658B (zh) 磁化率梯度绘图
CN101354433A (zh) 记录和处理时间上相继的图像数据组的序列的方法及设备
CN102914754A (zh) 用于产生磁共振照片的方法
CN106574954A (zh) 针对epi的具有奈奎斯特伪影校正的并行mr成像
CN104204838A (zh) 具有Dixon类型水/脂肪分离和关于主磁场的不均匀性的先验知识的MRI
CN105103001A (zh) 使用高snr同相图像和较低snr至少部分地异相图像的dixon型水/脂肪分离mri
Wright et al. Comparison of two ultrashort echo time sequences for the quantification of T1 within phantom and human Achilles tendon at 3 T
CN107076818B (zh) 零回波时间mr成像
CN103140167A (zh) 化学物类的磁共振成像
Brizi et al. Bone volume–to–total volume ratio measured in trabecular bone by single‐sided NMR devices
Zaitsev et al. Error reduction and parameter optimization of the TAPIR method for fast T1 mapping
US20140316251A1 (en) Mri method for investigating interstitial fluid velocity in a tissue using a nulling preparation
US11415652B2 (en) Multi-dimensional spectroscopic NMR and MRI using marginal distributions

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150819

Termination date: 20180310

CF01 Termination of patent right due to non-payment of annual fee