truemotion1.c 27.9 KB
Newer Older
1 2 3 4
/*
 * Duck TrueMotion 1.0 Decoder
 * Copyright (C) 2003 Alex Beregszaszi & Mike Melanson
 *
5 6 7
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
8 9
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
10
 * version 2.1 of the License, or (at your option) any later version.
11
 *
12
 * FFmpeg is distributed in the hope that it will be useful,
13 14 15 16 17
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
18
 * License along with FFmpeg; if not, write to the Free Software
19
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 21 22
 */

/**
23
 * @file
24
 * Duck TrueMotion v1 Video Decoder by
25
 * Alex Beregszaszi and
26 27 28
 * Mike Melanson (melanson@pcisys.net)
 *
 * The TrueMotion v1 decoder presently only decodes 16-bit TM1 data and
29
 * outputs RGB555 (or RGB565) data. 24-bit TM1 data is not supported yet.
30 31 32 33 34 35 36
 */

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include "avcodec.h"
37
#include "internal.h"
38
#include "libavutil/imgutils.h"
39
#include "libavutil/internal.h"
40
#include "libavutil/intreadwrite.h"
41
#include "libavutil/mem.h"
42 43 44 45 46

#include "truemotion1data.h"

typedef struct TrueMotion1Context {
    AVCodecContext *avctx;
47
    AVFrame *frame;
48

Michael Niedermayer's avatar
Michael Niedermayer committed
49
    const uint8_t *buf;
50 51
    int size;

Michael Niedermayer's avatar
Michael Niedermayer committed
52
    const uint8_t *mb_change_bits;
53
    int mb_change_bits_row_size;
Michael Niedermayer's avatar
Michael Niedermayer committed
54
    const uint8_t *index_stream;
55 56 57 58
    int index_stream_size;

    int flags;
    int x, y, w, h;
59

60 61
    uint32_t y_predictor_table[1024];
    uint32_t c_predictor_table[1024];
62 63
    uint32_t fat_y_predictor_table[1024];
    uint32_t fat_c_predictor_table[1024];
64

65 66 67 68 69
    int compression;
    int block_type;
    int block_width;
    int block_height;

70 71 72 73
    int16_t ydt[8];
    int16_t cdt[8];
    int16_t fat_ydt[8];
    int16_t fat_cdt[8];
74

75 76 77
    int last_deltaset, last_vectable;

    unsigned int *vert_pred;
78
    int vert_pred_size;
79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117

} TrueMotion1Context;

#define FLAG_SPRITE         32
#define FLAG_KEYFRAME       16
#define FLAG_INTERFRAME      8
#define FLAG_INTERPOLATED    4

struct frame_header {
    uint8_t header_size;
    uint8_t compression;
    uint8_t deltaset;
    uint8_t vectable;
    uint16_t ysize;
    uint16_t xsize;
    uint16_t checksum;
    uint8_t version;
    uint8_t header_type;
    uint8_t flags;
    uint8_t control;
    uint16_t xoffset;
    uint16_t yoffset;
    uint16_t width;
    uint16_t height;
};

#define ALGO_NOP        0
#define ALGO_RGB16V     1
#define ALGO_RGB16H     2
#define ALGO_RGB24H     3

/* these are the various block sizes that can occupy a 4x4 block */
#define BLOCK_2x2  0
#define BLOCK_2x4  1
#define BLOCK_4x2  2
#define BLOCK_4x4  3

typedef struct comp_types {
    int algorithm;
118 119
    int block_width; // vres
    int block_height; // hres
120 121 122
    int block_type;
} comp_types;

123
/* { valid for metatype }, algorithm, num of deltas, vert res, horiz res */
124
static const comp_types compression_types[17] = {
125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154
    { ALGO_NOP,    0, 0, 0 },

    { ALGO_RGB16V, 4, 4, BLOCK_4x4 },
    { ALGO_RGB16H, 4, 4, BLOCK_4x4 },
    { ALGO_RGB16V, 4, 2, BLOCK_4x2 },
    { ALGO_RGB16H, 4, 2, BLOCK_4x2 },

    { ALGO_RGB16V, 2, 4, BLOCK_2x4 },
    { ALGO_RGB16H, 2, 4, BLOCK_2x4 },
    { ALGO_RGB16V, 2, 2, BLOCK_2x2 },
    { ALGO_RGB16H, 2, 2, BLOCK_2x2 },

    { ALGO_NOP,    4, 4, BLOCK_4x4 },
    { ALGO_RGB24H, 4, 4, BLOCK_4x4 },
    { ALGO_NOP,    4, 2, BLOCK_4x2 },
    { ALGO_RGB24H, 4, 2, BLOCK_4x2 },

    { ALGO_NOP,    2, 4, BLOCK_2x4 },
    { ALGO_RGB24H, 2, 4, BLOCK_2x4 },
    { ALGO_NOP,    2, 2, BLOCK_2x2 },
    { ALGO_RGB24H, 2, 2, BLOCK_2x2 }
};

