--- branches/dev-api-4/xvidcore/src/motion/motion_est.h 2003/06/09 01:25:19 1053 +++ branches/dev-api-4/xvidcore/src/motion/motion_est.h 2003/08/03 10:10:54 1108 @@ -1,34 +1,29 @@ -/************************************************************************** +/***************************************************************************** * * XVID MPEG-4 VIDEO CODEC - * - Motion estimation header - + * - Motion Estimation related header - * - * This program is an implementation of a part of one or more MPEG-4 - * Video tools as specified in ISO/IEC 14496-2 standard. Those intending - * to use this software module in hardware or software products are - * advised that its use may infringe existing patents or copyrights, and - * any such use would be at such party's own risk. The original - * developer of this software module and his/her company, and subsequent - * editors and their companies, will have no liability for use of this - * software or modifications or derivatives thereof. + * Copyright(C) 2002 Christoph Lampert + * 2002 Michael Militzer + * 2002-2003 Radoslaw Czyz * - * This program is free software; you can redistribute it and/or modify + * This program is free software ; you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation ; either version 2 of the License, or * (at your option) any later version. * * This program is distributed in the hope that it will be useful, - * but WITHOUT ANY WARRANTY; without even the implied warranty of + * but WITHOUT ANY WARRANTY ; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU General Public License - * along with this program; if not, write to the Free Software + * along with this program ; if not, write to the Free Software * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA * - * $Id: motion_est.h,v 1.3.2.8 2003-06-09 01:22:18 edgomez Exp $ + * $Id: motion_est.h,v 1.3.2.15 2003-08-03 10:10:54 syskin Exp $ * - ***************************************************************************/ + ****************************************************************************/ #ifndef _MOTION_EST_H_ #define _MOTION_EST_H_ @@ -144,6 +139,7 @@ int16_t * dctSpace; uint32_t iQuant; uint32_t quant_type; + int * cbp; } SearchData; @@ -213,6 +209,7 @@ const IMAGE * const pRefH, const IMAGE * const pRefV, const IMAGE * const pRefHV, + const IMAGE * const pGMC, const uint32_t iLimit); static void @@ -232,8 +229,30 @@ const MACROBLOCK * const prevMBs, MACROBLOCK * const pMB); -static WARPPOINTS -GlobalMotionEst(const MACROBLOCK * const pMBs, +static __inline void +GMEanalyzeMB (const uint8_t * const pCur, + const uint8_t * const pRef, + const uint8_t * const pRefH, + const uint8_t * const pRefV, + const uint8_t * const pRefHV, + const int x, + const int y, + const MBParam * const pParam, + MACROBLOCK * const pMBs, + SearchData * const Data); + +void +GMEanalysis(const MBParam * const pParam, + const FRAMEINFO * const current, + const FRAMEINFO * const reference, + const IMAGE * const pRefH, + const IMAGE * const pRefV, + const IMAGE * const pRefHV); + + + +WARPPOINTS +GlobalMotionEst(MACROBLOCK * const pMBs, const MBParam * const pParam, const FRAMEINFO * const current, const FRAMEINFO * const reference, @@ -241,6 +260,31 @@ const IMAGE * const pRefV, const IMAGE * const pRefHV ); + +int +GlobalMotionEstRefine(WARPPOINTS *const startwp, + MACROBLOCK * const pMBs, + const MBParam * const pParam, + const FRAMEINFO * const current, + const FRAMEINFO * const reference, + const IMAGE * const pCurr, + const IMAGE * const pRef, + const IMAGE * const pRefH, + const IMAGE * const pRefV, + const IMAGE * const pRefHV); + + + + +int +globalSAD(const WARPPOINTS *const wp, + const MBParam * const pParam, + const MACROBLOCK * const pMBs, + const FRAMEINFO * const current, + const IMAGE * const pRef, + const IMAGE * const pCurr, + uint8_t *const GMCblock); + #define iDiamondSize 2 static __inline uint32_t @@ -248,39 +292,39 @@ { uint32_t Flags = MotionFlags; - if (!(VopFlags & XVID_VOP_MODEDECISION_BITS)) - Flags &= ~(XVID_ME_QUARTERPELREFINE16_BITS+XVID_ME_QUARTERPELREFINE8_BITS+XVID_ME_HALFPELREFINE16_BITS+XVID_ME_HALFPELREFINE8_BITS+XVID_ME_EXTSEARCH_BITS); + if (!(VopFlags & XVID_VOP_MODEDECISION_RD)) + Flags &= ~(XVID_ME_QUARTERPELREFINE16_RD+XVID_ME_QUARTERPELREFINE8_RD+XVID_ME_HALFPELREFINE16_RD+XVID_ME_HALFPELREFINE8_RD+XVID_ME_EXTSEARCH_RD); - if (Flags & XVID_ME_EXTSEARCH_BITS) - Flags |= XVID_ME_HALFPELREFINE16_BITS; + if (Flags & XVID_ME_EXTSEARCH_RD) + Flags |= XVID_ME_HALFPELREFINE16_RD; - if (Flags & XVID_ME_EXTSEARCH_BITS && MotionFlags & XVID_ME_EXTSEARCH8) - Flags |= XVID_ME_HALFPELREFINE8_BITS; + if (Flags & XVID_ME_EXTSEARCH_RD && MotionFlags & XVID_ME_EXTSEARCH8) + Flags |= XVID_ME_HALFPELREFINE8_RD; - if (Flags & XVID_ME_HALFPELREFINE16_BITS) - Flags |= XVID_ME_QUARTERPELREFINE16_BITS; + if (Flags & XVID_ME_HALFPELREFINE16_RD) + Flags |= XVID_ME_QUARTERPELREFINE16_RD; - if (Flags & XVID_ME_HALFPELREFINE8_BITS) { - Flags |= XVID_ME_QUARTERPELREFINE8_BITS; + if (Flags & XVID_ME_HALFPELREFINE8_RD) { + Flags |= XVID_ME_QUARTERPELREFINE8_RD; Flags &= ~XVID_ME_HALFPELREFINE8; } - if (Flags & XVID_ME_QUARTERPELREFINE8_BITS) + if (Flags & XVID_ME_QUARTERPELREFINE8_RD) Flags &= ~XVID_ME_QUARTERPELREFINE8; if (!(VolFlags & XVID_VOL_QUARTERPEL)) - Flags &= ~(XVID_ME_QUARTERPELREFINE16+XVID_ME_QUARTERPELREFINE8+XVID_ME_QUARTERPELREFINE16_BITS+XVID_ME_QUARTERPELREFINE8_BITS); + Flags &= ~(XVID_ME_QUARTERPELREFINE16+XVID_ME_QUARTERPELREFINE8+XVID_ME_QUARTERPELREFINE16_RD+XVID_ME_QUARTERPELREFINE8_RD); if (!(VopFlags & XVID_VOP_HALFPEL)) - Flags &= ~(XVID_ME_EXTSEARCH16+XVID_ME_HALFPELREFINE16+XVID_ME_HALFPELREFINE8+XVID_ME_HALFPELREFINE16_BITS+XVID_ME_HALFPELREFINE8_BITS); + Flags &= ~(XVID_ME_EXTSEARCH16+XVID_ME_HALFPELREFINE16+XVID_ME_HALFPELREFINE8+XVID_ME_HALFPELREFINE16_RD+XVID_ME_HALFPELREFINE8_RD); if ((VopFlags & XVID_VOP_GREYSCALE) || (VopFlags & XVID_VOP_REDUCED)) - Flags &= ~(XVID_ME_CHROMA16 + XVID_ME_CHROMA8); + Flags &= ~(XVID_ME_CHROMA_PVOP + XVID_ME_CHROMA_BVOP); return Flags; } -/* BITS mode decision and search */ +/* RD mode decision and search */ #include "../bitstream/zigzag.h" #include "../quant/quant_mpeg4.h" @@ -289,79 +333,85 @@ #include "../dct/fdct.h" static int -CountMBBitsInter(SearchData * const Data, - const MACROBLOCK * const pMBs, const int x, const int y, - const MBParam * const pParam, - const uint32_t MotionFlags); +findRDinter(SearchData * const Data, + const MACROBLOCK * const pMBs, const int x, const int y, + const MBParam * const pParam, + const uint32_t MotionFlags); static int -CountMBBitsInter4v(const SearchData * const Data, - MACROBLOCK * const pMB, const MACROBLOCK * const pMBs, - const int x, const int y, - const MBParam * const pParam, const uint32_t MotionFlags, - const VECTOR * const backup); +findRDinter4v(const SearchData * const Data, + MACROBLOCK * const pMB, const MACROBLOCK * const pMBs, + const int x, const int y, + const MBParam * const pParam, const uint32_t MotionFlags, + const VECTOR * const backup); static int -CountMBBitsIntra(const SearchData * const Data); +findRDintra(const SearchData * const Data); + +static int +findRDgmc(const SearchData * const Data, const IMAGE * const vGMC, const int x, const int y); int CodeCoeffIntra_CalcBits(const int16_t qcoeff[64], const uint16_t * zigzag); int CodeCoeffInter_CalcBits(const int16_t qcoeff[64], const uint16_t * zigzag); -#define LAMBDA ( (int)(1.0*BITS_MULT) ) +/* one over lambda for R-D mode decision and motion search */ +#define LAMBDA ( (int)(BITS_MULT/1.0) ) static __inline unsigned int Block_CalcBits( int16_t * const coeff, int16_t * const data, + int16_t * const dqcoeff, const uint32_t quant, const int quant_type, uint32_t * cbp, const int block) { int sum; int bits; - const int lambda = LAMBDA*quant*quant; int distortion = 0; int i; fdct(data); - if (quant_type == 0) sum = quant_inter(coeff, data, quant); + if (quant_type) sum = quant_inter(coeff, data, quant); else sum = quant4_inter(coeff, data, quant); if (sum > 0) { *cbp |= 1 << (5 - block); bits = BITS_MULT * CodeCoeffInter_CalcBits(coeff, scan_tables[0]); - } else bits = 0; - if (quant_type == 0) dequant_inter(coeff, coeff, quant); - else dequant4_inter(coeff, coeff, quant); + if (quant_type) dequant_inter(dqcoeff, coeff, quant); + else dequant4_inter(dqcoeff, coeff, quant); - for (i = 0; i < 64; i++) { - distortion += (data[i] - coeff[i])*(data[i] - coeff[i]); - } - bits += (BITS_MULT*BITS_MULT*distortion)/lambda; + for (i = 0; i < 64; i++) + distortion += (data[i] - dqcoeff[i])*(data[i] - dqcoeff[i]); + } else { + bits = 0; + for (i = 0; i < 64; i++) + distortion += data[i]*data[i]; + } - return bits; + return bits + (LAMBDA*distortion)/(quant*quant); } static __inline unsigned int Block_CalcBitsIntra(int16_t * const coeff, int16_t * const data, + int16_t * const dqcoeff, const uint32_t quant, const int quant_type, uint32_t * cbp, const int block, int * dcpred) { int bits, i; - const int lambda = LAMBDA*quant*quant; int distortion = 0; - uint32_t iDcScaler = get_dc_scaler(quant, block > 3); + uint32_t iDcScaler = get_dc_scaler(quant, block < 4); int b_dc; fdct(data); data[0] -= 1024; - if (quant_type == 0) quant_intra(coeff, data, quant, iDcScaler); + if (quant_type) quant_intra(coeff, data, quant, iDcScaler); else quant4_intra(coeff, data, quant, iDcScaler); b_dc = coeff[0]; @@ -370,22 +420,20 @@ *dcpred = b_dc; } - *cbp |= 1 << (5 - block); bits = BITS_MULT*CodeCoeffIntra_CalcBits(coeff, scan_tables[0]); - bits += BITS_MULT*dcy_tab[coeff[0] + 255].len; if (bits != 0) *cbp |= 1 << (5 - block); - coeff[0] = b_dc; - if (quant_type == 0) dequant_intra(coeff, coeff, quant, iDcScaler); - else dequant4_intra(coeff, coeff, quant, iDcScaler); + if (block < 4) bits += BITS_MULT*dcy_tab[coeff[0] + 255].len; + else bits += BITS_MULT*dcc_tab[coeff[0] + 255].len; - for (i = 0; i < 64; i++) { - distortion += (data[i] - coeff[i])*(data[i] - coeff[i]); - } + coeff[0] = b_dc; + if (quant_type) dequant_intra(dqcoeff, coeff, quant, iDcScaler); + else dequant4_intra(dqcoeff, coeff, quant, iDcScaler); - bits += (BITS_MULT*BITS_MULT*distortion)/lambda; + for (i = 0; i < 64; i++) + distortion += (data[i] - dqcoeff[i])*(data[i] - dqcoeff[i]); - return bits; + return bits + (LAMBDA*distortion)/(quant*quant); } #endif /* _MOTION_EST_H_ */