diff nonExposed/timeDomainVisualization.m @ 34:39399de892ef

better handling of displays with pleasantness colors

--- a/nonExposed/timeDomainVisualization.m	Mon May 22 10:46:10 2017 +0200
+++ b/nonExposed/timeDomainVisualization.m	Fri Jun 09 16:37:50 2017 +0200
@@ -20,7 +20,7 @@
 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
@@ -29,6 +29,7 @@
 end
 
 powers(powers<max(powers(:))/500) = 0;
+powers = gaussianSmoothing(powers', 15)';
 
 switch figuresOption
     case 1
@@ -39,19 +40,46 @@
 
 
 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
+
+