static void select_delta_tables(TrueMotion1Context *s, int delta_table_index)
{
    int i;

    if (delta_table_index > 3)
        return;

155 156 157 158
    memcpy(s->ydt, ydts[delta_table_index], 8 * sizeof(int16_t));
    memcpy(s->cdt, cdts[delta_table_index], 8 * sizeof(int16_t));
    memcpy(s->fat_ydt, fat_ydts[delta_table_index], 8 * sizeof(int16_t));
    memcpy(s->fat_cdt, fat_cdts[delta_table_index], 8 * sizeof(int16_t));
159 160 161 162 163 164 165 166 167 168 169 170

    /* Y skinny deltas need to be halved for some reason; maybe the
     * skinny Y deltas should be modified */
    for (i = 0; i < 8; i++)
    {
        /* drop the lsb before dividing by 2-- net effect: round down
         * when dividing a negative number (e.g., -3/2 = -2, not -1) */
        s->ydt[i] &= 0xFFFE;
        s->ydt[i] /= 2;
    }
}

171
#if HAVE_BIGENDIAN
172
static int make_ydt15_entry(int p2, int p1, int16_t *ydt)
173
#else
174
static int make_ydt15_entry(int p1, int p2, int16_t *ydt)
175
#endif
176 177
{
    int lo, hi;
178

179 180 181 182
    lo = ydt[p1];
    lo += (lo << 5) + (lo << 10);
    hi = ydt[p2];
    hi += (hi << 5) + (hi << 10);
183
    return (lo + (hi << 16)) << 1;
184 185
}

186
static int make_cdt15_entry(int p1, int p2, int16_t *cdt)
187 188
{
    int r, b, lo;
189

190 191 192
    b = cdt[p2];
    r = cdt[p1] << 10;
    lo = b + r;
193
    return (lo + (lo << 16)) << 1;
194 195
}

196
#if HAVE_BIGENDIAN
197 198 199 200 201 202
static int make_ydt16_entry(int p2, int p1, int16_t *ydt)
#else
static int make_ydt16_entry(int p1, int p2, int16_t *ydt)
#endif
{
    int lo, hi;
203

204 205 206 207
    lo = ydt[p1];
    lo += (lo << 6) + (lo << 11);
    hi = ydt[p2];
    hi += (hi << 6) + (hi << 11);
208
    return (lo + (hi << 16)) << 1;
209 210 211 212 213
}

static int make_cdt16_entry(int p1, int p2, int16_t *cdt)
{
    int r, b, lo;
214

215 216 217
    b = cdt[p2];
    r = cdt[p1] << 11;
    lo = b + r;
218
    return (lo + (lo << 16)) << 1;
219 220 221 222 223
}

static int make_ydt24_entry(int p1, int p2, int16_t *ydt)
{
    int lo, hi;
224

225 226
    lo = ydt[p1];
    hi = ydt[p2];
227
    return (lo + (hi << 8) + (hi << 16)) << 1;
228 229 230 231 232
}

static int make_cdt24_entry(int p1, int p2, int16_t *cdt)
{
    int r, b;
233

234 235
    b = cdt[p2];
    r = cdt[p1]<<16;
236
    return (b+r) << 1;
237 238
}

