Nilearn教程系列(3)-ICA静息功能磁共振成像的分组分析:CanICA

简介

将CanICA应用于静息状态数据的示例。此示例将其应用于ADHD200数据集的30个主题。
然后，它绘制一个地图，所有的组件一起和轴向切割的每个组件分开。
CanICA是一种用于fMRI数据组级分析的ICA方法。与其他策略相比，它带来了一个良好控制的组模型，以及一个阈值算法控制特异性和敏感性的显式模型的信号。

案例

1.加载ADHD200数据


from nilearn import datasetsadhd_dataset = datasets.fetch_adhd(n_subjects=30)
func_filenames = adhd_dataset.func  # list of 4D nifti files for each subject# 数据集基本信息
print('First functional nifti image (4D) is at: %s' %func_filenames[0])  # 4D data

2.将CanICA应用在核磁共振图像


from nilearn.decomposition import CanICAcanica = CanICA(n_components=20, smoothing_fwhm=6.,memory="nilearn_cache", memory_level=2,threshold=3., verbose=10, random_state=0)
canica.fit(func_filenames)# 检索大脑空间中的独立成分。可通过属性components_img_直接访问。
components_img = canica.components_img_
#components_img是一个Nifti Image对象，可以使用以下行保存到文件中：
components_img.to_filename('canica_resting_state.nii.gz')

3.可视化

3.1将所有的组件都放在同一组图片上

from nilearn.plotting import plot_prob_atlas
import warnings
warnings.filterwarnings("ignore")
# Plot all ICA components together
plot_prob_atlas(components_img, title='All ICA components')

3.2将所有的ICA组件单独绘制

from nilearn.image import iter_img
from nilearn.plotting import plot_stat_map, show
import warnings
warnings.filterwarnings("ignore")
for i, cur_img in enumerate(iter_img(components_img)):plot_stat_map(cur_img, display_mode="z", title="IC %d" % i,cut_coords=1, colorbar=False)show()

文献参考:

G. Varoquaux et al. “ICA-based sparse features recovery from fMRI datasets”, IEEE ISBI 2010, p. 1177
G. Varoquaux et al. “A group model for stable multi-subject ICA on fMRI datasets”, NeuroImage Vol 51 (2010), p. 288-299

脑机学习者Rose笔记分享,QQ交流群：903290195
更多分享，请关注公众号

Nilearn教程系列(3)-ICA静息功能磁共振成像的分组分析:CanICA相关推荐

Nilearn教程系列(4)-fMRI解码入门教程(一)
这是使用Nilearn解码的教程,它以Haxby 2001研究中猫辨别任务的数据为基础. fMRI解码入门教程功能性磁共振成像(FMRI,functional magnetic resonance ...
小甲鱼 OllyDbg 教程系列 (三) ：PJ 软件功能限制
小甲鱼OllyDbg教程:https://www.bilibili.com/video/av6889190?p=8 https://www.freebuf.com/articles/system/87 ...
Nilearn教程系列(4)-脑部地图集绘制
目录脑部地图集绘制第一步:下载数据集第二步:可视化Harvard-Oxford 脑部地图集玻璃脑图(Glass brain)绘制本分享为脑机学习者Rose整理发表于公众号:脑机接口社区 .Q ...
Nilearn教程系列(2)-3D和4D niimgs:处理和可视化
目录 3D和4D niimgs:处理和可视化第一步:加载数据第二步:可视化可视化4D文件本分享为脑机学习者Rose整理发表于公众号:脑机接口社区 .QQ交流群:903290195 3D和4D ...
小甲鱼 OllyDbg 教程系列 (六) ：PJ 软件功能限制（不修改jnz的非爆破方法）
小甲鱼视频地址:https://www.bilibili.com/video/av6889190?p=13 小甲鱼OD学习第9讲:https://www.bbsmax.com/A/x9J23xEM56 ...
cocos creator麻将教程系列(八)—— 达达麻将语音聊天源码分析
达达麻将语音聊天源码分析达达麻将版图语音聊天 1:语音聊天只支持Native平台,iOS与android; 2: 语音聊天的音频格式为amr; 3: native平台实现了语音的录制和播放,可以移 ...
第四章脑功能磁共振成像
第四章脑功能磁共振成像 fMRI简介 HBP:人类脑计划,利用现代化信息工具探索大脑现代神经科学:多学科.多水平的综合研究领域脑功能成像fMRI及其特点脑功能成像fMRI: fMRI根据MRI ...
静息态功能磁共振成像:关于静息态功能连接和脑网络分析方法
自诞生之初,人类就对大脑中发生的事情充满好奇.功能磁共振成像是一种重要的工具,它有助于无创地检查.定位和探索大脑的语言.记忆等功能.近年来,神经科学研究的焦点明显转向了"静息态"下 ...
eeglab教程系列(14)-Working with ICA Components
在本案例前,先要进行进行ICA数据分解操作: Tools > Run ICA,具体操作参照eeglab教程系列(12)-使用ICA分解数据的第三步. 1.绘制component spectra ...