source: examples/aubiomfcc.c @ 6a2478c

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

examples: more header updates to GPLv3

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