SPM中文教程汇总已

1、、SPM的安装与启动先安装matlab,然后将SPM复制到matlab下的一个文件夹（SPM2需要matlab6.0或以上版本）。 启动matlab,首先setpath,然后在matlab命令窗口中输入SPM即可启动，然后选择fMRI，也可 以直接输入SPMfMRI二、SPM数据处理概要先将所得数据进行 空间预处理（对齐，平滑，标准化等），然后进行模型估计（将刺激的时间、 间隔与血流动力函数进行卷积，所得结果与全脑象素信号进行相关分析），最后察看结果。三、SPM8数据处理的一般步骤为方便后续的数据处理，如果数据分散处理后整合，建议所有处理数据路径保持一致， 要统一 路径。处理前首先要采用数据转

2、换软件将 dicom数据转换成SPM解析格式，然后进行数据预处理， 预处理结束后到matlab安装目录中备份spm*ps文件，其中包含了空间校正和标准化的信息，然 后进行建模分析。运行命令：spmfmri，打开spm8的操作界面我们称左上侧的窗口为按钮窗口 （buttonwindow），左下侧的窗口为输入窗口 （inputwindow），右 侧大窗口为树形结构窗口或图形窗口 （TreeBuildingWindoworthegraphicswindow）。在spm8和spm5中，每一步处理都采用了直观的 树形结构”的面板，如果一个分支项左面有“+” 号，你可以双击显示子分支项，？如果一个分支项右

3、面有“X”号，你必须为之指定选项（否则不能 运行该tree），分支项的选项在其右侧面板指定，而帮助信息则在下面的面板中显示。如果我们处 理数据没有特殊需求，我们只关心带有up):1:2: nslices2:2: nslices，如1:2:25,2:2,25 in terleaved（top-dow n）: nslices:-2:1, nslices-1:-2:1ReferenceSlice,我们输入“31。选择参考扫描层（一般可使用默认值），其它扫描层的起始时间 都将以此层的起始时间为标准来移动进行校正。通常选择nslice/2，如25层时选择13层作为参考层。a*img文件，这Filenam

4、ePrefix，是指新生成的图像前加何标记，一般采用默认设置。 最后点击面板上方的向右的绿色三角即开始运行。运行完后将会生成一系列就是时间校准后的数据。注意很多研究者容易将时间校准和空间校准顺序颠倒，一般的观点是如果图像获取是隔层（interleaved）进行的，如 1、3、5、7、9、2、4、6、8、10，则要先进行 slicetiming 再进行 realign，如果图像各层是连续（sequential获取的，女口 1、2、3、4、5、6、7、8、9、10,则要先进行realign再做slicetiming。（为什么？）3、Realign （相当于 AFNI 中的 registration

5、）分两步：1）coregister，将每个session的第一个scan与第一个session的第一个scan进行比较，然后将每 个session中的其他scan与本session中的第一个scan进行比较，得到每个 filename.img文件的转: I;换参数，生成 filename.mat文件，同时为每个 session生成一个对齐参数（realignmentparameters， 文件名为 realignment_params_*.txt2） reslice，用filename.mat文件对filename.img重新切片，生成rfilename.img文件。并可依选择生 成一个平均象

6、，名为meanfilename.img。以上两步具体解释如下一一即使我们对被试的头部做了很好的固定，在实验过程中，被试也会不由自主的有一些轻微的头动，这在fMRI实验中尤为明显。这一步就是把一个实验序列中的每一 帧图像都和这个序列的第一桢图像按照一定的算法做对齐，以矫正头动。做完这一步，能给出该序列中被试的头动情况，以作为是否放弃该数据的依据，如果头动超过1个voxel （功能图像扫描矩阵一般是64*64，则体素的大小为（FOV/64）*（ FOV/64）* （层厚+层间距），则要考虑放弃该 时间点数据。该程序利用最小二乘法（leastsquaresapproach原理和含6个参数（刚体模型）