Michael Niedermayer's avatar
Michael Niedermayer committed
239
static void gen_vector_table15(TrueMotion1Context *s, const uint8_t *sel_vector_table)
240 241 242
{
    int len, i, j;
    unsigned char delta_pair;
243

244 245 246 247 248 249
    for (i = 0; i < 1024; i += 4)
    {
        len = *sel_vector_table++ / 2;
        for (j = 0; j < len; j++)
        {
            delta_pair = *sel_vector_table++;
250
            s->y_predictor_table[i+j] = 0xfffffffe &
251
                make_ydt15_entry(delta_pair >> 4, delta_pair & 0xf, s->ydt);
252
            s->c_predictor_table[i+j] = 0xfffffffe &
253 254 255 256 257 258 259
                make_cdt15_entry(delta_pair >> 4, delta_pair & 0xf, s->cdt);
        }
        s->y_predictor_table[i+(j-1)] |= 1;
        s->c_predictor_table[i+(j-1)] |= 1;
    }
}

Michael Niedermayer's avatar
Michael Niedermayer committed
260
static void gen_vector_table16(TrueMotion1Context *s, const uint8_t *sel_vector_table)
261 262 263
{
    int len, i, j;
    unsigned char delta_pair;
264

265 266 267 268 269 270
    for (i = 0; i < 1024; i += 4)
    {
        len = *sel_vector_table++ / 2;
        for (j = 0; j < len; j++)
        {
            delta_pair = *sel_vector_table++;
271
            s->y_predictor_table[i+j] = 0xfffffffe &
272
                make_ydt16_entry(delta_pair >> 4, delta_pair & 0xf, s->ydt);
273
            s->c_predictor_table[i+j] = 0xfffffffe &
274 275 276 277 278 279 280
                make_cdt16_entry(delta_pair >> 4, delta_pair & 0xf, s->cdt);
        }
        s->y_predictor_table[i+(j-1)] |= 1;
        s->c_predictor_table[i+(j-1)] |= 1;
    }
}

Michael Niedermayer's avatar
Michael Niedermayer committed
281
static void gen_vector_table24(TrueMotion1Context *s, const uint8_t *sel_vector_table)
282 283 284
{
    int len, i, j;
    unsigned char delta_pair;
285

286 287 288 289 290 291
    for (i = 0; i < 1024; i += 4)
    {
        len = *sel_vector_table++ / 2;
        for (j = 0; j < len; j++)
        {
            delta_pair = *sel_vector_table++;
292
            s->y_predictor_table[i+j] = 0xfffffffe &
293
                make_ydt24_entry(delta_pair >> 4, delta_pair & 0xf, s->ydt);
294
            s->c_predictor_table[i+j] = 0xfffffffe &
295
                make_cdt24_entry(delta_pair >> 4, delta_pair & 0xf, s->cdt);
296
            s->fat_y_predictor_table[i+j] = 0xfffffffe &
297
                make_ydt24_entry(delta_pair >> 4, delta_pair & 0xf, s->fat_ydt);
298
            s->fat_c_predictor_table[i+j] = 0xfffffffe &
299
                make_cdt24_entry(delta_pair >> 4, delta_pair & 0xf, s->fat_cdt);
300 301 302
        }
        s->y_predictor_table[i+(j-1)] |= 1;
        s->c_predictor_table[i+(j-1)] |= 1;
303 304
        s->fat_y_predictor_table[i+(j-1)] |= 1;
        s->fat_c_predictor_table[i+(j-1)] |= 1;
305 306 307 308
    }
}

/* Returns the number of bytes consumed from the bytestream. Returns -1 if
309
 * there was an error while decoding the header */
