Mercurial > libavcodec.hg
annotate h263.c @ 358:2066dc543be4 libavcodec
move put_string() to common.{c,h}
author | michaelni |
---|---|
date | Fri, 03 May 2002 23:13:24 +0000 |
parents | ac6fc4c7aecb |
children | df0736462b9f |
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. | |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
20 * |
327 | 21 * ac prediction encoding & b-frame support by Michael Niedermayer <michaelni@gmx.at> |
0 | 22 */ |
355 | 23 |
24 //#define DEBUG | |
0 | 25 #include "common.h" |
26 #include "dsputil.h" | |
27 #include "avcodec.h" | |
28 #include "mpegvideo.h" | |
29 #include "h263data.h" | |
30 #include "mpeg4data.h" | |
31 | |
255 | 32 //rounded divison & shift |
33 #define RDIV(a,b) ((a) > 0 ? ((a)+((b)>>1))/(b) : ((a)-((b)>>1))/(b)) | |
34 #define RSHIFT(a,b) ((a) > 0 ? ((a) + (1<<((b)-1)))>>(b) : ((a) + (1<<((b)-1))-1)>>(b)) | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
35 #define ABS(a) (((a)>=0)?(a):(-(a))) |
290 | 36 #define MAX(a,b) ((a) > (b) ? (a) : (b)) |
37 #define MIN(a,b) ((a) < (b) ? (a) : (b)) | |
255 | 38 |
0 | 39 static void h263_encode_block(MpegEncContext * s, DCTELEM * block, |
40 int n); | |
324 | 41 static void h263_encode_motion(MpegEncContext * s, int val, int fcode); |
78 | 42 static void h263p_encode_umotion(MpegEncContext * s, int val); |
0 | 43 static void mpeg4_encode_block(MpegEncContext * s, DCTELEM * block, |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
44 int n, int dc, UINT8 *scan_table); |
262 | 45 static int h263_decode_motion(MpegEncContext * s, int pred, int fcode); |
78 | 46 static int h263p_decode_umotion(MpegEncContext * s, int pred); |
0 | 47 static int h263_decode_block(MpegEncContext * s, DCTELEM * block, |
48 int n, int coded); | |
49 static int mpeg4_decode_block(MpegEncContext * s, DCTELEM * block, | |
50 int n, int coded); | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
51 static int h263_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr); |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
52 static inline int mpeg4_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr, int *dir_ptr); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
53 static void mpeg4_inv_pred_ac(MpegEncContext * s, INT16 *block, int n, |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
54 int dir); |
290 | 55 static void mpeg4_decode_sprite_trajectory(MpegEncContext * s); |
56 | |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
57 extern UINT32 inverse[256]; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
58 |
292
73a9ce3d9715
fcode_tables where too small, found by Klaas-Pieter Vlieg <vlieg@eurescom.de>
michaelni
parents:
291
diff
changeset
|
59 static UINT16 mv_penalty[MAX_FCODE+1][MAX_MV*2+1]; |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
60 static UINT8 fcode_tab[MAX_MV*2+1]; |
287 | 61 static UINT8 umv_fcode_tab[MAX_MV*2+1]; |
0 | 62 |
293 | 63 static UINT16 uni_DCtab_lum [512][2]; |
64 static UINT16 uni_DCtab_chrom[512][2]; | |
65 | |
0 | 66 int h263_get_picture_format(int width, int height) |
67 { | |
68 int format; | |
69 | |
70 if (width == 128 && height == 96) | |
355 | 71 format = 1; |
0 | 72 else if (width == 176 && height == 144) |
355 | 73 format = 2; |
0 | 74 else if (width == 352 && height == 288) |
355 | 75 format = 3; |
0 | 76 else if (width == 704 && height == 576) |
355 | 77 format = 4; |
0 | 78 else if (width == 1408 && height == 1152) |
355 | 79 format = 5; |
0 | 80 else |
81 format = 7; | |
82 return format; | |
83 } | |
84 | |
85 void h263_encode_picture_header(MpegEncContext * s, int picture_number) | |
86 { | |
78 | 87 int format; |
0 | 88 |
89 align_put_bits(&s->pb); | |
231 | 90 |
91 /* 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
|
92 s->ptr_lastgob = pbBufPtr(&s->pb); |
231 | 93 s->gob_number = 0; |
94 | |
95 put_bits(&s->pb, 22, 0x20); /* PSC */ | |
241 | 96 put_bits(&s->pb, 8, (((INT64)s->picture_number * 30 * FRAME_RATE_BASE) / |
0 | 97 s->frame_rate) & 0xff); |
98 | |
99 put_bits(&s->pb, 1, 1); /* marker */ | |
100 put_bits(&s->pb, 1, 0); /* h263 id */ | |
101 put_bits(&s->pb, 1, 0); /* split screen off */ | |
102 put_bits(&s->pb, 1, 0); /* camera off */ | |
103 put_bits(&s->pb, 1, 0); /* freeze picture release off */ | |
78 | 104 |
105 format = h263_get_picture_format(s->width, s->height); | |
0 | 106 if (!s->h263_plus) { |
107 /* H.263v1 */ | |
108 put_bits(&s->pb, 3, format); | |
109 put_bits(&s->pb, 1, (s->pict_type == P_TYPE)); | |
110 /* By now UMV IS DISABLED ON H.263v1, since the restrictions | |
111 of H.263v1 UMV implies to check the predicted MV after | |
112 calculation of the current MB to see if we're on the limits */ | |
113 put_bits(&s->pb, 1, 0); /* unrestricted motion vector: off */ | |
114 put_bits(&s->pb, 1, 0); /* SAC: off */ | |
115 put_bits(&s->pb, 1, 0); /* advanced prediction mode: off */ | |
116 put_bits(&s->pb, 1, 0); /* not PB frame */ | |
117 put_bits(&s->pb, 5, s->qscale); | |
118 put_bits(&s->pb, 1, 0); /* Continuous Presence Multipoint mode: off */ | |
119 } else { | |
120 /* H.263v2 */ | |
121 /* H.263 Plus PTYPE */ | |
122 put_bits(&s->pb, 3, 7); | |
123 put_bits(&s->pb,3,1); /* Update Full Extended PTYPE */ | |
78 | 124 if (format == 7) |
125 put_bits(&s->pb,3,6); /* Custom Source Format */ | |
126 else | |
127 put_bits(&s->pb, 3, format); | |
128 | |
0 | 129 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
|
130 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
|
131 put_bits(&s->pb, 1, s->umvplus); /* Unrestricted Motion Vector */ |
0 | 132 put_bits(&s->pb,1,0); /* SAC: off */ |
133 put_bits(&s->pb,1,0); /* Advanced Prediction Mode: off */ | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
134 put_bits(&s->pb,1,s->h263_aic); /* Advanced Intra Coding */ |
0 | 135 put_bits(&s->pb,1,0); /* Deblocking Filter: off */ |
136 put_bits(&s->pb,1,0); /* Slice Structured: off */ | |
137 put_bits(&s->pb,1,0); /* Reference Picture Selection: off */ | |
138 put_bits(&s->pb,1,0); /* Independent Segment Decoding: off */ | |
139 put_bits(&s->pb,1,0); /* Alternative Inter VLC: off */ | |
140 put_bits(&s->pb,1,0); /* Modified Quantization: off */ | |
141 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */ | |
142 put_bits(&s->pb,3,0); /* Reserved */ | |
143 | |
144 put_bits(&s->pb, 3, s->pict_type == P_TYPE); | |
145 | |
146 put_bits(&s->pb,1,0); /* Reference Picture Resampling: off */ | |
147 put_bits(&s->pb,1,0); /* Reduced-Resolution Update: off */ | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
148 if (s->pict_type == I_TYPE) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
149 s->no_rounding = 0; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
150 else |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
151 s->no_rounding ^= 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
152 put_bits(&s->pb,1,s->no_rounding); /* Rounding Type */ |
0 | 153 put_bits(&s->pb,2,0); /* Reserved */ |
154 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */ | |
155 | |
156 /* This should be here if PLUSPTYPE */ | |
157 put_bits(&s->pb, 1, 0); /* Continuous Presence Multipoint mode: off */ | |
158 | |
78 | 159 if (format == 7) { |
160 /* Custom Picture Format (CPFMT) */ | |
0 | 161 |
355 | 162 if (s->aspect_ratio_info) |
163 put_bits(&s->pb,4,s->aspect_ratio_info); | |
164 else | |
78 | 165 put_bits(&s->pb,4,2); /* Aspect ratio: CIF 12:11 (4:3) picture */ |
166 put_bits(&s->pb,9,(s->width >> 2) - 1); | |
167 put_bits(&s->pb,1,1); /* "1" to prevent start code emulation */ | |
168 put_bits(&s->pb,9,(s->height >> 2)); | |
169 } | |
170 | |
0 | 171 /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */ |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
172 if (s->umvplus) |
0 | 173 put_bits(&s->pb,1,1); /* Limited according tables of Annex D */ |
174 put_bits(&s->pb, 5, s->qscale); | |
175 } | |
176 | |
177 put_bits(&s->pb, 1, 0); /* no PEI */ | |
178 } | |
179 | |
162 | 180 int h263_encode_gob_header(MpegEncContext * s, int mb_line) |
181 { | |
182 int pdif=0; | |
183 | |
184 /* Check to see if we need to put a new GBSC */ | |
185 /* for RTP packetization */ | |
186 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
|
187 pdif = pbBufPtr(&s->pb) - s->ptr_lastgob; |
162 | 188 if (pdif >= s->rtp_payload_size) { |
189 /* Bad luck, packet must be cut before */ | |
190 align_put_bits(&s->pb); | |
231 | 191 flush_put_bits(&s->pb); |
192 /* Call the RTP callback to send the last GOB */ | |
193 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
|
194 pdif = pbBufPtr(&s->pb) - s->ptr_lastgob; |
231 | 195 s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number); |
196 } | |
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
|
197 s->ptr_lastgob = pbBufPtr(&s->pb); |
162 | 198 put_bits(&s->pb, 17, 1); /* GBSC */ |
231 | 199 s->gob_number = mb_line / s->gob_index; |
162 | 200 put_bits(&s->pb, 5, s->gob_number); /* GN */ |
231 | 201 put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */ |
162 | 202 put_bits(&s->pb, 5, s->qscale); /* GQUANT */ |
231 | 203 //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif); |
162 | 204 return pdif; |
205 } else if (pdif + s->mb_line_avgsize >= s->rtp_payload_size) { | |
206 /* Cut the packet before we can't */ | |
207 align_put_bits(&s->pb); | |
231 | 208 flush_put_bits(&s->pb); |
209 /* Call the RTP callback to send the last GOB */ | |
210 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
|
211 pdif = pbBufPtr(&s->pb) - s->ptr_lastgob; |
231 | 212 s->rtp_callback(s->ptr_lastgob, pdif, s->gob_number); |
213 } | |
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
|
214 s->ptr_lastgob = pbBufPtr(&s->pb); |
162 | 215 put_bits(&s->pb, 17, 1); /* GBSC */ |
231 | 216 s->gob_number = mb_line / s->gob_index; |
162 | 217 put_bits(&s->pb, 5, s->gob_number); /* GN */ |
231 | 218 put_bits(&s->pb, 2, s->pict_type == I_TYPE); /* GFID */ |
162 | 219 put_bits(&s->pb, 5, s->qscale); /* GQUANT */ |
231 | 220 //fprintf(stderr,"\nGOB: %2d size: %d", s->gob_number - 1, pdif); |
162 | 221 return pdif; |
222 } | |
223 } | |
224 return 0; | |
225 } | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
226 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
227 static inline int decide_ac_pred(MpegEncContext * s, DCTELEM block[6][64], int dir[6]) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
228 { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
229 int score0=0, score1=0; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
230 int i, n; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
231 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
232 for(n=0; n<6; n++){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
233 INT16 *ac_val, *ac_val1; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
234 |
266 | 235 ac_val = s->ac_val[0][0] + s->block_index[n] * 16; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
236 ac_val1= ac_val; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
237 if(dir[n]){ |
266 | 238 ac_val-= s->block_wrap[n]*16; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
239 for(i=1; i<8; i++){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
240 const int level= block[n][block_permute_op(i )]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
241 score0+= ABS(level); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
242 score1+= ABS(level - ac_val[i+8]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
243 ac_val1[i ]= block[n][block_permute_op(i<<3)]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
244 ac_val1[i+8]= level; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
245 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
246 }else{ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
247 ac_val-= 16; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
248 for(i=1; i<8; i++){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
249 const int level= block[n][block_permute_op(i<<3)]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
250 score0+= ABS(level); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
251 score1+= ABS(level - ac_val[i]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
252 ac_val1[i ]= level; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
253 ac_val1[i+8]= block[n][block_permute_op(i )]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
254 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
255 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
256 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
257 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
258 return score0 > score1 ? 1 : 0; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
259 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
260 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
261 void mpeg4_encode_mb(MpegEncContext * s, |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
262 DCTELEM block[6][64], |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
263 int motion_x, int motion_y) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
264 { |
324 | 265 int cbpc, cbpy, i, pred_x, pred_y; |
286 | 266 int bits; |
162 | 267 |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
268 // printf("**mb x=%d y=%d\n", s->mb_x, s->mb_y); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
269 if (!s->mb_intra) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
270 /* compute cbp */ |
324 | 271 int cbp = 0; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
272 for (i = 0; i < 6; i++) { |
324 | 273 if (s->block_last_index[i] >= 0) |
274 cbp |= 1 << (5 - i); | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
275 } |
324 | 276 |
277 if(s->pict_type==B_TYPE){ | |
278 static const int mb_type_table[8]= {-1, 2, 3, 1,-1,-1,-1, 0}; /* convert from mv_dir to type */ | |
279 int mb_type= mb_type_table[s->mv_dir]; | |
280 | |
281 if(s->mb_x==0){ | |
282 s->last_mv[0][0][0]= | |
283 s->last_mv[0][0][1]= | |
284 s->last_mv[1][0][0]= | |
285 s->last_mv[1][0][1]= 0; | |
286 } | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
287 |
324 | 288 /* nothing to do if this MB was skiped in the next P Frame */ |
289 if(s->mbskip_table[s->mb_y * s->mb_width + s->mb_x]){ | |
290 s->skip_count++; | |
291 s->mv[0][0][0]= | |
292 s->mv[0][0][1]= | |
293 s->mv[1][0][0]= | |
294 s->mv[1][0][1]= 0; | |
327 | 295 s->mv_dir= MV_DIR_FORWARD; //doesnt matter |
324 | 296 return; |
297 } | |
298 | |
299 if ((cbp | motion_x | motion_y | mb_type) ==0) { | |
300 /* direct MB with MV={0,0} */ | |
301 put_bits(&s->pb, 1, 1); /* mb not coded modb1=1 */ | |
302 s->misc_bits++; | |
303 s->last_bits++; | |
304 s->skip_count++; | |
305 return; | |
306 } | |
307 put_bits(&s->pb, 1, 0); /* mb coded modb1=0 */ | |
308 put_bits(&s->pb, 1, cbp ? 0 : 1); /* modb2 */ //FIXME merge | |
309 put_bits(&s->pb, mb_type+1, 1); // this table is so simple that we dont need it :) | |
310 if(cbp) put_bits(&s->pb, 6, cbp); | |
311 | |
312 if(cbp && mb_type) | |
313 put_bits(&s->pb, 1, 0); /* no q-scale change */ | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
314 |
295 | 315 bits= get_bit_count(&s->pb); |
316 s->misc_bits+= bits - s->last_bits; | |
317 s->last_bits=bits; | |
318 | |
324 | 319 switch(mb_type) |
320 { | |
321 case 0: /* direct */ | |
322 h263_encode_motion(s, motion_x, 1); | |
323 h263_encode_motion(s, motion_y, 1); | |
324 break; | |
325 case 1: /* bidir */ | |
326 h263_encode_motion(s, s->mv[0][0][0] - s->last_mv[0][0][0], s->f_code); | |
327 h263_encode_motion(s, s->mv[0][0][1] - s->last_mv[0][0][1], s->f_code); | |
328 h263_encode_motion(s, s->mv[1][0][0] - s->last_mv[1][0][0], s->b_code); | |
329 h263_encode_motion(s, s->mv[1][0][1] - s->last_mv[1][0][1], s->b_code); | |
330 s->last_mv[0][0][0]= s->mv[0][0][0]; | |
331 s->last_mv[0][0][1]= s->mv[0][0][1]; | |
332 s->last_mv[1][0][0]= s->mv[1][0][0]; | |
333 s->last_mv[1][0][1]= s->mv[1][0][1]; | |
334 break; | |
335 case 2: /* backward */ | |
336 h263_encode_motion(s, motion_x - s->last_mv[1][0][0], s->b_code); | |
337 h263_encode_motion(s, motion_y - s->last_mv[1][0][1], s->b_code); | |
338 s->last_mv[1][0][0]= motion_x; | |
339 s->last_mv[1][0][1]= motion_y; | |
340 break; | |
341 case 3: /* forward */ | |
342 h263_encode_motion(s, motion_x - s->last_mv[0][0][0], s->f_code); | |
343 h263_encode_motion(s, motion_y - s->last_mv[0][0][1], s->f_code); | |
344 s->last_mv[0][0][0]= motion_x; | |
345 s->last_mv[0][0][1]= motion_y; | |
346 break; | |
327 | 347 default: |
348 printf("unknown mb type\n"); | |
324 | 349 return; |
350 } | |
351 bits= get_bit_count(&s->pb); | |
352 s->mv_bits+= bits - s->last_bits; | |
353 s->last_bits=bits; | |
295 | 354 |
324 | 355 /* encode each block */ |
356 for (i = 0; i < 6; i++) { | |
357 mpeg4_encode_block(s, block[i], i, 0, zigzag_direct); | |
358 } | |
359 bits= get_bit_count(&s->pb); | |
360 s->p_tex_bits+= bits - s->last_bits; | |
361 s->last_bits=bits; | |
362 }else{ /* s->pict_type==B_TYPE */ | |
363 if ((cbp | motion_x | motion_y) == 0 && s->mv_type==MV_TYPE_16X16) { | |
331 | 364 /* check if the B frames can skip it too, as we must skip it if we skip here |
365 why didnt they just compress the skip-mb bits instead of reusing them ?! */ | |
366 if(s->max_b_frames>0){ | |
367 int i; | |
368 const int offset= (s->mb_x + s->mb_y*s->linesize)*16; | |
369 uint8_t *p_pic= s->new_picture[0] + offset; | |
370 s->mb_skiped=1; | |
371 for(i=0; i<s->max_b_frames; i++){ | |
339 | 372 uint8_t *b_pic; |
373 int diff; | |
374 | |
375 if(s->coded_order[i+1].pict_type!=B_TYPE) break; | |
376 | |
377 b_pic= s->coded_order[i+1].picture[0] + offset; | |
378 diff= pix_abs16x16(p_pic, b_pic, s->linesize); | |
331 | 379 if(diff>s->qscale*70){ |
380 s->mb_skiped=0; | |
381 break; | |
382 } | |
383 } | |
384 }else | |
385 s->mb_skiped=1; | |
386 | |
387 if(s->mb_skiped==1){ | |
388 /* skip macroblock */ | |
389 put_bits(&s->pb, 1, 1); | |
390 s->misc_bits++; | |
391 s->last_bits++; | |
392 s->skip_count++; | |
393 return; | |
394 } | |
295 | 395 } |
324 | 396 put_bits(&s->pb, 1, 0); /* mb coded */ |
397 if(s->mv_type==MV_TYPE_16X16){ | |
398 cbpc = cbp & 3; | |
399 put_bits(&s->pb, | |
400 inter_MCBPC_bits[cbpc], | |
401 inter_MCBPC_code[cbpc]); | |
402 cbpy = cbp >> 2; | |
403 cbpy ^= 0xf; | |
404 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); | |
405 | |
406 bits= get_bit_count(&s->pb); | |
407 s->misc_bits+= bits - s->last_bits; | |
408 s->last_bits=bits; | |
409 | |
410 /* motion vectors: 16x16 mode */ | |
411 h263_pred_motion(s, 0, &pred_x, &pred_y); | |
412 | |
413 h263_encode_motion(s, motion_x - pred_x, s->f_code); | |
414 h263_encode_motion(s, motion_y - pred_y, s->f_code); | |
415 }else{ | |
416 cbpc = (cbp & 3)+16; | |
417 put_bits(&s->pb, | |
418 inter_MCBPC_bits[cbpc], | |
419 inter_MCBPC_code[cbpc]); | |
420 cbpy = cbp >> 2; | |
421 cbpy ^= 0xf; | |
422 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); | |
423 | |
424 bits= get_bit_count(&s->pb); | |
425 s->misc_bits+= bits - s->last_bits; | |
426 s->last_bits=bits; | |
427 | |
428 for(i=0; i<4; i++){ | |
429 /* motion vectors: 8x8 mode*/ | |
430 h263_pred_motion(s, i, &pred_x, &pred_y); | |
431 | |
432 h263_encode_motion(s, s->motion_val[ s->block_index[i] ][0] - pred_x, s->f_code); | |
433 h263_encode_motion(s, s->motion_val[ s->block_index[i] ][1] - pred_y, s->f_code); | |
434 } | |
435 } | |
436 bits= get_bit_count(&s->pb); | |
437 s->mv_bits+= bits - s->last_bits; | |
438 s->last_bits=bits; | |
439 | |
440 /* encode each block */ | |
441 for (i = 0; i < 6; i++) { | |
442 mpeg4_encode_block(s, block[i], i, 0, zigzag_direct); | |
443 } | |
444 bits= get_bit_count(&s->pb); | |
445 s->p_tex_bits+= bits - s->last_bits; | |
446 s->last_bits=bits; | |
447 s->p_count++; | |
295 | 448 } |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
449 } else { |
324 | 450 int cbp; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
451 int dc_diff[6]; //dc values with the dc prediction subtracted |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
452 int dir[6]; //prediction direction |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
453 int zigzag_last_index[6]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
454 UINT8 *scan_table[6]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
455 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
456 for(i=0; i<6; i++){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
457 const int level= block[i][0]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
458 UINT16 *dc_ptr; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
459 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
460 dc_diff[i]= level - mpeg4_pred_dc(s, i, &dc_ptr, &dir[i]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
461 if (i < 4) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
462 *dc_ptr = level * s->y_dc_scale; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
463 } else { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
464 *dc_ptr = level * s->c_dc_scale; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
465 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
466 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
467 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
468 s->ac_pred= decide_ac_pred(s, block, dir); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
469 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
470 if(s->ac_pred){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
471 for(i=0; i<6; i++){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
472 UINT8 *st; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
473 int last_index; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
474 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
475 mpeg4_inv_pred_ac(s, block[i], i, dir[i]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
476 if (dir[i]==0) st = ff_alternate_vertical_scan; /* left */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
477 else st = ff_alternate_horizontal_scan; /* top */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
478 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
479 for(last_index=63; last_index>=0; last_index--) //FIXME optimize |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
480 if(block[i][st[last_index]]) break; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
481 zigzag_last_index[i]= s->block_last_index[i]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
482 s->block_last_index[i]= last_index; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
483 scan_table[i]= st; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
484 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
485 }else{ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
486 for(i=0; i<6; i++) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
487 scan_table[i]= zigzag_direct; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
488 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
489 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
490 /* compute cbp */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
491 cbp = 0; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
492 for (i = 0; i < 6; i++) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
493 if (s->block_last_index[i] >= 1) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
494 cbp |= 1 << (5 - i); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
495 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
496 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
497 cbpc = cbp & 3; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
498 if (s->pict_type == I_TYPE) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
499 put_bits(&s->pb, |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
500 intra_MCBPC_bits[cbpc], |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
501 intra_MCBPC_code[cbpc]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
502 } else { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
503 put_bits(&s->pb, 1, 0); /* mb coded */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
504 put_bits(&s->pb, |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
505 inter_MCBPC_bits[cbpc + 4], |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
506 inter_MCBPC_code[cbpc + 4]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
507 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
508 put_bits(&s->pb, 1, s->ac_pred); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
509 cbpy = cbp >> 2; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
510 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
511 |
286 | 512 bits= get_bit_count(&s->pb); |
513 s->misc_bits+= bits - s->last_bits; | |
514 s->last_bits=bits; | |
515 | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
516 /* encode each block */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
517 for (i = 0; i < 6; i++) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
518 mpeg4_encode_block(s, block[i], i, dc_diff[i], scan_table[i]); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
519 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
520 |
286 | 521 bits= get_bit_count(&s->pb); |
522 s->i_tex_bits+= bits - s->last_bits; | |
523 s->last_bits=bits; | |
524 s->i_count++; | |
525 | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
526 /* restore ac coeffs & last_index stuff if we messed them up with the prediction */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
527 if(s->ac_pred){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
528 for(i=0; i<6; i++){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
529 int j; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
530 INT16 *ac_val; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
531 |
266 | 532 ac_val = s->ac_val[0][0] + s->block_index[i] * 16; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
533 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
534 if(dir[i]){ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
535 for(j=1; j<8; j++) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
536 block[i][block_permute_op(j )]= ac_val[j+8]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
537 }else{ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
538 for(j=1; j<8; j++) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
539 block[i][block_permute_op(j<<3)]= ac_val[j ]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
540 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
541 s->block_last_index[i]= zigzag_last_index[i]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
542 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
543 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
544 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
545 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
546 |
0 | 547 void h263_encode_mb(MpegEncContext * s, |
548 DCTELEM block[6][64], | |
549 int motion_x, int motion_y) | |
550 { | |
551 int cbpc, cbpy, i, cbp, pred_x, pred_y; | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
552 INT16 pred_dc; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
553 INT16 rec_intradc[6]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
554 UINT16 *dc_ptr[6]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
555 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
556 //printf("**mb x=%d y=%d\n", s->mb_x, s->mb_y); |
324 | 557 if (!s->mb_intra) { |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
558 /* compute cbp */ |
324 | 559 cbp = 0; |
560 for (i = 0; i < 6; i++) { | |
561 if (s->block_last_index[i] >= 0) | |
562 cbp |= 1 << (5 - i); | |
563 } | |
564 if ((cbp | motion_x | motion_y) == 0) { | |
565 /* skip macroblock */ | |
566 put_bits(&s->pb, 1, 1); | |
567 return; | |
568 } | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
569 put_bits(&s->pb, 1, 0); /* mb coded */ |
324 | 570 cbpc = cbp & 3; |
571 put_bits(&s->pb, | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
572 inter_MCBPC_bits[cbpc], |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
573 inter_MCBPC_code[cbpc]); |
324 | 574 cbpy = cbp >> 2; |
575 cbpy ^= 0xf; | |
576 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); | |
0 | 577 |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
578 /* motion vectors: 16x16 mode only now */ |
324 | 579 h263_pred_motion(s, 0, &pred_x, &pred_y); |
78 | 580 |
324 | 581 if (!s->umvplus) { |
582 h263_encode_motion(s, motion_x - pred_x, s->f_code); | |
583 h263_encode_motion(s, motion_y - pred_y, s->f_code); | |
584 } | |
585 else { | |
586 h263p_encode_umotion(s, motion_x - pred_x); | |
587 h263p_encode_umotion(s, motion_y - pred_y); | |
588 if (((motion_x - pred_x) == 1) && ((motion_y - pred_y) == 1)) | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
589 /* To prevent Start Code emulation */ |
324 | 590 put_bits(&s->pb,1,1); |
591 } | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
592 } else { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
593 int li = s->h263_aic ? 0 : 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
594 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
595 cbp = 0; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
596 for(i=0; i<6; i++) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
597 /* Predict DC */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
598 if (s->h263_aic && s->mb_intra) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
599 INT16 level = block[i][0]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
600 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
601 pred_dc = h263_pred_dc(s, i, &dc_ptr[i]); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
602 level -= pred_dc; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
603 /* Quant */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
604 if (level < 0) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
605 level = (level + (s->qscale >> 1))/(s->y_dc_scale); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
606 else |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
607 level = (level - (s->qscale >> 1))/(s->y_dc_scale); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
608 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
609 /* AIC can change CBP */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
610 if (level == 0 && s->block_last_index[i] == 0) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
611 s->block_last_index[i] = -1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
612 else if (level < -127) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
613 level = -127; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
614 else if (level > 127) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
615 level = 127; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
616 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
617 block[i][0] = level; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
618 /* Reconstruction */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
619 rec_intradc[i] = (s->y_dc_scale*level) + pred_dc; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
620 /* Oddify */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
621 rec_intradc[i] |= 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
622 //if ((rec_intradc[i] % 2) == 0) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
623 // rec_intradc[i]++; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
624 /* Clipping */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
625 if (rec_intradc[i] < 0) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
626 rec_intradc[i] = 0; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
627 else if (rec_intradc[i] > 2047) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
628 rec_intradc[i] = 2047; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
629 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
630 /* Update AC/DC tables */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
631 *dc_ptr[i] = rec_intradc[i]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
632 } |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
633 /* compute cbp */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
634 if (s->block_last_index[i] >= li) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
635 cbp |= 1 << (5 - i); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
636 } |
0 | 637 |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
638 cbpc = cbp & 3; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
639 if (s->pict_type == I_TYPE) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
640 put_bits(&s->pb, |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
641 intra_MCBPC_bits[cbpc], |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
642 intra_MCBPC_code[cbpc]); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
643 } else { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
644 put_bits(&s->pb, 1, 0); /* mb coded */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
645 put_bits(&s->pb, |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
646 inter_MCBPC_bits[cbpc + 4], |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
647 inter_MCBPC_code[cbpc + 4]); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
648 } |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
649 if (s->h263_aic) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
650 /* XXX: currently, we do not try to use ac prediction */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
651 put_bits(&s->pb, 1, 0); /* no AC prediction */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
652 } |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
653 cbpy = cbp >> 2; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
654 put_bits(&s->pb, cbpy_tab[cbpy][1], cbpy_tab[cbpy][0]); |
0 | 655 } |
656 | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
657 for(i=0; i<6; i++) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
658 /* encode each block */ |
294 | 659 h263_encode_block(s, block[i], i); |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
660 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
661 /* Update INTRADC for decoding */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
662 if (s->h263_aic && s->mb_intra) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
663 block[i][0] = rec_intradc[i]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
664 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
665 } |
0 | 666 } |
667 } | |
668 | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
669 static int h263_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
670 { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
671 int x, y, wrap, a, c, pred_dc, scale; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
672 INT16 *dc_val, *ac_val; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
673 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
674 /* find prediction */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
675 if (n < 4) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
676 x = 2 * s->mb_x + 1 + (n & 1); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
677 y = 2 * s->mb_y + 1 + ((n & 2) >> 1); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
678 wrap = s->mb_width * 2 + 2; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
679 dc_val = s->dc_val[0]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
680 ac_val = s->ac_val[0][0]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
681 scale = s->y_dc_scale; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
682 } else { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
683 x = s->mb_x + 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
684 y = s->mb_y + 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
685 wrap = s->mb_width + 2; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
686 dc_val = s->dc_val[n - 4 + 1]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
687 ac_val = s->ac_val[n - 4 + 1][0]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
688 scale = s->c_dc_scale; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
689 } |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
690 /* B C |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
691 * A X |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
692 */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
693 a = dc_val[(x - 1) + (y) * wrap]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
694 c = dc_val[(x) + (y - 1) * wrap]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
695 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
696 /* No prediction outside GOB boundary */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
697 if (s->first_gob_line && ((n < 2) || (n > 3))) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
698 c = 1024; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
699 pred_dc = 1024; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
700 /* just DC prediction */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
701 if (a != 1024 && c != 1024) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
702 pred_dc = (a + c) >> 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
703 else if (a != 1024) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
704 pred_dc = a; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
705 else |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
706 pred_dc = c; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
707 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
708 /* we assume pred is positive */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
709 //pred_dc = (pred_dc + (scale >> 1)) / scale; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
710 *dc_val_ptr = &dc_val[x + y * wrap]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
711 return pred_dc; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
712 } |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
713 |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
714 |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
715 void h263_pred_acdc(MpegEncContext * s, INT16 *block, int n) |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
716 { |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
717 int x, y, wrap, a, c, pred_dc, scale, i; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
718 INT16 *dc_val, *ac_val, *ac_val1; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
719 |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
720 /* find prediction */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
721 if (n < 4) { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
722 x = 2 * s->mb_x + 1 + (n & 1); |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
723 y = 2 * s->mb_y + 1 + ((n & 2) >> 1); |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
724 wrap = s->mb_width * 2 + 2; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
725 dc_val = s->dc_val[0]; |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
726 ac_val = s->ac_val[0][0]; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
727 scale = s->y_dc_scale; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
728 } else { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
729 x = s->mb_x + 1; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
730 y = s->mb_y + 1; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
731 wrap = s->mb_width + 2; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
732 dc_val = s->dc_val[n - 4 + 1]; |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
733 ac_val = s->ac_val[n - 4 + 1][0]; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
734 scale = s->c_dc_scale; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
735 } |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
736 |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
737 ac_val += ((y) * wrap + (x)) * 16; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
738 ac_val1 = ac_val; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
739 |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
740 /* B C |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
741 * A X |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
742 */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
743 a = dc_val[(x - 1) + (y) * wrap]; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
744 c = dc_val[(x) + (y - 1) * wrap]; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
745 |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
746 /* No prediction outside GOB boundary */ |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
747 if (s->first_gob_line && ((n < 2) || (n > 3))) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
748 c = 1024; |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
749 pred_dc = 1024; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
750 if (s->ac_pred) { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
751 if (s->h263_aic_dir) { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
752 /* left prediction */ |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
753 if (a != 1024) { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
754 ac_val -= 16; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
755 for(i=1;i<8;i++) { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
756 block[block_permute_op(i*8)] += ac_val[i]; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
757 } |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
758 pred_dc = a; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
759 } |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
760 } else { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
761 /* top prediction */ |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
762 if (c != 1024) { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
763 ac_val -= 16 * wrap; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
764 for(i=1;i<8;i++) { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
765 block[block_permute_op(i)] += ac_val[i + 8]; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
766 } |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
767 pred_dc = c; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
768 } |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
769 } |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
770 } else { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
771 /* just DC prediction */ |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
772 if (a != 1024 && c != 1024) |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
773 pred_dc = (a + c) >> 1; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
774 else if (a != 1024) |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
775 pred_dc = a; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
776 else |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
777 pred_dc = c; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
778 } |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
779 |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
780 /* we assume pred is positive */ |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
781 block[0]=block[0]*scale + pred_dc; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
782 |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
783 if (block[0] < 0) |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
784 block[0] = 0; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
785 else if (!(block[0] & 1)) |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
786 block[0]++; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
787 |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
788 /* Update AC/DC tables */ |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
789 dc_val[(x) + (y) * wrap] = block[0]; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
790 |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
791 /* left copy */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
792 for(i=1;i<8;i++) |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
793 ac_val1[i] = block[block_permute_op(i * 8)]; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
794 /* top copy */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
795 for(i=1;i<8;i++) |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
796 ac_val1[8 + i] = block[block_permute_op(i)]; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
797 } |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
798 |
0 | 799 INT16 *h263_pred_motion(MpegEncContext * s, int block, |
800 int *px, int *py) | |
801 { | |
266 | 802 int xy, wrap; |
0 | 803 INT16 *A, *B, *C, *mot_val; |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
804 static const int off[4]= {2, 1, 1, -1}; |
0 | 805 |
266 | 806 wrap = s->block_wrap[0]; |
807 xy = s->block_index[block]; | |
0 | 808 |
245 | 809 mot_val = s->motion_val[xy]; |
0 | 810 |
811 /* special case for first line */ | |
272 | 812 if ((s->mb_y == 0 || s->first_slice_line || s->first_gob_line) && block<2) { |
245 | 813 A = s->motion_val[xy - 1]; |
0 | 814 *px = A[0]; |
815 *py = A[1]; | |
816 } else { | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
817 A = s->motion_val[xy - 1]; |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
818 B = s->motion_val[xy - wrap]; |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
819 C = s->motion_val[xy + off[block] - wrap]; |
0 | 820 *px = mid_pred(A[0], B[0], C[0]); |
821 *py = mid_pred(A[1], B[1], C[1]); | |
822 } | |
823 return mot_val; | |
824 } | |
825 | |
324 | 826 static void h263_encode_motion(MpegEncContext * s, int val, int f_code) |
0 | 827 { |
828 int range, l, m, bit_size, sign, code, bits; | |
829 | |
830 if (val == 0) { | |
831 /* zero vector */ | |
832 code = 0; | |
833 put_bits(&s->pb, mvtab[code][1], mvtab[code][0]); | |
834 } else { | |
324 | 835 bit_size = f_code - 1; |
0 | 836 range = 1 << bit_size; |
837 /* modulo encoding */ | |
838 l = range * 32; | |
839 m = 2 * l; | |
840 if (val < -l) { | |
841 val += m; | |
842 } else if (val >= l) { | |
843 val -= m; | |
844 } | |
845 | |
846 if (val >= 0) { | |
847 sign = 0; | |
848 } else { | |
849 val = -val; | |
850 sign = 1; | |
851 } | |
312 | 852 val--; |
853 code = (val >> bit_size) + 1; | |
854 bits = val & (range - 1); | |
0 | 855 |
856 put_bits(&s->pb, mvtab[code][1] + 1, (mvtab[code][0] << 1) | sign); | |
857 if (bit_size > 0) { | |
858 put_bits(&s->pb, bit_size, bits); | |
859 } | |
860 } | |
861 } | |
862 | |
78 | 863 /* Encode MV differences on H.263+ with Unrestricted MV mode */ |
864 static void h263p_encode_umotion(MpegEncContext * s, int val) | |
865 { | |
866 short sval = 0; | |
867 short i = 0; | |
868 short n_bits = 0; | |
869 short temp_val; | |
870 int code = 0; | |
871 int tcode; | |
872 | |
873 if ( val == 0) | |
874 put_bits(&s->pb, 1, 1); | |
875 else if (val == 1) | |
876 put_bits(&s->pb, 3, 0); | |
877 else if (val == -1) | |
878 put_bits(&s->pb, 3, 2); | |
879 else { | |
880 | |
881 sval = ((val < 0) ? (short)(-val):(short)val); | |
882 temp_val = sval; | |
883 | |
884 while (temp_val != 0) { | |
885 temp_val = temp_val >> 1; | |
886 n_bits++; | |
887 } | |
888 | |
889 i = n_bits - 1; | |
890 while (i > 0) { | |
891 tcode = (sval & (1 << (i-1))) >> (i-1); | |
892 tcode = (tcode << 1) | 1; | |
893 code = (code << 2) | tcode; | |
894 i--; | |
895 } | |
896 code = ((code << 1) | (val < 0)) << 1; | |
897 put_bits(&s->pb, (2*n_bits)+1, code); | |
898 //printf("\nVal = %d\tCode = %d", sval, code); | |
899 } | |
900 } | |
901 | |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
902 static void init_mv_penalty_and_fcode(MpegEncContext *s) |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
903 { |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
904 int f_code; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
905 int mv; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
906 for(f_code=1; f_code<=MAX_FCODE; f_code++){ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
907 for(mv=-MAX_MV; mv<=MAX_MV; mv++){ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
908 int len; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
909 |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
910 if(mv==0) len= mvtab[0][1]; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
911 else{ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
912 int val, bit_size, range, code; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
913 |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
914 bit_size = s->f_code - 1; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
915 range = 1 << bit_size; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
916 |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
917 val=mv; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
918 if (val < 0) |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
919 val = -val; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
920 val--; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
921 code = (val >> bit_size) + 1; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
922 if(code<33){ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
923 len= mvtab[code][1] + 1 + bit_size; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
924 }else{ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
925 len= mvtab[32][1] + 2 + bit_size; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
926 } |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
927 } |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
928 |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
929 mv_penalty[f_code][mv+MAX_MV]= len; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
930 } |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
931 } |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
932 |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
933 for(f_code=MAX_FCODE; f_code>0; f_code--){ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
934 for(mv=-(16<<f_code); mv<(16<<f_code); mv++){ |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
935 fcode_tab[mv+MAX_MV]= f_code; |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
936 } |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
937 } |
287 | 938 |
939 for(mv=0; mv<MAX_MV*2+1; mv++){ | |
940 umv_fcode_tab[mv]= 1; | |
941 } | |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
942 } |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
943 |
293 | 944 static void init_uni_dc_tab() |
945 { | |
946 int level, uni_code, uni_len; | |
947 | |
312 | 948 for(level=-256; level<256; level++){ |
293 | 949 int size, v, l; |
950 /* find number of bits */ | |
951 size = 0; | |
952 v = abs(level); | |
953 while (v) { | |
954 v >>= 1; | |
955 size++; | |
956 } | |
957 | |
958 if (level < 0) | |
959 l= (-level) ^ ((1 << size) - 1); | |
960 else | |
961 l= level; | |
962 | |
963 /* luminance */ | |
964 uni_code= DCtab_lum[size][0]; | |
965 uni_len = DCtab_lum[size][1]; | |
966 | |
967 if (size > 0) { | |
968 uni_code<<=size; uni_code|=l; | |
969 uni_len+=size; | |
970 if (size > 8){ | |
971 uni_code<<=1; uni_code|=1; | |
972 uni_len++; | |
973 } | |
974 } | |
975 uni_DCtab_lum[level+256][0]= uni_code; | |
976 uni_DCtab_lum[level+256][1]= uni_len; | |
977 | |
978 /* chrominance */ | |
979 uni_code= DCtab_chrom[size][0]; | |
980 uni_len = DCtab_chrom[size][1]; | |
981 | |
982 if (size > 0) { | |
983 uni_code<<=size; uni_code|=l; | |
984 uni_len+=size; | |
985 if (size > 8){ | |
986 uni_code<<=1; uni_code|=1; | |
987 uni_len++; | |
988 } | |
989 } | |
990 uni_DCtab_chrom[level+256][0]= uni_code; | |
991 uni_DCtab_chrom[level+256][1]= uni_len; | |
992 | |
993 } | |
994 } | |
995 | |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
996 void h263_encode_init(MpegEncContext *s) |
0 | 997 { |
998 static int done = 0; | |
999 | |
1000 if (!done) { | |
1001 done = 1; | |
293 | 1002 |
1003 init_uni_dc_tab(); | |
1004 | |
0 | 1005 init_rl(&rl_inter); |
1006 init_rl(&rl_intra); | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1007 init_rl(&rl_intra_aic); |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
1008 |
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
1009 init_mv_penalty_and_fcode(s); |
0 | 1010 } |
287 | 1011 s->mv_penalty= mv_penalty; //FIXME exact table for msmpeg4 & h263p |
277
5cb2978e701f
new motion estimation (epzs) not complete yet but allready pretty good :)
michaelni
parents:
274
diff
changeset
|
1012 |
287 | 1013 // use fcodes >1 only for mpeg4 & h263 & h263p FIXME |
344 | 1014 switch(s->codec_id){ |
1015 case CODEC_ID_MPEG4: | |
1016 s->fcode_tab= fcode_tab; | |
1017 s->min_qcoeff= -2048; | |
1018 s->max_qcoeff= 2047; | |
1019 break; | |
1020 case CODEC_ID_H263P: | |
1021 s->fcode_tab= umv_fcode_tab; | |
1022 s->min_qcoeff= -128; | |
1023 s->max_qcoeff= 127; | |
1024 break; | |
1025 default: //nothing needed default table allready set in mpegvideo.c | |
1026 s->min_qcoeff= -128; | |
1027 s->max_qcoeff= 127; | |
1028 } | |
1029 | |
1030 /* h263 type bias */ | |
1031 //FIXME mpeg4 mpeg quantizer | |
1032 s->intra_quant_bias=0; | |
1033 s->inter_quant_bias=-(1<<(QUANT_BIAS_SHIFT-2)); //(a - x/4)/x | |
0 | 1034 } |
1035 | |
1036 static void h263_encode_block(MpegEncContext * s, DCTELEM * block, int n) | |
1037 { | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1038 int level, run, last, i, j, last_index, last_non_zero, sign, slevel, code; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1039 RLTable *rl; |
0 | 1040 |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1041 rl = &rl_inter; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1042 if (s->mb_intra && !s->h263_aic) { |
231 | 1043 /* DC coef */ |
1044 level = block[0]; | |
0 | 1045 /* 255 cannot be represented, so we clamp */ |
1046 if (level > 254) { | |
1047 level = 254; | |
1048 block[0] = 254; | |
1049 } | |
231 | 1050 /* 0 cannot be represented also */ |
1051 else if (!level) { | |
1052 level = 1; | |
1053 block[0] = 1; | |
1054 } | |
1055 if (level == 128) | |
1056 put_bits(&s->pb, 8, 0xff); | |
1057 else | |
1058 put_bits(&s->pb, 8, level & 0xff); | |
1059 i = 1; | |
0 | 1060 } else { |
231 | 1061 i = 0; |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1062 if (s->h263_aic && s->mb_intra) |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1063 rl = &rl_intra_aic; |
0 | 1064 } |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1065 |
0 | 1066 /* AC coefs */ |
1067 last_index = s->block_last_index[n]; | |
1068 last_non_zero = i - 1; | |
1069 for (; i <= last_index; i++) { | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1070 j = zigzag_direct[i]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1071 level = block[j]; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1072 if (level) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1073 run = i - last_non_zero - 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1074 last = (i == last_index); |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1075 sign = 0; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1076 slevel = level; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1077 if (level < 0) { |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1078 sign = 1; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1079 level = -level; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1080 } |
0 | 1081 code = get_rl_index(rl, last, run, level); |
1082 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
1083 if (code == rl->n) { | |
1084 put_bits(&s->pb, 1, last); | |
1085 put_bits(&s->pb, 6, run); | |
1086 put_bits(&s->pb, 8, slevel & 0xff); | |
1087 } else { | |
1088 put_bits(&s->pb, 1, sign); | |
1089 } | |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1090 last_non_zero = i; |
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
1091 } |
0 | 1092 } |
1093 } | |
1094 | |
1095 /***************************************************/ | |
1096 | |
262 | 1097 static void mpeg4_stuffing(PutBitContext * pbc) |
1098 { | |
1099 int length; | |
1100 put_bits(pbc, 1, 0); | |
1101 length= (-get_bit_count(pbc))&7; | |
1102 put_bits(pbc, length, (1<<length)-1); | |
1103 } | |
1104 | |
327 | 1105 /* must be called before writing the header */ |
1106 void ff_set_mpeg4_time(MpegEncContext * s, int picture_number){ | |
1107 int time_div, time_mod; | |
1108 | |
1109 if(s->pict_type==I_TYPE){ //we will encode a vol header | |
1110 s->time_increment_resolution= s->frame_rate/ff_gcd(s->frame_rate, FRAME_RATE_BASE); | |
1111 if(s->time_increment_resolution>=256*256) s->time_increment_resolution= 256*128; | |
1112 | |
1113 s->time_increment_bits = av_log2(s->time_increment_resolution - 1) + 1; | |
1114 } | |
1115 | |
1116 s->time= picture_number*(int64_t)FRAME_RATE_BASE*s->time_increment_resolution/s->frame_rate; | |
1117 time_div= s->time/s->time_increment_resolution; | |
1118 time_mod= s->time%s->time_increment_resolution; | |
1119 | |
1120 if(s->pict_type==B_TYPE){ | |
1121 s->bp_time= s->last_non_b_time - s->time; | |
1122 }else{ | |
1123 s->last_time_base= s->time_base; | |
1124 s->time_base= time_div; | |
1125 s->pp_time= s->time - s->last_non_b_time; | |
1126 s->last_non_b_time= s->time; | |
1127 } | |
1128 } | |
1129 | |
263 | 1130 static void mpeg4_encode_vol_header(MpegEncContext * s) |
0 | 1131 { |
263 | 1132 int vo_ver_id=1; //must be 2 if we want GMC or q-pel |
322 | 1133 char buf[255]; |
336 | 1134 |
1135 s->vo_type= s->has_b_frames ? CORE_VO_TYPE : SIMPLE_VO_TYPE; | |
1136 | |
263 | 1137 if(get_bit_count(&s->pb)!=0) mpeg4_stuffing(&s->pb); |
1138 put_bits(&s->pb, 16, 0); | |
1139 put_bits(&s->pb, 16, 0x100); /* video obj */ | |
1140 put_bits(&s->pb, 16, 0); | |
1141 put_bits(&s->pb, 16, 0x120); /* video obj layer */ | |
1142 | |
1143 put_bits(&s->pb, 1, 0); /* random access vol */ | |
336 | 1144 put_bits(&s->pb, 8, s->vo_type); /* video obj type indication */ |
263 | 1145 put_bits(&s->pb, 1, 1); /* is obj layer id= yes */ |
1146 put_bits(&s->pb, 4, vo_ver_id); /* is obj layer ver id */ | |
1147 put_bits(&s->pb, 3, 1); /* is obj layer priority */ | |
281
1fc96b02142e
mpeg4 aspect_ratio_info in AVCodecContext (requested by alex)
michaelni
parents:
278
diff
changeset
|
1148 if(s->aspect_ratio_info) |
1fc96b02142e
mpeg4 aspect_ratio_info in AVCodecContext (requested by alex)
michaelni
parents:
278
diff
changeset
|
1149 put_bits(&s->pb, 4, s->aspect_ratio_info);/* aspect ratio info */ |
1fc96b02142e
mpeg4 aspect_ratio_info in AVCodecContext (requested by alex)
michaelni
parents:
278
diff
changeset
|
1150 else |
1fc96b02142e
mpeg4 aspect_ratio_info in AVCodecContext (requested by alex)
michaelni
parents:
278
diff
changeset
|
1151 put_bits(&s->pb, 4, 1); /* aspect ratio info= sqare pixel */ |
336 | 1152 |
1153 if(s->low_delay){ | |
1154 put_bits(&s->pb, 1, 1); /* vol control parameters= yes */ | |
341 | 1155 put_bits(&s->pb, 2, 1); /* chroma format YUV 420/YV12 */ |
336 | 1156 put_bits(&s->pb, 1, s->low_delay); |
1157 put_bits(&s->pb, 1, 0); /* vbv parameters= no */ | |
1158 }else{ | |
1159 put_bits(&s->pb, 1, 0); /* vol control parameters= no */ | |
1160 } | |
1161 | |
263 | 1162 put_bits(&s->pb, 2, RECT_SHAPE); /* vol shape= rectangle */ |
1163 put_bits(&s->pb, 1, 1); /* marker bit */ | |
324 | 1164 |
1165 put_bits(&s->pb, 16, s->time_increment_resolution); | |
263 | 1166 if (s->time_increment_bits < 1) |
1167 s->time_increment_bits = 1; | |
1168 put_bits(&s->pb, 1, 1); /* marker bit */ | |
1169 put_bits(&s->pb, 1, 0); /* fixed vop rate=no */ | |
1170 put_bits(&s->pb, 1, 1); /* marker bit */ | |
1171 put_bits(&s->pb, 13, s->width); /* vol width */ | |
1172 put_bits(&s->pb, 1, 1); /* marker bit */ | |
1173 put_bits(&s->pb, 13, s->height); /* vol height */ | |
1174 put_bits(&s->pb, 1, 1); /* marker bit */ | |
1175 put_bits(&s->pb, 1, 0); /* interlace */ | |
1176 put_bits(&s->pb, 1, 1); /* obmc disable */ | |
1177 if (vo_ver_id == 1) { | |
1178 put_bits(&s->pb, 1, s->vol_sprite_usage=0); /* sprite enable */ | |
1179 }else{ /* vo_ver_id == 2 */ | |
1180 put_bits(&s->pb, 2, s->vol_sprite_usage=0); /* sprite enable */ | |
1181 } | |
1182 put_bits(&s->pb, 1, 0); /* not 8 bit */ | |
1183 put_bits(&s->pb, 1, 0); /* quant type= h263 style*/ | |
1184 if (vo_ver_id != 1) | |
1185 put_bits(&s->pb, 1, s->quarter_sample=0); | |
1186 put_bits(&s->pb, 1, 1); /* complexity estimation disable */ | |
1187 put_bits(&s->pb, 1, 1); /* resync marker disable */ | |
1188 put_bits(&s->pb, 1, 0); /* data partitioned */ | |
1189 if (vo_ver_id != 1){ | |
1190 put_bits(&s->pb, 1, 0); /* newpred */ | |
1191 put_bits(&s->pb, 1, 0); /* reduced res vop */ | |
1192 } | |
1193 put_bits(&s->pb, 1, 0); /* scalability */ | |
1194 | |
262 | 1195 mpeg4_stuffing(&s->pb); |
1196 put_bits(&s->pb, 16, 0); | |
1197 put_bits(&s->pb, 16, 0x1B2); /* user_data */ | |
322 | 1198 sprintf(buf, "FFmpeg v%s / libavcodec build: %s", FFMPEG_VERSION, LIBAVCODEC_BUILD_STR); |
1199 put_string(&s->pb, buf); | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
1200 |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
1201 s->no_rounding = 0; |
263 | 1202 } |
1203 | |
1204 /* write mpeg4 VOP header */ | |
1205 void mpeg4_encode_picture_header(MpegEncContext * s, int picture_number) | |
1206 { | |
324 | 1207 int time_incr; |
1208 int time_div, time_mod; | |
1209 | |
263 | 1210 if(s->pict_type==I_TYPE) mpeg4_encode_vol_header(s); |
324 | 1211 |
1212 //printf("num:%d rate:%d base:%d\n", s->picture_number, s->frame_rate, FRAME_RATE_BASE); | |
1213 | |
263 | 1214 if(get_bit_count(&s->pb)!=0) mpeg4_stuffing(&s->pb); |
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
|
1215 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
|
1216 put_bits(&s->pb, 16, 0x1B6); /* vop header */ |
0 | 1217 put_bits(&s->pb, 2, s->pict_type - 1); /* pict type: I = 0 , P = 1 */ |
324 | 1218 |
327 | 1219 time_div= s->time/s->time_increment_resolution; |
1220 time_mod= s->time%s->time_increment_resolution; | |
324 | 1221 time_incr= time_div - s->last_time_base; |
1222 while(time_incr--) | |
1223 put_bits(&s->pb, 1, 1); | |
1224 | |
0 | 1225 put_bits(&s->pb, 1, 0); |
1226 | |
1227 put_bits(&s->pb, 1, 1); /* marker */ | |
324 | 1228 put_bits(&s->pb, s->time_increment_bits, time_mod); /* time increment */ |
0 | 1229 put_bits(&s->pb, 1, 1); /* marker */ |
1230 put_bits(&s->pb, 1, 1); /* vop coded */ | |
263 | 1231 if ( s->pict_type == P_TYPE |
1232 || (s->pict_type == S_TYPE && s->vol_sprite_usage==GMC_SPRITE)) { | |
1233 s->no_rounding ^= 1; | |
0 | 1234 put_bits(&s->pb, 1, s->no_rounding); /* rounding type */ |
1235 } | |
1236 put_bits(&s->pb, 3, 0); /* intra dc VLC threshold */ | |
263 | 1237 //FIXME sprite stuff |
0 | 1238 |
1239 put_bits(&s->pb, 5, s->qscale); | |
1240 | |
1241 if (s->pict_type != I_TYPE) | |
1242 put_bits(&s->pb, 3, s->f_code); /* fcode_for */ | |
263 | 1243 if (s->pict_type == B_TYPE) |
1244 put_bits(&s->pb, 3, s->b_code); /* fcode_back */ | |
0 | 1245 // printf("****frame %d\n", picture_number); |
1246 } | |
1247 | |
1248 void h263_dc_scale(MpegEncContext * s) | |
1249 { | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1250 #if 1 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1251 const static UINT8 y_tab[32]={ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1252 // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1253 0, 8, 8, 8, 8,10,12,14,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,34,36,38,40,42,44,46 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1254 }; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1255 const static UINT8 c_tab[32]={ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1256 // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1257 0, 8, 8, 8, 8, 9, 9,10,10,11,11,12,12,13,13,14,14,15,15,16,16,17,17,18,18,19,20,21,22,23,24,25 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1258 }; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1259 s->y_dc_scale = y_tab[s->qscale]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1260 s->c_dc_scale = c_tab[s->qscale]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1261 #else |
0 | 1262 int quant; |
1263 quant = s->qscale; | |
1264 /* luminance */ | |
1265 if (quant < 5) | |
1266 s->y_dc_scale = 8; | |
1267 else if (quant > 4 && quant < 9) | |
1268 s->y_dc_scale = (2 * quant); | |
1269 else if (quant > 8 && quant < 25) | |
1270 s->y_dc_scale = (quant + 8); | |
1271 else | |
1272 s->y_dc_scale = (2 * quant - 16); | |
1273 /* chrominance */ | |
1274 if (quant < 5) | |
1275 s->c_dc_scale = 8; | |
1276 else if (quant > 4 && quant < 25) | |
1277 s->c_dc_scale = ((quant + 13) / 2); | |
1278 else | |
1279 s->c_dc_scale = (quant - 6); | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1280 #endif |
0 | 1281 } |
1282 | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1283 static inline int mpeg4_pred_dc(MpegEncContext * s, int n, UINT16 **dc_val_ptr, int *dir_ptr) |
0 | 1284 { |
266 | 1285 int a, b, c, wrap, pred, scale; |
0 | 1286 UINT16 *dc_val; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1287 int dummy; |
0 | 1288 |
1289 /* find prediction */ | |
1290 if (n < 4) { | |
1291 scale = s->y_dc_scale; | |
1292 } else { | |
1293 scale = s->c_dc_scale; | |
1294 } | |
266 | 1295 wrap= s->block_wrap[n]; |
1296 dc_val = s->dc_val[0] + s->block_index[n]; | |
0 | 1297 |
1298 /* B C | |
1299 * A X | |
1300 */ | |
266 | 1301 a = dc_val[ - 1]; |
1302 b = dc_val[ - 1 - wrap]; | |
1303 c = dc_val[ - wrap]; | |
0 | 1304 |
1305 if (abs(a - b) < abs(b - c)) { | |
1306 pred = c; | |
1307 *dir_ptr = 1; /* top */ | |
1308 } else { | |
1309 pred = a; | |
1310 *dir_ptr = 0; /* left */ | |
1311 } | |
1312 /* we assume pred is positive */ | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1313 #ifdef ARCH_X86 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1314 asm volatile ( |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1315 "xorl %%edx, %%edx \n\t" |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1316 "mul %%ecx \n\t" |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1317 : "=d" (pred), "=a"(dummy) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1318 : "a" (pred + (scale >> 1)), "c" (inverse[scale]) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1319 ); |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1320 #else |
0 | 1321 pred = (pred + (scale >> 1)) / scale; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1322 #endif |
0 | 1323 |
1324 /* prepare address for prediction update */ | |
266 | 1325 *dc_val_ptr = &dc_val[0]; |
0 | 1326 |
1327 return pred; | |
1328 } | |
1329 | |
174
ac5075a55488
new IDCT code by Michael Niedermayer (michaelni@gmx.at) - #define SIMPLE_IDCT to enable
arpi_esp
parents:
162
diff
changeset
|
1330 void mpeg4_pred_ac(MpegEncContext * s, INT16 *block, int n, |
0 | 1331 int dir) |
1332 { | |
266 | 1333 int i; |
0 | 1334 INT16 *ac_val, *ac_val1; |
1335 | |
1336 /* find prediction */ | |
266 | 1337 ac_val = s->ac_val[0][0] + s->block_index[n] * 16; |
0 | 1338 ac_val1 = ac_val; |
1339 if (s->ac_pred) { | |
1340 if (dir == 0) { | |
1341 /* left prediction */ | |
1342 ac_val -= 16; | |
1343 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
|
1344 block[block_permute_op(i*8)] += ac_val[i]; |
0 | 1345 } |
1346 } else { | |
1347 /* top prediction */ | |
266 | 1348 ac_val -= 16 * s->block_wrap[n]; |
0 | 1349 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
|
1350 block[block_permute_op(i)] += ac_val[i + 8]; |
0 | 1351 } |
1352 } | |
1353 } | |
1354 /* left copy */ | |
1355 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
|
1356 ac_val1[i] = block[block_permute_op(i * 8)]; |
0 | 1357 /* top copy */ |
1358 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
|
1359 ac_val1[8 + i] = block[block_permute_op(i)]; |
0 | 1360 } |
1361 | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1362 static void mpeg4_inv_pred_ac(MpegEncContext * s, INT16 *block, int n, |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1363 int dir) |
0 | 1364 { |
266 | 1365 int i; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1366 INT16 *ac_val; |
0 | 1367 |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1368 /* find prediction */ |
266 | 1369 ac_val = s->ac_val[0][0] + s->block_index[n] * 16; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1370 |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1371 if (dir == 0) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1372 /* left prediction */ |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1373 ac_val -= 16; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1374 for(i=1;i<8;i++) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1375 block[block_permute_op(i*8)] -= ac_val[i]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1376 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1377 } else { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1378 /* top prediction */ |
266 | 1379 ac_val -= 16 * s->block_wrap[n]; |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1380 for(i=1;i<8;i++) { |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1381 block[block_permute_op(i)] -= ac_val[i + 8]; |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1382 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1383 } |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1384 } |
0 | 1385 |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1386 static inline void mpeg4_encode_dc(MpegEncContext * s, int level, int n) |
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1387 { |
293 | 1388 #if 1 |
1389 level+=256; | |
1390 if (n < 4) { | |
1391 /* luminance */ | |
1392 put_bits(&s->pb, uni_DCtab_lum[level][1], uni_DCtab_lum[level][0]); | |
1393 } else { | |
1394 /* chrominance */ | |
1395 put_bits(&s->pb, uni_DCtab_chrom[level][1], uni_DCtab_chrom[level][0]); | |
1396 } | |
1397 #else | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1398 int size, v; |
0 | 1399 /* find number of bits */ |
1400 size = 0; | |
1401 v = abs(level); | |
1402 while (v) { | |
1403 v >>= 1; | |
1404 size++; | |
1405 } | |
1406 | |
1407 if (n < 4) { | |
1408 /* luminance */ | |
1409 put_bits(&s->pb, DCtab_lum[size][1], DCtab_lum[size][0]); | |
1410 } else { | |
1411 /* chrominance */ | |
1412 put_bits(&s->pb, DCtab_chrom[size][1], DCtab_chrom[size][0]); | |
1413 } | |
1414 | |
1415 /* encode remaining bits */ | |
1416 if (size > 0) { | |
1417 if (level < 0) | |
1418 level = (-level) ^ ((1 << size) - 1); | |
1419 put_bits(&s->pb, size, level); | |
1420 if (size > 8) | |
1421 put_bits(&s->pb, 1, 1); | |
1422 } | |
293 | 1423 #endif |
0 | 1424 } |
1425 | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1426 static void mpeg4_encode_block(MpegEncContext * s, DCTELEM * block, int n, int intra_dc, UINT8 *scan_table) |
0 | 1427 { |
1428 int level, run, last, i, j, last_index, last_non_zero, sign, slevel; | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1429 int code; |
0 | 1430 const RLTable *rl; |
1431 | |
1432 if (s->mb_intra) { | |
1433 /* mpeg4 based DC predictor */ | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1434 mpeg4_encode_dc(s, intra_dc, n); |
0 | 1435 i = 1; |
1436 rl = &rl_intra; | |
1437 } else { | |
1438 i = 0; | |
1439 rl = &rl_inter; | |
1440 } | |
1441 | |
1442 /* AC coefs */ | |
1443 last_index = s->block_last_index[n]; | |
1444 last_non_zero = i - 1; | |
1445 for (; i <= last_index; i++) { | |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
1446 j = scan_table[i]; |
0 | 1447 level = block[j]; |
1448 if (level) { | |
1449 run = i - last_non_zero - 1; | |
1450 last = (i == last_index); | |
1451 sign = 0; | |
1452 slevel = level; | |
1453 if (level < 0) { | |
1454 sign = 1; | |
1455 level = -level; | |
1456 } | |
1457 code = get_rl_index(rl, last, run, level); | |
1458 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
1459 if (code == rl->n) { | |
1460 int level1, run1; | |
1461 level1 = level - rl->max_level[last][run]; | |
1462 if (level1 < 1) | |
1463 goto esc2; | |
1464 code = get_rl_index(rl, last, run, level1); | |
1465 if (code == rl->n) { | |
1466 esc2: | |
1467 put_bits(&s->pb, 1, 1); | |
1468 if (level > MAX_LEVEL) | |
1469 goto esc3; | |
1470 run1 = run - rl->max_run[last][level] - 1; | |
1471 if (run1 < 0) | |
1472 goto esc3; | |
1473 code = get_rl_index(rl, last, run1, level); | |
1474 if (code == rl->n) { | |
1475 esc3: | |
1476 /* third escape */ | |
1477 put_bits(&s->pb, 1, 1); | |
1478 put_bits(&s->pb, 1, last); | |
1479 put_bits(&s->pb, 6, run); | |
1480 put_bits(&s->pb, 1, 1); | |
1481 put_bits(&s->pb, 12, slevel & 0xfff); | |
1482 put_bits(&s->pb, 1, 1); | |
1483 } else { | |
1484 /* second escape */ | |
1485 put_bits(&s->pb, 1, 0); | |
1486 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
1487 put_bits(&s->pb, 1, sign); | |
1488 } | |
1489 } else { | |
1490 /* first escape */ | |
1491 put_bits(&s->pb, 1, 0); | |
1492 put_bits(&s->pb, rl->table_vlc[code][1], rl->table_vlc[code][0]); | |
1493 put_bits(&s->pb, 1, sign); | |
1494 } | |
1495 } else { | |
1496 put_bits(&s->pb, 1, sign); | |
1497 } | |
1498 last_non_zero = i; | |
1499 } | |
1500 } | |
1501 } | |
1502 | |
1503 | |
1504 | |
1505 /***********************************************/ | |
1506 /* decoding */ | |
1507 | |
1508 static VLC intra_MCBPC_vlc; | |
1509 static VLC inter_MCBPC_vlc; | |
1510 static VLC cbpy_vlc; | |
1511 static VLC mv_vlc; | |
1512 static VLC dc_lum, dc_chrom; | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1513 static VLC sprite_trajectory; |
262 | 1514 static VLC mb_type_b_vlc; |
0 | 1515 |
1516 void init_rl(RLTable *rl) | |
1517 { | |
1518 INT8 max_level[MAX_RUN+1], max_run[MAX_LEVEL+1]; | |
1519 UINT8 index_run[MAX_RUN+1]; | |
1520 int last, run, level, start, end, i; | |
1521 | |
1522 /* compute max_level[], max_run[] and index_run[] */ | |
1523 for(last=0;last<2;last++) { | |
1524 if (last == 0) { | |
1525 start = 0; | |
1526 end = rl->last; | |
1527 } else { | |
1528 start = rl->last; | |
1529 end = rl->n; | |
1530 } | |
1531 | |
1532 memset(max_level, 0, MAX_RUN + 1); | |
1533 memset(max_run, 0, MAX_LEVEL + 1); | |
1534 memset(index_run, rl->n, MAX_RUN + 1); | |
1535 for(i=start;i<end;i++) { | |
1536 run = rl->table_run[i]; | |
1537 level = rl->table_level[i]; | |
1538 if (index_run[run] == rl->n) | |
1539 index_run[run] = i; | |
1540 if (level > max_level[run]) | |
1541 max_level[run] = level; | |
1542 if (run > max_run[level]) | |
1543 max_run[level] = run; | |
1544 } | |
1545 rl->max_level[last] = malloc(MAX_RUN + 1); | |
1546 memcpy(rl->max_level[last], max_level, MAX_RUN + 1); | |
1547 rl->max_run[last] = malloc(MAX_LEVEL + 1); | |
1548 memcpy(rl->max_run[last], max_run, MAX_LEVEL + 1); | |
1549 rl->index_run[last] = malloc(MAX_RUN + 1); | |
1550 memcpy(rl->index_run[last], index_run, MAX_RUN + 1); | |
1551 } | |
1552 } | |
1553 | |
1554 void init_vlc_rl(RLTable *rl) | |
1555 { | |
1556 init_vlc(&rl->vlc, 9, rl->n + 1, | |
1557 &rl->table_vlc[0][1], 4, 2, | |
1558 &rl->table_vlc[0][0], 4, 2); | |
1559 } | |
1560 | |
1561 /* init vlcs */ | |
1562 | |
1563 /* XXX: find a better solution to handle static init */ | |
1564 void h263_decode_init_vlc(MpegEncContext *s) | |
1565 { | |
1566 static int done = 0; | |
1567 | |
1568 if (!done) { | |
1569 done = 1; | |
1570 | |
1571 init_vlc(&intra_MCBPC_vlc, 6, 8, | |
1572 intra_MCBPC_bits, 1, 1, | |
1573 intra_MCBPC_code, 1, 1); | |
161 | 1574 init_vlc(&inter_MCBPC_vlc, 9, 25, |
0 | 1575 inter_MCBPC_bits, 1, 1, |
1576 inter_MCBPC_code, 1, 1); | |
1577 init_vlc(&cbpy_vlc, 6, 16, | |
1578 &cbpy_tab[0][1], 2, 1, | |
1579 &cbpy_tab[0][0], 2, 1); | |
1580 init_vlc(&mv_vlc, 9, 33, | |
1581 &mvtab[0][1], 2, 1, | |
1582 &mvtab[0][0], 2, 1); | |
1583 init_rl(&rl_inter); | |
1584 init_rl(&rl_intra); | |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
1585 init_rl(&rl_intra_aic); |
0 | 1586 init_vlc_rl(&rl_inter); |
1587 init_vlc_rl(&rl_intra); | |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
1588 init_vlc_rl(&rl_intra_aic); |
0 | 1589 init_vlc(&dc_lum, 9, 13, |
1590 &DCtab_lum[0][1], 2, 1, | |
1591 &DCtab_lum[0][0], 2, 1); | |
1592 init_vlc(&dc_chrom, 9, 13, | |
1593 &DCtab_chrom[0][1], 2, 1, | |
1594 &DCtab_chrom[0][0], 2, 1); | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1595 init_vlc(&sprite_trajectory, 9, 15, |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1596 &sprite_trajectory_tab[0][1], 4, 2, |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1597 &sprite_trajectory_tab[0][0], 4, 2); |
262 | 1598 init_vlc(&mb_type_b_vlc, 4, 4, |
1599 &mb_type_b_tab[0][1], 2, 1, | |
1600 &mb_type_b_tab[0][0], 2, 1); | |
0 | 1601 } |
1602 } | |
1603 | |
162 | 1604 int h263_decode_gob_header(MpegEncContext *s) |
1605 { | |
1606 unsigned int val, gfid; | |
1607 | |
1608 /* Check for GOB Start Code */ | |
1609 val = show_bits(&s->gb, 16); | |
1610 if (val == 0) { | |
1611 /* We have a GBSC probably with GSTUFF */ | |
1612 skip_bits(&s->gb, 16); /* Drop the zeros */ | |
1613 while (get_bits1(&s->gb) == 0); /* Seek the '1' bit */ | |
1614 #ifdef DEBUG | |
1615 fprintf(stderr,"\nGOB Start Code at MB %d\n", (s->mb_y * s->mb_width) + s->mb_x); | |
1616 #endif | |
1617 s->gob_number = get_bits(&s->gb, 5); /* GN */ | |
1618 gfid = get_bits(&s->gb, 2); /* GFID */ | |
1619 s->qscale = get_bits(&s->gb, 5); /* GQUANT */ | |
1620 #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
|
1621 fprintf(stderr, "\nGN: %u GFID: %u Quant: %u\n", s->gob_number, gfid, s->qscale); |
162 | 1622 #endif |
1623 return 1; | |
1624 } | |
1625 return 0; | |
1626 | |
1627 } | |
1628 | |
290 | 1629 static inline void memsetw(short *tab, int val, int n) |
1630 { | |
1631 int i; | |
1632 for(i=0;i<n;i++) | |
1633 tab[i] = val; | |
1634 } | |
1635 | |
1636 static int mpeg4_resync(MpegEncContext *s) | |
1637 { | |
1638 int state, v, bits; | |
1639 int mb_num_bits= av_log2(s->mb_num - 1) + 1; | |
1640 int header_extension=0, mb_num; | |
1641 int c_wrap, c_xy, l_wrap, l_xy; | |
324 | 1642 int time_increment; |
290 | 1643 //printf("resync at %d %d\n", s->mb_x, s->mb_y); |
1644 //printf("%X\n", show_bits(&s->gb, 24)); | |
1645 | |
1646 if( get_bits_count(&s->gb) > s->gb.size*8-32) | |
1647 return 0; | |
1648 | |
1649 align_get_bits(&s->gb); | |
1650 state = 0xff; | |
1651 for(;;) { | |
1652 v = get_bits(&s->gb, 8); | |
1653 //printf("%X ", v); | |
1654 state = ((state << 8) | v) & 0xffff; | |
1655 if (state == 0) break; | |
1656 if( get_bits_count(&s->gb) > s->gb.size*8-32){ | |
1657 printf("resync failed\n"); | |
1658 return -1; | |
1659 } | |
1660 } | |
1661 //printf("%X\n", show_bits(&s->gb, 24)); | |
1662 bits=0; | |
1663 while(!get_bits1(&s->gb) && bits<30) bits++; | |
1664 if(s->pict_type == P_TYPE && bits != s->f_code-1) | |
1665 printf("marker does not match f_code\n"); | |
1666 //FIXME check bits for B-framess | |
1667 //printf("%X\n", show_bits(&s->gb, 24)); | |
1668 | |
1669 if(s->shape != RECT_SHAPE){ | |
1670 header_extension= get_bits1(&s->gb); | |
1671 //FIXME more stuff here | |
1672 } | |
1673 | |
1674 mb_num= get_bits(&s->gb, mb_num_bits); | |
1675 if(mb_num != s->mb_x + s->mb_y*s->mb_width){ | |
1676 printf("MB-num change not supported %d %d\n", mb_num, s->mb_x + s->mb_y*s->mb_width); | |
1677 // s->mb_x= mb_num % s->mb_width; | |
1678 // s->mb_y= mb_num / s->mb_width; | |
1679 //FIXME many vars are wrong now | |
1680 } | |
1681 | |
1682 if(s->shape != BIN_ONLY_SHAPE){ | |
1683 s->qscale= get_bits(&s->gb, 5); | |
1684 h263_dc_scale(s); | |
1685 } | |
1686 | |
1687 if(s->shape == RECT_SHAPE){ | |
1688 header_extension= get_bits1(&s->gb); | |
1689 } | |
1690 if(header_extension){ | |
1691 int time_incr=0; | |
1692 printf("header extension not really supported\n"); | |
1693 while (get_bits1(&s->gb) != 0) | |
1694 time_incr++; | |
1695 | |
1696 check_marker(&s->gb, "before time_increment in video packed header"); | |
324 | 1697 time_increment= get_bits(&s->gb, s->time_increment_bits); |
290 | 1698 if(s->pict_type!=B_TYPE){ |
324 | 1699 s->last_time_base= s->time_base; |
290 | 1700 s->time_base+= time_incr; |
324 | 1701 s->time= s->time_base*s->time_increment_resolution + time_increment; |
1702 s->pp_time= s->time - s->last_non_b_time; | |
1703 s->last_non_b_time= s->time; | |
290 | 1704 }else{ |
324 | 1705 s->time= (s->last_time_base + time_incr)*s->time_increment_resolution + time_increment; |
1706 s->bp_time= s->last_non_b_time - s->time; | |
290 | 1707 } |
1708 check_marker(&s->gb, "before vop_coding_type in video packed header"); | |
1709 | |
1710 skip_bits(&s->gb, 2); /* vop coding type */ | |
1711 //FIXME not rect stuff here | |
1712 | |
1713 if(s->shape != BIN_ONLY_SHAPE){ | |
1714 skip_bits(&s->gb, 3); /* intra dc vlc threshold */ | |
1715 | |
1716 if(s->pict_type == S_TYPE && s->vol_sprite_usage==GMC_SPRITE && s->num_sprite_warping_points){ | |
1717 mpeg4_decode_sprite_trajectory(s); | |
1718 } | |
1719 | |
1720 //FIXME reduced res stuff here | |
1721 | |
1722 if (s->pict_type != I_TYPE) { | |
1723 s->f_code = get_bits(&s->gb, 3); /* fcode_for */ | |
1724 if(s->f_code==0){ | |
1725 printf("Error, video packet header damaged or not MPEG4 header (f_code=0)\n"); | |
1726 return -1; // makes no sense to continue, as the MV decoding will break very quickly | |
1727 } | |
1728 } | |
1729 if (s->pict_type == B_TYPE) { | |
1730 s->b_code = get_bits(&s->gb, 3); | |
1731 } | |
1732 } | |
1733 | |
1734 } | |
1735 //FIXME new-pred stuff | |
1736 | |
1737 l_wrap= s->block_wrap[0]; | |
1738 l_xy= s->mb_y*l_wrap*2; | |
1739 c_wrap= s->block_wrap[4]; | |
1740 c_xy= s->mb_y*c_wrap; | |
1741 | |
1742 /* clean DC */ | |
1743 memsetw(s->dc_val[0] + l_xy, 1024, l_wrap*3); | |
1744 memsetw(s->dc_val[1] + c_xy, 1024, c_wrap*2); | |
1745 memsetw(s->dc_val[2] + c_xy, 1024, c_wrap*2); | |
1746 | |
1747 /* clean AC */ | |
1748 memset(s->ac_val[0] + l_xy, 0, l_wrap*3*16*sizeof(INT16)); | |
1749 memset(s->ac_val[1] + c_xy, 0, c_wrap*2*16*sizeof(INT16)); | |
1750 memset(s->ac_val[2] + c_xy, 0, c_wrap*2*16*sizeof(INT16)); | |
1751 | |
1752 /* clean MV */ | |
1753 memset(s->motion_val + l_xy, 0, l_wrap*3*2*sizeof(INT16)); | |
1754 // memset(s->motion_val, 0, 2*sizeof(INT16)*(2 + s->mb_width*2)*(2 + s->mb_height*2)); | |
1755 s->resync_x_pos= s->mb_x; | |
1756 s->first_slice_line=1; | |
1757 | |
1758 return 0; | |
1759 } | |
1760 | |
0 | 1761 int h263_decode_mb(MpegEncContext *s, |
1762 DCTELEM block[6][64]) | |
1763 { | |
1764 int cbpc, cbpy, i, cbp, pred_x, pred_y, mx, my, dquant; | |
1765 INT16 *mot_val; | |
110 | 1766 static INT8 quant_tab[4] = { -1, -2, 1, 2 }; |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
1767 |
290 | 1768 if(s->resync_marker){ |
1769 if( s->resync_x_pos == s->mb_x+1 | |
1770 || s->resync_x_pos == s->mb_x){ | |
1771 /* f*ck mpeg4 | |
1772 this is here so we dont need to slowdown h263_pred_motion with it */ | |
1773 if(s->resync_x_pos == s->mb_x+1 && s->mb_x==0){ | |
1774 int xy= s->block_index[0] - s->block_wrap[0]; | |
1775 s->motion_val[xy][0]= s->motion_val[xy+2][0]; | |
1776 s->motion_val[xy][1]= s->motion_val[xy+2][1]; | |
1777 } | |
1778 | |
1779 s->first_slice_line=0; | |
1780 s->resync_x_pos=0; // isnt needed but for cleanness sake ;) | |
1781 } | |
1782 | |
1783 if(show_aligned_bits(&s->gb, 1, 16) == 0){ | |
1784 if( mpeg4_resync(s) < 0 ) return -1; | |
1785 | |
1786 } | |
1787 } | |
1788 | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1789 if (s->pict_type == P_TYPE || s->pict_type==S_TYPE) { |
21 | 1790 if (get_bits1(&s->gb)) { |
0 | 1791 /* skip mb */ |
1792 s->mb_intra = 0; | |
1793 for(i=0;i<6;i++) | |
1794 s->block_last_index[i] = -1; | |
1795 s->mv_dir = MV_DIR_FORWARD; | |
1796 s->mv_type = MV_TYPE_16X16; | |
255 | 1797 if(s->pict_type==S_TYPE && s->vol_sprite_usage==GMC_SPRITE){ |
1798 const int a= s->sprite_warping_accuracy; | |
1799 // int l = (1 << (s->f_code - 1)) * 32; | |
1800 | |
1801 s->mcsel=1; | |
301 | 1802 if(s->divx_version==500 && s->divx_build==413){ |
1803 s->mv[0][0][0] = s->sprite_offset[0][0] / (1<<(a-s->quarter_sample)); | |
1804 s->mv[0][0][1] = s->sprite_offset[0][1] / (1<<(a-s->quarter_sample)); | |
1805 }else{ | |
1806 s->mv[0][0][0] = RSHIFT(s->sprite_offset[0][0], a-s->quarter_sample); | |
1807 s->mv[0][0][1] = RSHIFT(s->sprite_offset[0][1], a-s->quarter_sample); | |
1808 } | |
255 | 1809 /* if (s->mv[0][0][0] < -l) s->mv[0][0][0]= -l; |
1810 else if (s->mv[0][0][0] >= l) s->mv[0][0][0]= l-1; | |
1811 if (s->mv[0][0][1] < -l) s->mv[0][0][1]= -l; | |
1812 else if (s->mv[0][0][1] >= l) s->mv[0][0][1]= l-1;*/ | |
1813 | |
1814 s->mb_skiped = 0; | |
1815 }else{ | |
1816 s->mcsel=0; | |
1817 s->mv[0][0][0] = 0; | |
1818 s->mv[0][0][1] = 0; | |
1819 s->mb_skiped = 1; | |
1820 } | |
0 | 1821 return 0; |
1822 } | |
1823 cbpc = get_vlc(&s->gb, &inter_MCBPC_vlc); | |
144 | 1824 //fprintf(stderr, "\tCBPC: %d", cbpc); |
0 | 1825 if (cbpc < 0) |
1826 return -1; | |
161 | 1827 if (cbpc > 20) |
1828 cbpc+=3; | |
1829 else if (cbpc == 20) | |
1830 fprintf(stderr, "Stuffing !"); | |
144 | 1831 |
0 | 1832 dquant = cbpc & 8; |
1833 s->mb_intra = ((cbpc & 4) != 0); | |
262 | 1834 if (s->mb_intra) goto intra; |
1835 | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1836 if(s->pict_type==S_TYPE && s->vol_sprite_usage==GMC_SPRITE && (cbpc & 16) == 0) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1837 s->mcsel= get_bits1(&s->gb); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1838 else s->mcsel= 0; |
0 | 1839 cbpy = get_vlc(&s->gb, &cbpy_vlc); |
1840 cbp = (cbpc & 3) | ((cbpy ^ 0xf) << 2); | |
1841 if (dquant) { | |
1842 s->qscale += quant_tab[get_bits(&s->gb, 2)]; | |
1843 if (s->qscale < 1) | |
1844 s->qscale = 1; | |
1845 else if (s->qscale > 31) | |
1846 s->qscale = 31; | |
290 | 1847 h263_dc_scale(s); |
0 | 1848 } |
1849 s->mv_dir = MV_DIR_FORWARD; | |
1850 if ((cbpc & 16) == 0) { | |
1851 /* 16x16 motion prediction */ | |
1852 s->mv_type = MV_TYPE_16X16; | |
1853 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
|
1854 if (s->umvplus_dec) |
78 | 1855 mx = h263p_decode_umotion(s, pred_x); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1856 else if(!s->mcsel) |
262 | 1857 mx = h263_decode_motion(s, pred_x, s->f_code); |
255 | 1858 else { |
1859 const int a= s->sprite_warping_accuracy; | |
1860 // int l = (1 << (s->f_code - 1)) * 32; | |
301 | 1861 if(s->divx_version==500 && s->divx_build==413){ |
1862 mx = s->sprite_offset[0][0] / (1<<(a-s->quarter_sample)); | |
1863 }else{ | |
1864 mx = RSHIFT(s->sprite_offset[0][0], a-s->quarter_sample); | |
1865 } | |
1866 // if (mx < -l) mx= -l, printf("C"); | |
1867 // else if (mx >= l) mx= l-1, printf("C"); | |
255 | 1868 } |
0 | 1869 if (mx >= 0xffff) |
1870 return -1; | |
78 | 1871 |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
1872 if (s->umvplus_dec) |
78 | 1873 my = h263p_decode_umotion(s, pred_y); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
1874 else if(!s->mcsel) |
262 | 1875 my = h263_decode_motion(s, pred_y, s->f_code); |
255 | 1876 else{ |
1877 const int a= s->sprite_warping_accuracy; | |
1878 // int l = (1 << (s->f_code - 1)) * 32; | |
301 | 1879 if(s->divx_version==500 && s->divx_build==413){ |
1880 my = s->sprite_offset[0][1] / (1<<(a-s->quarter_sample)); | |
1881 }else{ | |
1882 my = RSHIFT(s->sprite_offset[0][1], a-s->quarter_sample); | |
1883 } | |
1884 // if (my < -l) my= -l, printf("C"); | |
1885 // else if (my >= l) my= l-1, printf("C"); | |
255 | 1886 } |
0 | 1887 if (my >= 0xffff) |
1888 return -1; | |
1889 s->mv[0][0][0] = mx; | |
1890 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
|
1891 /*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
|
1892 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
|
1893 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
|
1894 if (s->umvplus_dec && (mx - pred_x) == 1 && (my - pred_y) == 1) |
78 | 1895 skip_bits1(&s->gb); /* Bit stuffing to prevent PSC */ |
1896 | |
0 | 1897 } else { |
1898 s->mv_type = MV_TYPE_8X8; | |
1899 for(i=0;i<4;i++) { | |
1900 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
|
1901 if (s->umvplus_dec) |
78 | 1902 mx = h263p_decode_umotion(s, pred_x); |
1903 else | |
262 | 1904 mx = h263_decode_motion(s, pred_x, s->f_code); |
0 | 1905 if (mx >= 0xffff) |
1906 return -1; | |
78 | 1907 |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
1908 if (s->umvplus_dec) |
78 | 1909 my = h263p_decode_umotion(s, pred_y); |
1910 else | |
262 | 1911 my = h263_decode_motion(s, pred_y, s->f_code); |
0 | 1912 if (my >= 0xffff) |
1913 return -1; | |
1914 s->mv[0][i][0] = mx; | |
1915 s->mv[0][i][1] = my; | |
79
82e579c37bc3
Moved some H.263+ variables to MpegEncContext to be thread-safe.
pulento
parents:
78
diff
changeset
|
1916 if (s->umvplus_dec && (mx - pred_x) == 1 && (my - pred_y) == 1) |
78 | 1917 skip_bits1(&s->gb); /* Bit stuffing to prevent PSC */ |
0 | 1918 mot_val[0] = mx; |
1919 mot_val[1] = my; | |
1920 } | |
1921 } | |
262 | 1922 } else if(s->pict_type==B_TYPE) { |
1923 int modb1; // first bit of modb | |
1924 int modb2; // second bit of modb | |
1925 int mb_type; | |
324 | 1926 uint16_t time_pp; |
1927 uint16_t time_pb; | |
262 | 1928 int xy; |
1929 | |
1930 s->mb_intra = 0; //B-frames never contain intra blocks | |
1931 s->mcsel=0; // ... true gmc blocks | |
1932 | |
1933 if(s->mb_x==0){ | |
1934 s->last_mv[0][0][0]= | |
1935 s->last_mv[0][0][1]= | |
1936 s->last_mv[1][0][0]= | |
1937 s->last_mv[1][0][1]= 0; | |
327 | 1938 // printf("\n"); |
262 | 1939 } |
1940 | |
1941 /* if we skipped it in the future P Frame than skip it now too */ | |
1942 s->mb_skiped= s->mbskip_table[s->mb_y * s->mb_width + s->mb_x]; // Note, skiptab=0 if last was GMC | |
1943 | |
1944 if(s->mb_skiped){ | |
1945 /* skip mb */ | |
1946 for(i=0;i<6;i++) | |
1947 s->block_last_index[i] = -1; | |
1948 | |
1949 s->mv_dir = MV_DIR_FORWARD; | |
1950 s->mv_type = MV_TYPE_16X16; | |
1951 s->mv[0][0][0] = 0; | |
1952 s->mv[0][0][1] = 0; | |
1953 s->mv[1][0][0] = 0; | |
1954 s->mv[1][0][1] = 0; | |
291 | 1955 //FIXME is this correct? |
1956 /* s->last_mv[0][0][0]= | |
1957 s->last_mv[0][0][1]=0;*/ | |
327 | 1958 // printf("S"); |
262 | 1959 return 0; |
1960 } | |
1961 | |
1962 modb1= get_bits1(&s->gb); | |
1963 if(modb1==0){ | |
1964 modb2= get_bits1(&s->gb); | |
1965 mb_type= get_vlc(&s->gb, &mb_type_b_vlc); | |
1966 if(modb2==0) cbp= get_bits(&s->gb, 6); | |
1967 else cbp=0; | |
1968 if (mb_type && cbp) { | |
1969 if(get_bits1(&s->gb)){ | |
1970 s->qscale +=get_bits1(&s->gb)*4 - 2; | |
1971 if (s->qscale < 1) | |
1972 s->qscale = 1; | |
1973 else if (s->qscale > 31) | |
1974 s->qscale = 31; | |
290 | 1975 h263_dc_scale(s); |
262 | 1976 } |
1977 } | |
1978 }else{ | |
1979 mb_type=4; //like 0 but no vectors coded | |
1980 cbp=0; | |
1981 } | |
1982 s->mv_type = MV_TYPE_16X16; // we'll switch to 8x8 only if the last P frame had 8x8 for this MB and mb_type=0 here | |
1983 mx=my=0; //for case 4, we could put this to the mb_type=4 but than gcc compains about uninitalized mx/my | |
1984 switch(mb_type) | |
1985 { | |
324 | 1986 case 0: /* direct */ |
262 | 1987 mx = h263_decode_motion(s, 0, 1); |
1988 my = h263_decode_motion(s, 0, 1); | |
324 | 1989 case 4: /* direct with mx=my=0 */ |
262 | 1990 s->mv_dir = MV_DIR_FORWARD | MV_DIR_BACKWARD | MV_DIRECT; |
266 | 1991 xy= s->block_index[0]; |
324 | 1992 time_pp= s->pp_time; |
1993 time_pb= time_pp - s->bp_time; | |
262 | 1994 //if(time_pp>3000 )printf("%d %d ", time_pp, time_pb); |
1995 //FIXME 4MV | |
1996 //FIXME avoid divides | |
1997 s->mv[0][0][0] = s->motion_val[xy][0]*time_pb/time_pp + mx; | |
1998 s->mv[0][0][1] = s->motion_val[xy][1]*time_pb/time_pp + my; | |
1999 s->mv[1][0][0] = mx ? s->mv[0][0][0] - s->motion_val[xy][0] | |
2000 : s->motion_val[xy][0]*(time_pb - time_pp)/time_pp + mx; | |
2001 s->mv[1][0][1] = my ? s->mv[0][0][1] - s->motion_val[xy][1] | |
2002 : s->motion_val[xy][1]*(time_pb - time_pp)/time_pp + my; | |
2003 /* s->mv[0][0][0] = | |
2004 s->mv[0][0][1] = | |
2005 s->mv[1][0][0] = | |
2006 s->mv[1][0][1] = 1000;*/ | |
327 | 2007 // printf("D"); |
262 | 2008 break; |
2009 case 1: | |
2010 s->mv_dir = MV_DIR_FORWARD | MV_DIR_BACKWARD; | |
2011 mx = h263_decode_motion(s, s->last_mv[0][0][0], s->f_code); | |
2012 my = h263_decode_motion(s, s->last_mv[0][0][1], s->f_code); | |
2013 s->last_mv[0][0][0]= s->mv[0][0][0] = mx; | |
2014 s->last_mv[0][0][1]= s->mv[0][0][1] = my; | |
2015 | |
2016 mx = h263_decode_motion(s, s->last_mv[1][0][0], s->b_code); | |
2017 my = h263_decode_motion(s, s->last_mv[1][0][1], s->b_code); | |
2018 s->last_mv[1][0][0]= s->mv[1][0][0] = mx; | |
2019 s->last_mv[1][0][1]= s->mv[1][0][1] = my; | |
327 | 2020 // printf("I"); |
262 | 2021 break; |
2022 case 2: | |
2023 s->mv_dir = MV_DIR_BACKWARD; | |
2024 mx = h263_decode_motion(s, s->last_mv[1][0][0], s->b_code); | |
2025 my = h263_decode_motion(s, s->last_mv[1][0][1], s->b_code); | |
2026 s->last_mv[1][0][0]= s->mv[1][0][0] = mx; | |
2027 s->last_mv[1][0][1]= s->mv[1][0][1] = my; | |
327 | 2028 // printf("B"); |
262 | 2029 break; |
2030 case 3: | |
2031 s->mv_dir = MV_DIR_FORWARD; | |
2032 mx = h263_decode_motion(s, s->last_mv[0][0][0], s->f_code); | |
2033 my = h263_decode_motion(s, s->last_mv[0][0][1], s->f_code); | |
2034 s->last_mv[0][0][0]= s->mv[0][0][0] = mx; | |
2035 s->last_mv[0][0][1]= s->mv[0][0][1] = my; | |
327 | 2036 // printf("F"); |
262 | 2037 break; |
2038 default: return -1; | |
2039 } | |
2040 } else { /* I-Frame */ | |
2041 cbpc = get_vlc(&s->gb, &intra_MCBPC_vlc); | |
2042 if (cbpc < 0) | |
2043 return -1; | |
2044 dquant = cbpc & 4; | |
2045 s->mb_intra = 1; | |
2046 intra: | |
0 | 2047 s->ac_pred = 0; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2048 if (s->h263_pred || s->h263_aic) { |
21 | 2049 s->ac_pred = get_bits1(&s->gb); |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2050 if (s->ac_pred && s->h263_aic) |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2051 s->h263_aic_dir = get_bits1(&s->gb); |
0 | 2052 } |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2053 if (s->h263_aic) { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2054 s->y_dc_scale = 2 * s->qscale; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2055 s->c_dc_scale = 2 * s->qscale; |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2056 } |
0 | 2057 cbpy = get_vlc(&s->gb, &cbpy_vlc); |
2058 cbp = (cbpc & 3) | (cbpy << 2); | |
2059 if (dquant) { | |
2060 s->qscale += quant_tab[get_bits(&s->gb, 2)]; | |
2061 if (s->qscale < 1) | |
2062 s->qscale = 1; | |
2063 else if (s->qscale > 31) | |
2064 s->qscale = 31; | |
290 | 2065 h263_dc_scale(s); |
0 | 2066 } |
2067 } | |
2068 | |
2069 /* decode each block */ | |
2070 if (s->h263_pred) { | |
2071 for (i = 0; i < 6; i++) { | |
2072 if (mpeg4_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1) < 0) | |
2073 return -1; | |
2074 } | |
2075 } else { | |
2076 for (i = 0; i < 6; i++) { | |
2077 if (h263_decode_block(s, block[i], i, (cbp >> (5 - i)) & 1) < 0) | |
2078 return -1; | |
2079 } | |
2080 } | |
2081 return 0; | |
2082 } | |
2083 | |
262 | 2084 static int h263_decode_motion(MpegEncContext * s, int pred, int f_code) |
0 | 2085 { |
2086 int code, val, sign, shift, l, m; | |
2087 | |
2088 code = get_vlc(&s->gb, &mv_vlc); | |
2089 if (code < 0) | |
2090 return 0xffff; | |
2091 | |
2092 if (code == 0) | |
2093 return pred; | |
21 | 2094 sign = get_bits1(&s->gb); |
262 | 2095 shift = f_code - 1; |
0 | 2096 val = (code - 1) << shift; |
2097 if (shift > 0) | |
2098 val |= get_bits(&s->gb, shift); | |
2099 val++; | |
2100 if (sign) | |
2101 val = -val; | |
2102 val += pred; | |
2103 | |
2104 /* modulo decoding */ | |
2105 if (!s->h263_long_vectors) { | |
262 | 2106 l = (1 << (f_code - 1)) * 32; |
0 | 2107 m = 2 * l; |
2108 if (val < -l) { | |
2109 val += m; | |
2110 } else if (val >= l) { | |
2111 val -= m; | |
2112 } | |
2113 } else { | |
2114 /* horrible h263 long vector mode */ | |
2115 if (pred < -31 && val < -63) | |
2116 val += 64; | |
2117 if (pred > 32 && val > 63) | |
2118 val -= 64; | |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
2119 |
0 | 2120 } |
2121 return val; | |
2122 } | |
2123 | |
78 | 2124 /* Decodes RVLC of H.263+ UMV */ |
2125 static int h263p_decode_umotion(MpegEncContext * s, int pred) | |
2126 { | |
2127 int code = 0, sign; | |
2128 | |
2129 if (get_bits1(&s->gb)) /* Motion difference = 0 */ | |
2130 return pred; | |
2131 | |
2132 code = 2 + get_bits1(&s->gb); | |
2133 | |
2134 while (get_bits1(&s->gb)) | |
2135 { | |
2136 code <<= 1; | |
2137 code += get_bits1(&s->gb); | |
2138 } | |
2139 sign = code & 1; | |
2140 code >>= 1; | |
2141 | |
2142 code = (sign) ? (pred - code) : (pred + code); | |
2143 #ifdef DEBUG | |
2144 fprintf(stderr,"H.263+ UMV Motion = %d\n", code); | |
2145 #endif | |
2146 return code; | |
2147 | |
2148 } | |
2149 | |
0 | 2150 static int h263_decode_block(MpegEncContext * s, DCTELEM * block, |
2151 int n, int coded) | |
2152 { | |
2153 int code, level, i, j, last, run; | |
2154 RLTable *rl = &rl_inter; | |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2155 const UINT8 *scan_table; |
0 | 2156 |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2157 scan_table = zigzag_direct; |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2158 if (s->h263_aic && s->mb_intra) { |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2159 rl = &rl_intra_aic; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2160 i = 0; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2161 if (s->ac_pred) { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2162 if (s->h263_aic_dir) |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2163 scan_table = ff_alternate_vertical_scan; /* left */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2164 else |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2165 scan_table = ff_alternate_horizontal_scan; /* top */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2166 } |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2167 } else if (s->mb_intra) { |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2168 /* DC coef */ |
0 | 2169 if (s->h263_rv10 && s->rv10_version == 3 && s->pict_type == I_TYPE) { |
2170 int component, diff; | |
2171 component = (n <= 3 ? 0 : n - 4 + 1); | |
2172 level = s->last_dc[component]; | |
2173 if (s->rv10_first_dc_coded[component]) { | |
2174 diff = rv_decode_dc(s, n); | |
2175 if (diff == 0xffff) | |
2176 return -1; | |
2177 level += diff; | |
2178 level = level & 0xff; /* handle wrap round */ | |
2179 s->last_dc[component] = level; | |
2180 } else { | |
2181 s->rv10_first_dc_coded[component] = 1; | |
2182 } | |
2183 } else { | |
2184 level = get_bits(&s->gb, 8); | |
2185 if (level == 255) | |
2186 level = 128; | |
2187 } | |
2188 block[0] = level; | |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2189 i = 1; |
0 | 2190 } else { |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2191 i = 0; |
0 | 2192 } |
2193 if (!coded) { | |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2194 if (s->mb_intra && s->h263_aic) |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2195 goto not_coded; |
0 | 2196 s->block_last_index[n] = i - 1; |
2197 return 0; | |
2198 } | |
2199 | |
2200 for(;;) { | |
2201 code = get_vlc(&s->gb, &rl->vlc); | |
2202 if (code < 0) | |
2203 return -1; | |
2204 if (code == rl->n) { | |
2205 /* escape */ | |
21 | 2206 last = get_bits1(&s->gb); |
0 | 2207 run = get_bits(&s->gb, 6); |
2208 level = (INT8)get_bits(&s->gb, 8); | |
2209 if (s->h263_rv10 && level == -128) { | |
2210 /* XXX: should patch encoder too */ | |
2211 level = get_bits(&s->gb, 12); | |
2212 level = (level << 20) >> 20; | |
2213 } | |
2214 } else { | |
2215 run = rl->table_run[code]; | |
2216 level = rl->table_level[code]; | |
2217 last = code >= rl->last; | |
21 | 2218 if (get_bits1(&s->gb)) |
0 | 2219 level = -level; |
2220 } | |
2221 i += run; | |
2222 if (i >= 64) | |
2223 return -1; | |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2224 j = scan_table[i]; |
0 | 2225 block[j] = level; |
2226 if (last) | |
2227 break; | |
2228 i++; | |
2229 } | |
249
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2230 not_coded: |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2231 if (s->mb_intra && s->h263_aic) { |
42a0b7b16738
- Bug fixes in H.263+ Advanced INTRA Coding decoder.
pulento
parents:
248
diff
changeset
|
2232 h263_pred_acdc(s, block, n); |
265
4e9e728021d8
use ac prediction in mpeg4 encoding (5% smaller intra-blocks/keyframes)
michaelni
parents:
264
diff
changeset
|
2233 i = 63; |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2234 } |
0 | 2235 s->block_last_index[n] = i; |
2236 return 0; | |
2237 } | |
2238 | |
2239 static int mpeg4_decode_dc(MpegEncContext * s, int n, int *dir_ptr) | |
2240 { | |
2241 int level, pred, code; | |
2242 UINT16 *dc_val; | |
2243 | |
2244 if (n < 4) | |
2245 code = get_vlc(&s->gb, &dc_lum); | |
2246 else | |
2247 code = get_vlc(&s->gb, &dc_chrom); | |
2248 if (code < 0) | |
2249 return -1; | |
2250 if (code == 0) { | |
2251 level = 0; | |
2252 } else { | |
2253 level = get_bits(&s->gb, code); | |
2254 if ((level >> (code - 1)) == 0) /* if MSB not set it is negative*/ | |
2255 level = - (level ^ ((1 << code) - 1)); | |
2256 if (code > 8) | |
21 | 2257 skip_bits1(&s->gb); /* marker */ |
0 | 2258 } |
2259 | |
2260 pred = mpeg4_pred_dc(s, n, &dc_val, dir_ptr); | |
2261 level += pred; | |
2262 if (level < 0) | |
2263 level = 0; | |
2264 if (n < 4) { | |
2265 *dc_val = level * s->y_dc_scale; | |
2266 } else { | |
2267 *dc_val = level * s->c_dc_scale; | |
2268 } | |
2269 return level; | |
2270 } | |
2271 | |
2272 static int mpeg4_decode_block(MpegEncContext * s, DCTELEM * block, | |
2273 int n, int coded) | |
2274 { | |
2275 int code, level, i, j, last, run; | |
2276 int dc_pred_dir; | |
2277 RLTable *rl; | |
2278 const UINT8 *scan_table; | |
2279 | |
2280 if (s->mb_intra) { | |
2281 /* DC coef */ | |
2282 level = mpeg4_decode_dc(s, n, &dc_pred_dir); | |
2283 if (level < 0) | |
2284 return -1; | |
2285 block[0] = level; | |
2286 i = 1; | |
2287 if (!coded) | |
2288 goto not_coded; | |
2289 rl = &rl_intra; | |
2290 if (s->ac_pred) { | |
2291 if (dc_pred_dir == 0) | |
2292 scan_table = ff_alternate_vertical_scan; /* left */ | |
2293 else | |
2294 scan_table = ff_alternate_horizontal_scan; /* top */ | |
2295 } else { | |
2296 scan_table = zigzag_direct; | |
2297 } | |
2298 } else { | |
2299 i = 0; | |
2300 if (!coded) { | |
2301 s->block_last_index[n] = i - 1; | |
2302 return 0; | |
2303 } | |
2304 rl = &rl_inter; | |
2305 scan_table = zigzag_direct; | |
2306 } | |
2307 | |
2308 for(;;) { | |
2309 code = get_vlc(&s->gb, &rl->vlc); | |
2310 if (code < 0) | |
2311 return -1; | |
2312 if (code == rl->n) { | |
2313 /* escape */ | |
21 | 2314 if (get_bits1(&s->gb) != 0) { |
2315 if (get_bits1(&s->gb) != 0) { | |
0 | 2316 /* third escape */ |
21 | 2317 last = get_bits1(&s->gb); |
0 | 2318 run = get_bits(&s->gb, 6); |
21 | 2319 get_bits1(&s->gb); /* marker */ |
0 | 2320 level = get_bits(&s->gb, 12); |
2321 level = (level << 20) >> 20; /* sign extend */ | |
21 | 2322 skip_bits1(&s->gb); /* marker */ |
0 | 2323 } else { |
2324 /* second escape */ | |
2325 code = get_vlc(&s->gb, &rl->vlc); | |
2326 if (code < 0 || code >= rl->n) | |
2327 return -1; | |
2328 run = rl->table_run[code]; | |
2329 level = rl->table_level[code]; | |
2330 last = code >= rl->last; | |
2331 run += rl->max_run[last][level] + 1; | |
21 | 2332 if (get_bits1(&s->gb)) |
0 | 2333 level = -level; |
2334 } | |
2335 } else { | |
2336 /* first escape */ | |
2337 code = get_vlc(&s->gb, &rl->vlc); | |
2338 if (code < 0 || code >= rl->n) | |
2339 return -1; | |
2340 run = rl->table_run[code]; | |
2341 level = rl->table_level[code]; | |
2342 last = code >= rl->last; | |
2343 level += rl->max_level[last][run]; | |
21 | 2344 if (get_bits1(&s->gb)) |
0 | 2345 level = -level; |
2346 } | |
2347 } else { | |
2348 run = rl->table_run[code]; | |
2349 level = rl->table_level[code]; | |
2350 last = code >= rl->last; | |
21 | 2351 if (get_bits1(&s->gb)) |
0 | 2352 level = -level; |
2353 } | |
2354 i += run; | |
2355 if (i >= 64) | |
2356 return -1; | |
2357 j = scan_table[i]; | |
2358 block[j] = level; | |
2359 i++; | |
2360 if (last) | |
2361 break; | |
2362 } | |
2363 not_coded: | |
2364 if (s->mb_intra) { | |
2365 mpeg4_pred_ac(s, block, n, dc_pred_dir); | |
2366 if (s->ac_pred) { | |
2367 i = 64; /* XXX: not optimal */ | |
2368 } | |
2369 } | |
2370 s->block_last_index[n] = i - 1; | |
2371 return 0; | |
2372 } | |
2373 | |
2374 /* most is hardcoded. should extend to handle all h263 streams */ | |
2375 int h263_decode_picture_header(MpegEncContext *s) | |
2376 { | |
2377 int format, width, height; | |
2378 | |
355 | 2379 /* picture start code */ |
2380 if (get_bits(&s->gb, 22) != 0x20) { | |
2381 fprintf(stderr, "Bad picture start code\n"); | |
0 | 2382 return -1; |
355 | 2383 } |
2384 /* temporal reference */ | |
231 | 2385 s->picture_number = get_bits(&s->gb, 8); /* picture timestamp */ |
355 | 2386 |
2387 /* PTYPE starts here */ | |
2388 if (get_bits1(&s->gb) != 1) { | |
2389 /* marker */ | |
2390 fprintf(stderr, "Bad marker\n"); | |
2391 return -1; | |
2392 } | |
2393 if (get_bits1(&s->gb) != 0) { | |
2394 fprintf(stderr, "Bad H263 id\n"); | |
0 | 2395 return -1; /* h263 id */ |
355 | 2396 } |
21 | 2397 skip_bits1(&s->gb); /* split screen off */ |
2398 skip_bits1(&s->gb); /* camera off */ | |
2399 skip_bits1(&s->gb); /* freeze picture release off */ | |
0 | 2400 |
155
3c3449bce692
- Bug fix on MV prediction for MPEG4 caused by new H.263 GOB code.
pulento
parents:
154
diff
changeset
|
2401 /* Reset GOB number */ |
154
f914f710b8d0
- Fixed a bug on H.263 MV prediction for MB on GOBs limits.
pulento
parents:
144
diff
changeset
|
2402 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
|
2403 |
0 | 2404 format = get_bits(&s->gb, 3); |
355 | 2405 /* |
2406 0 forbidden | |
2407 1 sub-QCIF | |
2408 10 QCIF | |
2409 7 extended PTYPE (PLUSPTYPE) | |
2410 */ | |
0 | 2411 |
161 | 2412 if (format != 7 && format != 6) { |
0 | 2413 s->h263_plus = 0; |
2414 /* H.263v1 */ | |
2415 width = h263_format[format][0]; | |
2416 height = h263_format[format][1]; | |
2417 if (!width) | |
2418 return -1; | |
161 | 2419 |
2420 s->width = width; | |
2421 s->height = height; | |
21 | 2422 s->pict_type = I_TYPE + get_bits1(&s->gb); |
0 | 2423 |
21 | 2424 s->unrestricted_mv = get_bits1(&s->gb); |
0 | 2425 s->h263_long_vectors = s->unrestricted_mv; |
2426 | |
355 | 2427 if (get_bits1(&s->gb) != 0) { |
2428 fprintf(stderr, "H263 SAC not supported\n"); | |
0 | 2429 return -1; /* SAC: off */ |
355 | 2430 } |
161 | 2431 if (get_bits1(&s->gb) != 0) { |
2432 s->mv_type = MV_TYPE_8X8; /* Advanced prediction mode */ | |
2433 } | |
2434 | |
355 | 2435 if (get_bits1(&s->gb) != 0) { |
2436 fprintf(stderr, "H263 PB frame not supported\n"); | |
0 | 2437 return -1; /* not PB frame */ |
355 | 2438 } |
0 | 2439 s->qscale = get_bits(&s->gb, 5); |
21 | 2440 skip_bits1(&s->gb); /* Continuous Presence Multipoint mode: off */ |
0 | 2441 } else { |
161 | 2442 int ufep; |
2443 | |
0 | 2444 /* H.263v2 */ |
161 | 2445 s->h263_plus = 1; |
2446 ufep = get_bits(&s->gb, 3); /* Update Full Extended PTYPE */ | |
355 | 2447 |
2448 /* ufep other than 0 and 1 are reserved */ | |
161 | 2449 if (ufep == 1) { |
2450 /* OPPTYPE */ | |
2451 format = get_bits(&s->gb, 3); | |
355 | 2452 dprintf("ufep=1, format: %d\n", format); |
161 | 2453 skip_bits(&s->gb,1); /* Custom PCF */ |
2454 s->umvplus_dec = get_bits(&s->gb, 1); /* Unrestricted Motion Vector */ | |
2455 skip_bits1(&s->gb); /* Syntax-based Arithmetic Coding (SAC) */ | |
2456 if (get_bits1(&s->gb) != 0) { | |
2457 s->mv_type = MV_TYPE_8X8; /* Advanced prediction mode */ | |
2458 } | |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2459 if (get_bits1(&s->gb) != 0) { /* Advanced Intra Coding (AIC) */ |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2460 s->h263_aic = 1; |
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2461 } |
355 | 2462 |
248
56ee684c48bb
- H.263+ decoder support for Advanded INTRA Coding (buggy)
pulento
parents:
245
diff
changeset
|
2463 skip_bits(&s->gb, 7); |
355 | 2464 /* these are the 7 bits: (in order of appearence */ |
2465 /* Deblocking Filter */ | |
2466 /* Slice Structured */ | |
2467 /* Reference Picture Selection */ | |
2468 /* Independent Segment Decoding */ | |
2469 /* Alternative Inter VLC */ | |
2470 /* Modified Quantization */ | |
2471 /* Prevent start code emulation */ | |
2472 | |
161 | 2473 skip_bits(&s->gb, 3); /* Reserved */ |
355 | 2474 } else if (ufep != 0) { |
2475 fprintf(stderr, "Bad UFEP type (%d)\n", ufep); | |
0 | 2476 return -1; |
355 | 2477 } |
161 | 2478 |
78 | 2479 /* MPPTYPE */ |
355 | 2480 s->pict_type = get_bits(&s->gb, 3) + I_TYPE; |
2481 dprintf("pict_type: %d\n", s->pict_type); | |
0 | 2482 if (s->pict_type != I_TYPE && |
2483 s->pict_type != P_TYPE) | |
2484 return -1; | |
250
3449316664b5
- Bug fix on RTYPE (rounding type) not being honoured by H.263+ decoder.
pulento
parents:
249
diff
changeset
|
2485 skip_bits(&s->gb, 2); |
3449316664b5
- Bug fix on RTYPE (rounding type) not being honoured by H.263+ decoder.
pulento
parents:
249
diff
changeset
|
2486 s->no_rounding = get_bits1(&s->gb); |
355 | 2487 dprintf("RTYPE: %d\n", s->no_rounding); |
250
3449316664b5
- Bug fix on RTYPE (rounding type) not being honoured by H.263+ decoder.
pulento
parents:
249
diff
changeset
|
2488 skip_bits(&s->gb, 4); |
78 | 2489 |
2490 /* Get the picture dimensions */ | |
161 | 2491 if (ufep) { |
2492 if (format == 6) { | |
2493 /* Custom Picture Format (CPFMT) */ | |
355 | 2494 s->aspect_ratio_info = get_bits(&s->gb, 4); |
2495 dprintf("aspect: %d\n", s->aspect_ratio_info); | |
2496 /* aspect ratios: | |
2497 0 - forbidden | |
2498 1 - 1:1 | |
2499 2 - 12:11 (CIF 4:3) | |
2500 3 - 10:11 (525-type 4:3) | |
2501 4 - 16:11 (CIF 16:9) | |
2502 5 - 40:33 (525-type 16:9) | |
2503 6-14 - reserved | |
2504 */ | |
161 | 2505 width = (get_bits(&s->gb, 9) + 1) * 4; |
2506 skip_bits1(&s->gb); | |
2507 height = get_bits(&s->gb, 9) * 4; | |
355 | 2508 dprintf("\nH.263+ Custom picture: %dx%d\n",width,height); |
2509 if (s->aspect_ratio_info == EXTENDED_PAR) { | |
2510 /* aspected dimensions */ | |
2511 skip_bits(&s->gb, 8); /* width */ | |
2512 skip_bits(&s->gb, 8); /* height */ | |
2513 } | |
2514 } else { | |
161 | 2515 width = h263_format[format][0]; |
2516 height = h263_format[format][1]; | |
2517 } | |
2518 if ((width == 0) || (height == 0)) | |
2519 return -1; | |
2520 s->width = width; | |
2521 s->height = height; | |
2522 if (s->umvplus_dec) { | |
2523 skip_bits1(&s->gb); /* Unlimited Unrestricted Motion Vectors Indicator (UUI) */ | |
2524 } | |
78 | 2525 } |
2526 | |
0 | 2527 s->qscale = get_bits(&s->gb, 5); |
2528 } | |
2529 /* PEI */ | |
21 | 2530 while (get_bits1(&s->gb) != 0) { |
2531 skip_bits(&s->gb, 8); | |
0 | 2532 } |
2533 s->f_code = 1; | |
2534 return 0; | |
2535 } | |
2536 | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2537 static void mpeg4_decode_sprite_trajectory(MpegEncContext * s) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2538 { |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2539 int i; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2540 int a= 2<<s->sprite_warping_accuracy; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2541 int rho= 3-s->sprite_warping_accuracy; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2542 int r=16/a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2543 const int vop_ref[4][2]= {{0,0}, {s->width,0}, {0, s->height}, {s->width, s->height}}; // only true for rectangle shapes |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2544 int d[4][2]={{0,0}, {0,0}, {0,0}, {0,0}}; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2545 int sprite_ref[4][2]; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2546 int virtual_ref[2][2]; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2547 int w2, h2; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2548 int alpha=0, beta=0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2549 int w= s->width; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2550 int h= s->height; |
255 | 2551 //printf("SP %d\n", s->sprite_warping_accuracy); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2552 for(i=0; i<s->num_sprite_warping_points; i++){ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2553 int length; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2554 int x=0, y=0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2555 |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2556 length= get_vlc(&s->gb, &sprite_trajectory); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2557 if(length){ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2558 x= get_bits(&s->gb, length); |
255 | 2559 //printf("lx %d %d\n", length, x); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2560 if ((x >> (length - 1)) == 0) /* if MSB not set it is negative*/ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2561 x = - (x ^ ((1 << length) - 1)); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2562 } |
255 | 2563 if(!(s->divx_version==500 && s->divx_build==413)) skip_bits1(&s->gb); /* marker bit */ |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2564 |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2565 length= get_vlc(&s->gb, &sprite_trajectory); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2566 if(length){ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2567 y=get_bits(&s->gb, length); |
255 | 2568 //printf("ly %d %d\n", length, y); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2569 if ((y >> (length - 1)) == 0) /* if MSB not set it is negative*/ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2570 y = - (y ^ ((1 << length) - 1)); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2571 } |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2572 skip_bits1(&s->gb); /* marker bit */ |
255 | 2573 //printf("%d %d %d %d\n", x, y, i, s->sprite_warping_accuracy); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2574 //if(i>0 && (x!=0 || y!=0)) printf("AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\n"); |
255 | 2575 //x=y=0; |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2576 d[i][0]= x; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2577 d[i][1]= y; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2578 } |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2579 |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2580 while((1<<alpha)<w) alpha++; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2581 while((1<<beta )<h) beta++; // there seems to be a typo in the mpeg4 std for the definition of w' and h' |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2582 w2= 1<<alpha; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2583 h2= 1<<beta; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2584 |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2585 // Note, the 4th point isnt used for GMC |
262 | 2586 if(s->divx_version==500 && s->divx_build==413){ |
2587 sprite_ref[0][0]= a*vop_ref[0][0] + d[0][0]; | |
2588 sprite_ref[0][1]= a*vop_ref[0][1] + d[0][1]; | |
2589 sprite_ref[1][0]= a*vop_ref[1][0] + d[0][0] + d[1][0]; | |
2590 sprite_ref[1][1]= a*vop_ref[1][1] + d[0][1] + d[1][1]; | |
2591 sprite_ref[2][0]= a*vop_ref[2][0] + d[0][0] + d[2][0]; | |
2592 sprite_ref[2][1]= a*vop_ref[2][1] + d[0][1] + d[2][1]; | |
2593 } else { | |
2594 sprite_ref[0][0]= (a>>1)*(2*vop_ref[0][0] + d[0][0]); | |
2595 sprite_ref[0][1]= (a>>1)*(2*vop_ref[0][1] + d[0][1]); | |
2596 sprite_ref[1][0]= (a>>1)*(2*vop_ref[1][0] + d[0][0] + d[1][0]); | |
2597 sprite_ref[1][1]= (a>>1)*(2*vop_ref[1][1] + d[0][1] + d[1][1]); | |
2598 sprite_ref[2][0]= (a>>1)*(2*vop_ref[2][0] + d[0][0] + d[2][0]); | |
2599 sprite_ref[2][1]= (a>>1)*(2*vop_ref[2][1] + d[0][1] + d[2][1]); | |
2600 } | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2601 /* sprite_ref[3][0]= (a>>1)*(2*vop_ref[3][0] + d[0][0] + d[1][0] + d[2][0] + d[3][0]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2602 sprite_ref[3][1]= (a>>1)*(2*vop_ref[3][1] + d[0][1] + d[1][1] + d[2][1] + d[3][1]); */ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2603 |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2604 // this is mostly identical to the mpeg4 std (and is totally unreadable because of that ...) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2605 // perhaps it should be reordered to be more readable ... |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2606 // the idea behind this virtual_ref mess is to be able to use shifts later per pixel instead of divides |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2607 // so the distance between points is converted from w&h based to w2&h2 based which are of the 2^x form |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2608 virtual_ref[0][0]= 16*(vop_ref[0][0] + w2) |
255 | 2609 + RDIV(((w - w2)*(r*sprite_ref[0][0] - 16*vop_ref[0][0]) + w2*(r*sprite_ref[1][0] - 16*vop_ref[1][0])),w); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2610 virtual_ref[0][1]= 16*vop_ref[0][1] |
255 | 2611 + RDIV(((w - w2)*(r*sprite_ref[0][1] - 16*vop_ref[0][1]) + w2*(r*sprite_ref[1][1] - 16*vop_ref[1][1])),w); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2612 virtual_ref[1][0]= 16*vop_ref[0][0] |
255 | 2613 + RDIV(((h - h2)*(r*sprite_ref[0][0] - 16*vop_ref[0][0]) + h2*(r*sprite_ref[2][0] - 16*vop_ref[2][0])),h); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2614 virtual_ref[1][1]= 16*(vop_ref[0][1] + h2) |
255 | 2615 + RDIV(((h - h2)*(r*sprite_ref[0][1] - 16*vop_ref[0][1]) + h2*(r*sprite_ref[2][1] - 16*vop_ref[2][1])),h); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2616 |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2617 switch(s->num_sprite_warping_points) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2618 { |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2619 case 0: |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2620 s->sprite_offset[0][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2621 s->sprite_offset[0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2622 s->sprite_offset[1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2623 s->sprite_offset[1][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2624 s->sprite_delta[0][0][0]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2625 s->sprite_delta[0][0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2626 s->sprite_delta[0][1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2627 s->sprite_delta[0][1][1]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2628 s->sprite_delta[1][0][0]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2629 s->sprite_delta[1][0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2630 s->sprite_delta[1][1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2631 s->sprite_delta[1][1][1]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2632 s->sprite_shift[0][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2633 s->sprite_shift[0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2634 s->sprite_shift[1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2635 s->sprite_shift[1][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2636 break; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2637 case 1: //GMC only |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2638 s->sprite_offset[0][0]= sprite_ref[0][0] - a*vop_ref[0][0]; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2639 s->sprite_offset[0][1]= sprite_ref[0][1] - a*vop_ref[0][1]; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2640 s->sprite_offset[1][0]= ((sprite_ref[0][0]>>1)|(sprite_ref[0][0]&1)) - a*(vop_ref[0][0]/2); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2641 s->sprite_offset[1][1]= ((sprite_ref[0][1]>>1)|(sprite_ref[0][1]&1)) - a*(vop_ref[0][1]/2); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2642 s->sprite_delta[0][0][0]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2643 s->sprite_delta[0][0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2644 s->sprite_delta[0][1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2645 s->sprite_delta[0][1][1]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2646 s->sprite_delta[1][0][0]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2647 s->sprite_delta[1][0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2648 s->sprite_delta[1][1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2649 s->sprite_delta[1][1][1]= a; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2650 s->sprite_shift[0][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2651 s->sprite_shift[0][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2652 s->sprite_shift[1][0]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2653 s->sprite_shift[1][1]= 0; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2654 break; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2655 case 2: |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2656 case 3: //FIXME |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2657 s->sprite_offset[0][0]= (sprite_ref[0][0]<<(alpha+rho)) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2658 + ((-r*sprite_ref[0][0] + virtual_ref[0][0])*(-vop_ref[0][0]) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2659 +( r*sprite_ref[0][1] - virtual_ref[0][1])*(-vop_ref[0][1])); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2660 s->sprite_offset[0][1]= (sprite_ref[0][1]<<(alpha+rho)) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2661 + ((-r*sprite_ref[0][1] + virtual_ref[0][1])*(-vop_ref[0][0]) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2662 +(-r*sprite_ref[0][0] + virtual_ref[0][0])*(-vop_ref[0][1])); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2663 s->sprite_offset[1][0]= ((-r*sprite_ref[0][0] + virtual_ref[0][0])*(-2*vop_ref[0][0] + 1) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2664 +( r*sprite_ref[0][1] - virtual_ref[0][1])*(-2*vop_ref[0][1] + 1) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2665 +2*w2*r*sprite_ref[0][0] - 16*w2); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2666 s->sprite_offset[1][1]= ((-r*sprite_ref[0][1] + virtual_ref[0][1])*(-2*vop_ref[0][0] + 1) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2667 +(-r*sprite_ref[0][0] + virtual_ref[0][0])*(-2*vop_ref[0][1] + 1) |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2668 +2*w2*r*sprite_ref[0][1] - 16*w2); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2669 s->sprite_delta[0][0][0]= (-r*sprite_ref[0][0] + virtual_ref[0][0]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2670 s->sprite_delta[0][0][1]= ( r*sprite_ref[0][1] - virtual_ref[0][1]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2671 s->sprite_delta[0][1][0]= (-r*sprite_ref[0][1] + virtual_ref[0][1]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2672 s->sprite_delta[0][1][1]= (-r*sprite_ref[0][0] + virtual_ref[0][0]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2673 s->sprite_delta[1][0][0]= 4*(-r*sprite_ref[0][0] + virtual_ref[0][0]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2674 s->sprite_delta[1][0][1]= 4*( r*sprite_ref[0][1] - virtual_ref[0][1]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2675 s->sprite_delta[1][1][0]= 4*(-r*sprite_ref[0][1] + virtual_ref[0][1]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2676 s->sprite_delta[1][1][1]= 4*(-r*sprite_ref[0][0] + virtual_ref[0][0]); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2677 s->sprite_shift[0][0]= alpha+rho; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2678 s->sprite_shift[0][1]= alpha+rho; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2679 s->sprite_shift[1][0]= alpha+rho+2; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2680 s->sprite_shift[1][1]= alpha+rho+2; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2681 break; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2682 // case 3: |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2683 break; |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2684 } |
255 | 2685 /*printf("%d %d\n", s->sprite_delta[0][0][0], a<<s->sprite_shift[0][0]); |
2686 printf("%d %d\n", s->sprite_delta[0][0][1], 0); | |
2687 printf("%d %d\n", s->sprite_delta[0][1][0], 0); | |
2688 printf("%d %d\n", s->sprite_delta[0][1][1], a<<s->sprite_shift[0][1]); | |
2689 printf("%d %d\n", s->sprite_delta[1][0][0], a<<s->sprite_shift[1][0]); | |
2690 printf("%d %d\n", s->sprite_delta[1][0][1], 0); | |
2691 printf("%d %d\n", s->sprite_delta[1][1][0], 0); | |
2692 printf("%d %d\n", s->sprite_delta[1][1][1], a<<s->sprite_shift[1][1]);*/ | |
2693 /* try to simplify the situation */ | |
2694 if( s->sprite_delta[0][0][0] == a<<s->sprite_shift[0][0] | |
2695 && s->sprite_delta[0][0][1] == 0 | |
2696 && s->sprite_delta[0][1][0] == 0 | |
2697 && s->sprite_delta[0][1][1] == a<<s->sprite_shift[0][1] | |
2698 && s->sprite_delta[1][0][0] == a<<s->sprite_shift[1][0] | |
2699 && s->sprite_delta[1][0][1] == 0 | |
2700 && s->sprite_delta[1][1][0] == 0 | |
2701 && s->sprite_delta[1][1][1] == a<<s->sprite_shift[1][1]) | |
2702 { | |
2703 s->sprite_offset[0][0]>>=s->sprite_shift[0][0]; | |
2704 s->sprite_offset[0][1]>>=s->sprite_shift[0][1]; | |
2705 s->sprite_offset[1][0]>>=s->sprite_shift[1][0]; | |
2706 s->sprite_offset[1][1]>>=s->sprite_shift[1][1]; | |
2707 s->sprite_delta[0][0][0]= a; | |
2708 s->sprite_delta[0][0][1]= 0; | |
2709 s->sprite_delta[0][1][0]= 0; | |
2710 s->sprite_delta[0][1][1]= a; | |
2711 s->sprite_delta[1][0][0]= a; | |
2712 s->sprite_delta[1][0][1]= 0; | |
2713 s->sprite_delta[1][1][0]= 0; | |
2714 s->sprite_delta[1][1][1]= a; | |
2715 s->sprite_shift[0][0]= 0; | |
2716 s->sprite_shift[0][1]= 0; | |
2717 s->sprite_shift[1][0]= 0; | |
2718 s->sprite_shift[1][1]= 0; | |
2719 s->real_sprite_warping_points=1; | |
2720 } | |
2721 else | |
2722 s->real_sprite_warping_points= s->num_sprite_warping_points; | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2723 |
301 | 2724 //printf("%d %d %d %d\n", d[0][0], d[0][1], s->sprite_offset[0][0], s->sprite_offset[0][1]); |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2725 } |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2726 |
0 | 2727 /* decode mpeg4 VOP header */ |
2728 int mpeg4_decode_picture_header(MpegEncContext * s) | |
2729 { | |
2730 int time_incr, startcode, state, v; | |
324 | 2731 int time_increment; |
346
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
2732 int vol_control=-1; |
0 | 2733 |
2734 redo: | |
2735 /* search next start code */ | |
2736 align_get_bits(&s->gb); | |
2737 state = 0xff; | |
2738 for(;;) { | |
2739 v = get_bits(&s->gb, 8); | |
2740 if (state == 0x000001) { | |
2741 state = ((state << 8) | v) & 0xffffff; | |
2742 startcode = state; | |
2743 break; | |
2744 } | |
2745 state = ((state << 8) | v) & 0xffffff; | |
278 | 2746 if( get_bits_count(&s->gb) > s->gb.size*8-32){ |
333 | 2747 if(s->gb.size>50){ |
2748 printf("no VOP startcode found, frame size was=%d\n", s->gb.size); | |
2749 return -1; | |
2750 }else{ | |
2751 printf("frame skip\n"); | |
2752 return FRAME_SKIPED; | |
2753 } | |
263 | 2754 } |
0 | 2755 } |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2756 //printf("startcode %X %d\n", startcode, get_bits_count(&s->gb)); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2757 if (startcode == 0x120) { // Video Object Layer |
262 | 2758 int width, height, vo_ver_id; |
0 | 2759 |
2760 /* vol header */ | |
21 | 2761 skip_bits(&s->gb, 1); /* random access */ |
336 | 2762 s->vo_type= get_bits(&s->gb, 8); |
63 | 2763 if (get_bits1(&s->gb) != 0) { /* is_ol_id */ |
2764 vo_ver_id = get_bits(&s->gb, 4); /* vo_ver_id */ | |
2765 skip_bits(&s->gb, 3); /* vo_priority */ | |
2766 } else { | |
2767 vo_ver_id = 1; | |
2768 } | |
341 | 2769 //printf("vo type:%d\n",s->vo_type); |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2770 s->aspect_ratio_info= get_bits(&s->gb, 4); |
355 | 2771 if(s->aspect_ratio_info == EXTENDED_PAR){ |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2772 skip_bits(&s->gb, 8); //par_width |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2773 skip_bits(&s->gb, 8); // par_height |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2774 } |
281
1fc96b02142e
mpeg4 aspect_ratio_info in AVCodecContext (requested by alex)
michaelni
parents:
278
diff
changeset
|
2775 |
350
6ebbecc10063
- Advanced Intra Coding (AIC) support for H.263+ encoder, just DC by now.
pulento
parents:
346
diff
changeset
|
2776 if ((vol_control=get_bits1(&s->gb))) { /* vol control parameter */ |
336 | 2777 int chroma_format= get_bits(&s->gb, 2); |
2778 if(chroma_format!=1){ | |
2779 printf("illegal chroma format\n"); | |
2780 } | |
2781 s->low_delay= get_bits1(&s->gb); | |
2782 if(get_bits1(&s->gb)){ /* vbv parameters */ | |
2783 printf("vbv parameters not supported\n"); | |
2784 return -1; | |
2785 } | |
2786 }else{ | |
2787 s->low_delay=0; | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2788 } |
336 | 2789 |
63 | 2790 s->shape = get_bits(&s->gb, 2); /* vol shape */ |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
2791 if(s->shape != RECT_SHAPE) printf("only rectangular vol supported\n"); |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2792 if(s->shape == GRAY_SHAPE && vo_ver_id != 1){ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2793 printf("Gray shape not supported\n"); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2794 skip_bits(&s->gb, 4); //video_object_layer_shape_extension |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2795 } |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2796 |
21 | 2797 skip_bits1(&s->gb); /* marker */ |
0 | 2798 |
262 | 2799 s->time_increment_resolution = get_bits(&s->gb, 16); |
2800 s->time_increment_bits = av_log2(s->time_increment_resolution - 1) + 1; | |
63 | 2801 if (s->time_increment_bits < 1) |
2802 s->time_increment_bits = 1; | |
21 | 2803 skip_bits1(&s->gb); /* marker */ |
0 | 2804 |
63 | 2805 if (get_bits1(&s->gb) != 0) { /* fixed_vop_rate */ |
2806 skip_bits(&s->gb, s->time_increment_bits); | |
2807 } | |
0 | 2808 |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2809 if (s->shape != BIN_ONLY_SHAPE) { |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2810 if (s->shape == RECT_SHAPE) { |
63 | 2811 skip_bits1(&s->gb); /* marker */ |
2812 width = get_bits(&s->gb, 13); | |
2813 skip_bits1(&s->gb); /* marker */ | |
2814 height = get_bits(&s->gb, 13); | |
2815 skip_bits1(&s->gb); /* marker */ | |
274
d0c186bcf075
use the width & height from the mpeg4 header ... in the case that its complete
michaelni
parents:
272
diff
changeset
|
2816 if(width && height){ /* they should be non zero but who knows ... */ |
d0c186bcf075
use the width & height from the mpeg4 header ... in the case that its complete
michaelni
parents:
272
diff
changeset
|
2817 s->width = width; |
d0c186bcf075
use the width & height from the mpeg4 header ... in the case that its complete
michaelni
parents:
272
diff
changeset
|
2818 s->height = height; |
331 | 2819 // printf("width/height: %d %d\n", width, height); |
274
d0c186bcf075
use the width & height from the mpeg4 header ... in the case that its complete
michaelni
parents:
272
diff
changeset
|
2820 } |
63 | 2821 } |
2822 | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
2823 if(get_bits1(&s->gb)) printf("interlaced not supported\n"); /* interlaced */ |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2824 if(!get_bits1(&s->gb)) printf("OBMC not supported\n"); /* OBMC Disable */ |
63 | 2825 if (vo_ver_id == 1) { |
2826 s->vol_sprite_usage = get_bits1(&s->gb); /* vol_sprite_usage */ | |
2827 } else { | |
2828 s->vol_sprite_usage = get_bits(&s->gb, 2); /* vol_sprite_usage */ | |
2829 } | |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2830 if(s->vol_sprite_usage==STATIC_SPRITE) printf("Static Sprites not supported\n"); |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2831 if(s->vol_sprite_usage==STATIC_SPRITE || s->vol_sprite_usage==GMC_SPRITE){ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2832 if(s->vol_sprite_usage==STATIC_SPRITE){ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2833 s->sprite_width = get_bits(&s->gb, 13); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2834 skip_bits1(&s->gb); /* marker */ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2835 s->sprite_height= get_bits(&s->gb, 13); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2836 skip_bits1(&s->gb); /* marker */ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2837 s->sprite_left = get_bits(&s->gb, 13); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2838 skip_bits1(&s->gb); /* marker */ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2839 s->sprite_top = get_bits(&s->gb, 13); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2840 skip_bits1(&s->gb); /* marker */ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2841 } |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
2842 s->num_sprite_warping_points= get_bits(&s->gb, 6); |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2843 s->sprite_warping_accuracy = get_bits(&s->gb, 2); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2844 s->sprite_brightness_change= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2845 if(s->vol_sprite_usage==STATIC_SPRITE) |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2846 s->low_latency_sprite= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2847 } |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2848 // FIXME sadct disable bit if verid!=1 && shape not rect |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2849 |
63 | 2850 if (get_bits1(&s->gb) == 1) { /* not_8_bit */ |
2851 s->quant_precision = get_bits(&s->gb, 4); /* quant_precision */ | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
2852 if(get_bits(&s->gb, 4)!=8) printf("N-bit not supported\n"); /* bits_per_pixel */ |
290 | 2853 if(s->quant_precision!=5) printf("quant precission %d\n", s->quant_precision); |
63 | 2854 } else { |
2855 s->quant_precision = 5; | |
2856 } | |
2857 | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2858 // FIXME a bunch of grayscale shape things |
312 | 2859 |
2860 if(get_bits1(&s->gb)){ /* vol_quant_type */ | |
2861 int i, j, v; | |
2862 /* load default matrixes */ | |
2863 for(i=0; i<64; i++){ | |
2864 v= ff_mpeg4_default_intra_matrix[i]; | |
2865 s->intra_matrix[i]= v; | |
2866 s->chroma_intra_matrix[i]= v; | |
2867 | |
2868 v= ff_mpeg4_default_non_intra_matrix[i]; | |
344 | 2869 s->inter_matrix[i]= v; |
2870 s->chroma_inter_matrix[i]= v; | |
312 | 2871 } |
2872 | |
2873 /* load custom intra matrix */ | |
2874 if(get_bits1(&s->gb)){ | |
2875 for(i=0; i<64; i++){ | |
2876 v= get_bits(&s->gb, 8); | |
2877 if(v==0) break; | |
2878 | |
2879 j= zigzag_direct[i]; | |
2880 s->intra_matrix[j]= v; | |
2881 s->chroma_intra_matrix[j]= v; | |
2882 } | |
2883 } | |
2884 | |
2885 /* load custom non intra matrix */ | |
2886 if(get_bits1(&s->gb)){ | |
2887 for(i=0; i<64; i++){ | |
2888 v= get_bits(&s->gb, 8); | |
2889 if(v==0) break; | |
2890 | |
2891 j= zigzag_direct[i]; | |
344 | 2892 s->inter_matrix[j]= v; |
2893 s->chroma_inter_matrix[j]= v; | |
312 | 2894 } |
2895 | |
2896 /* replicate last value */ | |
2897 for(; i<64; i++){ | |
2898 j= zigzag_direct[i]; | |
344 | 2899 s->inter_matrix[j]= v; |
2900 s->chroma_inter_matrix[j]= v; | |
312 | 2901 } |
2902 } | |
2903 | |
325 | 2904 s->dct_unquantize= s->dct_unquantize_mpeg2; |
312 | 2905 |
2906 // FIXME a bunch of grayscale shape things | |
2907 }else | |
2908 s->dct_unquantize= s->dct_unquantize_h263; | |
2909 | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2910 if(vo_ver_id != 1) |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2911 s->quarter_sample= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2912 else s->quarter_sample=0; |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
2913 |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
2914 if(!get_bits1(&s->gb)) printf("Complexity estimation not supported\n"); |
290 | 2915 |
2916 s->resync_marker= !get_bits1(&s->gb); /* resync_marker_disabled */ | |
2917 | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2918 s->data_partioning= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2919 if(s->data_partioning){ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2920 printf("data partitioning not supported\n"); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2921 skip_bits1(&s->gb); // reversible vlc |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2922 } |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2923 |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2924 if(vo_ver_id != 1) { |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2925 s->new_pred= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2926 if(s->new_pred){ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2927 printf("new pred not supported\n"); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2928 skip_bits(&s->gb, 2); /* requested upstream message type */ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2929 skip_bits1(&s->gb); /* newpred segment type */ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2930 } |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2931 s->reduced_res_vop= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2932 if(s->reduced_res_vop) printf("reduced resolution VOP not supported\n"); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2933 } |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2934 else{ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2935 s->new_pred=0; |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2936 s->reduced_res_vop= 0; |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2937 } |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2938 |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2939 s->scalability= get_bits1(&s->gb); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2940 if (s->scalability) { |
285 | 2941 printf("scalability not supported\n"); |
63 | 2942 } |
2943 } | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2944 //printf("end Data %X %d\n", show_bits(&s->gb, 32), get_bits_count(&s->gb)&0x7); |
0 | 2945 goto redo; |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2946 } else if (startcode == 0x1b2) { //userdata |
255 | 2947 char buf[256]; |
2948 int i; | |
2949 int e; | |
2950 int ver, build; | |
2951 | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2952 //printf("user Data %X\n", show_bits(&s->gb, 32)); |
255 | 2953 buf[0]= show_bits(&s->gb, 8); |
2954 for(i=1; i<256; i++){ | |
2955 buf[i]= show_bits(&s->gb, 16)&0xFF; | |
2956 if(buf[i]==0) break; | |
2957 skip_bits(&s->gb, 8); | |
2958 } | |
2959 buf[255]=0; | |
2960 e=sscanf(buf, "DivX%dBuild%d", &ver, &build); | |
333 | 2961 if(e!=2) |
2962 e=sscanf(buf, "DivX%db%d", &ver, &build); | |
255 | 2963 if(e==2){ |
2964 s->divx_version= ver; | |
2965 s->divx_build= build; | |
2966 if(s->picture_number==0){ | |
2967 printf("This file was encoded with DivX%d Build%d\n", ver, build); | |
333 | 2968 if(ver==500 && build==413){ |
255 | 2969 printf("WARNING: this version of DivX is not MPEG4 compatible, trying to workaround these bugs...\n"); |
333 | 2970 #if 0 |
255 | 2971 }else{ |
2972 printf("hmm, i havnt seen that version of divx yet, lets assume they fixed these bugs ...\n" | |
2973 "using mpeg4 decoder, if it fails contact the developers (of ffmpeg)\n"); | |
333 | 2974 #endif |
255 | 2975 } |
2976 } | |
2977 } | |
2978 //printf("User Data: %s\n", buf); | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2979 goto redo; |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
2980 } else if (startcode != 0x1b6) { //VOP |
0 | 2981 goto redo; |
2982 } | |
2983 | |
355 | 2984 s->pict_type = get_bits(&s->gb, 2) + I_TYPE; /* pict type: I = 0 , P = 1 */ |
324 | 2985 // printf("pic: %d, qpel:%d\n", s->pict_type, s->quarter_sample); |
262 | 2986 time_incr=0; |
21 | 2987 while (get_bits1(&s->gb) != 0) |
0 | 2988 time_incr++; |
2989 | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
2990 check_marker(&s->gb, "before time_increment"); |
324 | 2991 time_increment= get_bits(&s->gb, s->time_increment_bits); |
2992 //printf(" type:%d incr:%d increment:%d\n", s->pict_type, time_incr, time_increment); | |
262 | 2993 if(s->pict_type!=B_TYPE){ |
324 | 2994 s->last_time_base= s->time_base; |
262 | 2995 s->time_base+= time_incr; |
324 | 2996 s->time= s->time_base*s->time_increment_resolution + time_increment; |
2997 s->pp_time= s->time - s->last_non_b_time; | |
2998 s->last_non_b_time= s->time; | |
262 | 2999 }else{ |
324 | 3000 s->time= (s->last_time_base + time_incr)*s->time_increment_resolution + time_increment; |
3001 s->bp_time= s->last_non_b_time - s->time; | |
335 | 3002 if(s->pp_time <=s->bp_time){ |
3003 // printf("messed up order, seeking?, skiping current b frame\n"); | |
3004 return FRAME_SKIPED; | |
3005 } | |
262 | 3006 } |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3007 |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3008 if(check_marker(&s->gb, "before vop_coded")==0 && s->picture_number==0){ |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3009 printf("hmm, seems the headers arnt complete, trying to guess time_increment_bits\n"); |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3010 for(s->time_increment_bits++ ;s->time_increment_bits<16; s->time_increment_bits++){ |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3011 if(get_bits1(&s->gb)) break; |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3012 } |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3013 printf("my guess is %d bits ;)\n",s->time_increment_bits); |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3014 } |
0 | 3015 /* vop coded */ |
21 | 3016 if (get_bits1(&s->gb) != 1) |
63 | 3017 goto redo; |
262 | 3018 //printf("time %d %d %d || %d %d %d\n", s->time_increment_bits, s->time_increment, s->time_base, |
3019 //s->time, s->last_non_b_time[0], s->last_non_b_time[1]); | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3020 if (s->shape != BIN_ONLY_SHAPE && ( s->pict_type == P_TYPE |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3021 || (s->pict_type == S_TYPE && s->vol_sprite_usage==GMC_SPRITE))) { |
0 | 3022 /* rounding type for motion estimation */ |
21 | 3023 s->no_rounding = get_bits1(&s->gb); |
63 | 3024 } else { |
3025 s->no_rounding = 0; | |
0 | 3026 } |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3027 //FIXME reduced res stuff |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3028 |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3029 if (s->shape != RECT_SHAPE) { |
63 | 3030 if (s->vol_sprite_usage != 1 || s->pict_type != I_TYPE) { |
3031 int width, height, hor_spat_ref, ver_spat_ref; | |
3032 | |
3033 width = get_bits(&s->gb, 13); | |
3034 skip_bits1(&s->gb); /* marker */ | |
3035 height = get_bits(&s->gb, 13); | |
3036 skip_bits1(&s->gb); /* marker */ | |
3037 hor_spat_ref = get_bits(&s->gb, 13); /* hor_spat_ref */ | |
3038 skip_bits1(&s->gb); /* marker */ | |
3039 ver_spat_ref = get_bits(&s->gb, 13); /* ver_spat_ref */ | |
3040 } | |
3041 skip_bits1(&s->gb); /* change_CR_disable */ | |
3042 | |
3043 if (get_bits1(&s->gb) != 0) { | |
3044 skip_bits(&s->gb, 8); /* constant_alpha_value */ | |
3045 } | |
3046 } | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3047 //FIXME complexity estimation stuff |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3048 |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3049 if (s->shape != BIN_ONLY_SHAPE) { |
290 | 3050 int t; |
3051 t=get_bits(&s->gb, 3); /* intra dc VLC threshold */ | |
3052 //printf("threshold %d\n", t); | |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3053 //FIXME interlaced specific bits |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3054 } |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
3055 |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3056 if(s->pict_type == S_TYPE && (s->vol_sprite_usage==STATIC_SPRITE || s->vol_sprite_usage==GMC_SPRITE)){ |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
3057 if(s->num_sprite_warping_points){ |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
3058 mpeg4_decode_sprite_trajectory(s); |
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
3059 } |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3060 if(s->sprite_brightness_change) printf("sprite_brightness_change not supported\n"); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3061 if(s->vol_sprite_usage==STATIC_SPRITE) printf("static sprite not supported\n"); |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3062 } |
254
b4fed8b24e3a
gmc bitstream decoding support (the real motion compensation isnt implemnted yet)
michaelni
parents:
253
diff
changeset
|
3063 |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3064 if (s->shape != BIN_ONLY_SHAPE) { |
63 | 3065 /* note: we do not use quant_precision to avoid problem if no |
3066 MPEG4 vol header as it is found on some old opendivx | |
3067 movies */ | |
3068 s->qscale = get_bits(&s->gb, 5); | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3069 if(s->qscale==0){ |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3070 printf("Error, header damaged or not MPEG4 header (qscale=0)\n"); |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3071 return -1; // makes no sense to continue, as there is nothing left from the image then |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3072 } |
63 | 3073 |
3074 if (s->pict_type != I_TYPE) { | |
3075 s->f_code = get_bits(&s->gb, 3); /* fcode_for */ | |
264
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3076 if(s->f_code==0){ |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3077 printf("Error, header damaged or not MPEG4 header (f_code=0)\n"); |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3078 return -1; // makes no sense to continue, as the MV decoding will break very quickly |
28c5c62b1c4c
support decoding (with mplayer) of 3 .mp4 files from mphq
michaelni
parents:
263
diff
changeset
|
3079 } |
63 | 3080 } |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3081 if (s->pict_type == B_TYPE) { |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3082 s->b_code = get_bits(&s->gb, 3); |
262 | 3083 //printf("b-code %d\n", s->b_code); |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3084 } |
331 | 3085 //printf("quant:%d fcode:%d bcode:%d type:%d\n", s->qscale, s->f_code, s->b_code, s->pict_type); |
253
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3086 if(!s->scalability){ |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3087 if (s->shape!=RECT_SHAPE && s->pict_type!=I_TYPE) { |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3088 skip_bits1(&s->gb); // vop shape coding type |
4448dd55d415
parsing more of the mpeg4 header & print some "not supported" stuff
michaelni
parents:
250
diff
changeset
|
3089 } |
63 | 3090 } |
3091 } | |
346
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3092 /* detect buggy encoders which dont set the low_delay flag (divx4/xvid/opendivx)*/ |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3093 // note we cannot detect divx5 without b-frames easyly (allthough its buggy too) |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3094 if(s->vo_type==0 && vol_control==0 && s->divx_version==0){ |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3095 if(s->picture_number==0) |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3096 printf("looks like this file was encoded with (divx4/(old)xvid/opendivx) -> forcing low_delay flag\n"); |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3097 s->low_delay=1; |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3098 } |
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3099 |
255 | 3100 s->picture_number++; // better than pic number==0 allways ;) |
333 | 3101 //printf("done\n"); |
346
c2f789fe4945
detecting xvid/divx4/opendivx and set low_delay flag
michaelni
parents:
344
diff
changeset
|
3102 |
63 | 3103 return 0; |
0 | 3104 } |
3105 | |
3106 /* don't understand why they choose a different header ! */ | |
3107 int intel_h263_decode_picture_header(MpegEncContext *s) | |
3108 { | |
3109 int format; | |
3110 | |
3111 /* picture header */ | |
355 | 3112 if (get_bits(&s->gb, 22) != 0x20) { |
3113 fprintf(stderr, "Bad picture start code\n"); | |
0 | 3114 return -1; |
355 | 3115 } |
3116 s->picture_number = get_bits(&s->gb, 8); /* picture timestamp */ | |
0 | 3117 |
355 | 3118 if (get_bits1(&s->gb) != 1) { |
3119 fprintf(stderr, "Bad marker\n"); | |
0 | 3120 return -1; /* marker */ |
355 | 3121 } |
3122 if (get_bits1(&s->gb) != 0) { | |
3123 fprintf(stderr, "Bad H263 id\n"); | |
0 | 3124 return -1; /* h263 id */ |
355 | 3125 } |
21 | 3126 skip_bits1(&s->gb); /* split screen off */ |
3127 skip_bits1(&s->gb); /* camera off */ | |
3128 skip_bits1(&s->gb); /* freeze picture release off */ | |
0 | 3129 |
3130 format = get_bits(&s->gb, 3); | |
355 | 3131 if (format != 7) { |
3132 fprintf(stderr, "Intel H263 free format not supported\n"); | |
0 | 3133 return -1; |
355 | 3134 } |
0 | 3135 s->h263_plus = 0; |
3136 | |
21 | 3137 s->pict_type = I_TYPE + get_bits1(&s->gb); |
0 | 3138 |
21 | 3139 s->unrestricted_mv = get_bits1(&s->gb); |
0 | 3140 s->h263_long_vectors = s->unrestricted_mv; |
3141 | |
355 | 3142 if (get_bits1(&s->gb) != 0) { |
3143 fprintf(stderr, "SAC not supported\n"); | |
0 | 3144 return -1; /* SAC: off */ |
355 | 3145 } |
3146 if (get_bits1(&s->gb) != 0) { | |
3147 fprintf(stderr, "Advanced Prediction Mode not supported\n"); | |
0 | 3148 return -1; /* advanced prediction mode: off */ |
355 | 3149 } |
3150 if (get_bits1(&s->gb) != 0) { | |
3151 fprintf(stderr, "PB frame mode no supported\n"); | |
3152 return -1; /* PB frame mode */ | |
3153 } | |
0 | 3154 |
3155 /* skip unknown header garbage */ | |
21 | 3156 skip_bits(&s->gb, 41); |
0 | 3157 |
3158 s->qscale = get_bits(&s->gb, 5); | |
21 | 3159 skip_bits1(&s->gb); /* Continuous Presence Multipoint mode: off */ |
0 | 3160 |
3161 /* PEI */ | |
21 | 3162 while (get_bits1(&s->gb) != 0) { |
3163 skip_bits(&s->gb, 8); | |
0 | 3164 } |
3165 s->f_code = 1; | |
3166 return 0; | |
3167 } | |
144 | 3168 |