source: src/pitch/pitch.c @ 168337e

feature/autosinkfeature/cnnfeature/cnn_orgfeature/constantqfeature/crepefeature/crepe_orgfeature/pitchshiftfeature/pydocstringsfeature/timestretchfix/ffmpeg5pitchshiftsamplertimestretchyinfft+
Last change on this file since 168337e was 168337e, checked in by Paul Brossier <piem@piem.org>, 14 years ago

src/pitch: switch to mono

  • Property mode set to 100644
File size: 10.1 KB
Line 
1/*
2  Copyright (C) 2003-2009 Paul Brossier <piem@aubio.org>
3
4  This file is part of aubio.
5
6  aubio is free software: you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation, either version 3 of the License, or
9  (at your option) any later version.
10
11  aubio is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  GNU General Public License for more details.
15
16  You should have received a copy of the GNU General Public License
17  along with aubio.  If not, see <http://www.gnu.org/licenses/>.
18
19*/
20
21#include "aubio_priv.h"
22#include "fvec.h"
23#include "cvec.h"
24#include "lvec.h"
25#include "mathutils.h"
26#include "musicutils.h"
27#include "spectral/phasevoc.h"
28#include "temporal/filter.h"
29#include "temporal/c_weighting.h"
30#include "pitch/pitchmcomb.h"
31#include "pitch/pitchyin.h"
32#include "pitch/pitchfcomb.h"
33#include "pitch/pitchschmitt.h"
34#include "pitch/pitchyinfft.h"
35#include "pitch/pitch.h"
36
37/** pitch detection algorithm */
38typedef enum
39{
40  aubio_pitcht_yin,     /**< YIN algorithm */
41  aubio_pitcht_mcomb,   /**< Multi-comb filter */
42  aubio_pitcht_schmitt, /**< Schmitt trigger */
43  aubio_pitcht_fcomb,   /**< Fast comb filter */
44  aubio_pitcht_yinfft,   /**< Spectral YIN */
45  aubio_pitcht_default = aubio_pitcht_yinfft, /**< the one used when "default" is asked */
46} aubio_pitch_type;
47
48/** pitch detection output mode */
49typedef enum
50{
51  aubio_pitchm_freq,   /**< Frequency (Hz) */
52  aubio_pitchm_midi,   /**< MIDI note (0.,127) */
53  aubio_pitchm_cent,   /**< Cent */
54  aubio_pitchm_bin,    /**< Frequency bin (0,bufsize) */
55  aubio_pitchm_default = aubio_pitchm_freq, /**< the one used when "default" is asked */
56} aubio_pitch_mode;
57
58typedef void (*aubio_pitch_func_t)
59  (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf);
60typedef smpl_t (*aubio_pitch_conv_t)
61  (smpl_t value, uint_t srate, uint_t bufsize);
62
63void aubio_pitch_slideblock (aubio_pitch_t * p, fvec_t * ibuf);
64
65void aubio_pitch_do_mcomb (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf);
66void aubio_pitch_do_yin (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf);
67void aubio_pitch_do_schmitt (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf);
68void aubio_pitch_do_fcomb (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf);
69void aubio_pitch_do_yinfft (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf);
70
71/** generic pitch detection structure */
72struct _aubio_pitch_t
73{
74  aubio_pitch_type type; /**< pitch detection mode */
75  aubio_pitch_mode mode; /**< pitch detection output mode */
76  uint_t srate;                   /**< samplerate */
77  uint_t bufsize;                 /**< buffer size */
78  aubio_pitchmcomb_t *mcomb;      /**< mcomb object */
79  aubio_pitchfcomb_t *fcomb;      /**< fcomb object */
80  aubio_pitchschmitt_t *schmitt;  /**< schmitt object */
81  aubio_pitchyinfft_t *yinfft;    /**< yinfft object */
82  aubio_pitchyin_t *yin;    /**< yinfft object */
83  aubio_filter_t *filter;         /**< filter */
84  aubio_pvoc_t *pv;               /**< phase vocoder for mcomb */
85  cvec_t *fftgrain;               /**< spectral frame for mcomb */
86  fvec_t *buf;                    /**< temporary buffer for yin */
87  aubio_pitch_func_t callback; /**< pointer to current pitch detection method */
88  aubio_pitch_conv_t freqconv; /**< pointer to current pitch conversion method */
89};
90
91/* convenience wrapper function for frequency unit conversions
92 * should probably be rewritten with #defines */
93smpl_t freqconvbin (smpl_t f, uint_t srate, uint_t bufsize);
94smpl_t
95freqconvbin (smpl_t f, uint_t srate, uint_t bufsize)
96{
97  return aubio_freqtobin (f, srate, bufsize);
98}
99
100smpl_t freqconvmidi (smpl_t f, uint_t srate, uint_t bufsize);
101smpl_t
102freqconvmidi (smpl_t f, uint_t srate UNUSED, uint_t bufsize UNUSED)
103{
104  return aubio_freqtomidi (f);
105}
106
107smpl_t freqconvpass (smpl_t f, uint_t srate, uint_t bufsize);
108smpl_t
109freqconvpass (smpl_t f, uint_t srate UNUSED, uint_t bufsize UNUSED)
110{
111  return f;
112}
113
114aubio_pitch_t *
115new_aubio_pitch (char_t * pitch_mode,
116    uint_t bufsize, uint_t hopsize, uint_t samplerate)
117{
118  aubio_pitch_t *p = AUBIO_NEW (aubio_pitch_t);
119  aubio_pitch_type pitch_type;
120  if (strcmp (pitch_mode, "mcomb") == 0)
121    pitch_type = aubio_pitcht_mcomb;
122  else if (strcmp (pitch_mode, "yinfft") == 0)
123    pitch_type = aubio_pitcht_yin;
124  else if (strcmp (pitch_mode, "yin") == 0)
125    pitch_type = aubio_pitcht_yin;
126  else if (strcmp (pitch_mode, "schmitt") == 0)
127    pitch_type = aubio_pitcht_schmitt;
128  else if (strcmp (pitch_mode, "fcomb") == 0)
129    pitch_type = aubio_pitcht_fcomb;
130  else if (strcmp (pitch_mode, "default") == 0)
131    pitch_type = aubio_pitcht_default;
132  else {
133    AUBIO_ERR ("unknown pitch detection method %s, using default.\n",
134        pitch_mode);
135    pitch_type = aubio_pitcht_default;
136    return NULL;
137  }
138  p->srate = samplerate;
139  p->type = pitch_type;
140  aubio_pitch_set_unit (p, "default");
141  p->bufsize = bufsize;
142  switch (p->type) {
143    case aubio_pitcht_yin:
144      p->buf = new_fvec (bufsize);
145      p->yin = new_aubio_pitchyin (bufsize);
146      p->callback = aubio_pitch_do_yin;
147      aubio_pitchyin_set_tolerance (p->yin, 0.15);
148      break;
149    case aubio_pitcht_mcomb:
150      p->pv = new_aubio_pvoc (bufsize, hopsize);
151      p->fftgrain = new_cvec (bufsize);
152      p->mcomb = new_aubio_pitchmcomb (bufsize, hopsize);
153      p->filter = new_aubio_filter_c_weighting (samplerate);
154      p->callback = aubio_pitch_do_mcomb;
155      break;
156    case aubio_pitcht_fcomb:
157      p->buf = new_fvec (bufsize);
158      p->fcomb = new_aubio_pitchfcomb (bufsize, hopsize);
159      p->callback = aubio_pitch_do_fcomb;
160      break;
161    case aubio_pitcht_schmitt:
162      p->buf = new_fvec (bufsize);
163      p->schmitt = new_aubio_pitchschmitt (bufsize);
164      p->callback = aubio_pitch_do_schmitt;
165      break;
166    case aubio_pitcht_yinfft:
167      p->buf = new_fvec (bufsize);
168      p->yinfft = new_aubio_pitchyinfft (bufsize);
169      p->callback = aubio_pitch_do_yinfft;
170      aubio_pitchyinfft_set_tolerance (p->yinfft, 0.85);
171      break;
172    default:
173      break;
174  }
175  return p;
176}
177
178void
179del_aubio_pitch (aubio_pitch_t * p)
180{
181  switch (p->type) {
182    case aubio_pitcht_yin:
183      del_fvec (p->buf);
184      del_aubio_pitchyin (p->yin);
185      break;
186    case aubio_pitcht_mcomb:
187      del_aubio_pvoc (p->pv);
188      del_cvec (p->fftgrain);
189      del_aubio_filter (p->filter);
190      del_aubio_pitchmcomb (p->mcomb);
191      break;
192    case aubio_pitcht_schmitt:
193      del_fvec (p->buf);
194      del_aubio_pitchschmitt (p->schmitt);
195      break;
196    case aubio_pitcht_fcomb:
197      del_fvec (p->buf);
198      del_aubio_pitchfcomb (p->fcomb);
199      break;
200    case aubio_pitcht_yinfft:
201      del_fvec (p->buf);
202      del_aubio_pitchyinfft (p->yinfft);
203      break;
204    default:
205      break;
206  }
207  AUBIO_FREE (p);
208}
209
210void
211aubio_pitch_slideblock (aubio_pitch_t * p, fvec_t * ibuf)
212{
213  uint_t j = 0, overlap_size = 0;
214  overlap_size = p->buf->length - ibuf->length;
215  for (j = 0; j < overlap_size; j++) {
216    p->buf->data[j] = p->buf->data[j + ibuf->length];
217  }
218  for (j = 0; j < ibuf->length; j++) {
219    p->buf->data[j + overlap_size] = ibuf->data[j];
220  }
221}
222
223uint_t
224aubio_pitch_set_unit (aubio_pitch_t * p, char_t * pitch_unit)
225{
226  aubio_pitch_mode pitch_mode;
227  if (strcmp (pitch_unit, "freq") == 0)
228    pitch_mode = aubio_pitchm_freq;
229  else if (strcmp (pitch_unit, "midi") == 0)
230    pitch_mode = aubio_pitchm_midi;
231  else if (strcmp (pitch_unit, "cent") == 0)
232    pitch_mode = aubio_pitchm_cent;
233  else if (strcmp (pitch_unit, "bin") == 0)
234    pitch_mode = aubio_pitchm_bin;
235  else if (strcmp (pitch_unit, "default") == 0)
236    pitch_mode = aubio_pitchm_default;
237  else {
238    AUBIO_ERR ("unknown pitch detection unit %s, using default\n", pitch_unit);
239    pitch_mode = aubio_pitchm_default;
240  }
241  p->mode = pitch_mode;
242  switch (p->mode) {
243    case aubio_pitchm_freq:
244      p->freqconv = freqconvpass;
245      break;
246    case aubio_pitchm_midi:
247      p->freqconv = freqconvmidi;
248      break;
249    case aubio_pitchm_cent:
250      /* bug: not implemented */
251      p->freqconv = freqconvmidi;
252      break;
253    case aubio_pitchm_bin:
254      p->freqconv = freqconvbin;
255      break;
256    default:
257      break;
258  }
259  return AUBIO_OK;
260}
261
262uint_t
263aubio_pitch_set_tolerance (aubio_pitch_t * p, smpl_t tol)
264{
265  switch (p->type) {
266    case aubio_pitcht_yin:
267      aubio_pitchyin_set_tolerance (p->yin, tol);
268      break;
269    case aubio_pitcht_yinfft:
270      aubio_pitchyinfft_set_tolerance (p->yinfft, tol);
271      break;
272    default:
273      break;
274  }
275  return AUBIO_OK;
276}
277
278void
279aubio_pitch_do (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf)
280{
281  p->callback (p, ibuf, obuf);
282  obuf->data[0] = p->freqconv (obuf->data[0], p->srate, p->bufsize);
283}
284
285void
286aubio_pitch_do_mcomb (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf)
287{
288  aubio_filter_do (p->filter, ibuf);
289  aubio_pvoc_do (p->pv, ibuf, p->fftgrain);
290  aubio_pitchmcomb_do (p->mcomb, p->fftgrain, obuf);
291  obuf->data[0] = aubio_bintofreq (obuf->data[0], p->srate, p->bufsize);
292}
293
294void
295aubio_pitch_do_yin (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf)
296{
297  smpl_t pitch = 0.;
298  aubio_pitch_slideblock (p, ibuf);
299  aubio_pitchyin_do (p->yin, p->buf, obuf);
300  pitch = obuf->data[0];
301  if (pitch > 0) {
302    pitch = p->srate / (pitch + 0.);
303  } else {
304    pitch = 0.;
305  }
306  obuf->data[0] = pitch;
307}
308
309
310void
311aubio_pitch_do_yinfft (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * obuf)
312{
313  smpl_t pitch = 0.;
314  aubio_pitch_slideblock (p, ibuf);
315  aubio_pitchyinfft_do (p->yinfft, p->buf, obuf);
316  pitch = obuf->data[0];
317  if (pitch > 0) {
318    pitch = p->srate / (pitch + 0.);
319  } else {
320    pitch = 0.;
321  }
322  obuf->data[0] = pitch;
323}
324
325void
326aubio_pitch_do_fcomb (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * out)
327{
328  aubio_pitch_slideblock (p, ibuf);
329  aubio_pitchfcomb_do (p->fcomb, p->buf, out);
330  out->data[0] = aubio_bintofreq (out->data[0], p->srate, p->bufsize);
331}
332
333void
334aubio_pitch_do_schmitt (aubio_pitch_t * p, fvec_t * ibuf, fvec_t * out)
335{
336  smpl_t period, pitch = 0.;
337  aubio_pitch_slideblock (p, ibuf);
338  aubio_pitchschmitt_do (p->schmitt, p->buf, out);
339  period = out->data[0];
340  if (period > 0) {
341    pitch = p->srate / period;
342  } else {
343    pitch = 0.;
344  }
345  out->data[0] = pitch;
346}
Note: See TracBrowser for help on using the repository browser.