tomwalters@0: % method of class @frame tomwalters@0: % tomwalters@0: % Accummulates the signal between "fromchannel" and "tochannel" tomwalters@0: % This is the famous intervalprofile tomwalters@0: % INPUT VALUES: tomwalters@0: % tomwalters@0: % tomwalters@0: % RETURN VALUE: tomwalters@0: % tomwalters@0: % tomwalters@0: % (c) 2003, University of Cambridge, Medical Research Council tomwalters@0: % Stefan Bleeck (stefan@bleeck.de) tomwalters@0: % http://www.mrc-cbu.cam.ac.uk/cnbh/aimmanual tomwalters@0: % $Date: 2003/07/09 15:24:57 $ tomwalters@0: % $Revision: 1.1 $ tomwalters@0: tomwalters@0: function sumsig=getintervalprofile(fr,fromchannel,tochannel) tomwalters@0: tomwalters@0: tomwalters@0: if nargin<3 tomwalters@0: tochannel=getnrchannels(fr(1)); tomwalters@0: end tomwalters@0: if nargin < 2 tomwalters@0: fromchannel=1; tomwalters@0: end tomwalters@0: tomwalters@0: tomwalters@0: number=max(size(fr)); tomwalters@0: if number==1 tomwalters@0: val=fr.values; tomwalters@0: l=getlength(fr); tomwalters@0: sr=getsr(fr); tomwalters@0: if getnrchannels(fr)>1 tomwalters@0: summe=sum(val(fromchannel:tochannel,:)); tomwalters@0: sumsig=signal(summe); tomwalters@0: sumsig=setsr(sumsig,sr); tomwalters@0: sumsig=setname(sumsig,sprintf('Sum of Frame: %s',getname(fr))); tomwalters@0: sumsig=setstarttime(sumsig,getminimumtime(fr)); tomwalters@0: else tomwalters@0: sumsig=getsinglechannel(fr,1); tomwalters@0: end tomwalters@0: else tomwalters@0: for e=1:number tomwalters@0: tomwalters@0: val=fr(e).values; tomwalters@0: tomwalters@0: l=getlength(fr(e)); tomwalters@0: sr=getsr(fr(e)); tomwalters@0: tomwalters@0: summe(e)=sum(val(fromchannel:tochannel,:)); tomwalters@0: tomwalters@0: sumsig(e)=signal(summe); tomwalters@0: sumsig(e)=setsr(sumsig(e),sr); tomwalters@0: tomwalters@0: sumsig(e)=setstarttime(sumsig(e),getminimumtime(fr(e))) tomwalters@0: tomwalters@0: end tomwalters@0: end