4. 图像域与k空间域伪影机理

MRI伪影的根源几乎都可以追溯到两个域之间的映射关系：k空间是采集域，图像域是重建结果，二者通过傅里叶变换相互绑定。理解伪影，关键不是记住“伪影长什么样”，而是判断它发生在哪个域、由什么采样或相位错误引起，以及工程上能在哪一环切断它。本章按“k空间采样错误→图像域表现→抑制策略”的链条，逐一分析四类高频伪影。

一、截断伪影：k空间高频缺失导致的振铃

理想图像需要k空间从负无穷到正无穷的完整采样，实际只能在有限范围 [-k_max, +k_max] 内采集。这相当于对真实k空间乘上一个矩形窗，其傅里叶变换是sinc函数，因此图像中任何强边缘（如颅骨与脑脊液、脊髓与椎管）都会被sinc旁瓣“振铃”污染，表现为平行于边缘的明暗条纹，即截断伪影（Gibbs伪影）。

工程抑制有三条主线。第一是提高采样范围，即增加 k_max，直接提高分辨率并削弱振铃，但受梯度性能和SNR限制。第二是k空间滤波，在重建前对k空间边缘施加平滑窗（如Hamming、Hanning、Tukey窗），代价是分辨率轻微下降。第三是部分傅里叶与零填充的配合：零填充只改变像素矩阵大小，不增加真实高频信息，不能消除截断伪影，只能让振铃在更细的网格上呈现。实际序列中常把Tukey窗的过渡带设在k空间外侧10%~20%区域，在1.5T头颅T1序列上可明显减轻脑室-脑实质界面的条纹，同时保持边缘锐度可接受。

二、混叠：FOV与采样间隔不匹配

混叠（wrap-around）来自采样定理：当被成像物体在相位编码方向超出设定FOV时，超出部分会折叠到对侧。其数学本质是k空间相位编码步长 Δk_y 与FOV满足 FOV = 1/Δk_y，若真实物体宽度大于FOV，则高频相位信息被欠采样，图像域出现周期性折叠。

频率编码方向的混叠通常由模拟滤波器带宽限制，表现为图像两端卷折；相位编码方向混叠则更常见，尤其在矢状位脊柱、腹部冠状位等解剖超出FOV的场景。抑制策略包括：增大FOV（降低分辨率或增加相位编码步数）、在相位编码方向施加过采样（如相位FOV 100%→125%，多采集25%数据后裁掉）、使用饱和带（SAT）压制FOV外组织信号、以及并行成像中的参考行校准与展开算法（SENSE/GRAPPA）在重建时解折叠。工程上，腹部冠状位T2序列常设相位编码为前后方向并加前腹壁饱和带，可把呼吸运动导致的卷折与运动伪影一并压低。

三、运动伪影：k空间相位一致性被破坏

MRI采集是逐行（或逐段、逐回波）填充k空间的，任何在采集期间发生的运动都会让不同k空间行对应不同时刻的解剖状态，破坏相位一致性。图像域表现为沿相位编码方向的鬼影（周期性复制）和模糊。呼吸、心跳、肠蠕动、吞咽、眼球运动都是常见来源。

抑制策略按“减少运动发生—补偿运动影响—重建端校正”三层展开。减少发生：屏气、呼吸门控、心电门控、镇静或快速序列缩短采集时间。补偿影响：导航回波（navigator）在每次采集前测量膈肌位置，只接受落在预设窗口内的数据；径向（radial）或螺旋（spiral）轨迹因过采样中心k空间，对运动不敏感，配合自校准重建可显著降低鬼影。重建端校正：基于k空间中心相位差估计运动位移，对各行做相位校正；或使用压缩感知/低秩约束在迭代中分离运动分量。工程上，腹部T2呼吸触发序列常设接受窗口为呼吸周期的30%~50%，配合脂肪抑制和径向采集，可在自由呼吸下获得可诊断图像。

四、化学位移伪影：频率差在频率编码方向的空间错位

脂肪与水质子的化学位移约为3.5 ppm。在1.5T下，进动频率差约 Δf ≈ 3.5 × 10⁻⁶ × 63.87 MHz ≈ 224 Hz；在3.0T下约448 Hz。频率编码方向的空间位置由频率映射决定，因此脂肪信号会相对水信号沿频率编码方向偏移，偏移像素数 Δx = Δf / (BW_per_pixel)。例如1.5T、接收带宽 BW = 32 kHz、矩阵256时，每像素带宽约125 Hz，偏移约1.8像素；若带宽降到16 kHz，偏移约3.6像素，在脂肪-水界面（如肾周脂肪、椎体骨髓）出现明显黑白勾边。

抑制策略包括：提高接收带宽（缩短读出时间，降低偏移像素数，但SNR下降）；施加脂肪抑制（CHESS、SPIR、SPAIR）直接消除脂肪信号；使用Dixon水脂分离（两点或三点法）在重建中把水与脂肪分开，既消除错位又获得脂肪分数图；在频率编码方向选择脂肪-水界面较少的方位，或把频率编码方向设在左右方向以避开主要界面。工程上，脊柱矢状位T2常把频率编码设为前后方向并加SPAIR脂肪抑制，可同时压制化学位移和脂肪高信号。

五、工程抑制策略的协同与取舍

实际序列中，四类伪影往往同时存在，抑制手段之间也存在耦合。例如提高带宽可减轻化学位移，但会降低SNR并可能迫使增加平均次数，延长采集时间，从而加重运动伪影；脂肪抑制可消除化学位移，但SPAIR等脉冲会消耗SAR并可能影响TR。因此工程上需要按“伪影严重度—SNR—采集时间—SAR”四个维度做权衡。

一个可参考的腹部冠状位T2加权序列参数示例（仅作教学示例，非实测结果）：FOV 380 mm、矩阵320×256、相位编码前后方向、相位FOV 100%、接收带宽31.25 kHz、TR 2000 ms、TE 90 ms、回波链长度16、呼吸触发接受窗口40%、SPAIR脂肪抑制、径向采集配合自校准重建。该配置下，混叠由前后饱和带和相位过采样控制，运动由呼吸触发和径向轨迹控制，化学位移由SPAIR和高带宽控制，截断伪影由Tukey窗和足够k_max控制。实际使用时应根据设备梯度性能、SAR余量和患者配合度调整，并通过观察图像边缘、对侧折叠区、相位编码方向鬼影和脂肪-水界面来验证抑制效果。

总结来看，伪影抑制不是单一开关，而是从k空间采样设计、序列时序、接收链带宽到重建算法的系统工程。掌握每个伪影的域归属和数学根源，才能在参数表上做出有依据的取舍，而不是盲目试错。