try ai
科普
编辑
分享
反馈
  • 相位编码

相位编码

SciencePedia玻尔百科
核心要点
  • 相位编码使用一个短暂、受控的磁场梯度,为质子自旋赋予独特的相位偏移,从而有效地标记其空间位置。
  • 在磁共振成像中,通过以不同的梯度强度重复相位编码过程,逐行填充一个称为k空间的数据网格来构建图像。
  • 回波平面成像 (EPI) 等先进方法通过单次激发后捕获所有k空间数据,加速了这一过程,从而实现了功能性磁共振成像 (fMRI) 等动态应用。
  • 将信息编码到相位中的原理超越了空间定位,可用于测量运动,例如相位对比磁共振血管成像中的血流,并且在地球物理学和神经科学等领域具有概念上的相似之处。

引言

磁共振成像 (MRI) 为我们提供了一扇观察人体的无与伦比的窗口,但它究竟是如何将旋转质子的简单信号转化为精细的解剖学图像的呢?其根本挑战在于空间定位:如果磁场中所有质子都以几乎相同的频率共振,那么区分它们就像试图绘制一座每个居民地址都相同的城市地图。相位编码是解决这个问题的巧妙方案之一,它优雅地为质子标记上其位置的记忆。这是一个基础概念,不仅使成像成为可能,还开启了大量先进的诊断和研究能力。

本文将相位编码分解为其核心组成部分,以揭开其神秘面纱。在第一章“原理与机制”中,我们将探讨磁场梯度的工作原理、相位如何被用作空间标签,以及这个过程如何让我们在抽象的k空间域中系统地构建图像。随后,“应用与跨学科联系”一章将展示这一原理的多功能性,从利用功能性磁共振成像 (fMRI) 实时捕捉大脑活动,到测量血流、校正图像畸变,甚至在地球物理学和神经科学等截然不同的领域中发现其惊人的相似之处。

原理与机制

自旋的交响乐

想象你置身于一个巨大而黑暗的音乐厅。舞台上是一个规模超乎想象的管弦乐队,其中每个音乐家都是患者体内一个微小的旋转质子。在MRI扫描仪的均匀磁场中,所有这些质子都以几乎完全相同的频率进动——它们像微小的陀螺一样摆动——这个频率被称为​​拉莫尔频率​​。就好像整个管弦乐队正在演奏一个完美同步、永无止境的单音。声音虽然强大,但只是一个音符。我们如何能从这种单调的嗡鸣中创造出一幅精细的画面呢?要制作一幅地图,我们需要变化。我们需要一种方法来问每个质子:“你在哪里?”并得到一个独特的答案。

这就是磁共振成像的根本挑战。如果所有质子都唱着同一首歌,我们就无法分辨信号是来自头部还是脚部。为了创建图像,我们必须打破这种完美的均匀性。我们需要给管弦乐队一个更复杂的乐谱,一个每个音乐家的部分都取决于其在舞台上位置的乐谱。实现这一点的关键是一个非常优雅的工具:​​磁场梯度​​。

磁场中的扭转:梯度

磁场梯度是一个简单但深刻的想法。我们不再保持磁场完全均匀,而是有意地让它在身体内以可控的方式变化。假设我们想要编码垂直位置,我们称之为 yyy 轴。我们施加一个梯度 GyG_yGy​,使得磁场在顶部(正 yyy)稍强,在底部(负 yyy)稍弱。在任何位置 yyy 的总磁场不再仅仅是主磁场 B0B_0B0​,而是变得与位置相关:

B(y)=B0+GyyB(y) = B_0 + G_y yB(y)=B0​+Gy​y

由于进动频率与磁场成正比(ω=γB\omega = \gamma Bω=γB,其中 γ\gammaγ 是旋磁比),质子们不再齐声歌唱。它们的进动频率现在取决于它们的 yyy 轴位置:

ω(y)=γ(B0+Gyy)=ω0+γGyy\omega(y) = \gamma (B_0 + G_y y) = \omega_0 + \gamma G_y yω(y)=γ(B0​+Gy​y)=ω0​+γGy​y