7、的空间变 换，对从一个被试获取的时间序列进行校正。用户可指定某个volume作为随后volumes的参考。可以是第1个volume，也可选择比较有代表性的 volume（更明智的选择），例如选择磁场相对稳定的第 4个volume。校正信息(头动信息)将在结果窗口(GraphicsWindow)显示。每个Session的校正 信息将存储为rp*txt，其中*为Session数据集名称。另外，头动校正信息将以plot图形显示。如下 图：Translation表示被试头部在X，丫，Z三个方向的平移，分别用蓝，绿，红三种颜色表示。Rotation 表示被试头部在实验过程当中绕 X(L-R)，丫(A-P

8、)，Z(S-I)三条轴的转动角度。横坐标代表这个序列 所采集的所有图像，纵坐标表示的是偏移量和偏转角度，分别以毫米和度为单位。采用SPM8，头动信息和空间标准化的图形文件将以spm“data” .p的形式保存于matlab的工作目录下，如我们是2009年4月30日处理的数据，则将以 spm_2009Apr30.ps文件存于 matlab的work目录下。我们在预处理面板校准选项中选择“Realign(Est&Res),”出现如下对话框，我们按下面设置进行：选中“data选择“NewSession然后选中data下出现的“Sessior选项。点击 “SpecifyFiles,用spm文件选择器选

9、择刚做完时间校准的图像(a*.img )。其余选项采用默认设置，点击上方绿色的 三角开始运行。Realign这一步也有分开进行的，具体描述如下3.1Realign:Estimate (重排参数的估计)此步骤采用最小方差原理和六参数刚体空间变换来重排从同一个被试上采集到的图像数据。使用者 所选取的第一幅图像文件将被作为其它图像重排的参考标准。也就是说，你想要用哪一幅图像作为参考标准，就先选哪一幅图像的文件。参考图像文件不一定非用采集到的第一幅图像，使用最有代 表性”的一幅图像也许更好。本步骤的目的 主要是去除fMRI和PET数据中的运动伪影。图像数据 的头文件会被改写以反映数据相对空间位置的变化

10、。此过程的具体参数会在结果窗口中以平移I/ j(translation)和旋转(rotation)曲线图显示。每个session的重排参数会被存储到名为 rp*.txt的文 件中。这些参数可以在最后的一般线性模型统计估计中作为混淆因素考虑进去。3.1.1Data选择一个被试需要进行此步骤处理的所有sessions注：在coregistration这一步，首先是对所有的session进行重排，其具体做法是把所选每个session的第一个scan与所选第一个session的第一个scan对齐。然后再把每个session里的其它scan与该session的第一个scan进行对齐。使 用此方式进行重排

11、是因为各个 session的数据之间可能会有较大差异。Session选择session里所有 的 scan=3.1.2Estimati on Optio ns这里包括各种注册参数选择项，若对某一个选项不确定，使用软件默认值 即可。Quality质量与速度的权衡。选择高质量以最慢的速度给出最精确的结果，低质量以较快的 速度给出较不精确的结果。此参数的设定实际影响到的是参与参数估计的象元(voxel)的数目。其依据是有些象元(voxel)其实对重排参数的估计贡献不大，可以舍弃。Separation此参数以毫米为单位，表示对参考图像文件进行重采样时采样点之间的间隔。采样点之间间隔越小，结果越精确，运

12、算速度越慢。Smooth in g（FWHM）高斯平滑的半高宽值。在估计重排参数之前一般先进行高 斯平滑。PET数据一般使用 7mm。MRI数据一般使用 5mm。NumPassesRegistertofirst:所有图像 文件对齐注册到第一幅图像。Registertomean:使用two pass处理将所有图像文件对齐注册到所有 图像文件的平均图像。PET数据一般注册到平均图像。因为 PET数据相比fMRI数据噪音更大， 文件更少，所以时间的影响更小。MRI数据一般注册到第一幅图像。虽然使用two pass处理可能 更精确，但是其对效果的提高与其所损失的运行时间相比得不偿失。Interpola

