VBM全流程详解：从CAT12分割到双样本t检验的SPM实操指南

简介：这份资料面向神经影像领域的科研人员与研究生，提供VBM（基于体素形态学分析）软件及其配套学习教程。VBM作为SPM的插件，用于对MRI数据进行标准化、分割和灰质密度校正，以探究大脑结构的细微变化。资源共21个文件，以MATLAB的.m脚本为主，包含VBM5.1-Manual.pdf操作手册、txt说明文档等，压缩包仅226KB，轻量便携，便于快速部署。已有1771人学习下载。内容涵盖VBM2与VBM5.1两个版本的工具包及详细教程，从MNI空间配准、灰质/白质/脑脊液分割到统计参数图解读均有涉及，可帮助初学者理解基于体素的形态学分析流程，并在实战中掌握t值/z值结果的解释方法。适合具备基础MRI知识、希望系统入门VBM分析的医学影像与认知神经科学方向读者。

1. VBM到底是什么：从“肉眼量脑区”到“逐体素找差异”

拿到一批T1像，几十个被试分两组，想找脑结构差异。传统做法是手动勾画ROI，费时不说，一次只能验证一两个预设脑区，换个脑区就得重画一遍。VBM（Voxel-Based Morphometry，基于体素的形态学分析）把这些工作标准化了：把所有人的脑配准到同一个MNI空间，分割出灰质、白质、脑脊液，再逐体素做统计，全脑几十万个点一次比完。耗时从几周压到一天，结果还能用MNI坐标直接写进论文。这套方案特别适合神经影像、老年痴呆、精分等方向的研究者，也是审稿人最熟悉的脑结构分析路径之一。下面从软件下载开始，把SPM+CAT12这条链路完整跑一遍。

2. VBM软件选型：SPM+CAT12与FSL-VBM怎么选，以及下载安装的完整步骤

VBM不是一个单独的安装包，而是“模板配准 + 组织分割 + 统计检验”的组合流程。能干活的主流方案有两套：SPM12配CAT12工具箱，以及FSL的fslvbm命令行流程。选错方案，后面每一步都会觉得别扭：一个卡在图形界面和批处理之间，一个卡在脚本调试和结果可视化上。

2.1 两套主流方案的边界在哪里

SPM+CAT12跑在MATLAB环境里，图形界面和批处理脚本都有，自带数据质量报告（QA），配准用的是DARTEL高维非线性算法，在体素形态学场景下非常成熟。FSL-VBM是纯命令行的bash流程，适合Linux服务器、喜欢全脚本管理的团队，但它的分割和配准分散在多个脚本里，中间产物多，出了问题调试时像个黑匣子。

我的建议很直接：除非你的课题组已经有一套成熟FSL流程在跑，否则新项目一律上SPM+CAT12。原因很现实：CAT12的QA报告能直接回答“这批数据能不能用”，这一条就能帮你挡掉大量返工。

2.2 CAT12下载与MATLAB环境准备

下载要去SPM官网和CAT官网拿官方压缩包。网上有些二次打包的“一键版”，带着旧版本或改过的路径配置，出问题后排查起来非常痛苦。MATLAB版本不要用太老的，SPM12对版本兼容性有明确要求，装好后先确认能正常启动再继续。

注意：SPM12和CAT12的安装路径不要出现中文和空格，这是MATLAB环境下最容易被忽视的坑。

我见过太多人把工具包放进 D:\研究\工具\cat12 这种路径，结果某个子函数读不了路径，报错还看不出来原因。换成 D:\tools\spm12 和 D:\tools\cat12 之后一次通过。目录规划建议是这样的：

• 工具目录： D:\tools\spm12 、 D:\tools\cat12
• 数据目录： D:\data\adni_t1 、 D:\data\stats
• 输出目录： D:\data\output

% init_spm_cat.m —— 每次启动MATLAB先运行这一段
% 把SPM12和CAT12挂到MATLAB搜索路径
addpath(genpath('D:\tools\spm12'));
addpath(genpath('D:\tools\cat12'));

spm('defaults', 'PET');    % 借用SPM的默认参数框架，PET/EPI均可
spm_jobman('initcfg');     % 初始化作业管理器，批处理必需

which cat12                % 能打印出cat12.m的完整路径就说明加载成功

逻辑说明： addpath(genpath(...)) 把工具包下的所有子目录都加进搜索路径，CAT12依赖SPM的多个子函数，这一步漏了会在后面报 Undefined function 。 spm('defaults', 'PET') 是为了让SPM进入统计参数图流程的默认设置，用PET还是EPI在这个环节没有实质差别。 spm_jobman('initcfg') 是批处理前必须执行的一步，跳过它运行 spm_jobman('run', ...) 会直接报错。最后 which cat12 是一个验证动作：如果输出为空，多半是CAT12压缩包没解压完整，或解压后多了一层嵌套目录，检查一下文件结构即可。