突然间,我们的管弦乐队有了空间结构。顶部的质子进动得快一些,底部的质子进动得慢一些。我们给每一排音乐家都分配了自己独特的节拍。这个原理是MRI中所有空间编码的基石。然而,我们如何使用它,决定了具体的编码方法。

相位的记忆

现在来看​​相位编码​​这个极其聪明的技巧。我们不是在“聆听”信号时一直保持梯度开启,而是只在一个非常短暂、精确控制的持续时间 τ\tauτ 内打开它。

在这短暂的时间窗口内,不同 yyy 位置的自旋以其独特的、与位置相关的速度进动。一个在较高 yyy 值的自旋,进动速度更快,其摆动会比中心的自旋稍微超前一些。一个在较低 yyy 值的自旋则会稍微落后一些。在脉冲持续时间 τ\tauτ 结束时,我们关闭梯度。瞬间,所有自旋都回到以完全相同的基准频率 ω0\omega_0ω0​ 进动的状态。

但有些东西改变了。位于位置 yyy 的自旋相对于中心位置的自旋累积了一个额外的相位——一个角度。这个相位偏移 Δϕ(y)\Delta\phi(y)Δϕ(y) 是它刚刚经历的梯度脉冲的“记忆”。我们可以很容易地计算它。额外的频率是 γGyy\gamma G_y yγGy​y。在时间 τ\tauτ 内,累积的总额外相位就是额外频率乘以时间:

Δϕ(y)=(γGyy)τ\Delta\phi(y) = (\gamma G_y y) \tauΔϕ(y)=(γGy​y)τ

这就是问题的核心。尽管所有自旋现在都再次以相同的速率进动,但它们已不再同步。每个 yyy 位置现在都由一个独特的相位偏移来标记。这就像一群完全相同的时钟,都以相同的速率滴答作响。如果我们短暂地让一些时钟走得快些,另一些走得慢些,仅仅一分钟,然后让它们全部恢复正常速率,它们将永远显示不同的时间。它们显示的时间差异,就是它们滴答速度暂时改变的永久记录——一种相位记忆。

这与它的同类方法——​​频率编码​​和​​层面选择​​——有着根本的不同。在频率编码中,我们在信号采集期间保持梯度开启,因此位置由我们听到的不同频率来编码。在层面选择中,我们在初始射频脉冲期间施加梯度,以便从一开始就只激发身体的一个薄层。相位编码则是在主要测量开始之前用相位标记位置的艺术。

构成图像:k空间的语言

我们现在有了一种方法,可以让来自每个 yyy 位置的信号携带一个独特的标签:它的相位。我们如何读取这个标签并用它来构建一幅图像呢?答案在于科学和工程领域最强大的工具之一:​​傅里叶变换​​。

我们接收到的总信号是来自被激发层面中所有质子信号的总和。来自位置 yyy 的信号被其特殊的相位因子 exp⁡(−iΔϕ(y))\exp(-i\Delta\phi(y))exp(−iΔϕ(y)) 调制。当我们写下这个式子时,总信号的方程看起来正像一个傅里叶变换。傅里叶世界中的数学变量是​​空间频率​​,而所有空间频率构成的景观被亲切地称为​​k空间​​。

我们施加梯度脉冲的物理行为可以直接转化为k空间的数学语言。傅里叶变换的核看起来像 exp⁡(−i2πkyy)\exp(-i 2\pi k_y y)exp(−i2πky​y),其中 kyk_yky​ 是 yyy 方向的空间频率。我们可以简单地将其与我们的物理相位因子等同起来:

2πkyy=Δϕ(y)=γGyyτ2\pi k_y y = \Delta\phi(y) = \gamma G_y y \tau2πky​y=Δϕ(y)=γGy​yτ

位置变量 yyy 被消掉,留下了一座连接我们操作与最终图像数据的美丽而直接的桥梁:

ky=γGyτ2πk_y = \frac{\gamma G_y \tau}{2\pi}ky​=2πγGy​τ​

我们在k空间中采样的位置 (kyk_yky​) 与我们施加的梯度脉冲的强度和持续时间成正比。这个组合项,即梯度幅值与其持续时间的乘积,被称为梯度​​面积​​或​​矩​​。它是我们用来在k空间中导航的基本旋钮。

