2 Copyright (C) 2003-2009 Paul Brossier <piem@aubio.org>
4 This file is part of aubio.
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.
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.
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/>.
21 #include "aubio_priv.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/pitchspecacf.h"
36 #include "pitch/pitch.h"
38 #define DEFAULT_PITCH_SILENCE -50.
40 /** pitch detection algorithms */
43 aubio_pitcht_yin, /**< `yin`, YIN algorithm */
44 aubio_pitcht_mcomb, /**< `mcomb`, Multi-comb filter */
45 aubio_pitcht_schmitt, /**< `schmitt`, Schmitt trigger */
46 aubio_pitcht_fcomb, /**< `fcomb`, Fast comb filter */
47 aubio_pitcht_yinfft, /**< `yinfft`, Spectral YIN */
48 aubio_pitcht_specacf, /**< `specacf`, Spectral autocorrelation */
50 = aubio_pitcht_yinfft, /**< `default` */
53 /** pitch detection output modes */
56 aubio_pitchm_freq, /**< Frequency (Hz) */
57 aubio_pitchm_midi, /**< MIDI note (0.,127) */
58 aubio_pitchm_cent, /**< Cent */
59 aubio_pitchm_bin, /**< Frequency bin (0,bufsize) */
60 aubio_pitchm_default = aubio_pitchm_freq, /**< the one used when "default" is asked */
63 /** callback to get pitch candidate, defined below */
64 typedef void (*aubio_pitch_detect_t) (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
66 /** callback to convert pitch from one unit to another, defined below */
67 typedef smpl_t(*aubio_pitch_convert_t) (smpl_t value, uint_t samplerate, uint_t bufsize);
69 /** callback to fetch the confidence of the algorithm */
70 typedef smpl_t (*aubio_pitch_get_conf_t) (void * p);
72 /** generic pitch detection structure */
75 aubio_pitch_type type; /**< pitch detection mode */
76 aubio_pitch_mode mode; /**< pitch detection output mode */
77 uint_t samplerate; /**< samplerate */
78 uint_t bufsize; /**< buffer size */
79 void *p_object; /**< pointer to pitch object */
80 aubio_filter_t *filter; /**< filter */
81 fvec_t *filtered; /**< filtered input */
82 aubio_pvoc_t *pv; /**< phase vocoder for mcomb */
83 cvec_t *fftgrain; /**< spectral frame for mcomb */
84 fvec_t *buf; /**< temporary buffer for yin */
85 aubio_pitch_detect_t detect_cb; /**< callback to get the pitch candidates */
86 aubio_pitch_convert_t conv_cb; /**< callback to convert it to the desired unit */
87 aubio_pitch_get_conf_t conf_cb; /**< pointer to the current confidence callback */
88 smpl_t silence; /**< silence threshold */
91 /* callback functions for pitch detection */
92 static void aubio_pitch_do_mcomb (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
93 static void aubio_pitch_do_yin (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
94 static void aubio_pitch_do_schmitt (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
95 static void aubio_pitch_do_fcomb (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
96 static void aubio_pitch_do_yinfft (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
97 static void aubio_pitch_do_specacf (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf);
99 /* internal functions for frequency conversion */
100 static smpl_t freqconvbin (smpl_t f, uint_t samplerate, uint_t bufsize);
101 static smpl_t freqconvmidi (smpl_t f, uint_t samplerate, uint_t bufsize);
102 static smpl_t freqconvpass (smpl_t f, uint_t samplerate, uint_t bufsize);
104 /* adapter to stack ibuf new samples at the end of buf, and trim `buf` to `bufsize` */
105 void aubio_pitch_slideblock (aubio_pitch_t * p, const fvec_t * ibuf);
109 new_aubio_pitch (const char_t * pitch_mode,
110 uint_t bufsize, uint_t hopsize, uint_t samplerate)
112 aubio_pitch_t *p = AUBIO_NEW (aubio_pitch_t);
113 aubio_pitch_type pitch_type;
114 if (strcmp (pitch_mode, "mcomb") == 0)
115 pitch_type = aubio_pitcht_mcomb;
116 else if (strcmp (pitch_mode, "yinfft") == 0)
117 pitch_type = aubio_pitcht_yinfft;
118 else if (strcmp (pitch_mode, "yin") == 0)
119 pitch_type = aubio_pitcht_yin;
120 else if (strcmp (pitch_mode, "schmitt") == 0)
121 pitch_type = aubio_pitcht_schmitt;
122 else if (strcmp (pitch_mode, "fcomb") == 0)
123 pitch_type = aubio_pitcht_fcomb;
124 else if (strcmp (pitch_mode, "specacf") == 0)
125 pitch_type = aubio_pitcht_specacf;
126 else if (strcmp (pitch_mode, "default") == 0)
127 pitch_type = aubio_pitcht_default;
129 AUBIO_ERR ("pitch: unknown pitch detection method ā%sā\n", pitch_mode);
133 // check parameters are valid
134 if ((sint_t)hopsize < 1) {
135 AUBIO_ERR("pitch: got hopsize %d, but can not be < 1\n", hopsize);
137 } else if ((sint_t)bufsize < 1) {
138 AUBIO_ERR("pitch: got buffer_size %d, but can not be < 1\n", bufsize);
140 } else if (bufsize < hopsize) {
141 AUBIO_ERR("pitch: hop size (%d) is larger than win size (%d)\n", hopsize, bufsize);
143 } else if ((sint_t)samplerate < 1) {
144 AUBIO_ERR("pitch: samplerate (%d) can not be < 1\n", samplerate);
148 p->samplerate = samplerate;
149 p->type = pitch_type;
150 aubio_pitch_set_unit (p, "default");
151 p->bufsize = bufsize;
152 p->silence = DEFAULT_PITCH_SILENCE;
155 case aubio_pitcht_yin:
156 p->buf = new_fvec (bufsize);
157 p->p_object = new_aubio_pitchyin (bufsize);
158 if (!p->p_object) goto beach;
159 p->detect_cb = aubio_pitch_do_yin;
160 p->conf_cb = (aubio_pitch_get_conf_t)aubio_pitchyin_get_confidence;
161 aubio_pitchyin_set_tolerance (p->p_object, 0.15);
163 case aubio_pitcht_mcomb:
164 p->filtered = new_fvec (hopsize);
165 p->pv = new_aubio_pvoc (bufsize, hopsize);
166 if (!p->pv) goto beach;
167 p->fftgrain = new_cvec (bufsize);
168 p->p_object = new_aubio_pitchmcomb (bufsize, hopsize);
169 p->filter = new_aubio_filter_c_weighting (samplerate);
170 p->detect_cb = aubio_pitch_do_mcomb;
172 case aubio_pitcht_fcomb:
173 p->buf = new_fvec (bufsize);
174 p->p_object = new_aubio_pitchfcomb (bufsize, hopsize);
175 if (!p->p_object) goto beach;
176 p->detect_cb = aubio_pitch_do_fcomb;
178 case aubio_pitcht_schmitt:
179 p->buf = new_fvec (bufsize);
180 p->p_object = new_aubio_pitchschmitt (bufsize);
181 p->detect_cb = aubio_pitch_do_schmitt;
183 case aubio_pitcht_yinfft:
184 p->buf = new_fvec (bufsize);
185 p->p_object = new_aubio_pitchyinfft (samplerate, bufsize);
186 if (!p->p_object) goto beach;
187 p->detect_cb = aubio_pitch_do_yinfft;
188 p->conf_cb = (aubio_pitch_get_conf_t)aubio_pitchyinfft_get_confidence;
189 aubio_pitchyinfft_set_tolerance (p->p_object, 0.85);
191 case aubio_pitcht_specacf:
192 p->buf = new_fvec (bufsize);
193 p->p_object = new_aubio_pitchspecacf (bufsize);
194 if (!p->p_object) goto beach;
195 p->detect_cb = aubio_pitch_do_specacf;
196 p->conf_cb = (aubio_pitch_get_conf_t)aubio_pitchspecacf_get_tolerance;
197 aubio_pitchspecacf_set_tolerance (p->p_object, 0.85);
205 if (p->filtered) del_fvec(p->filtered);
206 if (p->buf) del_fvec(p->buf);
212 del_aubio_pitch (aubio_pitch_t * p)
215 case aubio_pitcht_yin:
217 del_aubio_pitchyin (p->p_object);
219 case aubio_pitcht_mcomb:
220 del_fvec (p->filtered);
221 del_aubio_pvoc (p->pv);
222 del_cvec (p->fftgrain);
223 del_aubio_filter (p->filter);
224 del_aubio_pitchmcomb (p->p_object);
226 case aubio_pitcht_schmitt:
228 del_aubio_pitchschmitt (p->p_object);
230 case aubio_pitcht_fcomb:
232 del_aubio_pitchfcomb (p->p_object);
234 case aubio_pitcht_yinfft:
236 del_aubio_pitchyinfft (p->p_object);
238 case aubio_pitcht_specacf:
240 del_aubio_pitchspecacf (p->p_object);
249 aubio_pitch_slideblock (aubio_pitch_t * p, const fvec_t * ibuf)
251 uint_t overlap_size = p->buf->length - ibuf->length;
252 #if 1 //!HAVE_MEMCPY_HACKS
254 for (j = 0; j < overlap_size; j++) {
255 p->buf->data[j] = p->buf->data[j + ibuf->length];
257 for (j = 0; j < ibuf->length; j++) {
258 p->buf->data[j + overlap_size] = ibuf->data[j];
261 smpl_t *data = p->buf->data;
262 smpl_t *newdata = ibuf->data;
263 memmove(data, data + ibuf->length, overlap_size);
264 memcpy(data + overlap_size, newdata, ibuf->length);
269 aubio_pitch_set_unit (aubio_pitch_t * p, const char_t * pitch_unit)
271 uint_t err = AUBIO_OK;
272 aubio_pitch_mode pitch_mode;
273 if (strcmp (pitch_unit, "freq") == 0)
274 pitch_mode = aubio_pitchm_freq;
275 else if (strcmp (pitch_unit, "hertz") == 0)
276 pitch_mode = aubio_pitchm_freq;
277 else if (strcmp (pitch_unit, "Hertz") == 0)
278 pitch_mode = aubio_pitchm_freq;
279 else if (strcmp (pitch_unit, "Hz") == 0)
280 pitch_mode = aubio_pitchm_freq;
281 else if (strcmp (pitch_unit, "f0") == 0)
282 pitch_mode = aubio_pitchm_freq;
283 else if (strcmp (pitch_unit, "midi") == 0)
284 pitch_mode = aubio_pitchm_midi;
285 else if (strcmp (pitch_unit, "cent") == 0)
286 pitch_mode = aubio_pitchm_cent;
287 else if (strcmp (pitch_unit, "bin") == 0)
288 pitch_mode = aubio_pitchm_bin;
289 else if (strcmp (pitch_unit, "default") == 0)
290 pitch_mode = aubio_pitchm_default;
292 AUBIO_WRN("pitch: unknown pitch detection unit ā%sā, using default\n",
294 pitch_mode = aubio_pitchm_default;
297 p->mode = pitch_mode;
299 case aubio_pitchm_freq:
300 p->conv_cb = freqconvpass;
302 case aubio_pitchm_midi:
303 p->conv_cb = freqconvmidi;
305 case aubio_pitchm_cent:
306 /* bug: not implemented */
307 p->conv_cb = freqconvmidi;
309 case aubio_pitchm_bin:
310 p->conv_cb = freqconvbin;
319 aubio_pitch_set_tolerance (aubio_pitch_t * p, smpl_t tol)
322 case aubio_pitcht_yin:
323 aubio_pitchyin_set_tolerance (p->p_object, tol);
325 case aubio_pitcht_yinfft:
326 aubio_pitchyinfft_set_tolerance (p->p_object, tol);
335 aubio_pitch_get_tolerance (aubio_pitch_t * p)
337 smpl_t tolerance = 1.;
339 case aubio_pitcht_yin:
340 tolerance = aubio_pitchyin_get_tolerance (p->p_object);
342 case aubio_pitcht_yinfft:
343 tolerance = aubio_pitchyinfft_get_tolerance (p->p_object);
352 aubio_pitch_set_silence (aubio_pitch_t * p, smpl_t silence)
354 if (silence <= 0 && silence >= -200) {
355 p->silence = silence;
358 AUBIO_WRN("pitch: could not set silence to %.2f\n", silence);
364 aubio_pitch_get_silence (aubio_pitch_t * p)
370 /* do method, calling the detection callback, then the conversion callback */
372 aubio_pitch_do (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf)
374 p->detect_cb (p, ibuf, obuf);
375 if (aubio_silence_detection(ibuf, p->silence) == 1) {
378 obuf->data[0] = p->conv_cb (obuf->data[0], p->samplerate, p->bufsize);
381 /* do method for each algorithm */
383 aubio_pitch_do_mcomb (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf)
385 aubio_filter_do_outplace (p->filter, ibuf, p->filtered);
386 aubio_pvoc_do (p->pv, ibuf, p->fftgrain);
387 aubio_pitchmcomb_do (p->p_object, p->fftgrain, obuf);
388 obuf->data[0] = aubio_bintofreq (obuf->data[0], p->samplerate, p->bufsize);
392 aubio_pitch_do_yin (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf)
395 aubio_pitch_slideblock (p, ibuf);
396 aubio_pitchyin_do (p->p_object, p->buf, obuf);
397 pitch = obuf->data[0];
399 pitch = p->samplerate / (pitch + 0.);
403 obuf->data[0] = pitch;
408 aubio_pitch_do_yinfft (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * obuf)
411 aubio_pitch_slideblock (p, ibuf);
412 aubio_pitchyinfft_do (p->p_object, p->buf, obuf);
413 pitch = obuf->data[0];
415 pitch = p->samplerate / (pitch + 0.);
419 obuf->data[0] = pitch;
423 aubio_pitch_do_specacf (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * out)
425 smpl_t pitch = 0., period;
426 aubio_pitch_slideblock (p, ibuf);
427 aubio_pitchspecacf_do (p->p_object, p->buf, out);
428 //out->data[0] = aubio_bintofreq (out->data[0], p->samplerate, p->bufsize);
429 period = out->data[0];
431 pitch = p->samplerate / period;
435 out->data[0] = pitch;
439 aubio_pitch_do_fcomb (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * out)
441 aubio_pitch_slideblock (p, ibuf);
442 aubio_pitchfcomb_do (p->p_object, p->buf, out);
443 out->data[0] = aubio_bintofreq (out->data[0], p->samplerate, p->bufsize);
447 aubio_pitch_do_schmitt (aubio_pitch_t * p, const fvec_t * ibuf, fvec_t * out)
449 smpl_t period, pitch = 0.;
450 aubio_pitch_slideblock (p, ibuf);
451 aubio_pitchschmitt_do (p->p_object, p->buf, out);
452 period = out->data[0];
454 pitch = p->samplerate / period;
458 out->data[0] = pitch;
461 /* conversion callbacks */
463 freqconvbin(smpl_t f, uint_t samplerate, uint_t bufsize)
465 return aubio_freqtobin(f, samplerate, bufsize);
469 freqconvmidi (smpl_t f, uint_t samplerate UNUSED, uint_t bufsize UNUSED)
471 return aubio_freqtomidi (f);
475 freqconvpass (smpl_t f, uint_t samplerate UNUSED, uint_t bufsize UNUSED)
480 /* confidence callbacks */
482 aubio_pitch_get_confidence (aubio_pitch_t * p)
485 return p->conf_cb(p->p_object);