home *** CD-ROM | disk | FTP | other *** search
/ PC Welt 2004 March / PCWELT_3_2004.ISO / pcwsoft / flaskmpeg_078_39_src.z.exe / flaskmpeg / Video / motion.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-10-28  |  7.3 KB  |  245 lines
  1. /* motion.c, motion vector decoding                                         */
  2.  
  3. /* Copyright (C) 1996, MPEG Software Simulation Group. All Rights Reserved. */
  4.  
  5. /*
  6.  * Disclaimer of Warranty
  7.  *
  8.  * These software programs are available to the user without any license fee or
  9.  * royalty on an "as is" basis.  The MPEG Software Simulation Group disclaims
  10.  * any and all warranties, whether express, implied, or statuary, including any
  11.  * implied warranties or merchantability or of fitness for a particular
  12.  * purpose.  In no event shall the copyright-holder be liable for any
  13.  * incidental, punitive, or consequential damages of any kind whatsoever
  14.  * arising from the use of these programs.
  15.  *
  16.  * This disclaimer of warranty extends to the user of these programs and user's
  17.  * customers, employees, agents, transferees, successors, and assigns.
  18.  *
  19.  * The MPEG Software Simulation Group does not represent or warrant that the
  20.  * programs furnished hereunder are free of infringement of any third-party
  21.  * patents.
  22.  *
  23.  * Commercial implementations of MPEG-1 and MPEG-2 video, including shareware,
  24.  * are subject to royalty fees to patent holders.  Many of these patents are
  25.  * general enough such that they are unavoidable regardless of implementation
  26.  * design.
  27.  *
  28.  */
  29.  
  30. #include <stdio.h>
  31.  
  32. #include "config.h"
  33. #include "global.h"
  34. //#include "getbits.h"
  35.  
  36. char szTemp[256];
  37. /* private prototypes */
  38. static void decode_motion_vector _ANSI_ARGS_((int *pred, int r_size, int motion_code,
  39.   int motion_residualesidual, int full_pel_vector));
  40.  
  41. /* ISO/IEC 13818-2 sections 6.2.5.2, 6.3.17.2, and 7.6.3: Motion vectors */
  42. void motion_vectors(PMV,dmvector,
  43.   motion_vertical_field_select,s,motion_vector_count,mv_format,h_r_size,v_r_size,dmv,mvscale)
  44. int PMV[2][2][2];
  45. int dmvector[2];
  46. int motion_vertical_field_select[2][2];
  47. int s, motion_vector_count, mv_format, h_r_size, v_r_size, dmv, mvscale;
  48. {
  49.   if (motion_vector_count==1)
  50.   {
  51.     if (mv_format==MV_FIELD && !dmv)
  52.     {
  53.       motion_vertical_field_select[1][s] = motion_vertical_field_select[0][s] = Get_Bits(1);
  54. #ifdef TRACE
  55.       if (Trace_Flag)
  56.       {
  57.         sprintf(szTemp, "motion_vertical_field_select[][%d] (%d): %d\n",s,
  58.           motion_vertical_field_select[0][s],motion_vertical_field_select[0][s]);
  59.     //    ShowInConsole( szTemp );
  60.       }
  61. #endif /* TRACE */
  62.     }
  63.  
  64.     motion_vector(PMV[0][s],dmvector,h_r_size,v_r_size,dmv,mvscale,0);
  65.  
  66.     /* update other motion vector predictors */
  67.     PMV[1][s][0] = PMV[0][s][0];
  68.     PMV[1][s][1] = PMV[0][s][1];
  69.   }
  70.   else
  71.   {
  72.     motion_vertical_field_select[0][s] = Get_Bits(1);
  73. #ifdef TRACE
  74.     if (Trace_Flag)
  75.     {
  76.       sprintf(szTemp, "motion_vertical_field_select[0][%d] (%d): %d\n",s,
  77.         motion_vertical_field_select[0][s],motion_vertical_field_select[0][s]);
  78.               //ShowInConsole( szTemp );
  79.     }
  80. #endif /* TRACE */
  81.     motion_vector(PMV[0][s],dmvector,h_r_size,v_r_size,dmv,mvscale,0);
  82.  
  83.     motion_vertical_field_select[1][s] = Get_Bits(1);
  84. #ifdef TRACE
  85.     if (Trace_Flag)
  86.     {
  87.       sprintf(szTemp, "motion_vertical_field_select[1][%d] (%d): %d\n",s,
  88.         motion_vertical_field_select[1][s],motion_vertical_field_select[1][s]);
  89.               //ShowInConsole( szTemp );
  90.     }
  91. #endif /* TRACE */
  92.     motion_vector(PMV[1][s],dmvector,h_r_size,v_r_size,dmv,mvscale,0);
  93.   }
  94. }
  95.  
  96. /* get and decode motion vector and differential motion vector 
  97.    for one prediction */
  98. void motion_vector(PMV,dmvector,
  99.   h_r_size,v_r_size,dmv,mvscale,full_pel_vector)
  100. int *PMV;
  101. int *dmvector;
  102. int h_r_size;
  103. int v_r_size;
  104. int dmv; /* MPEG-2 only: get differential motion vectors */
  105. int mvscale; /* MPEG-2 only: field vector in frame pic */
  106. int full_pel_vector; /* MPEG-1 only */
  107. {
  108.   int motion_code, motion_residual;
  109.  
  110.   /* horizontal component */
  111.   /* ISO/IEC 13818-2 Table B-10 */
  112.   motion_code = Get_motion_code();
  113.  
  114.   motion_residual = (h_r_size!=0 && motion_code!=0) ? Get_Bits(h_r_size) : 0;
  115.  
  116. #ifdef TRACE
  117.   if (Trace_Flag)
  118.   {
  119.     if (h_r_size!=0 && motion_code!=0)
  120.     {
  121.       printf("motion_residual (");
  122.       Print_Bits(motion_residual,h_r_size,h_r_size);
  123.       printf("): %d\n",motion_residual);
  124.     }
  125.   }
  126. #endif /* TRACE */
  127.  
  128.  
  129.   decode_motion_vector(&PMV[0],h_r_size,motion_code,motion_residual,full_pel_vector);
  130.  
  131.   if (dmv)
  132.     dmvector[0] = Get_dmvector();
  133.  
  134.  
  135.   /* vertical component */
  136.   motion_code     = Get_motion_code();
  137.   motion_residual = (v_r_size!=0 && motion_code!=0) ? Get_Bits(v_r_size) : 0;
  138.  
  139. #ifdef TRACE
  140.   if (Trace_Flag)
  141.   {
  142.     if (v_r_size!=0 && motion_code!=0)
  143.     {
  144.       printf("motion_residual (");
  145.       Print_Bits(motion_residual,v_r_size,v_r_size);
  146.       printf("): %d\n",motion_residual);
  147.     }
  148.   }
  149. #endif /* TRACE */
  150.  
  151.   if (mvscale)
  152.     PMV[1] >>= 1; /* DIV 2 */
  153.  
  154.   decode_motion_vector(&PMV[1],v_r_size,motion_code,motion_residual,full_pel_vector);
  155.  
  156.   if (mvscale)
  157.     PMV[1] <<= 1;
  158.  
  159.   if (dmv)
  160.     dmvector[1] = Get_dmvector();
  161.  
  162. #ifdef TRACE
  163.   if (Trace_Flag)
  164.     printf("PMV = %d,%d\n",PMV[0],PMV[1]);
  165. #endif /* TRACE */
  166. }
  167.  
  168. /* calculate motion vector component */
  169. /* ISO/IEC 13818-2 section 7.6.3.1: Decoding the motion vectors */
  170. /* Note: the arithmetic here is more elegant than that which is shown 
  171.    in 7.6.3.1.  The end results (PMV[][][]) should, however, be the same.  */
  172.  
  173. static void decode_motion_vector(pred,r_size,motion_code,motion_residual,full_pel_vector)
  174. int *pred;
  175. int r_size, motion_code, motion_residual;
  176. int full_pel_vector; /* MPEG-1 (ISO/IEC 11172-1) support */
  177. {
  178.   int lim, vec;
  179.  
  180.   lim = 16<<r_size;
  181.   vec = full_pel_vector ? (*pred >> 1) : (*pred);
  182.  
  183.   if (motion_code!=0)
  184.   {
  185.       if (motion_code>0)
  186.       {
  187.         vec+= ((motion_code-1)<<r_size) + motion_residual + 1;
  188.         if (vec>=lim)
  189.           vec-= lim + lim;
  190.       }
  191.       else
  192.       {
  193.         vec-= ((-motion_code-1)<<r_size) + motion_residual + 1;
  194.         if (vec<-lim)
  195.           vec+= lim + lim;
  196.       }
  197.   }
  198.   *pred = full_pel_vector ? (vec<<1) : vec;
  199. }
  200.  
  201.  
  202. /* ISO/IEC 13818-2 section 7.6.3.6: Dual prime additional arithmetic */
  203. void Dual_Prime_Arithmetic(DMV,dmvector,mvx,mvy)
  204. int DMV[][2];
  205. int *dmvector; /* differential motion vector */
  206. int mvx, mvy;  /* decoded mv components (always in field format) */
  207. {
  208.   if (picture_structure==FRAME_PICTURE)
  209.   {
  210.     if (top_field_first)
  211.     {
  212.       /* vector for prediction of top field from bottom field */
  213.       DMV[0][0] = ((mvx  +(mvx>0))>>1) + dmvector[0];
  214.       DMV[0][1] = ((mvy  +(mvy>0))>>1) + dmvector[1] - 1;
  215.  
  216.       /* vector for prediction of bottom field from top field */
  217.       DMV[1][0] = ((3*mvx+(mvx>0))>>1) + dmvector[0];
  218.       DMV[1][1] = ((3*mvy+(mvy>0))>>1) + dmvector[1] + 1;
  219.     }
  220.     else
  221.     {
  222.       /* vector for prediction of top field from bottom field */
  223.       DMV[0][0] = ((3*mvx+(mvx>0))>>1) + dmvector[0];
  224.       DMV[0][1] = ((3*mvy+(mvy>0))>>1) + dmvector[1] - 1;
  225.  
  226.       /* vector for prediction of bottom field from top field */
  227.       DMV[1][0] = ((mvx  +(mvx>0))>>1) + dmvector[0];
  228.       DMV[1][1] = ((mvy  +(mvy>0))>>1) + dmvector[1] + 1;
  229.     }
  230.   }
  231.   else
  232.   {
  233.     /* vector for prediction from field of opposite 'parity' */
  234.     DMV[0][0] = ((mvx+(mvx>0))>>1) + dmvector[0];
  235.     DMV[0][1] = ((mvy+(mvy>0))>>1) + dmvector[1];
  236.  
  237.     /* correct for vertical field shift */
  238.     if (picture_structure==TOP_FIELD)
  239.       DMV[0][1]--;
  240.     else
  241.       DMV[0][1]++;
  242.   }
  243. }
  244.  
  245.