13、tion在估计最佳变换时对数据进行重采样的方法。高的degree提供更好的结果，但是也更慢，因为会采样更多的相邻象元（voxel） 52,53,54。Wrapping此参数指示一个 volume中数据 wraparoundin的方向（此处具 体理解有待大家补充）。Nowrapping：适用于PET数据或者已进行过空间变换的数据。同时当你不 确定自己数据类型时，推荐使用此选项。WrapinY :适用于没有重排（resilce）过的在丫方向上进行相位编码的MRI数据。Weighting提供一个加权图像，在估计重排参数时对参考图像的每一个 象元进行加权。加权系数与标准差成反比。例如当有大量额外的头动

14、（如说话或者特定区域内的严 重伪影）时。（此处具体理解有待大家补充）。3.2Realign:Reslice （据已估计出的参数重排）此功能重排以上步骤中已进行参数估计和注册的图像文件，使之与参考图像文件达到象元级的匹配精确。重排后的数据被命名为：r +原文件名。3.2.1Images选择要重排的数据文件3.2.2ResliceOptions各种重排参数设定，若对某一个选项不确定，使用软件默认值即可。 ReslicedimagesAlllmages（1. n）重排所有数据，包括标准参考图像（重排后还是保持原位置不变） 。I ./ JImages2.n重排除了标准参考图像之外的所有数据。此选项用于

15、当你 以MRI结构像为标准重排PET 图像数据，而又不想在结果中再生成一个等同的MRI标准结构像时。AllImages+MeanImage:重排图像文件之外，另生成一个重排后的平均图像文件。Mea nI mageO nly:只生成重排后的平均图像文件。Interpolation图像文件重采样和重写入的方式。NearestNeighbour：最快，但不推荐使用。Bilinearlnterpolation :可用于 PET 数据，但不是太适用于 fMRI 数据。Fourierlnterpolation:此选项 仅适用于纯刚体变换，也就是说象元大小必须是相同，并且等方性（正方体）的17,14 o W

16、rapping 此参数指示一个 volume中数据wraparoundin的方向（此处具体理解有待大家补充）。Nowrapping： 适用于PET数据或者已进行过空间变换的数据。同时当你不确定自己数据类型时，推荐使用此选 项。WrapinY:适用于没有重排（resilce）过的在丫方向上进行相位编码的 MRI数据。Masking因 为扫描过程中被试总会或多或少有头动，造成同一个时间系列数据里所采集到的图像的边界不会完 全重合。在有些图像还有数据的地方（信号值大于0）,其它一些图像已经超出了图像边界（信号值为0）了。在这些信号为0的区域是无法采样数据的，因此SPM只要检测到某一幅图像在某个区域已

17、经超出了边界(即信号为0),就会将其它所有图像的此区域信号值均设为0。此做法相当于取了时间系列数据中所有图像的交集。3.3Realig n:Estimate&Reslice将上述参数估计与数据重排合到一起做。全部选项与参数原理均与3.1和3.2中对应项相同。4、Normailze选用realign步骤中得到的平均象与模板进行比较，获得进行标准化的参数，参数文件命名为 filename_sn3d.ma，然后依据此参数文件对每个img文件进行标准化，生成文件nfilename.img。具 体操作如下：在预处理面板标准化选项中选择“Normalise:Estimate&Write 出现如下对话框：我

18、们做如下设置：选中“data”“newsbjec在data下新出现的“subjects项中作如下设置，“ sourceimag选择空间校准步骤中生成的 mean文件，“ imagetowrite选择所有刚进行完校准的文件 “ra*img，” “templateimag我们选择“EPI.nii，其余采用默认设置，点绿三角运行。I ? xX y x !5、SmoothFWHM推荐为象素大小的两至三倍。在预处理面板标准化选项中选择“smooth,”出现如下对话框：我们在“imagetosmooth选项中选择所有刚进行完标准化的文件“wra*.img；然后点绿三角运行即可。这里FWHM我们采用默认设置

19、“888。6 fMRImodels依据提示填入刺激出现的间隔与时间，并选择实验涉及类型，然后进行估计。估计结果生成spm.mat等文件，保留在当前工作目录。即以前版本的 “fMRImodel ” spm5和spm8的分析选项有 所变化，使用 Specify1stlevel做单 个被试(singlesubject)分析； 使用 Specify2ndlevel做组分析:-I:(groupanalysis)。比如我们选择 “Specify1stlevel,出现如下对话框(见下一页)：我们选择Directory指定一个文件夹存放结果数据，其余做如下设置：“Unitsfordesign选择 “Scans

