Statistics
| Branch: | Revision:

ffmpeg / libavcodec / vp5.c @ 9106a698

History | View | Annotate | Download (10.2 KB)

1
/**
2
 * @file libavcodec/vp5.c
3
 * VP5 compatible video decoder
4
 *
5
 * Copyright (C) 2006  Aurelien Jacobs <aurel@gnuage.org>
6
 *
7
 * This file is part of FFmpeg.
8
 *
9
 * FFmpeg is free software; you can redistribute it and/or
10
 * modify it under the terms of the GNU Lesser General Public
11
 * License as published by the Free Software Foundation; either
12
 * version 2.1 of the License, or (at your option) any later version.
13
 *
14
 * FFmpeg is distributed in the hope that it will be useful,
15
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
17
 * Lesser General Public License for more details.
18
 *
19
 * You should have received a copy of the GNU Lesser General Public
20
 * License along with FFmpeg; if not, write to the Free Software
21
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22
 */
23

    
24
#include <stdlib.h>
25
#include <string.h>
26

    
27
#include "avcodec.h"
28
#include "dsputil.h"
29
#include "get_bits.h"
30

    
31
#include "vp56.h"
32
#include "vp56data.h"
33
#include "vp5data.h"
34

    
35

    
36
static int vp5_parse_header(VP56Context *s, const uint8_t *buf, int buf_size,
37
                            int *golden_frame)
