学习笔记|fmri_3 preprocessing+batch

还是bilibili课程笔记 https://b23.tv/FvgWGP
P6 dicom转档与pre 这部分笔记是视频1h10min的地方开始(这集共两个半小时)
到这里已经整理好数据了

打开spm

两种方法
1、spm —— 出现窗口,选择 fmri
2、命令行输入: spm fmri

设置路径

step1: 在Menu界面,选择cd
在这里插入图片描述
左边点点点找路径
说两点注意的
1、找到想要的路径了,怎么选中:点击右面 . (对就是那个点,文件夹最上面,点击代表选中左边的路径),然后点击Done

在这里插入图片描述
2、点击之后如下图,. 不见了,最底下出现一行路径,最后点击Done 路径选择完成。
如果不想选这个路径,在最底下的那行路径上点一下,这行路径被删除。这时候你会发现右侧并没有刷新,这是一种保护机制,如果其他操作,需要选取文件时,这样可以避免重复选择。
如何刷新呢,点左边到上一级,再重新进下一级,页面就会刷新——该操作spm里通用。【视频1h20min左右】
在这里插入图片描述

一 slice timing

点击Menu——slice timing,出现下面右图内容

1 先读图

1、位置1显示在进行的操作是 slice timing
2、位置2 带X的都表示是需要补充输入内容的,只有把标X的都操作后,位置3的run图标才会亮
在这里插入图片描述

2 填内容

视频1h40min有详细说明,具体参数要看扫数据的核磁参数的设置,用不同地方核磁扫描时,核磁仪器参数可能有所不同。跟扫描员确定(一般会有仪器参数的pdf版说明)。
在这里插入图片描述

3 结果显示

跑完文件夹会多一批被试,文件夹的选择这里只选了一个run/一个block的扫描文件夹,没有全选,先跟着练一下,后面再批量。
选择是做第一轮(run1)的文件夹,原文件夹 437个文件,做完slice timing 文件夹有437*2 个文件,多的一倍新数据前面都对应有a做前缀。 如原数据 007_0345.img ——slice timing操作后——新数据:a007_0345.img
在这里插入图片描述

二 realign

选第三个,点击
在这里插入图片描述
出现realign的batch editor 页面,只有一个有X,就是选数据,这里就选上面跑出来的新数据,找带a前缀的
在这里插入图片描述
注:可以用filter筛选带a的, 输入 ^a.*,回车,全选——done——数据页面自动关闭——run
在这里插入图片描述

数据结果 新数据带前缀r

1、可以看Graphics看数据质量,里面两个折线图,translation和rotation,可以看到translation在0.1到-0.2范围内波动,算是可接受范围
在这里插入图片描述
2、realign后的数据 前缀变为ra
在这里插入图片描述

batch

再记视频最后一个内容,什么是batch,虽然听报告听了很多回,但还是在实践中学习,才理解的更深刻。
先说为什么要batch,像上面说的步骤,slice timing和realign,如果对拿来的数据进行这两步的处理,要一步操作完,再做另一个步骤。
如果有写好或输出好的batch文件,比如视频中的(下图),那么这个导入的batch文件是包含步骤比较全的了,从change directory 到results report(这个是连后面的统计步骤都做了的)。
这样的操作步骤,就不用一步步来,直接选好要处理的数据,一键run就跑完了。
后面再有新数据,也可以想同步骤,简直不要太方便。
在这里插入图片描述

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
clear;clc parentdir = 'F:\data process\fMRI\fmrioutput'; % 定义储存各被试源文件的上级文件夹 cd(parentdir); % 进入这个上级文件夹 allsubjects = dir('sub*');%查找该文件夹下的所有被试 subinfos = numel(allsubjects); for i=1:numel(allsubjects) % 对每个被试进行循环 cursubject = allsubjects(i).name; % 找到当前被试的名字 matlabbatch=cell(1); curWPAT = fullfile(parentdir,cursubject,'WPAT'); curfucout=fullfile('F:\data process\fMRI\fmrioutput',cursubject,'WPAT') matlabbatch{1}.spm.stats.fmri_spec.dir = {curfucout}; matlabbatch{1}.spm.stats.fmri_spec.sess.scans = cellstr(spm_select('ExtFPList', curWPAT, '^sw*.nii', Inf)) matlabbatch{1}.spm.stats.fmri_spec.timing.units = 'scans'; matlabbatch{1}.spm.stats.fmri_spec.timing.RT = 2; matlabbatch{1}.spm.stats.fmri_spec.timing.fmri_t = 16; matlabbatch{1}.spm.stats.fmri_spec.timing.fmri_t0 = 8; %% matlabbatch{1}.spm.stats.fmri_spec.sess.cond = struct('name', {}, 'onset', {}, 'duration', {}, 'tmod', {}, 'pmod', {}, 'orth', {}); matlabbatch{1}.spm.stats.fmri_spec.sess.multi = {'D:\data process\fMRI\onsets\subject(i)_run1.mat'}; matlabbatch{1}.spm.stats.fmri_spec.sess.regress = struct('name', {}, 'val', {}); matlabbatch{1}.spm.stats.fmri_spec.sess.tempxx=dir(fullfile(curfucout,'rp*.txt')) matlabbatch{1}.spm.stats.fmri_spec.sess.hpf = 128; matlabbatch{1}.spm.stats.fmri_spec.fact = struct('name', {}, 'levels', {}); matlabbatch{1}.spm.stats.fmri_spec.bases.hrf.derivs = [0 0]; matlabbatch{1}.spm.stats.fmri_spec.volt = 1; matlabbatch{1}.spm.stats.fmri_spec.global = 'None'; matlabbatch{1}.spm.stats.fmri_spec.mthresh = 0.8; matlabbatch{1}.spm.stats.fmri_spec.mask = {''}; matlabbatch{1}.spm.stats.fmri_spec.cvi = 'AR(1)'; matlabbatch{2}.spm.stats.fmri_est.spmmat(1) = cfg_dep('fMRI model specification: SPM.mat File', substruct('.','val', '{}',{1}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('.','spmmat')); matlabbatch{2}.spm.stats.fmri_est.write_residuals = 0; matlabbatch{2}.spm.stats.fmri_est.method.Classical = 1; matlabbatch{3}.spm.stats.con.spmmat(1) = cfg_dep('Model estimation: SPM.mat File', substruct('.','val', '{}',{2}, '.','val', '{}',{1}, '.','val', '{}',{1}), substruct('.','spmmat')); matlabbatch{3}.spm.stats.con.consess{1}.tcon.name = 'Old'; matlabbatch{3}.spm.stats.con.consess{1}.tcon.weights = 1; matlabbatch{3}.spm.stats.con.consess{1}.tcon.sessrep = 'none'; matlabbatch{3}.spm.stats.con.delete = 0; end;怎么改
最新发布
05-24

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值