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