310 311
static int truemotion1_decode_header(TrueMotion1Context *s)
{
312
    int i, ret;
313
    int width_shift = 0;
314
    int new_pix_fmt;
315
    struct frame_header header;
316
    uint8_t header_buffer[128] = { 0 };  /* logical maximum size of the header */
Michael Niedermayer's avatar
Michael Niedermayer committed
317
    const uint8_t *sel_vector_table;
318 319 320 321

    header.header_size = ((s->buf[0] >> 5) | (s->buf[0] << 3)) & 0x7f;
    if (s->buf[0] < 0x10)
    {
322
        av_log(s->avctx, AV_LOG_ERROR, "invalid header size (%d)\n", s->buf[0]);
323
        return AVERROR_INVALIDDATA;
324 325
    }

326 327 328
    if (header.header_size + 1 > s->size) {
        av_log(s->avctx, AV_LOG_ERROR, "Input packet too small.\n");
        return AVERROR_INVALIDDATA;
329 330 331 332
    }

    /* unscramble the header bytes with a XOR operation */
    for (i = 1; i < header.header_size; i++)
333
        header_buffer[i - 1] = s->buf[i] ^ s->buf[i + 1];
334

335 336 337
    header.compression = header_buffer[0];
    header.deltaset = header_buffer[1];
    header.vectable = header_buffer[2];
338 339 340
    header.ysize = AV_RL16(&header_buffer[3]);
    header.xsize = AV_RL16(&header_buffer[5]);
    header.checksum = AV_RL16(&header_buffer[7]);
341 342 343 344 345 346 347 348 349 350
    header.version = header_buffer[9];
    header.header_type = header_buffer[10];
    header.flags = header_buffer[11];
    header.control = header_buffer[12];

    /* Version 2 */
    if (header.version >= 2)
    {
        if (header.header_type > 3)
        {
351
            av_log(s->avctx, AV_LOG_ERROR, "invalid header type (%d)\n", header.header_type);
352
            return AVERROR_INVALIDDATA;
353 354 355 356 357 358 359 360
        } else if ((header.header_type == 2) || (header.header_type == 3)) {
            s->flags = header.flags;
            if (!(s->flags & FLAG_INTERFRAME))
                s->flags |= FLAG_KEYFRAME;
        } else
            s->flags = FLAG_KEYFRAME;
    } else /* Version 1 */
        s->flags = FLAG_KEYFRAME;
361

362
    if (s->flags & FLAG_SPRITE) {
363
        avpriv_request_sample(s->avctx, "Frame with sprite");
364
        /* FIXME header.width, height, xoffset and yoffset aren't initialized */
365
        return AVERROR_PATCHWELCOME;
366 367 368 369 370
    } else {
        s->w = header.xsize;
        s->h = header.ysize;
        if (header.header_type < 2) {
            if ((s->w < 213) && (s->h >= 176))
371
            {
372
                s->flags |= FLAG_INTERPOLATED;
373
                avpriv_request_sample(s->avctx, "Interpolated frame");
374
            }
375 376 377
        }
    }

378
    if (header.compression >= 17) {
Alex Beregszaszi's avatar
Alex Beregszaszi committed
379
        av_log(s->avctx, AV_LOG_ERROR, "invalid compression type (%d)\n", header.compression);
380
        return AVERROR_INVALIDDATA;
381
    }
382 383

    if ((header.deltaset != s->last_deltaset) ||
384 385 386 387 388 389
        (header.vectable != s->last_vectable))
        select_delta_tables(s, header.deltaset);

    if ((header.compression & 1) && header.header_type)
        sel_vector_table = pc_tbl2;
    else {
390
        if (header.vectable > 0 && header.vectable < 4)
391 392
            sel_vector_table = tables[header.vectable - 1];
        else {
Alex Beregszaszi's avatar
Alex Beregszaszi committed
393
            av_log(s->avctx, AV_LOG_ERROR, "invalid vector table id (%d)\n", header.vectable);
394
            return AVERROR_INVALIDDATA;
395 396
        }
    }
397

398
    if (compression_types[header.compression].algorithm == ALGO_RGB24H) {
399
        new_pix_fmt = AV_PIX_FMT_RGB32;
400 401
        width_shift = 1;
    } else
402
        new_pix_fmt = AV_PIX_FMT_RGB555; // RGB565 is supported as well
403

404 405
    s->w >>= width_shift;

406 407
    if (s->w != s->avctx->width || s->h != s->avctx->height ||
        new_pix_fmt != s->avctx->pix_fmt) {
408
        av_frame_unref(s->frame);
409
        s->avctx->sample_aspect_ratio = (AVRational){ 1 << width_shift, 1 };
410
        s->avctx->pix_fmt = new_pix_fmt;
411 412 413 414

        if ((ret = ff_set_dimensions(s->avctx, s->w, s->h)) < 0)
            return ret;

415 416
        ff_set_sar(s->avctx, s->avctx->sample_aspect_ratio);

417
        av_fast_malloc(&s->vert_pred, &s->vert_pred_size, s->avctx->width * sizeof(unsigned int));
418 419
        if (!s->vert_pred)
            return AVERROR(ENOMEM);
420
    }
421

422 423 424
    /* There is 1 change bit per 4 pixels, so each change byte represents
     * 32 pixels; divide width by 4 to obtain the number of change bits and
     * then round up to the nearest byte. */
425
    s->mb_change_bits_row_size = ((s->avctx->width >> (2 - width_shift)) + 7) >> 3;
426

427 428 429
    if ((header.deltaset != s->last_deltaset) || (header.vectable != s->last_vectable))
    {
        if (compression_types[header.compression].algorithm == ALGO_RGB24H)
430
            gen_vector_table24(s, sel_vector_table);
431
        else
432
        if (s->avctx->pix_fmt == AV_PIX_FMT_RGB555)
433
            gen_vector_table15(s, sel_vector_table);
434
        else
435
            gen_vector_table16(s, sel_vector_table);
436 437 438 439 440 441 442 443 444
    }

    /* set up pointers to the other key data chunks */
    s->mb_change_bits = s->buf + header.header_size;
    if (s->flags & FLAG_KEYFRAME) {
        /* no change bits specified for a keyframe; only index bytes */
        s->index_stream = s->mb_change_bits;
    } else {
        /* one change bit per 4x4 block */
445
        s->index_stream = s->mb_change_bits +
446 447 448 449 450 451 452 453 454 455 456
            (s->mb_change_bits_row_size * (s->avctx->height >> 2));
    }
    s->index_stream_size = s->size - (s->index_stream - s->buf);

    s->last_deltaset = header.deltaset;
    s->last_vectable = header.vectable;
    s->compression = header.compression;
    s->block_width = compression_types[header.compression].block_width;
    s->block_height = compression_types[header.compression].block_height;
    s->block_type = compression_types[header.compression].block_type;

457
    if (s->avctx->debug & FF_DEBUG_PICT_INFO)
458 459 460 461 462 463 464
        av_log(s->avctx, AV_LOG_INFO, "tables: %d / %d c:%d %dx%d t:%d %s%s%s%s\n",
            s->last_deltaset, s->last_vectable, s->compression, s->block_width,
            s->block_height, s->block_type,
            s->flags & FLAG_KEYFRAME ? " KEY" : "",
            s->flags & FLAG_INTERFRAME ? " INTER" : "",
            s->flags & FLAG_SPRITE ? " SPRITE" : "",
            s->flags & FLAG_INTERPOLATED ? " INTERPOL" : "");
465

466
    return header.header_size;
467 468
}