2.3 跑通第一个VBM示例数据

CAT12压缩包里的 data 目录会带几幅示例T1像，专门用来验证环境。别一上来就处理自己的数据，先把示例跑通，能排除掉九成环境问题。

% run_first_vbm.m —— 用示例T1验证CAT12分割链路
ex_t1 = spm_select(1, 'image', '选择示例T1');   % 弹窗选一幅示例T1
matlabbatch{1}.spm.tools.cat.estwrite.data = {ex_t1};
matlabbatch{1}.spm.tools.cat.estwrite.output.mod = 1;  % 输出调制后灰质图
spm_jobman('run', matlabbatch);

跑完后数据目录里会多出 c1 、 c2 、 c3 、 mwc1 等前缀的NIfTI文件。 mwc1 是调制后的灰质图，是后续统计分析用的主料。这一条跑通，说明MATLAB、SPM12、CAT12三者兼容性没问题，可以放心上自己的批量数据。

如果你的机器上没有MATLAB图形界面，也可以改成从命令行传参的方式。批处理脚本本身不依赖GUI，只有 spm_select 这一步需要界面交互；换成代码里直接写死文件路径，整个流程就能丢到服务器上跑。

3. 从T1到灰质体积图：CAT12批处理流程与核心参数

VBM的中间链路可以拆成四步：格式转换、数据检查、分割配准、平滑。每一步都有参数要调，VBM项目翻车，十有八九是这四步里的参数组合出了问题。下面按顺序把每一步的落地细节讲清楚。

3.1 数据准备：格式、命名和第一批次检查

医院拿到的基本都是DICOM格式，VBM流程需要的是NIfTI，转换工具用dcm2niix，免费跨平台。转换之前把每个被试的DICOM单独放到一个文件夹，避免同一个被试的多个序列混在一起。

# 批量转换：把sub001的DICOM目录转成NIfTI
# -z y 输出nii.gz压缩格式，SPM12可以直接读
# -f sub%03d_T1w 按“sub+三位序号”命名
# -o 指定输出目录
dcm2niix -z y -f sub%03d_T1w -o /data/nifti /data/dicom/sub001

参数说明： -z y 把输出压成 .nii.gz ，SPM12和CAT12都能直接读，不要担心兼容性。 -f 控制输出文件名， sub%03d 会生成 sub001_T1w.nii 这种带补零的编号。 -o 后面是输出目录，建议全部被试统一放在同一个 nifti 目录下，后面批处理只需要指一个文件夹。

命名这件事值得单独拎出来强调：三个被试如果叫 sub1 、 sub01 、 sub001 ，输出前缀各不相同，统计建模时稍不注意就会对错位。更危险的是重名文件，两个被试的T1都叫 T1.nii ，分割输出直接互相覆盖，跑完才发现人数少了。统一用 sub001_T1w.nii 这种格式，全局唯一。

第一批次检查要做的是看每幅T1的原始质量：有没有头动重影、金属伪影、大面积灰白质对比度丢失。CAT12之后会出正式的QA报告，但开始前用肉眼粗筛一遍，能把明显坏掉的数据提前剔掉，省得它们在批处理中途报错中断整个队列。

3.2 批量分割的核心参数：bias、采样与调制

CAT12的分割在GUI里叫Segment，批处理对应 estwrite 这个模块。研究组数据多时，把文件名写进cell数组一次跑完。

% batch_segment_cat12.m —— 批量分割配准所有T1
t1_files = {
    '/data/nifti/sub001_T1w.nii'
    '/data/nifti/sub002_T1w.nii'
    '/data/nifti/sub003_T1w.nii'
    % 继续往下加被试
};

matlabbatch{1}.spm.tools.cat.estwrite.data = t1_files;
matlabbatch{1}.spm.tools.cat.estwrite.opts.biasreg = 0.001;  % 偏置场正则化
matlabbatch{1}.spm.tools.cat.estwrite.opts.biasfwhm = 60;    % 偏置场平滑核(mm)
matlabbatch{1}.spm.tools.cat.estwrite.opts.samp = 3;         % 采样距离(mm)
matlabbatch{1}.spm.tools.cat.estwrite.output.mod = 1;        % 1=调制(体积)，0=不调制(密度)
spm_jobman('run', matlabbatch);

