Chris@105: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
Chris@105: /*
Chris@105:     Vamp feature extraction plugin using the MATCH audio alignment
Chris@105:     algorithm.
Chris@105: 
Chris@105:     Centre for Digital Music, Queen Mary, University of London.
Chris@105:     This file copyright 2007 Simon Dixon, Chris Cannam and QMUL.
Chris@105:     
Chris@105:     This program is free software; you can redistribute it and/or
Chris@105:     modify it under the terms of the GNU General Public License as
Chris@105:     published by the Free Software Foundation; either version 2 of the
Chris@105:     License, or (at your option) any later version.  See the file
Chris@105:     COPYING included with this distribution for more information.
Chris@105: */
Chris@105: 
Chris@105: #include "MatchPipeline.h"
Chris@105: 
Chris@105: MatchPipeline::MatchPipeline(FeatureExtractor::Parameters feParams,
Chris@105: 			     FeatureConditioner::Parameters fcParams,
Chris@105: 			     Matcher::Parameters matchParams) :
Chris@105:     m_fe1(feParams),
Chris@105:     m_fe2(feParams),
Chris@105:     m_fc1(fcParams),
Chris@105:     m_fc2(fcParams),
Chris@105:     m_pm1(matchParams, 0),
Chris@105:     m_pm2(matchParams, &m_pm1),
Chris@105:     m_feeder(&m_pm1, &m_pm2),
Chris@105:     m_lastFrameIn1(0),
Chris@105:     m_lastFrameIn2(0),
Chris@105:     m_frameNo(0)
Chris@105: {
Chris@105:     m_pm1.setOtherMatcher(&m_pm2);
Chris@105: }
Chris@105: 
Chris@105: MatchPipeline::~MatchPipeline()
Chris@105: {
Chris@105: }
Chris@105: 
Chris@105: void
Chris@105: MatchPipeline::feedFrequencyDomainAudio(const float *arr1, const float *arr2)
Chris@105: {
Chris@105:     feedFeatures(m_fe1.process(arr1), m_fe2.process(arr2));
Chris@105: }
Chris@105: 
Chris@105: void
Chris@105: MatchPipeline::feedFeatures(const vector<double> &f1, const vector<double> &f2)
Chris@105: {
Chris@106:     m_f1 = f1;
Chris@106:     m_f2 = f2;
Chris@106: 
Chris@105:     feedConditionedFeatures(m_fc1.process(f1), m_fc2.process(f2));
Chris@105: }
Chris@105: 
Chris@105: void
Chris@106: MatchPipeline::feedConditionedFeatures(const vector<double> &c1, const vector<double> &c2)
Chris@105: {
Chris@106:     m_c1 = c1;
Chris@106:     m_c2 = c2;
Chris@106:     
Chris@106:     m_feeder.feed(c1, c2);
Chris@105: 
Chris@106:     if (aboveThreshold(c1)) m_lastFrameIn1 = m_frameNo;
Chris@106:     if (aboveThreshold(c2)) m_lastFrameIn2 = m_frameNo;
Chris@105: 
Chris@105:     ++m_frameNo;
Chris@105: }
Chris@105: 
Chris@106: void
Chris@106: MatchPipeline::extractFeatures(vector<double> &f1, vector<double> &f2)
Chris@106: {
Chris@106:     f1 = m_f1;
Chris@106:     f2 = m_f2;
Chris@106: }
Chris@106: 
Chris@106: void
Chris@106: MatchPipeline::extractConditionedFeatures(vector<double> &c1, vector<double> &c2)
Chris@106: {
Chris@106:     c1 = m_c1;
Chris@106:     c2 = m_c2;
Chris@106: }
Chris@106: 
Chris@105: bool
Chris@105: MatchPipeline::aboveThreshold(const vector<double> &f)
Chris@105: {
Chris@105:     double threshold = 1e-4f;
Chris@105:     double sum = 0.f;
Chris@105:     for (int i = 0; i < int(f.size()); ++i) {
Chris@105:         sum += f[i] * f[i];
Chris@105:     }
Chris@105:     return (sum >= threshold);
Chris@105: }
Chris@105: 
Chris@105: void
Chris@105: MatchPipeline::finish()
Chris@105: {
Chris@105:     m_feeder.finish();
Chris@105:     getFinder()->setDurations(m_lastFrameIn1, m_lastFrameIn2);
Chris@105: }
Chris@105: 
Chris@135: MatchFeatureFeeder *
Chris@135: MatchPipeline::getFeeder()
Chris@135: {
Chris@135:     return &m_feeder;
Chris@135: }
Chris@135: 
Chris@105: Finder *
Chris@105: MatchPipeline::getFinder()
Chris@105: {
Chris@105:     return m_feeder.getFinder();
Chris@105: }
Chris@105: 
Chris@105: 
Chris@105: 
Chris@105: