0001
0002
0003
0004 root = '/home/cbi/rhayashi/temp/dynamics_movie/data';
0005
0006
0007 dbayes_file = fullfile(root, '/dynamics_estimation/estimation_result.dbayes.mat');
0008 load(dbayes_file, 'Model');
0009 Zact = Model.Z{1};
0010 MAR = Model.MAR;
0011
0012
0013 Npairs = 100;
0014 current_threshold = 0.0207;
0015 Econ = MAR-diag(diag(MAR));
0016 [parcel_to_list,parcel_from_list] = dmri_select_pairs_for_movie(Zact, Econ, Npairs, current_threshold);
0017
0018 fprintf('Number of pairs : %d\n', length(parcel_to_list));
0019
0020 dmri_file = fullfile(root, '/dmri/connectivity/connectivity.dmri.mat');
0021 parcel_dir = fullfile(root, '/dmri/parcels');
0022 mif_file = fullfile(root, '/dmri/fibertrack/CSD6.mif');
0023 mask_file = fullfile(root, '/dmri/fibertrack/wm.nii.gz');
0024 fa_gz_file = fullfile(root, '/dmri/DWI/data_FA.nii.gz');
0025 trans_info_dir = fullfile(root, '/dmri/transwarp_info');
0026 fs_brain_gz_file = fullfile(root, '/freesurfer/Subject/mri/brain.nii.gz');
0027
0028
0029
0030 output_dir = '~/dynamics_movie';
0031
0032
0033
0034
0035
0036
0037 process_host = {'cbi-node01g', 'cbi-node02g','cbi-node03g'};
0038
0039
0040 Nworker = 20;
0041
0042 dmri_fiber_track_prob_specific_pair(parcel_dir, parcel_from_list, parcel_to_list, ...
0043 mif_file, mask_file, ...
0044 fa_gz_file, fs_brain_gz_file, trans_info_dir, dmri_file, ...
0045 output_dir, process_host, Nworker);
0046
0047
0048