参数说明： biasreg 控制偏置场校正的强度，MRI信号不均匀时分割会出错，CAT12用这个参数做正则化约束。数值越小校正越激进，但它是个“按下葫芦浮起瓢”的参数，调太小会把真实组织对比度也抹掉，所以优先用默认值，只有QA报告明确提示偏置场问题时才动。 biasfwhm 是偏置场估计时的平滑核，60mm是覆盖常见硬件不均匀尺度的经验值。 samp 是分割时对图像做各向同性采样的距离，3mm是精度和耗时的折中点，数据质量差或重点脑区小时可以用2，但耗时可能翻倍。 mod 决定是否做调制，这是VBM最核心的概念：配准到MNI空间时，局部体积会被拉伸或压缩，调制就是把这种形变信息乘回去。开调制得到的是灰质体积，不开得到的是灰质密度。绝大多数VBM论文说的是体积，所以 mod=1 是默认选择。

3.3 平滑核：什么时候用8mm，什么时候换6或12

分割调制后的图像还不能直接进统计，需要先平滑。平滑有三个目的：提高信噪比、让体素数据更接近正态分布、抵消配准后的残余个体差异。平滑核太小，噪声压不住，统计结果碎片化；平滑核太大，小的结构被直接抹掉，等于白扫。

% smooth_gm.m —— 对调制后灰质图做8mm平滑
gm_imgs = spm_select('FPList', '/data/nifti', '^mwc1.*\.nii$'); % 选所有调制灰质图
matlabbatch{1}.spm.spatial.smooth.data = cellstr(gm_imgs);
matlabbatch{1}.spm.spatial.smooth.fwhm = [8 8 8];   % 各向同性8mm高斯核
matlabbatch{1}.spm.spatial.smooth.dtype = 0;        % 保持原数据类型
spm_jobman('run', matlabbatch);

fwhm = [8 8 8] 表示三个方向都用8mm半高全宽的高斯核，这是VBM的经典参数。平滑后的文件前缀是 smwc1 ，统计建模时认准这个前缀就不会拿错文件。如果研究的是小核团，比如蓝斑、杏仁核亚区，可以把平滑核降到6mm甚至4mm，但前提是配准精度足够好，否则残留误差会被凸显。如果样本量小、个体差异大，10到12mm的平滑反而能帮你找到更稳健的宽泛差异区域，代价是定位精度下降。这个取舍要在论文方法里写清楚。

4. 组间比较怎么做：两种研究设计、协变量与统计阈值

分割和平滑做完后，每个被试得到一张 smwc1 灰质图，接下来的问题是：两组人的灰质到底差在哪。VBM最常用的是双样本t检验，配对设计用配对t或ANOVA也能走同一套SPM框架。这一章的决策直接影响结果能不能发表。

4.1 双样本t检验的基本假设与协变量选择

VBM的统计本质是逐体素的线性模型，核心自变量是组别，但T1衍生的体积数据几乎总受年龄、性别、颅内总体积（TIV）影响。年龄必须放进模型：全脑灰质随年龄自然下降，如果两组年龄有差异，不控制年龄直接就是把年龄效应误读成疾病效应。TIV也要放：头大的人灰质绝对体积自然大，两组如果平均头围不同，局部灰质萎缩的差异会被淹没在全局尺度差异里。性别同理。

协变量数量不是越多越好。样本量只有三四十人时，塞进五六个协变量，模型自由度会被吃光，统计功效大幅下降。我的经验是：年龄、性别、TIV三个是标配，其他如教育年限、病程等只在假设明确时加入，并且要用独立样本检验确认它们确实与因变量相关。

4.2 统计批处理脚本与结果读出

% stats_ttest.m —— 病人组 vs 对照组，含年龄性别TIV协变量
% 先读入临床信息表，保证行顺序与分割时的图像顺序一致
demo = readtable('/data/demo.csv');  % 列：id, age, gender, tiv, group
age = demo.age; gender = demo.gender; tiv = demo.tiv;

% 用文件名列表构建两组输入，顺序与协变量一一对应
all_files = cellstr(spm_select('FPList', '/data/nifti', '^smwc1.*\.nii$'));
patient_smwc = all_files(1:32);   % 假设前32个被试是病人组
control_smwc = all_files(33:64);  % 后32个是对照组

matlabbatch{1}.spm.stats.factorial_design.des.t2.scans1 = patient_smwc;
matlabbatch{1}.spm.stats.factorial_design.des.t2.scans2 = control_smwc;

matlabbatch{1}.spm.stats.factorial_design.cov(1).c = age;
matlabbatch{1}.spm.stats.factorial_design.cov(1).cname = 'age';
matlabbatch{1}.spm.stats.factorial_design.cov(1).iCC = 1;   % 全局中心化
matlabbatch{1}.spm.stats.factorial_design.dir = {'/data/stats/ttest_AD'};
spm_jobman('run', matlabbatch);

