ADNI静息态功能核磁共振成像数据预处理总流程共如下八步

目录

1. 下载DICOM格式数据

2. DICOM -> NIFTI格式 输出s开头文件

3.  时间层校正Slice Timing 输出as开头文件

4. 头动校正Realignment 输出ras开头文件

5. 归一化Normalize 输出w开头文件

6. 平滑Smooth 输出sw开头文件

7. 配准Co-register reslice

8. 提取ROI


word资源:https://download.csdn.net/download/m0_37852937/19053920

1. 下载DICOM格式数据

登陆地址:http://adni.loni.usc.edu/data-samples/access-data/

输入登录账号及密码->ADNI GO->DOWNLOAD IMAGE Collections->Advanced Search(beta)->确定Research Group: EMCI LCMI AD CN Image Description:Resting State fMRI-> Display in result->csv doenload (注:网页加载很慢)

建议单个下,如果all同时下,会出现断网,网中断,下载失败或文件不全。且不好统计哪些是下载过的。

已知下载好的文件格式解压如下(注:文件必须是6720个dicom,如果少了,说明网络中断导致文件没下载全,需要重新下载。):

事先在每个subject里存好每个步骤对应路径(可以用程序在每个人对应路径下创建以下文件夹):

Matlab 配好spm12 环境,输入spm(以下三个窗口都不要随意叉掉,后面要用到)

Dicom import

文件输出选中6720个文件

 

文件输出保存路径

(如果点错了,想要重选,点击红框以选中路径,文件路径就会消失,重选即可)

选好后点击done:

想要批量处理,可以先直接把第一个没加文件的Module List复制

想要批量处理,可以先直接把第一个Module List 复制

显示处理进度,耗时5分钟

有时候源文件生出结果只有一个文件,这个subject是不可用的,

2. DICOM -> NIFTI格式 输出s开头文件

原始的多个dicm文件得到140的nii文件

3.  时间层校正Slice Timing 输出as开头文件

Number pf slice: 48

TR: 3

TA: 3-3/48

Slice order:输入修改参数 1:2:47 2:2:48

Reference Slice:输入47

选择上一步骤的140个文件,设置参数后运行。

4. 头动校正Realignment 输出ras开头文件

输入上一步骤得到的140个文件,没有参数设置,直接运行

运行后会得到如下头动文件,自行记录和判定即可。

得到的中间图查看是否符合规范

5. 归一化Normalize 输出w开头文件

修改Voxel sizes—[改为3 3 3]

6. 平滑Smooth 输出sw开头文件

将默认的[8,8,8],改成[6,6,6]

7. 配准Co-register reslice

选择Coregister(Reslice)

然后点击运行。运行完成后,将得到的rswaf开头的130个新文件。

8. 提取ROI

在matlab里输入dpabi,运行dpabi。在出现的可视化操作界面中,按下图顺序点击按钮

在出现的界面中,点击Add image 添加rswaf开头的187个文件

点击Define ROI 选择AAL atlas,然后点OK

然后点击Output Dir 来修改输出路径,接着点击Extract

运行完成后,在定义的输出路径下有130*7=910个文件

其中ROISignals开头的(130*2个文件)就是需要的fmri时间序列文件,

共包含txt和mat两种格式,两种格式文件所存内容一致,只是文件格式不一样。

ROISignals开头的每个mat文件,储存了1×116的矩阵,表示该时间点下,116个脑区各自的数值。将130个mat文件按顺序组合成130*116的矩阵(可用matlab批处理完成),即是最终的数据矩阵。

AAL模板中 前90个为大脑脑区,研究中使用较多。后26个为小脑脑区,研究中使用较少。

