fMRI数据处理（随时更新）

要开始学习处理fMRI的数据了。

一、使用matlab工具包SPM读取fMRI数据

（1）首先得安装工具包SPM，我参考的是下面这篇博客：

在matlab下安装spm工具_spmas包matlab-CSDN博客

（2）使用SPM读取数据，举例如下：

inputDir='D:\2_chenfang_labs\ECG-fMRI\result_removePhysio_40sub';
num=40; % number of subject
subSample=dir([inputDir,filesep]);

for i=1:num
    sub=i+2;
    imgNameSample = dir([inputDir,filesep,subSample(sub,1).name,filesep,'*.nii']);
    imgnum=length(imgNameSample);
    for jj=1:imgnum
        Dir_img=[inputDir,filesep,subSample(sub,1).name,filesep,imgNameSample(jj).name];
        vSub=spm_vol(Dir_img);  dim=vSub.dim; mat=vSub.mat;
        imgSub=spm_read_vols(vSub);
        data=reshape(imgSub,dim(1)*dim(2)*dim(3),1);
        data(isnan(data))=0;
        data1(jj,:,:,:)=data;% 读取fMRI数据
    end
end