Daniel@0: function r = mirfeatures(x,varargin)
Daniel@0: %   f = mirfeatures(x) computes a large set of features from one or several
Daniel@0: %       audio files. x can be either the name of an audio file, or the
Daniel@0: %       'Folder' keyword.
Daniel@0: %   mirfeatures(...,'Stat') returns the statistics of the features instead
Daniel@0: %       of the complete features themselves.
Daniel@0: %   mirfeatures(...,'Segment',t) segments the audio sequence at the 
Daniel@0: %       temporal positions indicated in the array t (in s.), and analyzes
Daniel@0: %       each segment separately.
Daniel@0: 
Daniel@0: %(not available yet)
Daniel@0: %   mirfeatures(...,'Filterbank',nc) computes the analysis on each channel
Daniel@0: %       of a filterbank decomposition.
Daniel@0: %       Default value: nc = 5
Daniel@0: %   mirfeatures(...,'Frame',...) 
Daniel@0: %   mirfeatures(...,'Normal')
Daniel@0: %   mirfeatures(...,'Sampling',s)
Daniel@0: %   miraudio options (Extract, ...)
Daniel@0: 
Daniel@0: [stat,nchan,segm,feat] = scanargin(varargin);
Daniel@0: 
Daniel@0: if isa(x,'miraudio') || isa(x,'mirdesign')
Daniel@0:     a = miraudio(x,'Normal'); % normalize with respect to RMS energy 
Daniel@0:                               % in order to consider timbre independently of
Daniel@0:                              % energy
Daniel@0: else
Daniel@0:     a = miraudio('Design','Normal');
Daniel@0: end
Daniel@0: 
Daniel@0: if not(isempty(segm))
Daniel@0:     a = mirsegment(a,segm);
Daniel@0: end
Daniel@0: 
Daniel@0: 
Daniel@0: 
Daniel@0: % DYNAMICS
Daniel@0: % --------
Daniel@0: 
Daniel@0: r.dynamics.rms = mirrms(a,'Frame');
Daniel@0: % Perceived dynamics: spectral slope?
Daniel@0: 
Daniel@0: % RHYTHM
Daniel@0: % ------
Daniel@0: 
Daniel@0: r.fluctuation = mirstruct;
Daniel@0: r.fluctuation.tmp.f = mirfluctuation(a,'Summary');
Daniel@0: r.fluctuation.peak = mirpeaks(r.fluctuation.tmp.f,'Total',1);%only one?
Daniel@0: r.fluctuation.centroid = mircentroid(r.fluctuation.tmp.f);
Daniel@0: 
Daniel@0: r.rhythm = mirstruct;
Daniel@0: r.rhythm.tmp.onsets = mironsets(a);
Daniel@0: 
Daniel@0: %r.rhythm.eventdensity = ...
Daniel@0: 
Daniel@0: r.rhythm.tempo = mirtempo(r.rhythm.tmp.onsets,'Frame');
Daniel@0: %r.rhythm.pulseclarity = mirpulseclarity(r.tmp.onsets,'Frame');
Daniel@0:     % Should use the second output of mirtempo.
Daniel@0: 
Daniel@0: attacks = mironsets(r.rhythm.tmp.onsets,'Attacks');
Daniel@0: r.rhythm.attack.time = mirattacktime(attacks);
Daniel@0: r.rhythm.attack.slope = mirattackslope(attacks);
Daniel@0: 
Daniel@0: % TIMBRE
Daniel@0: % ------
Daniel@0: 
Daniel@0: f = mirframe(a,.05,.5);
Daniel@0: r.spectral = mirstruct;
Daniel@0: r.spectral.tmp.s = mirspectrum(f);
Daniel@0: %pitch = mirpitch(a,'Frame',.05,.5);
Daniel@0: 
Daniel@0: r.spectral.centroid = mircentroid(r.spectral.tmp.s);
Daniel@0: r.spectral.brightness = mirbrightness(r.spectral.tmp.s);
Daniel@0: r.spectral.spread = mirspread(r.spectral.tmp.s);
Daniel@0: r.spectral.skewness = mirskewness(r.spectral.tmp.s);
Daniel@0: r.spectral.kurtosis = mirkurtosis(r.spectral.tmp.s);
Daniel@0: r.spectral.rolloff95 = mirrolloff(r.spectral.tmp.s,95);
Daniel@0: r.spectral.rolloff85 = mirrolloff(r.spectral.tmp.s,85);
Daniel@0: r.spectral.spectentropy = mirentropy(r.spectral.tmp.s);
Daniel@0: r.spectral.flatness = mirflatness(r.spectral.tmp.s);
Daniel@0: 
Daniel@0: r.spectral.roughness = mirroughness(r.spectral.tmp.s);
Daniel@0: r.spectral.irregularity = mirregularity(r.spectral.tmp.s);
Daniel@0: %r.spectral.inharmonicity = mirinharmonicity(r.spectral.tmp.s,'f0',pitch);
Daniel@0: 
Daniel@0: r.spectral.mfcc = mirmfcc(r.spectral.tmp.s);
Daniel@0: r.spectral.dmfcc = mirmfcc(r.spectral.mfcc,'Delta');
Daniel@0: r.spectral.ddmfcc = mirmfcc(r.spectral.dmfcc,'Delta');
Daniel@0: 
Daniel@0: r.timbre.zerocross = mirzerocross(f);
Daniel@0: r.timbre.lowenergy = mirlowenergy(f);
Daniel@0: r.timbre.spectralflux = mirflux(f);
Daniel@0: 
Daniel@0: % PITCH
Daniel@0: % -----
Daniel@0: 
Daniel@0: r.tonal = mirstruct;
Daniel@0: r.tonal.tmp.chromagram = mirchromagram(a,'Frame','Wrap',0,'Pitch',0);
Daniel@0: r.tonal.chromagram.peak=mirpeaks(r.tonal.tmp.chromagram,'Total',1);
Daniel@0: r.tonal.chromagram.centroid=mircentroid(r.tonal.tmp.chromagram);
Daniel@0: 
Daniel@0: % TONALITY/HARMONY
Daniel@0: % ----------------
Daniel@0: 
Daniel@0: keystrengths = mirkeystrength(r.tonal.tmp.chromagram);
Daniel@0: [k r.tonal.keyclarity] = mirkey(keystrengths,'Total',1);
Daniel@0: %r.tonal.keyclarity = k{2};
Daniel@0: r.tonal.mode = mirmode(keystrengths);
Daniel@0: r.tonal.hcdf = mirhcdf(r.tonal.tmp.chromagram);
Daniel@0: 
Daniel@0: if stat
Daniel@0:     r = mirstat(r);
Daniel@0:     % SHOULD COMPUTE STAT OF CURVES FROM FRAMED_DECOMPOSED HIGH FEATURES
Daniel@0: end
Daniel@0:     
Daniel@0: if not(isa(x,'miraudio')) && not(isa(x,'mirdesign'))
Daniel@0:     r = mireval(r,x);
Daniel@0: end
Daniel@0: 
Daniel@0: 
Daniel@0: function [stat,nchan,segm,feat] = scanargin(v)
Daniel@0: stat = 0;
Daniel@0: nchan = 1;
Daniel@0: segm = [];
Daniel@0: feat = {};
Daniel@0: i = 1;
Daniel@0: while i <= length(v)
Daniel@0:     arg = v{i};
Daniel@0:     if ischar(arg) && strcmpi(arg,'Filterbank')
Daniel@0:         i = i+1;
Daniel@0:         if i <= length(v)
Daniel@0:             nchan = v{i};
Daniel@0:         else
Daniel@0:             nchan = 10;
Daniel@0:         end
Daniel@0:     elseif ischar(arg) && strcmpi(arg,'Stat')
Daniel@0:         i = i+1;
Daniel@0:         if i <= length(v)
Daniel@0:             stat = v{i};
Daniel@0:         else
Daniel@0:             stat = 1;
Daniel@0:         end
Daniel@0:     elseif ischar(arg) && strcmpi(arg,'Segment')
Daniel@0:         i = i+1;
Daniel@0:         if i <= length(v)
Daniel@0:             segm = v{i};
Daniel@0:         else
Daniel@0:             segm = 1;
Daniel@0:         end
Daniel@0:     else
Daniel@0:         feat{end+1} = arg;
Daniel@0:     end    
Daniel@0:     i = i+1;
Daniel@0: end