[svn] / branches / dev-api-4 / xvidcore / examples / xvid_stat.c Repository:
ViewVC logotype

Diff of /branches/dev-api-4/xvidcore/examples/xvid_stat.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 211, Fri Jun 14 15:36:22 2002 UTC revision 684, Wed Nov 27 21:09:10 2002 UTC
# Line 1  Line 1 
1  /**************************************************************************  /*****************************************************************************
2   *   *
3   *      XVID MPEG-4 VIDEO CODEC - Example for encoding and decoding   *  XVID MPEG-4 VIDEO CODEC
4     *  - Console based test application  -
5     *
6     *  Copyright(C) 2002 Christoph Lampert
7   *   *
8   *      This program is free software; you can redistribute it and/or modify   *      This program is free software; you can redistribute it and/or modify
9   *      it under the terms of the GNU General Public License as published by   *      it under the terms of the GNU General Public License as published by
# Line 14  Line 17 
17   *   *
18   *      You should have received a copy of the GNU General Public License   *      You should have received a copy of the GNU General Public License
19   *      along with this program; if not, write to the Free Software   *      along with this program; if not, write to the Free Software
20   *      Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.   *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
21   *   *
22   *************************************************************************/   * $Id: xvid_stat.c,v 1.16 2002-11-27 21:09:10 edgomez Exp $
   
 /************************************************************************  
23   *   *
24   *  PSNR and Speed test routine for XviD using the XviD-API   ****************************************************************************/
25   *  (C) Christoph Lampert, 2002/04/13  
26    /*****************************************************************************
27     *  Application notes :
28   *   *
29   *  A sequence of YUV pics in PGM file format is encoded and decoded   *  A sequence of YUV pics in PGM file format is encoded and decoded
30   *  The speed is measured and PSNR of decoded picture is calculated.   *  The speed is measured and PSNR of decoded picture is calculated.
31   *   *
32   *  The program is plain C and needs no libraries except for libxvidcore,   *  The program is plain C and needs no libraries except for libxvidcore,
33   *  and maths-lib ,so with UN*X you simply compile by   *  and maths-lib.
  *  
  *   gcc xvid_stat.c -lxvidcore -lm -o xvid_stat  
  *  
  *  Run without or with illegal parameters, then PGM input input is read  
  *  from stdin.  
  *  
  *  Parameters are: xvid_stat XDIM YDIM QUALITY BITRATE/QUANTIZER FRAMERATE  
  *  
  *  if XDIM or YDIM are illegal (e.g. 0), they are ignored and input is  
  *  considered to be PGM. Otherwise (X and Y both greater than 0) raw YUV  
  *  is expected, as e.g. the standard MPEG test-files, like "foreman"  
  *  
  *  0 <= QUALITY <= 6  (default 5)  
34   *   *
35   *  BITRATE is in kbps (default 900),   * Usage : xvid_stat [OPTIONS]
36   *      if BITRATE<32, then value is taken is fixed QUANTIZER   * Options :
37     *  -w integer     : frame width ([1.2048])
38     *  -h integer     : frame height ([1.2048])
39     *  -b integer     : target bitrate (>0 | default=900kbit)
40     *  -f float       : target framerate (>0)
41     *  -i string      : input filename (default=stdin)
42     *  -t integer     : input data type (yuv=0, pgm=1)
43     *  -n integer     : number of frames to encode
44     *  -q integer     : quality ([0..5])
45     *  -d boolean     : save decoder output (0 False*, !=0 True)
46     *  -m boolean     : save mpeg4 raw stream (0 False*, !=0 True)
47     *  -help          : prints this help message
48     *  -quant integer : fixed quantizer (disables -b setting)
49     *  (* means default)
50     *
51     *  An input file named "stdin" is treated as standard input stream.
52     *
53     *
54     *  PGM input must be in a very specific format, basically it pgm file must
55     *  contain Y plane first just like usual P5 pgm files, and then U and V
56     *  planes are stored just after Y plane so y dimension is y*3/2 in reality
57   *   *
58   *  FRAMERATE is a float (with or without decimal dot), default is 25.00   *  See read_pgmheader for more details.
59   *   *
60   *  input/output and m4v-output is saved, if corresponding flags are set   *  Such a PGM file can be generated from MPEG2 by # mpeg2dec -o pgmpipe
61   *   *
62   *  PGM input must in a very specific format, see read_pgmheader   ****************************************************************************/
  *  it can be generated e.g. from MPEG2 by    mpeg2dec -o pgmpipe  
  *  
  ************************************************************************/  
   
 /************************************************************************  
  *  
  *  For EXAMPLES how to use this, see the seperate file xvid_stat.examples  
  *  
  ************************************************************************/  
