fMRI处理:基于spm的dparsf

刚刚踏上研究生之旅,开启一段新的研究旅程,希望能记录下自己的脚印。
最近刚刚接触fMRI,对fMRI序列什么的还不是很了解,不过先上网上手学习了一点给予软件的预处理知识,先记录一下,下回能用上再回来看看。

dparsf安装

dparsf是基于matlab的spm,我是网上下载spm12和最新版本的dparsf的安装包,然后添加到matlab路径,就可以使用了。在matlab编辑区输入dpabi就可以出现界面了,如下

fMRI预处理

预处理步骤包括:可视化(Visualization)、去伪影(Artifact removal)、时间配准(Slice time correction)、头动校正(Motion correction)、生理噪音校正(Correction for physiological effect)、结构功能配准(Co-registration)、标准化(Normalization)和时空间滤波(Spatial and temporal filtering)。

1/2:把你的数据由dicom数据转成可以识别的nii形式,不过我的数据已经是nii的了,所以没有用到这个;
3:去除前面10个时间点,前面的数据可能存在一点扰动;
4:时间校正:在扫描一次完整大脑(Brain volume)的周期(TR,Repetition time)内,我们会扫好几片脑片(Slice),由于我们一个时间点只能扫描一张脑片,如果是按照顺序一张一张扫下去的话,每张脑片之间的扫面时间点都会有区别。比如最顶端的脑片的扫描时间会相对低端的脑片延迟2s;
5:时间矫正的时候我们选择的扫描顺序,如果数据是646432,那么扫描一般都是按照1:2:31,2:2:32的顺序扫描的;
6:参考层:一般就是选择5中排序中间的那一层作为参考层;
7:头动矫正:头动(Head motion)也是一个很严重的干扰。我们在做后续分析的时候需要假设每个体素时时刻刻对应的都是大脑的同一个位置。可是扫描的时候被试的头肯定会有微小的移动。这样的话一个体素可能在下一秒就对应了脑中另外一个截然不同的位置。因此我们会用刚性变换(Rigid body transformation)把所有的图像中的脑都固定在同一个靶位置(Target imagine)
8/9:这个地方就需要我们手动定义图像原点的位置了,一般选择作为配准模板的是MNI图像,我们手动定位的时候尽量把图像的原点定义在和MNI一样的位置,就是前联合的位置;
10:去头骨和配准:结构图像与功能图像配准(Co-registration)了。我们希望能把功能图像上的点定位在有着较高分辨率的结构图像上。一般来说我们会在实验的每轮扫描的开始扫一张结构图像。虽然说结构图像和功能图像是在同一论扫描里扫得的。但往往它们往往是没法重合覆盖在一起的。一是因为他们测量不是一个东西(编者注:一个是组织一个是血流信号),二是因为测量过程中会有头动;
11:分割和标准化,这一步选择了分割和DARTEL,就是说分割之后还会保存分割图和分割图到标准化空间的一个形变场。
12:消除一些头动的协方差变量;
13:标准化:dpabi里面标准化有三种方式。一种基于EPI,一种基于T1,一种基于DARTEL,可以看一下软件给出的不同标准化的原理;

14:滤波:对图像进行空间滤波时,我们经常使用高斯核(Gaussian kernel),该平滑核的尺寸(平滑范围)由半高宽(full width at half maximum, FWHM)决定;
15:平滑:在统计分析前对采集到的数据进行空间平滑(smoothing)是很普遍的,这可以提高信噪比(signal to noise ratio),消除伪影(artifacts),使数据更符合分布假设(distributional assumptions)。
16:你开始处理的文件夹,我的数据是nii的,就要把开始处理的文件夹名字定义为FunImg,FunImg子文件就是一个个subject,处理好的文件他也会自己生成对应名字的文件夹;一旦你改好这个名字,软件会自动识别里面的subject名字,不用自己添加;

以上就是我暂时用到dpabi功能的基本介绍了,
接下来我就要准备实操了,到时候再更微
有什么问题实操见。

