loadtxt.m

含有多种ICA算法的eeglab工具箱

% loadtxt() - load ascii text file into numeric or cell arrays
%
% Usage:
%   >> array = loadtxt( filename, 'key', 'val' ...);
%
% Inputs:
%    filename - name of the input file
%
% Optional inputs
%   'skipline' - number of lines to skip {default:0}. If this number is
%                negative the program will only skip non-empty lines 
%                (can be usefull for files transmitted from one platform
%                to an other, as CR may be inserted at every lines).
%   'convert'  - 'on' standard text conversion, see note 1
%                'off' no conversion, considers text only
%                'force' force conversion, NaN are returned 
%                for non-numeric inputs {default:'on'}
%   'delim'    - ascii character for delimiters. {default:[9 32]
%                i.e space and tab}. It is also possible to enter 
%                strings, Ex: [9 ' ' ','].
%   'verbose'  - ['on'|'off'] {default:'on'}
%   'nlines'   - [integer] number of lines to read {default: all file}
%
% Outputs:
%    array - cell array. If the option 'force' is given, the function
%            retrun a numeric array.
%
% Notes: 1) Since it uses cell arrays, the function can handle text input.
%        The function reads each token and then try to convert it to a 
%        number. If the conversion is unsucessfull, the string itself
%        is included in the array.
%        2) The function adds empty entries for rows that contains
%        fewer columns than others.
%
% Author: Arnaud Delorme, CNL / Salk Institute, 29 March 2002

%123456789012345678901234567890123456789012345678901234567890123456789012

% Copyright (C) Arnaud Delorme, CNL / Salk Institute, 29 March 2002
%
% This program is free software; you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation; either version 2 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program; if not, write to the Free Software
% Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

% $Log: loadtxt.m,v $
% Revision 1.6  2004/01/29 01:35:35  arno
% debug numerical conversion
%
% Revision 1.5  2003/11/19 19:28:16  arno
% now reading empty tabs
%
% Revision 1.4  2003/11/07 01:45:32  arno
% adding nline argument
%
% Revision 1.3  2003/01/10 17:28:56  arno
% debug last
%
% Revision 1.2  2003/01/10 17:27:13  arno
% str2num -> str2double
%
% Revision 1.1  2002/04/05 17:39:45  jorn
% Initial revision
%

function array = loadtxt( filename, varargin );

if nargin < 1
	help loadtxt;
	return;
end;	
if ~isempty(varargin)
   try, g = struct(varargin{:});
   catch, disp('Wrong syntax in function arguments'); return; end;
else
    g = [];
end;

try, g.convert;	 	  catch, g.convert = 'on'; end;
try, g.skipline;      catch, g.skipline = 0; end;
try, g.verbose;       catch, g.verbose = 'on'; end;
try, g.delim; 	      catch, g.delim = [9 32]; end;
try, g.nlines; 	      catch, g.nlines = Inf; end;
g.convert = lower(g.convert);
g.verbose = lower(g.verbose);
g.delim = char(g.delim);

% open the file
% -------------
if exist(filename) ~=2, error( ['file ' filename ' not found'] ); end;  
fid=fopen(filename,'r','ieee-le');
if fid<0, error( ['file ' filename ' found but error while opening file'] ); end;  

index = 0;
while index < abs(g.skipline)
    tmpline = fgetl(fid); 
    if g.skipline > 0 | ~isempty(tmpline)
        index = index + 1;
    end;    
end; % skip lines ---------

inputline = fgetl(fid);
linenb = 1;
if strcmp(g.verbose, 'on'), fprintf('Reading file (lines): '); end;
while isempty(inputline) | inputline~=-1
     colnb = 1;
     if ~isempty(inputline)
	     switch g.convert
	        case 'off',
			     while ~isempty(deblank(inputline))
			         [array{linenb, colnb} inputline] = strtok(inputline, g.delim);
			         colnb = colnb+1;
			     end;
	        case 'on',
			     while ~isempty(deblank(inputline))
			         [tmp inputline] = mystrtok(inputline, g.delim);
			         if ~isempty(tmp) & tmp(1) > 43 & tmp(1) < 59, tmp2 = str2num(tmp);
                     else tmp2 = []; end;
			         if isempty( tmp2 )  , array{linenb, colnb} = tmp;
			         else                  array{linenb, colnb} = tmp2;
			         end;
			         colnb = colnb+1;
			     end;
	        case 'force',
			     while ~isempty(deblank(inputline))
			         [tmp inputline] = mystrtok(inputline, g.delim);
			         array{linenb, colnb} = str2double( tmp );
			         colnb = colnb+1;
			     end;
	        otherwise, error('Unrecognised converting option');
	     end;   
	     linenb = linenb +1;
     end;
     inputline = fgetl(fid);
     if linenb > g.nlines
         inputline = -1;
     end;
     if ~mod(linenb,10) & strcmp(g.verbose, 'on'), fprintf('%d ', linenb); end;
end;        
if strcmp(g.verbose, 'on'),  fprintf('%d\n', linenb-1); end;
if strcmp(g.convert, 'force'), array = cell2mat(array); end;
fclose(fid); 

% problem strtok do not consider tabulation
% -----------------------------------------
function [str, strout] = mystrtok(strin, delim);
    if delim == 9 % tab
        if length(strin) > 1 & strin(1) == 9 & strin(2) == 9 
            str = '';
            strout = strin(2:end);
        else
            [str, strout] = strtok(strin, delim);
        end;
    else
        [str, strout] = strtok(strin, delim);
    end;