当前位置: 首页 > news >正文

磁共振成像原理(理论)30:基本自旋回波成像 (Basic Spin-Echo Imaging)

本专题主要参考《Principles of Magnetic Resonance Imaging A Signal Processing Perspective 》-Sec 7.4

饱和恢复自旋回波序列

在磁共振成像中,普通的饱和恢复和反转恢复序列产生的是T1加权的自由感应衰减信号。然而,为了实现k空间的对称覆盖,通常采用回波信号进行数据采集。本节将详细讨论将自旋回波技术引入这些成像序列后所产生的图像对比度特性。

image.png

上图展示了饱和恢复自旋回波序列的时序图。该序列结合了饱和恢复的磁化准备机制和自旋回波的信号采集方式。

第1次饱和恢复:0≤t<TR0 \leq t< T_R0t<TR

t=0−t = 0_-t=0时刻,即施加90°射频脉冲 (假设是x′x'x方向的)前,系统处于热平衡状态,纵向和横向磁化强度分别是:
Mz′(0−)=Mz0Mx′y′(0−)=0() \begin{aligned} M_z'(0_-) &= M_z^0 \\ M_{x'y'}(0_-) &=0 \end{aligned}\tag{} Mz(0)Mxy(0)=Mz0=0()

t=0+t = 0_+t=0+时刻,即施加90°射频脉冲后,纵向和横向磁化强度分别为:
Mz′(0+)=0Mx′y′(0+)=Mz0() \begin{aligned} M_z'(0_+) &= 0 \\ M_{x'y'}(0_+) &= M_z^0 \end{aligned}\tag{} Mz(0+)Mxy(0+)=0=Mz0()

根据弛豫方程(式3.122)
{Mx′y′(t)=Mx′y′(0+)e−t/T2Mz′(t)=Mz0(1−e−t/T1)+Mz′(0+)e−t/T1(3.122) \left\{ \begin{array}{l} M_{x'y'}(t) = M_{x'y'}(0_+) e^{-t/T_2} \\ M_{z'}(t) = M_z^0 (1 - e^{-t/T_1}) + M_{z'}(0_+) e^{-t/T_1} \end{array} \right. \tag {3.122} {Mxy(t)=Mxy(0+)et/T2Mz(t)=Mz0(1et/T1)+Mz(0+)et/T1(3.122)

0≤t<TI0 \leq t < T_I0t<TI 时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−t/T1)(0<t<TI)Mx′y′(t)=Mz0e−t/T2(0<t<TI)() \begin{aligned} M_{z'}(t) &=M_z^0 \left(1 - e^{-t/T_1}\right) &\quad (0 < t < T_I) \\ M_{x'y'}(t) &= M_z^0e^{-t/T_2} &\quad (0 < t < T_I) \end{aligned} \tag{} Mz(t)Mxy(t)=Mz0(1et/T1)=Mz0et/T2(0<t<TI)(0<t<TI)()

t=TI−t = T_{I-}t=TI时刻,即180°脉冲施加前,横向和纵向磁化强度分别为:
Mz′(TI−)=Mz0(1−e−TI/T1)Mx′y′(TI−)=Mz0e−TI/T2() \begin{aligned} M_{z'}(T_{I-}) &= M_z^0 \left(1 - e^{-T_I/T_1} \right) &\quad\\ M_{x'y'}(T_{I-}) &= M_z^0e^{-T_I/T_2} &\quad \end{aligned} \tag{} Mz(TI)Mxy(TI)=Mz0(1eTI/T1)=Mz0eTI/T2()

理想180°射频脉冲给横向和纵向磁化带来的影响如下:
Mz′→180∘x′−Mz′() M_{z'} \xrightarrow{{180^\circ}_{x'}} -M_{z'} \tag{} Mz180xMz()
Mx′y′→180∘x′−Mx′y′() M_{x'y'} \xrightarrow{{180^\circ}_{x'}} -M_{x'y'} \tag{} Mxy180xMxy()

t=TI+t = T_{I+}t=TI+时刻,即180°脉冲施加后,横向和纵向磁化强度分别为:
Mz′(TI+)=−Mz′(TI−)=−Mz0(1−e−TI/T1)Mx′y′(TI+)=−Mx′y′(TI−)=−Mz0e−TI/T2() \begin{aligned} M_{z'}({T_I}+) &= -{M_{z'}({T_I}_-)} = -M_z^0 \left(1 - e^{-T_I/T_1} \right)\\ M_{x'y'}({T_I}+) &= -{M_{x'y'}({T_I}_-)} = -M_z^0e^{-T_I/T_2} \end{aligned} \tag{} Mz(TI+)Mxy(TI+)=Mz(TI)=Mz0(1eTI/T1)=Mxy(TI)=Mz0eTI/T2()

t>TI+t > T_{I+}t>TI+时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−(t−TI)/T1)−Mz0(1−e−TI/T1)e−(t−TI)/T1=Mz0(1−e−(t−TI)/T1)−Mz0(e−(t−TI)/T1−e−t/T1)=Mz0(1−2e−(t−TI)/T1)+Mz0e−t/T1=Mz0(1−2e−(t−TI)/T1+e−t/T1)Mx′y′(t)=−Mz0e−t/T2() \begin{aligned} M_{z'}(t) &= M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(1 - e^{-T_I/T_1} \right) e^{-(t-T_I)/T_1} &\quad\\ &=M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(e^{-(t-T_I)/T_1} - e^{-t/T_1} \right) \\ &=M_z^0 (1 - 2e^{-(t-T_I)/T_1}) + M_z^0 e^{-t/T_1} \\ &=M_z^0 (1 - 2e^{-(t-T_I)/T_1}+e^{-t/T_1}) \\ \\ M_{x'y'}(t) &= -M_z^0e^{-t/T_2} &\quad \end{aligned} \tag{} Mz(t)Mxy(t)=Mz0(1e(tTI)/T1)Mz0(1eTI/T1)e(tTI)/T1=Mz0(1e(tTI)/T1)Mz0(e(tTI)/T1et/T1)=Mz0(12e(tTI)/T1)+Mz0et/T1=Mz0(12e(tTI)/T1+et/T1)=Mz0et/T2()

第2次饱和恢复:TR≤t<2TRT_R \leq t< 2T_RTRt<2TR

t=TR−t = {T_R}_-t=TR时刻,横向和纵向磁化强度分别为:
Mz′(TR−)=Mz0(1−2e−(TR−TI)/T1−e−TR/T1)Mx′y′(TR−)=−Mz0e−TR/T2() \begin{aligned} M_{z'}({T_R}_-) &=M_z^0 (1 - 2e^{-(T_R-T_I)/T_1}-e^{-T_R/T_1}) \\ \\ M_{x'y'}({T_R}_-) &= -M_z^0e^{-T_R/T_2} &\quad \end{aligned} \tag{} Mz(TR)Mxy(TR)=Mz0(12e(TRTI)/T1eTR/T1)=Mz0eTR/T2()

这里将TE=2TIT_E = 2 T_ITE=2TI代入上式,可得在时间t=TR−t = {T_R}_-t=TR时刻,横向和纵向磁化强度分别为:
Mz′(TR−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)Mx′y′(TR−)=−Mz0e−TR/T2() \begin{aligned} M_{z'}({T_R}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\ \\ M_{x'y'}({T_R}_-) &= -M_z^0e^{-T_R/T_2} &\quad \end{aligned} \tag{} Mz(TR)Mxy(TR)=Mz0(12e(TRTE/2)/T1eTR/T1)=Mz0eTR/T2()
为了符号方便阅读,这里将TR−→0−{T_R}_- \to 0_-TR0,可得:
Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)Mx′y′(0−)=−Mz0e−TR/T2() \begin{aligned} M_{z'}({0}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\ \\ M_{x'y'}({0}_-) &= -M_z^0e^{-T_R/T_2} &\quad \end{aligned} \tag{} Mz(0)Mxy(0)=Mz0(12e(TRTE/2)/T1eTR/T1)=Mz0eTR/T2()

由于T2<TE<<TRT_2 < T_E<<T_RT2<TE<<TR,所以上面的式子中横向磁化可以忽略不计。因此

t=0−t = 0_-t=0时刻,即施加90°射频脉冲前,纵向和横向磁化强度分别为
Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)Mx′y′(0−)~0() \begin{aligned} M_{z'}({0}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\ \\ M_{x'y'}({0}_-) &~ 0 &\quad \end{aligned} \tag{} Mz(0)Mxy(0)=Mz0(12e(TRTE/2)/T1eTR/T1)0()

t=0+t = 0_+t=0+时刻,即施加90°射频脉冲后,纵向和横向磁化强度分别为:
Mz′(0+)=0Mx′y′(0+)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)() \begin{aligned} M_z'(0_+) &= 0\\ M_{x'y'}(0_+) &= M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \end{aligned}\tag{} Mz(0+)Mxy(0+)=0=Mz0(12e(TRTE/2)/T1eTR/T1)()

根据弛豫方程(式3.122)
{Mx′y′(t)=Mx′y′(0+)e−t/T2Mz′(t)=Mz0(1−e−t/T1)+Mz′(0+)e−t/T1(3.122) \left\{ \begin{array}{l} M_{x'y'}(t) = M_{x'y'}(0_+) e^{-t/T_2} \\ M_{z'}(t) = M_z^0 (1 - e^{-t/T_1}) + M_{z'}(0_+) e^{-t/T_1} \end{array} \right. \tag {3.122} {Mxy(t)=Mxy(0+)et/T2Mz(t)=Mz0(1et/T1)+Mz(0+)et/T1(3.122)

0≤t<TI0 \leq t < T_I0t<TI 时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−t/T1)(0<t<TI)Mx′y′(t)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−t/T2(0<t<TI)() \begin{aligned} M_{z'}(t) &=M_z^0 \left(1 - e^{-t/T_1}\right) &\quad (0 < t < T_I) \\ M_{x'y'}(t) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) e^{-t/T_2} &\quad (0 < t < T_I) \end{aligned} \tag{} Mz(t)Mxy(t)=Mz0(1et/T1)=Mz0(12e(TRTE/2)/T1eTR/T1)et/T2(0<t<TI)(0<t<TI)()

t=TI−t = T_{I-}t=TI时刻,即180°脉冲施加前,横向和纵向磁化强度分别为:
Mz′(TI−)=Mz0(1−e−TI/T1)Mx′y′(TI−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−TI/T2() \begin{aligned} M_{z'}({T_I}_-) &=M_z^0 \left(1 - e^{-{T_I}/T_1}\right) &\quad \\ M_{x'y'}({T_I}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) e^{-{T_I}/T_2} &\quad \end{aligned} \tag{} Mz(TI)Mxy(TI)=Mz0(1eTI/T1)=Mz0(12e(TRTE/2)/T1eTR/T1)eTI/T2()

理想180°射频脉冲给横向和纵向磁化带来的影响如下:
Mz′→180∘x′−Mz′() M_{z'} \xrightarrow{{180^\circ}_{x'}} -M_{z'} \tag{} Mz180xMz()
Mx′y′→180∘x′−Mx′y′() M_{x'y'} \xrightarrow{{180^\circ}_{x'}} -M_{x'y'} \tag{} Mxy180xMxy()

t=TI+t = T_{I+}t=TI+时刻,即180°脉冲施加后,横向和纵向磁化强度分别为:
Mz′(TI+)=−Mz′(TI−)=−Mz0(1−e−TI/T1)Mx′y′(TI+)=−Mx′y′(TI−)=−Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−TI/T2() \begin{aligned} M_{z'}({T_I}+) &= -{M_{z'}({T_I}_-)} = -M_z^0 \left(1 - e^{-{T_I}/T_1}\right)\\ M_{x'y'}({T_I}+) &= -{M_{x'y'}({T_I}_-)} = -M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) e^{-{T_I}/T_2} \end{aligned} \tag{} Mz(TI+)Mxy(TI+)=Mz(TI)=Mz0(1eTI/T1)=Mxy(TI)=Mz0(12e(TRTE/2)/T1eTR/T1)eTI/T2()

t>TI+t > T_{I+}t>TI+时刻,横向和纵向磁化强度分别为:
Mz′(t)=Mz0(1−e−(t−TI)/T1)−Mz0(1−e−TI/T1)e−(t−TI)/T1=Mz0(1−e−(t−TI)/T1)−Mz0(e−(t−TI)/T1−e−t/T1)=Mz0(1−2e−(t−TI)/T1)+Mz0e−t/T1=Mz0(1−2e−(t−TI)/T1+e−t/T1)Mx′y′(t)=−Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)e−t/T2() \begin{aligned} M_{z'}(t) &= M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(1 - e^{-T_I/T_1} \right) e^{-(t-T_I)/T_1} &\quad\\ &=M_z^0 (1 - e^{-(t-T_I)/T_1}) -M_z^0 \left(e^{-(t-T_I)/T_1} - e^{-t/T_1} \right) \\ &=M_z^0 (1 - 2e^{-(t-T_I)/T_1}) + M_z^0 e^{-t/T_1} \\ &=M_z^0 (1 - 2e^{-(t-T_I)/T_1}+e^{-t/T_1}) \\ \\ M_{x'y'}(t) &= -M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1})e^{-t/T_2} &\quad \end{aligned} \tag{} Mz(t)Mxy(t)=Mz0(1e(tTI)/T1)Mz0(1eTI/T1)e(tTI)/T1=Mz0(1e(tTI)/T1)Mz0(e(tTI)/T1et/T1)=Mz0(12e(tTI)/T1)+Mz0et/T1=Mz0(12e(tTI)/T1+et/T1)=Mz0(12e(TRTE/2)/T1eTR/T1)et/T2()

由上式可得:

  • 第一次循环后系统就已经进入了稳态,即后续的循环的起始时刻t=0−t=0_-t=0时的纵向磁化强度都是:
    Mz′(0−)=Mz0(1−2e−(TR−TE/2)/T1−e−TR/T1)(7.24) \begin{aligned} M_{z'}({0}_-) &=M_z^0 (1 - 2e^{-(T_R-T_E/2)/T_1}-e^{-T_R/T_1}) \\ \end{aligned} \tag{7.24} Mz(0)=Mz0(12e(TRTE/2)/T1eTR/T1)(7.24)

  • 自旋回波信号的幅度可表示为:
    AE=Mz0(1−2e−(TR−TE/2)/T1+e−TR/T1)e−TE/T2(7.25) A_E = M_z^0 \left(1 - 2e^{-(T_R - T_E/2)/T_1} + e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.25} AE=Mz0(12e(TRTE/2)/T1+eTR/T1)eTE/T2(7.25)

在实际成像中,通常满足TE≪TRT_E \ll T_RTETR的条件(即回波时间远小于重复时间)。在此条件下,可对公式(7.25)进行合理简化:当TE≪TRT_E \ll T_RTETR时,有TR−TE/2≈TRT_R - T_E/2 \approx T_RTRTE/2TR,所以:
e−(TR−TE/2)/T1≈e−TR/T1 e^{-(T_R - T_E/2)/T_1} \approx e^{-T_R/T_1} e(TRTE/2)/T1eTR/T1

代入公式(7.25)得:
AE=Mz0(1−2e−TR/T1+e−TR/T1)e−TE/T2=Mz0(1−e−TR/T1)e−TE/T2(7.26) A_E = M_z^0 \left(1 - 2e^{-T_R/T_1} + e^{-T_R/T_1}\right)e^{-T_E/T_2} = M_z^0 \left(1 - e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.26} AE=Mz0(12eTR/T1+eTR/T1)eTE/T2=Mz0(1eTR/T1)eTE/T2(7.26)

选择性对比度强调

公式(7.26)表明,饱和恢复自旋回波序列的图像强度同时具有三种加权特性:

T1T_1T1加权对比度:采用短TET_ETE和适当的TRT_RTR

  • TET_ETE使e−TE/T2→1e^{-T_E/T_2} \rightarrow 1eTE/T21,消除T2T_2T2加权效应
  • 适当TRT_RTR确保(1−e−TR/T1)(1 - e^{-T_R/T_1})(1eTR/T1)产生明显的T1T_1T1对比

T2T_2T2加权对比度:采用长TET_ETE和长TRT_RTR

  • TRT_RTR使(1−e−TR/T1)→1(1 - e^{-T_R/T_1}) \rightarrow 1(1eTR/T1)1,消除T1T_1T1加权效应
  • TET_ETE使e−TE/T2e^{-T_E/T_2}eTE/T2产生明显的T2T_2T2对比

自旋密度加权对比度:采用短TET_ETE和长TRT_RTR

  • TET_ETE消除T2T_2T2加权
  • TRT_RTR消除T1T_1T1加权
  • 剩余信号正比于自旋密度

下表总结了不同序列参数值对图像对比度的影响。
image.png

反转恢复自旋回波序列

image.png

上展示了反转恢复自旋回波序列的时序图。该序列首先施加180°反转脉冲,随后在反转时间TIT_ITI后施加90°激发脉冲,最后通过180°重聚脉冲产生自旋回波。

对于反转恢复自旋回波序列,自旋回波信号的幅度可推导为:
AE=Mz0(1−2e−TI/T1+2e−(TR−TE/2)/T1−e−TR/T1)e−TE/T2(7.27) A_E = M_z^0 \left(1 - 2e^{-T_I/T_1} + 2e^{-(T_R - T_E/2)/T_1} - e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.27} AE=Mz0(12eTI/T1+2e(TRTE/2)/T1eTR/T1)eTE/T2(7.27)

同样基于TE≪TRT_E \ll T_RTETR的实用假设,公式(7.27)可简化为:
AE=Mz0(1−2e−TI/T1+e−TR/T1)e−TE/T2(7.28) A_E = M_z^0 \left(1 - 2e^{-T_I/T_1} + e^{-T_R/T_1}\right)e^{-T_E/T_2} \tag{7.28} AE=Mz0(12eTI/T1+eTR/T1)eTE/T2(7.28)

该表达式表明,反转恢复自旋回波序列能够产生与饱和恢复自旋回波序列相似的T1T_1T1T2T_2T2和自旋密度加权对比度,但具有更强的T1T_1T1对比度调控能力。

结论

基本自旋回波成像序列通过将自旋回波技术引入饱和恢复和反转恢复序列,实现了对图像对比度的精细调控。公式(7.24)至(7.28)完整描述了这些序列的信号特性,为临床MRI扫描参数优化提供了理论基础。

通过合理选择TRT_RTRTET_ETETIT_ITI等参数,可以获得满足特定诊断需求的图像对比度,这是自旋回波序列在临床磁共振成像中得以广泛应用的重要原因。

http://www.dtcms.com/a/610948.html

相关文章:

  • 为什么 AI 推理天然更偏向 ARM(而不是 x86)
  • 广州网站建设哪里有国家信息公示网查营业执照年报
  • 网站 需求文档益阳市网站建设
  • 基于数据的人工智能建模流程及源码示例
  • 站长之家爱站网f006网站建设
  • 中国建设银行网上银行官方网站wordpress 隐私设置
  • 4网站建设软件工程师证书报考条件
  • Navicat把sqlserver数据库转到mysql
  • 在 Linux 系统上安装 JDK 8u202 Linux x64 RPM 包
  • FreeBSD 14.3 轻量级Jail虚拟机:内存资源占用仅13MB的实战指南
  • swift 中 CGImage CIImage UIImage 等 有什么区别
  • 智能网站开发自己建网站的优势
  • 下载flash网站蓄电池回收网站建设
  • 高级边界扫描 --4-- 使用Silicon Nail测试
  • 商品历史价格数据采集:基于DrissionPage的智能爬虫实战
  • mysql锁相关
  • 网站设计网站维护网站备案资料 下载
  • 右半平面零点(RHPZ)的挑战与解决方案:Boost和Flyback变换器稳定性深度解析
  • 张家港做网站排名网站排名标准
  • C# 容器实例生命周期
  • OFD(Open File Description)
  • LLM思维链提示(CoT)
  • 重庆网站建设快速建站中山优秀网站建设
  • 【愚公系列】《MCP协议与AI Agent开发》011-MCP协议标准与规范体系(交互协议与状态码体系)
  • 算法学习--快速幂
  • 快速排序简要介绍
  • 阆中 网站建设html网页制作小刘在线课程
  • 无锡网站建设无锡速联科技做网站首页可以用传媒公司吗
  • 建站公司网站 phpwind自己做的网站怎么给别人访问
  • 百度世界2025|实时互动数字人智能体让搜索更懂人心