逐行构建图像

一个单一的梯度脉冲让我们能够测量单个 kyk_yky​ 值处的信号。这为我们提供了关于物体的一小片信息——就像知道一首复杂乐曲中某个单一低音频率的量。要重建图像的完整“乐曲”,我们需要对整个k空间进行采样。

我们通过一遍又一遍地重复整个过程——激发、相位编码、读出——来做到这一点。对于一个在相位编码方向上具有256像素分辨率的典型图像,我们可能需要重复该过程256次。

在每次重复中(每隔一个“重复时间”或 TRT_RTR​ 发生一次),我们施加一个幅度略有不同的相位编码梯度。在标准的“自旋绕”采集中,我们从一个大的负梯度开始,然后在下一次重复时将其步进到一个稍微不那么负的值,依此类推,线性地穿过零点,直到一个大的正梯度。

对于这些相位编码步骤中的每一步(它设定了一个特定的 kyk_yky​ 值),我们接着打开频率编码梯度(沿 xxx 轴),并读出信号。这次单一的读出过程会扫过所有的 kxk_xkx​ 值,填充我们k空间数据网格中的一条水平线。经过256次重复后,我们就在k空间中填充了一个完整的256x256网格。最终可识别的图像,仅仅是这个抽象的k空间数据的二维逆傅里叶变换。第 nnn 条线所需的精确梯度幅度 GnG_nGn​ 可以被精确计算,以确保我们正确地、对称地围绕原点绘制出k空间。

游戏规则:视野与伪影

这种对k空间的逐步采样必须小心进行。我们采样这个抽象空间的方式对我们看到的图像有直接的、物理上的影响。这种关系由傅里叶变换的美妙互易性所支配。

其中一个最关键的参数是k空间各行之间的步长 Δky\Delta k_yΔky​。这个步长与我们最终图像的​​视野​​ (FOV) 成反比:FOVy=1/ΔkyFOV_y = 1/\Delta k_yFOVy​=1/Δky​。如果我们的k空间步长取得太大(即一个大的 Δky\Delta k_yΔky​),我们得到的FOV就会太小。如果我们正在成像的身体部位大于这个FOV,就会出现一种奇异的伪影。位于FOV之外的物体部分会被“折叠”或“卷绕”回图像中,叠加在我们想要看的解剖结构上。这被称为​​混叠​​或​​卷褶伪影​​。为了避免这种情况,我们必须仔细选择我们的梯度增量,以确保FOV足够大,能够包含整个物体。这就为连续步骤之间梯度面积的变化设定了一个最大限制。

这种逐行采集也使得相位编码方向最容易受到运动的影响。整个过程可能需要几分钟。如果病人呼吸或心脏跳动,他们的身体在每个相位编码步骤中都会处于稍微不同的位置。这会引入一个与位置相关的相位误差,这个误差在k空间的每一行之间都不同。对于周期性运动,这种误差模式也是周期性的。当我们进行傅里叶变换得到图像时,k空间中的这种周期性误差会表现为物体的“鬼影”复制品,沿着相位编码方向在图像上涂抹开来。正是相位编码缓慢而审慎的特性,使其对这些动态变化如此敏感。

因此,相位编码是物理学和数学之间的一场精妙舞蹈。通过用一个精心控制的磁场短暂地推动旋转的质子,我们用它们位置的记忆来标记它们。通过以系统变化的推动重复这个过程,我们耐心地聆听自旋的交响乐,并逐个频率地构建出身体的完整傅里叶表示。从这些抽象的数据中,一幅我们内在世界的精细而美丽的图像便浮现出来。

应用与跨学科联系

在我们之前的讨论中,我们深入探究了相位编码的核心,发现它是一种极其巧妙的方法,能将空间信息烙印在旋转质子信号的相位上。这就像教一个合唱团的歌手根据他们在舞台上的站位来改变音高。通过聆听合唱的声音,一个聪明的指挥家就能重构出每个歌手的位置。但这个原理,这种将信息“写入”波的相位的能力,远不止是一个单一的技巧。它是一把万能钥匙,开启了数量惊人的大门,引领我们从自己身体的内部运作到地球的深处,甚至进入思想本身的机制。现在,让我们踏上一段旅程,探索其中一些卓越的应用。

