rtjpeg.c 5.75 KB
Newer Older
1 2 3 4
/*
 * RTJpeg decoding functions
 * Copyright (c) 2006 Reimar Doeffinger
 *
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 20
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */
21
#include "libavutil/common.h"
22
#include "get_bits.h"
23 24 25 26 27 28
#include "rtjpeg.h"

#define PUT_COEFF(c) \
    i = scan[coeff--]; \
    block[i] = (c) * quant[i];

29
/// aligns the bitstream to the given power of two
30 31 32 33 34
#define ALIGN(a) \
    n = (-get_bits_count(gb)) & (a - 1); \
    if (n) {skip_bits(gb, n);}

/**
35 36 37 38 39 40
 * @brief read one block from stream
 * @param gb contains stream data
 * @param block where data is written to
 * @param scan array containing the mapping stream address -> block position
 * @param quant quantization factors
 * @return 0 means the block is not coded, < 0 means an error occurred.
41 42 43
 *
 * Note: GetBitContext is used to make the code simpler, since all data is
 * aligned this could be done faster in a different way, e.g. as it is done
44
 * in MPlayer libmpcodecs/native/rtjpegn.c.
45
 */
Diego Biurrun's avatar
Diego Biurrun committed
46
static inline int get_block(GetBitContext *gb, int16_t *block, const uint8_t *scan,
47
                            const uint32_t *quant) {
48 49 50 51 52 53 54 55 56 57
    int coeff, i, n;
    int8_t ac;
    uint8_t dc = get_bits(gb, 8);

    // block not coded
    if (dc == 255)
       return 0;

    // number of non-zero coefficients
    coeff = get_bits(gb, 6);
58
    if (get_bits_left(gb) < (coeff << 1))
59
        return AVERROR_INVALIDDATA;
60

61 62
    // normally we would only need to clear the (63 - coeff) last values,
    // but since we do not know where they are we just clear the whole block
Diego Biurrun's avatar
Diego Biurrun committed
63
    memset(block, 0, 64 * sizeof(int16_t));
64 65 66 67 68 69 70 71 72 73 74

    // 2 bits per coefficient
    while (coeff) {
        ac = get_sbits(gb, 2);
        if (ac == -2)
            break; // continue with more bits
        PUT_COEFF(ac);
    }

    // 4 bits per coefficient
    ALIGN(4);
75
    if (get_bits_left(gb) < (coeff << 2))
76
        return AVERROR_INVALIDDATA;
77 78 79 80 81 82 83 84 85
    while (coeff) {
        ac = get_sbits(gb, 4);
        if (ac == -8)
            break; // continue with more bits
        PUT_COEFF(ac);
    }

    // 8 bits per coefficient
    ALIGN(8);
86
    if (get_bits_left(gb) < (coeff << 3))
87
        return AVERROR_INVALIDDATA;
88 89 90 91 92 93 94 95 96 97
    while (coeff) {
        ac = get_sbits(gb, 8);
        PUT_COEFF(ac);
    }

    PUT_COEFF(dc);
    return 1;
}

/**
98
 * @brief decode one rtjpeg YUV420 frame
99
 * @param c context, must be initialized via ff_rtjpeg_decode_init
100
 * @param f AVFrame to place decoded frame into. If parts of the frame
101
 *          are not coded they are left unchanged, so consider initializing it
102 103 104
 * @param buf buffer containing input data
 * @param buf_size length of input data in bytes
 * @return number of bytes consumed from the input buffer
105
 */
106 107
int ff_rtjpeg_decode_frame_yuv420(RTJpegContext *c, AVFrame *f,
                                  const uint8_t *buf, int buf_size) {
108 109
    GetBitContext gb;
    int w = c->w / 16, h = c->h / 16;
110
    int x, y, ret;
111 112
    uint8_t *y1 = f->data[0], *y2 = f->data[0] + 8 * f->linesize[0];
    uint8_t *u = f->data[1], *v = f->data[2];
113 114 115 116

    if ((ret = init_get_bits8(&gb, buf, buf_size)) < 0)
        return ret;

117 118
    for (y = 0; y < h; y++) {
        for (x = 0; x < w; x++) {
119 120 121 122 123
#define BLOCK(quant, dst, stride) do { \
    int res = get_block(&gb, block, c->scan, quant); \
    if (res < 0) \
        return res; \
    if (res > 0) \
124
        c->dsp.idct_put(dst, stride, block); \
125
} while (0)
Diego Biurrun's avatar
Diego Biurrun committed
126
            int16_t *block = c->block;
127
            BLOCK(c->lquant, y1, f->linesize[0]);
128
            y1 += 8;
129
            BLOCK(c->lquant, y1, f->linesize[0]);
130
            y1 += 8;
131
            BLOCK(c->lquant, y2, f->linesize[0]);
132
            y2 += 8;
133
            BLOCK(c->lquant, y2, f->linesize[0]);
134
            y2 += 8;
135
            BLOCK(c->cquant, u,  f->linesize[1]);
136
            u += 8;
137
            BLOCK(c->cquant, v,  f->linesize[2]);
138 139 140 141 142 143 144 145 146 147 148
            v += 8;
        }
        y1 += 2 * 8 * (f->linesize[0] - w);
        y2 += 2 * 8 * (f->linesize[0] - w);
        u += 8 * (f->linesize[1] - w);
        v += 8 * (f->linesize[2] - w);
    }
    return get_bits_count(&gb) / 8;
}

/**
149 150 151
 * @brief initialize an RTJpegContext, may be called multiple times
 * @param c context to initialize
 * @param width width of image, will be rounded down to the nearest multiple
152
 *              of 16 for decoding
153
 * @param height height of image, will be rounded down to the nearest multiple
154
 *              of 16 for decoding
155 156
 * @param lquant luma quantization table to use
 * @param cquant chroma quantization table to use
157
 */
158
void ff_rtjpeg_decode_init(RTJpegContext *c, int width, int height,
159
                           const uint32_t *lquant, const uint32_t *cquant) {
160 161
    int i;
    for (i = 0; i < 64; i++) {
162
        int p = c->dsp.idct_permutation[i];
163 164 165 166 167 168
        c->lquant[p] = lquant[i];
        c->cquant[p] = cquant[i];
    }
    c->w = width;
    c->h = height;
}
169 170 171 172 173 174 175 176 177 178 179 180 181 182 183

void ff_rtjpeg_init(RTJpegContext *c, AVCodecContext *avctx)
{
    int i;

    ff_dsputil_init(&c->dsp, avctx);

    for (i = 0; i < 64; i++) {
        int z = ff_zigzag_direct[i];
        z = ((z << 3) | (z >> 3)) & 63; // rtjpeg uses a transposed variant

        // permute the scan and quantization tables for the chosen idct
        c->scan[i] = c->dsp.idct_permutation[z];
    }
}