38
{
39
    VP56RangeCoder *c = &s->c;
40
    int rows, cols;
41

    
42
    vp56_init_range_decoder(&s->c, buf, buf_size);
43
    s->framep[VP56_FRAME_CURRENT]->key_frame = !vp56_rac_get(c);
44
    vp56_rac_get(c);
45
    vp56_init_dequant(s, vp56_rac_gets(c, 6));
46
    if (s->framep[VP56_FRAME_CURRENT]->key_frame)
47
    {
48
        vp56_rac_gets(c, 8);
49
        if(vp56_rac_gets(c, 5) > 5)
50
            return 0;
51
        vp56_rac_gets(c, 2);
52
        if (vp56_rac_get(c)) {
53
            av_log(s->avctx, AV_LOG_ERROR, "interlacing not supported\n");
54
            return 0;
55
        }
56
        rows = vp56_rac_gets(c, 8);  /* number of stored macroblock rows */
57
        cols = vp56_rac_gets(c, 8);  /* number of stored macroblock cols */
58
        vp56_rac_gets(c, 8);  /* number of displayed macroblock rows */
59
        vp56_rac_gets(c, 8);  /* number of displayed macroblock cols */
60
        vp56_rac_gets(c, 2);
61
        if (!s->macroblocks || /* first frame */
62
            16*cols != s->avctx->coded_width ||
63
            16*rows != s->avctx->coded_height) {
64
            avcodec_set_dimensions(s->avctx, 16*cols, 16*rows);
65
            return 2;
66
        }
67
    } else if (!s->macroblocks)
68
        return 0;
69
    return 1;
70
}
71

    
72
/* Gives very similar result than the vp6 version except in a few cases */
73
static int vp5_adjust(int v, int t)
74
{
75
    int s2, s1 = v >> 31;
76
    v ^= s1;
77
    v -= s1;
78
    v *= v < 2*t;
79
    v -= t;
80
    s2 = v >> 31;
81
    v ^= s2;
82
    v -= s2;
83
    v = t - v;
84
    v += s1;
85
    v ^= s1;
86
    return v;
87
}
88

    
89
static void vp5_parse_vector_adjustment(VP56Context *s, VP56mv *vect)
90
{
91
    VP56RangeCoder *c = &s->c;
92
    VP56Model *model = s->modelp;
93
    int comp, di;
94

    
95
    for (comp=0; comp<2; comp++) {
96
        int delta = 0;
97
        if (vp56_rac_get_prob(c, model->vector_dct[comp])) {
98
            int sign = vp56_rac_get_prob(c, model->vector_sig[comp]);
99
            di  = vp56_rac_get_prob(c, model->vector_pdi[comp][0]);
100
            di |= vp56_rac_get_prob(c, model->vector_pdi[comp][1]) << 1;
101
            delta = vp56_rac_get_tree(c, vp56_pva_tree,
102
                                      model->vector_pdv[comp]);
103
            delta = di | (delta << 2);
104
            delta = (delta ^ -sign) + sign;
105
        }
106
        if (!comp)
107
            vect->x = delta;
108
        else
109
            vect->y = delta;
110
    }
111
}
112

    
113
static void vp5_parse_vector_models(VP56Context *s)
114
{
115
    VP56RangeCoder *c = &s->c;
116
    VP56Model *model = s->modelp;
117
    int comp, node;
118

    
119
    for (comp=0; comp<2; comp++) {
120
        if (vp56_rac_get_prob(c, vp5_vmc_pct[comp][0]))
121
            model->vector_dct[comp] = vp56_rac_gets_nn(c, 7);
122
        if (vp56_rac_get_prob(c, vp5_vmc_pct[comp][1]))
123
            model->vector_sig[comp] = vp56_rac_gets_nn(c, 7);
124
        if (vp56_rac_get_prob(c, vp5_vmc_pct[comp][2]))
125
            model->vector_pdi[comp][0] = vp56_rac_gets_nn(c, 7);
126
        if (vp56_rac_get_prob(c, vp5_vmc_pct[comp][3]))
127
            model->vector_pdi[comp][1] = vp56_rac_gets_nn(c, 7);
128
    }
129

    
130
    for (comp=0; comp<2; comp++)
131
        for (node=0; node<7; node++)
132
            if (vp56_rac_get_prob(c, vp5_vmc_pct[comp][4 + node]))
133
                model->vector_pdv[comp][node] = vp56_rac_gets_nn(c, 7);
134
}
135

    
136
static void vp5_parse_coeff_models(VP56Context *s)
137
{
138
    VP56RangeCoder *c = &s->c;
139
    VP56Model *model = s->modelp;
140
    uint8_t def_prob[11];
141
    int node, cg, ctx;
142
    int ct;    /* code type */
143
    int pt;    /* plane type (0 for Y, 1 for U or V) */
144

    
145
    memset(def_prob, 0x80, sizeof(def_prob));
146

    
147
    for (pt=0; pt<2; pt++)
148
        for (node=0; node<11; node++)
149
            if (vp56_rac_get_prob(c, vp5_dccv_pct[pt][node])) {
150
                def_prob[node] = vp56_rac_gets_nn(c, 7);
151
                model->coeff_dccv[pt][node] = def_prob[node];
152
            } else if (s->framep[VP56_FRAME_CURRENT]->key_frame) {
153
                model->coeff_dccv[pt][node] = def_prob[node];
154
            }
155

    
156
    for (ct=0; ct<3; ct++)
157
        for (pt=0; pt<2; pt++)
158
            for (cg=0; cg<6; cg++)
159
                for (node=0; node<11; node++)
160
                    if (vp56_rac_get_prob(c, vp5_ract_pct[ct][pt][cg][node])) {
161
                        def_prob[node] = vp56_rac_gets_nn(c, 7);
162
                        model->coeff_ract[pt][ct][cg][node] = def_prob[node];
163
                    } else if (s->framep[VP56_FRAME_CURRENT]->key_frame) {
164
                        model->coeff_ract[pt][ct][cg][node] = def_prob[node];
165
                    }
166

    
167
    /* coeff_dcct is a linear combination of coeff_dccv */
168
    for (pt=0; pt<2; pt++)
169
        for (ctx=0; ctx<36; ctx++)
170
            for (node=0; node<5; node++)
171
                model->coeff_dcct[pt][ctx][node] = av_clip(((model->coeff_dccv[pt][node] * vp5_dccv_lc[node][ctx][0] + 128) >> 8) + vp5_dccv_lc[node][ctx][1], 1, 254);
172

    
173
    /* coeff_acct is a linear combination of coeff_ract */
174
    for (ct=0; ct<3; ct++)
175
        for (pt=0; pt<2; pt++)
176
            for (cg=0; cg<3; cg++)
177
                for (ctx=0; ctx<6; ctx++)
178
                    for (node=0; node<5; node++)
179
                        model->coeff_acct[pt][ct][cg][ctx][node] = av_clip(((model->coeff_ract[pt][ct][cg][node] * vp5_ract_lc[ct][cg][node][ctx][0] + 128) >> 8) + vp5_ract_lc[ct][cg][node][ctx][1], 1, 254);
180
}
181

    
182
static void vp5_parse_coeff(VP56Context *s)
183
{
184
    VP56RangeCoder *c = &s->c;
185
    VP56Model *model = s->modelp;
186
    uint8_t *permute = s->scantable.permutated;
187
    uint8_t *model1, *model2;
188
    int coeff, sign, coeff_idx;
189
    int b, i, cg, idx, ctx, ctx_last;
190
    int pt = 0;    /* plane type (0 for Y, 1 for U or V) */
191

    
192
    for (b=0; b<6; b++) {
193
        int ct = 1;    /* code type */
194

    
195
        if (b > 3) pt = 1;
196

    
197
        ctx = 6*s->coeff_ctx[vp56_b6to4[b]][0]
198
              + s->above_blocks[s->above_block_idx[b]].not_null_dc;
199
        model1 = model->coeff_dccv[pt];
200
        model2 = model->coeff_dcct[pt][ctx];
201

    
202
        for (coeff_idx=0; coeff_idx<64; ) {
203
            if (vp56_rac_get_prob(c, model2[0])) {
204
                if (vp56_rac_get_prob(c, model2[2])) {
205
                    if (vp56_rac_get_prob(c, model2[3])) {
206
                        s->coeff_ctx[vp56_b6to4[b]][coeff_idx] = 4;
207
                        idx = vp56_rac_get_tree(c, vp56_pc_tree, model1);
208
                        sign = vp56_rac_get(c);
209
                        coeff = vp56_coeff_bias[idx+5];
210
                        for (i=vp56_coeff_bit_length[idx]; i>=0; i--)
211
                            coeff += vp56_rac_get_prob(c, vp56_coeff_parse_table[idx][i]) << i;
212
                    } else {
213
                        if (vp56_rac_get_prob(c, model2[4])) {
214
                            coeff = 3 + vp56_rac_get_prob(c, model1[5]);
215
                            s->coeff_ctx[vp56_b6to4[b]][coeff_idx] = 3;
216
                        } else {
217
                            coeff = 2;
218
                            s->coeff_ctx[vp56_b6to4[b]][coeff_idx] = 2;
219
                        }
220
                        sign = vp56_rac_get(c);
221
                    }
222
                    ct = 2;
223
                } else {
224
                    ct = 1;
225
                    s->coeff_ctx[vp56_b6to4[b]][coeff_idx] = 1;
226
                    sign = vp56_rac_get(c);
227
                    coeff = 1;
228
                }
229
                coeff = (coeff ^ -sign) + sign;
230
                if (coeff_idx)
231
                    coeff *= s->dequant_ac;
232
                s->block_coeff[b][permute[coeff_idx]] = coeff;
233
            } else {
234
                if (ct && !vp56_rac_get_prob(c, model2[1]))
235
                    break;
236
                ct = 0;
237
                s->coeff_ctx[vp56_b6to4[b]][coeff_idx] = 0;
238
            }
239

    
240
            cg = vp5_coeff_groups[++coeff_idx];
241
            ctx = s->coeff_ctx[vp56_b6to4[b]][coeff_idx];
242
            model1 = model->coeff_ract[pt][ct][cg];
243
            model2 = cg > 2 ? model1 : model->coeff_acct[pt][ct][cg][ctx];
244
        }
245

    
246
        ctx_last = FFMIN(s->coeff_ctx_last[vp56_b6to4[b]], 24);
247
        s->coeff_ctx_last[vp56_b6to4[b]] = coeff_idx;
248
        if (coeff_idx < ctx_last)
249
            for (i=coeff_idx; i<=ctx_last; i++)
250
                s->coeff_ctx[vp56_b6to4[b]][i] = 5;
251
        s->above_blocks[s->above_block_idx[b]].not_null_dc = s->coeff_ctx[vp56_b6to4[b]][0];
252
    }
253
}
254

    
255
static void vp5_default_models_init(VP56Context *s)
256
{
257
    VP56Model *model = s->modelp;
258
    int i;
259

    
260
    for (i=0; i<2; i++) {
261
        model->vector_sig[i] = 0x80;
262
        model->vector_dct[i] = 0x80;
263
        model->vector_pdi[i][0] = 0x55;
264
        model->vector_pdi[i][1] = 0x80;
265
    }
266
    memcpy(model->mb_types_stats, vp56_def_mb_types_stats, sizeof(model->mb_types_stats));
267
    memset(model->vector_pdv, 0x80, sizeof(model->vector_pdv));
268
}
269

    
270
static av_cold int vp5_decode_init(AVCodecContext *avctx)
271
{
272
    VP56Context *s = avctx->priv_data;
273

    
274
    vp56_init(avctx, 1, 0);
275
    s->vp56_coord_div = vp5_coord_div;
276
    s->parse_vector_adjustment = vp5_parse_vector_adjustment;
277
    s->adjust = vp5_adjust;
278
    s->parse_coeff = vp5_parse_coeff;
279
    s->default_models_init = vp5_default_models_init;
280
    s->parse_vector_models = vp5_parse_vector_models;
281
    s->parse_coeff_models = vp5_parse_coeff_models;
282
    s->parse_header = vp5_parse_header;
283

    
284
    return 0;
285
}
286

    
287
AVCodec vp5_decoder = {
288
    "vp5",
289
    CODEC_TYPE_VIDEO,
290
    CODEC_ID_VP5,
291
    sizeof(VP56Context),
292
    vp5_decode_init,
293
    NULL,
294
    vp56_free,
295
    vp56_decode_frame,
296
    CODEC_CAP_DR1,
297
    .long_name = NULL_IF_CONFIG_SMALL("On2 VP5"),
298
};