www.pudn.com > JM.rar > ldecod.c
/*!
***********************************************************************
* \mainpage
* This is the H.264/AVC decoder reference software. For detailed documentation
* see the comments in each file.
*
* \author
* The main contributors are listed in contributors.h
*
* \version
* JM 9.0 (FRExt)
*
* \note
* tags are used for document system "doxygen"
* available at http://www.doxygen.org
*/
/*!
* \file
* ldecod.c
* \brief
* H.264/AVC reference decoder project main()
* \author
* Main contributors (see contributors.h for copyright, address and affiliation details)
* - Inge Lille-Langøy
* - Rickard Sjoberg
* - Stephan Wenger
* - Jani Lainema
* - Sebastian Purreiter
* - Byeong-Moon Jeon
* - Gabi Blaettermann
* - Ye-Kui Wang
* - Valeri George
* - Karsten Suehring
*
***********************************************************************
*/
#include "contributors.h"
#include
#include
#include
#include
#include
#if defined WIN32
#include
#else
#include
#endif
#include
#include
#include
#include "global.h"
#include "rtp.h"
#include "memalloc.h"
#include "mbuffer.h"
#include "leaky_bucket.h"
#include "fmo.h"
#include "annexb.h"
#include "output.h"
#include "cabac.h"
#include "erc_api.h"
#define JM "9 (FRExt)"
#define VERSION "9.0"
#define EXT_VERSION "(FRExt)"
#define LOGFILE "log.dec"
#define DATADECFILE "dataDec.txt"
#define TRACEFILE "trace_dec.txt"
extern objectBuffer_t *erc_object_list;
extern ercVariables_t *erc_errorVar;
extern ColocatedParams *Co_located;
// I have started to move the inp and img structures into global variables.
// They are declared in the following lines. Since inp is defined in conio.h
// and cannot be overridden globally, it is defined here as input
//
// Everywhere, input-> and img-> can now be used either globally or with
// the local override through the formal parameter mechanism
extern FILE* bits;
extern StorablePicture* dec_picture;
struct inp_par *input; //!< input parameters from input configuration file
struct snr_par *snr; //!< statistics
struct img_par *img; //!< image parameters
int global_init_done = 0;
/*!
***********************************************************************
* \brief
* main function for TML decoder
***********************************************************************
*/
int main(int argc, char **argv)
{
// allocate memory for the structures
if ((input = (struct inp_par *)calloc(1, sizeof(struct inp_par)))==NULL) no_mem_exit("main: input");
if ((snr = (struct snr_par *)calloc(1, sizeof(struct snr_par)))==NULL) no_mem_exit("main: snr");
if ((img = (struct img_par *)calloc(1, sizeof(struct img_par)))==NULL) no_mem_exit("main: img");
// Read Configuration File
if (argc != 2)
{
snprintf(errortext, ET_SIZE, "Usage: %s \n\t defines decoder parameters",argv[0]);
error(errortext, 300);
}
init_conf(input, argv[1]);
init_old_slice();
switch (input->FileFormat)
{
case 0:
OpenBitstreamFile (input->infile);
break;
case 1:
OpenRTPFile (input->infile);
break;
default:
printf ("Unsupported file format %d, exit\n", input->FileFormat);
}
// Allocate Slice data struct
malloc_slice(input,img);
init(img);
dec_picture = NULL;
dpb.init_done = 0;
g_nFrame = 0;
// init_dpb(input);
init_out_buffer();
img->idr_psnr_number=input->ref_offset;
img->psnr_number=0;
img->number=0;
img->type = I_SLICE;
img->dec_ref_pic_marking_buffer = NULL;
// B pictures
Bframe_ctr=0;
// time for total decoding session
tot_time = 0;
while (decode_one_frame(img, input, snr) != EOS)
;
report(input, img, snr);
free_slice(input,img);
FmoFinit();
free_global_buffers();
flush_dpb();
#ifdef PAIR_FIELDS_IN_OUTPUT
flush_pending_output(p_out);
#endif
CloseBitstreamFile();
close(p_out);
// fclose(p_out2);
if (p_ref!=-1)
close(p_ref);
#if TRACE
fclose(p_trace);
#endif
ercClose(erc_errorVar);
free_dpb();
uninit_out_buffer();
free_collocated(Co_located);
free (input);
free (snr);
free (img);
//while( !kbhit() );
return 0;
}
/*!
***********************************************************************
* \brief
* Initilize some arrays
***********************************************************************
*/
void init(struct img_par *img) //!< image parameters
{
img->oldFrameSizeInMbs = -1;
imgY_ref = NULL;
imgUV_ref = NULL;
}
/*!
***********************************************************************
* \brief
* Initilize FREXT variables
***********************************************************************
*/
void init_frext(struct img_par *img) //!< image parameters
{
//pel bitdepth init
img->bitdepth_luma_qp_scale = 6*(img->bitdepth_luma - 8);
if(img->bitdepth_luma > img->bitdepth_chroma || img->yuv_format == YUV400)
img->pic_unit_bitsize_on_disk = (img->bitdepth_luma > 8)? 16:8;
else
img->pic_unit_bitsize_on_disk = (img->bitdepth_chroma > 8)? 16:8;
img->dc_pred_value = 1<<(img->bitdepth_luma - 1);
img->max_imgpel_value = (1<bitdepth_luma) - 1;
if (img->yuv_format != YUV400)
{
//for chrominance part
img->bitdepth_chroma_qp_scale = 6*(img->bitdepth_chroma - 8);
img->max_imgpel_value_uv = (1<bitdepth_chroma) - 1;
img->num_blk8x8_uv = (1<yuv_format)&(~(0x1));
img->num_cdc_coeff = img->num_blk8x8_uv<<1;
img->mb_cr_size_x = (img->yuv_format==YUV420 || img->yuv_format==YUV422)? 8:16;
img->mb_cr_size_y = (img->yuv_format==YUV444 || img->yuv_format==YUV422)? 16:8;
// Residue Color Transform
if(img->residue_transform_flag)
img->bitdepth_chroma_qp_scale += 6;
}
else
{
img->bitdepth_chroma_qp_scale = 0;
img->max_imgpel_value_uv = 0;
img->num_blk8x8_uv = 0;
img->num_cdc_coeff = 0;
img->mb_cr_size_x = 0;
img->mb_cr_size_y = 0;
}
}
/*!
************************************************************************
* \brief
* Read input from configuration file
*
* \par Input:
* Name of configuration filename
*
* \par Output
* none
************************************************************************
*/
void init_conf(struct inp_par *inp,
char *config_filename)
{
FILE *fd;
int NAL_mode;
// read the decoder configuration file
if((fd=fopen(config_filename,"r")) == NULL)
{
snprintf(errortext, ET_SIZE, "Error: Control file %s not found\n",config_filename);
error(errortext, 300);
}
fscanf(fd,"%s",inp->infile); // H.264 compressed input bitsream
fscanf(fd,"%*[^\n]");
fscanf(fd,"%s",inp->outfile); // YUV 4:2:2 input format
fscanf(fd,"%*[^\n]");
fscanf(fd,"%s",inp->reffile); // reference file
fscanf(fd,"%*[^\n]");
// Frame buffer size
fscanf(fd,"%d,",&inp->dpb_size); // may be overwritten in case of RTP NAL
fscanf(fd,"%*[^\n]");
if (inp->dpb_size < 1)
{
snprintf(errortext, ET_SIZE, "Decoded Picture Buffer Size is %d. It has to be at least 1",inp->dpb_size);
error(errortext,1);
}
fscanf(fd,"%d",&(NAL_mode)); // NAL mode
fscanf(fd,"%*[^\n]");
switch(NAL_mode)
{
case 0:
inp->FileFormat = PAR_OF_ANNEXB;
break;
case 1:
inp->FileFormat = PAR_OF_RTP;
break;
default:
snprintf(errortext, ET_SIZE, "NAL mode %i is not supported", NAL_mode);
error(errortext,400);
}
fscanf(fd,"%d,",&inp->ref_offset); // offset used for SNR computation
fscanf(fd,"%*[^\n]");
fscanf(fd,"%d,",&inp->poc_scale); // offset used for SNR computation
fscanf(fd,"%*[^\n]");
if (inp->poc_scale < 1 || inp->poc_scale > 10)
{
snprintf(errortext, ET_SIZE, "Poc Scale is %d. It has to be within range 1 to 10",inp->poc_scale);
error(errortext,1);
}
#ifdef _LEAKYBUCKET_
fscanf(fd,"%ld,",&inp->R_decoder); // Decoder rate
fscanf(fd, "%*[^\n]");
fscanf(fd,"%ld,",&inp->B_decoder); // Decoder buffer size
fscanf(fd, "%*[^\n]");
fscanf(fd,"%ld,",&inp->F_decoder); // Decoder initial delay
fscanf(fd, "%*[^\n]");
fscanf(fd,"%s",inp->LeakyBucketParamFile); // file where Leaky Bucket params (computed by encoder) are stored
fscanf(fd,"%*[^\n]");
#endif
fclose (fd);
#if TRACE
if ((p_trace=fopen(TRACEFILE,"w"))==0) // append new statistic at the end
{
snprintf(errortext, ET_SIZE, "Error open file %s!",TRACEFILE);
error(errortext,500);
}
#endif
if ((p_out=open(inp->outfile, OPENFLAGS_WRITE, OPEN_PERMISSIONS))==-1)
{
snprintf(errortext, ET_SIZE, "Error open file %s ",inp->outfile);
error(errortext,500);
}
/* if ((p_out2=fopen("out.yuv","wb"))==0)
{
snprintf(errortext, ET_SIZE, "Error open file %s ",inp->outfile);
error(errortext,500);
}*/
fprintf(stdout,"----------------------------- JM %s %s -----------------------------\n", VERSION, EXT_VERSION);
fprintf(stdout," Decoder config file : %s \n",config_filename);
fprintf(stdout,"--------------------------------------------------------------------------\n");
fprintf(stdout," Input H.264 bitstream : %s \n",inp->infile);
fprintf(stdout," Output decoded YUV : %s \n",inp->outfile);
fprintf(stdout," Output status file : %s \n",LOGFILE);
if ((p_ref=open(inp->reffile,OPENFLAGS_READ))==-1)
{
fprintf(stdout," Input reference file : %s does not exist \n",inp->reffile);
fprintf(stdout," SNR values are not available\n");
}
else
fprintf(stdout," Input reference file : %s \n",inp->reffile);
fprintf(stdout,"--------------------------------------------------------------------------\n");
#ifdef _LEAKYBUCKET_
fprintf(stdout," Rate_decoder : %8ld \n",inp->R_decoder);
fprintf(stdout," B_decoder : %8ld \n",inp->B_decoder);
fprintf(stdout," F_decoder : %8ld \n",inp->F_decoder);
fprintf(stdout," LeakyBucketParamFile: %s \n",inp->LeakyBucketParamFile); // Leaky Bucket Param file
calc_buffer(inp);
fprintf(stdout,"--------------------------------------------------------------------------\n");
#endif
fprintf(stdout,"POC must = frame# or field# for SNRs to be correct\n");
fprintf(stdout,"Frame POC Pic# QP SnrY SnrU SnrV Y:U:V Time(ms)\n");
}
/*!
************************************************************************
* \brief
* Reports the gathered information to appropriate outputs
*
* \par Input:
* struct inp_par *inp,
* struct img_par *img,
* struct snr_par *stat
*
* \par Output:
* None
************************************************************************
*/
void report(struct inp_par *inp, struct img_par *img, struct snr_par *snr)
{
#define OUTSTRING_SIZE 255
char string[OUTSTRING_SIZE];
FILE *p_log;
char yuv_formats[4][4]= { {"400"}, {"420"}, {"422"}, {"444"} };
#ifndef WIN32
time_t now;
struct tm *l_time;
#else
char timebuf[128];
#endif
fprintf(stdout,"-------------------- Average SNR all frames ------------------------------\n");
fprintf(stdout," SNR Y(dB) : %5.2f\n",snr->snr_ya);
fprintf(stdout," SNR U(dB) : %5.2f\n",snr->snr_ua);
fprintf(stdout," SNR V(dB) : %5.2f\n",snr->snr_va);
fprintf(stdout," Total decoding time : %.3f sec \n",tot_time*0.001);
fprintf(stdout,"--------------------------------------------------------------------------\n");
fprintf(stdout," Exit JM %s decoder, ver %s ",JM, VERSION);
fprintf(stdout,"\n");
// write to log file
snprintf(string, OUTSTRING_SIZE, "%s", LOGFILE);
if ((p_log=fopen(string,"r"))==0) // check if file exist
{
if ((p_log=fopen(string,"a"))==0)
{
snprintf(errortext, ET_SIZE, "Error open file %s for appending",string);
error(errortext, 500);
}
else // Create header to new file
{
fprintf(p_log," -------------------------------------------------------------------------------------------------------------------\n");
fprintf(p_log,"| Decoder statistics. This file is made first time, later runs are appended |\n");
fprintf(p_log," ------------------------------------------------------------------------------------------------------------------- \n");
fprintf(p_log,"| ver | Date | Time | Sequence |#Img| Format | YUV |Coding|SNRY 1|SNRU 1|SNRV 1|SNRY N|SNRU N|SNRV N|\n");
fprintf(p_log," -------------------------------------------------------------------------------------------------------------------\n");
}
}
else
{
fclose(p_log);
p_log=fopen(string,"a"); // File exist,just open for appending
}
fprintf(p_log,"|%s/%-4s", VERSION, EXT_VERSION);
#ifdef WIN32
_strdate( timebuf );
fprintf(p_log,"| %1.5s |",timebuf );
_strtime( timebuf);
fprintf(p_log," % 1.5s |",timebuf);
#else
now = time ((time_t *) NULL); // Get the system time and put it into 'now' as 'calender time'
time (&now);
l_time = localtime (&now);
strftime (string, sizeof string, "%d-%b-%Y", l_time);
fprintf(p_log,"| %1.5s |",string );
strftime (string, sizeof string, "%H:%M:%S", l_time);
fprintf(p_log,"| %1.5s |",string );
#endif
fprintf(p_log,"%20.20s|",inp->infile);
fprintf(p_log,"%3d |",img->number);
fprintf(p_log,"%4dx%-4d|", img->width, img->height);
fprintf(p_log," %s |", &(yuv_formats[img->yuv_format][0]));
if (active_pps->entropy_coding_mode_flag == UVLC)
fprintf(p_log," CAVLC|");
else
fprintf(p_log," CABAC|");
fprintf(p_log,"%6.3f|",snr->snr_y1);
fprintf(p_log,"%6.3f|",snr->snr_u1);
fprintf(p_log,"%6.3f|",snr->snr_v1);
fprintf(p_log,"%6.3f|",snr->snr_ya);
fprintf(p_log,"%6.3f|",snr->snr_ua);
fprintf(p_log,"%6.3f|\n",snr->snr_va);
fclose(p_log);
snprintf(string, OUTSTRING_SIZE,"%s", DATADECFILE);
p_log=fopen(string,"a");
if(Bframe_ctr != 0) // B picture used
{
fprintf(p_log, "%3d %2d %2d %2.2f %2.2f %2.2f %5d "
"%2.2f %2.2f %2.2f %5d "
"%2.2f %2.2f %2.2f %5d %.3f\n",
img->number, 0, img->qp,
snr->snr_y1,
snr->snr_u1,
snr->snr_v1,
0,
0.0,
0.0,
0.0,
0,
snr->snr_ya,
snr->snr_ua,
snr->snr_va,
0,
(double)0.001*tot_time/(img->number+Bframe_ctr-1));
}
else
{
fprintf(p_log, "%3d %2d %2d %2.2f %2.2f %2.2f %5d "
"%2.2f %2.2f %2.2f %5d "
"%2.2f %2.2f %2.2f %5d %.3f\n",
img->number, 0, img->qp,
snr->snr_y1,
snr->snr_u1,
snr->snr_v1,
0,
0.0,
0.0,
0.0,
0,
snr->snr_ya,
snr->snr_ua,
snr->snr_va,
0,
(double)0.001*tot_time/img->number);
}
fclose(p_log);
}
/*!
************************************************************************
* \brief
* Allocates a stand-alone partition structure. Structure should
* be freed by FreePartition();
* data structures
*
* \par Input:
* n: number of partitions in the array
* \par return
* pointer to DataPartition Structure, zero-initialized
************************************************************************
*/
DataPartition *AllocPartition(int n)
{
DataPartition *partArr, *dataPart;
int i;
partArr = (DataPartition *) calloc(n, sizeof(DataPartition));
if (partArr == NULL)
{
snprintf(errortext, ET_SIZE, "AllocPartition: Memory allocation for Data Partition failed");
error(errortext, 100);
}
for (i=0; ibitstream = (Bitstream *) calloc(1, sizeof(Bitstream));
if (dataPart->bitstream == NULL)
{
snprintf(errortext, ET_SIZE, "AllocPartition: Memory allocation for Bitstream failed");
error(errortext, 100);
}
dataPart->bitstream->streamBuffer = (byte *) calloc(MAX_CODED_FRAME_SIZE, sizeof(byte));
if (dataPart->bitstream->streamBuffer == NULL)
{
snprintf(errortext, ET_SIZE, "AllocPartition: Memory allocation for streamBuffer failed");
error(errortext, 100);
}
}
return partArr;
}
/*!
************************************************************************
* \brief
* Frees a partition structure (array).
*
* \par Input:
* Partition to be freed, size of partition Array (Number of Partitions)
*
* \par return
* None
*
* \note
* n must be the same as for the corresponding call of AllocPartition
************************************************************************
*/
void FreePartition (DataPartition *dp, int n)
{
int i;
assert (dp != NULL);
assert (dp->bitstream != NULL);
assert (dp->bitstream->streamBuffer != NULL);
for (i=0; istreamBuffer);
free (dp[i].bitstream);
}
free (dp);
}
/*!
************************************************************************
* \brief
* Allocates the slice structure along with its dependent
* data structures
*
* \par Input:
* Input Parameters struct inp_par *inp, struct img_par *img
************************************************************************
*/
void malloc_slice(struct inp_par *inp, struct img_par *img)
{
Slice *currSlice;
img->currentSlice = (Slice *) calloc(1, sizeof(Slice));
if ( (currSlice = img->currentSlice) == NULL)
{
snprintf(errortext, ET_SIZE, "Memory allocation for Slice datastruct in NAL-mode %d failed", inp->FileFormat);
error(errortext,100);
}
// img->currentSlice->rmpni_buffer=NULL;
//! you don't know whether we do CABAC hre, hence initialize CABAC anyway
// if (inp->symbol_mode == CABAC)
if (1)
{
// create all context models
currSlice->mot_ctx = create_contexts_MotionInfo();
currSlice->tex_ctx = create_contexts_TextureInfo();
}
currSlice->max_part_nr = 3; //! assume data partitioning (worst case) for the following mallocs()
currSlice->partArr = AllocPartition(currSlice->max_part_nr);
}
/*!
************************************************************************
* \brief
* Memory frees of the Slice structure and of its dependent
* data structures
*
* \par Input:
* Input Parameters struct inp_par *inp, struct img_par *img
************************************************************************
*/
void free_slice(struct inp_par *inp, struct img_par *img)
{
Slice *currSlice = img->currentSlice;
FreePartition (currSlice->partArr, 3);
// if (inp->symbol_mode == CABAC)
if (1)
{
// delete all context models
delete_contexts_MotionInfo(currSlice->mot_ctx);
delete_contexts_TextureInfo(currSlice->tex_ctx);
}
free(img->currentSlice);
currSlice = NULL;
}
/*!
************************************************************************
* \brief
* Dynamic memory allocation of frame size related global buffers
* buffers are defined in global.h, allocated memory must be freed in
* void free_global_buffers()
*
* \par Input:
* Input Parameters struct inp_par *inp, Image Parameters struct img_par *img
*
* \par Output:
* Number of allocated bytes
***********************************************************************
*/
int init_global_buffers()
{
int memory_size=0;
int quad_range, i;
if (global_init_done)
{
free_global_buffers();
}
// allocate memory for reference frame in find_snr
memory_size += get_mem2Dpel(&imgY_ref, img->height, img->width);
if (img->yuv_format != YUV400)
memory_size += get_mem3Dpel(&imgUV_ref, 2, img->height_cr, img->width_cr);
// allocate memory in structure img
if(((img->mb_data) = (Macroblock *) calloc(img->FrameSizeInMbs, sizeof(Macroblock))) == NULL)
no_mem_exit("init_global_buffers: img->mb_data");
if(((img->intra_block) = (int*)calloc(img->FrameSizeInMbs, sizeof(int))) == NULL)
no_mem_exit("init_global_buffers: img->intra_block");
memory_size += get_mem2Dint(&(img->ipredmode), 4*img->PicWidthInMbs , 4*img->FrameHeightInMbs);
memory_size += get_mem2Dint(&(img->field_anchor),4*img->FrameHeightInMbs, 4*img->PicWidthInMbs);
memory_size += get_mem3Dint(&(img->wp_weight), 2, MAX_REFERENCE_PICTURES, 3);
memory_size += get_mem3Dint(&(img->wp_offset), 6, MAX_REFERENCE_PICTURES, 3);
memory_size += get_mem4Dint(&(img->wbp_weight), 6, MAX_REFERENCE_PICTURES, MAX_REFERENCE_PICTURES, 3);
// CAVLC mem
memory_size += get_mem3Dint(&(img->nz_coeff), img->FrameSizeInMbs, 4, 4 + img->num_blk8x8_uv);
memory_size += get_mem2Dint(&(img->siblock),img->PicWidthInMbs , img->FrameHeightInMbs);
if(img->max_imgpel_value > img->max_imgpel_value_uv || img->yuv_format == YUV400)
quad_range = (img->max_imgpel_value + 1) * 2;
else
quad_range = (img->max_imgpel_value_uv + 1) * 2;
if ((img->quad = (int*)calloc (quad_range, sizeof(int))) == NULL)
no_mem_exit ("init_img: img->quad");
for (i=0; i < quad_range/2; ++i)
{
img->quad[i]=i*i;
}
global_init_done = 1;
img->oldFrameSizeInMbs = img->FrameSizeInMbs;
return (memory_size);
}
/*!
************************************************************************
* \brief
* Free allocated memory of frame size related global buffers
* buffers are defined in global.h, allocated memory is allocated in
* int init_global_buffers()
*
* \par Input:
* Input Parameters struct inp_par *inp, Image Parameters struct img_par *img
*
* \par Output:
* none
*
************************************************************************
*/
void free_global_buffers()
{
free_mem2Dpel (imgY_ref);
if (img->yuv_format != YUV400)
free_mem3Dpel (imgUV_ref,2);
// CAVLC free mem
free_mem3Dint(img->nz_coeff, img->oldFrameSizeInMbs);
free_mem2Dint(img->siblock);
// free mem, allocated for structure img
if (img->mb_data != NULL) free(img->mb_data);
free (img->intra_block);
free_mem2Dint (img->ipredmode);
free_mem2Dint(img->field_anchor);
free_mem3Dint(img->wp_weight, 2);
free_mem3Dint(img->wp_offset, 6);
free_mem4Dint(img->wbp_weight, 6, MAX_REFERENCE_PICTURES);
free (img->quad);
global_init_done = 0;
}