并行磁共振成像——SENSE 敏感度编码成像
SENSitivity Enconding(SENSE) 敏感度编码重建是一种基于图像域的并行成像重技术,是在1999年由 Klaas P. Pruessmann 提出的。
SENSE 重建的核心是从欠采样具有伪影的图像中,多个线圈采集得到的有伪影图像中恢复出没有伪影完整的图像。
MRI 成像过程中每一个线圈对于靠近线圈部分的成像物体敏感度最高,随着距离增加敏感度下降,即线圈的成像具有局部敏感性(local sensitivity)。线圈对于不同空间位置磁共振信号的敏感程度即为线圈敏感度。
在并行成像中,根据不同线圈对图像的敏感度对磁共振图像进行加权即可以得到一张完整的图像,但是在上图所示的情况下 k 空间中并没有进行欠采样,也就是没有实现加速成像。而在 k 空间欠采样会导致图像出现重叠伪影(Phase wrapping),这是因为
K 空间为频域空间,在频域中欠采样,即对频域信号进行离散化,在时域中对应的有限长序列为:原序列以采样点数为周期进行周期延拓后的主值序列
k空间欠采样后重建得到具有伪影的图像可以看成是原图像的超空间混叠,利用先验的线圈敏感度和FOV知识便可以分开这种空间混叠便可以得到没有伪影的图像。
SENSE就是利用敏感度编码,从而实现可以基于伪影图像的重建。
SENSE 主要分为三个步骤
一、线圈敏感度估计
线圈敏感度估计是SENSE算法中重要的一部分,在实际采集之前先计算各线圈的敏感度。首先在获取每个线圈的低分辨率全视场图像,对于图像进行归一化,然后平滑滤波,阈值处理以及点估计从而生成平滑的敏感度图像,即每个线圈接收区域内不同源点的信号的相对权重。
二、获取 k 空间欠采样数据
在 k 空间进行欠采样,从而实现加速成像,每个线圈获得了具有混叠的欠FOV图像。欠采样倍数受到采集线圈数量的限制,即欠采样倍数要小于采集线圈数量。
三、混叠图像重建
目前,已知了线圈的敏感度和混叠得到的信号,混叠信号可以看作全FOV内信号的信号叠加。以2倍 k 空间欠采样为例来举例。
线圈 1 和 2 用于检测该信号的灵敏度将分别表示为 S 1 A S_{1A} S1A 和 S 2 A S_{2A} S2A。类似地,任何其他点 B 的线圈灵敏度也是已知的,并将表示为 S1B 和 S2B。
用 P 1 P_1 P1 和 P 2 P_2 P2 表示来自线圈 1 和 2 的 重建后的结果,为2倍混叠的1/2 FOV的混叠图像。
P
1
=
A
S
1
A
+
B
S
1
B
P_1 =A S_{1A}+BS_{1B}
P1=AS1A+BS1B
P
2
=
A
S
2
A
+
B
S
2
B
P_2 = AS_{2A}+BS_{2B}
P2=AS2A+BS2B
其中,S 都已知,P都已知,所以只需要解线性方程就可以得到 A 和 B 两点的图像信号了。
参考资料
- SENSE/ASSET? - Questions and Answers in MRI (mriquestions.com)
- 快速磁共振成像 张明辉 刘且根 徐晓玲