63    
64  #include <stdio.h>  #include <stdio.h>
65  #include <stdlib.h>  #include <stdlib.h>
66  #include <math.h>               // needed for log10  #include <string.h>
67  #include <sys/time.h>           // only needed for gettimeofday  #include <math.h>
68    #ifndef _MSC_VER
69    #include <sys/time.h>
70    #else
71    #include <time.h>
72    #endif
73    
74  #include "../src/xvid.h"                /* comes with XviD */  #include "xvid.h"
75    
76    /****************************************************************************
77     *                               Prototypes
78     ***************************************************************************/
79    
80    /* Prints program usage message */
81    static void usage();
82    
83    /* Statistical functions */
84    static double msecond();
85    static double absdistq(int x, int y,
86                                               unsigned char* buf1, int stride1,
87                                               unsigned char* buf2, int stride2);
88    static double PSNR(int x, int y,
89                                       unsigned char* buf1, int stride1,
90                                       unsigned char* buf2, int stride2);
91    
92    /* PGM related functions */
93    static int read_pgmheader(FILE* handle);
94    static int read_pgmdata(FILE* handle, unsigned char *image);
95    static int read_yuvdata(FILE* handle, unsigned char *image);
96    static int write_pgm(char *filename, unsigned char *image);
97    
98    /* Encoder related functions */
99    static int enc_init(int use_assembler);
100    static int enc_stop();
101    static int enc_main(unsigned char* image, unsigned char* bitstream,
102                                            int *streamlength, int* frametype);
103    
104    /* Decoder related functions */
105    static int dec_stop();
106    static int dec_main(unsigned char *m4v_buffer, unsigned char *out_buffer,
107                                            int m4v_size);
108    static int dec_init(int use_assembler);
109    
110    /*****************************************************************************
111     *                            Quality presets
112     ****************************************************************************/
113    
114    static int const motion_presets[7] = {
115            0,                                                        /* Q 0 */
116            PMV_EARLYSTOP16,                                          /* Q 1 */
117            PMV_EARLYSTOP16,                                          /* Q 2 */
118            PMV_EARLYSTOP16 | PMV_HALFPELREFINE16,                    /* Q 3 */
119            PMV_EARLYSTOP16 | PMV_HALFPELREFINE16,                    /* Q 4 */
120            PMV_EARLYSTOP16 | PMV_HALFPELREFINE16 | PMV_EARLYSTOP8 |  /* Q 5 */
121            PMV_HALFPELREFINE8,
122            PMV_EARLYSTOP16 | PMV_HALFPELREFINE16 | PMV_EXTSEARCH16 | /* Q 6 */
123            PMV_USESQUARES16 | PMV_EARLYSTOP8 | PMV_HALFPELREFINE8
124    };
125    
126  int motion_presets[7] = {  static int const general_presets[7] = {
127          0,                                                              // Q 0          XVID_H263QUANT,                               /* Q 0 */
128          PMV_EARLYSTOP16,                                                // Q 1          XVID_MPEGQUANT,                               /* Q 1 */
129          PMV_EARLYSTOP16,                                                // Q 2          XVID_H263QUANT,                               /* Q 2 */
130          PMV_EARLYSTOP16 | PMV_HALFPELREFINE16,                          // Q 3          XVID_H263QUANT | XVID_HALFPEL,                /* Q 3 */
131          PMV_EARLYSTOP16 | PMV_HALFPELREFINE16,                          // Q 4          XVID_H263QUANT | XVID_HALFPEL | XVID_INTER4V, /* Q 4 */
132          PMV_EARLYSTOP16 | PMV_HALFPELREFINE16 | PMV_EARLYSTOP8          // Q 5          XVID_H263QUANT | XVID_HALFPEL | XVID_INTER4V, /* Q 5 */
133                          | PMV_HALFPELREFINE8,          XVID_H263QUANT | XVID_HALFPEL | XVID_INTER4V  /* Q 6 */
         PMV_EARLYSTOP16 | PMV_HALFPELREFINE16 | PMV_EXTSEARCH16         // Q 6  
                         | PMV_USESQUARES16 | PMV_EARLYSTOP8 | PMV_HALFPELREFINE8  
134          };          };
135    
 int general_presets[7] = {  
         XVID_H263QUANT, /* or use XVID_MPEGQUANT */             // Q 0  
         XVID_MPEGQUANT,                                         // Q 1  
         XVID_H263QUANT,                                // Q 2  
         XVID_H263QUANT | XVID_HALFPEL,                          // Q 3  
         XVID_H263QUANT | XVID_HALFPEL | XVID_INTER4V,           // Q 4  
         XVID_H263QUANT | XVID_HALFPEL | XVID_INTER4V,           // Q 5  
         XVID_H263QUANT | XVID_HALFPEL | XVID_INTER4V };         // Q 6  
