0001 function construct_brain_model(p)
0002
0003
0004
0005
0006
0007 disp(mfilename);
0008
0009
0010 freesurfer_dir = fullfile(p.proj_root, p.fs_dirname, p.struct_name);
0011 brain_parm = vb_set_brain_parm();
0012 brain_parm.FS_left_file = fullfile(freesurfer_dir, 'bem', 'lh.smoothwm.asc');
0013 brain_parm.FS_right_file = fullfile(freesurfer_dir, 'bem', 'rh.smoothwm.asc');
0014 brain_parm.FS_left_infl_file = fullfile(freesurfer_dir, 'bem', 'lh.inflated.asc');
0015 brain_parm.FS_right_infl_file = fullfile(freesurfer_dir, 'bem', 'rh.inflated.asc');
0016 brain_parm.FS_left_curv_file = fullfile(freesurfer_dir, 'bem', 'lh.curv.asc');
0017 brain_parm.FS_right_curv_file = fullfile(freesurfer_dir, 'bem', 'rh.curv.asc');
0018 brain_parm.FS_left_sphere_file = fullfile(freesurfer_dir, 'bem', 'lh.sphere.reg.asc');
0019 brain_parm.FS_right_sphere_file = fullfile(freesurfer_dir, 'bem', 'rh.sphere.reg.asc');
0020 brain_parm.FS_sphere_key = 'sphere.reg';
0021 brain_parm.FS_left_label_file = fullfile(freesurfer_dir, 'label', 'lh.cortex.label');
0022 brain_parm.FS_right_label_file = fullfile(freesurfer_dir, 'label', 'rh.cortex.label');
0023
0024
0025
0026 brain_parm.registration_mode = 'FS';
0027
0028 brain_dir = fullfile(p.proj_root, p.brain_dirname);
0029 if exist(brain_dir, 'dir') ~= 7
0030 mkdir(brain_dir);
0031 end
0032 brain_parm.analyze_file = fullfile(p.proj_root, p.t1_dirname, ['m' p.struct_name '.nii']);
0033 brain_parm.brain_file = [p.struct_name, '.brain.mat'];
0034 brain_parm.area_file = [p.struct_name, '.area.mat'];
0035 brain_parm.act_file = [p.struct_name, '.act.mat'];
0036
0037
0038 vb_job_brain(brain_dir, brain_parm);
0039
0040