内部观察的艺术:构建图像

相位编码最直接,或许也是最著名的应用,是在磁共振成像 (MRI) 中。我们如何将一个三维物体,比如人脑,变成一幅精细的二维层面图像?

想象一个艺术家的画布,我们称之为“kkk-空间”。这不是一个普通的由位置构成的画布;它是一个由空间频率构成的画布,一本最终图像的“配方书”。要创建我们的图像,我们不是逐个像素地绘画。相反,我们必须填满这整本配方书。相位编码就是我们用来选择要书写配方书哪一行的工具。

对于每一行,一个短暂的磁场梯度被施加在,比如说,垂直 (yyy) 方向上。对于每个连续的行,这个梯度的强度都以离散的步长变化。这就是相位编码步骤。对于弱梯度,我们准备在kkk-空间画布中心附近填充一行;对于强梯度,我们准备在顶部或底部附近填充一行。在画布上的这个垂直“位置”设定好之后,另一个梯度,即频率编码梯度,会沿着水平 (xxx) 方向开启。当这个梯度开启时,我们通过聆听演变的信号,在我们的kkk-空间画布上“绘制”水平线。通过为许多不同的相位编码步骤重复此过程,我们有条不紊地逐行填满整个kkk-空间配方书。一旦画布被填满,一个简单的数学运算——傅里叶变换——就能立即将这份配方转换成我们所熟悉的精美、细致的解剖图像。我们选择的相位编码步数直接决定了图像在该方向上的分辨率;步数越多,图像就越精细、越详细。

对速度的需求:捕捉运动中的思想

逐行扫描的方法虽然细致且能产生精美的图像,但速度很慢。如果我们想对快速变化的事物进行成像,比如思考时的大脑,或者跳动时的心脏,该怎么办呢?等待采集数百条独立的线,就像试图用19世纪的相机拍摄蜂鸟;等到照片拍好,鸟儿早已飞远了。

这正是相位编码在一种名为回波平面成像 (EPI) 的技术中大放异彩的地方。与庄重、逐行推进的方式不同,EPI是一次性在整个kkk-空间画布上疯狂、闪电般地涂鸦。单次激发后,一个强的频率编码梯度被迅速地来回切换,使“画笔”在kkk-空间左右扫动。在这些扫动之间的微小瞬间,一个非常短暂、尖锐的相位编码梯度“脉冲”被施加。这个脉冲的强度恰好足以将画笔踢到画布的下一行。之字形扫动,脉冲,反向扫动,脉冲,再之字形扫动... 在几十毫秒内,整本配方书被填满,一幅完整的图像被捕获。

这种令人难以置信的速度是功能性磁共振成像 (fMRI) 的基础,这项技术使我们能够观察与神经活动相关的大脑血流变化。通过使用EPI每隔一两秒捕获整个大脑的图像,我们可以制作一部大脑工作的电影。当被推向极限时,相位编码为我们提供了一扇窥探认知本身的窗户。

当“颜料”溢出时:伪影与巧妙的修正

当然,任何强大的技术都不是没有代价的。我们用来编码的相位敏感性本身也是麻烦的来源。MRI扫描仪内部的磁场永远不会完美均匀,尤其是在不同组织之间的边界处,比如充满空气的鼻窦和脑组织。这些场的不完美性会导致一种与所施加梯度无关的“离共振”相位偏移。

在缓慢的常规扫描中,这些影响通常很小。但在快速的EPI扫描中,长采集窗口给了这种不必要的相位足够的时间来累积。扫描仪对离共振一无所知,将这种额外的相位误解为空间编码的一部分。结果是几何畸变——图像的扭曲、拉伸和挤压,在相位编码方向上最为严重。图像就像哈哈镜中的反射。