136    
137    /*****************************************************************************
138     *                     Command line global variables
139     ****************************************************************************/
140    
141    /* Maximum number of frames to encode */
142    #define ABS_MAXFRAMENR 9999
143    
144    static int   ARG_BITRATE = 900;
145    static int   ARG_QUANTI = 0;
146    static int   ARG_QUALITY = 6;
147    static int   ARG_MINQUANT = 1;
148    static int   ARG_MAXQUANT = 31;
149    static float ARG_FRAMERATE = 25.00f;
150    static int   ARG_MAXFRAMENR = ABS_MAXFRAMENR;
151    static char *ARG_INPUTFILE = NULL;
152    static int   ARG_INPUTTYPE = 0;
153    static int   ARG_SAVEDECOUTPUT = 0;
154    static int   ARG_SAVEMPEGSTREAM = 0;
155    static int   XDIM = 0;
156    static int   YDIM = 0;
157    #define IMAGE_SIZE(x,y) ((x)*(y)*3/2)
158    
159  /* my default values for encoding */  #define MAX(A,B) ( ((A)>(B)) ? (A) : (B) )
160    #define SMALL_EPS 1e-10
161    
162  #define ABS_MAXFRAMENR 9999               // max number of frames  /****************************************************************************
163     *                     Nasty global vars ;-)
164     ***************************************************************************/
165    
166    static int i,filenr = 0;
167    static int save_ref_flag = 0;
168    
169    /* the path where to save output */
170    static char filepath[256] = "./";
171    
172    /* Internal structures (handles) for encoding and decoding */
173    static void *enc_handle = NULL;
174    static void *dec_handle = NULL;
175    
176    /*****************************************************************************
177     *                               Main program
178     ****************************************************************************/
179    
180  int ARG_BITRATE=900;  int main(int argc, char *argv[])
181  int ARG_QUANTI=0;  {
182    
183  int ARG_QUALITY =6;          unsigned char *divx_buffer = NULL;
184  int ARG_MINQUANT=1;          unsigned char *in_buffer = NULL;
185  int ARG_MAXQUANT=31;          unsigned char *out_buffer = NULL;
 float ARG_FRAMERATE=25.00;  
186    
187  int ARG_MAXFRAMENR=ABS_MAXFRAMENR;          double enctime,dectime;
188            double totalenctime=0.;
189            double totaldectime=0.;
190    
191            long totalsize=0;
192            int status;
193    
194  #define MAX(A,B) ( ((A)>(B)) ? (A) : (B) )          int m4v_size;
195  #define SMALL_EPS 1e-10          int frame_type[ABS_MAXFRAMENR];
196            int Iframes=0, Pframes=0, use_assembler=0;
197            double framepsnr[ABS_MAXFRAMENR];
198    
199            double Ipsnr=0.,Imaxpsnr=0.,Iminpsnr=999.,Ivarpsnr=0.;
200            double Ppsnr=0.,Pmaxpsnr=0.,Pminpsnr=999.,Pvarpsnr=0.;
201    
202            char filename[256];
203    
204            FILE *filehandle;
205            FILE *in_file = stdin;
206    
207            printf("xvid_stat - XviD core library test program ");
208            printf("written by Christoph Lampert 2002\n\n");
209    
210    /*****************************************************************************
211     *                            Command line parsing
212     ****************************************************************************/
213    
214            for (i=1; i< argc; i++) {
215    
216                    if (strcmp("-asm", argv[i]) == 0 ) {
217                            use_assembler = 1;
218                    }
219                    else if (strcmp("-w", argv[i]) == 0 && i < argc - 1 ) {
220                            i++;
221                            XDIM = atoi(argv[i]);
222                    }
223                    else if (strcmp("-h", argv[i]) == 0 && i < argc - 1 ) {
224                            i++;
225                            YDIM = atoi(argv[i]);
226                    }
227                    else if (strcmp("-b", argv[i]) == 0 && i < argc - 1 ) {
228                            i++;
229                            ARG_BITRATE = atoi(argv[i]);
230                    }
231                    else if (strcmp("-q", argv[i]) == 0 && i < argc - 1 ) {
232                            i++;
233                            ARG_QUALITY = atoi(argv[i]);
234                    }
235                    else if (strcmp("-f", argv[i]) == 0 && i < argc - 1 ) {
236                            i++;
237                            ARG_FRAMERATE = (float)atof(argv[i]);
238                    }
239                    else if (strcmp("-i", argv[i]) == 0 && i < argc - 1 ) {
240                            i++;
241                            ARG_INPUTFILE = argv[i];
242                    }
243                    else if (strcmp("-t", argv[i]) == 0 && i < argc - 1 ) {
244                            i++;
245                            ARG_INPUTTYPE = atoi(argv[i]);
246                    }
247                    else if(strcmp("-n", argv[i]) == 0 && i < argc - 1 ) {
248                            i++;
249                            ARG_MAXFRAMENR  = atoi(argv[i]);
250                    }
251                    else if (strcmp("-quant", argv[i]) == 0 && i < argc - 1 ) {
252                            i++;
253                            ARG_QUANTI = atoi(argv[i]);
254                    }
255                    else if (strcmp("-d", argv[i]) == 0 && i < argc - 1 ) {
256                            i++;
257                            ARG_SAVEDECOUTPUT = atoi(argv[i]);
258                    }
259                    else if (strcmp("-m", argv[i]) == 0 && i < argc - 1 ) {
260                            i++;
261                            ARG_SAVEMPEGSTREAM = atoi(argv[i]);
262                    }
263                    else if (strcmp("-help", argv[i])) {
264                            usage();
265                            return(0);
266                    }
267                    else {
268                            usage();
269                            exit(-1);
270                    }
271    
272            }
273    
274    /*****************************************************************************
275     *                            Arguments checking
276     ****************************************************************************/
277    
278            if (XDIM <= 0 || XDIM >= 2048 || YDIM <=0 || YDIM >= 2048 ) {
279                    fprintf(stderr, "Trying to retreive width and height from PGM header\n");
280                    ARG_INPUTTYPE = 1; /* pgm */
281            }
282    
283            if ( ARG_QUALITY < 0 || ARG_QUALITY > 6) {
284                    fprintf(stderr,"Wrong Quality\n");
285                    return -1;
286            }
287    
288            if ( ARG_BITRATE <= 0 && ARG_QUANTI == 0) {
289                    fprintf(stderr,"Wrong Bitrate\n");
290                    return -1;
291            }
292    
293            if ( ARG_FRAMERATE <= 0) {
294                    fprintf(stderr,"Wrong Framerate %s \n",argv[5]);
295                    return -1;
296            }
297    
298            if ( ARG_MAXFRAMENR <= 0) {
299                    fprintf(stderr,"Wrong number of frames\n");
300                    return -1;
301            }
302    
303            if ( ARG_INPUTFILE == NULL || strcmp(ARG_INPUTFILE, "stdin") == 0) {
304                    in_file = stdin;
305            }
306            else {
307    
308                    in_file = fopen(ARG_INPUTFILE, "rb");
309                    if (in_file == NULL) {
310                            fprintf(stderr, "Error opening input file %s\n", ARG_INPUTFILE);
311                            return -1;
312                    }
313            }
314    
315            if (ARG_INPUTTYPE) {
316                    if (read_pgmheader(in_file)) {
317                            fprintf(stderr, "Wrong input format, I want YUV encapsulated in PGM\n");
318                            return -1;
319                    }
320            }
321    
322            /* now we know the sizes, so allocate memory */
323    
324  /* these are global variables. Not very elegant, but easy, and this is an easy program */          in_buffer = (unsigned char *) malloc(IMAGE_SIZE(XDIM,YDIM));
325            if (!in_buffer)
326                    goto free_all_memory;
327    
328  int XDIM=0;          /* this should really be enough memory ! */
329  int YDIM=0;     // will be set when reading first image          divx_buffer = (unsigned char *) malloc(IMAGE_SIZE(XDIM,YDIM)*2);
330  int i,filenr = 0;          if (!divx_buffer)
331                    goto free_all_memory;
332    
333  int save_m4v_flag = 0;          // save MPEG4-bytestream?          out_buffer = (unsigned char *) malloc(IMAGE_SIZE(XDIM,YDIM)*4);
334  int save_dec_flag = 1;          // save decompressed bytestream?          if (!out_buffer)
335  int save_ref_flag = 0;          //                  goto free_all_memory;
336    
 int pgmflag = 0;                // a flag, if input is in PGM format, overwritten in init-phase  
 char filepath[256] = "./";      // the path where to save output  
337    
338  void *enc_handle = NULL;                // internal structures (handles) for encoding  /*****************************************************************************
339  void *dec_handle = NULL;                // and decoding   *                            XviD PART  Start
340     ****************************************************************************/
341    
342    
343  /*********************************************************************/          status = enc_init(use_assembler);
344  /*                     "statistical" functions                       */          if (status)
345  /*                                                                   */          {
346  /*  these are not needed for encoding or decoding, but for measuring */                  fprintf(stderr, "Encore INIT problem, return value %d\n", status);
347  /*  time and quality, there in nothing specific to XviD in these     */                  goto release_all;
348  /*                                                                   */          }
 /*********************************************************************/  
349    
350  double msecond()          status = dec_init(use_assembler);
351  /* return the current time in seconds(!)  */          if (status)
352  {  {
353                    fprintf(stderr, "Decore INIT problem, return value %d\n", status);
354                    goto release_all;
355            }
356    
357    
358    /*****************************************************************************
359     *                            Main loop
360     ****************************************************************************/
361    
362            do {
363    
364                    if (ARG_INPUTTYPE)
365                            status = read_pgmdata(in_file, in_buffer);      /* read PGM data (YUV-format) */
366                    else
367                            status = read_yuvdata(in_file, in_buffer);      /* read raw data (YUV-format) */
368    
369                    if (status)
370                    {
371                            /* Couldn't read image, most likely end-of-file */
372                            continue;
373                    }
374    
375    
376                    if (save_ref_flag)
377                    {
378                            sprintf(filename, "%s%05d.pgm", filepath, filenr);
379                            write_pgm(filename,in_buffer);
380                    }
381    
382    
383    /*****************************************************************************
384     *                       Analyse this frame before encoding
385     ****************************************************************************/
386    
387    /*
388     *      nothing is done here at the moment, but you could e.g. create
389     *      histograms or measure entropy or apply preprocessing filters...
390     */
391    
392    /*****************************************************************************
393     *                       Encode and decode this frame
394     ****************************************************************************/
395    
396                    enctime = msecond();
397                    status = enc_main(in_buffer, divx_buffer, &m4v_size, &frame_type[filenr]);
398                    enctime = msecond() - enctime;
399    
400                    totalenctime += enctime;
401                    totalsize += m4v_size;
402    
403                    printf("Frame %5d: intra %1d, enctime=%6.1f ms, size=%6d bytes ",
404                               (int)filenr, (int)frame_type[filenr], (float)enctime, (int)m4v_size);
405    
406                    if (ARG_SAVEMPEGSTREAM)
407                    {
408                            sprintf(filename, "%sframe%05d.m4v", filepath, filenr);
409                            filehandle = fopen(filename, "wb");
410                            fwrite(divx_buffer, m4v_size, 1, filehandle);
411                            fclose(filehandle);
412                    }
413    
414                    dectime = msecond();
415                    status = dec_main(divx_buffer, out_buffer, m4v_size);
416                    dectime = msecond() - dectime;
417    
418                    totaldectime += dectime;
419    
420    
421    /*****************************************************************************
422     *             Analyse the decoded frame and compare to original
423     ****************************************************************************/
424    
425                    framepsnr[filenr] = PSNR(XDIM,YDIM*3/2, in_buffer, XDIM, out_buffer, XDIM);
426    
427                    printf("dectime =%6.1f ms PSNR %5.2f\n",dectime, framepsnr[filenr]);
428    
429                    if (ARG_SAVEDECOUTPUT)
430                    {
431                            sprintf(filename, "%sdec%05d.pgm", filepath, filenr);
432                            write_pgm(filename, out_buffer);
433                    }
434    
435                    /* Read the header if it's pgm stream */
436                    if (ARG_INPUTTYPE)
437                            status = read_pgmheader(in_file);
438    
439                    filenr++;
440    
441            } while ( (!status) && (filenr<ARG_MAXFRAMENR) );
442    
443    
444    
445    /*****************************************************************************
446     *         Calculate totals and averages for output, print results
447     ****************************************************************************/
448    
449            totalsize    /= filenr;
450            totalenctime /= filenr;
451            totaldectime /= filenr;
452    
453            for (i=0;i<filenr;i++)
454            {
455                    switch (frame_type[i])
456                    {
457                    case 0:
458                            Pframes++;
459                            Ppsnr += framepsnr[i];
460                            break;
461                    case 1:
462                            Iframes++;
463                            Ipsnr += framepsnr[i];
464                            break;
465                    default:
466                            break;
467                    }
468            }
469    
470            if (Pframes)
471                    Ppsnr /= Pframes;
472            if (Iframes)
473                    Ipsnr /= Iframes;
474    
475            /* calculate statistics for every frametype: P,I */
476            for (i=0;i<filenr;i++)
477            {
478                    switch (frame_type[i])
479                    {
480                    case 0:
481                            if (framepsnr[i] > Pmaxpsnr)
482                                    Pmaxpsnr = framepsnr[i];
483                            if (framepsnr[i] < Pminpsnr)
484                                    Pminpsnr = framepsnr[i];
485                            Pvarpsnr += (framepsnr[i] - Ppsnr)*(framepsnr[i] - Ppsnr) /Pframes;
486                            break;
487                    case 1:
488                            if (framepsnr[i] > Imaxpsnr)
489                                    Imaxpsnr = framepsnr[i];
490                            if (framepsnr[i] < Pminpsnr)
491                                    Iminpsnr = framepsnr[i];
492                            Ivarpsnr += (framepsnr[i] - Ipsnr)*(framepsnr[i] - Ipsnr) /Iframes;
493                    default:
494                            break;
495                    }
496            }
497    
498            /* Print all statistics */
499            printf("Avg. Q%1d %2s ",ARG_QUALITY, (ARG_QUANTI ? " q" : "br"));
500            printf("%04d ",(ARG_QUANTI)?ARG_QUANTI:ARG_BITRATE);
501            printf("( %.2f bpp) ", (double)ARG_BITRATE*1000/XDIM/YDIM/ARG_FRAMERATE);
502            printf("size %6d ", (int)totalsize);
503            printf("( %4d kbps ",(int)(totalsize*8*ARG_FRAMERATE/1000));
504            printf("/ %.2f bpp) ",(double)totalsize*8/XDIM/YDIM);
505            printf("enc: %6.1f fps, dec: %6.1f fps \n",1000/totalenctime, 1000/totaldectime);
506            printf("PSNR P(%d): %5.2f ( %5.2f , %5.2f ; %5.4f ) ",Pframes,Ppsnr,Pminpsnr,Pmaxpsnr,sqrt(Pvarpsnr/filenr));
507            printf("I(%d): %5.2f ( %5.2f , %5.2f ; %5.4f ) ",Iframes,Ipsnr,Iminpsnr,Imaxpsnr,sqrt(Ivarpsnr/filenr));
508            printf("\n");
509    
510    /*****************************************************************************
511     *                            XviD PART  Stop
512     ****************************************************************************/
513    
514     release_all:
515    
516            if (enc_handle)
517            {
518                    status = enc_stop();
519                    if (status)
520                            fprintf(stderr, "Encore RELEASE problem return value %d\n", status);
521            }
522    
523            if (dec_handle)
524            {
525                    status = dec_stop();
526                    if (status)
527                            fprintf(stderr, "Decore RELEASE problem return value %d\n", status);
528            }
529    
530            fclose(in_file);
531    
532     free_all_memory:
533            free(out_buffer);
534            free(divx_buffer);
535            free(in_buffer);
536    
537            return 0;
538    
539    }
540    
541    /*****************************************************************************
542     *                        "statistical" functions
543     *
544     *  these are not needed for encoding or decoding, but for measuring
545     *  time and quality, there in nothing specific to XviD in these
546     *
547     *****************************************************************************/
548    
549    
550    
551    /* Return time elapsed time in miliseconds since the program started */
552    static double msecond()
553    {
554    #ifndef _MSC_VER
555          struct timeval  tv;          struct timeval  tv;
556          gettimeofday(&tv, 0);          gettimeofday(&tv, 0);
557          return tv.tv_sec + tv.tv_usec * 1.0e-6;          return tv.tv_sec*1.0e3 + tv.tv_usec * 1.0e-3;
558    #else
559            clock_t clk;
560            clk = clock();
561            return clk * 1000 / CLOCKS_PER_SEC;
562    #endif
563  }  }
564    
565    
566    /*
567  double absdistq(int x,int y, unsigned char* buf1, int stride1, unsigned char* buf2, int stride2)   * Returns the sum of squared distances (SSD) between two images of dimensions
568  /* returns the sum of squared distances (SSD) between two images of dimensions x times y */   * x times y
569     */
570    static double absdistq(int x, int y,
571                                               unsigned char* buf1, int stride1,
572                                               unsigned char* buf2, int stride2)
573  {  {
574          double dist=0.;          double dist=0.;
575          int i,j,val;          int i,j,val;
# Line 163  Line 588 
588  }  }
589    
590    
591  double PSNR(int x,int y, unsigned char* buf1, int stride1, unsigned char* buf2, int stride2 )  /*
592  /* return the PSNR between to images                                               */   * Returns the PSNR between to images.
593  /* this is a logarithmic measure for "quality" from the world of signal processing */   *
594  /* if you don't know what it is, simply accept that higher values are better       */   * This is a common logarithmic measure for "quality" from the world of signal
595     * processing if you don't know what it is, simply accept that higher values
596     * are better.
597     *
598     * PSNR represents the ratio of useful signal over noise signal. In our case,
599     * useful signal is refernce image, noise signal is the difference between
600     * reference and decoded frame from encoded bitstream.
601     *
602     * The problem is this type of value is dependant of image source and so, is
603     * not reliable as a common "quality" indicator.
604     * So PSNR computes the ratio of maximum/noise. Maximum being set to 2^bpp/channel
605     * This way, PSNR is not dependant anymore of image data type.
606     *
607     */
608    static double PSNR(int x, int y,
609                                       unsigned char* buf1, int stride1,
610                                       unsigned char* buf2, int stride2)
611  {  {
612     return 10*(log10(255*255)-log10( absdistq(x, y, buf1, stride1, buf2, stride2) ));     return 10*(log10(255*255)-log10( absdistq(x, y, buf1, stride1, buf2, stride2) ));
613  }  }
614    
615    /*****************************************************************************
616     *                             Usage message
617     *****************************************************************************/
618    
619    static void usage()
620    {
621    
622  /*********************************************************************/          fprintf(stderr, "Usage : xvid_stat [OPTIONS]\n");
623  /*                    input and output functions                     */          fprintf(stderr, "Options :\n");
624  /*                                                                   */          fprintf(stderr, " -w integer     : frame width ([1.2048])\n");
625  /* the are small and simple routines to read and write PGM and YUV   */          fprintf(stderr, " -h integer     : frame height ([1.2048])\n");
626  /* image. It's just for convenience, again nothing specific to XviD  */          fprintf(stderr, " -b integer     : target bitrate (>0 | default=900kbit)\n");
627  /*                                                                   */          fprintf(stderr, " -f float       : target framerate (>0)\n");
628  /*********************************************************************/          fprintf(stderr, " -i string      : input filename (default=stdin)\n");
629            fprintf(stderr, " -t integer     : input data type (yuv=0, pgm=1)\n");
630            fprintf(stderr, " -n integer     : number of frames to encode\n");
631            fprintf(stderr, " -q integer     : quality ([0..5])\n");
632            fprintf(stderr, " -d boolean     : save decoder output (0 False*, !=0 True)\n");
633            fprintf(stderr, " -m boolean     : save mpeg4 raw stream (0 False*, !=0 True)\n");
634            fprintf(stderr, " -help          : prints this help message\n");
635            fprintf(stderr, " -quant integer : fixed quantizer (disables -b setting)\n");
636            fprintf(stderr, " (* means default)\n");
637    
638  int read_pgmheader(FILE* handle)  }
639    
640    /*****************************************************************************
641     *                       Input and output functions
642     *
643     *      the are small and simple routines to read and write PGM and YUV
644     *      image. It's just for convenience, again nothing specific to XviD
645     *
646     *****************************************************************************/
647    
648    static int read_pgmheader(FILE* handle)
649  {  {
650          int bytes,xsize,ysize,depth;          int bytes,xsize,ysize,depth;
651          char dummy[2];          char dummy[2];
# Line 189  Line 654 
654    
655          if ( (bytes < 2) || (dummy[0] != 'P') || (dummy[1] != '5' ))          if ( (bytes < 2) || (dummy[0] != 'P') || (dummy[1] != '5' ))
656                  return 1;                  return 1;
657    
658          fscanf(handle,"%d %d %d",&xsize,&ysize,&depth);          fscanf(handle,"%d %d %d",&xsize,&ysize,&depth);
659          if ( (xsize > 1440) || (ysize > 2880 ) || (depth != 255) )          if ( (xsize > 1440) || (ysize > 2880 ) || (depth != 255) )
660          {          {
# Line 196  Line 662 
662                  return 2;                  return 2;
663          }          }
664          if ( (XDIM==0) || (YDIM==0) )          if ( (XDIM==0) || (YDIM==0) )
665          {       XDIM=xsize;          {
666                  YDIM=ysize;                  XDIM=xsize;
667                    YDIM=ysize*2/3;
668          }          }
669    
670          return 0;          return 0;
671  }  }
672    
673  int read_pgmdata(FILE* handle, unsigned char *image)  static int read_pgmdata(FILE* handle, unsigned char *image)
674  {  {
675          int i,status;          int i;
676          char dummy;          char dummy;
677    
678          unsigned char* buff1_ptr2 = image + XDIM*YDIM;          unsigned char *y = image;
679          unsigned char* buff1_ptr3 = image + XDIM*YDIM + XDIM/2*YDIM/2;          unsigned char *u = image + XDIM*YDIM;
680            unsigned char *v = image + XDIM*YDIM + XDIM/2*YDIM/2;
681    
682          fread(image,XDIM*YDIM,1,stdin); // read Y component of picture          /* read Y component of picture */
683            fread(y, 1, XDIM*YDIM, handle);
684    
685          for (i=0;i<YDIM/2;i++)          for (i=0;i<YDIM/2;i++)
686          {          {
687                  fread(buff1_ptr2,XDIM/2,1,stdin);        // read U                  /* read U */
688                  buff1_ptr2 += XDIM/2;                  fread(u, 1, XDIM/2, handle);
689                  fread(buff1_ptr3,XDIM/2,1,stdin);        // read V  
690                  buff1_ptr3 += XDIM/2;                  /* read V */
691                    fread(v, 1, XDIM/2, handle);
692    
693                    /* Update pointers */
694                    u += XDIM/2;
695                    v += XDIM/2;
696          }          }
697          fread(&dummy,1,1,handle);       //  I don't know why, but this seems needed  
698        /*  I don't know why, but this seems needed */
699            fread(&dummy, 1, 1, handle);
700    
701          return 0;          return 0;
702  }  }
703    
704  int read_yuvdata(FILE* handle, unsigned char *image)  static int read_yuvdata(FILE* handle, unsigned char *image)
705  {       int i;  {
         char dummy;  
   
         unsigned char* buff1_ptr2 = image + XDIM*YDIM;  
         unsigned char* buff1_ptr3 = image + XDIM*YDIM + XDIM/2*YDIM/2;  
706    
707          if (fread(image,XDIM,YDIM*3/2,stdin) != YDIM*3/2)          if (fread(image, 1, IMAGE_SIZE(XDIM, YDIM), handle) != (unsigned int)IMAGE_SIZE(XDIM, YDIM))
708                  return 1;                  return 1;
709          else          else
710                  return 0;                  return 0;
711  }  }
712    
713  int write_pgm(char *filename, unsigned char *image)  static int write_pgm(char *filename, unsigned char *image)
714  {  {
715            int loop;
716    
717            unsigned char *y = image;
718            unsigned char *u = image + XDIM*YDIM;
719            unsigned char *v = image + XDIM*YDIM + XDIM/2*YDIM/2;
720    
721          FILE *filehandle;          FILE *filehandle;
722          filehandle=fopen(filename,"wb");          filehandle=fopen(filename,"w+b");
723          if (filehandle)          if (filehandle)
724          {          {
725                  fprintf(filehandle,"P5\n\n");           //                  /* Write header */
726                  fprintf(filehandle,"%d %d 255\n",XDIM,YDIM*3/2);                  fprintf(filehandle,"P5\n\n%d %d 255\n", XDIM,YDIM*3/2);
727                  fwrite(image,XDIM,YDIM*3/2,filehandle);  
728                    /* Write Y data */
729                    fwrite(y, 1, XDIM*YDIM, filehandle);
730    
731                    for(loop=0; loop<YDIM/2; loop++)
732                    {
733                            /* Write U scanline */
734                            fwrite(u, 1, XDIM/2, filehandle);
735    
736                            /* Write V scanline */
737                            fwrite(v, 1, XDIM/2, filehandle);
738    
739                            /* Update pointers */
740                            u += XDIM/2;
741                            v += XDIM/2;
742    
743                    }
744    
745                    /* Close file */
746                  fclose(filehandle);                  fclose(filehandle);
747    
748                  return 0;                  return 0;
749          }          }
750          else          else
751                  return 1;                  return 1;
752  }  }
753    
754    /*****************************************************************************
755     *     Routines for encoding: init encoder, frame step, release encoder
756  /*********************************************************************/   ****************************************************************************/
 /* Routines for encoding: init encoder, frame step, release encoder  */  
 /*********************************************************************/  