469
static av_cold int truemotion1_decode_init(AVCodecContext *avctx)
470
{
471
    TrueMotion1Context *s = avctx->priv_data;
472 473 474

    s->avctx = avctx;

475 476
    // FIXME: it may change ?
//    if (avctx->bits_per_sample == 24)
477
//        avctx->pix_fmt = AV_PIX_FMT_RGB24;
478
//    else
479
//        avctx->pix_fmt = AV_PIX_FMT_RGB555;
480

481 482 483
    s->frame = av_frame_alloc();
    if (!s->frame)
        return AVERROR(ENOMEM);
484 485 486

    /* there is a vertical predictor for each pixel in a line; each vertical
     * predictor is 0 to start with */
487
    av_fast_malloc(&s->vert_pred, &s->vert_pred_size, s->avctx->width * sizeof(unsigned int));
488 489
    if (!s->vert_pred)
        return AVERROR(ENOMEM);
490 491 492 493

    return 0;
}

494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519
/*
Block decoding order:

dxi: Y-Y
dxic: Y-C-Y
dxic2: Y-C-Y-C

hres,vres,i,i%vres (0 < i < 4)
2x2 0: 0 dxic2
2x2 1: 1 dxi
2x2 2: 0 dxic2
2x2 3: 1 dxi
2x4 0: 0 dxic2
2x4 1: 1 dxi
2x4 2: 2 dxi
2x4 3: 3 dxi
4x2 0: 0 dxic
4x2 1: 1 dxi
4x2 2: 0 dxic
4x2 3: 1 dxi
4x4 0: 0 dxic
4x4 1: 1 dxi
4x4 2: 2 dxi
4x4 3: 3 dxi
*/

520 521 522
#define GET_NEXT_INDEX() \
{\
    if (index_stream_index >= s->index_stream_size) { \
Alex Beregszaszi's avatar
Alex Beregszaszi committed
523
        av_log(s->avctx, AV_LOG_INFO, " help! truemotion1 decoder went out of bounds\n"); \
524 525 526 527 528
        return; \
    } \
    index = s->index_stream[index_stream_index++] * 4; \
}

