c1fdf7250fe76dcfed30b081083586cc56900943
[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 pf;
20
21    for (i=0;i<FRAME_SIZE;i++)
22       bak2[i]=0;
23    st = speex_encoder_init(&speex_wb_mode);
24    dec = speex_decoder_init(&speex_wb_mode);
25
26    pf=0;
27    speex_decoder_ctl(dec, SPEEX_SET_PF, &pf);
28    pf=10;
29    speex_encoder_ctl(st, SPEEX_SET_QUALITY, &pf);
30
31    if (argc != 4 && argc != 3)
32    {
33       fprintf (stderr, "Usage: encode [in file] [out file] [bits file]\nargc = %d", argc);
34       exit(1);
35    }
36    inFile = argv[1];
37    fin = fopen(inFile, "r");
38    outFile = argv[2];
39    fout = fopen(outFile, "w");
40    if (argc==4)
41    {
42       bitsFile = argv[3];
43       fbits = fopen(bitsFile, "w");
44    }
45    speex_bits_init(&bits);
46    while (!feof(fin))
47    {
48       fread(in, sizeof(short), FRAME_SIZE, fin);
49       if (feof(fin))
50          break;
51       for (i=0;i<FRAME_SIZE;i++)
52          bak[i]=input[i]=in[i];
53       speex_bits_reset(&bits);
54       speex_encode(st, input, &bits);
55       nbBits = speex_bits_write(&bits, cbits, 200);
56       printf ("Encoding frame in %d bits\n", nbBits*8);
57       if (argc==4)
58          fwrite(cbits, 1, nbBits, fbits);
59       {
60          float enoise=0, esig=0, snr;
61          for (i=0;i<FRAME_SIZE;i++)
62          {
63             enoise+=(bak2[i]-input[i])*(bak2[i]-input[i]);
64             esig += bak2[i]*bak2[i];
65          }
66          snr = 10*log10((esig+1)/(enoise+1));
67          printf ("real SNR = %f\n", snr);
68       }
69       speex_bits_rewind(&bits);
70       
71       speex_decode(dec, &bits, input, 0);
72
73       /* Save the bits here */
74       for (i=0;i<FRAME_SIZE;i++)
75       {
76          if (input[i]>32000)
77             input[i]=32000;
78          else if (input[i]<-32000)
79             input[i]=-32000;
80       }
81       speex_bits_reset(&bits);
82       for (i=0;i<FRAME_SIZE;i++)
83          in[i]=input[i];
84       for (i=0;i<FRAME_SIZE;i++)
85          bak2[i]=bak[i];
86       fwrite(in, sizeof(short), FRAME_SIZE, fout);
87    }
88    
89    speex_encoder_destroy(st);
90    speex_decoder_destroy(dec);
91    return 1;
92 }