Mauch MIREX 2010

function song = song_syncchroma_special(song, param)

[song.syncchroma.treble, treblenrm, beats] = syncchroma(song.chroma.treble, song.beattimes, song.t);

[song.syncchroma.bass, bassnrm] = syncchroma(song.chroma.bass, song.beattimes, song.t(:));

song.syncchroma.bass(isnan(sum(song.syncchroma.bass,2)),:) = 1;

song.syncchroma.treble(isnan(sum(song.syncchroma.treble,2)),:) = 1;

song.syncchroma.bass(sum(song.syncchroma.bass,2)<.01,:) = 1;

song.syncchroma.treble(sum(song.syncchroma.treble,2)<.01,:) = 1;

bcnrm = qnorm(song.syncchroma.bass,1,2);

if param.dsp.bassbin13 % if the model requires a thirteenth bass bin

    bc = qnormalise(song.syncchroma.bass,1,2);

    bc13 = (12 * max(bc,[],2)).^(-2) .* bcnrm;

    bc13(isnan(bc13)) = 1;

    song.syncchroma.bass = [song.syncchroma.bass, bc13];

end

switch param.dsp.norm

    case 'max'

        song.syncchroma.bass = qnormalise(song.syncchroma.bass,inf,2);

        song.syncchroma.treble = qnormalise(song.syncchroma.treble,inf,2);

    case 'L1'

        song.syncchroma.bass = qnormalise(song.syncchroma.bass,1,2);

        song.syncchroma.treble = qnormalise(song.syncchroma.treble,1,2);

    case 'L2'

        song.syncchroma.bass = qnormalise(song.syncchroma.bass,2,2);

        song.syncchroma.treble = qnormalise(song.syncchroma.treble,2,2);

    case 'standardize'

        song.syncchroma.bass = standardize(song.syncchroma.bass);

        song.syncchroma.treble = standardize(song.syncchroma.treble);

    case 'meta'

        song.syncchroma.bass = qnormalise(song.syncchroma.bass,inf,2);

        disp('no standardization for treble chroma; max normalisation for bass')

    otherwise

        'blub'

end

% make all NaNs to 1

song.syncchroma.bass(isnan(sum(song.syncchroma.bass,2)),:) = 1;

song.syncchroma.treble(isnan(sum(song.syncchroma.treble,2)),:) = 1;

% aah, what does this do???

% song.syncchroma.bass(sum(song.syncchroma.bass,2) < 1,:) = 1;

% song.syncchroma.treble(sum(song.syncchroma.treble,2) < 1,:) = 1;

% below "wide" chroma processing, that's neglected ...

try

    [song.syncchroma.wide, widenrm] = syncchroma(song.chroma.wide, song.beattimes, song.t);

    for iBeat = 1:song.nBeat

        song.syncchroma.wide(iBeat,:) = song.syncchroma.wide(iBeat,:) / widenrm(iBeat);

    end

catch

    warning('could not find wide chroma');

end

function [bc, nrm, beats] = syncchroma(chromag,beats,frame_t)

% assigning beats to frames

nBeat = size(beats,1);

nFrame = length(frame_t);

bc = ones(nBeat,size(chromag,1));

chromag(:,isnan(sum(chromag))) = 1;

% first beat

currind = find(frame_t >= beats(1,1) & frame_t < beats(1,2));

bc(1,:) = median(chromag(:,currind),2)';

for iBeat = 2:nBeat

%     prestarttimes = [beats(iBeat - 1,1), beats(iBeat,1), beats(iBeat,2)];

%     starttimes = (prestarttimes + beats(iBeat,1))/2;

    starttimes = beats(iBeat,1) + [-.1 -.05 0 .05 .1];

    entro = 0;

    for iST = 1:5

        currind1 = find(frame_t >= beats(iBeat - 1,1) & frame_t < starttimes(iST));

        currind2 = find(frame_t >= starttimes(iST) & frame_t < beats(iBeat,2));

        vec = [ mean(chromag(:,currind2),2)'];

        vec = vec/sum(vec);

        vec = vec(vec>0);

        ent = -vec * log2(vec');

        if  ent > entro 

            entro = ent;

            beats(iBeat - 1, 2) = starttimes(iST);

            beats(iBeat, 1) = starttimes(iST);

        end

    end

    currind = find(frame_t >= beats(iBeat,1) & frame_t < beats(iBeat,2));

    bc(iBeat,:) = median(chromag(:,currind),2)';

end

nrm = qnorm(bc,inf,2);

% bc = qnormalise(bc,inf,2);