annotate ffmpeg/libavcodec/eamad.c @ 13:844d341cf643 tip

Back up before ISMIR
author Yading Song <yading.song@eecs.qmul.ac.uk>
date Thu, 31 Oct 2013 13:17:06 +0000
parents 6840f77b83aa
children
rev   line source
yading@10 1 /*
yading@10 2 * Electronic Arts Madcow Video Decoder
yading@10 3 * Copyright (c) 2007-2009 Peter Ross
yading@10 4 *
yading@10 5 * This file is part of FFmpeg.
yading@10 6 *
yading@10 7 * FFmpeg is free software; you can redistribute it and/or
yading@10 8 * modify it under the terms of the GNU Lesser General Public
yading@10 9 * License as published by the Free Software Foundation; either
yading@10 10 * version 2.1 of the License, or (at your option) any later version.
yading@10 11 *
yading@10 12 * FFmpeg is distributed in the hope that it will be useful,
yading@10 13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
yading@10 15 * Lesser General Public License for more details.
yading@10 16 *
yading@10 17 * You should have received a copy of the GNU Lesser General Public
yading@10 18 * License along with FFmpeg; if not, write to the Free Software
yading@10 19 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
yading@10 20 */
yading@10 21
yading@10 22 /**
yading@10 23 * @file
yading@10 24 * Electronic Arts Madcow Video Decoder
yading@10 25 * @author Peter Ross <pross@xvid.org>
yading@10 26 *
yading@10 27 * @see technical details at
yading@10 28 * http://wiki.multimedia.cx/index.php?title=Electronic_Arts_MAD
yading@10 29 */
yading@10 30
yading@10 31 #include "avcodec.h"
yading@10 32 #include "get_bits.h"
yading@10 33 #include "aandcttab.h"
yading@10 34 #include "eaidct.h"
yading@10 35 #include "internal.h"
yading@10 36 #include "mpeg12.h"
yading@10 37 #include "mpeg12data.h"
yading@10 38 #include "libavutil/imgutils.h"
yading@10 39
yading@10 40 #define EA_PREAMBLE_SIZE 8
yading@10 41 #define MADk_TAG MKTAG('M', 'A', 'D', 'k') /* MAD i-frame */
yading@10 42 #define MADm_TAG MKTAG('M', 'A', 'D', 'm') /* MAD p-frame */
yading@10 43 #define MADe_TAG MKTAG('M', 'A', 'D', 'e') /* MAD lqp-frame */
yading@10 44
yading@10 45 typedef struct MadContext {
yading@10 46 AVCodecContext *avctx;
yading@10 47 DSPContext dsp;
yading@10 48 AVFrame last_frame;
yading@10 49 GetBitContext gb;
yading@10 50 void *bitstream_buf;
yading@10 51 unsigned int bitstream_buf_size;
yading@10 52 DECLARE_ALIGNED(16, int16_t, block)[64];
yading@10 53 ScanTable scantable;
yading@10 54 uint16_t quant_matrix[64];
yading@10 55 int mb_x;
yading@10 56 int mb_y;
yading@10 57 } MadContext;
yading@10 58
yading@10 59 static av_cold int decode_init(AVCodecContext *avctx)
yading@10 60 {
yading@10 61 MadContext *s = avctx->priv_data;
yading@10 62 s->avctx = avctx;
yading@10 63 avctx->pix_fmt = AV_PIX_FMT_YUV420P;
yading@10 64 ff_dsputil_init(&s->dsp, avctx);
yading@10 65 ff_init_scantable_permutation(s->dsp.idct_permutation, FF_NO_IDCT_PERM);
yading@10 66 ff_init_scantable(s->dsp.idct_permutation, &s->scantable, ff_zigzag_direct);
yading@10 67 ff_mpeg12_init_vlcs();
yading@10 68 return 0;
yading@10 69 }
yading@10 70
yading@10 71 static inline void comp(unsigned char *dst, int dst_stride,
yading@10 72 unsigned char *src, int src_stride, int add)
yading@10 73 {
yading@10 74 int j, i;
yading@10 75 for (j=0; j<8; j++)
yading@10 76 for (i=0; i<8; i++)
yading@10 77 dst[j*dst_stride + i] = av_clip_uint8(src[j*src_stride + i] + add);
yading@10 78 }
yading@10 79
yading@10 80 static inline void comp_block(MadContext *t, AVFrame *frame,
yading@10 81 int mb_x, int mb_y,
yading@10 82 int j, int mv_x, int mv_y, int add)
yading@10 83 {
yading@10 84 if (j < 4) {
yading@10 85 unsigned offset = (mb_y*16 + ((j&2)<<2) + mv_y)*t->last_frame.linesize[0] + mb_x*16 + ((j&1)<<3) + mv_x;
yading@10 86 if (offset >= (t->avctx->height - 7) * t->last_frame.linesize[0] - 7)
yading@10 87 return;
yading@10 88 comp(frame->data[0] + (mb_y*16 + ((j&2)<<2))*frame->linesize[0] + mb_x*16 + ((j&1)<<3),
yading@10 89 frame->linesize[0],
yading@10 90 t->last_frame.data[0] + offset,
yading@10 91 t->last_frame.linesize[0], add);
yading@10 92 } else if (!(t->avctx->flags & CODEC_FLAG_GRAY)) {
yading@10 93 int index = j - 3;
yading@10 94 unsigned offset = (mb_y * 8 + (mv_y/2))*t->last_frame.linesize[index] + mb_x * 8 + (mv_x/2);
yading@10 95 if (offset >= (t->avctx->height/2 - 7) * t->last_frame.linesize[index] - 7)
yading@10 96 return;
yading@10 97 comp(frame->data[index] + (mb_y*8)*frame->linesize[index] + mb_x * 8,
yading@10 98 frame->linesize[index],
yading@10 99 t->last_frame.data[index] + offset,
yading@10 100 t->last_frame.linesize[index], add);
yading@10 101 }
yading@10 102 }
yading@10 103
yading@10 104 static inline void idct_put(MadContext *t, AVFrame *frame, int16_t *block,
yading@10 105 int mb_x, int mb_y, int j)
yading@10 106 {
yading@10 107 if (j < 4) {
yading@10 108 ff_ea_idct_put_c(
yading@10 109 frame->data[0] + (mb_y*16 + ((j&2)<<2))*frame->linesize[0] + mb_x*16 + ((j&1)<<3),
yading@10 110 frame->linesize[0], block);
yading@10 111 } else if (!(t->avctx->flags & CODEC_FLAG_GRAY)) {
yading@10 112 int index = j - 3;
yading@10 113 ff_ea_idct_put_c(
yading@10 114 frame->data[index] + (mb_y*8)*frame->linesize[index] + mb_x*8,
yading@10 115 frame->linesize[index], block);
yading@10 116 }
yading@10 117 }
yading@10 118
yading@10 119 static inline int decode_block_intra(MadContext *s, int16_t * block)
yading@10 120 {
yading@10 121 int level, i, j, run;
yading@10 122 RLTable *rl = &ff_rl_mpeg1;
yading@10 123 const uint8_t *scantable = s->scantable.permutated;
yading@10 124 int16_t *quant_matrix = s->quant_matrix;
yading@10 125
yading@10 126 block[0] = (128 + get_sbits(&s->gb, 8)) * quant_matrix[0];
yading@10 127
yading@10 128 /* The RL decoder is derived from mpeg1_decode_block_intra;
yading@10 129 Escaped level and run values a decoded differently */
yading@10 130 i = 0;
yading@10 131 {
yading@10 132 OPEN_READER(re, &s->gb);
yading@10 133 /* now quantify & encode AC coefficients */
yading@10 134 for (;;) {
yading@10 135 UPDATE_CACHE(re, &s->gb);
yading@10 136 GET_RL_VLC(level, run, re, &s->gb, rl->rl_vlc[0], TEX_VLC_BITS, 2, 0);
yading@10 137
yading@10 138 if (level == 127) {
yading@10 139 break;
yading@10 140 } else if (level != 0) {
yading@10 141 i += run;
yading@10 142 j = scantable[i];
yading@10 143 level = (level*quant_matrix[j]) >> 4;
yading@10 144 level = (level-1)|1;
yading@10 145 level = (level ^ SHOW_SBITS(re, &s->gb, 1)) - SHOW_SBITS(re, &s->gb, 1);
yading@10 146 LAST_SKIP_BITS(re, &s->gb, 1);
yading@10 147 } else {
yading@10 148 /* escape */
yading@10 149 UPDATE_CACHE(re, &s->gb);
yading@10 150 level = SHOW_SBITS(re, &s->gb, 10); SKIP_BITS(re, &s->gb, 10);
yading@10 151
yading@10 152 UPDATE_CACHE(re, &s->gb);
yading@10 153 run = SHOW_UBITS(re, &s->gb, 6)+1; LAST_SKIP_BITS(re, &s->gb, 6);
yading@10 154
yading@10 155 i += run;
yading@10 156 j = scantable[i];
yading@10 157 if (level < 0) {
yading@10 158 level = -level;
yading@10 159 level = (level*quant_matrix[j]) >> 4;
yading@10 160 level = (level-1)|1;
yading@10 161 level = -level;
yading@10 162 } else {
yading@10 163 level = (level*quant_matrix[j]) >> 4;
yading@10 164 level = (level-1)|1;
yading@10 165 }
yading@10 166 }
yading@10 167 if (i > 63) {
yading@10 168 av_log(s->avctx, AV_LOG_ERROR, "ac-tex damaged at %d %d\n", s->mb_x, s->mb_y);
yading@10 169 return -1;
yading@10 170 }
yading@10 171
yading@10 172 block[j] = level;
yading@10 173 }
yading@10 174 CLOSE_READER(re, &s->gb);
yading@10 175 }
yading@10 176 return 0;
yading@10 177 }
yading@10 178
yading@10 179 static int decode_motion(GetBitContext *gb)
yading@10 180 {
yading@10 181 int value = 0;
yading@10 182 if (get_bits1(gb)) {
yading@10 183 if (get_bits1(gb))
yading@10 184 value = -17;
yading@10 185 value += get_bits(gb, 4) + 1;
yading@10 186 }
yading@10 187 return value;
yading@10 188 }
yading@10 189
yading@10 190 static int decode_mb(MadContext *s, AVFrame *frame, int inter)
yading@10 191 {
yading@10 192 int mv_map = 0;
yading@10 193 int mv_x, mv_y;
yading@10 194 int j;
yading@10 195
yading@10 196 if (inter) {
yading@10 197 int v = decode210(&s->gb);
yading@10 198 if (v < 2) {
yading@10 199 mv_map = v ? get_bits(&s->gb, 6) : 63;
yading@10 200 mv_x = decode_motion(&s->gb);
yading@10 201 mv_y = decode_motion(&s->gb);
yading@10 202 }
yading@10 203 }
yading@10 204
yading@10 205 for (j=0; j<6; j++) {
yading@10 206 if (mv_map & (1<<j)) { // mv_x and mv_y are guarded by mv_map
yading@10 207 int add = 2*decode_motion(&s->gb);
yading@10 208 if (s->last_frame.data[0])
yading@10 209 comp_block(s, frame, s->mb_x, s->mb_y, j, mv_x, mv_y, add);
yading@10 210 } else {
yading@10 211 s->dsp.clear_block(s->block);
yading@10 212 if(decode_block_intra(s, s->block) < 0)
yading@10 213 return -1;
yading@10 214 idct_put(s, frame, s->block, s->mb_x, s->mb_y, j);
yading@10 215 }
yading@10 216 }
yading@10 217 return 0;
yading@10 218 }
yading@10 219
yading@10 220 static void calc_quant_matrix(MadContext *s, int qscale)
yading@10 221 {
yading@10 222 int i;
yading@10 223
yading@10 224 s->quant_matrix[0] = (ff_inv_aanscales[0]*ff_mpeg1_default_intra_matrix[0]) >> 11;
yading@10 225 for (i=1; i<64; i++)
yading@10 226 s->quant_matrix[i] = (ff_inv_aanscales[i]*ff_mpeg1_default_intra_matrix[i]*qscale + 32) >> 10;
yading@10 227 }
yading@10 228
yading@10 229 static int decode_frame(AVCodecContext *avctx,
yading@10 230 void *data, int *got_frame,
yading@10 231 AVPacket *avpkt)
yading@10 232 {
yading@10 233 const uint8_t *buf = avpkt->data;
yading@10 234 int buf_size = avpkt->size;
yading@10 235 const uint8_t *buf_end = buf+buf_size;
yading@10 236 MadContext *s = avctx->priv_data;
yading@10 237 AVFrame *frame = data;
yading@10 238 int width, height;
yading@10 239 int chunk_type;
yading@10 240 int inter, ret;
yading@10 241
yading@10 242 if (buf_size < 26) {
yading@10 243 av_log(avctx, AV_LOG_ERROR, "Input buffer too small\n");
yading@10 244 *got_frame = 0;
yading@10 245 return AVERROR_INVALIDDATA;
yading@10 246 }
yading@10 247
yading@10 248 chunk_type = AV_RL32(&buf[0]);
yading@10 249 inter = (chunk_type == MADm_TAG || chunk_type == MADe_TAG);
yading@10 250 buf += 8;
yading@10 251
yading@10 252 av_reduce(&avctx->time_base.num, &avctx->time_base.den,
yading@10 253 AV_RL16(&buf[6]), 1000, 1<<30);
yading@10 254
yading@10 255 width = AV_RL16(&buf[8]);
yading@10 256 height = AV_RL16(&buf[10]);
yading@10 257 calc_quant_matrix(s, buf[13]);
yading@10 258 buf += 16;
yading@10 259
yading@10 260 if (avctx->width != width || avctx->height != height) {
yading@10 261 if((width * height)/2048*7 > buf_end-buf)
yading@10 262 return AVERROR_INVALIDDATA;
yading@10 263 if ((ret = av_image_check_size(width, height, 0, avctx)) < 0)
yading@10 264 return ret;
yading@10 265 avcodec_set_dimensions(avctx, width, height);
yading@10 266 av_frame_unref(&s->last_frame);
yading@10 267 }
yading@10 268
yading@10 269 if ((ret = ff_get_buffer(avctx, frame, AV_GET_BUFFER_FLAG_REF)) < 0)
yading@10 270 return ret;
yading@10 271
yading@10 272 if (inter && !s->last_frame.data[0]) {
yading@10 273 av_log(avctx, AV_LOG_WARNING, "Missing reference frame.\n");
yading@10 274 ret = ff_get_buffer(avctx, &s->last_frame, AV_GET_BUFFER_FLAG_REF);
yading@10 275 if (ret < 0)
yading@10 276 return ret;
yading@10 277 memset(s->last_frame.data[0], 0, s->last_frame.height *
yading@10 278 s->last_frame.linesize[0]);
yading@10 279 memset(s->last_frame.data[1], 0x80, s->last_frame.height / 2 *
yading@10 280 s->last_frame.linesize[1]);
yading@10 281 memset(s->last_frame.data[2], 0x80, s->last_frame.height / 2 *
yading@10 282 s->last_frame.linesize[2]);
yading@10 283 }
yading@10 284
yading@10 285 av_fast_padded_malloc(&s->bitstream_buf, &s->bitstream_buf_size,
yading@10 286 buf_end - buf);
yading@10 287 if (!s->bitstream_buf)
yading@10 288 return AVERROR(ENOMEM);
yading@10 289 s->dsp.bswap16_buf(s->bitstream_buf, (const uint16_t*)buf, (buf_end-buf)/2);
yading@10 290 memset((uint8_t*)s->bitstream_buf + (buf_end-buf), 0, FF_INPUT_BUFFER_PADDING_SIZE);
yading@10 291 init_get_bits(&s->gb, s->bitstream_buf, 8*(buf_end-buf));
yading@10 292
yading@10 293 for (s->mb_y=0; s->mb_y < (avctx->height+15)/16; s->mb_y++)
yading@10 294 for (s->mb_x=0; s->mb_x < (avctx->width +15)/16; s->mb_x++)
yading@10 295 if(decode_mb(s, frame, inter) < 0)
yading@10 296 return AVERROR_INVALIDDATA;
yading@10 297
yading@10 298 *got_frame = 1;
yading@10 299
yading@10 300 if (chunk_type != MADe_TAG) {
yading@10 301 av_frame_unref(&s->last_frame);
yading@10 302 if ((ret = av_frame_ref(&s->last_frame, frame)) < 0)
yading@10 303 return ret;
yading@10 304 }
yading@10 305
yading@10 306 return buf_size;
yading@10 307 }
yading@10 308
yading@10 309 static av_cold int decode_end(AVCodecContext *avctx)
yading@10 310 {
yading@10 311 MadContext *t = avctx->priv_data;
yading@10 312 av_frame_unref(&t->last_frame);
yading@10 313 av_free(t->bitstream_buf);
yading@10 314 return 0;
yading@10 315 }
yading@10 316
yading@10 317 AVCodec ff_eamad_decoder = {
yading@10 318 .name = "eamad",
yading@10 319 .type = AVMEDIA_TYPE_VIDEO,
yading@10 320 .id = AV_CODEC_ID_MAD,
yading@10 321 .priv_data_size = sizeof(MadContext),
yading@10 322 .init = decode_init,
yading@10 323 .close = decode_end,
yading@10 324 .decode = decode_frame,
yading@10 325 .capabilities = CODEC_CAP_DR1,
yading@10 326 .long_name = NULL_IF_CONFIG_SMALL("Electronic Arts Madcow Video")
yading@10 327 };