SPM12实战：手把手教你搞定fMRI数据预处理（从时间矫正到空间平滑）
SPM12实战：零基础入门fMRI数据预处理全流程解析
第一次接触功能磁共振成像（fMRI）数据分析时，面对SPM12复杂的界面和晦涩的术语，很多新手都会感到无从下手。这篇文章将带你从零开始，用最直观的方式掌握fMRI数据预处理的完整流程。不同于简单的步骤罗列，我们将重点关注每个环节背后的原理、常见问题解决方案以及实际案例演示，让你不仅能按部就班完成操作，更能理解为什么要这样做。
1. 环境准备与数据导入
在开始预处理之前，确保你的电脑已经安装了MATLAB（建议R2018b或更新版本）和SPM12工具箱。SPM12的安装非常简单：
- 从官网下载最新版SPM12压缩包
- 解压到MATLAB工具箱目录（如
C:\toolbox\） - 在MATLAB命令窗口输入：
addpath('C:\toolbox\spm12') savepath spm fmri
注意：如果遇到路径错误，检查MATLAB是否以管理员权限运行
准备测试数据时，建议使用公开数据集如OpenNeuro上的ds000117，它包含了完整的fMRI和结构像数据。下载后按照以下结构组织文件：
project_folder/ ├── func/ # 功能像数据 ├── anat/ # 结构像数据 └── deriv/ # 处理后数据重要提示：原始DICOM数据需先用
spm_dicom_convert转换为NIfTI格式，转换时保持默认参数即可，不要修改体素尺寸。
2. 时间层校正：解决扫描时序问题
fMRI采集时并非同时获取所有切片，而是逐层扫描。时间层校正（Slice Timing）就是补偿这种时间差异的关键步骤。在SPM12界面选择Slice Timing模块后，需要特别注意以下参数：
| 参数 | 说明 | 典型值 |
|---|---|---|
| Number of Slices | 图像切片数量 | 36 |
| TR | 重复时间 | 2s |
| TA | 采集时间 | TR-(TR/NSlices) |
| Slice Order | 扫描顺序 | 隔层升序(1:2:36) |
| Reference Slice | 参考层 | 中间层(18) |
计算TA的MATLAB示例：
TR = 2; % 秒 nSlices = 36; TA = TR - (TR/nSlices)常见问题排查：
- 图像显示异常：检查
Data中是否选择了正确的4D文件（通常以.nii结尾） - 结果文件缺失：确认输出目录有写入权限，文件名前缀不要包含特殊字符
- 时间参数错误：使用
spm_dicom_header查看原始DICOM中的TR和切片数
3. 头动校正：消除被试移动影响
即使被试尽量保持不动，微小的头部移动也会严重影响fMRI数据质量。Realign模块通过估计每个时间点的位移来校正这种影响。操作要点：
- 选择
Realign (Est & Res)模块 - 在
Data选项卡添加时间校正后的文件（通常以a开头） - 保持默认参数运行
校正后会生成rp_*.txt文件，包含6列运动参数（3平移+3旋转）。用以下代码可视化头动情况：
motion = load('rp_func_001.txt'); plot(motion(:,1:3)); % 平移 legend('X','Y','Z'); title('头部平移运动');临界值参考：平移>2mm或旋转>2°的数据建议剔除
对于亚洲被试，可能需要调整以下参数：
Quality调至0.95以提高配准精度Separation设为4mm以适应较小的脑结构
4. 空间标准化：统一到标准脑空间
为了使不同被试的数据能够比较，需要将图像配准到标准模板。SPM12提供了针对不同人群的模板选项：
4.1 功能像与结构像配准
- 选择
Coregister模块 Reference Image选择高分辨率的T1结构像Source Image选择头动校正后的平均功能像（mean*文件）
4.2 结构像标准化
- 选择
Normalise模块 Image to Align选择T1结构像- 关键参数设置：
Template Image: 选择TPM.niiAffine Regularisation: 亚洲被试选East AsianVoxel sizes: 设为[2 2 2]
注意：标准化后的文件会以w开头，如wT1.nii
5. 空间平滑：提高信噪比
空间平滑通过高斯核卷积减少噪声，但过度平滑会损失细节。推荐设置：
- FWHM值：通常为体素大小的2-3倍（如3×3×3mm体素用6mm）
- 核类型：保持默认的
Gaussian - 图像选择：标准化后的功能像（
swr*文件）
平滑效果检查方法：
spm_check_registration('smoothed.nii'); spm_orthviews('context_menu','show_smooth',true);6. 预处理质量检查
完整的预处理流程结束后，必须进行质量评估：
时间层校正检查：
- 用
spm_check_registration比较校正前后图像 - 确保没有明显的层间错位
- 用
头动校正验证：
- 检查
rp_*.txt中最大位移值 - 用
spm_plot_motion可视化运动曲线
- 检查
标准化评估：
- 将
wT1.nii与标准模板叠加显示 - 确保关键脑区（如海马、杏仁核）对齐良好
- 将
常见问题解决方案：
- 配准失败：尝试改用
DARTEL方法 - 图像变形：检查是否选择了正确的模板（亚洲/欧洲）
- 信噪比低：考虑增加平滑核大小或重新检查原始数据
预处理后的数据现在可以进行统计分析。建议保存完整的处理日志，包括所有参数设置和中间文件路径，这对后续分析和结果复现至关重要。