Mercurial > libavcodec.hg
annotate h263.c @ 245:a023f78e0cc2 libavcodec
* simplified indexing
author | kabi |
---|---|
date | Mon, 18 Feb 2002 09:32:11 +0000 |
parents | 4bb6289eff93 |
children | 56ee684c48bb |
rev | line source |
---|---|
0 | 1 /* |
2 * H263/MPEG4 backend for ffmpeg encoder and decoder | |
3 * Copyright (c) 2000,2001 Gerard Lantau. | |
78 | 4 * H263+ support. |
0 | 5 * Copyright (c) 2001 Juan J. Sierralta P. |
6 * | |
7 * This program is free software; you can redistribute it and/or modify | |
8 * it under the terms of the GNU General Public License as published by | |
9 * the Free Software Foundation; either version 2 of the License, or | |
10 * (at your option) any later version. | |
11 * | |
12 * This program is distributed in the hope that it will be useful, | |
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
15 * GNU General Public License for more details. | |
16 * | |
17 * You should have received a copy of the GNU General Public License | |
18 * along with this program; if not, write to the Free Software | |
19 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. | |
20 */ | |
21 #include "common.h" | |
22 #include "dsputil.h" | |
23 #include "avcodec.h" | |
24 #include "mpegvideo.h" | |
25 #include "h263data.h" | |
26 #include "mpeg4data.h" | |
27 | |
28 static void h263_encode_block(MpegEncContext * s, DCTELEM * block, | |
29 int n); | |
30 static void h263_encode_motion(MpegEncContext * s, int val); | |
78 | 31 static void h263p_encode_umotion(MpegEncContext * s, int val); |
0 | 32 static void mpeg4_encode_block(MpegEncContext * s, DCTELEM * block, |
33 int n); | |
34 static int h263_decode_motion(MpegEncContext * s, int pred); | |
78 | 35 static int h263p_decode_umotion(MpegEncContext * s, int pred); |
0 | 36 static int h263_decode_block(MpegEncContext * s, DCTELEM * block, |
37 int n, int coded); | |
38 static int mpeg4_decode_block(MpegEncContext * s, DCTELEM * block, | |
39 int n, int coded); | |
40 | |
41 int h263_get_picture_format(int width, int height) | |
42 { | |
43 int format; | |
44 | |
45 if (width == 128 && height == 96) | |
46 format = 1; | |
47 else if (width == 176 && height == 144) | |
48 format = 2; | |
49 else if (width == 352 && height == 288) | |
50 format = 3; | |
51 else if (width == 704 && height == 576) | |
52 format = 4; | |
53 else if (width == 1408 && height == 1152) | |
54 format = 5; | |
55 else | |
56 format = 7; | |
57 return format; | |
58 } | |
59 | |
60 void h263_encode_picture_header(MpegEncContext * s, int picture_number) | |
61 { | |
78 | 62 int format; |
0 | 63 |
64 align_put_bits(&s->pb); | |
231 | 65 |
66 /* Update the pointer to last GOB */ | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
67 s->ptr_lastgob = pbBufPtr(&s->pb); |
231 | 68 s->gob_number = 0; |
69 | |
70 put_bits(&s->pb, 22, 0x20); /* PSC */ | |
241 | 71 put_bits(&s->pb, 8, (((INT64)s->picture_number * 30 * FRAME_RATE_BASE) / |
0 | 72 s->frame_rate) & 0xff); |
73 | |
74 put_bits(&s->pb, 1, 1); /* marker */ | |
75 put_bits(&s->pb, 1, 0); /* h263 id */ | |
76 put_bits(&s->pb, 1, 0); /* split screen off */ | |
77 put_bits(&s->pb, 1, 0); /* camera off */ | |
78 put_bits(&s->pb, 1, 0); /* freeze picture release off */ | |
78 | 79 |
80 format = h263_get_picture_format(s->width, s->height); | |
0 | 81 if (!s->h263_plus) { |
82 /* H.263v1 */ | |
83 put_bits(&s->pb, 3, format); | |
84 put_bits(&s->pb, 1, (s->pict_type == P_TYPE)); | |
85 /* By now UMV IS DISABLED ON H.263v1, since the restrictions | |
86 of H.263v1 UMV implies to check the predicted MV after | |
87 calculation of the current MB to see if we're on the limits */ | |
88 put_bits(&s->pb, 1, 0); /* unrestricted motion vector: off */ | |
89 put_bits(&s->pb, 1, 0); /* SAC: off */ | |
90 put_bits(&s->pb, 1, 0); /* advanced prediction mode: off */ | |
91 put_bits(&s->pb, 1, 0); /* not PB frame */ | |
92 put_bits(&s->pb, 5, s->qscale); | |
93 put_bits(&s->pb, 1, 0); /* Continuous Presence Multipoint mode: off */ | |
94 } else { | |
95 /* H.263v2 */ | |
96 /* H.263 Plus PTYPE */ | |
97 put_bits(&s->pb, 3, 7); | |
98 put_bits(&s->pb,3,1); /* Update Full Extended PTYPE */ | |
78 | 99 if (format == 7) |
100 put_bits(&s->pb,3,6); /* Custom Source Format */ | |
101 else | |
102 put_bits(&s->pb, 3, format); | |
103 | |
0 | 104 put_bits(&s->pb,1,0); /* Custom PCF: off */ |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
105 s->umvplus = (s->pict_type == P_TYPE) && s->unrestricted_mv; |
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
106 put_bits(&s->pb, 1, s->umvplus); /* Unrestricted Motion Vector */ |
0 | 107 put_bits(&s->pb,1,0); /* SAC: off */ |
108 put_bits(&s->pb,1,0); /* Advanced Prediction Mode: off */ | |
109 put_bits(&s->pb,1,0); /* Advanced Intra Coding: off */ | |
110 put_bits(&s->pb,1,0); /* Deblocking Filter: off */ | |
111 put_bits(&s->pb,1,0); /* Slice Structured: off */ | |
112 put_bits(&s->pb,1,0); /* Reference Picture Selection: off */ | |
113 put_bits(&s->pb,1,0); /* Independent Segment Decoding: off */ | |
114 put_bits(&s->pb,1,0); /* Alternative Inter VLC: off */ | |
115 put_bits(&s->pb,1,0); /* Modified Quantization: off */ | |
116 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */ | |
117 put_bits(&s->pb,3,0); /* Reserved */ | |
118 | |
119 put_bits(&s->pb, 3, s->pict_type == P_TYPE); | |
120 | |
121 put_bits(&s->pb,1,0); /* Reference Picture Resampling: off */ | |
122 put_bits(&s->pb,1,0); /* Reduced-Resolution Update: off */ | |
123 put_bits(&s->pb,1,0); /* Rounding Type */ | |
124 put_bits(&s->pb,2,0); /* Reserved */ | |
125 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */ | |
126 | |
127 /* This should be here if PLUSPTYPE */ | |
128 put_bits(&s->pb, 1, 0); /* Continuous Presence Multipoint mode: off */ | |
129 | |
78 | 130 if (format == 7) { |
131 /* Custom Picture Format (CPFMT) */ | |
0 | 132 |
78 | 133 put_bits(&s->pb,4,2); /* Aspect ratio: CIF 12:11 (4:3) picture */ |
134 put_bits(&s->pb,9,(s->width >> 2) - 1); | |
135 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */ | |
136 put_bits(&s->pb,9,(s->height >> 2)); | |
137 } | |
138 | |
0 | 139 /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */ |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
140 if (s->umvplus) |
0 | 141 put_bits(&s->pb,1,1); /* Limited according tables of Annex D */ |
142 put_bits(&s->pb, 5, s->qscale); | |
143 } | |
144 | |
145 put_bits(&s->pb, 1, 0); /* no PEI */ | |
146 } | |
147 | |
162 | 148 int h263_encode_gob_header(MpegEncContext * s, int mb_line) |
149 { | |
150 int pdif=0; | |
151 | |
152 /* Check to see if we need to put a new GBSC */ | |
153 /* for RTP packetization */ | |
154 if (s->rtp_mode) { | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
155 pdif = pbBufPtr(&s->pb) - s->ptr_lastgob; |
162 | 156 if (pdif >= s->rtp_payload_size) { |
157 /* Bad luck, packet must be cut before */ | |
158 align_put_bits(&s->pb); | |
231 | 159 flush_put_bits(&s->pb); |
160 /* Call the RTP callback to send the last GOB */ | |
161 if (s->rtp_callback) { | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
162 pdif = pbBufPtr(&s->pb) - s->ptr_lastgob; |
231 | 163 s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number); |
164 } | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
165 s->ptr_lastgob = pbBufPtr(&s->pb); |
162 | 166 put_bits(&s->pb, 17, 1); /* GBSC */ |
231 | 167 s->gob_number = mb_line / s->gob_index; |
162 | 168 put_bits(&s->pb, 5, s->gob_number); /* GN */ |
231 | 169 put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */ |
162 | 170 put_bits(&s->pb, 5, s->qscale); /* GQUANT */ |
231 | 171 //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif); |
162 | 172 return pdif; |
173 } else if (pdif + s->mb_line_avgsize >= s->rtp_payload_size) { | |
174 /* Cut the packet before we can't */ | |
175 align_put_bits(&s->pb); | |
231 | 176 flush_put_bits(&s->pb); |
177 /* Call the RTP callback to send the last GOB */ | |
178 if (s->rtp_callback) { | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
179 pdif = pbBufPtr(&s->pb) - s->ptr_lastgob; |
231 | 180 s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number); |
181 } | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
182 s->ptr_lastgob = pbBufPtr(&s->pb); |
162 | 183 put_bits(&s->pb, 17, 1); /* GBSC */ |
231 | 184 s->gob_number = mb_line / s->gob_index; |
162 | 185 put_bits(&s->pb, 5, s->gob_number); /* GN */ |
231 | 186 put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */ |
162 | 187 put_bits(&s->pb, 5, s->qscale); /* GQUANT */ |
231 | 188 //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif); |
162 | 189 return pdif; |
190 } | |
191 } | |
192 return 0; | |
193 } | |
194 | |
0 | 195 void h263_encode_mb(MpegEncContext * s, |
196 DCTELEM block[6][64], | |
197 int motion_x, int motion_y) | |
198 { | |
199 int cbpc, cbpy, i, cbp, pred_x, pred_y; | |
162 | 200 |
0 | 201 // printf("**mb x=%d y=%d\n", s->mb_x, s->mb_y); |
78 | 202 if (!s->mb_intra) { |
203 /* compute cbp */ | |
204 cbp = 0; | |
205 for (i = 0; i < 6; i++) { | |
206 if (s->block_last_index[i] >= 0) | |
207 cbp |= 1 << (5 - i); | |
208 } | |
209 if ((cbp | motion_x | motion_y) == 0) { | |
210 /* skip macroblock */ | |
211 put_bits(&s->pb, 1, 1); | |
212 return; | |
213 } | |
214 put_bits(&s->pb, 1, 0); /* mb coded */ | |
215 cbpc = cbp & 3; | |
216 put_bits(&s->pb, | |
217 inter_MCBPC_bits[cbpc], | |
218 inter_MCBPC_code[cbpc]); | |
219 cbpy = cbp >> 2; | |
220 cbpy ^= 0xf; | |
221 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); | |
0 | 222 |
78 | 223 /* motion vectors: 16x16 mode only now */ |
224 h263_pred_motion(s, 0, &pred_x, &pred_y); | |
225 | |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
226 if (!s->umvplus) { |
78 | 227 h263_encode_motion(s, motion_x - pred_x); |
228 h263_encode_motion(s, motion_y - pred_y); | |
229 } | |
230 else { | |
231 h263p_encode_umotion(s, motion_x - pred_x); | |
232 h263p_encode_umotion(s, motion_y - pred_y); | |
233 if (((motion_x - pred_x) == 1) && ((motion_y - pred_y) == 1)) | |
234 /* To prevent Start Code emulation */ | |
235 put_bits(&s->pb,1,1); | |
236 } | |
237 } else { | |
0 | 238 /* compute cbp */ |
239 cbp = 0; | |
240 for (i = 0; i < 6; i++) { | |
241 if (s->block_last_index[i] >= 1) | |
242 cbp |= 1 << (5 - i); | |
243 } | |
244 | |
245 cbpc = cbp & 3; | |
246 if (s->pict_type == I_TYPE) { | |
247 put_bits(&s->pb, | |
248 intra_MCBPC_bits[cbpc], | |
249 intra_MCBPC_code[cbpc]); | |
250 } else { | |
251 put_bits(&s->pb, 1, 0); /* mb coded */ | |
252 put_bits(&s->pb, | |
253 inter_MCBPC_bits[cbpc + 4], | |
254 inter_MCBPC_code[cbpc + 4]); | |
255 } | |
256 if (s->h263_pred) { | |
257 /* XXX: currently, we do not try to use ac prediction */ | |
258 put_bits(&s->pb, 1, 0); /* no ac prediction */ | |
259 } | |
260 cbpy = cbp >> 2; | |
261 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); | |
262 } | |
263 | |
264 /* encode each block */ | |
265 if (s->h263_pred) { | |
266 for (i = 0; i < 6; i++) { | |
267 mpeg4_encode_block(s, block[i], i); | |
268 } | |
269 } else { | |
270 for (i = 0; i < 6; i++) { | |
271 h263_encode_block(s, block[i], i); | |
272 } | |
273 } | |
274 } | |
275 | |
276 static inline int mid_pred(int a, int b, int c) | |
277 { | |
278 int vmin, vmax; | |
245 | 279 vmax = vmin = a; |
0 | 280 if (b < vmin) |
281 vmin = b; | |
245 | 282 else |
283 vmax = b; | |
284 | |
0 | 285 if (c < vmin) |
286 vmin = c; | |
245 | 287 else if (c > vmax) |
0 | 288 vmax = c; |
289 | |
290 return a + b + c - vmin - vmax; | |
291 } | |
292 | |
293 INT16 *h263_pred_motion(MpegEncContext * s, int block, | |
294 int *px, int *py) | |
295 { | |
245 | 296 int xy, y, wrap; |
0 | 297 INT16 *A, *B, *C, *mot_val; |
298 | |
299 wrap = 2 * s->mb_width + 2; | |
245 | 300 y = xy = 2 * s->mb_y + 1 + ((block >> 1) & 1); // y |
301 xy *= wrap; // y * wrap | |
302 xy += 2 * s->mb_x + 1 + (block & 1); // x + y * wrap | |
0 | 303 |
245 | 304 mot_val = s->motion_val[xy]; |
0 | 305 |
306 /* special case for first line */ | |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
307 if (y == 1 || s->first_slice_line || s->first_gob_line) { |
245 | 308 A = s->motion_val[xy - 1]; |
0 | 309 *px = A[0]; |
310 *py = A[1]; | |
311 } else { | |
312 switch(block) { | |
313 default: | |
314 case 0: | |
245 | 315 A = s->motion_val[xy - 1]; |
316 B = s->motion_val[xy - wrap]; | |
317 C = s->motion_val[xy + 2 - wrap]; | |
0 | 318 break; |
319 case 1: | |
320 case 2: | |
245 | 321 A = s->motion_val[xy - 1]; |
322 B = s->motion_val[xy - wrap]; | |
323 C = s->motion_val[xy + 1 - wrap]; | |
0 | 324 break; |
325 case 3: | |
245 | 326 A = s->motion_val[xy - 1]; |
327 B = s->motion_val[xy - 1 - wrap]; | |
328 C = s->motion_val[xy - wrap]; | |
0 | 329 break; |
330 } | |
331 *px = mid_pred(A[0], B[0], C[0]); | |
332 *py = mid_pred(A[1], B[1], C[1]); | |
333 } | |
334 return mot_val; | |
335 } | |
336 | |
337 | |
338 static void h263_encode_motion(MpegEncContext * s, int val) | |
339 { | |
340 int range, l, m, bit_size, sign, code, bits; | |
341 | |
342 if (val == 0) { | |
343 /* zero vector */ | |
344 code = 0; | |
345 put_bits(&s->pb, mvtab[code][1], mvtab[code][0]); | |
346 } else { | |
347 bit_size = s->f_code - 1; | |
348 range = 1 << bit_size; | |
349 /* modulo encoding */ | |
350 l = range * 32; | |
351 m = 2 * l; | |
352 if (val < -l) { | |
353 val += m; | |
354 } else if (val >= l) { | |
355 val -= m; | |
356 } | |
357 | |
358 if (val >= 0) { | |
359 val--; | |
360 code = (val >> bit_size) + 1; | |
361 bits = val & (range - 1); | |
362 sign = 0; | |
363 } else { | |
364 val = -val; | |
365 val--; | |
366 code = (val >> bit_size) + 1; | |
367 bits = val & (range - 1); | |
368 sign = 1; | |
369 } | |
370 | |
371 put_bits(&s->pb, mvtab[code][1] + 1, (mvtab[code][0] << 1) | sign); | |
372 if (bit_size > 0) { | |
373 put_bits(&s->pb, bit_size, bits); | |
374 } | |
375 } | |
376 } | |
377 | |
78 | 378 /* Encode MV differences on H.263+ with Unrestricted MV mode */ |
379 static void h263p_encode_umotion(MpegEncContext * s, int val) | |
380 { | |
381 short sval = 0; | |
382 short i = 0; | |
383 short n_bits = 0; | |
384 short temp_val; | |
385 int code = 0; | |
386 int tcode; | |
387 | |
388 if ( val == 0) | |
389 put_bits(&s->pb, 1, 1); | |
390 else if (val == 1) | |
391 put_bits(&s->pb, 3, 0); | |
392 else if (val == -1) | |
393 put_bits(&s->pb, 3, 2); | |
394 else { | |
395 | |
396 sval = ((val < 0) ? (short)(-val):(short)val); | |
397 temp_val = sval; | |
398 | |
399 while (temp_val != 0) { | |
400 temp_val = temp_val >> 1; | |
401 n_bits++; | |
402 } | |
403 | |
404 i = n_bits - 1; | |
405 while (i > 0) { | |
406 tcode = (sval & (1 << (i-1))) >> (i-1); | |
407 tcode = (tcode << 1) | 1; | |
408 code = (code << 2) | tcode; | |
409 i--; | |
410 } | |
411 code = ((code << 1) | (val < 0)) << 1; | |
412 put_bits(&s->pb, (2*n_bits)+1, code); | |
413 //printf("\nVal = %d\tCode = %d", sval, code); | |
414 } | |
415 } | |
416 | |
0 | 417 void h263_encode_init_vlc(MpegEncContext *s) |
418 { | |
419 static int done = 0; | |
420 | |
421 if (!done) { | |
422 done = 1; | |
423 init_rl(&rl_inter); | |
424 init_rl(&rl_intra); | |
425 } | |
426 } | |
427 | |
428 static void h263_encode_block(MpegEncContext * s, DCTELEM * block, int n) | |
429 { | |
430 int level, run, last, i, j, last_index, last_non_zero, sign, slevel; | |
431 int code; | |
432 RLTable *rl = &rl_inter; | |
433 | |
434 if (s->mb_intra) { | |
231 | 435 /* DC coef */ |
436 level = block[0]; | |
0 | 437 /* 255 cannot be represented, so we clamp */ |
438 if (level > 254) { | |
439 level = 254; | |
440 block[0] = 254; | |
441 } | |
231 | 442 /* 0 cannot be represented also */ |
443 else if (!level) { | |
444 level = 1; | |
445 block[0] = 1; | |
446 } | |
447 if (level == 128) | |
448 put_bits(&s->pb, 8, 0xff); | |
449 else | |
450 put_bits(&s->pb, 8, level & 0xff); | |
451 i = 1; | |
0 | 452 } else { |
231 | 453 i = 0; |
0 | 454 } |
455 | |
456 /* AC coefs */ | |
457 last_index = s->block_last_index[n]; | |
458 last_non_zero = i - 1; | |
459 for (; i <= last_index; i++) { | |
460 j = zigzag_direct[i]; | |
461 level = block[j]; | |
462 if (level) { | |
463 run = i - last_non_zero - 1; | |
464 last = (i == last_index); | |
465 sign = 0; | |
466 slevel = level; | |
467 if (level < 0) { | |
468 sign = 1; | |
469 level = -level; | |
470 } | |
471 code = get_rl_index(rl, last, run, level); | |
472 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
473 if (code == rl->n) { | |
474 put_bits(&s->pb, 1, last); | |
475 put_bits(&s->pb, 6, run); | |
476 put_bits(&s->pb, 8, slevel & 0xff); | |
477 } else { | |
478 put_bits(&s->pb, 1, sign); | |
479 } | |
480 last_non_zero = i; | |
481 } | |
482 } | |
483 } | |
484 | |
485 /***************************************************/ | |
486 | |
487 /* write mpeg4 VOP header */ | |
488 void mpeg4_encode_picture_header(MpegEncContext * s, int picture_number) | |
489 { | |
490 align_put_bits(&s->pb); | |
491 | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
492 put_bits(&s->pb, 16, 0); /* vop header */ |
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
493 put_bits(&s->pb, 16, 0x1B6); /* vop header */ |
0 | 494 put_bits(&s->pb, 2, s->pict_type - 1); /* pict type: I = 0 , P = 1 */ |
495 /* XXX: time base + 1 not always correct */ | |
496 put_bits(&s->pb, 1, 1); | |
497 put_bits(&s->pb, 1, 0); | |
498 | |
499 put_bits(&s->pb, 1, 1); /* marker */ | |
500 put_bits(&s->pb, 4, 1); /* XXX: correct time increment */ | |
501 put_bits(&s->pb, 1, 1); /* marker */ | |
502 put_bits(&s->pb, 1, 1); /* vop coded */ | |
503 if (s->pict_type == P_TYPE) { | |
504 s->no_rounding = 0; | |
505 put_bits(&s->pb, 1, s->no_rounding); /* rounding type */ | |
506 } | |
507 put_bits(&s->pb, 3, 0); /* intra dc VLC threshold */ | |
508 | |
509 put_bits(&s->pb, 5, s->qscale); | |
510 | |
511 if (s->pict_type != I_TYPE) | |
512 put_bits(&s->pb, 3, s->f_code); /* fcode_for */ | |
513 // printf("****frame %d\n", picture_number); | |
514 } | |
515 | |
516 void h263_dc_scale(MpegEncContext * s) | |
517 { | |
518 int quant; | |
519 | |
520 quant = s->qscale; | |
521 /* luminance */ | |
522 if (quant < 5) | |
523 s->y_dc_scale = 8; | |
524 else if (quant > 4 && quant < 9) | |
525 s->y_dc_scale = (2 * quant); | |
526 else if (quant > 8 && quant < 25) | |
527 s->y_dc_scale = (quant + 8); | |
528 else | |
529 s->y_dc_scale = (2 * quant - 16); | |
530 /* chrominance */ | |
531 if (quant < 5) | |
532 s->c_dc_scale = 8; | |
533 else if (quant > 4 && quant < 25) | |
534 s->c_dc_scale = ((quant + 13) / 2); | |
535 else | |
536 s->c_dc_scale = (quant - 6); | |
537 } | |
538 | |
539 static int mpeg4_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr, int *dir_ptr) | |
540 { | |
245 | 541 int a, b, c, xy, wrap, pred, scale; |
0 | 542 UINT16 *dc_val; |
543 | |
544 /* find prediction */ | |
545 if (n < 4) { | |
546 wrap = s->mb_width * 2 + 2; | |
245 | 547 xy = 2 * s->mb_y + 1 + ((n & 2) >> 1); |
548 xy *= wrap; | |
549 xy += 2 * s->mb_x + 1 + (n & 1); | |
0 | 550 dc_val = s->dc_val[0]; |
551 scale = s->y_dc_scale; | |
552 } else { | |
553 wrap = s->mb_width + 2; | |
245 | 554 xy = s->mb_y + 1; |
555 xy *= wrap; | |
556 xy += s->mb_x + 1; | |
0 | 557 dc_val = s->dc_val[n - 4 + 1]; |
558 scale = s->c_dc_scale; | |
559 } | |
560 | |
561 /* B C | |
562 * A X | |
563 */ | |
245 | 564 a = dc_val[xy - 1]; |
565 b = dc_val[xy - 1 - wrap]; | |
566 c = dc_val[xy - wrap]; | |
0 | 567 |
568 if (abs(a - b) < abs(b - c)) { | |
569 pred = c; | |
570 *dir_ptr = 1; /* top */ | |
571 } else { | |
572 pred = a; | |
573 *dir_ptr = 0; /* left */ | |
574 } | |
575 /* we assume pred is positive */ | |
576 pred = (pred + (scale >> 1)) / scale; | |
577 | |
578 /* prepare address for prediction update */ | |
245 | 579 *dc_val_ptr = &dc_val[xy]; |
0 | 580 |
581 return pred; | |
582 } | |
583 | |
174
ac5075a55488
new IDCT code by Michael Niedermayer (michaelni@gmx.at) - #define SIMPLE_IDCT to enable
arpi_esp
parents:
162
diff
changeset
|
584 void mpeg4_pred_ac(MpegEncContext * s, INT16 *block, int n, |
0 | 585 int dir) |
586 { | |
587 int x, y, wrap, i; | |
588 INT16 *ac_val, *ac_val1; | |
589 | |
590 /* find prediction */ | |
591 if (n < 4) { | |
592 x = 2 * s->mb_x + 1 + (n & 1); | |
593 y = 2 * s->mb_y + 1 + ((n & 2) >> 1); | |
594 wrap = s->mb_width * 2 + 2; | |
595 ac_val = s->ac_val[0][0]; | |
596 } else { | |
597 x = s->mb_x + 1; | |
598 y = s->mb_y + 1; | |
599 wrap = s->mb_width + 2; | |
600 ac_val = s->ac_val[n - 4 + 1][0]; | |
601 } | |
602 ac_val += ((y) * wrap + (x)) * 16; | |
603 ac_val1 = ac_val; | |
604 if (s->ac_pred) { | |
605 if (dir == 0) { | |
606 /* left prediction */ | |
607 ac_val -= 16; | |
608 for(i=1;i<8;i++) { | |
174
ac5075a55488
new IDCT code by Michael Niedermayer (michaelni@gmx.at) - #define SIMPLE_IDCT to enable
arpi_esp
parents:
162
diff
changeset
|
609 block[block_permute_op(i*8)] += ac_val[i]; |
0 | 610 } |
611 } else { | |
612 /* top prediction */ | |
613 ac_val -= 16 * wrap; | |
614 for(i=1;i<8;i++) { | |
174
ac5075a55488
new IDCT code by Michael Niedermayer (michaelni@gmx.at) - #define SIMPLE_IDCT to enable
arpi_esp
parents:
162
diff
changeset
|
615 block[block_permute_op(i)] += ac_val[i + 8]; |
0 | 616 } |
617 } | |
618 } | |
619 /* left copy */ | |
620 for(i=1;i<8;i++) | |
174
ac5075a55488
new IDCT code by Michael Niedermayer (michaelni@gmx.at) - #define SIMPLE_IDCT to enable
arpi_esp
parents:
162
diff
changeset
|
621 ac_val1[i] = block[block_permute_op(i * 8)]; |
0 | 622 /* top copy */ |
623 for(i=1;i<8;i++) | |
174
ac5075a55488
new IDCT code by Michael Niedermayer (michaelni@gmx.at) - #define SIMPLE_IDCT to enable
arpi_esp
parents:
162
diff
changeset
|
624 ac_val1[8 + i] = block[block_permute_op(i)]; |
0 | 625 } |
626 | |
627 static inline void mpeg4_encode_dc(MpegEncContext * s, int level, int n, int *dir_ptr) | |
628 { | |
629 int size, v, pred; | |
630 UINT16 *dc_val; | |
631 | |
632 pred = mpeg4_pred_dc(s, n, &dc_val, dir_ptr); | |
633 if (n < 4) { | |
634 *dc_val = level * s->y_dc_scale; | |
635 } else { | |
636 *dc_val = level * s->c_dc_scale; | |
637 } | |
638 | |
639 /* do the prediction */ | |
640 level -= pred; | |
641 /* find number of bits */ | |
642 size = 0; | |
643 v = abs(level); | |
644 while (v) { | |
645 v >>= 1; | |
646 size++; | |
647 } | |
648 | |
649 if (n < 4) { | |
650 /* luminance */ | |
651 put_bits(&s->pb, DCtab_lum[size][1], DCtab_lum[size][0]); | |
652 } else { | |
653 /* chrominance */ | |
654 put_bits(&s->pb, DCtab_chrom[size][1], DCtab_chrom[size][0]); | |
655 } | |
656 | |
657 /* encode remaining bits */ | |
658 if (size > 0) { | |
659 if (level < 0) | |
660 level = (-level) ^ ((1 << size) - 1); | |
661 put_bits(&s->pb, size, level); | |
662 if (size > 8) | |
663 put_bits(&s->pb, 1, 1); | |
664 } | |
665 } | |
666 | |
667 static void mpeg4_encode_block(MpegEncContext * s, DCTELEM * block, int n) | |
668 { | |
669 int level, run, last, i, j, last_index, last_non_zero, sign, slevel; | |
670 int code, dc_pred_dir; | |
671 const RLTable *rl; | |
672 | |
673 if (s->mb_intra) { | |
674 /* mpeg4 based DC predictor */ | |
675 mpeg4_encode_dc(s, block[0], n, &dc_pred_dir); | |
676 i = 1; | |
677 rl = &rl_intra; | |
678 } else { | |
679 i = 0; | |
680 rl = &rl_inter; | |
681 } | |
682 | |
683 /* AC coefs */ | |
684 last_index = s->block_last_index[n]; | |
685 last_non_zero = i - 1; | |
686 for (; i <= last_index; i++) { | |
687 j = zigzag_direct[i]; | |
688 level = block[j]; | |
689 if (level) { | |
690 run = i - last_non_zero - 1; | |
691 last = (i == last_index); | |
692 sign = 0; | |
693 slevel = level; | |
694 if (level < 0) { | |
695 sign = 1; | |
696 level = -level; | |
697 } | |
698 code = get_rl_index(rl, last, run, level); | |
699 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
700 if (code == rl->n) { | |
701 int level1, run1; | |
702 level1 = level - rl->max_level[last][run]; | |
703 if (level1 < 1) | |
704 goto esc2; | |
705 code = get_rl_index(rl, last, run, level1); | |
706 if (code == rl->n) { | |
707 esc2: | |
708 put_bits(&s->pb, 1, 1); | |
709 if (level > MAX_LEVEL) | |
710 goto esc3; | |
711 run1 = run - rl->max_run[last][level] - 1; | |
712 if (run1 < 0) | |
713 goto esc3; | |
714 code = get_rl_index(rl, last, run1, level); | |
715 if (code == rl->n) { | |
716 esc3: | |
717 /* third escape */ | |
718 put_bits(&s->pb, 1, 1); | |
719 put_bits(&s->pb, 1, last); | |
720 put_bits(&s->pb, 6, run); | |
721 put_bits(&s->pb, 1, 1); | |
722 put_bits(&s->pb, 12, slevel & 0xfff); | |
723 put_bits(&s->pb, 1, 1); | |
724 } else { | |
725 /* second escape */ | |
726 put_bits(&s->pb, 1, 0); | |
727 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
728 put_bits(&s->pb, 1, sign); | |
729 } | |
730 } else { | |
731 /* first escape */ | |
732 put_bits(&s->pb, 1, 0); | |
733 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
734 put_bits(&s->pb, 1, sign); | |
735 } | |
736 } else { | |
737 put_bits(&s->pb, 1, sign); | |
738 } | |
739 last_non_zero = i; | |
740 } | |
741 } | |
742 } | |
743 | |
744 | |
745 | |
746 /***********************************************/ | |
747 /* decoding */ | |
748 | |
749 static VLC intra_MCBPC_vlc; | |
750 static VLC inter_MCBPC_vlc; | |
751 static VLC cbpy_vlc; | |
752 static VLC mv_vlc; | |
753 static VLC dc_lum, dc_chrom; | |
754 | |
755 void init_rl(RLTable *rl) | |
756 { | |
757 INT8 max_level[MAX_RUN+1], max_run[MAX_LEVEL+1]; | |
758 UINT8 index_run[MAX_RUN+1]; | |
759 int last, run, level, start, end, i; | |
760 | |
761 /* compute max_level[], max_run[] and index_run[] */ | |
762 for(last=0;last<2;last++) { | |
763 if (last == 0) { | |
764 start = 0; | |
765 end = rl->last; | |
766 } else { | |
767 start = rl->last; | |
768 end = rl->n; | |
769 } | |
770 | |
771 memset(max_level, 0, MAX_RUN + 1); | |
772 memset(max_run, 0, MAX_LEVEL + 1); | |
773 memset(index_run, rl->n, MAX_RUN + 1); | |
774 for(i=start;i<end;i++) { | |
775 run = rl->table_run[i]; | |
776 level = rl->table_level[i]; | |
777 if (index_run[run] == rl->n) | |
778 index_run[run] = i; | |
779 if (level > max_level[run]) | |
780 max_level[run] = level; | |
781 if (run > max_run[level]) | |
782 max_run[level] = run; | |
783 } | |
784 rl->max_level[last] = malloc(MAX_RUN + 1); | |
785 memcpy(rl->max_level[last], max_level, MAX_RUN + 1); | |
786 rl->max_run[last] = malloc(MAX_LEVEL + 1); | |
787 memcpy(rl->max_run[last], max_run, MAX_LEVEL + 1); | |
788 rl->index_run[last] = malloc(MAX_RUN + 1); | |
789 memcpy(rl->index_run[last], index_run, MAX_RUN + 1); | |
790 } | |
791 } | |
792 | |
793 void init_vlc_rl(RLTable *rl) | |
794 { | |
795 init_vlc(&rl->vlc, 9, rl->n + 1, | |
796 &rl->table_vlc[0][1], 4, 2, | |
797 &rl->table_vlc[0][0], 4, 2); | |
798 } | |
799 | |
800 /* init vlcs */ | |
801 | |
802 /* XXX: find a better solution to handle static init */ | |
803 void h263_decode_init_vlc(MpegEncContext *s) | |
804 { | |
805 static int done = 0; | |
806 | |
807 if (!done) { | |
808 done = 1; | |
809 | |
810 init_vlc(&intra_MCBPC_vlc, 6, 8, | |
811 intra_MCBPC_bits, 1, 1, | |
812 intra_MCBPC_code, 1, 1); | |
161 | 813 init_vlc(&inter_MCBPC_vlc, 9, 25, |
0 | 814 inter_MCBPC_bits, 1, 1, |
815 inter_MCBPC_code, 1, 1); | |
816 init_vlc(&cbpy_vlc, 6, 16, | |
817 &cbpy_tab[0][1], 2, 1, | |
818 &cbpy_tab[0][0], 2, 1); | |
819 init_vlc(&mv_vlc, 9, 33, | |
820 &mvtab[0][1], 2, 1, | |
821 &mvtab[0][0], 2, 1); | |
822 init_rl(&rl_inter); | |
823 init_rl(&rl_intra); | |
824 init_vlc_rl(&rl_inter); | |
825 init_vlc_rl(&rl_intra); | |
826 init_vlc(&dc_lum, 9, 13, | |
827 &DCtab_lum[0][1], 2, 1, | |
828 &DCtab_lum[0][0], 2, 1); | |
829 init_vlc(&dc_chrom, 9, 13, | |
830 &DCtab_chrom[0][1], 2, 1, | |
831 &DCtab_chrom[0][0], 2, 1); | |
832 } | |
833 } | |
834 | |
162 | 835 int h263_decode_gob_header(MpegEncContext *s) |
836 { | |
837 unsigned int val, gfid; | |
838 | |
839 /* Check for GOB Start Code */ | |
840 val = show_bits(&s->gb, 16); | |
841 if (val == 0) { | |
842 /* We have a GBSC probably with GSTUFF */ | |
843 skip_bits(&s->gb, 16); /* Drop the zeros */ | |
844 while (get_bits1(&s->gb) == 0); /* Seek the '1' bit */ | |
845 #ifdef DEBUG | |
846 fprintf(stderr,"\nGOB Start Code at MB %d\n", (s->mb_y * s->mb_width) + s->mb_x); | |
847 #endif | |
848 s->gob_number = get_bits(&s->gb, 5); /* GN */ | |
849 gfid = get_bits(&s->gb, 2); /* GFID */ | |
850 s->qscale = get_bits(&s->gb, 5); /* GQUANT */ | |
851 #ifdef DEBUG | |
234
5fc0c3af3fe4
alternative bitstream writer (disabled by default, uncomment #define ALT_BISTREAM_WRITER in common.h if u want to try it)
michaelni
parents:
231
diff
changeset
|
852 fprintf(stderr, "\nGN: %u GFID: %u Quant: %u\n", s->gob_number, gfid, s->qscale); |
162 | 853 #endif |
854 return 1; | |
855 } | |
856 return 0; | |
857 | |
858 } | |
859 | |
0 | 860 int h263_decode_mb(MpegEncContext *s, |
861 DCTELEM block[6][64]) | |
862 { | |
863 int cbpc, cbpy, i, cbp, pred_x, pred_y, mx, my, dquant; | |
864 INT16 *mot_val; | |
110 | 865 static INT8 quant_tab[4] = { -1, -2, 1, 2 }; |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
866 |
0 | 867 if (s->pict_type == P_TYPE) { |
21 | 868 if (get_bits1(&s->gb)) { |
0 | 869 /* skip mb */ |
870 s->mb_intra = 0; | |
871 for(i=0;i<6;i++) | |
872 s->block_last_index[i] = -1; | |
873 s->mv_dir = MV_DIR_FORWARD; | |
874 s->mv_type = MV_TYPE_16X16; | |
875 s->mv[0][0][0] = 0; | |
876 s->mv[0][0][1] = 0; | |
7
1d3ac9654178
added skip macroblock optimization (big perf win on black regions for example)
glantau
parents:
0
diff
changeset
|
877 s->mb_skiped = 1; |
0 | 878 return 0; |
879 } | |
880 cbpc = get_vlc(&s->gb, &inter_MCBPC_vlc); | |
144 | 881 //fprintf(stderr, "\tCBPC: %d", cbpc); |
0 | 882 if (cbpc < 0) |
883 return -1; | |
161 | 884 if (cbpc > 20) |
885 cbpc+=3; | |
886 else if (cbpc == 20) | |
887 fprintf(stderr, "Stuffing !"); | |
144 | 888 |
0 | 889 dquant = cbpc & 8; |
890 s->mb_intra = ((cbpc & 4) != 0); | |
891 } else { | |
892 cbpc = get_vlc(&s->gb, &intra_MCBPC_vlc); | |
893 if (cbpc < 0) | |
894 return -1; | |
895 dquant = cbpc & 4; | |
896 s->mb_intra = 1; | |
897 } | |
898 | |
899 if (!s->mb_intra) { | |
900 cbpy = get_vlc(&s->gb, &cbpy_vlc); | |
901 cbp = (cbpc & 3) | ((cbpy ^ 0xf) << 2); | |
902 if (dquant) { | |
903 s->qscale += quant_tab[get_bits(&s->gb, 2)]; | |
904 if (s->qscale < 1) | |
905 s->qscale = 1; | |
906 else if (s->qscale > 31) | |
907 s->qscale = 31; | |
908 } | |
909 s->mv_dir = MV_DIR_FORWARD; | |
910 if ((cbpc & 16) == 0) { | |
911 /* 16x16 motion prediction */ | |
912 s->mv_type = MV_TYPE_16X16; | |
913 h263_pred_motion(s, 0, &pred_x, &pred_y); | |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
914 if (s->umvplus_dec) |
78 | 915 mx = h263p_decode_umotion(s, pred_x); |
916 else | |
917 mx = h263_decode_motion(s, pred_x); | |
0 | 918 if (mx >= 0xffff) |
919 return -1; | |
78 | 920 |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
921 if (s->umvplus_dec) |
78 | 922 my = h263p_decode_umotion(s, pred_y); |
923 else | |
924 my = h263_decode_motion(s, pred_y); | |
0 | 925 if (my >= 0xffff) |
926 return -1; | |
927 s->mv[0][0][0] = mx; | |
928 s->mv[0][0][1] = my; | |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
929 /*fprintf(stderr, "\n MB %d", (s->mb_y * s->mb_width) + s->mb_x); |
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
930 fprintf(stderr, "\n\tmvx: %d\t\tpredx: %d", mx, pred_x); |
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
931 fprintf(stderr, "\n\tmvy: %d\t\tpredy: %d", my, pred_y);*/ |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
932 if (s->umvplus_dec && (mx - pred_x) == 1 && (my - pred_y) == 1) |
78 | 933 skip_bits1(&s->gb); /* Bit stuffing to prevent PSC */ |
934 | |
0 | 935 } else { |
936 s->mv_type = MV_TYPE_8X8; | |
937 for(i=0;i<4;i++) { | |
938 mot_val = h263_pred_motion(s, i, &pred_x, &pred_y); | |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
939 if (s->umvplus_dec) |
78 | 940 mx = h263p_decode_umotion(s, pred_x); |
941 else | |
942 mx = h263_decode_motion(s, pred_x); | |
0 | 943 if (mx >= 0xffff) |
944 return -1; | |
78 | 945 |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
946 if (s->umvplus_dec) |
78 | 947 my = h263p_decode_umotion(s, pred_y); |
948 else | |
949 my = h263_decode_motion(s, pred_y); | |
0 | 950 if (my >= 0xffff) |
951 return -1; | |
952 s->mv[0][i][0] = mx; | |
953 s->mv[0][i][1] = my; | |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
954 if (s->umvplus_dec && (mx - pred_x) == 1 && (my - pred_y) == 1) |
78 | 955 skip_bits1(&s->gb); /* Bit stuffing to prevent PSC */ |
0 | 956 mot_val[0] = mx; |
957 mot_val[1] = my; | |
958 } | |
959 } | |
960 } else { | |
961 s->ac_pred = 0; | |
144 | 962 if (s->h263_pred) { |
21 | 963 s->ac_pred = get_bits1(&s->gb); |
0 | 964 } |
965 cbpy = get_vlc(&s->gb, &cbpy_vlc); | |
966 cbp = (cbpc & 3) | (cbpy << 2); | |
967 if (dquant) { | |
968 s->qscale += quant_tab[get_bits(&s->gb, 2)]; | |
969 if (s->qscale < 1) | |
970 s->qscale = 1; | |
971 else if (s->qscale > 31) | |
972 s->qscale = 31; | |
973 } | |
974 } | |
975 | |
976 /* decode each block */ | |
977 if (s->h263_pred) { | |
978 for (i = 0; i < 6; i++) { | |
979 if (mpeg4_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1) < 0) | |
980 return -1; | |
981 } | |
982 } else { | |
983 for (i = 0; i < 6; i++) { | |
984 if (h263_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1) < 0) | |
985 return -1; | |
986 } | |
987 } | |
988 return 0; | |
989 } | |
990 | |
991 static int h263_decode_motion(MpegEncContext * s, int pred) | |
992 { | |
993 int code, val, sign, shift, l, m; | |
994 | |
995 code = get_vlc(&s->gb, &mv_vlc); | |
996 if (code < 0) | |
997 return 0xffff; | |
998 | |
999 if (code == 0) | |
1000 return pred; | |
21 | 1001 sign = get_bits1(&s->gb); |
0 | 1002 shift = s->f_code - 1; |
1003 val = (code - 1) << shift; | |
1004 if (shift > 0) | |
1005 val |= get_bits(&s->gb, shift); | |
1006 val++; | |
1007 if (sign) | |
1008 val = -val; | |
1009 val += pred; | |
1010 | |
1011 /* modulo decoding */ | |
1012 if (!s->h263_long_vectors) { | |
1013 l = (1 << (s->f_code - 1)) * 32; | |
1014 m = 2 * l; | |
1015 if (val < -l) { | |
1016 val += m; | |
1017 } else if (val >= l) { | |
1018 val -= m; | |
1019 } | |
1020 } else { | |
1021 /* horrible h263 long vector mode */ | |
1022 if (pred < -31 && val < -63) | |
1023 val += 64; | |
1024 if (pred > 32 && val > 63) | |
1025 val -= 64; | |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
1026 |
0 | 1027 } |
1028 return val; | |
1029 } | |
1030 | |
78 | 1031 /* Decodes RVLC of H.263+ UMV */ |
1032 static int h263p_decode_umotion(MpegEncContext * s, int pred) | |
1033 { | |
1034 int code = 0, sign; | |
1035 | |
1036 if (get_bits1(&s->gb)) /* Motion difference = 0 */ | |
1037 return pred; | |
1038 | |
1039 code = 2 + get_bits1(&s->gb); | |
1040 | |
1041 while (get_bits1(&s->gb)) | |
1042 { | |
1043 code <<= 1; | |
1044 code += get_bits1(&s->gb); | |
1045 } | |
1046 sign = code & 1; | |
1047 code >>= 1; | |
1048 | |
1049 code = (sign) ? (pred - code) : (pred + code); | |
1050 #ifdef DEBUG | |
1051 fprintf(stderr,"H.263+ UMV Motion = %d\n", code); | |
1052 #endif | |
1053 return code; | |
1054 | |
1055 } | |
1056 | |
0 | 1057 static int h263_decode_block(MpegEncContext * s, DCTELEM * block, |
1058 int n, int coded) | |
1059 { | |
1060 int code, level, i, j, last, run; | |
1061 RLTable *rl = &rl_inter; | |
1062 | |
1063 if (s->mb_intra) { | |
1064 /* DC coef */ | |
1065 if (s->h263_rv10 && s->rv10_version == 3 && s->pict_type == I_TYPE) { | |
1066 int component, diff; | |
1067 component = (n <= 3 ? 0 : n - 4 + 1); | |
1068 level = s->last_dc[component]; | |
1069 if (s->rv10_first_dc_coded[component]) { | |
1070 diff = rv_decode_dc(s, n); | |
1071 if (diff == 0xffff) | |
1072 return -1; | |
1073 level += diff; | |
1074 level = level & 0xff; /* handle wrap round */ | |
1075 s->last_dc[component] = level; | |
1076 } else { | |
1077 s->rv10_first_dc_coded[component] = 1; | |
1078 } | |
1079 } else { | |
1080 level = get_bits(&s->gb, 8); | |
1081 if (level == 255) | |
1082 level = 128; | |
1083 } | |
1084 block[0] = level; | |
1085 i = 1; | |
1086 } else { | |
1087 i = 0; | |
1088 } | |
1089 if (!coded) { | |
1090 s->block_last_index[n] = i - 1; | |
1091 return 0; | |
1092 } | |
1093 | |
1094 for(;;) { | |
1095 code = get_vlc(&s->gb, &rl->vlc); | |
1096 if (code < 0) | |
1097 return -1; | |
1098 if (code == rl->n) { | |
1099 /* escape */ | |
21 | 1100 last = get_bits1(&s->gb); |
0 | 1101 run = get_bits(&s->gb, 6); |
1102 level = (INT8)get_bits(&s->gb, 8); | |
1103 if (s->h263_rv10 && level == -128) { | |
1104 /* XXX: should patch encoder too */ | |
1105 level = get_bits(&s->gb, 12); | |
1106 level = (level << 20) >> 20; | |
1107 } | |
1108 } else { | |
1109 run = rl->table_run[code]; | |
1110 level = rl->table_level[code]; | |
1111 last = code >= rl->last; | |
21 | 1112 if (get_bits1(&s->gb)) |
0 | 1113 level = -level; |
1114 } | |
1115 i += run; | |
1116 if (i >= 64) | |
1117 return -1; | |
1118 j = zigzag_direct[i]; | |
1119 block[j] = level; | |
1120 if (last) | |
1121 break; | |
1122 i++; | |
1123 } | |
1124 s->block_last_index[n] = i; | |
1125 return 0; | |
1126 } | |
1127 | |
1128 static int mpeg4_decode_dc(MpegEncContext * s, int n, int *dir_ptr) | |
1129 { | |
1130 int level, pred, code; | |
1131 UINT16 *dc_val; | |
1132 | |
1133 if (n < 4) | |
1134 code = get_vlc(&s->gb, &dc_lum); | |
1135 else | |
1136 code = get_vlc(&s->gb, &dc_chrom); | |
1137 if (code < 0) | |
1138 return -1; | |
1139 if (code == 0) { | |
1140 level = 0; | |
1141 } else { | |
1142 level = get_bits(&s->gb, code); | |
1143 if ((level >> (code - 1)) == 0) /* if MSB not set it is negative*/ | |
1144 level = - (level ^ ((1 << code) - 1)); | |
1145 if (code > 8) | |
21 | 1146 skip_bits1(&s->gb); /* marker */ |
0 | 1147 } |
1148 | |
1149 pred = mpeg4_pred_dc(s, n, &dc_val, dir_ptr); | |
1150 level += pred; | |
1151 if (level < 0) | |
1152 level = 0; | |
1153 if (n < 4) { | |
1154 *dc_val = level * s->y_dc_scale; | |
1155 } else { | |
1156 *dc_val = level * s->c_dc_scale; | |
1157 } | |
1158 return level; | |
1159 } | |
1160 | |
1161 static int mpeg4_decode_block(MpegEncContext * s, DCTELEM * block, | |
1162 int n, int coded) | |
1163 { | |
1164 int code, level, i, j, last, run; | |
1165 int dc_pred_dir; | |
1166 RLTable *rl; | |
1167 const UINT8 *scan_table; | |
1168 | |
1169 if (s->mb_intra) { | |
1170 /* DC coef */ | |
1171 level = mpeg4_decode_dc(s, n, &dc_pred_dir); | |
1172 if (level < 0) | |
1173 return -1; | |
1174 block[0] = level; | |
1175 i = 1; | |
1176 if (!coded) | |
1177 goto not_coded; | |
1178 rl = &rl_intra; | |
1179 if (s->ac_pred) { | |
1180 if (dc_pred_dir == 0) | |
1181 scan_table = ff_alternate_vertical_scan; /* left */ | |
1182 else | |
1183 scan_table = ff_alternate_horizontal_scan; /* top */ | |
1184 } else { | |
1185 scan_table = zigzag_direct; | |
1186 } | |
1187 } else { | |
1188 i = 0; | |
1189 if (!coded) { | |
1190 s->block_last_index[n] = i - 1; | |
1191 return 0; | |
1192 } | |
1193 rl = &rl_inter; | |
1194 scan_table = zigzag_direct; | |
1195 } | |
1196 | |
1197 for(;;) { | |
1198 code = get_vlc(&s->gb, &rl->vlc); | |
1199 if (code < 0) | |
1200 return -1; | |
1201 if (code == rl->n) { | |
1202 /* escape */ | |
21 | 1203 if (get_bits1(&s->gb) != 0) { |
1204 if (get_bits1(&s->gb) != 0) { | |
0 | 1205 /* third escape */ |
21 | 1206 last = get_bits1(&s->gb); |
0 | 1207 run = get_bits(&s->gb, 6); |
21 | 1208 get_bits1(&s->gb); /* marker */ |
0 | 1209 level = get_bits(&s->gb, 12); |
1210 level = (level << 20) >> 20; /* sign extend */ | |
21 | 1211 skip_bits1(&s->gb); /* marker */ |
0 | 1212 } else { |
1213 /* second escape */ | |
1214 code = get_vlc(&s->gb, &rl->vlc); | |
1215 if (code < 0 || code >= rl->n) | |
1216 return -1; | |
1217 run = rl->table_run[code]; | |
1218 level = rl->table_level[code]; | |
1219 last = code >= rl->last; | |
1220 run += rl->max_run[last][level] + 1; | |
21 | 1221 if (get_bits1(&s->gb)) |
0 | 1222 level = -level; |
1223 } | |
1224 } else { | |
1225 /* first escape */ | |
1226 code = get_vlc(&s->gb, &rl->vlc); | |
1227 if (code < 0 || code >= rl->n) | |
1228 return -1; | |
1229 run = rl->table_run[code]; | |
1230 level = rl->table_level[code]; | |
1231 last = code >= rl->last; | |
1232 level += rl->max_level[last][run]; | |
21 | 1233 if (get_bits1(&s->gb)) |
0 | 1234 level = -level; |
1235 } | |
1236 } else { | |
1237 run = rl->table_run[code]; | |
1238 level = rl->table_level[code]; | |
1239 last = code >= rl->last; | |
21 | 1240 if (get_bits1(&s->gb)) |
0 | 1241 level = -level; |
1242 } | |
1243 i += run; | |
1244 if (i >= 64) | |
1245 return -1; | |
1246 j = scan_table[i]; | |
1247 block[j] = level; | |
1248 i++; | |
1249 if (last) | |
1250 break; | |
1251 } | |
1252 not_coded: | |
1253 if (s->mb_intra) { | |
1254 mpeg4_pred_ac(s, block, n, dc_pred_dir); | |
1255 if (s->ac_pred) { | |
1256 i = 64; /* XXX: not optimal */ | |
1257 } | |
1258 } | |
1259 s->block_last_index[n] = i - 1; | |
1260 return 0; | |
1261 } | |
1262 | |
1263 /* most is hardcoded. should extend to handle all h263 streams */ | |
1264 int h263_decode_picture_header(MpegEncContext *s) | |
1265 { | |
1266 int format, width, height; | |
1267 | |
1268 /* picture header */ | |
1269 if (get_bits(&s->gb, 22) != 0x20) | |
1270 return -1; | |
231 | 1271 s->picture_number = get_bits(&s->gb, 8); /* picture timestamp */ |
1272 | |
21 | 1273 if (get_bits1(&s->gb) != 1) |
0 | 1274 return -1; /* marker */ |
21 | 1275 if (get_bits1(&s->gb) != 0) |
0 | 1276 return -1; /* h263 id */ |
21 | 1277 skip_bits1(&s->gb); /* split screen off */ |
1278 skip_bits1(&s->gb); /* camera off */ | |
1279 skip_bits1(&s->gb); /* freeze picture release off */ | |
0 | 1280 |
155
3c3449bce692
- Bug fix on MV prediction for MPEG4 caused by new H.263 GOB code.
pulento
parents:
154
diff
changeset
|
1281 /* Reset GOB number */ |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
1282 s->gob_number = 0; |
155
3c3449bce692
- Bug fix on MV prediction for MPEG4 caused by new H.263 GOB code.
pulento
parents:
154
diff
changeset
|
1283 |
0 | 1284 format = get_bits(&s->gb, 3); |
1285 | |
161 | 1286 if (format != 7 && format != 6) { |
0 | 1287 s->h263_plus = 0; |
1288 /* H.263v1 */ | |
1289 width = h263_format[format][0]; | |
1290 height = h263_format[format][1]; | |
1291 if (!width) | |
1292 return -1; | |
161 | 1293 |
1294 s->width = width; | |
1295 s->height = height; | |
21 | 1296 s->pict_type = I_TYPE + get_bits1(&s->gb); |
0 | 1297 |
21 | 1298 s->unrestricted_mv = get_bits1(&s->gb); |
0 | 1299 s->h263_long_vectors = s->unrestricted_mv; |
1300 | |
21 | 1301 if (get_bits1(&s->gb) != 0) |
0 | 1302 return -1; /* SAC: off */ |
161 | 1303 if (get_bits1(&s->gb) != 0) { |
1304 s->mv_type = MV_TYPE_8X8; /* Advanced prediction mode */ | |
1305 } | |
1306 | |
21 | 1307 if (get_bits1(&s->gb) != 0) |
0 | 1308 return -1; /* not PB frame */ |
1309 | |
1310 s->qscale = get_bits(&s->gb, 5); | |
21 | 1311 skip_bits1(&s->gb); /* Continuous Presence Multipoint mode: off */ |
0 | 1312 } else { |
161 | 1313 int ufep; |
1314 | |
0 | 1315 /* H.263v2 */ |
161 | 1316 s->h263_plus = 1; |
1317 ufep = get_bits(&s->gb, 3); /* Update Full Extended PTYPE */ | |
1318 | |
1319 if (ufep == 1) { | |
1320 /* OPPTYPE */ | |
1321 format = get_bits(&s->gb, 3); | |
1322 skip_bits(&s->gb,1); /* Custom PCF */ | |
1323 s->umvplus_dec = get_bits(&s->gb, 1); /* Unrestricted Motion Vector */ | |
1324 skip_bits1(&s->gb); /* Syntax-based Arithmetic Coding (SAC) */ | |
1325 if (get_bits1(&s->gb) != 0) { | |
1326 s->mv_type = MV_TYPE_8X8; /* Advanced prediction mode */ | |
1327 } | |
1328 skip_bits(&s->gb, 8); | |
1329 skip_bits(&s->gb, 3); /* Reserved */ | |
1330 } else if (ufep != 0) | |
0 | 1331 return -1; |
161 | 1332 |
78 | 1333 /* MPPTYPE */ |
0 | 1334 s->pict_type = get_bits(&s->gb, 3) + 1; |
1335 if (s->pict_type != I_TYPE && | |
1336 s->pict_type != P_TYPE) | |
1337 return -1; | |
21 | 1338 skip_bits(&s->gb, 7); |
78 | 1339 |
1340 /* Get the picture dimensions */ | |
161 | 1341 if (ufep) { |
1342 if (format == 6) { | |
1343 /* Custom Picture Format (CPFMT) */ | |
1344 skip_bits(&s->gb, 4); /* aspect ratio */ | |
1345 width = (get_bits(&s->gb, 9) + 1) * 4; | |
1346 skip_bits1(&s->gb); | |
1347 height = get_bits(&s->gb, 9) * 4; | |
78 | 1348 #ifdef DEBUG |
161 | 1349 fprintf(stderr,"\nH.263+ Custom picture: %dx%d\n",width,height); |
78 | 1350 #endif |
161 | 1351 } |
1352 else { | |
1353 width = h263_format[format][0]; | |
1354 height = h263_format[format][1]; | |
1355 } | |
1356 if ((width == 0) || (height == 0)) | |
1357 return -1; | |
1358 s->width = width; | |
1359 s->height = height; | |
1360 if (s->umvplus_dec) { | |
1361 skip_bits1(&s->gb); /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */ | |
1362 } | |
78 | 1363 } |
1364 | |
0 | 1365 s->qscale = get_bits(&s->gb, 5); |
1366 } | |
1367 /* PEI */ | |
21 | 1368 while (get_bits1(&s->gb) != 0) { |
1369 skip_bits(&s->gb, 8); | |
0 | 1370 } |
1371 s->f_code = 1; | |
1372 return 0; | |
1373 } | |
1374 | |
1375 /* decode mpeg4 VOP header */ | |
1376 int mpeg4_decode_picture_header(MpegEncContext * s) | |
1377 { | |
1378 int time_incr, startcode, state, v; | |
1379 | |
1380 redo: | |
1381 /* search next start code */ | |
1382 align_get_bits(&s->gb); | |
1383 state = 0xff; | |
1384 for(;;) { | |
1385 v = get_bits(&s->gb, 8); | |
1386 if (state == 0x000001) { | |
1387 state = ((state << 8) | v) & 0xffffff; | |
1388 startcode = state; | |
1389 break; | |
1390 } | |
1391 state = ((state << 8) | v) & 0xffffff; | |
1392 /* XXX: really detect end of frame */ | |
1393 if (state == 0) | |
1394 return -1; | |
1395 } | |
1396 | |
1397 if (startcode == 0x120) { | |
63 | 1398 int time_increment_resolution, width, height, vo_ver_id; |
0 | 1399 |
1400 /* vol header */ | |
21 | 1401 skip_bits(&s->gb, 1); /* random access */ |
1402 skip_bits(&s->gb, 8); /* vo_type */ | |
63 | 1403 if (get_bits1(&s->gb) != 0) { /* is_ol_id */ |
1404 vo_ver_id = get_bits(&s->gb, 4); /* vo_ver_id */ | |
1405 skip_bits(&s->gb, 3); /* vo_priority */ | |
1406 } else { | |
1407 vo_ver_id = 1; | |
1408 } | |
0 | 1409 |
21 | 1410 skip_bits(&s->gb, 4); /* aspect_ratio_info */ |
63 | 1411 skip_bits1(&s->gb); /* vol control parameter */ |
1412 s->shape = get_bits(&s->gb, 2); /* vol shape */ | |
21 | 1413 skip_bits1(&s->gb); /* marker */ |
0 | 1414 |
1415 time_increment_resolution = get_bits(&s->gb, 16); | |
70 | 1416 s->time_increment_bits = av_log2(time_increment_resolution - 1) + 1; |
63 | 1417 if (s->time_increment_bits < 1) |
1418 s->time_increment_bits = 1; | |
21 | 1419 skip_bits1(&s->gb); /* marker */ |
0 | 1420 |
63 | 1421 if (get_bits1(&s->gb) != 0) { /* fixed_vop_rate */ |
1422 skip_bits(&s->gb, s->time_increment_bits); | |
1423 } | |
0 | 1424 |
63 | 1425 if (s->shape != 2) { |
1426 if (s->shape == 0) { | |
1427 skip_bits1(&s->gb); /* marker */ | |
1428 width = get_bits(&s->gb, 13); | |
1429 skip_bits1(&s->gb); /* marker */ | |
1430 height = get_bits(&s->gb, 13); | |
1431 skip_bits1(&s->gb); /* marker */ | |
1432 } | |
1433 | |
1434 skip_bits1(&s->gb); /* interlaced */ | |
1435 skip_bits1(&s->gb); /* OBMC */ | |
1436 if (vo_ver_id == 1) { | |
1437 s->vol_sprite_usage = get_bits1(&s->gb); /* vol_sprite_usage */ | |
1438 } else { | |
1439 s->vol_sprite_usage = get_bits(&s->gb, 2); /* vol_sprite_usage */ | |
1440 } | |
1441 if (get_bits1(&s->gb) == 1) { /* not_8_bit */ | |
1442 s->quant_precision = get_bits(&s->gb, 4); /* quant_precision */ | |
1443 skip_bits(&s->gb, 4); /* bits_per_pixel */ | |
1444 } else { | |
1445 s->quant_precision = 5; | |
1446 } | |
1447 | |
1448 skip_bits1(&s->gb); /* vol_quant_type */ | |
1449 skip_bits1(&s->gb); /* vol_quarter_pixel */ | |
1450 skip_bits1(&s->gb); /* complexity_estimation_disabled */ | |
1451 skip_bits1(&s->gb); /* resync_marker_disabled */ | |
1452 skip_bits1(&s->gb); /* data_partioning_enabled */ | |
1453 if (get_bits1(&s->gb) != 0) { /* scalability */ | |
1454 printf("bad scalability!!!\n"); | |
1455 return -1; | |
1456 } | |
1457 } | |
0 | 1458 goto redo; |
1459 } else if (startcode != 0x1b6) { | |
1460 goto redo; | |
1461 } | |
1462 | |
1463 s->pict_type = get_bits(&s->gb, 2) + 1; /* pict type: I = 0 , P = 1 */ | |
1464 if (s->pict_type != I_TYPE && | |
1465 s->pict_type != P_TYPE) | |
1466 return -1; | |
1467 | |
1468 /* XXX: parse time base */ | |
1469 time_incr = 0; | |
21 | 1470 while (get_bits1(&s->gb) != 0) |
0 | 1471 time_incr++; |
1472 | |
21 | 1473 skip_bits1(&s->gb); /* marker */ |
1474 skip_bits(&s->gb, s->time_increment_bits); | |
1475 skip_bits1(&s->gb); /* marker */ | |
0 | 1476 /* vop coded */ |
21 | 1477 if (get_bits1(&s->gb) != 1) |
63 | 1478 goto redo; |
0 | 1479 |
63 | 1480 if (s->shape != 2 && s->pict_type == P_TYPE) { |
0 | 1481 /* rounding type for motion estimation */ |
21 | 1482 s->no_rounding = get_bits1(&s->gb); |
63 | 1483 } else { |
1484 s->no_rounding = 0; | |
0 | 1485 } |
63 | 1486 |
1487 if (s->shape != 0) { | |
1488 if (s->vol_sprite_usage != 1 || s->pict_type != I_TYPE) { | |
1489 int width, height, hor_spat_ref, ver_spat_ref; | |
1490 | |
1491 width = get_bits(&s->gb, 13); | |
1492 skip_bits1(&s->gb); /* marker */ | |
1493 height = get_bits(&s->gb, 13); | |
1494 skip_bits1(&s->gb); /* marker */ | |
1495 hor_spat_ref = get_bits(&s->gb, 13); /* hor_spat_ref */ | |
1496 skip_bits1(&s->gb); /* marker */ | |
1497 ver_spat_ref = get_bits(&s->gb, 13); /* ver_spat_ref */ | |
1498 } | |
1499 skip_bits1(&s->gb); /* change_CR_disable */ | |
1500 | |
1501 if (get_bits1(&s->gb) != 0) { | |
1502 skip_bits(&s->gb, 8); /* constant_alpha_value */ | |
1503 } | |
1504 } | |
0 | 1505 |
63 | 1506 if (s->shape != 2) { |
1507 skip_bits(&s->gb, 3); /* intra dc VLC threshold */ | |
1508 | |
1509 /* note: we do not use quant_precision to avoid problem if no | |
1510 MPEG4 vol header as it is found on some old opendivx | |
1511 movies */ | |
1512 s->qscale = get_bits(&s->gb, 5); | |
1513 | |
1514 if (s->pict_type != I_TYPE) { | |
1515 s->f_code = get_bits(&s->gb, 3); /* fcode_for */ | |
1516 } | |
1517 if (s->shape && (s->pict_type != I_TYPE)) { | |
1518 skip_bits1(&s->gb); // vop shape coding type | |
1519 } | |
1520 } | |
1521 return 0; | |
0 | 1522 } |
1523 | |
1524 /* don't understand why they choose a different header ! */ | |
1525 int intel_h263_decode_picture_header(MpegEncContext *s) | |
1526 { | |
1527 int format; | |
1528 | |
1529 /* picture header */ | |
1530 if (get_bits(&s->gb, 22) != 0x20) | |
1531 return -1; | |
21 | 1532 skip_bits(&s->gb, 8); /* picture timestamp */ |
0 | 1533 |
21 | 1534 if (get_bits1(&s->gb) != 1) |
0 | 1535 return -1; /* marker */ |
21 | 1536 if (get_bits1(&s->gb) != 0) |
0 | 1537 return -1; /* h263 id */ |
21 | 1538 skip_bits1(&s->gb); /* split screen off */ |
1539 skip_bits1(&s->gb); /* camera off */ | |
1540 skip_bits1(&s->gb); /* freeze picture release off */ | |
0 | 1541 |
1542 format = get_bits(&s->gb, 3); | |
1543 if (format != 7) | |
1544 return -1; | |
1545 | |
1546 s->h263_plus = 0; | |
1547 | |
21 | 1548 s->pict_type = I_TYPE + get_bits1(&s->gb); |
0 | 1549 |
21 | 1550 s->unrestricted_mv = get_bits1(&s->gb); |
0 | 1551 s->h263_long_vectors = s->unrestricted_mv; |
1552 | |
21 | 1553 if (get_bits1(&s->gb) != 0) |
0 | 1554 return -1; /* SAC: off */ |
21 | 1555 if (get_bits1(&s->gb) != 0) |
0 | 1556 return -1; /* advanced prediction mode: off */ |
21 | 1557 if (get_bits1(&s->gb) != 0) |
0 | 1558 return -1; /* not PB frame */ |
1559 | |
1560 /* skip unknown header garbage */ | |
21 | 1561 skip_bits(&s->gb, 41); |
0 | 1562 |
1563 s->qscale = get_bits(&s->gb, 5); | |
21 | 1564 skip_bits1(&s->gb); /* Continuous Presence Multipoint mode: off */ |
0 | 1565 |
1566 /* PEI */ | |
21 | 1567 while (get_bits1(&s->gb) != 0) { |
1568 skip_bits(&s->gb, 8); | |
0 | 1569 } |
1570 s->f_code = 1; | |
1571 return 0; | |
1572 } | |
144 | 1573 |