% 估计参数：SPM把基于体素的统计都归到fmri_est，VBM也走这个
matlabbatch{2}.spm.stats.fmri_est.spmmat = {'/data/stats/ttest_AD/SPM.mat'};
spm_jobman('run', matlabbatch);

% 定义对比：病人组 < 对照组
matlabbatch{3}.spm.stats.con.spmmat = {'/data/stats/ttest_AD/SPM.mat'};
matlabbatch{3}.spm.stats.con.consess{1}.tcon.name = 'AD < NC';
matlabbatch{3}.spm.stats.con.consess{1}.tcon.convec = [1 -1 0 0 0];
spm_jobman('run', matlabbatch);

convec 的长度是“两组均值 + 协变量数量”，这里放了年龄、性别、TIV三个协变量，所以是 [1 -1 0 0 0] 。前两个数字 1 -1 表示病人组均值减对照组均值，后三个0表示协变量不参与对比。如果要看反方向，写成 [-1 1 0 0 0] 即可。

结果读出在SPM的Results按钮里，选中SPM.mat，指定对比和阈值后就能看到叠加在模板上的显著簇。值得注意的是SPM把VBM统计归在 fmri_est 模型里，新手看到 fmri 字段会疑惑，其实这里只是SPM内部的模型估计入口，与功能磁共振无关。

4.3 多重比较校正：FWE、FDR与TFCE怎么选

几十万个体素同时检验，不做多重比较校正等于主动生产假阳性。SPM里常用的有三档：体素级FWE要求peak的p<0.05，最严格，样本量不够时往往一个点也剩不下；FDR温和一些，适合探索性研究；cluster级FWE是现在论文里的主流做法——先设一个体素水平的筛选阈值，通常用未校正p<0.001，再对连成片的体素簇做FWE校正。

FSL的randomise里还有TFCE（基于阈值无簇增强），不需要预先设定体素阈值，置换检验直接给每个体素打分，对小效应且分布弥散的差异更敏感。务实建议是：新研究、样本量不大，用体素未校正p<0.001筛簇并报cluster FWE；如果结果一个簇都过不了，再看FDR，但要在论文里明确写“本结果经FDR校正”。同一批数据反复换阈值直到显著，这是审稿人最反感的p-hacking，千万别干。

5. VBM避坑指南：这五个问题我几乎在每批数据里都见过

VBM流程链长，前面任何环节的错误都会在最后一环爆发。以下五条是真实数据上反复出现的问题，每条按现象、原因、解决的顺序写，遇到类似情况可以直接对照排查。

5.1 现象：QA报告大面积红叉、分割严重错位

CAT12跑完会生成QA报告（check data quality），用颜色标出每例的分割质量，红色代表严重问题。最常见的诱因是头动伪影：被试扫描时轻微移动几毫米，T1边缘就出现重影，分割时脑脊液被分进灰质，概率图到处是斑块。另一类是偏置场严重，信号不均匀导致灰质概率图呈奇怪的梯度分布。

原因：头动和线圈不均匀是硬件层面的问题，参数调整只能掩盖，不能根治。解决：先看QA给出的exclude建议，如果超过10%的被试被标记，就要在论文里报告排除标准，并重新采集这些被试。轻微偏置场问题可以把 biasreg 从0.001调到0.0005重跑，但大的头动伪影必须走重采流程，不要指望算法能修复。

5.2 现象：跑完发现输出文件被覆盖，统计时数量对不上

批处理几十个被试后建模，报错说某个smwc文件不存在，或者统计人数比预想少了好几个。原因十有八九是输入命名不唯一：两个被试的T1都叫 T1.nii ，或者 sub001_T1w.nii 与 sub01_T1w.nii 生成的前缀极易混淆。

解决：DICOM转换阶段就用全局唯一数字编号，并加一道清点防线：

# 检查smwc1文件数量是否与被试数一致
ls /data/nifti/smwc1*.nii | wc -l

# 如果少了，列出已经跑完的被试编号，和原始名单比对
ls /data/nifti/smwc1*.nii | sed 's/smwc1//; s/.nii//' | sort > /tmp/done.txt

第一条命令输出文件总数，第二条生成已完成编号清单。数量对不上时，先查命名，再翻批处理日志看是否有被试分割到一半报错退出。CAT12跑长队列时偶尔会因为内存不足中断某个被试，日志里会有 Out of memory 记录，这种重跑单例就行。