757    
758  #define FRAMERATE_INCR 1001  #define FRAMERATE_INCR 1001
759    
760    /* Initialize encoder for first use, pass all needed parameters to the codec */
761  int enc_init(int use_assembler)  static int enc_init(int use_assembler)
762  {       /* initialize encoder for first use, pass all needed parameters to the codec */  {
763          int xerr;          int xerr;
764    
765          XVID_INIT_PARAM xinit;          XVID_INIT_PARAM xinit;
766          XVID_ENC_PARAM xparam;          XVID_ENC_PARAM xparam;
767    
768          if(use_assembler)          if(use_assembler) {
769    
770  #ifdef ARCH_IA64  #ifdef ARCH_IA64
771                  xinit.cpu_flags = XVID_CPU_FORCE | XVID_CPU_IA64;                  xinit.cpu_flags = XVID_CPU_FORCE | XVID_CPU_IA64;
772  #else  #else
773                  xinit.cpu_flags = 0;                  xinit.cpu_flags = 0;
774  #endif  #endif
775            }
776          else          else {
777                  xinit.cpu_flags = XVID_CPU_FORCE;                  xinit.cpu_flags = XVID_CPU_FORCE;
778            }
779    
780          xvid_init(NULL, 0, &xinit, NULL);          xvid_init(NULL, 0, &xinit, NULL);
781    
# Line 298  Line 795 
795          xparam.rc_averaging_period = 100;          xparam.rc_averaging_period = 100;
796          xparam.rc_buffer = 10;          xparam.rc_buffer = 10;
797          xparam.rc_bitrate = ARG_BITRATE*1000;          xparam.rc_bitrate = ARG_BITRATE*1000;
798          xparam.min_quantizer = 1;          xparam.min_quantizer = ARG_MINQUANT;
799          xparam.max_quantizer = 31;          xparam.max_quantizer = ARG_MAXQUANT;
800          xparam.max_key_interval = (int)ARG_FRAMERATE*10;          xparam.max_key_interval = (int)ARG_FRAMERATE*10;
801    
802                  /* I use a small value here, since will not encode whole movies, but short clips */                  /* I use a small value here, since will not encode whole movies, but short clips */
# Line 310  Line 807 
807          return xerr;          return xerr;
808  }  }
809    
810  int  enc_stop()  static int enc_stop()
811  {       int xerr;  {
812            int xerr;
813    
814          xerr = xvid_encore(enc_handle, XVID_ENC_DESTROY, NULL, NULL);          xerr = xvid_encore(enc_handle, XVID_ENC_DESTROY, NULL, NULL);
815          return xerr;          return xerr;
816    
817  }  }
818    
819  int  enc_main(unsigned char* image, unsigned char* bitstream, int *streamlength, int* frametype)  static int enc_main(unsigned char* image, unsigned char* bitstream,
820  {       int xerr;                                          int *streamlength, int* frametype)
821    {
822            int xerr;
823    
824          XVID_ENC_FRAME xframe;          XVID_ENC_FRAME xframe;
825          XVID_ENC_STATS xstats;          XVID_ENC_STATS xstats;
826    
827          xframe.bitstream = bitstream;          xframe.bitstream = bitstream;
828          xframe.length = -1;     // this is written by the routine          xframe.length = -1;     /* this is written by the routine */
829    
830          xframe.image = image;          xframe.image = image;
831          xframe.colorspace = XVID_CSP_YV12;      // defined in <xvid.h>          xframe.colorspace = XVID_CSP_YV12;      /* defined in <xvid.h> */
832    
833          xframe.intra = -1; // let the codec decide between I-frame (1) and P-frame (0)          xframe.intra = -1; /* let the codec decide between I-frame (1) and P-frame (0) */
834    
835          xframe.quant = ARG_QUANTI;      // is quant != 0, use a fixed quant (and ignore bitrate)          xframe.quant = ARG_QUANTI;      /* is quant != 0, use a fixed quant (and ignore bitrate) */
836    
837          xframe.motion = motion_presets[ARG_QUALITY];          xframe.motion = motion_presets[ARG_QUALITY];
838          xframe.general = general_presets[ARG_QUALITY];          xframe.general = general_presets[ARG_QUALITY];
# Line 339  Line 840 
840    
841          xerr = xvid_encore(enc_handle, XVID_ENC_ENCODE, &xframe, &xstats);          xerr = xvid_encore(enc_handle, XVID_ENC_ENCODE, &xframe, &xstats);
842    
843  /*              enc_result->is_key_frame = xframe.intra;          /*
844                  enc_result->quantizer = xframe.quant;           * This is statictical data, e.g. for 2-pass. If you are not
845                  enc_result->total_bits = xframe.length * 8;           * interested in any of this, you can use NULL instead of &xstats
                 enc_result->motion_bits = xstats.hlength * 8;  
                 enc_result->texture_bits = enc_result->total_bits - enc_result->motion_bits;  
 */  
   
 /*  This is statictical data, e.g. for 2-pass.  
     If you are not interested in any of this, you can use NULL instead of &xstats  
846  */  */
847          *frametype = xframe.intra;          *frametype = xframe.intra;
848          *streamlength = xframe.length;          *streamlength = xframe.length;
# Line 355  Line 850 
850          return xerr;          return xerr;
851  }  }
852    
853    /*****************************************************************************
854     * Routines for decoding: init encoder, frame step, release encoder
855     ****************************************************************************/
856    
857  /*********************************************************************/  /* init decoder before first run */
858  /* Routines for decoding: init encoder, frame step, release encoder  */  static int dec_init(int use_assembler)
 /*********************************************************************/  
   
 int dec_init(int use_assembler) /* init decoder before first run */  
