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