parcel_fa_volume

 Parcel FA volume based on parcelled cortical surface

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

0001 function parcel_fa_volume(p)
0002 % Parcel FA volume based on parcelled cortical surface
0003 %
0004 % Copyright (C) 2011, ATR All Rights Reserved.
0005 % License : New BSD License(see VBMEG_LICENSE.txt)
0006 
0007 disp(mfilename);
0008 
0009 % Convert parcels to label
0010 dmri_root = fullfile(p.proj_root, p.dmri_dirname);
0011 parcel_dir = fullfile(dmri_root, p.parcels_dirname);
0012 parcel_file = fullfile(parcel_dir, p.parcels_file);
0013 fs_info_file = fullfile(parcel_dir, 'fs_info', 'fs_info.mat');
0014 label_dir = fullfile(parcel_dir, 'parcels_label');
0015 
0016 dmri_parcel_to_label_file(parcel_file, fs_info_file, label_dir);
0017 
0018 % Convert label to FA volume.
0019 fs_dir = fullfile(p.proj_root, p.fs_dirname, p.struct_name);
0020 ref_brain_file = fullfile(dmri_root, p.t1_dirname,  [p.struct_name '_brain.nii.gz']);
0021 ref_FA_file = fullfile(dmri_root, p.dmri_dirname, [p.FA_basename '.nii.gz']);
0022 trans_info_dir = fullfile(dmri_root, p.transform_info_dirname);
0023 Nprocess = 10; % The number of Linux processes to make volume files.
0024 
0025 dmri_label_file_to_FA_volume(label_dir, ...
0026                              fs_dir,    ...
0027                              ref_brain_file, ref_FA_file, ...
0028                              trans_info_dir, ...
0029                              parcel_dir, ...
0030                              Nprocess);
0031 
0032 % Check parcellation result
0033 dmri_parcel_view(fs_dir, parcel_dir, 'l');% 'l' : left hemisphere
0034 dmri_parcel_view(fs_dir, parcel_dir, 'r');% 'r' : right hemisphere