529 530 531 532 533 534 535 536 537
#define INC_INDEX                                                   \
do {                                                                \
    if (index >= 1023) {                                            \
        av_log(s->avctx, AV_LOG_ERROR, "Invalid index value.\n");   \
        return;                                                     \
    }                                                               \
    index++;                                                        \
} while (0)

538 539 540 541 542 543 544 545 546 547 548 549
#define APPLY_C_PREDICTOR() \
    predictor_pair = s->c_predictor_table[index]; \
    horiz_pred += (predictor_pair >> 1); \
    if (predictor_pair & 1) { \
        GET_NEXT_INDEX() \
        if (!index) { \
            GET_NEXT_INDEX() \
            predictor_pair = s->c_predictor_table[index]; \
            horiz_pred += ((predictor_pair >> 1) * 5); \
            if (predictor_pair & 1) \
                GET_NEXT_INDEX() \
            else \
550
                INC_INDEX; \
551 552
        } \
    } else \
553
        INC_INDEX;
554

555 556
#define APPLY_C_PREDICTOR_24() \
    predictor_pair = s->c_predictor_table[index]; \
557
    horiz_pred += (predictor_pair >> 1); \
558 559 560 561 562
    if (predictor_pair & 1) { \
        GET_NEXT_INDEX() \
        if (!index) { \
            GET_NEXT_INDEX() \
            predictor_pair = s->fat_c_predictor_table[index]; \
563
            horiz_pred += (predictor_pair >> 1); \
564 565 566
            if (predictor_pair & 1) \
                GET_NEXT_INDEX() \
            else \
567
                INC_INDEX; \
568 569
        } \
    } else \
570
        INC_INDEX;
571 572


573 574 575 576 577 578 579 580 581 582 583 584
#define APPLY_Y_PREDICTOR() \
    predictor_pair = s->y_predictor_table[index]; \
    horiz_pred += (predictor_pair >> 1); \
    if (predictor_pair & 1) { \
        GET_NEXT_INDEX() \
        if (!index) { \
            GET_NEXT_INDEX() \
            predictor_pair = s->y_predictor_table[index]; \
            horiz_pred += ((predictor_pair >> 1) * 5); \
            if (predictor_pair & 1) \
                GET_NEXT_INDEX() \
            else \
585
                INC_INDEX; \
586 587
        } \
    } else \
588
        INC_INDEX;
589

590 591 592 593 594 595 596 597 598 599 600 601
#define APPLY_Y_PREDICTOR_24() \
    predictor_pair = s->y_predictor_table[index]; \
    horiz_pred += (predictor_pair >> 1); \
    if (predictor_pair & 1) { \
        GET_NEXT_INDEX() \
        if (!index) { \
            GET_NEXT_INDEX() \
            predictor_pair = s->fat_y_predictor_table[index]; \
            horiz_pred += (predictor_pair >> 1); \
            if (predictor_pair & 1) \
                GET_NEXT_INDEX() \
            else \
602
                INC_INDEX; \
603 604
        } \
    } else \
605
        INC_INDEX;
606

607 608
#define OUTPUT_PIXEL_PAIR() \
    *current_pixel_pair = *vert_pred + horiz_pred; \
609
    *vert_pred++ = *current_pixel_pair++;
610 611 612 613 614 615 616 617 618

