forked from netstim/leaddbs
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathea_createCM_fmri.m
22 lines (18 loc) · 968 Bytes
/
ea_createCM_fmri.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
function [fMRI_CM,gmtc] = ea_createCM_fmri(options)
expfolder=[options.root,options.patientname,filesep,'connectomics',filesep,options.lc.general.parcellation,filesep];
% check if rest_preprocessing has been performed:
[~,restfname]=fileparts(options.prefs.rest);
[~,anatfname]=fileparts(options.prefs.prenii_unnormalized);
if ~ea_coreglocked(options,['r',restfname,'_',anatfname]) || ...
~exist([options.root,options.patientname,filesep,'r',restfname,'_',anatfname,'.nii'],'file') % preprocessing needs to be performed
disp('No (approved) preprocessed fMRI-images found, processing...');
ea_preprocess_fmri(options);
disp('Done preprocessing fMRI data.');
end
if ~ea_coreglocked(options,['r',restfname,'_',anatfname]) || ~exist([expfolder,options.prefs.rest,'_tc.mat'],'file')
disp('No timecourses found, processing...');
gmtc=ea_extract_timecourses(options);
else
load([expfolder,options.prefs.gmtc]);
end
fMRI_CM=corrcoef(gmtc);