dmri_process_FA_parcel_volume_files_create

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

0001 function dmri_process_FA_parcel_volume_files_create(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 %
0009 % --- Convert parcels to label
0010 %
0011 parcel_file    = fullfile(p.output_dir, p.vbmeg_dirname, p.parcels_file);
0012 fs_info_file   = fullfile(p.output_dir, p.vbmeg_dirname, 'fs_info/fs_info.mat');
0013 label_dir      = fullfile(p.output_dir, p.parcels_dirname, 'parcels_label');
0014 
0015 dmri_parcel_to_label_file(parcel_file, fs_info_file, label_dir);
0016 
0017 
0018 %
0019 % --- Convert label to FA volume.
0020 %
0021 if isfield(p, 'ref_fs_dir') && ~isempty(p.ref_fs_dir)
0022     fs_dir = p.ref_fs_dir;
0023 else
0024     fs_dir = fullfile(p.output_dir, p.fs_dirname, p.fs_sbj_name);
0025 end
0026 ref_struct_file  = fullfile(p.output_dir, p.t1_dirname,  [p.t1_brain_basename, '_brain.nii.gz']);
0027 ref_FA_file      = fullfile(p.output_dir, p.dwi_dirname, [p.FA_basename '.nii.gz']);
0028 trans_info_dir   = fullfile(p.output_dir, p.transform_info_dirname);
0029 output_dir       = fullfile(p.output_dir, p.parcels_dirname);
0030 Nprocess         = 10; % The number of Linux processes to make volume files.
0031 
0032 dmri_label_file_to_FA_volume(label_dir, ...
0033                              fs_dir,    ...
0034                              ref_struct_file, ref_FA_file, ...
0035                              trans_info_dir, ...
0036                              output_dir, ...
0037                              Nprocess);
0038 
0039 %
0040 % --- Check parcellation result
0041 %
0042 dmri_parcel_view(fs_dir, output_dir, 'l'); % 'l' : left hemisphere
0043                                           % 'r' : right hemisphere