SPM8处理流程

1、SPM8处理磁共振实验数据说明一、数据准备为方便后续的数据处理，如果数据分散处理后整合，建议所有处理数据路径 保持一致，要统一路径。处理前首先要采用数据转换软件将dicom数据转换成SPM解析格式，转换时格式请选择 NIfTI，可用SPM俞入面板中的DiCOMImport 模块转换，也可以采用专门的转换软件，如MRIcovert。然后进行数据预处理，预处理结束后到matlab安装目录中备份spm*ps文件，其中包含了空间校正和 标准化的信息，然后进行建模分析。二、数据处理流程Spm8和以往版本一样，先安装 matlab7.1以上版本，然后将spm8文件夹放 在matlab安装目录中，我们一般放

2、在toolbox文件夹中，然后启动matlab用菜 单中设置路径命令增加spm8的路径并保存。之后我们运行命令：spm fmri，这 样将打开spm8的操作界面，左上侧的窗口为按钮窗口(button win dow)，左下侧的窗口为输入窗口(input window)，右侧大窗口为树形结构窗口或图形窗口 (Tree Building Window or the graphics win dow) 。在spm8和spm5中，每一步处理都采用了直观的“树形结构”的面板，如果一个分支项左面有“ +”号，你可以双击显示子分支项，如果一个分支项右面有“<X”号，你必须为之指定选项(否则不能运行该t

3、ree)，分支项的选项在其右侧面板指定，而帮助信息则在下面的面板中显示。如果我们处理数据没有特殊 需求，我们只关心带有“ <X”项目并完成输入即可，其余均可采用默认设置。另外注意在Tree Building Window 的顶部菜单，新增了一个菜单项“ TASKS，在使用批处理分析时非常重要。对上图右侧选项我们做如下设置：1: Slice TimingSlice Timimg用来校正1个volume中层与层之间获取(采集)时间的差异， 对事件相关设计的实验尤为重要。我们在按钮窗口中的预处理面板中点击“ SliceTimimg”，将出现如下对话框:Data，点击data并在下面的面板中点击

4、“ newsession ” ,这样在data下会 出现“ session ”的分支项，选中该项并点击面板下方的“select files ”，然后用spm文件选择器选择你要处理的数据，最后点击“ down”。选择数据时可以 把静息态、数值任务和物理大小任务分为三个session来选(datanewsessionsession )，也可以作为一个 session 来选，结果是一样的。Number of Slices ，我们输入每祯图像的层数，如“32”TR输入重复时间，一般为2秒，输入“ 2”TA是每祯图像获取第一层开始到获取最后一层图像的时间间隔，我们输入TR-TR/nslice，可直接输入

5、公式，如我们输入“ 2-2/32 ”Slice order ，输入“ 1:2:31,2:2:32”。指定层获取顺序的层次序参数是一个含N个数的向量，这里N是每个volume所含的层数。每一个数表示该层在图像(volume)中的位置。向量内的数字排列顺序是这些层的获取时间顺序。 如行 向量1 3 5 7 9 1113 15 17 19 21 23 25 2 4 6 8 10 12 14 16 18 20 22 24(在Matlab 中可表示为1:2:25,2:2,25)Referenee Slice,我们输入“ 31”。选择参考层，通常选择 n slice/2 :如25层时选择13层作为参考层。