5.3 现象：统计结果一个显著簇都没有

双样本t检验全脑空荡荡，这是VBM里最常见的挫败时刻。先查三件事：调制开没开、协变量有没有放TIV、平滑核是不是太大。

如果 estwrite 里设了 mod = 0 ，实际上比的是灰质密度而不是体积，配准形变会把真实的局部体积差异吸收掉，结果自然一片空白。TIV不放进去，组间头围差异会消耗统计功效，等于在模型里留了一个巨大的噪声源。平滑核从8mm改到16mm，局部萎缩效应被抹平，也找不到簇。还有一个容易踩的坑：有人把未调制的 c1 图当成输入，正确输入是 smwc1 ，也就是“平滑+调制+灰质”产物的组合。改完这些再跑，先用未校正p<0.005看效应分布，确认差异脑区符合解剖直觉，再把阈值收紧到p<0.001。

5.4 现象：SPM估计时设计矩阵提示数据不足

跑 fmri_est 时报“设计矩阵不可逆”或“InSufficient data”。这类报错通常和样本量、协变量数量、数据完整性三个因素有关。总人数少于“组数+协变量数”时，模型数学上根本估计不出来；协变量里混入缺失值NaN，矩阵直接变成奇异矩阵。

另一个隐蔽原因是图像列表与协变量行序错位：统计脚本里写的第5个smwc文件和demo表格第5行的年龄不是同一个人，模型拟合出的方差会异常大，间接导致同样的报错。

解决：确认被试数至少是协变量数加组数的三倍以上；把图像列表与协变量矩阵按被试编号排序后打印出来逐一核对，笨办法最可靠。如果个别被试分割失败没有输出smwc文件，记得从两个列表同时删掉，不要让长度错位。协变量里绝对不能有NaN，用 demo = rmmissing(demo) 提前清理。

5.5 现象：peak坐标落在脑室或白质，和文献对不上

结果图上读出的peak坐标在脑室或白质里，和文献已知的灰质萎缩区域对不上。原因通常是配准失败：个别被试有脑畸形或大片白质病变，DARTEL把结构拉错了位；另一个原因是阈值太低，出来的小簇本身就是噪声簇。

解决：把统计图叠加到组平均模板上，肉眼检查每个报告的簇。SPM的Results界面里打开“明显示簇范围”，确认peak体素确实落在灰质结构内。写论文时把坐标旁标注模板来源（MNI152或ICBM），审稿人追问时才能说清楚。配准失败的个例，直接考虑从分析中剔除，并在排除流程里写明依据。

6. 收尾：结果可信不可信，先过这三关再投稿

6.1 先过质量关：DQI值和exclude名单

CAT12的QA报告不只看彩色图，要盯三个数字：IQR（图像质量等级）、DQI平均值、exclude标记。我自己的习惯是跑完批处理就把QA报告全部导出，凡是被标记 exclude 或两个以上指标亮红的被试，直接列进论文排除名单，并在方法里写一句“排除标准基于CAT12的DQI异常值”。这样审稿人不会在这个环节纠缠，因为你有明确、可复现的排除依据。

6.2 再过配准关：肉眼检查每一例

配准质量是VBM最大的黑匣子。批处理跑完绝不意味着每颗脑袋都对齐得很好。我一般会把每例的 smwc1 图叠加在模板上，逐层看矢状位和轴状位，重点观察侧脑室边缘、扣带回、小脑扁桃体三个结构有没有错位。三十个被试大约二十分钟，这二十分钟能挡住绝大多数翻车。如果发核心期刊，这一步的检查截图甚至可以放进补充材料，审稿观感会好很多。

6.3 最后过效应量关：别只报p值

显著结果至少要补一个效应量，否则审稿人一句“效果是否小到无意义”就能让返修来回折腾。两样本t检验在两组样本量接近时，可以用峰值t值近似估算Cohen's d：

t_val = 4.32;         % 从SPM结果表抄的峰值t
n1 = 32; n2 = 32;     % 两组样本量
d = 2 * t_val / sqrt(n1 + n2 - 2);  % 近似Cohen's d
fprintf('Cohen''s d ≈ %.2f\n', d);

这个公式只适合两组样本量相等或接近的场景，正式报告建议用G*Power或R的 effectsize 包算精确值。最后一个习惯是：每次跑完一批数据，把QA报告、配准检查记录、统计结果表归档在同一个项目文件夹里，文件名带上日期。VBM这个流程我已经跑过几十批数据，下载、分割、统计每个环节都有各自的玄学，但把QA、配准、效应量这三关守住，结果基本站得住。希望帮到你。

本文还有配套的精品资源，点击获取