Home > vbmeg > demo > tutorial_for_vbmeg2 > advanced > construct_brain_model.m

construct_brain_model

PURPOSE ^

Construct brain model

SYNOPSIS ^

function construct_brain_model(p)

DESCRIPTION ^

 Construct brain model

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function construct_brain_model(p)
0002 % Construct brain model
0003 %
0004 % Copyright (C) 2011, ATR All Rights Reserved.
0005 % License : New BSD License(see VBMEG_LICENSE.txt)
0006 
0007 disp(mfilename);
0008 
0009 % Set parameters
0010 freesurfer_dir = fullfile(p.proj_root, p.fs_dirname, p.struct_name);
0011 brain_parm = vb_set_brain_parm(); % get default parameter
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 % registration method to standard brain model.
0025 %     = 'FS' : Freesurfer sphere file is used.
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 % Construct brain model
0038 vb_job_brain(brain_dir, brain_parm);
0039 
0040

Generated on Mon 22-May-2023 06:53:56 by m2html © 2005