tomwalters@0: % method of class @signal
tomwalters@0: % 
tomwalters@0: %   INPUT VALUES:
tomwalters@0: %  
tomwalters@0: %   RETURN VALUE:
tomwalters@0: %
bleeck@3: % This external file is included as part of the 'aim-mat' distribution package
bleeck@3: % (c) 2011, University of Southampton
bleeck@3: % Maintained by Stefan Bleeck (bleeck@gmail.com)
bleeck@3: % download of current version is on the soundsoftware site: 
bleeck@3: % http://code.soundsoftware.ac.uk/projects/aimmat
bleeck@3: % documentation and everything is on http://www.acousticscale.org
tomwalters@0: 
tomwalters@0: function rethandle=fill(sig,border,stil)
tomwalters@0: % usage: handle=fill(sig,border,stil)
tomwalters@0: % works as plot, but fills the content with the color of the signal
tomwalters@0: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: if nargin<3
tomwalters@0: 	stil='b';
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: % if we are dealing with milliseconds (very often), then the units are in
tomwalters@0: % ms as well
tomwalters@0: if strfind(sig.unit_x,'(ms)')>0
tomwalters@0:     multiplier=1000;
tomwalters@0: else
tomwalters@0:     multiplier=1;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: 
tomwalters@0: sr=getsr(sig);
tomwalters@0: durationbin=size(sig.werte,1);
tomwalters@0: duration=durationbin/sr;
tomwalters@0: if nargin<2
tomwalters@0: 	t1=1;
tomwalters@0: 	t2=durationbin;
tomwalters@0: 	border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isstruct(border) || isempty(border)
tomwalters@0: 	options=border;
tomwalters@0: 	t1=1;
tomwalters@0: 	t2=durationbin;
tomwalters@0: 	border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0: else
tomwalters@0: 	options=[];
tomwalters@0: 	if ischar(border)
tomwalters@0: 		stil=border;
tomwalters@0: 		t1=1;
tomwalters@0: 		t2=durationbin;
tomwalters@0: 		border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0: 	else
tomwalters@0: 		if nargin <2
tomwalters@0: 			t1=1;
tomwalters@0: 			t2=durationbin;
tomwalters@0: 			border=[t1 t2 min(sig.werte(t1:t2)) max(sig.werte(t1:t2))];
tomwalters@0: 		else
tomwalters@0: 			t1=border(1);
tomwalters@0: 			t2=border(2);
tomwalters@0: 			x1=time2bin(sig,t1);
tomwalters@0: 			x2=time2bin(sig,t2);
tomwalters@0: 			nr=size(border,2);
tomwalters@0: 			if nr==2    % wenn nur die x-Werte angegeben werden
tomwalters@0: 				border=[x1 x2 min(sig.werte(x1+1:x2)) max(sig.werte(x1+1:x2))];
tomwalters@0: 			else
tomwalters@0: 				border(1)=x1;
tomwalters@0: 				border(2)=x2;
tomwalters@0: 			end
tomwalters@0: 		end
tomwalters@0: 	end
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if ~isfield(options,'is_log');
tomwalters@0: 	is_log=0;
tomwalters@0: else
tomwalters@0: 	is_log=options.is_log;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isfield(options,'minimum_time');
tomwalters@0: 	border(1)=time2bin(sig,options.minimum_time);
tomwalters@0: 	minimum_time=options.minimum_time;
tomwalters@0: else
tomwalters@0: 	minimum_time=bin2time(sig,border(1));
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if isfield(options,'maximum_time');
tomwalters@0: 	border(2)=time2bin(sig,options.maximum_time);
tomwalters@0: 	maximum_time=options.maximum_time;
tomwalters@0: else
tomwalters@0: 	maximum_time=bin2time(sig,border(2));
tomwalters@0: end
tomwalters@0: 
tomwalters@0: 
tomwalters@0: % if the time scale is reversed (time from left to right)
tomwalters@0: if isfield(options,'time_reversed');
tomwalters@0: 	time_reversed=options.time_reversed;
tomwalters@0: else
tomwalters@0: 	time_reversed=0;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: start_time=getminimumtime(sig);
tomwalters@0: 
tomwalters@0: min_x_screen=border(1);   % einer wird abgezogen damit wir bei Null beginnen, nicht beim ersten bin, was blöde aussieht
tomwalters@0: max_x_screen=border(2);
tomwalters@0: minshowy=border(3);
tomwalters@0: maxy=border(4);
tomwalters@0: 
tomwalters@0: xvals=getxvalues(sig);
tomwalters@0: 
tomwalters@0: 
tomwalters@0: xvals=xvals.*multiplier;
tomwalters@0: yvals=getvalues(sig);
tomwalters@0: 
tomwalters@0: 
tomwalters@0: % fill up the corners
tomwalters@0: yvals=[0 yvals' 0];
tomwalters@0: tstart=getminimumtime(sig);
tomwalters@0: tstop=getmaximumtime(sig);
tomwalters@0: % xvals=[tstart:1/sr*multiplier:tstop*multiplier+1/sr*multiplier];
tomwalters@0: xvals=[tstart*multiplier:1/sr*multiplier:tstop*multiplier+1/sr*multiplier];
tomwalters@0: % xvals=1:length(yvals);
tomwalters@0: 
tomwalters@0: % this is the plotting command:
tomwalters@0: handle=fill(xvals,yvals,stil,'facecolor',stil,'EdgeColor',stil);
tomwalters@0: 
tomwalters@0: if time_reversed
tomwalters@0:     set(gca,'XDir','reverse')   % turn them around, because the higher values shell end on the right
tomwalters@0: else
tomwalters@0:     set(gca,'XDir','normal')   % normale ausrichtung
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if is_log
tomwalters@0:     set(gca,'XScale','log')
tomwalters@0:     t=minimum_time*multiplier;
tomwalters@0:     tix=[t 2*t 4*t 8*t 16*t 32*t 64*t 128*t 256*t 512*t 1024*t];
tomwalters@0:     if ~isempty(sig.x_tick_labels)
tomwalters@0:         ti=sig.x_tick_labels;
tomwalters@0:         set(gca,'XTicklabel',ti);
tomwalters@0:     end
tomwalters@0:     set(gca,'XTick',tix);
tomwalters@0: else % its linear
tomwalters@0:     set(gca,'XScale','linear')
tomwalters@0: end
tomwalters@0: 
tomwalters@0: miny=border(3);
tomwalters@0: maxy=border(4);
tomwalters@0: y=[miny maxy];
tomwalters@0: if miny==maxy
tomwalters@0: 	maxy=miny+1;
tomwalters@0: 	miny=miny-1;
tomwalters@0: end
tomwalters@0: if isnan(miny) || isnan(maxy)
tomwalters@0: else
tomwalters@0:     axis([minimum_time*multiplier maximum_time*multiplier miny*1.05 maxy*1.05]);
tomwalters@0: end
tomwalters@0: 
tomwalters@0: xlabel(sig.unit_x);
tomwalters@0: ylabel(sig.unit_y);
tomwalters@0: title(sig.name,'Interpreter','none');
tomwalters@0: 
tomwalters@0: if nargout==1
tomwalters@0: 	rethandle=handle;
tomwalters@0: end
tomwalters@0: 
tomwalters@0: return