comparison nonExposed/timeDomainVisualization.m @ 34:39399de892ef

better handling of displays with pleasantness colors

 
 % Producing a time domain visualization
 step=0.5; % seconds
 wStep = round(step*settingFigure.sr);
 wSize = wStep*2;
-for t=1:size(tracks,1);
+for t=1:size(tracks,1)
     for i=0:floor((size(tracks,2)-wSize)/wStep)
         % Not really power, but more nicely additive, better suited for
         % this representation I think
         powers(i+1,t) = norm(tracks(t,i*wStep+1:min(size(tracks,2), i*wStep+wSize)));
     end
 end
 
 powers(powers<max(powers(:))/500) = 0;
+powers = gaussianSmoothing(powers', 15)';
 
 switch figuresOption
     case 1
        f=figure('Visible', 'off');
     case 2
        f=figure(figNum);
 end
 
 
 clf;
-area(powers,'LineWidth', 1/1.6, 'EdgeColor', [.3, .3, .3]);
+h = area(powers,'LineWidth', 2, 'EdgeColor', [.2, .2, .2]);
+for t=1:size(tracks,1)
+    h(t).FaceColor = settingFigure.cmap(t, :);
+end
+
 xlim([0 size(powers,1)])
 xtick=0:round(20/step):size(powers,1); % every 20 sec
 set(gca,'YTick', [],'YTicklabel', [],'xtick',xtick,'xticklabel',xtick*step);
 xlabel('time (sec)')
+% ylabel('amplitude')
 
-set(f,'PaperUnits','centimeters')
-set(f,'PaperPositionMode','manual')
-set(f,'papersize',[settingFigure.width,settingFigure.height])
-set(f,'paperposition',[0,0,settingFigure.width,settingFigure.height])
 set(findall(f,'-property','FontSize'),'FontSize',settingFigure.FontSize)
 set(findall(f,'-property','FontName'),'FontName','Arial')
 
-print(f,fileName,'-dpng')
 
+box off %// remove outer border
+set(gca,'ycolor', [1 1 1]);
+ax1 = gca;
+yruler = ax1.YRuler;
+yruler.Axle.Visible = 'off';
+
+if figuresOption == 1
+    set(f,'PaperUnits','centimeters')
+    set(f,'PaperPositionMode','manual')
+    set(f,'papersize',[settingFigure.width,settingFigure.height])
+    set(f,'paperposition',[0,0,settingFigure.width,settingFigure.height])
+    set(gca,'ycolor', [1 1 1]);
+    print(f,fileName,'-dpng')
 end
+
+
+function as = gaussianSmoothing(as, factor)
+
+if ~exist('factor', 'var'), factor = 200; end
+
+% gaussian filtering for smooth display
+g = gausswin(ceil(size(as, 2)/factor)); % shall adapt the size to the length of the file
+g = g/sum(g);
+for k=1:size(as, 1)
+    as(k, :) = conv(as(k, :), g, 'same');
+end
+
+