matthiasm@0: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
Chris@9: 
matthiasm@0: /*
Chris@9:     pYIN - A fundamental frequency estimator for monophonic audio
Chris@9:     Centre for Digital Music, Queen Mary, University of London.
Chris@9:     
Chris@9:     This program is free software; you can redistribute it and/or
Chris@9:     modify it under the terms of the GNU General Public License as
Chris@9:     published by the Free Software Foundation; either version 2 of the
Chris@9:     License, or (at your option) any later version.  See the file
Chris@9:     COPYING included with this distribution for more information.
matthiasm@0: */
matthiasm@0: 
matthiasm@0: #ifndef _YIN_H_
matthiasm@0: #define _YIN_H_
matthiasm@0: 
matthiasm@0: #include "vamp-sdk/FFT.h"
matthiasm@0: #include "MeanFilter.h"
matthiasm@0: 
matthiasm@0: #include <cmath>
matthiasm@0: 
matthiasm@0: #include <iostream>
matthiasm@0: #include <vector>
matthiasm@0: #include <exception>
matthiasm@0: 
matthiasm@0: using std::vector;
matthiasm@0: using std::pair;
matthiasm@0: 
Chris@136: class YinUtil;
Chris@136: 
matthiasm@0: class Yin
matthiasm@0: {
matthiasm@0: public:
matthiasm@117:     Yin(size_t frameSize, size_t inputSampleRate, double thresh = 0.2, bool fast = true);
matthiasm@0:     virtual ~Yin();
matthiasm@0: 
matthiasm@0:     struct YinOutput {
matthiasm@0:         double f0;
matthiasm@0:         double periodicity;
matthiasm@0:         double rms;
matthiasm@0:         vector<double> salience;
matthiasm@0:         vector<pair<double, double> > freqProb;
matthiasm@0:         YinOutput() :  f0(0), periodicity(0), rms(0), 
matthiasm@0:             salience(vector<double>(0)), freqProb(vector<pair<double, double> >(0)) { }
matthiasm@0:         YinOutput(double _f, double _p, double _r) :
matthiasm@0:             f0(_f), periodicity(_p), rms(_r), 
matthiasm@0:             salience(vector<double>(0)), freqProb(vector<pair<double, double> >(0)) { }
matthiasm@0:         YinOutput(double _f, double _p, double _r, vector<double> _salience) :
matthiasm@0:             f0(_f), periodicity(_p), rms(_r), salience(_salience), 
matthiasm@0:             freqProb(vector<pair<double, double> >(0)) { }
matthiasm@0:     };
matthiasm@0:     
matthiasm@0:     int setThreshold(double parameter);
matthiasm@0:     int setThresholdDistr(float parameter);
matthiasm@0:     int setFrameSize(size_t frameSize);
matthiasm@117:     int setFast(bool fast);
matthiasm@0:     // int setRemoveUnvoiced(bool frameSize);
matthiasm@0:     YinOutput process(const double *in) const;
matthiasm@0:     YinOutput processProbabilisticYin(const double *in) const;
matthiasm@0: 
matthiasm@0: private:
matthiasm@0:     mutable size_t m_frameSize;
matthiasm@0:     mutable size_t m_inputSampleRate;
matthiasm@0:     mutable double m_thresh;
matthiasm@0:     mutable size_t m_threshDistr;
matthiasm@0:     mutable size_t m_yinBufferSize;
matthiasm@117:     mutable bool   m_fast;
matthiasm@0:     // mutable bool m_removeUnvoiced;
Chris@136:     YinUtil *m_yinUtil;
matthiasm@0: };
matthiasm@0: 
matthiasm@0: #endif