ref: 7abfc82f117bda633cae850d6c7d57bc0c7d3067
dir: /examples/aubiomfcc.c/
/* Copyright (C) 2007 Amaury Hazan <ahazan@iua.upf.edu> and Paul Brossier <piem@piem.org> 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. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. */ #include "utils.h" /* mfcc objects */ fvec_t * mfcc_out; aubio_mfcc_t * mfcc; uint_t n_filters = 40; uint_t n_coefs = 11; unsigned int pos = 0; /*frames%dspblocksize*/ uint_t usepitch = 0; int aubio_process(float **input, float **output, int nframes); int aubio_process(float **input, float **output, int nframes) { unsigned int i; /*channels*/ unsigned int j; /*frames*/ for (j=0;j<(unsigned)nframes;j++) { if(usejack) { for (i=0;i<channels;i++) { /* write input to datanew */ fvec_write_sample(ibuf, input[i][j], i, pos); /* put synthnew in output */ output[i][j] = fvec_read_sample(obuf, i, pos); } } /*time for fft*/ if (pos == overlap_size-1) { /* block loop */ //compute mag spectrum aubio_pvoc_do (pv,ibuf, fftgrain); //compute mfccs aubio_mfcc_do(mfcc, fftgrain, mfcc_out); /* end of block loop */ pos = -1; /* so it will be zero next j loop */ } pos++; } return 1; } void process_print (void); void process_print (void) { /* output times in seconds write extracted mfccs */ uint_t coef_cnt; if (output_filename == NULL) { if(frames >= 4) { outmsg("%f\t",(frames-4)*overlap_size/(float)samplerate); } else if (frames < 4) { outmsg("%f\t",0.); } for (coef_cnt = 0; coef_cnt < n_coefs; coef_cnt++) { outmsg("%f ",mfcc_out->data[0][coef_cnt]); } outmsg("\n"); } } int main(int argc, char **argv) { // params buffer_size = 1024; overlap_size = 512; examples_common_init(argc,argv); smpl_t lowfreq = 133.333f; smpl_t highfreq = 44100.f; mfcc_out = new_fvec(n_coefs,channels); //populating the filter mfcc = new_aubio_mfcc(buffer_size, samplerate, n_filters, n_coefs , lowfreq, highfreq, channels); //process examples_common_process(aubio_process,process_print); //destroying mfcc del_aubio_mfcc(mfcc); del_fvec(mfcc_out); examples_common_del(); debug("End of program.\n"); fflush(stderr); return 0; }