source: examples/aubiomfcc.c @ 1b25a70

feature/autosinkfeature/constantqfeature/pitchshiftfeature/pydocstringsfeature/timestretchpitchshiftsamplertimestretchyinfft+
Last change on this file since 1b25a70 was 1b25a70, checked in by Paul Brossier <piem@piem.org>, 6 years ago

examples/: move parse_args to parse_args.h, clean up, remove lash and old frames_delay

  • Property mode set to 100644
File size: 2.4 KB
Line 
1/*
2  Copyright (C) 2007-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 "utils.h"
22#include "parse_args.h"
23
24/* mfcc objects */
25fvec_t * mfcc_out;
26aubio_mfcc_t * mfcc;
27aubio_pvoc_t *pv;
28cvec_t *fftgrain;
29
30uint_t n_filters = 40;
31uint_t n_coefs = 13;
32
33unsigned int pos = 0; /*frames%dspblocksize*/
34
35static int aubio_process(smpl_t **input, smpl_t **output, int nframes) {
36  unsigned int j;       /*frames*/
37 
38  for (j=0;j<(unsigned)nframes;j++) {
39    if(usejack) {
40      /* write input to datanew */
41      fvec_write_sample(ibuf, input[0][j], pos);
42      /* put synthnew in output */
43      output[0][j] = fvec_read_sample(obuf, pos);
44    }
45    /*time for fft*/
46    if (pos == overlap_size-1) {         
47      /* block loop */
48     
49      //compute mag spectrum
50      aubio_pvoc_do (pv, ibuf, fftgrain);
51     
52      //compute mfccs
53      aubio_mfcc_do(mfcc, fftgrain, mfcc_out);
54     
55      /* end of block loop */
56      pos = -1; /* so it will be zero next j loop */
57    }
58    pos++;
59  }
60  return 1;
61}
62
63static void process_print (void) {
64  /* output times in seconds and extracted mfccs */
65  if (sink_uri == NULL) {
66    outmsg("%f\t",frames*overlap_size/(float)samplerate);
67    fvec_print(mfcc_out);
68  }
69}
70
71int main(int argc, char **argv) {
72  // params
73  buffer_size  = 512;
74  overlap_size = 256;
75 
76  examples_common_init(argc,argv);
77
78  /* phase vocoder */
79  pv = new_aubio_pvoc (buffer_size, overlap_size);
80
81  fftgrain = new_cvec (buffer_size);
82
83  //populating the filter
84  mfcc = new_aubio_mfcc(buffer_size, n_filters, n_coefs, samplerate);
85 
86  mfcc_out = new_fvec(n_coefs);
87 
88  //process
89  examples_common_process(aubio_process,process_print);
90 
91  //destroying mfcc
92  del_aubio_pvoc (pv);
93  del_cvec (fftgrain);
94  del_aubio_mfcc(mfcc);
95  del_fvec(mfcc_out);
96
97  examples_common_del();
98  debug("End of program.\n");
99  fflush(stderr);
100 
101  return 0;
102}
103
Note: See TracBrowser for help on using the repository browser.