static void truemotion1_decode_16bit(TrueMotion1Context *s)
{
    int y;
    int pixels_left;  /* remaining pixels on this line */
    unsigned int predictor_pair;
    unsigned int horiz_pred;
    unsigned int *vert_pred;
    unsigned int *current_pixel_pair;
619
    unsigned char *current_line = s->frame->data[0];
620 621 622
    int keyframe = s->flags & FLAG_KEYFRAME;

    /* these variables are for managing the stream of macroblock change bits */
Michael Niedermayer's avatar
Michael Niedermayer committed
623
    const unsigned char *mb_change_bits = s->mb_change_bits;
624 625 626 627 628 629 630 631 632
    unsigned char mb_change_byte;
    unsigned char mb_change_byte_mask;
    int mb_change_index;

    /* these variables are for managing the main index stream */
    int index_stream_index = 0;  /* yes, the index into the index stream */
    int index;

    /* clean out the line buffer */
633
    memset(s->vert_pred, 0, s->avctx->width * sizeof(unsigned int));
634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653

    GET_NEXT_INDEX();

    for (y = 0; y < s->avctx->height; y++) {

        /* re-init variables for the next line iteration */
        horiz_pred = 0;
        current_pixel_pair = (unsigned int *)current_line;
        vert_pred = s->vert_pred;
        mb_change_index = 0;
        mb_change_byte = mb_change_bits[mb_change_index++];
        mb_change_byte_mask = 0x01;
        pixels_left = s->avctx->width;

        while (pixels_left > 0) {

            if (keyframe || ((mb_change_byte & mb_change_byte_mask) == 0)) {

                switch (y & 3) {
                case 0:
654
                    /* if macroblock width is 2, apply C-Y-C-Y; else
655
                     * apply C-Y-Y */
656
                    if (s->block_width == 2) {
657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681
                        APPLY_C_PREDICTOR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_C_PREDICTOR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                    } else {
                        APPLY_C_PREDICTOR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                    }
                    break;

                case 1:
                case 3:
                    /* always apply 2 Y predictors on these iterations */
                    APPLY_Y_PREDICTOR();
                    OUTPUT_PIXEL_PAIR();
                    APPLY_Y_PREDICTOR();
                    OUTPUT_PIXEL_PAIR();
                    break;

                case 2:
682
                    /* this iteration might be C-Y-C-Y, Y-Y, or C-Y-Y
683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707
                     * depending on the macroblock type */
                    if (s->block_type == BLOCK_2x2) {
                        APPLY_C_PREDICTOR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_C_PREDICTOR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                    } else if (s->block_type == BLOCK_4x2) {
                        APPLY_C_PREDICTOR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                    } else {
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_Y_PREDICTOR();
                        OUTPUT_PIXEL_PAIR();
                    }
                    break;
                }

            } else {

708
                /* skip (copy) four pixels, but reassign the horizontal
709 710 711 712
                 * predictor */
                *vert_pred++ = *current_pixel_pair++;
                horiz_pred = *current_pixel_pair - *vert_pred;
                *vert_pred++ = *current_pixel_pair++;
713

714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732
            }

            if (!keyframe) {
                mb_change_byte_mask <<= 1;

                /* next byte */
                if (!mb_change_byte_mask) {
                    mb_change_byte = mb_change_bits[mb_change_index++];
                    mb_change_byte_mask = 0x01;
                }
            }

            pixels_left -= 4;
        }

        /* next change row */
        if (((y + 1) & 3) == 0)
            mb_change_bits += s->mb_change_bits_row_size;

733
        current_line += s->frame->linesize[0];
734 735 736
    }
}

