dmri_process_fodf_calc

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

0001 function dmri_process_fodf_calc(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 % --- Get area for response function estimation
0010 %
0011 fa_file        = fullfile(p.output_dir, p.dwi_dirname, [p.FA_basename, '.nii.gz']);
0012 resp_area_file = fullfile(p.output_dir, p.fiber_track_dirname, p.fiber_resp_area_file);
0013 
0014 dmri_response_estimation_prepare(fa_file, resp_area_file);
0015 
0016 
0017 %
0018 % --- Calculate fiber orientation density function.
0019 %
0020 dwi_file                = fullfile(p.output_dir, p.dwi_dirname, [p.dwi_basename, '_m.nii.gz']);
0021 bval_file               = fullfile(p.output_dir, p.dwi_dirname, [p.dwi_basename, '_m.bval']);
0022 bvec_file               = fullfile(p.output_dir, p.dwi_dirname, [p.dwi_basename, '_m.bvec']);
0023 struct_brain_on_fa_file = fullfile(p.output_dir, p.dwi_dirname, [p.t1_brain_basename, '_brain_FA.nii.gz']);
0024 output_mif_file         = fullfile(p.output_dir, p.fiber_track_dirname, p.fiber_mif_file);
0025 
0026 dmri_fodf_calc(dwi_file, bval_file, bvec_file, ...
0027                resp_area_file, struct_brain_on_fa_file, output_mif_file);
0028 
0029 
0030 %
0031 % --- Check result
0032 %
0033 
0034 % Check1
0035 cmd = ['fslview ', struct_brain_on_fa_file, ' ', resp_area_file, ' -l Blue'];
0036 dmri_system(cmd);
0037 
0038 % Check2
0039 % dmri_mrview(fa_file);