H263编解码器，使用纯C语言实现

/************************************************************************ * * getpic.c, picture decoding for tmndecode (H.263 decoder) * Copyright (C) 1995, 1996 Telenor R&D, Norway * * Contacts: * Robert Danielsen <Robert.Danielsen@nta.no> * * Telenor Research and Development http://www.nta.no/brukere/DVC/ * P.O.Box 83 tel.: +47 63 84 84 00 * N-2007 Kjeller, Norway fax.: +47 63 81 00 76 * * Copyright (C) 1997 University of BC, Canada * Modified by: Mike Gallant <mikeg@ee.ubc.ca> * Guy Cote <guyc@ee.ubc.ca> * Berna Erol <bernae@ee.ubc.ca> * * Contacts: * Micheal Gallant <mikeg@ee.ubc.ca> * * UBC Image Processing Laboratory http://www.ee.ubc.ca/image * 2356 Main Mall tel.: +1 604 822 4051 * Vancouver BC Canada V6T1Z4 fax.: +1 604 822 5949 * ************************************************************************/ /* Disclaimer of Warranty * * These software programs are available to the user without any license fee * or royalty on an "as is" basis. The University of British Columbia * disclaims any and all warranties, whether express, implied, or * statuary, including any implied warranties or merchantability or of * fitness for a particular purpose. In no event shall the * copyright-holder be liable for any incidental, punitive, or * consequential damages of any kind whatsoever arising from the use of * these programs. * * This disclaimer of warranty extends to the user of these programs and * user's customers, employees, agents, transferees, successors, and * assigns. * * The University of British Columbia does not represent or warrant that the * programs furnished hereunder are free of infringement of any * third-party patents. * * Commercial implementations of H.263, including shareware, are subject to * royalty fees to patent holders. Many of these patents are general * enough such that they are unavoidable regardless of implementation * design. * */ /* modified to support annex O true B frames, mikeg@ee.ubc.ca * * modified by Wayne Ellis BT Labs to run Annex E Arithmetic Decoding * <ellis_w_wayne@bt-web.bt.co.uk> * * based on mpeg2decode, (C) 1994, MPEG Software Simulation Group and * mpeg2play, (C) 1994 Stefan Eckart <stefan@lis.e-technik.tu-muenchen.de> */ #include <stdio.h> #include <stdlib.h> #include <string.h> #include "config.h" #include "tmndec.h" #include "global.h" #include "indices.h" #include "sactbls.h" static int coded_map[MBR + 1][MBC + 1]; static int quant_map[MBR + 1][MBC + 1]; static int STRENGTH[] = {1, 1, 2, 2, 3, 3, 4, 4, 4, 5, 5, 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 11, 11, 11, 12, 12, 12}; static int STRENGTH1[] = {1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4}; static int STRENGTH2[] = {1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3}; /* private prototypes */ static void get_I_P_MBs _ANSI_ARGS_ ((int framenum, int gob)); static void get_B_MBs _ANSI_ARGS_ ((int framenum)); static void get_EI_EP_MBs _ANSI_ARGS_ ((int framenum)); static void clearblock _ANSI_ARGS_ ((int comp)); static int motion_decode _ANSI_ARGS_ ((int vec, int pmv)); static int find_pmv _ANSI_ARGS_ ((int x, int y, int block, int comp)); static void addblock _ANSI_ARGS_ ((int comp, int bx, int by, int addflag)); static void reconblock_b _ANSI_ARGS_ ((int comp, int bx, int by, int mode, int bdx, int bdy)); static void find_bidir_limits _ANSI_ARGS_ ((int vec, int *start, int *stop, int nhv)); static void find_bidir_chroma_limits _ANSI_ARGS_ ((int vec, int *start, int *stop)); static void make_edge_image _ANSI_ARGS_ ((unsigned char *src, unsigned char *dst, int width, int height, int edge)); static void init_enhancement_layer _ANSI_ARGS_((int layer)); void edge_filter _ANSI_ARGS_ ((unsigned char *lum, unsigned char *Cb, unsigned char *Cr, int width, int height)); void horiz_edge_filter _ANSI_ARGS_ ((unsigned char *rec, int width, int height, int chr)); void vert_edge_filter _ANSI_ARGS_ ((unsigned char *rec, int width, int height, int chr)); void vert_post_filter (unsigned char *rec, int width, int height, int chr); void horiz_post_filter (unsigned char *rec, int width, int height, int chr); void PostFilter (unsigned char *lum, unsigned char *Cb, unsigned char *Cr, int width, int height); /* reference picture selection */ int get_reference_picture(void); void store_picture(int); /* error concealment */ void conceal_missing_gobs(int start_mb_row_missing, int number_of_mb_rows_missing); /********************************************************************** * * Name: getpicture * Description: decode and display one picture * Input: frame number * Returns: * Side effects: * * Date: 971102 Author: mikeg@ee.ubc.ca * * ***********************************************************************/ void getpicture (int *framenum, int gob) { unsigned char *tmp; int i, store_pb, quality=0; static int absolute_temp_ref = 0; static int prev_temp_ref = 0; FILE *fout; switch (pict_type) { case PCT_INTRA: case PCT_INTER: case PCT_PB: case PCT_IPB: enhance_pict = NO; if (0 == UFEP && plus_type) { mv_outside_frame = prev_mv_outside_frame; deblocking_filter_mode = prev_df; adv_pred_mode = prev_adv_pred; overlapping_MC = prev_obmc; use_4mv = prev_4mv; long_vectors = prev_long_vectors; syntax_arith_coding = prev_sac; advanced_intra_coding = prev_aic; slice_structured_mode = prev_slice_struct; reference_picture_selection_mode = prev_rps; independently_segmented_decoding_mode = prev_isd; alternative_inter_VLC_mode = prev_aivlc; modified_quantization_mode = prev_mq; } for (i = 0; i < 3; i++) { tmp = prev_frame[i]; prev_I_P_frame[i] = prev_frame[i] = next_I_P_frame[i]; next_I_P_frame[i] = current_frame[i] = tmp; } /* if the picture header is missing, assume temp_ref increments */ /* not done yet*/ /* Annex N: get the correct reference picture for prediction */ if (reference_picture_selection_mode) { if (-1 == (quality = get_reference_picture())) break; if (quality > 256) { fprintf(stderr,"completely out of sync -- waiting for I-frame\n"); stop_decoder = 1; break; } } #ifdef DEBUG if (temp_ref < prev_temp_ref) prev_temp_ref -= 256; absolute_temp_ref += temp_ref - prev_temp_ref; prev_temp_ref = temp_ref; fprintf(stdout, "%d\n", absolute_temp_ref); #endif if ((mv_outside_frame) && (*framenum > 0)) { make_edge_image (prev_I_P_frame[0], edgeframe[0], coded_picture_width, coded_picture_height, 32); make_edge_image (prev_I_P_frame[1], edgeframe[1], chrom_width, chrom_height, 16); make_edge_image (prev_I_P_frame[2], edgeframe[2], chrom_width, chrom_height, 16); } get_I_P_MBs (*framenum, gob); store_pb = pb_frame; pb_frame = 0; if (deblocking_filter_mode) { edge_filter (current_frame[0], current_frame[1], current_frame[2], coded_picture_width, coded_picture_height); printf("do deblocking\n"); } #if 1 fout = fopen("origen2.yuv", "wb"); if (fout == NULL) { printf("open output file fail\n"); exit(-1); } for(i=0;i<288;i++) { if (fwrite((current_frame[0]+coded_picture_width*i), 352, 1, fout) < 0) { printf("error/n"); break; } } for(i=0;i<144;i++) { if(fwrite((current_frame[1]+(coded_picture_width>>1)*i),