annotate AgentList.h @ 11:7169eb47b1bc

formatting only
author Chris Cannam
date Fri, 30 Sep 2011 23:11:14 +0100
parents 1c1e98cd1b2e
children 59520cd6abac
rev   line source
Chris@8 1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
Chris@8 2
Chris@8 3 /*
Chris@8 4 Vamp feature extraction plugin for the BeatRoot beat tracker.
Chris@8 5
Chris@8 6 Centre for Digital Music, Queen Mary, University of London.
Chris@8 7 This file copyright 2011 Simon Dixon, Chris Cannam and QMUL.
Chris@8 8
Chris@8 9 This program is free software; you can redistribute it and/or
Chris@8 10 modify it under the terms of the GNU General Public License as
Chris@8 11 published by the Free Software Foundation; either version 2 of the
Chris@8 12 License, or (at your option) any later version. See the file
Chris@8 13 COPYING included with this distribution for more information.
Chris@8 14 */
Chris@8 15
Chris@8 16 #ifndef _AGENT_LIST_H_
Chris@8 17 #define _AGENT_LIST_H_
Chris@8 18
Chris@8 19 #include "Agent.h"
Chris@8 20 #include "Event.h"
Chris@8 21
Chris@8 22 #include <algorithm>
Chris@8 23
Chris@8 24 /** Class for maintaining the set of all Agents involved in beat tracking a piece of music.
Chris@8 25 */
Chris@8 26 class AgentList
Chris@8 27 {
Chris@8 28 public:
Chris@8 29 typedef std::vector<Agent> Container;
Chris@8 30 typedef Container::iterator iterator;
Chris@8 31
Chris@8 32 protected:
Chris@8 33 Container list;
Chris@8 34
Chris@8 35 public:
Chris@8 36 // expose some vector methods
Chris@8 37 //!!! can we remove these again once the rest of AgentList is implemented?
Chris@8 38 bool empty() const { return list.empty(); }
Chris@8 39 Container::iterator begin() { return list.begin(); }
Chris@8 40 Container::iterator end() { return list.end(); }
Chris@8 41 void push_back(const Agent &a) { list.push_back(a); }
Chris@8 42
Chris@11 43 /** Flag for choice between sum and average beat salience values for Agent scores.
Chris@11 44 * The use of summed saliences favours faster tempi or lower metrical levels. */
Chris@11 45 static bool useAverageSalience;
Chris@8 46
Chris@11 47 /** For the purpose of removing duplicate agents, the default JND of IBI */
Chris@11 48 static const double DEFAULT_BI;
Chris@8 49
Chris@11 50 /** For the purpose of removing duplicate agents, the default JND of phase */
Chris@11 51 static const double DEFAULT_BT;
Chris@8 52
Chris@11 53 /** Inserts newAgent into the list in ascending order of beatInterval */
Chris@11 54 void add(Agent a) {
Chris@11 55 add(a, true);
Chris@11 56 } // add()/1
Chris@8 57
Chris@11 58 /** Appends newAgent to list (sort==false), or inserts newAgent into the list
Chris@11 59 * in ascending order of beatInterval
Chris@11 60 * @param newAgent The agent to be added to the list
Chris@11 61 * @param sort Flag indicating whether the list is sorted or not
Chris@11 62 */
Chris@11 63 void add(Agent newAgent, bool sort){
Chris@11 64 list.push_back(newAgent);
Chris@11 65 if (sort) this->sort();
Chris@11 66 } // add()/2
Chris@8 67
Chris@11 68 /** Sorts the AgentList by increasing beatInterval. */
Chris@11 69 void sort() {
Chris@11 70 std::sort(list.begin(), list.end());
Chris@11 71 } // sort()
Chris@8 72
Chris@11 73 /** Removes the given item from the list.
Chris@11 74 * @param ptr Points to the Agent which is removed from the list
Chris@11 75 */
Chris@11 76 void remove(iterator itr) {
Chris@11 77 list.erase(itr);
Chris@11 78 } // remove()
Chris@8 79
Chris@8 80 protected:
Chris@11 81 /** Removes Agents from the list which are duplicates of other Agents.
Chris@11 82 * A duplicate is defined by the tempo and phase thresholds
Chris@11 83 * thresholdBI and thresholdBT respectively.
Chris@11 84 */
Chris@11 85 void removeDuplicates() {
Chris@11 86 sort();
Chris@11 87 for (iterator itr = begin(); itr != end(); ++itr) {
Chris@11 88 if (itr->phaseScore < 0.0) // already flagged for deletion
Chris@11 89 continue;
Chris@11 90 iterator itr2 = itr;
Chris@11 91 for (++itr2; itr2 != end(); ++itr2) {
Chris@11 92 if (itr2->beatInterval - itr->beatInterval > DEFAULT_BI)
Chris@11 93 break;
Chris@11 94 if (fabs(itr->beatTime - itr2->beatTime) > DEFAULT_BT)
Chris@11 95 continue;
Chris@11 96 if (itr->phaseScore < itr2->phaseScore) {
Chris@11 97 itr->phaseScore = -1.0; // flag for deletion
Chris@11 98 if (itr2->topScoreTime < itr->topScoreTime)
Chris@11 99 itr2->topScoreTime = itr->topScoreTime;
Chris@11 100 break;
Chris@11 101 } else {
Chris@11 102 itr2->phaseScore = -1.0; // flag for deletion
Chris@11 103 if (itr->topScoreTime < itr2->topScoreTime)
Chris@11 104 itr->topScoreTime = itr2->topScoreTime;
Chris@11 105 }
Chris@11 106 }
Chris@11 107 }
Chris@11 108 for (iterator itr = begin(); itr != end(); ) {
Chris@11 109 if (itr->phaseScore < 0.0) {
Chris@11 110 list.erase(itr);
Chris@11 111 } else {
Chris@11 112 ++itr;
Chris@11 113 }
Chris@11 114 }
Chris@11 115 } // removeDuplicates()
Chris@8 116
Chris@8 117 public:
Chris@11 118 /** Perform beat tracking on a list of events (onsets).
Chris@11 119 * @param el The list of onsets (or events or peaks) to beat track
Chris@11 120 */
Chris@11 121 void beatTrack(EventList el) {
Chris@11 122 beatTrack(el, -1.0);
Chris@11 123 } // beatTrack()/1
Chris@8 124
Chris@11 125 /** Perform beat tracking on a list of events (onsets).
Chris@11 126 * @param el The list of onsets (or events or peaks) to beat track.
Chris@11 127 * @param stop Do not find beats after <code>stop</code> seconds.
Chris@11 128 */
Chris@11 129 void beatTrack(EventList el, double stop) {
Chris@11 130 EventList::iterator ei = el.begin();
Chris@11 131 bool phaseGiven = !empty() && (begin()->beatTime >= 0); // if given for one, assume given for others
Chris@11 132 while (ei != el.end()) {
Chris@11 133 Event ev = *ei;
Chris@11 134 ++ei;
Chris@11 135 if ((stop > 0) && (ev.time > stop))
Chris@11 136 break;
Chris@11 137 bool created = phaseGiven;
Chris@11 138 double prevBeatInterval = -1.0;
Chris@11 139 // cc: Duplicate our list of agents, and scan through
Chris@11 140 // the (now immutable) copy. This means we can safely
Chris@11 141 // add agents to our own list while scanning without
Chris@11 142 // disrupting our scan.
Chris@11 143 Container currentAgents = list;
Chris@11 144 for (Container::iterator ai = currentAgents.begin();
Chris@11 145 ai != currentAgents.end(); ++ai) {
Chris@11 146 Agent currentAgent = *ai;
Chris@11 147 if (currentAgent.beatInterval != prevBeatInterval) {
Chris@11 148 if ((prevBeatInterval>=0) && !created && (ev.time<5.0)) {
Chris@11 149 // Create new agent with different phase
Chris@11 150 Agent newAgent(prevBeatInterval);
Chris@11 151 // This may add an agent to our list
Chris@11 152 newAgent.considerAsBeat(ev, *this);
Chris@8 153 }
Chris@11 154 prevBeatInterval = currentAgent.beatInterval;
Chris@11 155 created = phaseGiven;
Chris@11 156 }
Chris@11 157 if (currentAgent.considerAsBeat(ev, *this))
Chris@11 158 created = true;
Chris@11 159 } // loop for each agent
Chris@11 160 removeDuplicates();
Chris@11 161 } // loop for each event
Chris@11 162 } // beatTrack()
Chris@8 163
Chris@11 164 /** Finds the Agent with the highest score in the list, or NULL if beat tracking has failed.
Chris@11 165 * @return The Agent with the highest score
Chris@11 166 */
Chris@11 167 Agent *bestAgent() {
Chris@11 168 double best = -1.0;
Chris@11 169 Agent *bestAg = 0;
Chris@11 170 for (iterator itr = begin(); itr != end(); ++itr) {
Chris@11 171 if (itr->events.empty()) continue;
Chris@11 172 double startTime = itr->events.begin()->time;
Chris@11 173 double conf = (itr->phaseScore + itr->tempoScore) /
Chris@11 174 (useAverageSalience? (double)itr->beatCount: 1.0);
Chris@11 175 if (conf > best) {
Chris@11 176 bestAg = &(*itr);
Chris@11 177 best = conf;
Chris@8 178 }
Chris@11 179 }
Chris@11 180 return bestAg;
Chris@11 181 } // bestAgent()
Chris@8 182
Chris@8 183
Chris@8 184 }; // class AgentList
Chris@8 185
Chris@8 186 #endif
Chris@8 187