Home > functions > job > vb_job_flat.m

vb_job_flat

PURPOSE ^

---

SYNOPSIS ^

function vb_job_flat(proj_root,flat_parm)

DESCRIPTION ^

 ---
 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)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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');

Generated on Tue 27-Aug-2013 11:46:04 by m2html © 2005