vb_load_current_j_tr

 load estimated current

 USAGE:
  [Jall,ix_act] = vb_load_current_j_tr(currfile,trial,ix_area,verbose)

 --- INPUT:
  currfile  - current file name
 --- Optional INPUT:
  curr_type - current type 
            = 0 : J-current : original vertex current [Default]
            = 1 : Z-current : internal current 
                              correspond to locally spread source
     If currfile is old format, curr_type is neglected
     and Jact is returned
  trial   : trial number to load 
            if empty, all trials are loaded
  ix_area : (absolute) vertex index to load
            if empty, all vertices are loaded
  verbose : 0 (No print) , 1 ('.') , 2 ('session , trial')
 --- OUTPUT:
 ix_act  : Vertex index corresponding to current
 Jact    : J-current in focal region
   Jact(Nact,Nsample,Ntrials)  
     Nact     = Lact * Nvact,  Nvact = length(ix_act)
     Nsample  : # of time sample, 
     Ntrials  : # of trials in all session]

   Jact( n + Nvact*(i-1), t, :) 
    current at the vertex 'ix_act(n)', 'i-th' direction, time index 't'

   Time sample in MEG data : t_meg = Tsample(t)
   Time (ms)   in MEG data : t_ms  = (t_meg - Pretrigger)*(1000/SampleFreq)


 2007-3-5 Masa-aki Sato
 2008-7-9 Masa-aki Sato

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

0001 function [Jall,ix_act] = vb_load_current_j_tr(currfile,trial,ix_area,verbose)
0002 % load estimated current
0003 %
0004 % USAGE:
0005 %  [Jall,ix_act] = vb_load_current_j_tr(currfile,trial,ix_area,verbose)
0006 %
0007 % --- INPUT:
0008 %  currfile  - current file name
0009 % --- Optional INPUT:
0010 %  curr_type - current type
0011 %            = 0 : J-current : original vertex current [Default]
0012 %            = 1 : Z-current : internal current
0013 %                              correspond to locally spread source
0014 %     If currfile is old format, curr_type is neglected
0015 %     and Jact is returned
0016 %  trial   : trial number to load
0017 %            if empty, all trials are loaded
0018 %  ix_area : (absolute) vertex index to load
0019 %            if empty, all vertices are loaded
0020 %  verbose : 0 (No print) , 1 ('.') , 2 ('session , trial')
0021 % --- OUTPUT:
0022 % ix_act  : Vertex index corresponding to current
0023 % Jact    : J-current in focal region
0024 %   Jact(Nact,Nsample,Ntrials)
0025 %     Nact     = Lact * Nvact,  Nvact = length(ix_act)
0026 %     Nsample  : # of time sample,
0027 %     Ntrials  : # of trials in all session]
0028 %
0029 %   Jact( n + Nvact*(i-1), t, :)
0030 %    current at the vertex 'ix_act(n)', 'i-th' direction, time index 't'
0031 %
0032 %   Time sample in MEG data : t_meg = Tsample(t)
0033 %   Time (ms)   in MEG data : t_ms  = (t_meg - Pretrigger)*(1000/SampleFreq)
0034 %
0035 %
0036 % 2007-3-5 Masa-aki Sato
0037 % 2008-7-9 Masa-aki Sato
0038 %
0039 % Copyright (C) 2011, ATR All Rights Reserved.
0040 % License : New BSD License(see VBMEG_LICENSE.txt)
0041 
0042 if ~exist('ix_area','var'), ix_area = []; end;
0043 if ~exist('trial','var'),   trial = []; end;
0044 if ~exist('verbose','var'), verbose=2; end;
0045 
0046 load(currfile)
0047 
0048 [Jinfo] = vb_load_current_info(currfile);
0049 
0050 if ~isfield(Jinfo,'jactdir') 
0051     error(['No trial data directory in' currfile])
0052 end
0053 
0054 jactdir = [fileparts(currfile)  '/' Jinfo.jactdir];
0055 
0056 if exist(jactdir,'dir')==0,
0057     error('No trial data directory')
0058 end
0059 
0060 if ~exist('Jact','var') 
0061     error(['No J-current in ' currfile])
0062 end
0063 
0064 [NJ, T] = size(Jact);
0065 
0066 % Select vertex index 'ix_area' within the active current region
0067 [ix_act, jx_act] = vb_current_area_info(Jinfo, ix_area);
0068 
0069 Ntrial  = sum(Jinfo.Ntrial);
0070 
0071 if isempty(trial), trial=1:Ntrial; end;
0072 
0073 Ntrial = length(trial);
0074 NJ = length(jx_act);
0075 
0076 Jall = zeros(NJ, T, Ntrial);
0077 
0078 session = [];
0079 ntrial  = [];
0080 
0081 % Session & Trial number list
0082 for n=1:Jinfo.Nsession
0083     Ntry = Jinfo.Ntrial(n);
0084     session = [session, repmat(n,1,Ntry)];
0085     ntrial  = [ntrial, 1:Ntry];
0086 end
0087 
0088 % Trial loop
0089 for n=1:Ntrial
0090     nt = trial(n);
0091     
0092     if verbose==1,
0093         fprintf('.')
0094     elseif verbose==2,
0095         fprintf('Session: %04d , Trial: %04d\n',session(nt), ntrial(nt))
0096     end
0097     
0098     fname = sprintf('data_s%04dt%04d', session(nt), ntrial(nt));
0099     load([jactdir '/' fname], 'Jact');
0100     
0101     Jall(:,:,n) = Jact(jx_act,:);
0102 end
0103 
0104 if verbose==1,
0105     fprintf('\n')
0106 end