function PlotSgnMrkNc(sgn,pqrs,freq, name, nband, nsecr, normModef, swYTickLabel)
% --------------------------------------------------------------------------------------------
% PlotSgnMrkNc.m: Plot signals (on multiple horizontal strips) with vertical markers.
% --> With signal normalization
% sgn = signal matrix
% pqrs = cell array (or matrix) of time indexes (markers position)
% freq = sampling frequency
% name = figure title
% nbands = number of strips per page
% nsecr = number of seconds per row
% normMode = normalization mode
% (0 -> on all signals,
% 1 -> for each channel
% 2 -> on each page
% 3 -> on each row)
% swYTickLabel = 0 -> no yticklabel; 1 -> yticklabel
%
% Version allowing marker positions in cell array
% --------------------------------------------------------------------------------------------
% Maurizio Varanini, Clinical Physiology Institute, CNR, Pisa, Italy
% For any comment or bug report, please send e-mail to: maurizio.varanini@ifc.cnr.it
% --------------------------------------------------------------------------------------------
% convert column to row if necessary
[m,n] = size(sgn); if n<m, sgn = sgn.'; end
if(iscell(pqrs))
[m,n] = size(pqrs); if n>m, pqrs = pqrs.'; end
for ia=1:length(pqrs)
[m,n] = size(pqrs{ia}); if n<m, pqrs{ia} = pqrs{ia}.'; end
end
else
[m,n] = size(pqrs); if n<m, pqrs = pqrs.'; end
end
nsgn=size(sgn,1);
nsampt=size(sgn,2);
if(nargin<3 || isempty(freq)), freq=1; end
nsect=nsampt/freq;
if(nargin<4 || isempty(name)), name=''; end
if(nargin<5 || isempty(nband)), nband=6; end
if(nargin<6 || isempty(nsecr)), nsecr=nsect/nband; end
if(nargin<7 || isempty(normModef)), normModef=1; end
if(nargin<8), swYTickLabel=1; end
normMode= floor(normModef);
ymargc = normModef-normMode;
ymint=min(min(sgn));
ymaxt=max(max(sgn));
nsampr=nsecr*freq; % number of samples per row
isampi=1;
isampf=0;
isampf_i=floor(isampf);
vmrkcol=['m', 'b', 'r', 'g', 'k'];
if(normMode==0)
yminc=ymint*ones(1,nsgn);
ymaxc=ymaxt*ones(1,nsgn);
ymaxc= ymaxc + (ymaxc==yminc);
yrangec=ymaxc-yminc;
yminc=yminc-ymargc*yrangec;
ymaxc=ymaxc+ymargc*yrangec;
elseif(normMode==1)
yminc=min(sgn,[],2);
ymaxc=max(sgn,[],2);
ymaxc= ymaxc + (ymaxc==yminc);
yrangec=ymaxc-yminc;
yminc=yminc-ymargc*yrangec;
ymaxc=ymaxc+ymargc*yrangec;
end
while isampf_i < nsampt
figure, set(gcf,'Color','white');
if(normMode==2)
yminc=min(sgn(:,isampi:isampf_i),[],2);
ymaxc=max(sgn(:,isampi:isampf_i),[],2);
ymaxc= ymaxc + (ymaxc==yminc);
yrangec=ymaxc-yminc;
yminc=yminc-ymargc*yrangec;
ymaxc=ymaxc+ymargc*yrangec;
end
for i=0:nband-1
subplot(nband,1,i+1);
isampf = isampi+nsampr-1;
isampf_i=ceil(isampf);
if(isampf_i > nsampt), isampf_i=nsampt; end;
if(normMode>=3)
yminc=min(sgn(:,isampi:isampf_i),[],2);
ymaxc=max(sgn(:,isampi:isampf_i),[],2);
ymaxc= ymaxc + (ymaxc==yminc);
yrangec=ymaxc-yminc;
yminc=yminc-ymargc*yrangec;
ymaxc=ymaxc+ymargc*yrangec;
end
ymin=-1; ymax=1;
yn=(ymax-ymin)./(ymaxc-yminc);
yoffset=(1*(ymax-ymin)).*((nsgn:-1:1)-1)';
clear ysgn;
ysgn=zeros(nsgn,isampf_i-isampi+1);
for is=1:nsgn
ysgn(is,:)= yoffset(is) +ymin + yn(is)*(sgn(is,isampi:isampf_i)-yminc(is));
end
yming=ymin;
ymaxg=ymax+max(yoffset);
xmin=isampi/freq;
xmax=(isampi+nsampr)/freq;
plot([isampi:isampf_i]/freq ,ysgn);
if(~swYTickLabel), set(gca,'YTickLabel',[]); end
axis tight;
if(i==0), ht=title(name); set(ht,'Interpreter','none'); end;
xlim([xmin, xmax]);
ylim([yming, ymaxg]);
if(~isempty(pqrs))
hold on;
for im=1:size(pqrs,1); % loop on marker type
if(iscell(pqrs))
vmrk=pqrs{im};
else
vmrk=pqrs(im,:);
end
mrkcol=vmrkcol(1+mod(im-1,size(vmrkcol,2)));
iqrsint=find((isampi < vmrk)&(isampf_i > vmrk));
if(~isempty(iqrsint))
pqrsint=vmrk(iqrsint);
fxpqrs=[pqrsint; pqrsint]/freq;
fypqrs = ones(size(fxpqrs)); fypqrs(1,:)=yming; fypqrs(2,:)=ymaxg;
plot(fxpqrs,fypqrs, [mrkcol ':+']);
end
end
end
if(isampf_i == nsampt), break; end;
isampi = isampf_i+1;
end
end
end %== function ================================================================
%