我们该如何修正这个问题?解决方案是一招漂亮的科学“柔道”。我们采集第二幅EPI图像,但这一次,我们施加极性相反的相位编码脉冲。这会产生反转畸变方向的效果。如果一个脑区在第一幅图像中向右拉伸,那么在第二幅图像中它将被向左挤压。现在,我们有了两幅不同但相关的哈哈镜图像。一个聪明的计算机算法可以同时观察两者,计算出“镜子”(离共振场)的确切形状,并计算出大脑在被扭曲之前必定是什么样子。这种方法,通常被称为“反转相位编码”法,是现代神经影像学的基石,使我们能够获得几何上忠实的大脑功能和结构图。同样这种操纵相位以校正误差的精神,也体现在更先进的技术中,用于金属植入物附近的成像 或解开同时采集的多个层面信号。

超越位置:编码运动与化学信息

到目前为止,我们已经用相位来编码物体的位置。但这个原理更具普适性。我们可以编码任何随时间影响自旋感受磁场的物理量。例如,如果我们不仅想对血管的解剖结构成像,还想对其中流动的血液成像,该怎么办?

这是通过相位对比磁共振血管成像 (PC-MRA) 实现的。其技巧是使用一个“双极”梯度脉冲——一个正叶紧跟着一个相同的负叶。对于静止的自旋,第一叶期间获得的相位被第二叶期间失去的相位完美抵消。净相位偏移为零。但对于一个正在移动的自旋,比如流经动脉的血细胞,它在第二叶时的位置将与第一叶时不同。它经历了不同的磁场历史,相位抵消不再完美。它会留下一个与它的速度成正比的残余相位!通过从这个依赖于速度的相位创建图像,我们可以实时可视化血流,诊断堵塞和动脉瘤,而无需接触患者。

这还不止。在磁共振波谱成像 (MRSI) 中,我们可以使用相位编码来创建空间图谱,然后,对于图谱上的每个位置,我们使用另一种技术,根据不同化学物质独特的共振频率来分离它们的信号。这使我们能够创建器官的化学图谱,观察代谢物或神经递质的分布。这是研究肿瘤和神经系统疾病中新陈代谢的有力工具。当然,就像任何采样系统一样,如果存在来自我们预期视野之外的化学信号,它的信号可能会“卷绕”并作为鬼影出现在我们的图像中,这正是支配所有傅里叶成像的相同采样原理的直接后果。

地球与心智中的回响:相位的普适语言

一个深刻物理原理的真正美在于其普适性。利用相位编码信息的想法并不仅限于MRI扫描仪的纯净环境。

考虑计算地球物理学领域。为了绘制地下深处的岩层,地球物理学家会引爆受控的爆炸(震源)并记录反射回来的地震波。为了提高效率,他们希望同时引爆多个震源。但这会在接收器处产生一堆混乱的重叠信号。如何将它们解开呢?答案出人意料,是随机相位编码。通过在每个震源被引爆前对其施加一个独特的、随机的时间延迟(相位偏移),所产生的“非相干”信号叠加可以通过压缩感知的数学工具成功地分离出来。这假设了潜在的地质结构相对简单,或称“稀疏”。这正是在一个完全不同的背景下,用来窥探地球的完全相同的原理——用相位标记信息以便后续分离。

也许所有联系中最深刻的,不是我们工程设计的,而是我们发现的。在我们自己的大脑内部,在海马体——记忆和导航的中心——进化已经发展出了它自己的相位编码形式。被称为“位置细胞”的特化神经元会在动物处于特定位置,即其“位置野”时放电。但大脑所知晓的,不仅仅是动物在那个位置野里。海马体产生一种稳定、有节奏的电波,称为θ节律,这是一个内部时钟。当一只老鼠跑过一个位置野时,它的位置细胞在θ波的越来越早的相位上放电。这种现象被称为“相位进动”,它编码了动物在位置野内的精确位置。这是一种相位编码的生物学实现,一种关于空间神经代码,与我们为成像机器开发的技术惊人地相似。

从构建图像,到观察其运动,到修正其缺陷,到看到流动和化学成分,最后,到在地球地壳和我们自己心智的机制中发现其回响——相位编码的旅程是科学统一性与优雅的有力例证。一个单一、简单的想法,在好奇心和创造力的驱使下,揭示了一个充满隐藏联系和可能性的世界。