fMRI处理:基于spm的dparsf相关推荐

  1. fmri数据分析图像格式及转换问题——基于spm讨论

    1.几大常用格式 这张ppt介绍了现在常用的三种格式:dicom主要由扫描器产生,一般是一个slice一个文件:analyze 格式= .img/.hdr spm下用的最多:nifti ,由fsl & ...

  2. fmri|SPM contrast manager

    整理一下遇到的基础问题 问题1:batch建立的contrast manager 和 results建的contrast是一样的么? 是的,一样的 建立contrast的两种方法,一种点点点,一种通过 ...

  3. fMRI基础理论知识学习

    前言 时隔多年,再次上线,重新经营csdn.自读研以来,不是干饭就是摆烂,实在颓废,能重新开始写博客,已然不易.在这里立下flag,争取以后每周都能写点什么东西,一来锻炼文笔,二来记录学习历程 我的研 ...

  4. spm,afni的安装和使用心得

    原文地址: http://hi.baidu.com/luo_xiang/blog/item/8866662990d5f0fb99250afb.html 一.Spm的安装与启动: 先安装matlab,然 ...

  5. NeoRS:新生儿静息态fMRI数据预处理流程

    大家好,这里是 "茗创科技" .茗创科技专注于脑科学数据处理,涵盖(EEG/ERP, fMRI,结构像,DTI,ASL, ,FNIRS)等,欢迎留言讨论及转发推荐,也欢迎了解茗创科 ...

  6. Science Advances:胶质细胞功能障碍导致抑郁症静息态fMRI功能连接异常

    1. 摘要 重度抑郁症(MDD)是一种毁灭性的精神障碍,影响着全球高达17%的人口.尽管通过静息状态功能磁共振成像(rsfMRI)检测的MDD患者存在全脑网络水平的异常,但这些网络变化的机制尚不清楚, ...

  7. 致远SPM解决方案之库存管理

    卷首语 基于SPM系统实现库存管理,帮助企业了解设材库存现状,灵活调配库存设材,提高效率,节约资源. 库存管理应用概述 库存管理包括设备档案.材料档案.物资到货单.物资入库单.物资领用单.物资调拨申请 ...

  8. 致远SPM解决方案之安全管理

    卷首语 基于SPM系统实现安全管理,帮助企业构建安全管理要素.梳理安全管理流程.设计安全管理文件.安全风险辨识及风险评估方法,降低企业安全风险. 安全管理应用概述 安全管理包含危险源辨识.安全培训.安 ...

  9. 致远SPM解决方案之招标管理

    卷首语 基于SPM系统实现招标管理,帮助企业解决招标可能面对的各种痛点难题,提高企业招标效率,降低企业招标管理成本. 招标管理应用概述 招标管理包括专家库.评审小组组建.评审方案审批.招标文件审批.资 ...

最新文章

  1. 环形熔断器设计与gobreaker源码分析
  2. Python 学习第十八天 js 正则及其它前端知识
  3. TI DSP入门学习
  4. js组装知识(待续……)
  5. unsigned char与char的区别
  6. r语言回归分析_R语言之cox回归分析
  7. Java设计模式学习总结(7)——结构型模式之适配器模式
  8. JavaScript语言精粹:对象——读书笔记
  9. PyQt5学习--基本窗口控件--QButton类控件
  10. 拓端tecdat|R语言用综合信息准则比较随机波动率(SV)模型对股票价格时间序列建模
  11. 小白都能了解的聚类算法之一(Kmeans与GMM)
  12. python信息安全书籍_GitHub - riusksk/secbook: 信息安全从业者书单推荐
  13. C#使用itextsharp生成PDF文件含页脚页眉
  14. 计算机房维修保养记录表,机房设备系统运行及维护记录学习表格.docx
  15. Shenyu网关本地打docker镜像包
  16. 断层约束的等值线追踪
  17. python控制小爱同学_神秘鸭,用Siri小爱同学语音助手控制你的电脑
  18. Codeforces Round #617 (Div. 3)
  19. 第 3-5 课:Flutter 调试及 Android 和 iOS 打包
  20. 《社会调查数据管理——基于Stata 14管理CGSS数据》一2.2 数据管理的工作标准...

热门文章

  1. 中国人身材测试小软件编码,测试 || 我怀疑你是假中国人,不信你来试试
  2. 微信小程序获取当前位置并根据经纬度跳转地图导航
  3. 做为一个中国的ITer,你感到耻辱吗?看CSDN的“软件中国2006风云榜之10大最具创新性技术”...
  4. 全国行政区域代码/编码
  5. Confidence-guided roadside individual tree extraction for ecological benefit estimation
  6. ESP8266给服务器发送数据
  7. sinx/x的极限为什么是1_在sinx/x中当x→0时极限为什么为1?
  8. 码云生成公钥 公钥的作用是什么?
  9. ZStack实践汇 | ZStack云平台应用堡垒机教程
  10. 联想拯救者Y7000如何重装系统