20、， “Interscaninterval输入” “2， 选择 “DataandDesigri后选择 “NewSubject/Session再选择新出现的“Subject/Session选择“Sean并用文件选择器选择相应任 务的所有平滑后的功能图像(swra*.img)然后点击“dow”选择“Condition后选择“Newcondition， 然后选中新出现的“Condition，” “nam”选项输入任务条件的名称，“onse输入任务条件的启动向量， 代表任务刺激启动的扫描数，选中“Durations输入任务组块的持续时长，如是事件相关设计请输入“ 0”如还有其他任务，要再次选中 “ Co

21、ndition ”“ Newcondition ”“定义其他任务条件。设置完毕后点击绿三角运行。这样将会在开始选择的目录中生成文件spm.mato下面要估计我们刚建立的模型，在模型设置面板中点击“estimate”将打开如下对话框：很简单，我们只须选择刚生成的“spm.ma文件点击“dow6然后点击绿三角运行即可。估计完成后，我们选择“results,将打开如下对话框:上图设计矩阵表明一个扫描序列中我们有三个任务条件，均为事件相关设计，选中-contrasts，”点击 “ definenewcontras第一个任务 “ contras定义为 “ 1,第二个任务定义为 “ 01,” 第 三个任务

22、nn定义为“001 ”第一个任务减第三个任务定义为“101”，第二个任务减第三个任务定义为“01-1”，第二个任务减第一个任务定义为-11”，其余操作以及激活图显示和以前版本都是一致的。值得我们注意的是，我们指定或输入一系列的参数来进行每一步的处理， 完成后我们可以通过Save 按钮将每一步存为一个*.mat文件。以后我们可以通过Load按钮重新加载并使用这些*.mat文件，我 们适当修改后（例如改变所运行的数据集）再选择 Run按钮运行。另外我们可以使用TASKS菜单项指定一系列的操作（预处理和/或分析），在TASKS-Batch菜单 项下，你可以在一个大文件中指定数据处理所有步骤。 批处理

23、交互界面非常灵巧。它知道根据指定 的步骤将产生什么文件。例如，在 Tasks-Batch菜单项下，选择NewSpatial。在你的SPM任务树.1（Jobstree）中选择（highlight）-Spatial，并从选项面板中选择NewRealig n选项指定我们前述的 realig nme nt的详细步骤。现在，当我们向任务树中添加标准化步骤时，我们将看见为标准化操作选f n.择r*文件的选项，尽管事实上，我们还没有真的生成r*文件。所以，批处理非常智能，可以预测我们的需要。如果我们已经生成了几个批处理任务文件，我们可以使用 TASKS-Util-ExecuteBatchJobs选项来运行它

24、们。该工具允许我们选择一系列的 mat文件（不一定是 batchjobs文件）来运行。如果你怀念旧版本SPM操作界面，在TASKS菜单下，选择Sequential,这样 就不会出现树（tree），选项将会出现在SPM输入窗口。7、Result选中刚才生成的spm.mat文件，定义con strast,看结果。在上面的“fMRImodels ”部分有讲述。四、SPM的多种数据转换方法使用SPM进行数据处理前，必须先将其它档案格式转换成 SPM可以读取的Analyze档案格式，包 含.img档和.hdr标头档，相关的转档软件有 XMedCon和MRIcro。1利用AFNI数据转换首先使用 AFNI

25、的三维数据重建：to3d-time:tz177202saltplus*生成+orig文件，然后:3dAFNItoAnalyze-4D-orientLPI*epi+orig，将生成*.hdr与*.img文件.（*代表你所用的任意文件名）， 然后打开MRIcro软件：选择要转换的hdr,img文件,processing,ok择面板上部 File-Saveas4dto3d-Saveasintel-Save这里有点忘了 :）（这一段我压根没有看懂）2直接用MRIcro转换:（此部分转自核医学论坛）注:SPM的dicomimport是可以导入dicom文件，但导入的这些图像还是 2D的，不是3D的，不能 用