pop_plottopo.m

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

% pop_plottopo() - plot concatenated multichannel data epochs in a
%                  topographic array format using plottopo()
% Usage:
%   >> pop_plottopo( EEG ); % pop-up
%   >> pop_plottopo( EEG, channels );
%   >> pop_plottopo( EEG, channels, title, singletrials);
%   >> pop_plottopo( EEG, channels, title, singletrials, axsize, ...
%                         'color', ydir, vert);
%
% Inputs:
%   EEG          - input dataset
%   channels     - indices of channels to plot
%   title        - plot title. Default is none.
%   singletrials - [0|1], 0 plot average, 1 plot individual
%                  single trials. Default is 0.
%   others...    - additional plottopo arguments {'axsize', 'color', 'ydir'
%                  'vert'} (see >> help plottopo)
%
% Author: Arnaud Delorme, CNL / Salk Institute, 10 March 2002
%
% See also: plottopo()

%123456789012345678901234567890123456789012345678901234567890123456789012

% Copyright (C) 10 March 2002 Arnaud Delorme, Salk Institute, [e-mail]
%
% 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: pop_plottopo.m,v $
% Revision 1.14  2004/04/06 17:31:33  arno
% adding edit box for options
%
% Revision 1.13  2003/05/10 02:30:01  arno
% continuous array
%
% Revision 1.12  2003/03/12 03:18:58  arno
% help button
%
% Revision 1.11  2003/02/21 00:33:09  scott
% edit header
%
% Revision 1.10  2002/08/29 18:34:29  arno
% changing default background
%
% Revision 1.9  2002/08/14 01:39:30  scott
% figure name plottopo()
%
% Revision 1.8  2002/08/12 22:27:52  arno
% change title
%
% Revision 1.7  2002/08/12 16:34:05  arno
% same
%
% Revision 1.6  2002/08/12 02:35:29  arno
% [6~[6~inputdlg2
%
% Revision 1.5  2002/08/12 01:39:10  arno
% color
%
% Revision 1.4  2002/08/11 22:18:23  arno
% color
%
% Revision 1.3  2002/04/24 18:56:29  arno
% adding vert and others
%
% Revision 1.2  2002/04/18 15:45:11  scott
% editted msgs -sm
%
% Revision 1.1  2002/04/05 17:32:13  jorn
% Initial revision
%

% 01-25-02 reformated help & license -ad 
% 02-16-02 text interface editing -sm & ad 
% 03-18-02 added title -ad & sm
% 03-30-02 added single trial capacities -ad

function com = pop_plottopo( EEG, channels, plottitle, singletrials, varargin);

com = '';
if nargin < 1
	help pop_plottopo;
	return;
end;	

if isempty(EEG.chanlocs)
	fprintf('Cannot plot without knowing channel locations. Use Edit/Dataset info\n');
	return;
end;

if nargin < 2
	promptstr    = { 'Channels to plot:' ...
					 'Plot title:' ...
					 'Plot single trials (yes|no)' ... 
                     'Vertical lines (ms)' ...
                     'Plotting options (see help)' };
	inistr       = { [ '1:' num2str( EEG.nbchan ) ] ...
					 fastif(isempty(EEG.setname), '',EEG.setname) ...
					 'no' '' '''ydir'', -1'};
	result       = inputdlg2( promptstr, 'Topographic ERP plot - pop_plottopo()', 1, inistr, 'pop_plottopo');
	if size(result,1) == 0 return; end;
	channels     = eval( [ '[' result{1} ']' ] );
	plottitle    = result{2};
	singletrials = strcmp( lower(result{3}), 'yes');
	vert         = eval( [ '[' result{4} ']' ] );
    addoptions   = eval( [ '{' result{5} '}' ] );
    figure('name', ' plottopo()');
	if ~isempty(vert)
		addoptions ={ addoptions{:} 'vert' vert };
	end;
    options ={ 'chanlocs' EEG.chanlocs 'frames' EEG.pnts 'limits' [EEG.xmin EEG.xmax 0 0]*1000 ...
               'title' plottitle 'chans' channels addoptions{:} };
else 
	options ={ 'chanlocs' EEG.chanlocs 'frames' EEG.pnts 'limits' [EEG.xmin EEG.xmax 0 0]*1000 ...
               'title' plottitle 'chans' channels varargin{:}};
    addoptions = {};
end;
try, icadefs; set(gcf, 'color', BACKCOLOR); catch, end;
	
if exist('plottitle') ~= 1
    plottitle = '';
end;    
if exist('singletrials') ~= 1
    singletrials = 0;
end;    

if singletrials 
    plottopo( EEG.data, options{:} );
else
    plottopo( mean(EEG.data,3), options{:} );
end;

if ~isempty(addoptions)
	com = sprintf('figure; pop_plottopo(%s, %s, ''%s'', %d, %s);', ...
				  inputname(1), vararg2str(channels), plottitle, singletrials, vararg2str(addoptions));
else
	com = sprintf('figure; pop_plottopo(%s, %s, ''%s'', %d);', ...
				  inputname(1), vararg2str(channels), plottitle, singletrials);
end;

return;