--- function vb_job_flat(flat_parm) flat_parm.brain_file : .brain.mat file flat_parm.FS_left_file : (for counting vertex number in original model) flat_parm.FS_left_flat : FreeSurfer patch file (left cortex) flat_parm.FS_right_flat: FreeSurfer patch file (right cortex) flat_parm.FS_left_curv : FreeSurfer curvature file flat_parm.FS_right_curv: FreeSurfer curvature file flat_parm.flat_file : Flatmap file 2006-02-20 Taku Yoshioka --- Copyright (C) 2011, ATR All Rights Reserved. License : New BSD License(see VBMEG_LICENSE.txt)
0001 function vb_job_flat(proj_root,flat_parm) 0002 % --- 0003 % function vb_job_flat(flat_parm) 0004 % 0005 % flat_parm.brain_file : .brain.mat file 0006 % flat_parm.FS_left_file : (for counting vertex number in original model) 0007 % flat_parm.FS_left_flat : FreeSurfer patch file (left cortex) 0008 % flat_parm.FS_right_flat: FreeSurfer patch file (right cortex) 0009 % flat_parm.FS_left_curv : FreeSurfer curvature file 0010 % flat_parm.FS_right_curv: FreeSurfer curvature file 0011 % flat_parm.flat_file : Flatmap file 0012 % 0013 % 2006-02-20 Taku Yoshioka 0014 % --- 0015 % 0016 % Copyright (C) 2011, ATR All Rights Reserved. 0017 % License : New BSD License(see VBMEG_LICENSE.txt) 0018 0019 proj_root = vb_rm_trailing_slash(proj_root); 0020 0021 % Vertex indices of reduced model 0022 load(flat_parm.brain_file,'BV_index','F'); 0023 0024 % Count number of vertices in the original model 0025 fid = fopen(flat_parm.FS_left_file,'r'); 0026 fgetl(fid); % comment line 0027 tmp = fscanf(fid,'%d %d',2); % # of vertices and patches 0028 NdipoleL0 = tmp(1); 0029 fclose(fid); 0030 0031 % Load flat map 0032 tic; 0033 fprintf('--- Load FreeSurfer patch data '); 0034 [V0L,F0L,C0L,IndxL0] ... 0035 = vb_fs_load_patch(flat_parm.FS_left_flat,... 0036 flat_parm.FS_left_curv); 0037 [V0R,F0R,C0R,IndxR0] ... 0038 = vb_fs_load_patch(flat_parm.FS_right_flat,... 0039 flat_parm.FS_right_curv); 0040 fprintf('%f[sec]\n',toc); 0041 0042 Vinfo.Ndipole = size(V0L,1)+size(V0R,1); 0043 Vinfo.NdipoleL = size(V0L,1); 0044 Vinfo.Npatch = size(F0L,1)+size(F0R,1); 0045 Vinfo.Coord = 'flat'; 0046 0047 Vflat = [V0L; V0R]; 0048 Cflat = [C0L; C0R]; 0049 Fflat.F3L = F0L; 0050 Fflat.F3R = F0R+Vinfo.NdipoleL; 0051 BV_index.Left = IndxL; % prefix 'BV' for compatibility 0052 BV_index.Right = IndxR+NdipoleL0; 0053 0054 vb_fsave([proj_root '/' flat_parm.flat_file,],'Vflat','Fflat',... 0055 'Cflat','BV_index','Vinfo');