dmri_process_cortex_parcel

 
 Copyright (C) 2011, ATR All Rights Reserved.
 License : New BSD License(see VBMEG_LICENSE.txt)

0001 function dmri_process_cortex_parcel(p)
0002 %
0003 % Copyright (C) 2011, ATR All Rights Reserved.
0004 % License : New BSD License(see VBMEG_LICENSE.txt)
0005 
0006 disp(mfilename);
0007 
0008 if isfield(p, 'parcel_method') && strcmpi(p.parcel_method, 'area')
0009    if isfield(p,'ref_parcel_area_file') && ~isempty(p.ref_parcel_area_file)
0010        parcel_area_file = p.ref_parcel_area_file;
0011    else
0012        parcel_area_file = fullfile(p.output_dir, p.vbmeg_dirname, [p.vbmeg_cortex_basename, '_aal.area.mat']);
0013    end
0014 else
0015     parcel_area_file = '';
0016     Nvertex = 2000;
0017 end
0018 
0019 if isfield(p, 'ref_brain_file') && ~isempty(p.ref_brain_file)
0020     brain_file = p.ref_brain_file;
0021 else
0022     brain_file = fullfile(p.output_dir, p.vbmeg_dirname, [p.vbmeg_cortex_basename, '.brain.mat']);
0023 end
0024 
0025 if isfield(p, 'ref_fs_dir') && ~isempty(p.ref_fs_dir)
0026     fs_dir = p.ref_fs_dir;
0027 else
0028     fs_dir = fullfile(p.output_dir, p.fs_dirname, p.fs_sbj_name);
0029 end
0030 
0031 output_dir   = fullfile(p.output_dir, p.vbmeg_dirname);
0032 
0033 %
0034 % --- Read FreeSurfer information
0035 %
0036 fs_info_file = fullfile(output_dir, 'fs_info', 'fs_info.mat');
0037 dmri_fs_info_read(fs_dir, fs_info_file);
0038 
0039 %
0040 % --- Create parcels
0041 %
0042 parcel_file  = fullfile(output_dir, p.parcels_file);
0043 if isempty(parcel_area_file)
0044     dmri_cortex_parcel(brain_file, fs_info_file, Nvertex, parcel_file);
0045 else
0046     dmri_cortex_parcel_from_area(brain_file, fs_info_file, parcel_area_file, parcel_file);
0047 end