libspeex/testenc_wb.c

   1 #ifdef HAVE_CONFIG_H
   2 #include "config.h"
   3 #endif
   4
   5 #include <speex/speex.h>
   6 #include <stdio.h>
   7 #include <stdlib.h>
   8 #include <speex/speex_callbacks.h>
   9
  10 #ifdef FIXED_DEBUG
  11 extern long long spx_mips;
  12 #endif
  13
  14 #define FRAME_SIZE 320
  15 #include <math.h>
  16 int main(int argc, char **argv)
  17 {
  18    char *inFile, *outFile, *bitsFile;
  19    FILE *fin, *fout, *fbits=NULL;
  20    short in_short[FRAME_SIZE];
  21    short out_short[FRAME_SIZE];
  22    float in_float[FRAME_SIZE];
  23    float sigpow,errpow,snr, seg_snr=0;
  24    int snr_frames = 0;
  25    char cbits[200];
  26    int nbBits;
  27    int i;
  28    void *st;
  29    void *dec;
  30    SpeexBits bits;
  31    int tmp;
  32    int bitCount=0;
  33    int skip_group_delay;
  34    SpeexCallback callback;
  35
  36    sigpow = 0;
  37    errpow = 0;
  38
  39    st = speex_encoder_init(&speex_wb_mode);
  40    dec = speex_decoder_init(&speex_wb_mode);
  41
  42    callback.callback_id = SPEEX_INBAND_CHAR;
  43    callback.func = speex_std_char_handler;
  44    callback.data = stderr;
  45    speex_decoder_ctl(dec, SPEEX_SET_HANDLER, &callback);
  46
  47    callback.callback_id = SPEEX_INBAND_MODE_REQUEST;
  48    callback.func = speex_std_mode_request_handler;
  49    callback.data = st;
  50    speex_decoder_ctl(dec, SPEEX_SET_HANDLER, &callback);
  51
  52    tmp=1;
  53    speex_decoder_ctl(dec, SPEEX_SET_ENH, &tmp);
  54    tmp=0;
  55    speex_encoder_ctl(st, SPEEX_SET_VBR, &tmp);
  56    tmp=8;
  57    speex_encoder_ctl(st, SPEEX_SET_QUALITY, &tmp);
  58    tmp=3;
  59    speex_encoder_ctl(st, SPEEX_SET_COMPLEXITY, &tmp);
  60    /*tmp=3;
  61    speex_encoder_ctl(st, SPEEX_SET_HIGH_MODE, &tmp);
  62    tmp=6;
  63    speex_encoder_ctl(st, SPEEX_SET_LOW_MODE, &tmp);
  64 */
  65
  66    speex_encoder_ctl(st, SPEEX_GET_LOOKAHEAD, &skip_group_delay);
  67    speex_decoder_ctl(dec, SPEEX_GET_LOOKAHEAD, &tmp);
  68    skip_group_delay += tmp;
  69
  70
  71    if (argc != 4 && argc != 3)
  72    {
  73       fprintf (stderr, "Usage: encode [in file] [out file] [bits file]\nargc = %d", argc);
  74       exit(1);
  75    }
  76    inFile = argv[1];
  77    fin = fopen(inFile, "r");
  78    outFile = argv[2];
  79    fout = fopen(outFile, "w+");
  80    if (argc==4)
  81    {
  82       bitsFile = argv[3];
  83       fbits = fopen(bitsFile, "w");
  84    }
  85    speex_bits_init(&bits);
  86    while (!feof(fin))
  87    {
  88       fread(in_short, sizeof(short), FRAME_SIZE, fin);
  89       if (feof(fin))
  90          break;
  91       for (i=0;i<FRAME_SIZE;i++)
  92          in_float[i]=in_short[i];
  93       speex_bits_reset(&bits);
  94
  95       speex_encode_int(st, in_short, &bits);
  96       nbBits = speex_bits_write(&bits, cbits, 200);
  97       bitCount+=bits.nbBits;
  98
  99       if (argc==4)
 100          fwrite(cbits, 1, nbBits, fbits);
 101       speex_bits_rewind(&bits);
 102
 103       speex_decode_int(dec, &bits, out_short);
 104       speex_bits_reset(&bits);
 105
 106       fwrite(&out_short[skip_group_delay], sizeof(short), FRAME_SIZE-skip_group_delay, fout);
 107       skip_group_delay = 0;
 108    }
 109    fprintf (stderr, "Total encoded size: %d bits\n", bitCount);
 110    speex_encoder_destroy(st);
 111    speex_decoder_destroy(dec);
 112    speex_bits_destroy(&bits);
 113
 114    rewind(fin);
 115    rewind(fout);
 116
 117    while ( FRAME_SIZE == fread(in_short, sizeof(short), FRAME_SIZE, fin)
 118            &&
 119            FRAME_SIZE ==  fread(out_short, sizeof(short), FRAME_SIZE,fout) )
 120    {
 121         float s=0, e=0;
 122         for (i=0;i<FRAME_SIZE;++i) {
 123             s += (float)in_short[i] * in_short[i];
 124             e += ((float)in_short[i]-out_short[i]) * ((float)in_short[i]-out_short[i]);
 125         }
 126         seg_snr += 10*log10((s+160)/(e+160));
 127         sigpow += s;
 128         errpow += e;
 129         snr_frames++;
 130    }
 131    fclose(fin);
 132    fclose(fout);
 133
 134    snr = 10 * log10( sigpow / errpow );
 135    seg_snr /= snr_frames;
 136    fprintf(stderr,"SNR = %f\nsegmental SNR = %f\n",snr, seg_snr);
 137
 138 #ifdef FIXED_DEBUG
 139    printf ("Total: %f MIPS\n", (float)(1e-6*50*spx_mips/snr_frames));
 140 #endif
 141
 142    return 1;
 143 }