tomwalters@0: % demonstration file for 'aim-mat'
tomwalters@0: %
tomwalters@0: % (c) 2006-2008, University of Cambridge, Medical Research Council 
tomwalters@0: % Written by Tom Walters (tcw24@cam.ac.uk)
tomwalters@0: % http://www.pdn.cam.ac.uk/cnbh/aim2006
tomwalters@0: % $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
tomwalters@0: % $Revision: 585 $
tomwalters@0: 
tomwalters@0: function plot_strobes(input, options, current_scale, titlestr)
tomwalters@0: 
tomwalters@0: start_time=options.minimum_time;
tomwalters@0: stop_time=options.maximum_time;
tomwalters@0: 
tomwalters@0: sig=input.data.signal;
tomwalters@0: current_frame=input.data.nap;
tomwalters@0: strobes=input.data.strobes;
tomwalters@0: input.info.current_plot=4;
tomwalters@0: str=get_graphics_options(input,input.info.current_strobes_module);
tomwalters@0: str.minimum_time_interval=start_time;
tomwalters@0: str.maximum_time_interval=stop_time;
tomwalters@0: duration=stop_time-start_time;
tomwalters@0: sig=getpart(sig,start_time,stop_time);
tomwalters@0: nrchan=getnrchannels(current_frame);
tomwalters@0: temp_scale=1;   % for debugging
tomwalters@0: % make the dots of a size of constant ratio to the window size
tomwalters@0: axpos=get(gcf,'Position');
tomwalters@0: marker_dot_size=axpos(3)/168;
tomwalters@0: 
tomwalters@0:   %several channels
tomwalters@0:     hand=plot(current_frame/temp_scale,str);hold on
tomwalters@0:     xlabel('time (ms)');ylabel('Frequency (kHz)');title(titlestr);
tomwalters@0:     colscale=length(hsv)/nrchan;
tomwalters@0:     for i=1:nrchan
tomwalters@0:         if isfield(strobes,'grouped')
tomwalters@0:             nr_sources=size(strobes.cross_strobes{1}.source_cross_channel_value,2);
tomwalters@0:             cols=hsv;
tomwalters@0:             % first plot the originals
tomwalters@0:             %                     herestrobes=strobes.original{i};
tomwalters@0:             %                     nr_here=length(herestrobes.strobes);
tomwalters@0:             %                     col=cols(round(i*colscale),:);
tomwalters@0:             %                     xoffs=-0.001;
tomwalters@0:             %                     for j=1:nr_here
tomwalters@0:             %                         time=herestrobes.strobes(j);
tomwalters@0:             %                         if time>start_time & time<start_time+duration
tomwalters@0:             %                             val=herestrobes.strobe_vals(j);%/current_scale*max(current_frame);
tomwalters@0:             %                             plot3(time2bin(sig,time+xoffs),i,val,'Marker','o','MarkerFaceColor',col,'MarkerEdgeColor',col,'MarkerSize',3);
tomwalters@0:             %                         end
tomwalters@0:             %                     end
tomwalters@0: 
tomwalters@0:             % then plot all the different connections
tomwalters@0:             herestrobesproces=strobes.cross_strobes{i};
tomwalters@0:             herestrobes=herestrobesproces.strobe_times;
tomwalters@0:             col=cols(round(i*colscale),:);
tomwalters@0:             cursize=1;
tomwalters@0:             nr_here=length(herestrobes);
tomwalters@0:             for j=1:nr_here
tomwalters@0:                 count=0;
tomwalters@0:                 for k=1:nr_sources
tomwalters@0:                     target_chan_act=herestrobesproces.source_cross_channel_value(j,k);
tomwalters@0:                     if target_chan_act>1
tomwalters@0:                         colnr=round(k*length(hsv)/nr_sources);
tomwalters@0:                         colnr=min(colnr,length(hsv));
tomwalters@0:                         colnr=max(colnr,1);
tomwalters@0:                         col=cols(colnr,:);
tomwalters@0:                         cursize=target_chan_act/100;
tomwalters@0:                     else
tomwalters@0:                         continue
tomwalters@0:                     end
tomwalters@0:                     time=herestrobes(j);
tomwalters@0:                     if time>start_time & time<start_time+duration
tomwalters@0:                         offx=(stop_time-start_time)/80*count;
tomwalters@0:                         count=count+1;
tomwalters@0:                         val=herestrobesproces.strobe_vals(j)/temp_scale;%/current_scale*max(current_frame);
tomwalters@0:                         if strcmp(handles.screen_modus,'paper')
tomwalters@0:                             p=plot3(time2bin(sig,time+offx),i,val,'Marker','o','MarkerSize',4,'MarkerFaceColor','k','MarkerEdgeColor','k','LineWidth',1);
tomwalters@0:                         else
tomwalters@0:                             plot3(time2bin(sig,time+offx),i,val,'Marker','o','MarkerFaceColor','r','MarkerEdgeColor','w','MarkerSize',cursize);
tomwalters@0:                         end
tomwalters@0:                     end
tomwalters@0:                 end
tomwalters@0:             end
tomwalters@0:         else    % only one set of strobes
tomwalters@0:             herestrobes=strobes{i};
tomwalters@0:             nr_here=length(herestrobes.strobes);
tomwalters@0:             %                     marker_dot_size=4;
tomwalters@0:             for j=1:nr_here
tomwalters@0:                 time=herestrobes.strobes(j);
tomwalters@0:                 if time>start_time & time<stop_time
tomwalters@0:                     % 						if j<=length(herestrobes.strobe_vals)
tomwalters@0:                     val=herestrobes.strobe_vals(j);%/current_scale*max(current_frame);
tomwalters@0:                     %if strcmp(handles.screen_modus,'paper')
tomwalters@0:                     %    p=plot3(time2bin(sig,time),i,val/temp_scale,'Marker','o','MarkerSize',4,'MarkerFaceColor','k','MarkerEdgeColor','k','LineWidth',1);
tomwalters@0:                     %else
tomwalters@0:                         plot3(time2bin(sig,time),i,val,'Marker','o','MarkerSize',marker_dot_size,'MarkerFaceColor','r','MarkerEdgeColor','r','LineWidth',1);
tomwalters@0:                     %end
tomwalters@0:                     % 						end
tomwalters@0:                 end
tomwalters@0:             end
tomwalters@0:         end % one or two strobe sources
tomwalters@0:     end
tomwalters@0:     zmin=0;zmax=max(current_frame)/current_scale;
tomwalters@0: 
tomwalters@0:     % TCW AIM 2006
tomwalters@0: %    if hand_scaling == 1
tomwalters@0:         set(gca,'Zlim',[zmin,1/current_scale]);
tomwalters@0: %    else
tomwalters@0: %        set(gca,'Zlim',[zmin,zmax]);
tomwalters@0: %    end
tomwalters@0: 
tomwalters@0:     %set(gca,'Zlim',[zmin,zmax]);
tomwalters@0:  %   if strcmp(handles.screen_modus,'paper')
tomwalters@0:  %       par=get(hand,'parent');
tomwalters@0:  %       set(par,'FontSize',12);xlab=get(par,'xlabel');set(xlab,'FontSize',12);ylab=get(par,'ylabel');set(ylab,'FontSize',12);
tomwalters@0:  %   end
tomwalters@0: 
tomwalters@0: hold off
tomwalters@0: