磁共振实验数据SPM8处理流程
MRI数据预处理流程资料讲解
数据处理基本流程由于MRI是断层扫描,耗费时间较长,患者在进行MRI扫描的时候不可避免的会头部挪动,导致照射出来的图像不能一一映射;不同人的头颅,脑部大小,形状都会有所差异,获得的MRI图像也千差万别,无法对其进行对比。
所以我们就必须用一种算法将所有的MRI图像进行空间转换到一个比较标准的空间(目前使用较多的是被神经学家广泛认可的Talairach坐标系)将各个解剖结构一一对应后,再与标准化图谱或者不同个体之间相互比较(目前使用的是Talairach-Tournoux图谱)本文使用的是SPM软件和MRIcro软件处理图像数据,将MRI图像进行数据分析。
数据分析的基本流程:(1)数据预处理:○1图像格式转换○2slice timing获取时间校正○3realign头动校正○4Coregister不同成像方法间的图像融合○5nomalize 不同被试之间的图像标准化(归一化)○6smooth空间平滑《2 3 4统称图像的空间变换》(2)模型构建与参数估计:○:1建立统计模型○2将数据应用于统计模型○3进行参数统计得到单个被试的结果,多个被试的组分析数据预处理SPM是一款以MATLAB为平台的软件,所以使用SPM前一定要安装MATLAB。
打开MATLAB软件,界面如下:1.图像格式转换。
在进行数据预处理第一步要先将图像格式转换成SPM可以识别的ANALYZE格式。
转换之前先将原始数据放在MATLAB下面的mri image文件夹下,将路径设置成D:\MATLAB\work\mri image\ 设置过程如下:点击红色方块所指的按钮,在弹出的窗口中选择工作路径,按确定按钮即可。
设置完工作路径后,利用如下方法,将SPM2及其所有子文件夹添加到MATLAB的搜索途径中(1.点击file按钮,在下拉菜单选择set path2.在弹出的路径设置窗口点击"Add Folder"浏览并选择目标文件夹,eg:D:\spm2\3.点击save按钮4.点击close按钮,完成添加)在打开SPM之前,应先确定默认变量的设置是否准确,具体做法如下:1.在matlab命令窗口输入“edit spm_defaults"打开spm_defaults.m文件2.查看defaults.analyze.flip条目,确认defaults.analyze.fip值是否为1,若不是,改成1打开SPM:在matlab命令窗口输入“spm"回车后出现下面窗口,按黄色长方形覆盖的按钮,方可打开SPM软件(或者直接输入spm fmri即可打开)左上角的窗口代表按钮窗口,用以对数据进行处理分析;左下方的窗口为输入窗口,右边的窗口为图形窗口,用以显示图像结果图像格式转换:1.点选按钮窗口中Tool boxes菜单中的Dicom选项2.在弹出对话框中选择所有*IMA文件,点done 3.转换完毕,删除原文件夹下的*IMA数据 4.结果将在该文件夹下生成三个*.img/hdr/mat 文件5.结束图像格式转换方法二:用MRIcro软件进行格式转换(1)单击convert foreign toanalyze(2)出现第二个窗口,将左边的数据填完后,单击箭头所指“选择”选择所要转换的文件夹后,单击确定即可生成两个文件"*i.img/hdr" 。
核磁数据处理方法
核磁数据处理方法一、引言核磁共振(Nuclear Magnetic Resonance,NMR)是一种重要的分析技术,广泛应用于化学、生物学、医学等领域。
核磁数据处理方法是对从核磁共振仪器中获得的原始数据进行处理和分析的过程。
本文将介绍一种常用的核磁数据处理方法,包括数据预处理、谱线拟合和数据解释等步骤。
二、数据预处理1. 数据导入将从核磁共振仪器中获得的原始数据导入到数据处理软件中。
通常,原始数据以多个文件的形式存在,每一个文件对应一个谱图或者一个谱段。
2. 数据校正对导入的原始数据进行校正,以消除仪器本身的误差和噪声。
校正方法包括零点校正、线性校正和非线性校正等。
3. 数据滤波对校正后的数据进行滤波处理,以去除高频噪声和干扰信号。
滤波方法可以采用低通滤波、高通滤波或者带通滤波等。
4. 数据平滑为了减小数据的噪声和波动性,可以对滤波后的数据进行平滑处理。
平滑方法包括挪移平均、Savitzky-Golay平滑和傅里叶平滑等。
三、谱线拟合1. 谱线拟合模型的选择根据实际问题和数据特点,选择合适的谱线拟合模型。
常见的谱线拟合模型包括高斯模型、洛伦兹模型、Voigt模型等。
2. 初始参数的设定根据经验或者先验知识,设定谱线拟合模型的初始参数。
初始参数的设定对拟合结果的准确性和稳定性有重要影响。
3. 谱线拟合算法的选择选择适合的谱线拟合算法进行计算。
常用的谱线拟合算法包括最小二乘法、非线性最小二乘法、Levenberg-Marquardt算法等。
4. 拟合效果评估评估谱线拟合的效果,包括拟合优度(Goodness of Fit)指标、残差分析、参数可信区间等。
通过评估拟合效果,可以判断拟合模型的合理性和可靠性。
四、数据解释1. 谱峰分析对谱线拟合的结果进行分析,提取谱峰的位置、强度和宽度等信息。
谱峰分析可以匡助确定样品的成份和结构。
2. 化学位移计算根据谱线拟合结果和参考物质,计算样品中各组分的化学位移。
核磁数据处理方法
核磁数据处理方法核磁共振(Nuclear Magnetic Resonance, NMR)是一种重要的分析技术,广泛应用于化学、生物学、医学等领域。
核磁数据处理方法是对核磁共振实验所得到的数据进行处理和分析的过程,旨在从原始数据中提取有用的信息,并进行解释和解析。
一、数据采集与预处理1. 核磁共振实验数据的采集:核磁共振实验通常通过激励样品中的核自旋,利用核自旋与外加磁场的相互作用来获取信号。
实验过程中需要记录样品的信号强度和频率等参数。
2. 数据预处理:核磁共振实验所得到的原始数据通常包含噪声和干扰。
在进行进一步分析之前,需要对数据进行预处理,包括去除噪声、基线校正、频率校准等。
常用的预处理方法有傅里叶变换、高斯滤波、峰值对准等。
二、谱图处理与解析1. 谱图处理:核磁共振实验所得到的数据通常以谱图的形式呈现。
谱图处理是对谱图进行分析和解读的过程。
常见的处理方法包括峰识别、峰积分、峰拟合等。
2. 化学位移和耦合常数的测定:核磁共振谱图中的峰代表了不同化学物质的核自旋共振信号。
通过测定峰的化学位移和耦合常数,可以确定分子的结构和化学环境。
3. 数据解析与结构推断:通过对谱图中峰的位置、形状和强度等进行分析,可以推断出分子的结构和化学性质。
常见的解析方法包括化学位移对应关系、峰形分析、耦合常数分析等。
三、数据模拟与拟合1. 数据模拟:通过对已知结构的分子进行核磁共振数据的模拟,可以得到理论上的谱图。
数据模拟可以帮助解释实验数据中的峰和信号,并验证分析结果的准确性。
2. 数据拟合:对实验数据进行拟合,可以得到与实验数据最佳拟合的理论谱图。
拟合过程中需要调整分子结构参数和谱图参数,以寻找最佳拟合结果。
四、数据可视化与报告1. 数据可视化:对处理和分析后的数据进行可视化呈现,可以更直观地展示分析结果。
常见的可视化方法包括谱图绘制、峰图绘制、化学结构图绘制等。
2. 报告撰写:将数据处理和分析的结果整理成报告,包括实验方法、数据处理过程、分析结果和结论等。
SPM8数据预处理
SPM8数据预处理流程一、数据准备处理前首先通过SPM DICOM Import工具将dicom数据转换成NIfTI格式;转完后删除前10个volumes文件二、数据处理流程数据处理包括slice timing,realignment,normalization和smoothing四步注意:如果图像获取是隔层(interleaved)进行的,如1、3、5、7、9、2、4、6、8、10,则要先进性slice timing 再进行realign,如果图像各层是连续(sequential)获取的,如1、2、3、4、5、6、7、8、9、10,则要先进行realign 再做slice timing。
1:Slice TimingData:双击Data,双击Session,添加数据Number of Slices:输入每祯图像的层数,如“32”TR:输入TR时间,一般为2秒,我们输入“2”TA:是每祯图像获取第一层开始到获取最后一层图像的时间间隔,我们输入TR-TR/nslice,可直接输入公式,如输入“2-2/32”Slice order:我们输入“1:2:31, 2:2:32”。
指定层获取顺序的层次序参数是一个含N个数的向量,这里N是每个volume所含的层数。
每一个数表示该层在图像(volume)中的位置。
向量内的数字排列顺序是这些层的获取时间顺序。
如行向量 [1 3 5 7 9 11 13 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])各种扫描类型和输入的层顺序如下:ascending 升序扫描(first slice=bottom): [1:1:nslices];descending 降序扫描(first slice=top): [nslices:-1:1];interleaved 间隔扫描(middle-top):for k = 1:nslices,round((nslices-k)/2 + (rem((nslices-k),2) * (nslices - 1)/2)) + 1,end; interleaved (bottom -> up): [1:2:nslices 2:2:nslices],如[1:2:25,2:2,25];interleaved (top -> down): [nslices:-2:1, nslices-1:-2:1]Reference Slice:我们输入“31”。
磁共振实验数据SPM8处理流程
磁共振实验数据SPM8处理流程磁共振实验数据的处理流程主要包括原始数据预处理、功能图像处理和统计分析三个步骤。
其中,SPM8是一款常用的功能磁共振成像软件,可以用于实现磁共振实验数据的处理与分析。
以下是磁共振实验数据SPM8处理流程的详细步骤:1.原始数据预处理:(1)数据导入:将原始数据导入SPM8软件,通常数据的格式为DICOM 或NIfTI格式。
(2)切片时间校正:根据数据采集过程中的切片顺序和时间信息,对数据进行切片时间校正,以获得各个时间点上的切片数据。
(3)运动校正:校正数据中由于受试者运动引起的图像运动,通常使用刚体变换或流形形变方法进行运动校正,以消除不同时间点之间的运动偏移。
(4)标准空间转换:将运动校正后的数据转换到参考空间中,通常使用统一分割方法,将每个受试者的数据映射到一个标准的空间模板上。
(5)降噪处理:对数据进行降噪处理,一般使用高斯平滑方法,平滑数据以减小噪声干扰。
2.功能图像处理:(1)激活检测:通过对原始数据进行激活检测,识别可能激活区域。
常用的方法包括广义线性模型(GLM)和独立成分分析(ICA)等。
(2) 功能图像标准化:将激活检测得到的功能图像进行标准化处理,以便进行后续的统计分析。
常用的方法包括z-score标准化和百分位标准化等。
(4)时间序列提取:从功能图像中提取感兴趣脑区的时间序列数据,以便进行统计分析。
3.统计分析:(1)组间比较:使用统计学方法对不同组别的数据进行比较,常用的方法包括单样本t检验、双样本t检验和方差分析等,可以得到激活区域和激活程度的差异信息。
(2) 相关性分析:分析不同脑区之间的功能连接关系,通常使用Pearson相关系数或互信息等方法进行分析。
(3)多变量分析:使用多变量模式识别方法(如支持向量机、随机森林等)对脑图像数据进行分类或预测,以探索脑区之间的潜在模式。
以上是磁共振实验数据SPM8处理流程的主要步骤,这些步骤将原始数据进行预处理,提取功能图像,以及进行统计分析,从而得出与研究目的相关的结论。
核磁数据处理方法
核磁数据处理方法核磁共振(NMR)是一种重要的分析技术,广泛应用于化学、生物学和医学等领域。
在进行核磁实验时,我们需要对得到的数据进行处理和分析,以获得实用的信息。
本文将详细介绍核磁数据处理的方法和步骤。
一、数据获取和预处理1. 数据获取:通过核磁共振仪器获取样品的核磁共振谱图。
谱图通常包含两个维度的数据,即横轴表示化学位移,纵轴表示信号强度。
2. 数据预处理:对获得的谱图进行预处理,包括去除噪声、基线校正和相位校正等。
去除噪声可以提高信噪比,基线校正可以消除谱图中的基线漂移,相位校正可以调整信号的相位。
二、数据处理和分析1. 峰识别:通过峰识别算法找出谱图中的峰,并确定其化学位移和峰面积。
常用的峰识别算法包括峰拟合法和峰积分法。
2. 化学位移校正:根据已知的内部标准品或者参考物质的峰位置,对谱图进行化学位移校正,以提高数据的准确性和可比性。
3. 峰归属:通过与数据库或者文献比对,确定峰的归属,即确定峰所代表的化合物或者功能基团。
4. 峰积分和积分校正:对峰进行积分,计算峰的面积。
如果需要,还可以进行积分校正,以消除不同峰的积分因子差异带来的影响。
5. 峰定量:根据峰的面积和浓度关系,进行峰定量分析,计算样品中目标物质的含量。
6. 数据可视化:将处理后的数据进行可视化展示,通常使用谱图或者曲线图来展示数据的特征和变化趋势。
三、常用的数据处理软件和工具1. NMRPipe:是一个常用的核磁数据处理软件,提供了丰富的数据处理和分析功能。
2. MestReNova:是一款功能强大的核磁数据处理软件,可以进行多维数据处理、峰识别和峰归属等操作。
3. TopSpin:是一款广泛应用于核磁实验室的数据处理软件,提供了直观的用户界面和丰富的数据处理功能。
4. MATLAB:是一种常用的科学计算和数据处理工具,可以编写自定义的数据处理算法和脚本。
四、注意事项和常见问题1. 数据质量:核磁数据处理的结果受到数据质量的影响,因此在进行数据处理前,需要确保获得的数据质量良好。
核磁数据处理方法
核磁数据处理方法核磁共振(Nuclear Magnetic Resonance,NMR)是一种重要的分析技术,广泛应用于化学、生物、医学等领域。
核磁共振实验产生的数据需要进行处理和分析,以提取有用的信息。
在本文中,将介绍一种常用的核磁数据处理方法,包括数据预处理、峰识别、峰积分和谱图解析等步骤。
1. 数据预处理核磁共振实验得到的原始数据通常包含噪声和基线漂移等干扰。
为了提高信噪比和准确性,需要对数据进行预处理。
常见的预处理方法包括零填充、傅里叶变换、基线校正和峰对齐等。
1.1 零填充零填充是指在原始数据的两端添加零值,以增加数据点的数量。
这样可以提高频谱分辨率,使峰更加清晰。
零填充后的数据可以通过快速傅里叶变换(Fast Fourier Transform,FFT)进行频谱分析。
1.2 傅里叶变换傅里叶变换是将时域数据转换为频域数据的数学方法。
通过傅里叶变换,可以将核磁共振信号从时间域转换为频率域,得到谱图。
傅里叶变换可以通过离散傅里叶变换(Discrete Fourier Transform,DFT)或快速傅里叶变换来实现。
1.3 基线校正基线校正是指对数据中的基线漂移进行校正,以消除基线对峰的影响。
常用的基线校正方法包括多项式拟合、小波变换和自适应基线校正等。
1.4 峰对齐峰对齐是指将不同谱图中的峰位置对齐,以便进行比较和分析。
峰对齐可以通过寻找共有峰或使用内部参考物质来实现。
2. 峰识别峰识别是指从处理后的数据中找出峰的位置和强度。
常用的峰识别方法包括阈值法、波形拟合法和小波变换等。
2.1 阈值法阈值法是一种简单直观的峰识别方法,通过设置一个阈值来判断哪些数据点属于峰。
超过阈值的数据点被认为是峰的一部分。
2.2 波形拟合法波形拟合法是一种更精确的峰识别方法,通过拟合峰的形状来确定峰的位置和强度。
常用的拟合函数包括高斯函数、洛伦兹函数和Voigt函数等。
2.3 小波变换小波变换是一种时频分析方法,可以同时提取峰的位置和强度信息。
大物实验~~核磁共振实验数据处理
大物实验~~核磁共振实验数据处理核磁共振(NMR)实验是物理学和化学领域中常用的一种实验方法,其数据处理过程包括多个步骤,包括数据采集、数据预处理、数据分析和数据可视化等。
以下是对这些步骤的详细描述。
一、数据采集在核磁共振实验中,数据采集是实验的核心部分。
实验人员需要设置适当的实验条件,如磁场强度、射频脉冲频率和脉冲宽度等,以获取清晰的核磁共振信号。
在实验过程中,通常使用核磁共振谱仪来收集数据。
核磁共振谱仪可以产生高精度的射频脉冲,并测量它们与原子核之间的相互作用。
二、数据预处理数据预处理是去除噪声和干扰,提高数据质量的过程。
在核磁共振实验中,数据预处理包括对数据进行平滑处理、基线校正、相位调整等操作。
这些操作可以改善数据的信噪比,并使后续的数据分析和可视化更加准确。
三、数据分析数据分析是通过对预处理后的数据进行处理和分析,提取有关样品中原子核分布的信息。
在核磁共振实验中,数据分析包括对谱峰的识别、峰面积的测量、化学位移的计算等操作。
这些操作可以得出原子核在不同磁场下的分布情况,从而了解样品的分子结构和化学性质。
四、数据可视化数据可视化是将数据分析得到的结果以图表的形式呈现出来。
在核磁共振实验中,数据可视化包括绘制核磁共振谱图、制作三维图像等操作。
这些图像可以直观地展示样品中原子核的分布情况,帮助实验人员更好地理解实验结果。
除了以上四个步骤外,核磁共振实验的数据处理还包括其他一些步骤,如实验设计、实验操作、数据处理和结果解释等。
这些步骤需要实验人员具备一定的物理学和化学知识,以及对数据处理方法的了解和应用能力。
在核磁共振实验中,数据处理是一个非常重要的环节。
通过对数据的采集、预处理、分析和可视化,实验人员可以得出有关样品中原子核分布的信息,并了解样品的分子结构和化学性质。
这些信息对于科学研究、化学分析、材料开发等领域都具有重要的意义。
需要注意的是,核磁共振实验的数据处理过程具有一定的复杂性和专业性,需要实验人员具备一定的技能和经验。
核磁数据处理方法
核磁数据处理方法概述:核磁共振(NMR)是一种常用的分析技术,在化学、生物学和医学等领域具有广泛应用。
核磁数据处理是指对采集到的核磁共振数据进行处理和分析的过程。
本文将介绍一种标准的核磁数据处理方法,包括数据预处理、峰识别、峰积分、峰归属和数据可视化等步骤。
1. 数据预处理:在进行核磁数据处理之前,需要对采集到的原始数据进行预处理,以去除噪声和不必要的信号。
常见的数据预处理方法包括基线校正、峰对齐和峰平滑等。
基线校正是通过拟合基线曲线并将其从原始数据中减去来消除背景噪声。
峰对齐是将不同样品的峰位置进行校正,以便后续的峰识别和峰积分。
峰平滑是通过应用平滑算法来减少数据中的噪声。
2. 峰识别:峰识别是核磁数据处理的关键步骤,其目的是自动识别出数据中的峰,并确定其峰位置和峰强度。
常见的峰识别算法包括阈值法、拟合法和小波变换法等。
阈值法是根据峰的信号强度与背景噪声的差异来确定峰的位置。
拟合法是通过拟合峰形曲线来确定峰的位置和峰强度。
小波变换法是利用小波变换对数据进行频谱分析,从而确定峰的位置和峰强度。
3. 峰积分:峰积分是核磁数据处理的重要步骤,其目的是计算每个峰的面积,以获取样品中各组分的相对含量。
常见的峰积分方法包括直接积分法、拟合积分法和峰高法等。
直接积分法是将峰的面积直接计算为峰下的积分值。
拟合积分法是通过拟合峰形曲线来计算峰的面积。
峰高法是将峰的面积计算为峰高与峰宽的乘积。
4. 峰归属:峰归属是核磁数据处理的关键步骤,其目的是将每个峰与相应的化合物进行关联,以确定样品中各组分的化学结构和含量。
常见的峰归属方法包括峰数据库比对、化学位移计算和二维谱图分析等。
峰数据库比对是将实验测得的峰与已知化合物的峰数据库进行比对,以确定峰的归属。
化学位移计算是通过计算峰的化学位移与已知化合物的化学位移进行比对,以确定峰的归属。
二维谱图分析是通过绘制二维谱图并进行交叉比对,以确定峰的归属。
5. 数据可视化:数据可视化是核磁数据处理的最后一步,其目的是将处理后的数据以图形的形式进行展示,以便于结果的分析和解释。
核磁数据处理步骤
核磁数据处理步骤介绍核磁共振(Nuclear Magnetic Resonance, NMR)是一种重要的分析技术,可用于分析和研究物质的结构和性质。
核磁共振数据处理是核磁共振实验中非常重要的一步,它涉及到数据的预处理、数据解析和数据解释等多个环节。
本文将详细介绍核磁数据处理步骤,并探讨每个步骤的具体内容和作用。
核磁数据处理步骤核磁数据处理通常包括以下几个步骤:1. 数据获取核磁共振实验需要通过核磁共振仪获取原始数据。
在这一步骤中,需要设置实验参数,如磁场强度、扫描方式和扫描范围等,以获取合适的核磁共振谱图。
同时,还需要进行系统校准,以保证数据的准确性和可靠性。
2. 数据预处理数据预处理是核磁数据处理的关键一步,它主要包括去噪、基线校正和谱图平滑处理。
去噪是指去除谱图中的噪声信号,常用的方法有傅里叶变换滤波和小波变换等。
基线校正是指对谱图中的基线进行修正,常用的方法有多项式拟合和简单直线法等。
谱图平滑处理是为了提高数据的信噪比,常用的方法有移动平均法和高斯平滑法等。
3. 数据解析数据解析是将核磁共振谱图中的峰进行定量分析,以确定样品中各组分的相对含量和结构。
这一步骤主要包括峰识别、峰集成和峰归属等。
峰识别是指在谱图中找出所有的峰,并对其进行编号和标记。
峰集成是指对每个峰进行积分,以得到峰面积和峰高等定量信息。
峰归属是指将每个峰与相应的化学位移和耦合常数进行关联,以确定相应的化学结构和相互作用类型。
4. 数据解释数据解释是核磁共振谱图中各峰的化学解释,以确定各个峰的来源和物质的结构。
这一步骤主要包括化学位移解释、耦合常数解释和化学结构解释等。
化学位移解释是指将峰的化学位移与特定化学官能团和化学键联系起来,以确定它们的存在和相对含量。
耦合常数解释是指通过峰之间的耦合常数和相对强度,推断出化学键的取向和相互作用类型。
化学结构解释是将所有的峰归属进行整合,以得到最终的化学结构和分子式。
核磁数据处理步骤的应用核磁数据处理步骤在化学、药学、生物学和材料科学等领域具有广泛的应用。
功能磁共振数据处理一般流程
功能磁共振数据处理⼀般流程⼀般基于MATLAB平台使⽤SPM⼯具包进⾏处理。
由于SPM操作较为复杂,不适合批处理,因⽽有很多实验室开发了⼀系列基于SPM的⼯具包,也即开发界⾯,调⽤SPM功能实现操作计算。
具体的⼯具包会再另⼀篇⾥详述。
本部分主要进⾏流程简述。
⼀、预处理0.删除Slice:为了防⽌初期设备不稳定,删除最初的⼏张slice(4-10);不过现在机器都有预热时间,开始试验后都已经进⼊稳定⼯作状态,这步也可不做。
1.时间校正 (slice timing)由于磁共振图像采⽤层层扫描,层与层之间保留⼀定的间隔,但是这也很难避免相邻两个层⾯之间的影响。
为保证采集到数据的准确性及数据的空间分辨率,常采⽤间隔扫描的⽅法,即⾸先采集 1,3,5……层,然后对2,4,6……层进⾏采集。
但⽆论是隔层扫描还是连续扫描,任意两层的采集时间是不相同的,通常采集⼀个体(volume)需要⼀个TR 的时间。
尤其是间隔扫描,这使得连续两层之间相差TR/2的时间。
在进⾏静态磁共振数据处理分析的时候,这样⼀个时间段会产⽣很⼤的影响,需要对每个体的时间进⾏校正,也就是将每个体的所有时间有同⼀个时间点。
||注意:在设置的时候,要检查⼀共扫了多少层,如果是奇数层(如33)顺序应为1:2:33,2:2:32,中间层为33;若是偶数层(如28),则为 2:2:28,1:2:27,中间层为28.2.头动校正(head motion correction)spm 软件对数据的头动处理主要是将第⼀幅图像(?存疑:是否应该是中间层?)作为基准图像,然后通过旋转或平移等刚性变换将时间序列上的其它图像与第⼀幅图像的位置匹配,然后⽤内插值算法对这系列的图像进⾏重新采样。
此步会⽣成头动的位置变化和⾓度变化图,以及头动⽂件,包含六列参数。
3.图像配准(coregister)将⼀个⼈的图像配准到⼀起,如功能像配到本⼈的结构像上。
要求将图像进⾏配准,使得格式⼤⼩统⼀。
高原脱习服期大脑结构变化
3. SPM8的处理流程及算法原理介绍本文处理数据用的是SPM8这个软件,SPM也就是我们所知道的统计参数图, 也是本文所使用的主要软件,经过这个软件处理的磁共振成像的数据最终得出的就是统计参数图,SPM这个软件是由Friston 教授等人在Matlab的基础上做出的一款数据处理以及统计功能非常强大的软件,他们设计这个软件最初是为了用统计的方法对SPECT,PET 和fMRI的数据进行处理和分析,最后得到统计参数图。
在使用SPM 对数据进行处理和分析的时候主要分为两大部分:一、对数据的预处理过程;二、对数据的统计分析过程。
那么下面就会比较详细的将SPM8的数学原理及操作步骤呈现出来。
3.1 SPM8的处理流程及算法原理1. 头动校正(realignment)因为脑成像所持续的时间比较长,而且需要测量的次数比较的多,由于被试者的呼吸和血流脉动等本身内部原因,所以头动是无法避免的,因此进行头动校正。
头动校正指对每一幅图像以一副图像(一般的为第一副图像)为标准进行的变换从而来抵偿未知的差异。
简要说来,对象图像向目标图像的配准是我们需要处理的问题。
头动校正只针对同一被试者的图像数据,通过刚性变换来描述两张不同的图像之间所存在的位移关系。
刚性变换分为平移和旋转这两种坐标变换过程,表征变换的未知量是使用“自动图像匹配法”求出来的。
这种算法认为假如有两副图像能完全重合,那么这两幅图像相对应的像素灰度值的比将会在任意的位置上都会呈现出接近于一个常量,所以它是按照最小化的两幅图灰度值比方差来求取的位移参数。
采用刚性变换时,头部运动可以分解成头部旋转和平移两个部分[5]。
如果是在三维空间内,则可以用6 个参数对这种变换进行表示。
原始的图像和目标图像分别用f 和g 来表示,图像的变换第利用矩阵相乘来计算,以便于计算机实现。
设在变换之前的空间坐标是,经过变换之后的空间坐标是,则可以得出:(2)其中是坐标原点平移之后的变换,,,分别是绕x,y 和z 轴旋转了角之后的变换。
核磁数据处理方法
核磁数据处理方法一、引言核磁共振成像(MRI)是一种非侵入性的医学成像技术,广泛应用于临床诊断和科学研究领域。
核磁共振成像通过对人体组织中的核自旋进行激发和检测,产生高质量的图象。
然而,MRI数据的处理是实现准确诊断和分析的关键步骤之一。
本文旨在介绍一种标准的核磁数据处理方法,以提供详细的步骤和指导。
二、数据预处理1. 数据获取:首先,需要获取原始的核磁共振成像数据。
这些数据通常以DICOM格式保存,并包含多个图象切片。
2. 数据校准:对于每一个图象切片,需要进行校准以消除可能存在的伪影和运动伪影。
这可以通过运动校准算法和伪影校正算法实现。
3. 数据重采样:如果需要将数据与其他成像模态进行配准或者比较,可以对数据进行重采样,以使其具有相同的空间分辨率和尺寸。
三、图象增强1. 噪声去除:核磁共振成像数据通常受到噪声的影响,因此需要进行噪声去除处理。
可以使用滤波算法,如高斯滤波或者中值滤波,来减少噪声。
2. 对照度增强:为了更好地显示组织结构和病变区域,可以使用直方图均衡化或者对照度拉伸等方法来增强图象的对照度。
3. 空间滤波:如果图象存在含糊或者伪影,可以使用空间滤波算法,如锐化滤波或者边缘增强滤波,来提高图象的清晰度和细节。
四、图象分割1. 区域生长:区域生长算法是一种常用的图象分割方法,可根据像素之间的相似性将图象划分为不同的区域。
可以根据特定的阈值或者生长准则来执行区域生长。
2. 边缘检测:边缘检测算法可用于检测图象中的边缘和轮廓。
常用的边缘检测算法包括Sobel算子、Canny算子等。
3. 阈值分割:阈值分割是一种简单而有效的分割方法,可根据像素的灰度值将图象分为不同的区域。
可以使用全局阈值或者自适应阈值来执行阈值分割。
五、特征提取1. 形态学操作:形态学操作可用于提取图象中的结构特征,如边缘、孔洞和连通区域等。
常用的形态学操作包括膨胀、腐蚀、开运算和闭运算等。
2. 纹理分析:纹理分析可以提取图象中的纹理特征,用于描述组织的细微变化。
核磁数据处理方法
核磁数据处理方法核磁共振(Nuclear Magnetic Resonance, NMR)是一种重要的分析技术,在化学、生物化学、药物研发等领域得到广泛应用。
核磁数据处理方法是指对核磁共振实验所得到的数据进行处理和分析的一系列步骤和算法。
一、数据预处理1. 数据校正:对核磁共振仪器所得到的原始数据进行校正,包括信号强度校正和基线校正。
信号强度校正可以通过内部参考物质的信号进行,基线校正可以通过拟合多项式或者使用特定算法进行。
2. 数据滤波:对原始数据进行滤波处理,以去除噪声和干扰信号。
常用的滤波方法包括高斯滤波、中值滤波等。
二、数据处理1. 峰识别:通过寻觅数据中的峰,确定样品中的化学成份及其相对浓度。
常用的峰识别方法包括阈值法、导数法等。
2. 峰积分:对识别出的峰进行积分,计算峰的面积,从而得到样品中各化学成份的绝对浓度。
常用的峰积分方法包括峰高积分、峰面积积分等。
3. 谱图拟合:对核磁共振谱进行拟合,以确定样品中的化学结构和相对浓度。
常用的谱图拟合方法包括高斯拟合、洛伦兹拟合等。
三、数据分析1. 化学位移分析:通过分析核磁共振谱中的化学位移,确定样品中的化学环境和化学结构。
化学位移分析可以通过与数据库对照或者使用特定算法进行。
2. 耦合常数分析:通过分析核磁共振谱中的耦合常数,确定样品中的化学键和份子结构。
耦合常数分析可以通过与已知标准对照或者使用特定算法进行。
3. 动力学分析:通过分析核磁共振谱中的峰形和峰位随时间的变化,研究样品中的反应动力学和份子运动。
四、数据可视化1. 谱图绘制:将处理后的核磁共振谱数据进行绘图,以直观展示样品中的化学成份和结构。
2. 谱图解释:对绘制的核磁共振谱进行解释,标注峰的化学位移和相对强度,解析样品中的化学结构和组成。
以上是核磁数据处理方法的普通步骤和算法,具体的处理方法和参数设置可以根据实际需求和样品特性进行调整。
核磁数据处理方法的准确性和可靠性对于得到准确的化学信息和结构解析非常重要,因此在处理过程中应严格按照标准格式进行操作,并进行必要的验证和校正。
如何进行核磁原始数据处理
目前,脑成像数据主要有DTI、fmri、3D三种模态。
这些数据在分析前都要进行格式转换,不同公司的扫描仪存储格式也不尽相同。
脑成像处理软件也很多,不同软件使用的格式也不一样,所以数据转换是脑成像数据处理的第一步,必须非常清楚。
这里主要以siemens的机器为准,介绍在windowx下的MRIcron 的dcm2nii转换和MRIConvert转换.从扫描中心下载的原始数据是以dicom数据格式存在的压缩文件,解压后,得到原始文件。
来自siemens的扫描仪的原始文件以“IMA”下为后缀。
对于功能像(fMRI)的数据,有多少个TR就有多少个IMA图像文件,即每个IMA 文件就是一个完整的volume;对于DTI数据,有n个方向,有m个b0像,就有n+m张IMA图片,即n+m个完整的volume。
当然有的DTI数据有的只有一个b0像,有的有6个b0像之多。
对于3D结构像数据,如果扫描了128层,就会有128张IMA图像,每张图像就是一张slice,不是volume。
数据转换后,主要有spm2之前使用的Analyze格式,以及fsl和spm5和spm8使用的NifTI_1格式。
Analyze格式是成对的hdr和img文件表示一个3D的volume,而NifTI_1格式可以是3D也可以是4D的,同时可以是hdr和img 成对文件,也可以是NifTI_1的nii一个文件。
如下:Spm2使用3D Analyze hdr/img;spm5和spm8使用3D NifTI hdr/img.fsl使用NifTI_1的4D的nii 格式。
目前数据转换主要有MRIcron的dcm2nii转换和MRIConvert转换。
现在一一介绍一下:在MRIcron的安装目录下,有一个dcm2nii.exe和dcm2niigui.exe,并且分别有:dcm2nii.nii和dcm2niigui.nii两个配置文件。
dcm2nii.exe是Dos的命令行操作,而dcm2niigui.exe是图形界面。
核磁数据处理方法
核磁数据处理方法一、引言核磁共振(Nuclear Magnetic Resonance,简称NMR)是一种重要的分析技术,广泛应用于化学、生物医学、材料科学等领域。
核磁共振技术通过观察核自旋在外加磁场下的行为,获取样品的结构和性质信息。
核磁共振实验产生的数据需要进行处理和分析,以提取实用的信息。
本文将介绍一种标准的核磁数据处理方法,以匡助研究人员准确、高效地处理核磁数据。
二、数据预处理1. 数据导入首先,将核磁共振实验得到的数据导入数据处理软件中。
常见的数据格式包括Bruker格式、Varian格式等。
根据实际情况选择合适的导入方式,并确保数据导入的准确性。
2. 数据校正在数据导入后,通常需要进行数据校正,以消除仪器噪声和仪器漂移对数据的影响。
校正方法包括零点校正、基线校正等。
校正后的数据更加准确可靠,有利于后续处理和分析。
三、数据处理1. 数据平滑为了降低数据的噪声,往往需要对数据进行平滑处理。
平滑方法包括挪移平均、高斯平滑等。
选择合适的平滑方法,并根据实际情况调整平滑参数,以获得平滑后的数据。
2. 数据相位校正在核磁共振实验中,样品的信号可能受到相位偏移的影响,导致数据失真。
因此,需要进行相位校正以恢复数据的准确性。
相位校正方法包括零次相位校正、一次相位校正等。
选择合适的相位校正方法,并进行必要的调整,以获得准确的数据。
3. 数据峰识别核磁共振实验中,往往需要识别出样品中的峰,以确定样品的组成和结构。
数据峰识别方法包括阈值法、峰拟合法等。
根据实际情况选择合适的峰识别方法,并进行必要的参数调整,以获得准确的峰信息。
4. 数据积分核磁共振实验中,峰的面积与样品中的物质浓度相关。
因此,需要对数据进行积分,以计算峰的面积。
积分方法包括峰高积分、峰面积积分等。
选择合适的积分方法,并进行必要的参数调整,以获得准确的积分结果。
四、数据分析1. 数据拟合核磁共振实验中,往往需要对数据进行拟合,以获得样品的参数信息。
SPM中文教程汇总已修订版
S P M中文教程汇总已 Document number:PBGCG-0857-BTDO-0089-PTT1998一、SPM的安装与启动先安装matlab,然后将SPM复制到matlab下的一个文件夹(SPM2需要或以上版本)。
启动matlab,首先set path,然后在matlab命令窗口中输入SPM即可启动,然后选择fMRI,也可以直接输入SPM fMRI二、SPM数据处理概要先将所得数据进行空间预处理(对齐,平滑,标准化等),然后进行模型估计(将刺激的时间、间隔与血流动力函数进行卷积,所得结果与全脑象素信号进行相关分析),最后察看结果。
三、SPM8数据处理的一般步骤为方便后续的数据处理,如果数据分散处理后整合,建议所有处理数据路径保持一致,要统一路径。
处理前首先要采用数据转换软件将dicom数据转换成SPM解析格式,然后进行数据预处理,预处理结束后到matlab安装目录中备份spm*.ps文件,其中包含了空间校正和标准化的信息,然后进行建模分析。
运行命令:spm fmri,打开spm8的操作界面我们称左上侧的窗口为按钮窗口(button window),左下侧的窗口为输入窗口(input window),右侧大窗口为树形结构窗口或图形窗口(Tree Building Window or the graphics window)。
在spm8和spm5中,每一步处理都采用了直观的“树形结构”的面板,如果一个分支项左面有“+”号,你可以双击显示子分支项,如果一个分支项右面有“<-X”号,你必须为之指定选项(否则不能运行该tree),分支项的选项在其右侧面板指定,而帮助信息则在下面的面板中显示。
如果我们处理数据没有特殊需求,我们只关心带有“<-X”项目并完成输入即可,其余均可采用默认设置。
另外注意在Tree Building Window的顶部菜单,新增了一个菜单项“TASKS”,在使用批处理分析时非常重要。
功能磁共振数据处理
Outline
知识准备:
FMRI原理
自旋 进动 磁化矢量 共振 豫弛
T1像 T2像
自旋回波
TR & TE T1像&T2像
BOLD 实验设计
Block
Event-related
SPM:
SPM介绍 数据准备 预处理
结果呈现
实验设计
Outline
知识准备:
FMRI原理
自旋 进动 磁化矢量 共振 豫弛
T1像 T2像
自旋回波
TR & TE T1像&T2像
BOLD 实验设计
Block
Event-related
SPM:
SPM介绍 数据准备 预处理
分析时需要
Normalize
Why?
T2像和模板配准
How? 备注:
模板T1.nii在 spm8\templates 目录下
Voxel sizes: 333
← mean*.img ← r*.img ← T1.nii
Normalize
批处理
Smooth
Why?
平滑配准后的图像
脱氧血红蛋白比氧合血红蛋白更具有顺磁性 血流量越大氧合血红蛋白的比例越高
与基线相比顺磁的脱氧血红蛋白浓度越小 横向弛豫速度减慢,信号衰减更慢 活动脑区信号相对较强
Outline
知识准备:
FMRI原理
自旋 进动 磁化矢量 共振 豫弛
T1像 T2像
自旋回波
TR & TE T1像&T2像
SE序列结构
MRI数据预处理流程资料讲解
Command Window ? xTo fet ft art ed, select HATLAB HCporDgmowThe element type "namr" must be terminated byCouilri not parse the ■file: d.: \matlab7\t oolboK\c »数据处理基本流程由于MRI是断层扫描,耗费时间较长,患者在进行MRI扫描的时候不可避免的会头部挪动,导致照射出来的图像不能一一映射;不同人的头颅,脑部大小, 形状都会有所差异,获得的MRI图像也千差万别,无法对其进行对比。
所以我们就必须用一种算法将所有的MRI图像进行空间转换到一个比较标准的空间(目前使用较多的是被神经学家广泛认可的Talairach坐标系)将各个解剖结构一一对应后,再与标准化图谱或者不同个体之间相互比较(目前使用的是Talairach-To urnoux 图谱)本文使用的是SPM软件和MRIcro软件处理图像数据,将MRI图像进行数据分析。
数据分析的基本流程:(1)数据预处理:。
图像格式转换C2 slice timing获取时间校正◎ realign头动校正C Coregister不同成像方法间的图像融合C5 nomalize不同被试之间的图像标准化(归一化)C smooth空间平滑《2 3 4统称图像的空间变换》(2)模型构建与参数估计:C建立统计模型C2将数据应用于统计模型C3进行参数统计得到单个被试的结果,多个被试的组分析数据预处理SPM是一款以MATLAB为平台的软件,所以使用SPM前一定要安装MATLAB。
File Edit View Graphics Debug Desktop Windlow Help□ R 晶■宦EC 體当I D:lmatlab?\work 二二]阖Shortcuts E How to Add 0 What's NewStvt1•图像格式转换在进行数据预处理第一步要先将图像格式转换成SPM 可以识别的ANALYZE 格式。
磁共振实验数据SPM8处理流程
磁共振实验数据SPM8处理说明一、数据准备为方便后续的数据处理,如果数据分散处理后整合,建议所有处理数据路径保持一致,要统一路径。
处理前首先要采用数据转换软件将dicom数据转换成SPM 解析格式,转换时格式请选择NIfTI,可用SPM输入面板中的DiCOM Import模块转换,也可以采用专门的转换软件,如MRIcovert。
然后进行数据预处理,预处理结束后到matlab安装目录中备份spm*.ps文件,其中包含了空间校正和标准化的信息,然后进行建模分析。
二、数据处理流程1:Slice TimingSpm8和以往版本一样,先安装matlab7.1以上版本,然后将spm8文件夹放在matlab安装目录中,我们一般放在toolbox文件夹中,然后启动matlab用菜单中设置路径命令增加spm8的路径并保存。
之后我们运行命令:spm fmri,这样将打开spm8的操作界面,我们称左上侧的窗口为按钮窗口(button window),左下侧的窗口为输入窗口(input window),右侧大窗口为树形结构窗口或图形窗口(Tree Building Window or the graphics window)。
Slice Timimg用来校正1个volume中层与层之间获取(采集)时间的差异,对事件相关设计的实验尤为重要。
我们在按钮窗口中的预处理面板中点击“Slice Timimg”,将出现如下对话框:在spm8和spm5中,每一步处理都采用了直观的“树形结构”的面板,如果一个分支项左面有“+”号,你可以双击显示子分支项,如果一个分支项右面有“<—X”号,你必须为之指定选项(否则不能运行该tree),分支项的选项在其右侧面板指定,而帮助信息则在下面的面板中显示。
如果我们处理数据没有特殊需求,我们只关心带有“<—X”项目并完成输入即可,其余均可采用默认设置。
另外注意在Tree Building Window的顶部菜单,新增了一个菜单项“TASKS”,在使用批处理分析时非常重要。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
磁共振实验数据SPM8处理说明一、数据准备为方便后续的数据处理,如果数据分散处理后整合,建议所有处理数据路径保持一致,要统一路径。
处理前首先要采用数据转换软件将dicom数据转换成SPM 解析格式,转换时格式请选择NIfTI,可用SPM输入面板中的DiCOM Import模块转换,也可以采用专门的转换软件,如MRIcovert。
然后进行数据预处理,预处理结束后到matlab安装目录中备份spm*.ps文件,其中包含了空间校正和标准化的信息,然后进行建模分析。
二、数据处理流程1:Slice TimingSpm8和以往版本一样,先安装matlab7.1以上版本,然后将spm8文件夹放在matlab安装目录中,我们一般放在toolbox文件夹中,然后启动matlab用菜单中设置路径命令增加spm8的路径并保存。
之后我们运行命令:spm fmri,这样将打开spm8的操作界面,我们称左上侧的窗口为按钮窗口(button window),左下侧的窗口为输入窗口(input window),右侧大窗口为树形结构窗口或图形窗口(Tree Building Window or the graphics window)。
Slice Timimg用来校正1个volume中层与层之间获取(采集)时间的差异,对事件相关设计的实验尤为重要。
我们在按钮窗口中的预处理面板中点击“Slice Timimg”,将出现如下对话框:在spm8和spm5中,每一步处理都采用了直观的“树形结构”的面板,如果一个分支项左面有“+”号,你可以双击显示子分支项,如果一个分支项右面有“<—X”号,你必须为之指定选项(否则不能运行该tree),分支项的选项在其右侧面板指定,而帮助信息则在下面的面板中显示。
如果我们处理数据没有特殊需求,我们只关心带有“<—X”项目并完成输入即可,其余均可采用默认设置。
另外注意在Tree Building Window的顶部菜单,新增了一个菜单项“TASKS”,在使用批处理分析时非常重要。
对上图右侧选项我们做如下设置:Data,点击data并在下面的面板中点击“new session”,这样在data下会出现“session”的分支项,选中该项并点击面板下方的“select files”,然后用spm文件选择器选择你要处理的数据,最后点击“down”。
选择数据时可以把静息态、数值任务和物理大小任务分为三个session来选(data——new session——session),也可以作为一个session来选,结果是一样的。
Number of Slices,我们输入每祯图像的层数,如“32”TR,我们输入重复时间,一般为2秒,我们输入“2”TA,是每祯图像获取第一层开始到获取最后一层图像的时间间隔,我们输入TR-TR/nslice,可直接输入公式,如我们输入“2-2/32”Slice order,我们输入“1:2:31,2:2:32”。
指定层获取顺序的层次序参数是一个含N个数的向量,这里N是每个volume所含的层数。
每一个数表示该层在图像(volume)中的位置。
向量内的数字排列顺序是这些层的获取时间顺序。
如行向量[135791113151719212325246810121416182022 24](在Matlab中可表示为[1:2:25,2:2,25])各种扫描类型和输入的层顺序如下:ascending升序扫描(first slice=bottom):[1:1:nslices];descending降序扫描(first slice=top):[nslices:-1:1];interleaved间隔扫描(middle-top):for k=1:nslices,round((nslices-k)/2+(rem((nslices-k),2) *(nslices-1)/2))+1,end;interleaved(bottom->up):[1:2:nslices 2:2:nslices],如[1:2:25,2:2,25];interleaved(top->down): [nslices:-2:1,nslices-1:-2:1]Reference Slice,我们输入“31”。
选择参考层,通常选择n slice/2,如25层时选择13层作为参考层。
Filename Prefix,是指新生成的图像前加何标记,一般采用默认设置。
最后点击面板上方的向右的绿色三角即开始运行。
运行完后将会生成一系列a*.img文件,这就是时间校准后的数据。
注意很多研究者容易将时间校准和空间校准顺序颠倒,一般的观点是如果图像获取是隔层(interleaved)进行的,如1、3、5、7、9、2、4、6、8、10,则要先进性slice timing再进行realign,如果图像各层是连续(sequential)获取的,如1、2、3、4、5、6、7、8、9、10,则要先进行realign再做slice timing。
2:Realignment即使我们对被试的头部做了很好的固定,在实验过程中,被试也会不由自主的有一些轻微的头动,这在fMRI实验中尤为明显。
这一步就是把一个实验序列中的每一帧图像都和这个序列的第一桢图像按照一定的算法做对齐,以矫正头动。
做完这一步,能给出该序列中被试的头动情况,以作为是否放弃该数据的依据,如果头动超过1个voxel(功能图像扫描矩阵一般是64*64,则体素的大小为(FOV/64)*(FOV/64)*(层厚+层间距)),则要考虑放弃该时间点数据。
该程序利用最小二乘法(least squares approach)原理和含6个参数(刚体模型)的空间变换,对从一个被试获取的时间序列进行校正。
用户可指定某个volume 作为随后volumes的参考。
可以是第1个volume,也可选择比较有代表性的volume(更明智的选择),例如选择磁场相对稳定的第4个volume。
校正信息(头动信息)将在结果窗口(Graphics Window)显示。
每个Session的校正信息将存储为rp*.txt,其中*为Session数据集名称。
另外,头动校正信息将以plot图形显示。
如下图:Translation表示被试头部在X,Y,Z三个方向的平移,分别用蓝,绿,红三种颜色表示。
Rotation表示被试头部在实验过程当中绕X(L-R),Y(A-P),Z(S-I)三条轴的转动角度。
横坐标代表这个序列所采集的所有图像,纵坐标表示的是偏移量和偏转角度,分别以毫米和度为单位。
采用SPM8,头动信息和空间标准化的图形文件将以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)。
其余选项采用默认设置,点击上方绿色的三角开始运行。
3:Normalise在预处理面板标准化选项中选择“Normalise:Estimate&Write”,出现如下对话框:我们做如下设置:选中“data”——“new sbject”,在data下新出现的“subject”选项中作如下设置,“source image”选择空间校准步骤中生成的mean文件,“image to write”选择所有刚进行完校准的文件“ra*.img”,“template image”我们选择“EPI.nii”,其余采用默认设置,点绿三角运行。
4:Smoothing在预处理面板标准化选项中选择“smooth”,出现如下对话框:我们在“image to smooth”选项中选择所有刚进行完标准化的文件“wra*.img”然后点绿三角运行即可。
这里FWHM我们采用默认设置“888”。
5:Specify1st level即以前版本的“fMRI model”,spm5和spm8种的分析选项有所变化,使用Specify1st level做单个被试(single subject)分析;使用Specify2nd level 做组分析(group analysis)。
我们选择“Specify1st level”,出现如下对话框:我们选择Directory指定一个文件夹存放结果数据,其余做如下设置:“Units for design”选择“Scans”,“Interscan interval”输入“2”,选择“Data and Design”后选择“New Subject/Session”,再选择新出现的“Subject/Session”,选择“Scans”并用文件选择器选择相应任务的所有平滑后的功能图像(swra*.img)然后点击“down”,选择“Condition”后选择“New condition”,然后选中新出现的“Condition”,“name”选项输入任务条件的名称,“onset”输入任务条件的启动向量,代表任务刺激启动的扫描数,选中“Durations”输入任务组块的持续时长,如是事件相关设计请输入“0”。
如还有其他任务,要再次选中“Condition”——“New condition”——“Condition”定义其他任务条件。
设置完毕后点击绿三角运行。
这样将会在开始选择的目录中生成文件spm.mat.下面要估计我们刚建立的模型,在模型设置面板中点击“estimate”,将打开如下对话框:很简单,我们只须选择刚生成的“spm.mat”文件点击“down”然后点击绿三角运行即可。
估计完成后,我们选择“results”,将打开如下对话框:上图设计矩阵表明一个扫描序列中我们有三个任务条件,均为事件相关设计,选中“t—contrasts”,点击“define new contrast”,第一个任务“contrast”定义为“1”,第二个任务定义为“01”,第三个任务nn定义为“001”,第一个任务减第三个任务定义为“10-1”,第二个任务减第三个任务定义为“0 1-1”,第二个任务减第一个任务定义为“-11”,其余操作以及激活图显示和以前版本都是一致的。
值得我们注意的是,我们指定或输入一系列的参数来进行每一步的处理,完成后我们可以通过Save按钮将每一步存为一个*.mat文件。
以后我们可以通过Load按钮重新加载并使用这些*.mat文件,我们适当修改后(例如改变所运行的数据集)再选择Run按钮运行。