% TFTOPO - Generate a figure showing a selected or representative image (e.g.,
% an ERSP, ITC or ERP-image) from a supplied set of images, one for each
% scalp channel. Then, plot TOPOPLOT scalp maps of value distributions
% at specified (time, frequency) image points. Else, image the signed
% (selected) between-channel STD. Inputs may be outputs of
% TIMEF, CROSSF, or ERPIMAGE.
% Usage:
% >> tftopo(tfdata,times,freqs, 'key1', 'val1', 'key2', val2' ...)
% Inputs:
% tfdata = Set of time/freq images, one for each channel. Matrix dims:
% (time,freq),(time,freq,chans). Else, (time,freq,chans,subjects) for grand mean
% RMS plotting.
% times = Vector of image (x-value) times in msec, from TIMEF).
% freqs = Vector of image (y-value) frequencies in Hz, from TIMEF).
%
% Optional inputs:
% 'timefreqs' = Array of time/frequency points at which to plot TOPOPLOT maps.
% Size: (nrows,2), each row given the [ms Hz] location
% of one point. Or size (nrows,4), each row given [min_ms
% max_ms min_hz max_hz].
% 'showchan' = [integer] Channel number of the tfdata to image. Else 0 to image
% the (median-signed) RMS values across channels. {default: 0}
% 'chanlocs' = ['string'|structure] Electrode locations file (for format, see
% >> topoplot example) or EEG.chanlocs structure {default: none}
% 'limits' = Vector of plotting limits [minms maxms minhz maxhz mincaxis maxcaxis]
% May omit final vales, or use NaN's to use the input data limits.
% Ex: [nan nan -100 400];
% 'signifs' = (times,freqs) Matrix of significance level(s) (e.g., from TIMEF)
% to zero out non-signif. tfdata points. Matrix size must be
% ([1|2], freqs, chans, subjects)
% if using the same threshold for all time points at each frequency, or
% ([1|2], freqs, times, chans, subjects).
% If first dimension is of size 1, data are assumed to contain
% positive values only {default: none}
% 'sigthresh' = [K L] After masking time-frequency decomposition using the 'signifs'
% array (above), concatenate (time,freq) values for which no more than
% K electrodes have non-0 (significant) values. If several subjects,
% the second value L is used to concatenate subjects in the same way.
% {default: [1 1]}
% 'selchans' = Channels to include in the TOPOPLOT scalp maps (and image values)
% {default: all}
% 'smooth' = [pow2] magnification and smoothing factor. power of 2 (default: 1}.
% 'mode' = ['rms'|'ave'] ('rms') return root-mean-square, else ('ave') average
% power {default: 'rms' }
% 'logfreq' = ['on'|'off'|'native'] plot log frequencies {default: 'off'}
% 'native' means that the input is already in log frequencies
% 'vert' = [times vector] (in msec) plot vertical dashed lines at specified times
% {default: 0}
% 'ylabel' = [string] label for the ordinate axis. Default is
% "Frequency (Hz)"
% 'shiftimgs' = [response_times_vector] shift time/frequency images from several
% subjects by each subject's response time {default: no shift}
% 'title' = [quoted_string] plot title (default: provided_string).
% 'cbar' = ['on'|'off'] plot color bar {default: 'on'}
% 'cmode' = ['common'|'separate'] 'common' or 'separate' color axis for each
% topoplot {default: 'common'}
% 'plotscalponly' = [x,y] location (e.g. msec,hz). Plot one scalp map only; no
% time-frequency image.
% 'events' = [real array] plot event latencies. The number of event
% must be the same as the number of "frequecies".
% 'verbose' = ['on'|'off'] comment on operations on command line {default: 'on'}.
% 'axcopy' = ['on'|'off'] creates a copy of the figure axis and its graphic objects in a new pop-up window
% using the left mouse button {default: 'on'}..
% 'denseLogTicks' = ['on'|'off'] creates denser labels on log freuqncy axis {default: 'off'}
%
% Notes:
% 1) Additional TOPOPLOT optional arguments can be used.
% 2) In the topoplot maps, average power (not masked by significance) is used
% instead of the (signed and masked) root-mean-square (RMS) values used in the image.
% 3) If tfdata from several subjects is used (via a 4-D tfdata input), RMS power is first
% computed across electrodes, then across the subjects.
%
% Authors: Scott Makeig, Arnaud Delorme & Marissa Westerfield, SCCN/INC/UCSD, La Jolla, 3/01
%
% See also: TIMEF, TOPOPLOT, SPECTOPO, TIMTOPO, ENVTOPO, CHANGEUNITS
% hidden parameter: 'shiftimgs' = array with one value per subject for shifting in time the
% time/freq images. Had to be inserted in tftopo because
% the shift happen after the smoothing
% Copyright (C) Scott Makeig, Arnaud Delorme & Marissa Westerfield, SCCN/INC/UCSD,
% La Jolla, 3/01
%
% This file is part of EEGLAB, see http://www.eeglab.org
% for the documentation and details.
%
% Redistribution and use in source and binary forms, with or without
% modification, are permitted provided that the following conditions are met:
%
% 1. Redistributions of source code must retain the above copyright notice,
% this list of conditions and the following disclaimer.
%
% 2. Redistributions in binary form must reproduce the above copyright notice,
% this list of conditions and the following disclaimer in the documentation
% and/or other materials provided with the distribution.
%
% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
% ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
% LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
% CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
% SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
% INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
% CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
% ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
% THE POSSIBILITY OF SUCH DAMAGE.
% 01-25-02 reformated help & license -ad
function tfave = tftopo(tfdata,times,freqs,varargin)
%timefreqs,showchan,chanlocs,limits,signifs,selchans)
LINECOLOR= 'k';
LINEWIDTH = 2;
ZEROLINEWIDTH = 2.8;
if nargin<3
help tftopo
return
end
icadefs_flag = 1;
try
icadefs;
catch
warning('icadefs.m can not be located in the path');
icadefs_flag = 0 ;
end
if ~icadefs_flag
AXES_FONTSIZE = 10;
PLOT_LINEWIDTH = 2;
end
% reshape tfdata
% --------------
if length(size(tfdata))==2
if size(tfdata,1) ~= length(freqs), tfdata = tfdata'; end
nchans = round(size(tfdata,2)/length(times));
tfdata = reshape(tfdata, size(tfdata,1), length(times), nchans);
elseif length(size(tfdata))>=3
nchans = size(tfdata,3);
else
help tftopo
return
end
% for topoplot
if length(size(tfdata)) >= 4
tfdataori = mean(tfdata,4);
else
tfdataori = tfdata;
end
% test inputs
% -----------
% 'key' 'val' sequence
fieldlist = { 'chanlocs' { 'string','struct' } [] '' ;
'limits' 'real' [] [nan nan nan nan nan nan];
'logfreq' 'string' {'on','off','native'} 'off';
'cbar' 'string' {'on','off' } 'on';
'mode' 'string' { 'ave','rms' } 'rms';
'title' 'string' [] '';
'verbose' 'string' {'on','off' } 'on';
'axcopy' 'string' {'on','off' } 'on';
'cmode' 'string' {'common','separate' } 'common';
'selchans' 'integer' [1 nchans] [1:nchans];
'shiftimgs' 'real' [] [];
'plotscalponly' 'real' [] [];
'events' 'real' [] [];
'showchan' 'integer' [0 nchans] 0 ;
'signifs' 'real' [] [];
'sigthresh' 'integer' [1 Inf] [1 1];
'smooth' 'real' [0 Inf] 1;
'timefreqs' 'real' [] [];
'ylabel' 'string' {} 'Frequency (Hz)';
'vert' 'real' [times(1) times(end)] [min(max(0, times(1)), times(end))];
'denseLogTicks' 'string' {'on','off'} 'off'
};
[g, varargin] = finputcheck( varargin, fieldlist, 'tftopo', 'ignore');
if ischar(g), error(g); end
% setting more defaults
% ---------------------
if length(times) ~= size(tfdata,2)
fprintf('tftopo(): tfdata columns must be a multiple of the length of times (%d)\n',...
length(times));
return
end
if length(g.showchan) > 1
error('tftopo(): showchan must be a single number');
end
if length(g.limits)<1 || isnan(g.limits(1))
g.limits(1) = times(1);
end
if length(g.limits)<2 || isnan(g.limits(2))
g.limits(2) = times(end);
end
if length(g.limits)<3 || isnan(g.limits(3))
g.limits(3) = freqs(1);
end
if length(g.limits)<4 || isnan(g.limits(4))
g.limits(4) = freqs(end);
end
if length(g.limits)<5 || isnan(g.limits(5)) % default caxis plotting limits
g.limits(5) = -max(abs(tfdata(:)));
mincax = g.limits(5);
end
if length(g.limits)<6 || isnan(g.limits(6))
defaultlim = 1;
if exist('mincax')
g.limits(6) = -mincax; % avoid recalculation
else
g.limits(6) = max(abs(tfdata(:)));
end
else
defaultlim = 0;
end
if length(g.sigthresh) == 1
g.sigthresh(2) = 1;
end
if g.sigthresh(1) > nchans
error('tftopo(): ''sigthresh'' first number must be less than or equal to the number of channels');
end
if g.sigthresh(2) > size(tfdata,4)
error('tftopo(): ''sigthresh'' second number must be less than or equal to the number of subjects');
end
if ~isempty(g.signifs)
if size(g.signifs,1) > 2 || size(g.signifs,2) ~= size(tfdata,1)|| ...
(size(g.signifs,3) ~= size(tfdata,3) && size(g.signifs,4) ~= size(tfdata,3))
fprintf('tftopo(): error in ''signifs'' array size not compatible with data size, trying to transpose.\n');
g.signifs = permute(g.signifs, [2 1 3 4]);
if size(g.signifs,1) > 2 || size(g.signifs,2) ~= size(tfdata,1)|| ...
(size(g.signifs,3) ~= size(tfdata,3) && size(g.signifs,4) ~= size(tfdata,3))
fprintf('tftopo(): ''signifs'' still the wrong size.\n');
return
end
end
end
if length(g.selchans) ~= nchans
selchans_opt = { 'plotchans' g.selchans };
else selchans_opt = { };
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% process time/freq data points
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ~isempty(g.timefreqs)
if size(g.timefreqs,2) == 2
g.timefreqs(:,3) = g.timefreqs(:,2);
g.timefreqs(:,4) = g.timefreqs(:,2);
g.timefreqs(:,2) = g.timefreqs(:,1);
end
if isempty(g.chanlocs)
error('tftopo(): ''chanlocs'' must be defined to plot time/freq points');
end
if min(min(g.timefreqs(:,[3 4])))<min(freqs)
fprintf('tftopo(): selected plotting frequency %g out of range.\n',min(min(g.timefreqs(:,[3 4]))));
return
end
if max(max(g.timefreqs(:,[3 4])))>max(freqs)
fprintf('tftopo(): selected plotting frequency %g out of range.\n',max(max(g.timefreqs(:,[3 4]))));
return
end
if min(min(g.timefreqs(:,[1 2])))<min(times)
fprintf('tftopo(): selected plotting time %g out of range.\n',min(min(g.timefreqs(:,[1 2]))));
return
end
if max(max(g.timefreqs(:,[1 2])))>max(times)
fprintf('tftopo(): selected plotting time %g out of range.\n',max(max(g.timefreqs(:,[1 2]))));
return
end
if 0 % USE USER-SUPPLIED SCALP MAP ORDER. A GOOD ALGORITHM FOR SELECTING
% g.timefreqs POINT ORDER GIVING MAX UNCROSSED LINES IS DIFFICULT!
[tmp tfi] = sort(g.timefreqs(:,1)); % sort on times
tmp = g.timefreqs;
for t=1:size(g.timefreqs,1)
g.timefreqs(t,:) = tmp(tfi(t),:);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Compute timefreqs point indices
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tfpoints = size(g.timefreqs,1);
for f=1:tfpoints
[~, fi1] = min(abs(freqs-g.timefreqs(f,3)));
[~, fi2] = min(abs(freqs-g.timefreqs(f,4)));
freqidx{f}=[fi1:fi2];
end
for f=1:tfpoints
[~, fi1] = min(abs(times-g.timefreqs(f,1)));
[~, fi2] = min(abs(times-g.timefreqs(f,2)));
timeidx{f}=[fi1:fi2];
end
else
tfpoints = 0;
end
% only plot one scalp map
% -----------------------
if ~isempty(g.plotscalponly)
[~, fi] = min(abs(freqs-g.plotscalponly(2)));
[~, ti] = min(abs(times-g.plotscalponly(1)));
scalpmap = squeeze(tfdataori(fi, ti, :));
if ~isempty(varargin)
topoplot(scalpmap,g.chanlocs,'electrodes','on', selchans_opt{:}, varargin{:});
else
topoplot(scalpmap,g.chanlocs,'electrodes','on', selchans_opt{:});
end
% 'interlimits','electrodes')
axis square;
hold on
tl=title([int2str(g.plotscalponly(2)),' ms, ',int2str(g.plotscalponly(1)),' Hz'],'interpreter','none');
set(tl,'fontsize',AXES_FONTSIZE+3); % 13
return;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Zero out non-significant image features
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
range = g.limits(6)-g.limits(5);
cc = feval(DEFAULT_COLORMAP, 256);
if ~isempty(g.signifs)
if strcmpi(g.verbose, 'on')
fprintf('Applying ''signifs'' mask by zeroing non-significant values\n');
end
for subject = 1:size(tfdata,4)
for elec = 1:size(tfdata,3)
if size(g.signifs,1) == 2
if ndims(g.signifs) > ndims(tfdata)
tmpfilt = (tfdata(:,:,elec,subject) >= squeeze(g.signifs(2,:,:,elec, subject))') | ...
(tfdata(:,:,elec,subject) <= squeeze(g.signifs(1,:,:,elec, subject))');
else
tmpfilt = (tfdata(:,:,elec,subject) >= repmat(g.signifs(2,:,elec, subject)', [1 size(tfdata,2)])) | ...
(tfdata(:,:,elec,subject) <= repmat(g.signifs(1,:,elec, subject)', [1 size(tfdata,2)]));
end
else
if ndims(g.signifs) > ndims(tfdata)
tmpfilt = (tfdata(:,:,elec,subject) >= squeeze(g.signifs(1,:,:,elec, subject))');
else
tmpfilt = (tfdata(:,:,elec,subject) >= repmat(g.signifs(1,:,elec, subject)', [1 size(tfdata,2)]));
end
end
tfdata(:,:,elec,subject) = tfdata(:,:,elec,subject) .* tmpfilt;
end
end
end
%%%%%%%%%%%%%%%%
% magnify inputs
%%%%%%%%%%%%%%%%
if g.smooth ~= 1
if strcmpi(g.verbose, 'on')
fprintf('Smoothing...\n');
end
for index = 1:round(log2(g.smooth))
[tfdata, times, freqs] = magnifytwice(tfdata, times, freqs);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%
% Shift time/freq images
%%%%%%%%%%%%%%%%%%%%%%%%
if ~isempty(g.shiftimgs)
timestep = times(2) - times(1);
for S = 1:size(tfdata,4)
nbsteps = round(g.shiftimgs(S)/timestep);
if strcmpi(g.verbose, 'on')
fprintf('Shifing images of subect %d by %3.3f ms or %d time steps\n', S, g.shiftimgs(S), nbsteps);
end
if nbsteps < 0, tfdata(:,-nbsteps+1:end,:,S) = tfdata(:,1:end+nbsteps,:,S);
else tfdata(:,1:end-nbsteps,:,S) = tfdata(:,nbsteps+1:end,:,S);
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%
% Adjust plotting limits
%%%%%%%%%%%%%%%%%%%%%%%%%
[~, minfreqidx] = min(abs(g.limits(3)-freqs)); % adjust min frequency
g.limits(3) = freqs(minfreqidx);
[~, maxfreqidx] = min(abs(g.limits(4)-freqs)); % adjust max frequency
g.limits(4) = freqs(maxfreqidx);
[~, mintimeidx] = min(abs(g.limits(1)-times)); % adjust min time
g.limits(1) = times(mintimeidx);
[~, maxtimeidx] = min(abs(g.limits(2)-times)); % adjust max time
g.limits(2) = times(maxtimeidx);
mmidx = [mintimeidx maxtimeidx minfreqidx maxfreqidx];
%colormap('jet');
%c = colormap;
%cc = zeros(256,3);
%if size(c,1)==64
% for i=1:3
% cc(:,i) = interp(c(:,i),4);
% end
%else
% cc=c;
%nd
%cc(find(cc<0))=0;
%cc(find(cc>1))=1;
%if exist('g.signif')
% minnull = round(256*(g.signif(1)-g.limits(5))/range);
% if minnull<1
% minnull = 1;
% end
% maxnull = round(256*(g.signif(2)-g.limits(5))/range);
% if maxnull>256
% maxnull = 256;
% end
% nullrange = minnull:maxnull;
% cc(nullrange,:) = repmat(cc(128,:),length(nullrange),1);
%end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot tfdata image for specified channel or selchans STD
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
axis off;
colormap(cc);
curax = gca; % current plot axes to plot into
if tfpoints ~= 0
plotdim = max(1+floor(tfpoints/2),4); % number of topoplots on top of image
imgax = sbplot(plotdim,plotdim,[plotdim*(plotdim-1)+1,2*plotdim-1],'ax',curax);
else
imgax = curax;
end
tftimes = mmidx(1):mmidx(2);
tffreqs = mmidx(3):mmidx(4);
if g.showchan>0 % -> image showchan data
tfave = tfdata(tffreqs, tftimes,g.showchan);
else % g.showchan==0 -> image std() of selchans
tfdat = tfdata(tffreqs,tftimes,g.selchans,:);
% average across electrodes
if strcmpi(g.verbose, 'on')
if strcmpi( g.mode, 'rms')
fprintf('Applying RMS across channels (mask for at least %d non-zeros values at each time/freq)\n', g.sigthresh(1));
else
fprintf('Computing average across channels (mask for at least %d non-zeros values at each time/freq)\n', g.sigthresh(1));
end
end
tfdat = avedata(tfdat, 3, g.sigthresh(1), g.mode);
% if several subject, first (RMS) averaging across subjects
if size(tfdata,4) > 1
if strcmpi(g.verbose, 'on')
if strcmpi( g.mode, 'rms')
fprintf('Applying RMS across subjects (mask for at least %d non-zeros values at each time/freq)\n', g.sigthresh(2));
else
fprintf('Computing average across channels (mask for at least %d non-zeros values at each time/freq)\n', g.sigthresh(2));
end
end
tfdat = avedata(tfdat, 4, g.sigthresh(2), g.mode);
end
tfave = tfdat;
if defaultlim
g.limits(6) = max(max(abs(tfave)));
g.limits(5) = -g.limits(6); % make symmetrical
end
end
if ~isreal(tfave(1)), tfave = abs(tfave); end
if strcmpi(g.logfreq, 'on')
logimagesc(times(tftimes),freqs(tffreqs),tfave);
axis([g.limits(1) g.limits(2) log(g.limits(3)), log(g.limits(4))]);
elseif strcmpi(g.logfreq, 'native')
imagesc(times(tftimes),log(freqs(tffreqs)),tfave);
axis([g.limits(1:2) log(g.limits(3:4))]);
if g.denseLogTicks
minTick = min(ylim);
maxTick = max(ylim);
set(gca,'ytick',linspace(minTick, maxTick,50));
end
tmpval = get(gca,'yticklabel');
if iscell(tmpval)
ft = cellfun(@str2double, tmpval(:));
else
% MATLAB version < 8.04
ft = str2num(tmpval);
end
ft = exp(1).^ft;
ft = unique_bc(round(ft));
ftick = get(gca,'ytick');
ftick = exp(1).^ftick;
ftick = unique_bc(round(ftick));
ftick = log(ftick);
inds = unique_bc(round(exp(linspace(log(1), log(length(ft))))));
set(gca,'ytick',ftick(inds(1:2:end)),'yticklabel', num2str(ft(inds(1:2:end))));
else
imagesc(times(tftimes),freqs(tffreqs),tfave);
axis([g.limits(1:4)]);
end
if length(unique(g.limits(5:6))) > 1
caxis([g.limits(5:6)]);
end
hold on;
%%%%%%%%%%%%%%%%%%%%%%%%%%
% Title and vertical lines
%%%%%%%%%%%%%%%%%%%%%%%%%%
axes(imgax)
xl=xlabel('Time (ms)');
set(xl,'fontsize',AXES_FONTSIZE+2);%12
set(gca,'yaxislocation','left')
if g.showchan>0
% tl=title(['Channel ',int2str(g.showchan)]);
% set(tl,'fontsize',14);
else
if isempty(g.title)
if strcmpi(g.mode, 'rms')
tl=title(['Signed channel rms']);
else
tl=title(['Signed channel average']);
end
else
tl = title(g.title,'interpreter','none');
end
set(tl,'fontsize',AXES_FONTSIZE + 2); %12
set(tl,'fontweigh','normal');
end
yl=ylabel(g.ylabel);
set(yl,'fontsize',AXES_FONTSIZE + 2); %12
set(gca,'fontsize',AXES_FONTSIZE + 2); %12
set(gca,'ydir','normal');
set(gca, 'userdata', { 'Time (ms)' g.ylabel });
for indtime = g.vert
tmpy = ylim;
htmp = plot([indtime indtime],tmpy,[LINECOLOR ':']);
set(htmp,'linewidth',PLOT_LINEWIDTH)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot topoplot maps at specified timefreqs points
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ~isempty(g.events)
tmpy = ylim;
yvals = linspace(tmpy(1), tmpy(2), length(g.events));
plot(g.events, yvals, 'k', 'linewidth', 2);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot topoplot maps at specified timefreqs points
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ~isempty(g.timefreqs)
wholeax = sbplot(1,1,1,'ax',curax);
topoaxes = zeros(1,tfpoints);
for n=1:tfpoints
if n<=plotdim
topoaxes(n)=sbplot(plotdim,plotdim,n,'ax',curax);
else
topoaxes(n)=sbplot(plotdim,plotdim,plotdim*(n+1-plotdim),'ax',curax);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot connecting lines using changeunits()
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tmptimefreq = [ mean(g.timefreqs(n,[1 2])) mean(g.timefreqs(n,[3 4])) ];
if strcmpi(g.logfreq, 'off')
from = changeunits(tmptimefreq,imgax,wholeax);
else
from = changeunits([tmptimefreq(1) log(tmptimefreq(2))],imgax,wholeax);
end
to = changeunits([0.5,0.5],topoaxes(n),wholeax);
axes(wholeax);
plot([from(1) to(1)],[from(2) to(2)],LINECOLOR,'linewidth',LINEWIDTH);
hold on
mk=plot(from(1),from(2),[LINECOLOR 'o'],'markersize',9);
set(mk,'markerfacecolor',LINECOLOR);
axis([0 1 0 1]);
axis off;
end
endcaxis = 0;
for n=1:tfpoints
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Plot scalp map using topoplot()
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
axes(topoaxes(n));
scalpmap = squeeze(mean(mean(tfdataori(freqidx{n},timeidx{n},:),1),2));
%topoplot(scalpmap,g.chanlocs,'maplimits',[g.limits(5) g.limits(6)],...
% 'electrodes','on');
if ~isempty(varargin)
topoplot(scalpmap,g.chanlocs,'electrodes','on', selchans_opt{:}, varargin{:});
else
topoplot(scalpmap,g.chanlocs,'electrodes','on', selchans_opt{:});
end
% 'interlimits','electrodes')
axis square;
hold on
if g.timefreqs(n,1) == g.timefreqs(n,2)
tl=title([int2str(g.timefreqs(n,1)),' ms, ',int2str(g.timefreqs(n,3)),' Hz'],'interpreter','none');
else
tl=title([int2str(g.timefreqs(n,1)) '-' int2str(g.timefreqs(n,2)) 'ms, ' ...
int2str(g.timefreqs(n,3)) '-' int2str(g.timefreqs(n,4)) ' Hz'],'interpreter','none');
end
set(tl,'fontsize',AXES_FONTSIZE + 3); %13
endcaxis = max(endcaxis,max(abs(caxis)));
%caxis([g.limits(5:6)]);
end
if strcmpi(g.cmode, 'common')
for n=1:tfpoints
axes(topoaxes(n));
caxis([-endcaxis endcaxis]);
if n==tfpoints && strcmpi(g.cbar, 'on') % && (mod(tfpoints,2)~=0) % image color bar by last map
cb=cbar;
pos = get(cb,'position');
set(cb,'position',[pos(1:2) 0.023 pos(4)]);
end
drawnow
end
end
end
if g.showchan>0 && ~isempty(g.chanlocs)
sbplot(4,4,1,'ax',imgax);
topoplot(g.showchan,g.chanlocs,'electrodes','off', ...
'style', 'blank', 'emarkersize1chan', 10 );
axis('square');
end
if strcmpi(g.axcopy, 'on')
com = 'set(gcf, ''''units'''', ''''pixels''''); tmppos = get(gcf, ''''position''''); set(gcf, ''''position'''', [tmppos(1) tmppos(2) 560 440]); clear tmppos; axis on;';
if strcmpi(g.logfreq, 'native')
% com = [ 'lb = get(gca,''''yticklabel'''');' ...
% 'if iscell(lb) lb = strvcat(lb); end;' ...
% 'ft = str2num(lb);' ...
% 'ft = exp(1).^ft;' ...
% 'ft = unique_bc(round(ft));' ...
% 'ftick = get(gca,''''ytick'''');' ...
% 'ftick = exp(1).^ftick;' ...
% 'ftick = unique_bc(round(log(ftick)));' ...
% 'set(gca,''''ytick'''',ftick);' ...
% 'set(gca,''''yticklabel'''', num2str(ft)); axis on;' ];
axcopy(gcf, com); % turn on axis copying on mouse click
else
axcopy(gca, com); % turn on axis copying on mouse click
end
end
%%%%%%%%%%%%%%%%%%%%%%%% embedded functions %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function tfdat = avedata(tfdat, dim, thresh, mode)
tfsign = sign(mean(tfdat,dim));
tfmask = sum(tfdat ~= 0,dim) >= thresh;
if strcmpi(mode, 'rms')
tfdat = tfmask.*tfsign.*sqrt(mean(tfdat.*tfdat,dim)); % std of all channels
else
tfdat = tfmask.*mean(tfdat,dim); % std of all channels
end
function [tfdatnew, times, freqs] = magnifytwice(tfdat, times, freqs);
indicetimes = [floor(1:0.5:size(tfdat,1)) size(tfdat,1)];
indicefreqs = [floor(1:0.5:size(tfdat,2)) size(tfdat,2)];
tfdatnew = tfdat(indicetimes, indicefreqs, :, :);
times = linspace(times(1), times(end), size(tfdat,2)*2);
freqs = linspace(freqs(1), freqs(end), size(tfdat,1)*2);
% smoothing
gauss2 = gauss2d(3,3);
for S = 1:size(tfdat,4)
for elec = 1:size(tfdat,3)
tfdatnew(:,:,elec,S) = conv2(tfdatnew(:,:,elec,S), gauss2, 'same');
end
end
%tfdatnew = convn(tfdatnew, gauss2, 'same'); % is equivalent to the loop for slowlier
Datasets
Models
