SPM中文教程 汇总(已整理) 联系客服

发布时间 : 星期六 文章SPM中文教程 汇总(已整理)更新完毕开始阅读4e5795aa844769eae109ed2e

件将以spm_“data”.ps的形式保存于matlab的工作目录下,如我们是2009年4月30日处理的数据,则将以spm_2009Apr30.ps文件存于matlab的work目录下。

我们在预处理面板校准选项中选择“Realign (Est & Res)”,出现如下对话框,我们按下面设置进行:

选中“data”, 选择“New Session”, 然后选中data下出现的 “Session”选项。点击“Specify Files”,用spm文件选择器选择刚做完时间校准的图像(a*.img)。其余选项采用默认设置,点击上方绿色的三角开始运行。

Realign这一步也有分开进行的,具体描述如下——

5

3.1 Realign: Estimate(重排参数的估计)

此步骤采用最小方差原理和六参数刚体空间变换来重排从同一个被试上采集到的图像数据。使用者所选取的第一幅图像文件将被作为其它图像重排的参考标准。也就是说,你想要用哪一幅图像作为参考标准,就先选哪一幅图像的文件。参考图像文件不一定非用采集到的第一幅图像,使用最有“代表性”的一幅图像也许更好。本步骤的目的主要是去除fMRI和PET数据中的运动伪影。图像数据的头文件会被改写以反映数据相对空间位置的变化。此过程的具体参数会在结果窗口中以平移(translation)和旋转(rotation)曲线图显示。每个session的重排参数会被存储到名为 rp*.txt 的文件中。这些参数可以在最后的一般线性模型统计估计中作为混淆因素考虑进去。

3.1.1Data 选择一个被试需要进行此步骤处理的所有sessions。注:在 coregistration 这一步,首先是对所有的session进行重排,其具体做法是把所选每个session的第一个scan与所选第一个session的第一个scan对齐。然后再把每个session里的其它scan与该session的第一个scan进行对齐。使用此方式进行重排是因为各个session的数据之间可能会有较大差异。 Session选择session里所有的scan。 3.1.2 Estimation Options 这里包括各种注册参数选择项,若对某一个选项不确定,使用软件默认值即可。 Quality质量与速度的权衡。选择高质量以最慢的速度给出最精确的结果,低质量以较快的速度给出较不精确的结果。此参数的设定实际影响到的是参与参数估计的象元(voxel)的数目。其依据是有些象元(voxel)其实对重排参数的估计贡献不大,可以舍弃。 Separation 此参数以毫米为单位,表示对参考图像文件进行重采样时采样点之间的间隔。采样点之间间隔越小,结果越精确,运算速度越慢。 Smoothing (FWHM)高斯平滑的半高宽值。在估计重排参数之前一般先进行高斯平滑。PET数据一般使用7mm。MRI数据一般使用5mm。 Num Passes Register to first: 所有图像文件对齐注册到第一幅图像。 Register to mean: 使用two-pass处理将所有图像文件对齐注册到所有图像文件的平均图像。PET 数据一般注册到平均图像。因为PET数据相比fMRI数据噪音更大,文件更少,所以时间的影响更小。MRI数据一般注册到第一幅图像。虽然使用two-pass处理可能更精确,但是其对效果的提高与其所损失的运行时间相比得不偿失。 Interpolation在估计最佳变换时对数据进行重采样的方法。高的degree提供更好的结果,但是也更慢,因为会采样更多的相邻象元(voxel)[52, 53, 54]。 Wrapping此参数指示一个volume中数据wrap around in的方向(此处具体理解有待大家补充)。No wrapping:适用于PET数据或者已进行过空间变换的数据。同时当你不确定自己数据类型时,推荐使用此选项。Wrap in Y:适用于没有重排(resilce)过的在Y方向上进行相位编码的MRI数据。 Weighting提供一个加权图像,在估计重排参数时对参考图像的每一个象元进行加权。加权系数与标准差成反比。例如当有大量额外的头动(如说话或者特定区域内的严重伪影)时。(此处具体理解有待大家补充)。

6

3.2 Realign: Reslice(据已估计出的参数重排)

此功能重排以上步骤中已进行参数估计和注册的图像文件,使之与参考图像文件达到象元级的匹配精确。重排后的数据被命名为:r+原文件名。 3.2.1 Images选择要重排的数据文件

3.2.2 Reslice Options 各种重排参数设定,若对某一个选项不确定,使用软件默认值即可。Resliced images All Images (1..n) : 重排所有数据,包括标准参考图像(重排后还是保持原位置不变)。Images 2..n :重排除了标准参考图像之外的所有数据。此选项用于当你以MRI结构像为标准重排PET图像数据,而又不想在结果中再生成一个等同的MRI标准结构像时。All Images + Mean Image : 重排图像文件之外,另生成一个重排后的平均图像文件。Mean Image Only : 只生成重排后的平均图像文件。 Interpolation图像文件重采样和重写入的方式。Nearest Neighbour :最快,但不推荐使用。Bilinear Interpolation:可用于PET数据,但不是太适用于fMRI数据。Fourier Interpolation:此选项仅适用于纯刚体变换,也就是说象元大小必须是相同,并且等方性(正方体)的[17, 14]。 Wrapping 此参数指示一个volume中数据wrap around in的方向(此处具体理解有待大家补充)。No wrapping:适用于PET数据或者已进行过空间变换的数据。同时当你不确定自己数据类型时,推荐使用此选项。Wrap in Y:适用于没有重排(resilce)过的在Y方向上进行相位编码的MRI数据。 Masking因为扫描过程中被试总会或多或少有头动,造成同一个时间系列数据里所采集到的图像的边界不会完全重合。在有些图像还有数据的地方(信号值大于0),其它一些图像已经超出了图像边界(信号值为0)了。在这些信号为0的区域是无法采样数据的,因此SPM只要检测到某一幅图像在某个区域已经超出了边界(即信号为0),就会将其它所有图像的此区域信号值均设为0。此做法相当于取了时间系列数据中所有图像的交集。

3.3 Realign: Estimate & Reslice

将上述参数估计与数据重排合到一起做。全部选项与参数原理均与3.1和3.2中对应项相同。

4、Normailze

选用realign步骤中得到的平均象与模板进行比较,获得进行标准化的参数,参数文件命名为filename_sn3d.mat,然后依据此参数文件对每个img文件进行标准化,生成文件nfilename.img。具体操作如下:在预处理面板标准化选项中选择“Normalise: Estimate & Write”,出现如下对话框:

7

我们做如下设置:选中“data”——“new sbject”,在data下新出现的 “subject” 选项中作如下设置, “source image”选择空间校准步骤中生成的mean文件, “image to write”选择所有刚进行完校准的文件“ra*.img”,“template image”我们选择“EPI.nii”,其余采用默认设置,点绿三角运行。

5、Smooth

FWHM推荐为象素大小的两至三倍。在预处理面板标准化选项中选择“smooth”,出现如下对话框:

8