737 738 739 740 741 742 743 744
static void truemotion1_decode_24bit(TrueMotion1Context *s)
{
    int y;
    int pixels_left;  /* remaining pixels on this line */
    unsigned int predictor_pair;
    unsigned int horiz_pred;
    unsigned int *vert_pred;
    unsigned int *current_pixel_pair;
745
    unsigned char *current_line = s->frame->data[0];
746 747 748
    int keyframe = s->flags & FLAG_KEYFRAME;

    /* these variables are for managing the stream of macroblock change bits */
Michael Niedermayer's avatar
Michael Niedermayer committed
749
    const unsigned char *mb_change_bits = s->mb_change_bits;
750 751 752 753 754 755 756 757 758
    unsigned char mb_change_byte;
    unsigned char mb_change_byte_mask;
    int mb_change_index;

    /* these variables are for managing the main index stream */
    int index_stream_index = 0;  /* yes, the index into the index stream */
    int index;

    /* clean out the line buffer */
759
    memset(s->vert_pred, 0, s->avctx->width * sizeof(unsigned int));
760 761 762 763 764 765

    GET_NEXT_INDEX();

    for (y = 0; y < s->avctx->height; y++) {

        /* re-init variables for the next line iteration */
766
        horiz_pred = 0;
767 768 769 770 771 772 773 774 775 776 777 778 779
        current_pixel_pair = (unsigned int *)current_line;
        vert_pred = s->vert_pred;
        mb_change_index = 0;
        mb_change_byte = mb_change_bits[mb_change_index++];
        mb_change_byte_mask = 0x01;
        pixels_left = s->avctx->width;

        while (pixels_left > 0) {

            if (keyframe || ((mb_change_byte & mb_change_byte_mask) == 0)) {

                switch (y & 3) {
                case 0:
780
                    /* if macroblock width is 2, apply C-Y-C-Y; else
781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807
                     * apply C-Y-Y */
                    if (s->block_width == 2) {
                        APPLY_C_PREDICTOR_24();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_C_PREDICTOR_24();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                    } else {
                        APPLY_C_PREDICTOR_24();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                    }
                    break;

                case 1:
                case 3:
                    /* always apply 2 Y predictors on these iterations */
                    APPLY_Y_PREDICTOR_24();
                    OUTPUT_PIXEL_PAIR();
                    APPLY_Y_PREDICTOR_24();
                    OUTPUT_PIXEL_PAIR();
                    break;

                case 2:
808
                    /* this iteration might be C-Y-C-Y, Y-Y, or C-Y-Y
809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833
                     * depending on the macroblock type */
                    if (s->block_type == BLOCK_2x2) {
                        APPLY_C_PREDICTOR_24();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_C_PREDICTOR_24();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                    } else if (s->block_type == BLOCK_4x2) {
                        APPLY_C_PREDICTOR_24();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                    } else {
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                        APPLY_Y_PREDICTOR_24();
                        OUTPUT_PIXEL_PAIR();
                    }
                    break;
                }

            } else {

834
                /* skip (copy) four pixels, but reassign the horizontal
835 836 837 838
                 * predictor */
                *vert_pred++ = *current_pixel_pair++;
                horiz_pred = *current_pixel_pair - *vert_pred;
                *vert_pred++ = *current_pixel_pair++;
839

840 841 842 843 844 845 846 847 848 849 850 851
            }

            if (!keyframe) {
                mb_change_byte_mask <<= 1;

                /* next byte */
                if (!mb_change_byte_mask) {
                    mb_change_byte = mb_change_bits[mb_change_index++];
                    mb_change_byte_mask = 0x01;
                }
            }

852
            pixels_left -= 2;
853 854 855 856 857 858
        }

        /* next change row */
        if (((y + 1) & 3) == 0)
            mb_change_bits += s->mb_change_bits_row_size;

859
        current_line += s->frame->linesize[0];
860 861 862 863
    }
}


864
static int truemotion1_decode_frame(AVCodecContext *avctx,
865
                                    void *data, int *got_frame,
866
                                    AVPacket *avpkt)
867
{
868
    const uint8_t *buf = avpkt->data;
869
    int ret, buf_size = avpkt->size;
870
    TrueMotion1Context *s = avctx->priv_data;
871 872 873 874

    s->buf = buf;
    s->size = buf_size;

875 876
    if ((ret = truemotion1_decode_header(s)) < 0)
        return ret;
877

878
    if ((ret = ff_reget_buffer(avctx, s->frame)) < 0)
879
        return ret;
880

881
    if (compression_types[s->compression].algorithm == ALGO_RGB24H) {
882
        truemotion1_decode_24bit(s);
883
    } else if (compression_types[s->compression].algorithm != ALGO_NOP) {
884 885 886
        truemotion1_decode_16bit(s);
    }

887
    if ((ret = av_frame_ref(data, s->frame)) < 0)
888 889
        return ret;

890
    *got_frame      = 1;
891 892 893 894 895

    /* report that the buffer was completely consumed */
    return buf_size;
}

896
static av_cold int truemotion1_decode_end(AVCodecContext *avctx)
897
{
898
    TrueMotion1Context *s = avctx->priv_data;
899

900
    av_frame_free(&s->frame);
901
    av_freep(&s->vert_pred);
902 903 904 905

    return 0;
}

906
AVCodec ff_truemotion1_decoder = {
907
    .name           = "truemotion1",
908
    .long_name      = NULL_IF_CONFIG_SMALL("Duck TrueMotion 1.0"),
909
    .type           = AVMEDIA_TYPE_VIDEO,
910
    .id             = AV_CODEC_ID_TRUEMOTION1,
911 912 913 914 915
    .priv_data_size = sizeof(TrueMotion1Context),
    .init           = truemotion1_decode_init,
    .close          = truemotion1_decode_end,
    .decode         = truemotion1_decode_frame,
    .capabilities   = CODEC_CAP_DR1,
916
};