compplot.m

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

% compplot() - plot a data epoch and maps its scalp topography at a given time
%
% Usage: To plot the projection of an ICA component onto the scalp
%        >> projdata = icaproj(data,weights,compindex);
%
% then   >> compplot(projdata);
%
% else to plot an EEG epoch with a topoplot at one selected time point
%        >> compplot(data,plotframe,chan_file,xstart,srate,title, splinefile);
%
% Inputs:
%  data        = data returned from icaproj() *ELSE* any EEG/ERP data epoch
%  plotframe   = frame to plot topographically {default|0 -> frame of max(var())}
%  'chan_file' = chan file, see >> topoplot example {def|0 -> 'chan_file'}
%  xstart      = start time in seconds {default|0 -> 0}
%  srate       = data sampling rate in Hz {default|0 -> 256 Hz}
%  'title'     = plot title {default|0 -> none}
%  splinefile  = headplot spline file (optional) for 3-d head image {default: none}
%
% Authors: Colin Humphries & Scott Makeig, SCCN/INC/UCSD, CNL / Salk Institute, 1997 
%
% See also: icaproj()

%123456789012345678901234567890123456789012345678901234567890123456789012

% Copyright (C) 2000 Colin Humphries & Scott Makeig, SCCN/INC/UCSD
%
% 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: compplot.m,v $
% Revision 1.1  2002/04/05 17:36:45  jorn
% Initial revision
%

% 6-07-97 changed order of args in icaproj above -sm
% Revised for Matlab 5.0.0.4064 11/97 -ch & sm
% 2-18-98 added undocumented headplot() option -sm
% 3-09-98 added check for chan_file -sm
% 3-18-98 use new eegplot('noui') -sm
% 7-25-98 make sure length(hkids) >= 3 on line 122 ff -sm
% 12-19-00 updated icaproj() call in help msg -sm
% 1-24-01 corrected input error in header -ad
% 01-25-02 reformated help & license, added links -ad 

function M = compplot(data,plotframe,chan_file,xstart,srate,titl,spline_file)

if nargin < 1
   help compplot
   return
end

[chans,frames] = size(data);
icadefs;   % read DEFAULT_SRATE

if nargin < 7
    spline_file = nan;
end
if nargin < 6,
   titl = '';     % DEFAULT TITLE
end
if nargin < 5
    srate = 0;  
end
if nargin < 4
    xstart = 0;   % DEFAULT START TIME
end
if nargin < 3
    chan_file = 'chan_file';  % DEFAULT CHAN_FILE
end
if chan_file == 0,
    chan_file = 'chan_file';  % DEFAULT CHAN_FILE
end
if nargin < 2
    plotframe = 0;
end
if plotframe == 0
    [mx plotframe] = max(mean(data.*data)); 
                         % default plotting frame has max variance
end
if plotframe > frames
    fprintf('Plot frame %d is > frames in data (%d)\n',plotframe,frames);
    return
end
fprintf('Topoplot will show frame %d\n',plotframe);
if nargin < 1
    help compplot
end

if srate==0,
	srate = DEFAULT_SRATE;
end

clf % clear the current figure
set(gcf,'Color',BACKCOLOR); % set the background color
axe = axes('Units','Normalized','Position',[.75 .12 .2 .8]);
set(gca,'Color','none');

% >> eegplot('noui',data,srate,spacing,eloc_file,startsec,color)
eegplot('noui',-data(:,xstart*srate+1:end),srate,0,chan_file,0,'r')

plottime = xstart + (plotframe-1)/srate;
timetext = num2str(plottime,'%4.3f');
set(gca,'XTick',plotframe,'XtickLabel',timetext)       %%CJH
set(gca,'Color','none');

limits = get(gca,'Ylim');
set(gca,'GridLineStyle',':')
set(gca,'Xgrid','off')
set(gca,'Ygrid','on')

% axes(axe)
plottime = xstart + (plotframe-1)/srate;
l1 = line([plotframe plotframe],limits,'color','b'); % was [plottime plottime]

fid = fopen(chan_file); % check whether topoplot will find chan_fil
if fid<1,
  fprintf('topoplot()^G: cannot open eloc_file (%s).\n',chan_file);
  return
end

axcolor = get(gcf,'Color');
axt = axes('Units','Normalized','Position',[.00 .10 .75 .80],'Color',axcolor); 
axes(axt)                                                       % topoplot axes
cla
if isnan(spline_file)
   topoplot(data(:,plotframe),chan_file,'style','both');
else
   headplot(data(:,plotframe),spline_file);
end
text(0.00,0.70,titl,'FontSize',14,'HorizontalAlignment','Center');

axt = axes('Units','Normalized','Position',[.05 .05 .055 .18]); % colorbar axes
h=cbar(axt);  

% hkids=get(h,'children');
% delno = 3;
% if delno > length(hkids)
%    delno = length(hkids);
% end
% for i=1:delno
%    delete(hkids(i));
% end

axis(axis)
set(h,'XTickLabel',[]);
HYLIM = get(h,'YLim');              %%CJH
set(h,'YTick',[HYLIM(1) (HYLIM(1)+HYLIM(2))/2 HYLIM(2)])  %%CJH
set(h,'YTickLabel',['-';'0';'+']);     %%CJH
set(h,'YTickMode','manual');

axt = axes('Units','Normalized','Position',[0 0 1 1],...
            'Visible','Off','Fontsize',16);                     % topoplot axes
axes(axt)
timetext = [ 't = ',num2str(plottime) ' s'];
text(0.16,0.09,timetext,'FontSize',14);
text(0.85,0.04,'Time (sec)','FontSize',14,'HorizontalAlignment','Center');