1 /*
2 * Copyright (C) 2006 Aurelien Jacobs <aurel@gnuage.org>
3 *
4 * This file is part of FFmpeg.
5 *
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU Lesser General Public
8 * License as published by the Free Software Foundation; either
9 * version 2.1 of the License, or (at your option) any later version.
10 *
11 * FFmpeg is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 * Lesser General Public License for more details.
15 *
16 * You should have received a copy of the GNU Lesser General Public
17 * License along with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
21 /**
22 * @file
23 * VP6 compatible video decoder
24 *
25 * The VP6F decoder accepts an optional 1 byte extradata. It is composed of:
26 * - upper 4 bits: difference between encoded width and visible width
27 * - lower 4 bits: difference between encoded height and visible height
28 */
29
30 #include <stdlib.h>
31
36
40
41 #define VP6_MAX_HUFF_SIZE 12
42
45
47 {
49 int parse_filter_info = 0;
50 int coeff_offset = 0;
51 int vrt_shift = 0;
52 int sub_version;
53 int rows, cols;
54 int res = 0;
55 int separated_coeff = buf[0] & 1;
56
59
61 sub_version = buf[1] >> 3;
62 if (sub_version > 8)
64 s->filter_header = buf[1] & 0x06;
65 if (buf[1] & 1) {
68 }
69 if (separated_coeff || !s->filter_header) {
70 coeff_offset =
AV_RB16(buf+2) - 2;
71 buf += 2;
72 buf_size -= 2;
73 }
74
75 rows = buf[2]; /* number of stored macroblock rows */
76 cols = buf[3]; /* number of stored macroblock cols */
77 /* buf[4] is number of displayed macroblock rows */
78 /* buf[5] is number of displayed macroblock cols */
79 if (!rows || !cols) {
82 }
83
84 if (!s->macroblocks || /* first frame */
85 16*cols != s->avctx->coded_width ||
86 16*rows != s->avctx->coded_height) {
87 if (s->avctx->extradata_size == 0 &&
88 FFALIGN(s->avctx->width, 16) == 16 * cols &&
89 FFALIGN(s->avctx->height, 16) == 16 * rows) {
90 // We assume this is properly signalled container cropping,
91 // in an F4V file. Just set the coded_width/height, don't
92 // touch the cropped ones.
93 s->avctx->coded_width = 16 * cols;
94 s->avctx->coded_height = 16 * rows;
95 } else {
97 if (ret < 0)
99
100 if (s->avctx->extradata_size == 1) {
101 s->avctx->width -= s->avctx->extradata[0] >> 4;
102 s->avctx->height -= s->avctx->extradata[0] & 0x0F;
103 }
104 }
106 }
107
110
111 parse_filter_info = s->filter_header;
112 if (sub_version < 8)
113 vrt_shift = 5;
114 s->sub_version = sub_version;
115 s->golden_frame = 0;
116 } else {
117 if (!s->sub_version || !s->avctx->coded_width || !s->avctx->coded_height)
119
120 if (separated_coeff || !s->filter_header) {
121 coeff_offset =
AV_RB16(buf+1) - 2;
122 buf += 2;
123 buf_size -= 2;
124 }
126
128 if (s->filter_header) {
130 if (s->deblock_filtering)
132 if (s->sub_version > 7)
134 }
135 }
136
137 if (parse_filter_info) {
139 s->filter_mode = 2;
140 s->sample_variance_threshold =
vp56_rac_gets(c, 5) << vrt_shift;
143 s->filter_mode = 1;
144 } else {
145 s->filter_mode = 0;
146 }
147 if (s->sub_version > 7)
149 else
150 s->filter_selection = 16;
151 }
152
154
156 if (coeff_offset) {
157 buf += coeff_offset;
158 buf_size -= coeff_offset;
159 if (buf_size < 0) {
163 }
164 if (s->use_huffman) {
167 } else {
169 s->ccp = &s->cc;
170 }
171 } else {
172 s->ccp = &s->c;
173 }
174
175 return res;
176 }
177
179 {
180 int i, pos, idx = 1;
181
182 s->modelp->coeff_index_to_pos[0] = 0;
183 for (i=0; i<16; i++)
184 for (pos=1; pos<64; pos++)
185 if (s->modelp->coeff_reorder[pos] == i)
186 s->modelp->coeff_index_to_pos[idx++] = pos;
187 }
188
190 {
192
197
203
205 }
206
208 {
212
213 for (comp=0; comp<2; comp++) {
218 }
219
220 for (comp=0; comp<2; comp++)
221 for (node=0; node<7; node++)
224
225 for (comp=0; comp<2; comp++)
226 for (node=0; node<8; node++)
229 }
230
231 /* nodes must ascend by count, but with descending symbol order */
233 {
234 const Node *
a = va, *
b = vb;
235 return (a->
count - b->count)*16 + (b->sym - a->
sym);
236 }
237
240 {
243
244 /* first compute probabilities from model */
246 for (i=0; i<size-1; i++) {
247 a = tmp[i].
count * coeff_model[i] >> 8;
248 b = tmp[i].
count * (255 - coeff_model[i]) >> 8;
249 nodes[map[2*i ]].
count = a + !
a;
250 nodes[map[2*i+1]].
count = b + !
b;
251 }
252
254 /* then build the huffman tree according to probabilities */
258 }
259
261 {
264 int def_prob[11];
265 int node, cg, ctx, pos;
266 int ct; /* code type */
267 int pt;
/* plane type (0 for Y, 1 for U or V) */
268
269 memset(def_prob, 0x80, sizeof(def_prob));
270
271 for (pt=0; pt<2; pt++)
272 for (node=0; node<11; node++)
278 }
279
281 for (pos=1; pos<64; pos++)
285 }
286
287 for (cg=0; cg<2; cg++)
288 for (node=0; node<14; node++)
291
292 for (ct=0; ct<3; ct++)
293 for (pt=0; pt<2; pt++)
294 for (cg=0; cg<6; cg++)
295 for (node=0; node<11; node++)
301 }
302
303 if (s->use_huffman) {
304 for (pt=0; pt<2; pt++) {
307 return -1;
310 return -1;
311 for (ct=0; ct<3; ct++)
312 for (cg = 0; cg < 6; cg++)
315 &s->ract_vlc[pt][ct][cg]))
316 return -1;
317 }
318 memset(s->nb_null, 0, sizeof(s->nb_null));
319 } else {
320 /* coeff_dcct is a linear combination of coeff_dccv */
321 for (pt=0; pt<2; pt++)
322 for (ctx=0; ctx<3; ctx++)
323 for (node=0; node<5; node++)
325 }
326 return 0;
327 }
328
330 {
334
336 if (s->vector_candidate_pos < 2)
337 *vect = s->vector_candidate[0];
338
339 for (comp=0; comp<2; comp++) {
341
343 static const uint8_t prob_order[] = {0, 1, 2, 7, 6, 5, 4};
344 for (i=0; i<sizeof(prob_order); i++) {
345 int j = prob_order[i];
347 }
348 if (delta & 0xF0)
350 else
351 delta |= 8;
352 } else {
355 }
356
359
360 if (!comp)
362 else
364 }
365 }
366
367 /**
368 * Read number of consecutive blocks with null DC or AC.
369 * This value is < 74.
370 */
372 {
374 if (val == 2)
376 else if (val == 3) {
378 val = 6+val +
get_bits(&s->gb, 2+val);
379 }
381 }
382
384 {
388 int coeff, sign, coeff_idx;
390 int pt = 0;
/* plane type (0 for Y, 1 for U or V) */
391
392 for (b=0; b<6; b++) {
393 int ct = 0; /* code type */
394 if (b > 3) pt = 1;
395 vlc_coeff = &s->dccv_vlc[
pt];
396
397 for (coeff_idx = 0;;) {
399 if (coeff_idx<2 && s->nb_null[coeff_idx][pt]) {
400 s->nb_null[coeff_idx][
pt]--;
401 if (coeff_idx)
402 break;
403 } else {
405 return;
407 if (coeff == 0) {
408 if (coeff_idx) {
409 int pt = (coeff_idx >= 6);
411 if (run >= 9)
413 } else
415 ct = 0;
416 } else if (coeff == 11) { /* end of block */
417 if (coeff_idx == 1) /* first AC coeff ? */
419 break;
420 } else {
422 if (coeff > 4)
423 coeff2 +=
get_bits(&s->gb, coeff <= 9 ? coeff - 4 : 11);
424 ct = 1 + (coeff2 > 1);
426 coeff2 = (coeff2 ^ -sign) + sign;
427 if (coeff_idx)
428 coeff2 *= s->dequant_ac;
430 s->block_coeff[
b][permute[idx]] = coeff2;
431 }
432 }
434 if (coeff_idx >= 64)
435 break;
437 vlc_coeff = &s->ract_vlc[
pt][ct][cg];
438 }
439 }
440 }
441
443 {
447 uint8_t *model1, *model2, *model3;
448 int coeff, sign, coeff_idx;
449 int b, i, cg, idx, ctx;
450 int pt = 0;
/* plane type (0 for Y, 1 for U or V) */
451
452 for (b=0; b<6; b++) {
453 int ct = 1; /* code type */
455
456 if (b > 3) pt = 1;
457
459 + s->above_blocks[s->above_block_idx[
b]].not_null_dc;
462
463 coeff_idx = 0;
464 for (;;) {
466 /* parse a coeff */
473 } else {
476 else
477 coeff = 2;
478 }
479 ct = 2;
480 } else {
481 ct = 1;
482 coeff = 1;
483 }
485 coeff = (coeff ^ -sign) + sign;
486 if (coeff_idx)
487 coeff *= s->dequant_ac;
489 s->block_coeff[
b][permute[idx]] =
coeff;
490 run = 1;
491 } else {
492 /* parse a run */
493 ct = 0;
494 if (coeff_idx > 0) {
496 break;
497
500 if (!run)
501 for (run=9, i=0; i<6; i++)
503 }
504 }
506 if (coeff_idx >= 64)
507 break;
510 }
511
513 s->above_blocks[s->above_block_idx[
b]].not_null_dc = !!s->block_coeff[
b][0];
514 }
515 }
516
518 {
519 int sum = 0, square_sum = 0;
521
522 for (y=0; y<8; y+=2) {
523 for (x=0; x<8; x+=2) {
524 sum += src[x];
525 square_sum += src[x]*src[x];
526 }
527 src += 2*stride;
528 }
529 return (16*square_sum - sum*sum) >> 8;
530 }
531
533 int delta,
const int16_t *weights)
534 {
536
537 for (y=0; y<8; y++) {
538 for (x=0; x<8; x++) {
539 dst[x] = av_clip_uint8(( src[x-delta ] * weights[0]
540 + src[x ] * weights[1]
541 + src[x+delta ] * weights[2]
542 + src[x+2*delta] * weights[3] + 64) >> 7);
543 }
544 src += stride;
545 dst += stride;
546 }
547 }
548
550 int stride,
int h_weight,
int v_weight)
551 {
552 uint8_t *tmp = s->edge_emu_buffer+16;
553 s->h264chroma.put_h264_chroma_pixels_tab[0](tmp,
src, stride, 9, h_weight, 0);
554 s->h264chroma.put_h264_chroma_pixels_tab[0](dst, tmp, stride, 8, 0, v_weight);
555 }
556
558 int offset1,
int offset2,
int stride,
560 {
561 int filter4 = 0;
562 int x8 = mv.
x &
mask;
563 int y8 = mv.
y &
mask;
564
565 if (luma) {
566 x8 *= 2;
567 y8 *= 2;
568 filter4 = s->filter_mode;
569 if (filter4 == 2) {
570 if (s->max_vector_length &&
571 (
FFABS(mv.
x) > s->max_vector_length ||
572 FFABS(mv.
y) > s->max_vector_length)) {
573 filter4 = 0;
574 } else if (s->sample_variance_threshold
576 < s->sample_variance_threshold)) {
577 filter4 = 0;
578 }
579 }
580 }
581
582 if ((y8 && (offset2-offset1)*s->flip<0) || (!y8 && offset1 > offset2)) {
583 offset1 = offset2;
584 }
585
586 if (filter4) {
587 if (!y8) { /* left or right combine */
590 } else if (!x8) { /* above or below combine */
593 } else {
594 s->vp56dsp.vp6_filter_diag4(dst, src+offset1+((mv.
x^mv.
y)>>31), stride,
597 }
598 } else {
599 if (!x8 || !y8) {
600 s->h264chroma.put_h264_chroma_pixels_tab[0](dst, src + offset1, stride, 8, x8, y8);
601 } else {
603 }
604 }
605 }
606
608
610 {
613
617
619
620 if (s->has_alpha) {
621 s->alpha_context =
av_mallocz(
sizeof(VP56Context));
623 s->flip == -1, s->has_alpha);
625 }
626
627 return 0;
628 }
629
631 {
632 s->deblock_filtering = 0;
640 }
641
643
645 {
647
650
651 if (s->alpha_context) {
655 }
656
657 return 0;
658 }
659
661 {
663
664 for (pt=0; pt<2; pt++) {
667 for (ct=0; ct<3; ct++)
668 for (cg=0; cg<6; cg++)
670 }
671 }
672
678 .priv_data_size = sizeof(VP56Context),
683 };
684
685 /* flash version, not flipped upside-down */
691 .priv_data_size = sizeof(VP56Context),
696 };
697
698 /* flash version, not flipped upside-down, with alpha channel */
704 .priv_data_size = sizeof(VP56Context),
709 };