tutorial_fiber_tracking_for_movie_for_dbayes

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

0001 %
0002 % Copyright (C) 2011, ATR All Rights Reserved.
0003 % License : New BSD License(see VBMEG_LICENSE.txt)
0004 
0005 root = '/home/cbi/rhayashi/temp/dynamics_movie/data';
0006 
0007 % Load current & MAR matrix
0008 dbayes_file  = fullfile(root, '/dynamics_estimation/estimation_result.dbayes.mat');
0009 load(dbayes_file,  'Model');
0010 Zact = Model.Z{1};
0011 MAR = Model.MAR;
0012 
0013 % Choose strong pairs
0014 Npairs = 100;
0015 current_threshold = 0.0207;
0016 Econ = MAR-diag(diag(MAR));
0017 [parcel_to_list,parcel_from_list] = dmri_select_pairs_for_movie(Zact, Econ, Npairs, current_threshold);
0018 
0019 fprintf('Number of pairs : %d\n', length(parcel_to_list));
0020 
0021 dmri_file        = fullfile(root, '/dmri/connectivity/connectivity.dmri.mat');
0022 parcel_dir       = fullfile(root, '/dmri/parcels');
0023 mif_file         = fullfile(root, '/dmri/fibertrack/CSD6.mif');
0024 mask_file        = fullfile(root, '/dmri/fibertrack/parcel_wm.nii.gz');
0025 fa_gz_file       = fullfile(root, '/dmri/DWI/data_FA.nii.gz');
0026 trans_info_dir   = fullfile(root, '/dmri/transwarp_info');
0027 fs_brain_gz_file = fullfile(root, '/freesurfer/Subject/mri/brain.nii.gz');
0028 
0029 
0030 % Output directory(filename is fixed : Vmni_connect.mat)
0031 output_dir       = '~/dynamics_movie';
0032 
0033 % Fiber tracking hosts(Do not specify a host that runs this script.)
0034 % Available hosts
0035 % - cbi-node20-33
0036 % - cbi-node01g-08g
0037 % - cbi-node20g-23g
0038 process_host     = {'cbi-node01g', 'cbi-node02g','cbi-node03g'};
0039 
0040 % The number of MATLAB processes which are used to process fiber tracking result files.
0041 Nworker          = 20;
0042 
0043 dmri_fiber_track_prob_specific_pair(parcel_dir, parcel_from_list, parcel_to_list, ...
0044                               mif_file, mask_file, ...
0045                               fa_gz_file, fs_brain_gz_file, trans_info_dir, dmri_file, ...
0046                               output_dir, process_host, Nworker);
0047 
0048 
0049