ADNI静息态功能核磁共振成像数据预处理总流程相关推荐

  1. 利用matlab的SPM12功能核磁共振成像数据处理_Auditory fMRI data(二)

    2.3 Segmentation图像分割 第一步:Highlight the "Volumes" field and then select the subject's regis ...

  2. 数据分析问题(异常值识别)中数据预处理部分流程(含2022年全国服务外包大赛实例)

      博主个人理解的数据预处理主要包括 个方面:读取文件 => 数据概览 => 缺失值填补 => 数据分布预览 => 衍生特征设计.这套流程在完成异常值识别时作为数据预处理时没有 ...

  3. 利用matlab的SPM12进行功能核磁共振成像数据处理_Auditory fMRI data(一)

    大脑成像的每次采样包括64个连续的切片.大脑的三维图像大小为64X64X64,每个像素的大小为3X3X3立方毫米.TR时间设置为7s.(TR是指第一個RF到下一個相同RF出現時,中間間隔的時間.RF, ...

  4. 利用matlab的SPM12进行功能核磁共振成像数据处理_Auditory fMRI data(四)

    懒癌患者可以直接利用SPM12 batch script的代码对数据处理进行操作,连数据都不需要下载,直接打开matlab,运行以下代码即可. data_path = fileparts(mfilen ...

  5. 功能核磁共振成像资料分享汇总

    一本书 Functional Magnetic Resonance Imaging Author(s): Scott A. Huettel; Allen W. Song; Gregory McCart ...

  6. 【教程】GF-1和GF-2PMS数据预处理全流程

    GF-1和GF-2的PMS数据全部的预处理流程主要包括5个部分:(1)辐射定标(2)大气校正(3)正射校正(4)图像融合(5)影像配准.由于全色图像无法进行大气校正,在定量遥感中一般仅使用多光谱数据. ...

  7. 数据预处理常用方法流程

    文章目录 数据读取和选择 缺失数据 找出缺失数据 处理缺失数据 丢弃观察值 丢弃特征 填充缺失数据 数值变量 分类变量 替换缺失值 异常值 如何找出异常值 数值变量 直方图 箱形图 分类变量 条形图 ...

  8. 从诱发反应中解码动态脑模式:应用于时间序列神经成像数据的多元模式分析教程

    多变量模式分析(MVPA)或大脑解码方法已经成为分析功能磁共振数据的标准做法.虽然解码方法已广泛应用于脑机接口,但其应用于时间序列神经成像数据(如脑磁图.脑电图)以解决认知神经科学中的实验问题是最近的 ...

  9. Principles of fMRI 1课程笔记8--fMRI的数据预处理

    刚采集的原始图像数据会经过一系列的预处理步骤.这些步骤主要是分辨并去除伪影(Artifact),或者是检验一些模型所需的假设是否成立.具体来说,有三个目的: 尽量减少因为数据采集(Data acqui ...

  10. 2022-2028全球及中国先进核磁共振成像可视化系统行业研究及十四五规划分析报告

    受新冠肺炎疫情等影响,2021年全球先进核磁共振成像可视化系统市场规模大约为 亿元(人民币),预计2028年将达到 亿元,2022-2028期间年复合增长率(CAGR)为 %.未来几年,本行业具有很大 ...

最新文章

  1. 功能测试常用6种方法_建筑管道常用的8种连接方法
  2. IOS开发学习记录第5天之C语言学习
  3. 常见的排序算法(1)
  4. 【java学习笔记-io流 文件读写和键盘读写】带缓存的输入/输出流和数据输入/输出流
  5. Go程序:变量声明、赋值与输出
  6. 【数据结构----笔记1】查找算法之【顺序查找和折半查找】
  7. windows API 菜鸟学习之路(一)
  8. 计算机技能比赛范文,计算机技能大赛策划书范文.doc
  9. cocos+kbe问题记录
  10. 我的世界java出生蘑菇岛,《我的世界》出生就有蘑菇岛和海底神庙的地图,附近还有村庄!...
  11. 使用3dmax渲染出二维画面的方法
  12. 如何卸载twincat3_TwinCAT软件安装FAQ
  13. 如何绕过校园网认证,不付费就可以使用自己的网络。
  14. 线性代数的本质(Essense Of Linear Algebra)[1]
  15. GNSS常用术语(缩略语)总结------------持续更新
  16. 疏忽了给该内页做外链
  17. 超详细的Android APP 集成 PayPal 境外支付 信用卡 借记卡 支付,支付宝 微信海外 境外支付
  18. DataGrip for Mac破解步骤详解,期限到2099年
  19. 解线性方程组的python实现(1)——高斯主元消去法
  20. 基础30讲 第六讲 中值定理

热门文章

  1. Photoshop脚本 图层的栅格化
  2. c++ gdal 矢量转栅格_gdal的矢量栅格化接口GDALRasterizeLayers使用(一)
  3. 外贸自建站优势渐明,你还要观望多久!2018加油
  4. Docker配置的Redis集群报错CLUSTERDOWN Hash slot not served的解决办法
  5. excel两个表格数据对比_表格技巧—Excel中如何查找相同的数据
  6. 因严重 OpenSSL 漏洞,Fedora 37 推迟至 11 月中旬发布
  7. 在Word2019中不能插入公式的解决办法
  8. 415错误及解决方法
  9. python获取图像灰度极值点_【图像处理】灰度图、亮度峰值极值查找
  10. 2021年中国程序员薪资和生活现状调查:年薪5-25万之间占比66.3%