0001 function dmri_process_connection_calc(p)
0002
0003
0004
0005
0006 disp(mfilename);
0007
0008 parcels_file = fullfile(p.output_dir, p.vbmeg_dirname, p.parcels_file);
0009 parcel_dir = fullfile(p.output_dir, p.parcels_dirname);
0010 anat_matrix_file = fullfile(p.output_dir, p.fiber_track_dirname, p.anat_matrix_file);
0011 dmri_file = fullfile(p.output_dir, p.connectivity_dirname, p.dmri_file);
0012
0013
0014
0015 dmri_connect_parm_calc(parcels_file, parcel_dir, anat_matrix_file, dmri_file);
0016
0017
0018
0019
0020 if isfield(p, 'ref_brain_file') && ~isempty(p.ref_brain_file)
0021 brain_file = p.ref_brain_file;
0022 else
0023 brain_file = fullfile(p.output_dir, p.vbmeg_dirname, [p.vbmeg_cortex_basename, '.brain.mat']);
0024 end
0025 if isfield(p, 'ref_aal_area_file') && ~isempty(p.ref_aal_area_file)
0026 area_file = p.ref_aal_area_file;
0027 else
0028 area_file = fullfile(p.output_dir, p.vbmeg_dirname, [p.vbmeg_cortex_basename, '_aal.area.mat']);
0029 end
0030 dmri_plot_dmri_connections_on_brain(dmri_file, brain_file, area_file);
0031 dmri_plot_dmri_connections_matrix(dmri_file, brain_file, area_file);