静息态功能磁共振数据分析工具包使用手册

静息态功能磁共振数据分析工具包使用手册宋晓伟（[e-mail]）文档版本:1.3文档修订日期:2008-2-25北京师范大学认知神经科学与学习国家重点实验室目录一、开发背景介绍....................................................................................................1二、软件用途和技术特点.........................................................................................4三、设计与实现.......................................................................................................4四、测试..................................................................................................................5五、使用要求...........................................................................................................5六、使用方法演示....................................................................................................6(一)计算功能连接......................................................................................7(二)计算局部一致性..................................................................................9(三)计算低频振幅....................................................................................11七、详细使用说明..................................................................................................13(一)安装REST..........................................................................................13(二)卸载REST........................................................................................13(三)启动REST..........................................................................................13(四)在REST中设置待处理的数据目录.....................................................16(五)Mask的设定.....................................................................................16(六)在REST中设定输出参数...................................................................17(七)可选项：去线性漂移.........................................................................18(八)可选项：滤波....................................................................................19(九)局部一致性计算参数的设定..............................................................20(十)低频振幅计算参数的设定.................................................................21(十一)功能连接参数的设定.........................................................................21(十二)点击“Doall”开始计算.......................................................................23(十三)耗时估计...........................................................................................24(十四)其它工具...........................................................................................24八、附注说明.........................................................................................................26九、参考文献.........................................................................................................28静息态功能磁共振数据分析软件（REST）1一、开发背景介绍大脑是人体中最迷人也是人们了解最少的部分，科学家哲学家们一直在寻找大脑与行为、情感、记忆、思想、意识等的关系，却缺少一个非侵入性的高分辨率的技术方法来直接观察并确立这种联系，直到上世纪末功能磁共振成像（functionalmagneticresonanceimaging,fMRI）的出现（Ogawaetal.,1990），既能让人们观察到大脑结构又能让人们观察大脑结构的某一部分所具有的特定功能（Clare,1997）。fMRI机制是血氧水平依赖性（Bloodoxygenleveldependent,BOLD）信号的变化。目前认识到的大多数的脑功能都是通过对任务或刺激的控制，并同时记录与任务或刺激相应的行为学上的变化和神经活动的变化来得到的。从Hubel和Wiesel电生理学上的实验，到现在神经影像学上的认知激活实验范式，都说明这种方法是很成功的。如图1被试睁眼或闭眼交替进行，这种简单的任务刺激范式所带来的BOLD信号的变化可以清楚地在大脑的特定区域看到（图1是在视觉区），从而把大脑的功能和解剖结构联系了起来（Foxetal.,2007）。这种基于任务刺激的实验范式一般都使用广义线性模型（Generallinearmodel,GLM）计算刺激或控制变量的效应，检测相应于刺激的大脑激活区，从而认识大脑的功能。图1、传统fMRI任务激活范式的分析：睁眼闭眼任务范式和初级视觉皮层的某个体素的BOLD信号。（引自Foxetal.,2007）对任务状态fMRI数据的分析和处理，研究者现在一般都使用软件SPM（Friston,1995）或AFNI（Cox,1996），这两个软件都可以使研究者很方便地基于GLM模型来分析和处理任务状态的fMRI数据。如图2是包括2个控制变量的GLM模型，研究者需要提供给软件的是设计矩阵，即研究者的控制变量，然后使用软件SPM（Friston,1995）或AFNI（Cox,1996）就可以很方便地估计出控制变量的效应大小，进而找到受控制变量影响的脑区，即和任务刺激相对应而激活的脑区。静息态功能磁共振数据分析软件（REST）2图2、包括2个控制变量的GLM模型（引自）但是GLM模型局限于对控制变量效应的估计，对已知刺激之外的自发BOLD信号的波动无法解释，按噪声来处理。而Foxetal.(2006)发现自发的BOLD信号可以显著地解释一部分在任务之上的BOLD信号的波动，并认为任务引起的和自发的BOLD信号的波动在人脑上是线性叠加的。Biswaletal.（1995）发现在没有明显的运动任务状态下，大脑左侧的感觉运动区和右侧感觉运动区、运动区内侧的BOLD信号波动却非常相关，说明静息状态下这种持续的、自发的BOLD信号波动并不是随机噪声，而是具有一定的模式。使用fMRI任务激活范式在很多任务条件下都发现，后扣带回、楔前叶、下顶叶、前额叶内侧等脑区总是表现出任务非特异性的活动减低（Raichleetal.,2001）。相比于任务状态下，在被动注视或静息状态下这些脑区都表现得更加活跃。Greiciusetal.（2003）研究发现静息态下后扣带回、前扣带回腹侧、下顶叶、前额叶内侧等脑区的BOLD信号高度相关，说明这些脑区在静息状态下的自发活动有很强的一致性，有很强的功能连接，构成一个功能网络。这些静息态下自发的脑神经活动或许与巩固记忆、预测未来、保持警觉等神经功能有关（Raichleetal.,2007）。静息态fMRI已经被用于某些神经精神疾病的研究，比如儿童注意缺陷多动障碍（Attention-DeficitHyperactivityDisorder，ADHD）、阿尔兹海默氏病（Alzheimer’sDisease，AD）等，Zangetal.（2007）发现静息态下ADHD的右侧额下回、左侧感觉运动区、双侧小脑、右侧扣带回前部、双侧脑干活动异常。Heetal.（2007）发现在静息态AD的后扣带回的自发活动与正常对照组相比明显降低。综上所述，研究静息态时的自发神经活动对理解神经功能有重要意义。另外，静息态fMRI一般连续数分钟扫描，没有任何特定的任务，不需要复杂精细的实验设计，也不需要被试做出反应，只要求被试安静地躺在扫描仪里，故研究者容易操作和控制扫描任务，被试也容易配合，易于多中心、大样本的研究，这是静息态fMRI脑成像的优点，非常适合临床上的研究和应用。目前静息态数据的计算方法中线性相关的方法被使用得最多，另外还有局部一致性（RegionalHomogeneity,ReHo）、低频振幅（AmplitudeofLowFrequencyFluctuation,ALFF）、独立成分分析（independentcomponentanalysis,ICA）等方法。REST目前的版本中首先加入静息态功能磁共振数据分析软件（REST）3了线性相关功能连接方法、ReHo方法以及ALFF方法。线性相关是静息态fMRI研究中应用最广泛的方法，它考察的是感兴趣区域（Regionofinterest,ROI）之间或者ROI与全脑所有的体素之间的线性相关程度，由此判断是否与ROI在功能上有较高的相似性，即有无功能连接。这种方法简单、灵敏、易于解释（Foxetal.,2007）。“局部一致性”(Zangetal.,2004)假设在一定条