Fieldtrip 和 spm 文件读取

原文地址：Fieldtrip 和 spm 文件读取作者：yyigeren

FIELDTRIP官方网 http://fieldtrip.fcdonders.nl/start
SPM  官方网 http://www.fil.ion.ucl.ac.uk/spm/

做脑影像的人应该都知道Fieldtrip 和 SPM这两个Matlab toolbox。用他们分析数据，事半功倍，如有神助。
目前我主要关注EEG信号的处理。在SPM上下载了EEG的example 数据，把流程走了一遍，感觉这两个工具箱各有所长：SPM的图形化界面非常，Fieldtrip的函数开放性强。

目前用到Fieldtrip集成的ICA，感觉非常好用。有时候需要用fieldtrip，有时候需要用spm，下面就把数据读取、格式相互转换的函数总结一下。

一、关于数据读取
1、Fieldtrip可以直接读Biosemi仪器采集的数据，还包含其他很多种格式的数据（见http://fieldtrip.fcdonders.nl/dataformat）
读取程序如下：data = ft_read_data(filename); % 读取的是eeg数据。还可以加参数控制，读取head之类
或者整个读取，得到1*1的struct
filename_raweeg = 'subject1.bdf';
cfg = [];
cfg.dataset = filename_raweeg;
cfg.bpfilter = [0.5 40]; % Band pass filter
cfg.detrend = 'yes'; % Detrend
[data] = ft_preprocessing (cfg);

2、SPM数据读取
1）SPM读取Biosemi采集的数据，需要先convert
SPM的图形界面太傻瓜了，直接在matlab中输入spm eeg，手动convert *.bdf数据，就保存成*.mat数据了。
SPM有个很好的地方：图形化界面的操作可以save as script，变为matlab code。 不过程序没有直接用fieldtrip函数写的来的整洁。
2）SPM读取mat数据
D = spm_eeg_load('spm8_subject1.mat');

二、关于数据转换
data = spm2fieldtrip(D);
D = spm_eeg_ft2spm(data, 'filedata.mat'); % data为fieldtrip数据，filedata.mat为存放data的matlab文件。