6、File name Prefix ，是指新生成的图像前加何标记，一般采用默认设置。最后点击面板上方的向右的绿色三角即开始运行。运行完后将会生成一系列a*img文件，这就是时间校准后的数据2： Realignment即使我们对被试的头部做了很好的固定，在实验过程中，被试也会不由自主的有一些轻微的头动，这在fMRI实验中尤为明显。这一步就是把一个实验序列 中的每一帧图像都和这个序列的第一桢图像按照一定的算法做对齐，以矫正头 动。做完这一步，能给出该序列中被试的头动情况，以作为是否放弃该数据的依据，如果头动超过1个voxel (功能图像扫描矩阵一般是 64*64，则体素的大小 为(FOV/64)*(

7、 FOV/64)* (层厚+层间距)，则要考虑放弃该时间点数据。 另外，头动校正信息将以plot图形显示。如下图：Tran slation 表示被试头部 在X，丫，Z三个方向的平移，分别用蓝，绿，红三种颜色表示。Rotation表示被试头部在实验过程当中绕X(L-R)，丫(A-P)，Z(S-I)三条轴的转动角度。横坐标代表这个序列所采集的所有图像，纵坐标表示的是偏移量和偏转角度，分别以毫米和度为单位。采用SPM8头动信息和空间标准化的图形文件将以spm_“ data ” .ps的形式保存于matlab的工作目录下，如我们是2009年4月30日处理的数据，则将以spm_2009Apr30.ps

8、文件存于 matlab的work目录下。translationrotation在预处理面板校准选项中选择“ Realign (Est & Res)”，出现如下对话框, 我们按下面设置进行：Bat elk. Kdi iorgdj t SPIA Ea=icIO Toy ExanplaUd U HModule ListData|Help on. Realiqn: Estimat白 & ResliceJbData<-XEstimation OptionsQuality0.9Separation4Smoothing (FWHM)5Num Passes.r to meanInterp

9、olation. Et-Spl neVWappingNo jrapWeightingRe slice OptionsRediced images.an ImageInterp elation. B-SplineWrappingNo vvrapzl匚urrenE Module: IHeaMgn: Eslimate 芯 R&s-lieeCurncnl. Ilcm: DislaEdWHAdd new sessions for this subject. In the co registration step, tfie sessions are tirst reali gned to eac

10、h ottier, by aligning ttie first scan from each session to the tirst scan or ttie irst session. Then ttie images within each sess on are aligned to the first image of ttie session The oarameter estimation >s performed this vw because it is assumed选中“data” ,选择“NewSession” ,然后选中data下出现的 “Session 选项

11、。点击“ Specify Files ”，用spm文件选择器选择刚做完时间校准的图像(a*img )。其余选项采用默认设置，点击上方绿色的三角开始运行。3： Normalise在预处理面板标准化选项中选择“ Normalise: Estimate & Write ”，出现如下对话框：|Help on- Normalise: Estimate & WriteData<-XEstimation OptionsTemplate Image<-XTemplate Weighting Image0 filesSource Image Smoothing8Template Im

12、age Smoothing0Affine Regularisation.templateNonlinear Frequency Cutoff25Nonlinear Iterations16Nonlinear Regularisation1Writing OptionsPreserve.entrationsBounding box2x3 doubled 口1刘Le £dit ( Sfi £asitIO 0jr Erampl? Td a；ra>Cuirent Module: Normalise: Estimate I WiiteNormalise' Estifl*

13、J- Batch Edi torModule ListCurrent Item; DataEdit 畑 e |DataList of subjects. Images of each subject should be warped differently1 or more options must be selected from:* Subject_d如下设置：选中“ data ”“ newsbject ”，在 data 下新出现的 “subject ” 选项中作如下设置，“source image”选择空间校准步骤中生成的 mean文件，“image to write ”选择所有刚进行完

14、校准的文件“ ra*.img ”，“ template image ”我 们选择“ EPI. nii ”，其余采用默认设置，点绿三角运行。4： Smoothing在预处理面板标准化选项中选择“ smooth”，出现如下对话框：在“ image to smooth”选项中选择所有刚进行完标准化的文件“ wra*img 然后点绿三角运行即可。这里 FWH我们采用默认设置“ 8 8 8 ”。5: Specify 1st level即以前版本的“ fMRI model ”，spm5和spm8种的分析选项有所变化,使用Specify 1st level 做单个被试(single subject)分析；使

15、用 Specify 2nd level 做组分析(group analysis)。选择“ Specify 1st level ”，出现如下对话框：Ea.tch Idi tar-|n| x|Edi t SPM fi asi eliQ Toy Exople：Module ListCuPfertt Module: IMRI moddl specificalionfMRI model spefion fMRI model specificationDirectoryTiming parametersUnih for desigrInterScan intervalMicrotime resolutio

16、nMicrotime onsetData & DesignFactorial designBasis FunctionsCanonical HRFModel derivativesModel Interactions (Volterra) Global norrnalisationItEHn: Directoru:116<1No derivatives el InteractionsNoneDirectory3Select a directory where the SPM.mat file design matrix will be 轲社en.containing the sp

17、ecifieddSelect Fies选择Directory指定一个文件夹存放结果数据，其余做如下设置：“ Un its for desig n ” 选择“ Sea ns”， “In tersca n in terval ” 输入“ 2”， 选择“ Data and Design ”后选择“ New Subject/Session ”，再选择新出现的 “Subject/Session ”，选择“Scans”并用文件选择器选择相应任务的所有平滑 后的功能图像(swra*.img )然后点击“ down'，选择“ Con diti on ”后选择“ Newcondition ”，然后选中新

18、出现的“ Condition ”，“namd'选项输入任务条件的 名称，“ on set ”输入任务条件的启动向量，代表任务刺激启动的扫描数，选中Durations ”输入任务组块的持续时长，如是事件相关设计请输入“0”。如还有其他任务，要再次选中 “ Condition ”“Newcondition ”“Condition定义其他任务条件。设置完毕后点击绿三角运行。这样将会在开始选择的目录中生成文件spm.mat.下面要估计刚建立的模型，在模型设置面板中点击“estimate ”，将打开如下对话框：Jnl >dFileJ”®矽 B&sidO Toy E:caj

19、iipl«Module ListId q Help on: Model estimationJSelect SPMmat<-xMethodClassicalCurrent Item: Select SPM.matModel estimaticlCurrent Module: Model estimationSelMtFile* |Select SPM.matSelect the SPM.mat file that contains the design specification.The directory containing this file is known as th

20、e input directory只须选择刚生成的“ spm.mat”文件点击“ down"然后点击绿三角运行即可。 估计完成后，选择“results ”，将打开如下对话框：上图设计矩阵表明一个扫描序列中我们有三个任务条件，均为事件相关设 计，选中“ t con trasts ”，点击“ defi ne n ewe on trast ”，第一个任务“ con trast ” 定义为“ T,第二个任务定义为“ 0 1 ”，第三个任务nn定义为“ 0 0 1 ”， 第一个任务减第三个任务定义为“ 1 0 -1 ”，第二个任务减第三个任务定义为“ 0 1 -1 ”，第二个任务减第一个任务定义为“ -1 1”，其余操作以及激活图显示和 以前版本