Implemented packet loss interpolation for wideband... not optimal yet...
[speexdsp.git] / libspeex / testenc.c
1 #include "speex.h"
2 #include <stdio.h>
3 #include <stdlib.h>
4
5 #define FRAME_SIZE 160
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_nb_mode);
25    dec = speex_decoder_init(&speex_nb_mode);
26
27    tmp=1;
28    speex_decoder_ctl(dec, SPEEX_SET_ENH, &tmp);
29    tmp=0;
30    speex_encoder_ctl(st, SPEEX_SET_VBR, &tmp);
31    tmp=6;
32    speex_encoder_ctl(st, SPEEX_SET_QUALITY, &tmp);
33    tmp=4;
34    speex_encoder_ctl(st, SPEEX_SET_COMPLEXITY, &tmp);
35
36    if (argc != 4 && argc != 3)
37    {
38       fprintf (stderr, "Usage: encode [in file] [out file] [bits file]\nargc = %d", argc);
39       exit(1);
40    }
41    inFile = argv[1];
42    fin = fopen(inFile, "r");
43    outFile = argv[2];
44    fout = fopen(outFile, "w");
45    if (argc==4)
46    {
47       bitsFile = argv[3];
48       fbits = fopen(bitsFile, "w");
49    }
50    speex_bits_init(&bits);
51    while (!feof(fin))
52    {
53       fread(in, sizeof(short), FRAME_SIZE, fin);
54       if (feof(fin))
55          break;
56       for (i=0;i<FRAME_SIZE;i++)
57          bak[i]=input[i]=in[i];
58       speex_bits_reset(&bits);
59       speex_encode(st, input, &bits);
60       nbBits = speex_bits_write(&bits, cbits, 200);
61       bitCount+=bits.nbBits;
62       printf ("Encoding frame in %d bits\n", nbBits*8);
63       if (argc==4)
64          fwrite(cbits, 1, nbBits, fbits);
65       {
66          float enoise=0, esig=0, snr;
67          for (i=0;i<FRAME_SIZE;i++)
68          {
69             enoise+=(bak2[i]-input[i])*(bak2[i]-input[i]);
70             esig += bak2[i]*bak2[i];
71          }
72          snr = 10*log10((esig+1)/(enoise+1));
73          printf ("real SNR = %f\n", snr);
74       }
75       speex_bits_rewind(&bits);
76       
77       speex_decode(dec, &bits, input, 0);
78       
79       /* Save the bits here */
80       for (i=0;i<FRAME_SIZE;i++)
81       {
82          if (input[i]>32000)
83             input[i]=32000;
84          else if (input[i]<-32000)
85             input[i]=-32000;
86       }
87       speex_bits_reset(&bits);
88       for (i=0;i<FRAME_SIZE;i++)
89          in[i]=input[i];
90       for (i=0;i<FRAME_SIZE;i++)
91          bak2[i]=bak[i];
92       fwrite(in, sizeof(short), FRAME_SIZE, fout);
93    }
94    fprintf (stderr, "Total encoded size: %d bits\n", bitCount);
95    
96    speex_encoder_destroy(st);
97    speex_decoder_destroy(dec);
98    return 1;
99 }