859  {  {
860          int xerr;          int xerr;
861    
# Line 388  Line 883 
883          return xerr;          return xerr;
884  }  }
885    
886  int dec_main(unsigned char *m4v_buffer, unsigned char *out_buffer, int m4v_size)  /* decode one frame  */
887  {       /* decode one frame  */  static int dec_main(unsigned char *m4v_buffer, unsigned char *out_buffer,
888                                            int m4v_size)
889    {
890          int xerr;          int xerr;
891          XVID_DEC_FRAME xframe;          XVID_DEC_FRAME xframe;
892    
# Line 398  Line 894 
894          xframe.length = m4v_size;          xframe.length = m4v_size;
895          xframe.image = out_buffer;          xframe.image = out_buffer;
896          xframe.stride = XDIM;          xframe.stride = XDIM;
897          xframe.colorspace = XVID_CSP_YV12;             // XVID_CSP_USER is fastest (no memcopy involved)          xframe.colorspace = XVID_CSP_YV12;             /* XVID_CSP_USER is fastest (no memcopy involved) */
898    
899          xerr = xvid_decore(dec_handle, XVID_DEC_DECODE, &xframe, NULL);          xerr = xvid_decore(dec_handle, XVID_DEC_DECODE, &xframe, NULL);
900    
901          return xerr;          return xerr;
902  }  }
903    
904  int dec_stop()  /* close decoder to release resources */  /* close decoder to release resources */
905    static int dec_stop()
906  {  {
907          int xerr;          int xerr;
908          xerr = xvid_decore(dec_handle, XVID_DEC_DESTROY, NULL, NULL);          xerr = xvid_decore(dec_handle, XVID_DEC_DESTROY, NULL, NULL);
# Line 413  Line 910 
910          return xerr;          return xerr;
911  }  }
912    
913    /* EOF */
 /*********************************************************************/  
 /*                          Main program                             */  
 /*********************************************************************/  
   
 int main(int argc, char *argv[])  
 {  
   unsigned char *divx_buffer = NULL;  
   unsigned char *in_buffer = NULL;  
   unsigned char *out_buffer = NULL;  
   
   double enctime,dectime;  
   double totalenctime=0.;  
   double totaldectime=0.;  
   
   long totalsize=0;  
   int status;  
   
   int m4v_size;  
   int frame_type[ABS_MAXFRAMENR];  
   int Iframes=0, Pframes=0, Bframes=0, use_assembler=0;  
   double framepsnr[ABS_MAXFRAMENR];  
   
   double Ipsnr=0.,Imaxpsnr=0.,Iminpsnr=999.,Ivarpsnr=0.;  
   double Ppsnr=0.,Pmaxpsnr=0.,Pminpsnr=999.,Pvarpsnr=0.;  
   double Bpsnr=0.,Bmaxpsnr=0.,Bminpsnr=999.,Bvarpsnr=0.;  
   
   char filename[256];  
   
   FILE *filehandle;  
   
 /* read YUV in pgm format from stdin */  
   if (!pgmflag)  
   {  
         pgmflag = 1;  
   
         if (argc==2 && !strcmp(argv[1],"-asm"))  
           use_assembler = 1;  
         if (argc>=3)  
         {       XDIM = atoi(argv[1]);  
                 YDIM = atoi(argv[2]);  
                 if ( (XDIM <= 0) || (XDIM >= 2048) || (YDIM <=0) || (YDIM >= 2048) )  
                 {       fprintf(stderr,"Wrong frames size %d %d, trying PGM \n",XDIM, YDIM);  
                 }  
                 else  
                 {  
                         YDIM = YDIM*3/2; /* for YUV */  
                         pgmflag = 0;  
                 }  
         }  
   }  
   
   if (pgmflag)  
   {     if (read_pgmheader(stdin))  
             {  
               printf("Wrong input format, I want YUV encapsulated in PGM\n");  
               return 1;  
             }  
   }  
   if (argc>=4)  
   {     ARG_QUALITY = atoi(argv[3]);  
         if ( (ARG_QUALITY < 0) || (ARG_QUALITY > 6) )  
                 { fprintf(stderr,"Wrong Quality\n"); return -1; }  
         else  
                   printf("Quality %d\n",ARG_QUALITY);  
   }  
   if (argc>=5)  
   {     ARG_BITRATE = atoi(argv[4]);  
         if ( (ARG_BITRATE <= 0) )  
                 { fprintf(stderr,"Wrong Bitrate\n"); return -1; }  
         if ( (ARG_BITRATE <= 32) )  
                 { ARG_QUANTI = ARG_BITRATE;  
                   ARG_BITRATE=0;  
                   printf("Quantizer %d\n",ARG_QUANTI);  
                 }  
         else  
                   printf("Bitrate %d kbps\n",ARG_BITRATE);  
   }  
   if (argc>=6)  
   {     ARG_FRAMERATE = (float)atof(argv[5]);  
         if ( (ARG_FRAMERATE <= 0) )  
                 { fprintf(stderr,"Wrong Fraterate %s \n",argv[5]); return -1; }  
         printf("Framerate %6.3f fps\n",ARG_FRAMERATE);  
   }  
   
   if (argc>=7)  
   {     ARG_MAXFRAMENR = atoi(argv[6]);  
         if ( (ARG_MAXFRAMENR <= 0) )  
          { fprintf(stderr,"Wrong number of frames\n"); return -1; }  
         printf("max. Framenr. %d\n",ARG_MAXFRAMENR);  
   }  
   
 /* now we know the sizes, so allocate memory */  
   
   in_buffer = (unsigned char *) malloc(XDIM*YDIM);  
   if (!in_buffer)  
     goto free_all_memory;       // goto is one of the most underestimated instructions in C !!!  
   
   divx_buffer = (unsigned char *) malloc(XDIM*YDIM*2);  // this should really be enough memory!  
   if (!divx_buffer)  
     goto free_all_memory;  
   
   YDIM = YDIM*2/3; // PGM is YUV 4:2:0 format, so real image height is *2/3 of PGM picture  
   
   out_buffer = (unsigned char *) malloc(XDIM*YDIM*4);  
   if (!out_buffer)  
     goto free_all_memory;  
   
   
 /*********************************************************************/  
 /*                         XviD PART  Start                          */  
 /*********************************************************************/  
   
   
   status = enc_init(use_assembler);  
         if (status)  
         {  
                 printf("Encore INIT problem, return value %d\n", status);  
                 goto release_all;  
         }  
   
         status = dec_init(use_assembler);  
         if (status)  
         {  
                 printf("Decore INIT problem, return value %d\n", status);  
                 goto release_all;  
         }  
   
   
 /*********************************************************************/  
 /*                               Main loop                           */  
 /*********************************************************************/  
   
   do  
     {  
         if (pgmflag)  
               status = read_pgmdata(stdin, in_buffer);  // read PGM data (YUV-format)  
         else  
               status = read_yuvdata(stdin, in_buffer);  // read raw data (YUV-format)  
   
       if (status)  
         {  
           // Couldn't read image, most likely end-of-file  
           continue;  
         }  
   
   
       if (save_ref_flag)  
         {  
                 sprintf(filename, "%s%05d.pgm", filepath, filenr);  
                 write_pgm(filename,in_buffer);  
         }  
   
   
 /*********************************************************************/  
 /*               analyse this frame before encoding                  */  
 /*********************************************************************/  
   
 //      nothing is done here at the moment, but you could e.g. create  
 //      histograms or measure entropy or apply preprocessing filters...  
   
 /*********************************************************************/  
 /*               encode and decode this frame                        */  
 /*********************************************************************/  
   
         enctime = -msecond();  
         status = enc_main(in_buffer, divx_buffer, &m4v_size, &frame_type[filenr]);  
         enctime += msecond();  
   
         totalenctime += enctime;  
         totalsize += m4v_size;  
   
         printf("Frame %5d: intra %d, enctime =%6.1f ms length=%7d bytes ",  
                  filenr, frame_type[filenr], enctime*1000, m4v_size);  
   
         if (save_m4v_flag)  
         {  
                 sprintf(filename, "%sframe%05d.m4v", filepath, filenr);  
                 filehandle = fopen(filename, "wb");  
                 fwrite(divx_buffer, m4v_size, 1, filehandle);  
                 fclose(filehandle);  
         }  
   
         dectime = -msecond();  
         status = dec_main(divx_buffer, out_buffer, m4v_size);  
         dectime += msecond();  
   
         totaldectime += dectime;  
   
   
 /*********************************************************************/  
 /*        analyse the decoded frame and compare to original          */  
 /*********************************************************************/  
   
         framepsnr[filenr] = PSNR(XDIM,YDIM, in_buffer, XDIM, out_buffer, XDIM );  
   
         printf("dectime =%6.1f ms PSNR %5.2f\n",dectime*1000, framepsnr[filenr]);  
   
         if (save_dec_flag)  
         {  
                 sprintf(filename, "%sdec%05d.pgm", filepath, filenr);  
                 write_pgm(filename,out_buffer);  
         }  
   
         if (pgmflag)  
                 status = read_pgmheader(stdin);         // because if this was the last PGM, stop now  
   
         filenr++;  
   
    } while ( (!status) && (filenr<ARG_MAXFRAMENR) );  
   
   
   
 /*********************************************************************/  
 /*     calculate totals and averages for output, print results       */  
 /*********************************************************************/  
   
         totalsize    /= filenr;  
         totalenctime /= filenr;  
         totaldectime /= filenr;  
   
         for (i=0;i<filenr;i++)  
         {  
                 switch (frame_type[i])  
                 {  
                 case 0:  
                         Pframes++;  
                         Ppsnr += framepsnr[i];  
                         break;  
                 case 1:  
                         Iframes++;  
                         Ipsnr += framepsnr[i];  
                         break;  
                 case 2:  
                 default:  
                         Bframes++;  
                         Bpsnr += framepsnr[i];  
                         break;  
                 }  
         }  
   
         if (Pframes)  
                 Ppsnr /= Pframes;  
         if (Iframes)  
                 Ipsnr /= Iframes;  
         if (Bframes)  
                 Bpsnr /= Bframes;  
   
   
         for (i=0;i<filenr;i++)  // calculate statistics for every frametype: P,I (and B)  
         {  
                 switch (frame_type[i])  
                 {  
                 case 0:  
                         if (framepsnr[i] > Pmaxpsnr)  
                                 Pmaxpsnr = framepsnr[i];  
                         if (framepsnr[i] < Pminpsnr)  
                                 Pminpsnr = framepsnr[i];  
                         Pvarpsnr += (framepsnr[i] - Ppsnr)*(framepsnr[i] - Ppsnr) /Pframes;  
                         break;  
                 case 1:  
                         if (framepsnr[i] > Imaxpsnr)  
                                 Imaxpsnr = framepsnr[i];  
                         if (framepsnr[i] < Pminpsnr)  
                                 Iminpsnr = framepsnr[i];  
                         Ivarpsnr += (framepsnr[i] - Ipsnr)*(framepsnr[i] - Ipsnr) /Iframes;  
                         break;  
                 case 2:  
                         if (framepsnr[i] > Bmaxpsnr)  
                                 Bmaxpsnr = framepsnr[i];  
                         if (framepsnr[i] < Pminpsnr)  
                                 Bminpsnr = framepsnr[i];  
                         Bvarpsnr += (framepsnr[i] - Bpsnr)*(framepsnr[i] - Bpsnr) /Bframes;  
                         break;  
                 }  
         }  
   
         printf("Avg. Q%1d %2s ",ARG_QUALITY, (ARG_QUANTI ? " q" : "br"));  
         printf("%04d ",MAX(ARG_QUANTI,ARG_BITRATE));  
         printf("( %.2f bpp) ", (double)ARG_BITRATE*1000/XDIM/YDIM/ARG_FRAMERATE);  
         printf("size %6d ",totalsize);  
         printf("( %4d kbps ",(int)(totalsize*8*ARG_FRAMERATE/1000));  
         printf("/ %.2f bpp) ",(double)totalsize*8/XDIM/YDIM);  
         printf("enc: %6.1f fps, dec: %6.1f fps \n",1/totalenctime, 1/totaldectime);  
         printf("PSNR P(%d): %5.2f ( %5.2f , %5.2f ; %5.4f ) ",Pframes,Ppsnr,Pminpsnr,Pmaxpsnr,sqrt(Pvarpsnr/filenr));  
         printf("I(%d): %5.2f ( %5.2f , %5.2f ; %5.4f ) ",Iframes,Ipsnr,Iminpsnr,Imaxpsnr,sqrt(Ivarpsnr/filenr));  
         if (Bframes)  
                 printf("B(%d): %5.2f ( %5.2f , %5.2f ; %5.4f ) ",Bframes,Bpsnr,Bminpsnr,Bmaxpsnr,sqrt(Bvarpsnr/filenr));  
         printf("\n");  
   
 /*********************************************************************/  
 /*                         XviD PART  Stop                           */  
 /*********************************************************************/  
   
 release_all:  
   
         if (enc_handle)  
         {  
                 status = enc_stop();  
                 if (status)  
                         printf("Encore RELEASE problem return value %d\n", status);  
         }  
   
         if (dec_handle)  
         {  
                 status = dec_stop();  
                 if (status)  
                         printf("Decore RELEASE problem return value %d\n", status);  
         }  
   
   
 free_all_memory:  
         free(out_buffer);  
         free(divx_buffer);  
         free(in_buffer);  
   
   return 0;  
 }  
Legend:
Removed from v.211  
changed lines
  Added in v.684
No admin address has been configured
 ViewVC Help
Powered by ViewVC 1.0.4