Mercurial > hg > pyin
view MonoNote.cpp @ 132:926c292fa3ff fixedlag
fixed lag smoothing for pitch track working
author | Matthias Mauch <mail@matthiasmauch.net> |
---|---|
date | Fri, 03 Jul 2015 17:34:38 +0100 |
parents | 080fe18f5ebf |
children | 83978b93aac1 |
line wrap: on
line source
/* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */ /* pYIN - A fundamental frequency estimator for monophonic audio Centre for Digital Music, Queen Mary, University of London. This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. See the file COPYING included with this distribution for more information. */ #include "MonoNote.h" #include <vector> #include <cstdio> #include <cmath> #include <complex> using std::vector; using std::pair; MonoNote::MonoNote() : hmm(0) { } MonoNote::~MonoNote() { } const vector<MonoNote::FrameOutput> MonoNote::process(const vector<vector<pair<double, double> > > pitchProb) { vector<vector<double> > obsProb; for (size_t iFrame = 0; iFrame < pitchProb.size(); ++iFrame) { obsProb.push_back(hmm.calculateObsProb(pitchProb[iFrame])); } vector<MonoNote::FrameOutput> out; vector<int> path = hmm.decodeViterbi(obsProb); for (size_t iFrame = 0; iFrame < path.size(); ++iFrame) { double currPitch = -1; int stateKind = 0; currPitch = hmm.par.minPitch + (path[iFrame]/hmm.par.nSPP) * 1.0/hmm.par.nPPS; stateKind = (path[iFrame]) % hmm.par.nSPP + 1; out.push_back(FrameOutput(iFrame, currPitch, stateKind)); // std::cerr << path[iFrame] << " -- "<< pitchProb[iFrame][0].first << " -- "<< currPitch << " -- " << stateKind << std::endl; } return(out); }