gregoire@16: function status = timeDomainVisualization(tracks,figNum,settingFigure,figuresOption,fileName)
mathieu@14: 
mathieu@14: % This program was written by Mathias Rossignol & Grégoire Lafay
mathieu@14: % is Copyright (C) 2015 IRCAM <http://www.ircam.fr>
mathieu@14: %
mathieu@14: % This program is free software: you can redistribute it and/or modify it
mathieu@14: % under the terms of the GNU General Public License as published by the Free
mathieu@14: % Software Foundation, either version 3 of the License, or (at your option)
mathieu@14: % any later version.
mathieu@14: %
mathieu@14: % This program is distributed in the hope that it will be useful, but
mathieu@14: % WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
mathieu@14: % or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
mathieu@14: % for more details.
mathieu@14: %
mathieu@14: % You should have received a copy of the GNU General Public License along
mathieu@14: % with this program.  If not, see <http://www.gnu.org/licenses/>.
mathieu@14: 
mathieu@14: % Producing a time domain visualization
gregoire@16: step=0.5; % seconds
gregoire@16: wStep = round(step*settingFigure.sr);
gregoire@16: wSize = wStep*2;
mathieu@34: for t=1:size(tracks,1)
gregoire@16:     for i=0:floor((size(tracks,2)-wSize)/wStep)
mathieu@14:         % Not really power, but more nicely additive, better suited for
mathieu@14:         % this representation I think
gregoire@16:         powers(i+1,t) = norm(tracks(t,i*wStep+1:min(size(tracks,2), i*wStep+wSize)));
mathieu@14:     end
mathieu@14: end
mathieu@14: 
mathieu@14: powers(powers<max(powers(:))/500) = 0;
mathieu@34: powers = gaussianSmoothing(powers', 15)';
mathieu@14: 
gregoire@16: switch figuresOption
gregoire@16:     case 1
gregoire@16:        f=figure('Visible', 'off');
gregoire@16:     case 2
gregoire@16:        f=figure(figNum);
gregoire@16: end
mathieu@14: 
mathieu@14: 
mathieu@14: clf;
mathieu@34: h = area(powers,'LineWidth', 2, 'EdgeColor', [.2, .2, .2]);
mathieu@34: for t=1:size(tracks,1)
mathieu@34:     h(t).FaceColor = settingFigure.cmap(t, :);
mathieu@34: end
mathieu@34: 
gregoire@16: xlim([0 size(powers,1)])
gregoire@16: xtick=0:round(20/step):size(powers,1); % every 20 sec
gregoire@16: set(gca,'YTick', [],'YTicklabel', [],'xtick',xtick,'xticklabel',xtick*step);
gregoire@16: xlabel('time (sec)')
mathieu@34: % ylabel('amplitude')
mathieu@14: 
gregoire@16: set(findall(f,'-property','FontSize'),'FontSize',settingFigure.FontSize)
gregoire@16: set(findall(f,'-property','FontName'),'FontName','Arial')
mathieu@14: 
mathieu@14: 
mathieu@34: box off %// remove outer border
mathieu@34: set(gca,'ycolor', [1 1 1]);
mathieu@34: ax1 = gca;
mathieu@34: yruler = ax1.YRuler;
mathieu@34: yruler.Axle.Visible = 'off';
mathieu@34: 
mathieu@34: if figuresOption == 1
mathieu@34:     set(f,'PaperUnits','centimeters')
mathieu@34:     set(f,'PaperPositionMode','manual')
mathieu@34:     set(f,'papersize',[settingFigure.width,settingFigure.height])
mathieu@34:     set(f,'paperposition',[0,0,settingFigure.width,settingFigure.height])
mathieu@34:     set(gca,'ycolor', [1 1 1]);
mathieu@34:     print(f,fileName,'-dpng')
mathieu@14: end
mathieu@34: 
mathieu@34: 
mathieu@34: function as = gaussianSmoothing(as, factor)
mathieu@34: 
mathieu@34: if ~exist('factor', 'var'), factor = 200; end
mathieu@34: 
mathieu@34: % gaussian filtering for smooth display
mathieu@34: g = gausswin(ceil(size(as, 2)/factor)); % shall adapt the size to the length of the file
mathieu@34: g = g/sum(g);
mathieu@34: for k=1:size(as, 1)
mathieu@34:     as(k, :) = conv(as(k, :), g, 'same');
mathieu@34: end
mathieu@34: 
mathieu@34: