cannam@108: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
cannam@108: 
cannam@108: /*
cannam@108:     Vamp
cannam@108: 
cannam@108:     An API for audio analysis and feature extraction plugins.
cannam@108: 
cannam@108:     Centre for Digital Music, Queen Mary, University of London.
cannam@108:     Copyright 2006 Chris Cannam.
cannam@108:   
cannam@108:     Permission is hereby granted, free of charge, to any person
cannam@108:     obtaining a copy of this software and associated documentation
cannam@108:     files (the "Software"), to deal in the Software without
cannam@108:     restriction, including without limitation the rights to use, copy,
cannam@108:     modify, merge, publish, distribute, sublicense, and/or sell copies
cannam@108:     of the Software, and to permit persons to whom the Software is
cannam@108:     furnished to do so, subject to the following conditions:
cannam@108: 
cannam@108:     The above copyright notice and this permission notice shall be
cannam@108:     included in all copies or substantial portions of the Software.
cannam@108: 
cannam@108:     THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
cannam@108:     EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
cannam@108:     MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
cannam@108:     NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS BE LIABLE FOR
cannam@108:     ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF
cannam@108:     CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
cannam@108:     WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
cannam@108: 
cannam@108:     Except as contained in this notice, the names of the Centre for
cannam@108:     Digital Music; Queen Mary, University of London; and Chris Cannam
cannam@108:     shall not be used in advertising or otherwise to promote the sale,
cannam@108:     use or other dealings in this Software without prior written
cannam@108:     authorization.
cannam@108: */
cannam@108: 
cannam@108: #include "PercussionOnsetDetector.h"
cannam@108: 
cannam@108: using std::string;
cannam@108: using std::vector;
cannam@108: using std::cerr;
cannam@108: using std::endl;
cannam@108: 
cannam@108: #include <cmath>
cannam@108: 
cannam@108: 
cannam@108: PercussionOnsetDetector::PercussionOnsetDetector(float inputSampleRate) :
cannam@108:     Plugin(inputSampleRate),
cannam@108:     m_stepSize(0),
cannam@108:     m_blockSize(0),
cannam@108:     m_threshold(3),
cannam@108:     m_sensitivity(40),
cannam@108:     m_priorMagnitudes(0),
cannam@108:     m_dfMinus1(0),
cannam@108:     m_dfMinus2(0)
cannam@108: {
cannam@108: }
cannam@108: 
cannam@108: PercussionOnsetDetector::~PercussionOnsetDetector()
cannam@108: {
cannam@108:     delete[] m_priorMagnitudes;
cannam@108: }
cannam@108: 
cannam@108: string
cannam@108: PercussionOnsetDetector::getIdentifier() const
cannam@108: {
cannam@108:     return "percussiononsets";
cannam@108: }
cannam@108: 
cannam@108: string
cannam@108: PercussionOnsetDetector::getName() const
cannam@108: {
cannam@108:     return "Simple Percussion Onset Detector";
cannam@108: }
cannam@108: 
cannam@108: string
cannam@108: PercussionOnsetDetector::getDescription() const
cannam@108: {
cannam@108:     return "Detect percussive note onsets by identifying broadband energy rises";
cannam@108: }
cannam@108: 
cannam@108: string
cannam@108: PercussionOnsetDetector::getMaker() const
cannam@108: {
cannam@108:     return "Vamp SDK Example Plugins";
cannam@108: }
cannam@108: 
cannam@108: int
cannam@108: PercussionOnsetDetector::getPluginVersion() const
cannam@108: {
cannam@108:     return 2;
cannam@108: }
cannam@108: 
cannam@108: string
cannam@108: PercussionOnsetDetector::getCopyright() const
cannam@108: {
cannam@108:     return "Code copyright 2006 Queen Mary, University of London, after Dan Barry et al 2005.  Freely redistributable (BSD license)";
cannam@108: }
cannam@108: 
cannam@108: size_t
cannam@108: PercussionOnsetDetector::getPreferredStepSize() const
cannam@108: {
cannam@108:     return 0;
cannam@108: }
cannam@108: 
cannam@108: size_t
cannam@108: PercussionOnsetDetector::getPreferredBlockSize() const
cannam@108: {
cannam@108:     return 1024;
cannam@108: }
cannam@108: 
cannam@108: bool
cannam@108: PercussionOnsetDetector::initialise(size_t channels, size_t stepSize, size_t blockSize)
cannam@108: {
cannam@108:     if (channels < getMinChannelCount() ||
cannam@108: 	channels > getMaxChannelCount()) return false;
cannam@108: 
cannam@108:     m_stepSize = stepSize;
cannam@108:     m_blockSize = blockSize;
cannam@108: 
cannam@108:     m_priorMagnitudes = new float[m_blockSize/2];
cannam@108: 
cannam@108:     for (size_t i = 0; i < m_blockSize/2; ++i) {
cannam@108:         m_priorMagnitudes[i] = 0.f;
cannam@108:     }
cannam@108: 
cannam@108:     m_dfMinus1 = 0.f;
cannam@108:     m_dfMinus2 = 0.f;
cannam@108: 
cannam@108:     return true;
cannam@108: }
cannam@108: 
cannam@108: void
cannam@108: PercussionOnsetDetector::reset()
cannam@108: {
cannam@108:     for (size_t i = 0; i < m_blockSize/2; ++i) {
cannam@108:         m_priorMagnitudes[i] = 0.f;
cannam@108:     }
cannam@108: 
cannam@108:     m_dfMinus1 = 0.f;
cannam@108:     m_dfMinus2 = 0.f;
cannam@108: }
cannam@108: 
cannam@108: PercussionOnsetDetector::ParameterList
cannam@108: PercussionOnsetDetector::getParameterDescriptors() const
cannam@108: {
cannam@108:     ParameterList list;
cannam@108: 
cannam@108:     ParameterDescriptor d;
cannam@108:     d.identifier = "threshold";
cannam@108:     d.name = "Energy rise threshold";
cannam@108:     d.description = "Energy rise within a frequency bin necessary to count toward broadband total";
cannam@108:     d.unit = "dB";
cannam@108:     d.minValue = 0;
cannam@108:     d.maxValue = 20;
cannam@108:     d.defaultValue = 3;
cannam@108:     d.isQuantized = false;
cannam@108:     list.push_back(d);
cannam@108: 
cannam@108:     d.identifier = "sensitivity";
cannam@108:     d.name = "Sensitivity";
cannam@108:     d.description = "Sensitivity of peak detector applied to broadband detection function";
cannam@108:     d.unit = "%";
cannam@108:     d.minValue = 0;
cannam@108:     d.maxValue = 100;
cannam@108:     d.defaultValue = 40;
cannam@108:     d.isQuantized = false;
cannam@108:     list.push_back(d);
cannam@108: 
cannam@108:     return list;
cannam@108: }
cannam@108: 
cannam@108: float
cannam@108: PercussionOnsetDetector::getParameter(std::string id) const
cannam@108: {
cannam@108:     if (id == "threshold") return m_threshold;
cannam@108:     if (id == "sensitivity") return m_sensitivity;
cannam@108:     return 0.f;
cannam@108: }
cannam@108: 
cannam@108: void
cannam@108: PercussionOnsetDetector::setParameter(std::string id, float value)
cannam@108: {
cannam@108:     if (id == "threshold") {
cannam@108:         if (value < 0) value = 0;
cannam@108:         if (value > 20) value = 20;
cannam@108:         m_threshold = value;
cannam@108:     } else if (id == "sensitivity") {
cannam@108:         if (value < 0) value = 0;
cannam@108:         if (value > 100) value = 100;
cannam@108:         m_sensitivity = value;
cannam@108:     }
cannam@108: }
cannam@108: 
cannam@108: PercussionOnsetDetector::OutputList
cannam@108: PercussionOnsetDetector::getOutputDescriptors() const
cannam@108: {
cannam@108:     OutputList list;
cannam@108: 
cannam@108:     OutputDescriptor d;
cannam@108:     d.identifier = "onsets";
cannam@108:     d.name = "Onsets";
cannam@108:     d.description = "Percussive note onset locations";
cannam@108:     d.unit = "";
cannam@108:     d.hasFixedBinCount = true;
cannam@108:     d.binCount = 0;
cannam@108:     d.hasKnownExtents = false;
cannam@108:     d.isQuantized = false;
cannam@108:     d.sampleType = OutputDescriptor::VariableSampleRate;
cannam@108:     d.sampleRate = m_inputSampleRate;
cannam@108:     list.push_back(d);
cannam@108: 
cannam@108:     d.identifier = "detectionfunction";
cannam@108:     d.name = "Detection Function";
cannam@108:     d.description = "Broadband energy rise detection function";
cannam@108:     d.binCount = 1;
cannam@108:     d.isQuantized = true;
cannam@108:     d.quantizeStep = 1.0;
cannam@108:     d.sampleType = OutputDescriptor::OneSamplePerStep;
cannam@108:     list.push_back(d);
cannam@108: 
cannam@108:     return list;
cannam@108: }
cannam@108: 
cannam@108: PercussionOnsetDetector::FeatureSet
cannam@108: PercussionOnsetDetector::process(const float *const *inputBuffers,
cannam@108:                                  Vamp::RealTime ts)
cannam@108: {
cannam@108:     if (m_stepSize == 0) {
cannam@108: 	cerr << "ERROR: PercussionOnsetDetector::process: "
cannam@108: 	     << "PercussionOnsetDetector has not been initialised"
cannam@108: 	     << endl;
cannam@108: 	return FeatureSet();
cannam@108:     }
cannam@108: 
cannam@108:     int count = 0;
cannam@108: 
cannam@108:     for (size_t i = 1; i < m_blockSize/2; ++i) {
cannam@108: 
cannam@108:         float real = inputBuffers[0][i*2];
cannam@108:         float imag = inputBuffers[0][i*2 + 1];
cannam@108: 
cannam@108:         float sqrmag = real * real + imag * imag;
cannam@108: 
cannam@108:         if (m_priorMagnitudes[i] > 0.f) {
cannam@108:             float diff = 10.f * log10f(sqrmag / m_priorMagnitudes[i]);
cannam@108: 
cannam@108: //        std::cout << "i=" << i << ", sqrmag=" << sqrmag << ", prior=" << m_priorMagnitudes[i] << ", diff=" << diff << ", threshold=" << m_threshold << " " << (diff >= m_threshold ? "[*]" : "") << std::endl;
cannam@108: 
cannam@108:             if (diff >= m_threshold) ++count;
cannam@108:         }
cannam@108: 
cannam@108:         m_priorMagnitudes[i] = sqrmag;
cannam@108:     }
cannam@108: 
cannam@108:     FeatureSet returnFeatures;
cannam@108: 
cannam@108:     Feature detectionFunction;
cannam@108:     detectionFunction.hasTimestamp = false;
cannam@108:     detectionFunction.values.push_back(count);
cannam@108:     returnFeatures[1].push_back(detectionFunction);
cannam@108: 
cannam@108:     if (m_dfMinus2 < m_dfMinus1 &&
cannam@108:         m_dfMinus1 >= count &&
cannam@108:         m_dfMinus1 > ((100 - m_sensitivity) * m_blockSize) / 200) {
cannam@108: 
cannam@108: //std::cout << "result at " << ts << "! (count == " << count << ", prev == " << m_dfMinus1 << ")" << std::endl;
cannam@108: 
cannam@108:         Feature onset;
cannam@108:         onset.hasTimestamp = true;
cannam@108:         onset.timestamp = ts - Vamp::RealTime::frame2RealTime
cannam@108:             (m_stepSize, int(m_inputSampleRate + 0.5));
cannam@108:         returnFeatures[0].push_back(onset);
cannam@108:     }
cannam@108: 
cannam@108:     m_dfMinus2 = m_dfMinus1;
cannam@108:     m_dfMinus1 = count;
cannam@108: 
cannam@108:     return returnFeatures;
cannam@108: }
cannam@108: 
cannam@108: PercussionOnsetDetector::FeatureSet
cannam@108: PercussionOnsetDetector::getRemainingFeatures()
cannam@108: {
cannam@108:     return FeatureSet();
cannam@108: }
cannam@108: