Improved spectral folding wideband mode
[speexdsp.git] / libspeex / testenc_wb.c
1 #include "speex.h"
2 #include <stdio.h>
3 #include <stdlib.h>
4
5 #define FRAME_SIZE 320
6 #include <math.h>
7 int main(int argc, char **argv)
8 {
9    char *inFile, *outFile, *bitsFile;
10    FILE *fin, *fout, *fbits=NULL;
11    short in[FRAME_SIZE];
12    float input[FRAME_SIZE], bak[FRAME_SIZE], bak2[FRAME_SIZE];
13    char cbits[200];
14    int nbBits;
15    int i;
16    void *st;
17    void *dec;
18    SpeexBits bits;
19    int tmp;
20    int bitCount=0;
21
22    for (i=0;i<FRAME_SIZE;i++)
23       bak2[i]=0;
24    st = speex_encoder_init(&speex_wb_mode);
25    dec = speex_decoder_init(&speex_wb_mode);
26    
27    tmp=0;
28    speex_decoder_ctl(dec, SPEEX_SET_ENH, &tmp);
29    tmp=0;
30    speex_encoder_ctl(st, SPEEX_SET_VBR, &tmp);
31    tmp=10;
32    speex_encoder_ctl(st, SPEEX_SET_QUALITY, &tmp);
33    tmp=3;
34    speex_encoder_ctl(st, SPEEX_SET_COMPLEXITY, &tmp);
35    tmp=3;
36    speex_encoder_ctl(st, SPEEX_SET_HIGH_MODE, &tmp);
37    tmp=4;
38    speex_encoder_ctl(st, SPEEX_SET_LOW_MODE, &tmp);
39    
40    if (argc != 4 && argc != 3)
41    { 
42       fprintf (stderr, "Usage: encode [in file] [out file] [bits file]\nargc = %d", argc);
43       exit(1);
44    }
45    inFile = argv[1];
46    fin = fopen(inFile, "r");
47    outFile = argv[2];
48    fout = fopen(outFile, "w");
49    if (argc==4)
50    {
51       bitsFile = argv[3];
52       fbits = fopen(bitsFile, "w");
53    }
54    speex_bits_init(&bits);
55    while (!feof(fin))
56    {
57       fread(in, sizeof(short), FRAME_SIZE, fin);
58       if (feof(fin))
59          break;
60       for (i=0;i<FRAME_SIZE;i++)
61          bak[i]=input[i]=in[i];
62       speex_bits_reset(&bits);
63       speex_encode(st, input, &bits);
64       nbBits = speex_bits_write(&bits, cbits, 200);
65       bitCount+=bits.nbBits;
66       printf ("Encoding frame in %d bits\n", nbBits*8);
67       if (argc==4)
68          fwrite(cbits, 1, nbBits, fbits);
69       {
70          float enoise=0, esig=0, snr;
71          for (i=0;i<FRAME_SIZE;i++)
72          {
73             enoise+=(bak2[i]-input[i])*(bak2[i]-input[i]);
74             esig += bak2[i]*bak2[i];
75          }
76          snr = 10*log10((esig+1)/(enoise+1));
77          printf ("real SNR = %f\n", snr);
78       }
79       speex_bits_rewind(&bits);
80       
81       speex_decode(dec, &bits, input, 0);
82
83       /* Save the bits here */
84       for (i=0;i<FRAME_SIZE;i++)
85       {
86          if (input[i]>32000)
87             input[i]=32000;
88          else if (input[i]<-32000)
89             input[i]=-32000;
90       }
91       speex_bits_reset(&bits);
92       for (i=0;i<FRAME_SIZE;i++)
93          in[i]=input[i];
94       for (i=0;i<FRAME_SIZE;i++)
95          bak2[i]=bak[i];
96       fwrite(in, sizeof(short), FRAME_SIZE, fout);
97    }
98    
99    fprintf (stderr, "Total encoded size: %d bits\n", bitCount);
100    speex_encoder_destroy(st);
101    speex_decoder_destroy(dec);
102    fclose(fin);
103    fclose(fout);
104    return 1;
105 }