00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00027 #include "avcodec.h"
00028 #include "bytestream.h"
00029 #include "get_bits.h"
00030 #include "dsputil.h"
00031
00032 #define TM2_ESCAPE 0x80000000
00033 #define TM2_DELTAS 64
00034
00035 enum TM2_STREAMS{ TM2_C_HI = 0, TM2_C_LO, TM2_L_HI, TM2_L_LO,
00036 TM2_UPD, TM2_MOT, TM2_TYPE, TM2_NUM_STREAMS};
00037
00038 enum TM2_BLOCKS{ TM2_HI_RES = 0, TM2_MED_RES, TM2_LOW_RES, TM2_NULL_RES,
00039 TM2_UPDATE, TM2_STILL, TM2_MOTION};
00040
00041 typedef struct TM2Context{
00042 AVCodecContext *avctx;
00043 AVFrame pic;
00044
00045 GetBitContext gb;
00046 DSPContext dsp;
00047
00048 uint8_t *buffer;
00049 int buffer_size;
00050
00051
00052 int *tokens[TM2_NUM_STREAMS];
00053 int tok_lens[TM2_NUM_STREAMS];
00054 int tok_ptrs[TM2_NUM_STREAMS];
00055 int deltas[TM2_NUM_STREAMS][TM2_DELTAS];
00056
00057 int D[4];
00058 int CD[4];
00059 int *last;
00060 int *clast;
00061
00062
00063 int *Y1_base, *U1_base, *V1_base, *Y2_base, *U2_base, *V2_base;
00064 int *Y1, *U1, *V1, *Y2, *U2, *V2;
00065 int y_stride, uv_stride;
00066 int cur;
00067 } TM2Context;
00068
00072 typedef struct TM2Codes{
00073 VLC vlc;
00074 int bits;
00075 int *recode;
00076 int length;
00077 } TM2Codes;
00078
00082 typedef struct TM2Huff{
00083 int val_bits;
00084 int max_bits;
00085 int min_bits;
00086 int nodes;
00087 int num;
00088 int max_num;
00089 int *nums;
00090 uint32_t *bits;
00091 int *lens;
00092 } TM2Huff;
00093
00094 static int tm2_read_tree(TM2Context *ctx, uint32_t prefix, int length, TM2Huff *huff)
00095 {
00096 if(length > huff->max_bits) {
00097 av_log(ctx->avctx, AV_LOG_ERROR, "Tree exceeded its given depth (%i)\n", huff->max_bits);
00098 return -1;
00099 }
00100
00101 if(!get_bits1(&ctx->gb)) {
00102 if (length == 0) {
00103 length = 1;
00104 }
00105 if(huff->num >= huff->max_num) {
00106 av_log(ctx->avctx, AV_LOG_DEBUG, "Too many literals\n");
00107 return -1;
00108 }
00109 huff->nums[huff->num] = get_bits_long(&ctx->gb, huff->val_bits);
00110 huff->bits[huff->num] = prefix;
00111 huff->lens[huff->num] = length;
00112 huff->num++;
00113 return 0;
00114 } else {
00115 if(tm2_read_tree(ctx, prefix << 1, length + 1, huff) == -1)
00116 return -1;
00117 if(tm2_read_tree(ctx, (prefix << 1) | 1, length + 1, huff) == -1)
00118 return -1;
00119 }
00120 return 0;
00121 }
00122
00123 static int tm2_build_huff_table(TM2Context *ctx, TM2Codes *code)
00124 {
00125 TM2Huff huff;
00126 int res = 0;
00127
00128 huff.val_bits = get_bits(&ctx->gb, 5);
00129 huff.max_bits = get_bits(&ctx->gb, 5);
00130 huff.min_bits = get_bits(&ctx->gb, 5);
00131 huff.nodes = get_bits_long(&ctx->gb, 17);
00132 huff.num = 0;
00133
00134
00135 if((huff.val_bits < 1) || (huff.val_bits > 32) ||
00136 (huff.max_bits < 0) || (huff.max_bits > 25)) {
00137 av_log(ctx->avctx, AV_LOG_ERROR, "Incorrect tree parameters - literal length: %i, max code length: %i\n",
00138 huff.val_bits, huff.max_bits);
00139 return -1;
00140 }
00141 if((huff.nodes <= 0) || (huff.nodes > 0x10000)) {
00142 av_log(ctx->avctx, AV_LOG_ERROR, "Incorrect number of Huffman tree nodes: %i\n", huff.nodes);
00143 return -1;
00144 }
00145
00146 if(huff.max_bits == 0)
00147 huff.max_bits = 1;
00148
00149
00150 huff.max_num = (huff.nodes + 1) >> 1;
00151 huff.nums = av_mallocz(huff.max_num * sizeof(int));
00152 huff.bits = av_mallocz(huff.max_num * sizeof(uint32_t));
00153 huff.lens = av_mallocz(huff.max_num * sizeof(int));
00154
00155 if(tm2_read_tree(ctx, 0, 0, &huff) == -1)
00156 res = -1;
00157
00158 if(huff.num != huff.max_num) {
00159 av_log(ctx->avctx, AV_LOG_ERROR, "Got less codes than expected: %i of %i\n",
00160 huff.num, huff.max_num);
00161 res = -1;
00162 }
00163
00164
00165 if(res != -1) {
00166 int i;
00167
00168 res = init_vlc(&code->vlc, huff.max_bits, huff.max_num,
00169 huff.lens, sizeof(int), sizeof(int),
00170 huff.bits, sizeof(uint32_t), sizeof(uint32_t), 0);
00171 if(res < 0) {
00172 av_log(ctx->avctx, AV_LOG_ERROR, "Cannot build VLC table\n");
00173 res = -1;
00174 } else
00175 res = 0;
00176 if(res != -1) {
00177 code->bits = huff.max_bits;
00178 code->length = huff.max_num;
00179 code->recode = av_malloc(code->length * sizeof(int));
00180 for(i = 0; i < code->length; i++)
00181 code->recode[i] = huff.nums[i];
00182 }
00183 }
00184
00185 av_free(huff.nums);
00186 av_free(huff.bits);
00187 av_free(huff.lens);
00188
00189 return res;
00190 }
00191
00192 static void tm2_free_codes(TM2Codes *code)
00193 {
00194 av_free(code->recode);
00195 if(code->vlc.table)
00196 ff_free_vlc(&code->vlc);
00197 }
00198
00199 static inline int tm2_get_token(GetBitContext *gb, TM2Codes *code)
00200 {
00201 int val;
00202 val = get_vlc2(gb, code->vlc.table, code->bits, 1);
00203 if(val<0)
00204 return -1;
00205 return code->recode[val];
00206 }
00207
00208 static inline int tm2_read_header(TM2Context *ctx, const uint8_t *buf)
00209 {
00210 uint32_t magic;
00211
00212 magic = AV_RL32(buf);
00213 buf += 4;
00214
00215 if(magic == 0x00000100) {
00216 av_log_missing_feature(ctx->avctx, "TM2 old header", 1);
00217 return 40;
00218 } else if(magic == 0x00000101) {
00219 return 40;
00220 } else {
00221 av_log (ctx->avctx, AV_LOG_ERROR, "Not a TM2 header: 0x%08X\n", magic);
00222 return -1;
00223 }
00224 }
00225
00226 static int tm2_read_deltas(TM2Context *ctx, int stream_id) {
00227 int d, mb;
00228 int i, v;
00229
00230 d = get_bits(&ctx->gb, 9);
00231 mb = get_bits(&ctx->gb, 5);
00232
00233 if((d < 1) || (d > TM2_DELTAS) || (mb < 1) || (mb > 32)) {
00234 av_log(ctx->avctx, AV_LOG_ERROR, "Incorrect delta table: %i deltas x %i bits\n", d, mb);
00235 return -1;
00236 }
00237
00238 for(i = 0; i < d; i++) {
00239 v = get_bits_long(&ctx->gb, mb);
00240 if(v & (1 << (mb - 1)))
00241 ctx->deltas[stream_id][i] = v - (1 << mb);
00242 else
00243 ctx->deltas[stream_id][i] = v;
00244 }
00245 for(; i < TM2_DELTAS; i++)
00246 ctx->deltas[stream_id][i] = 0;
00247
00248 return 0;
00249 }
00250
00251 static int tm2_read_stream(TM2Context *ctx, const uint8_t *buf, int stream_id, int buf_size)
00252 {
00253 int i;
00254 int skip = 0;
00255 int len, toks, pos;
00256 TM2Codes codes;
00257 GetByteContext gb;
00258
00259 if (buf_size < 4) {
00260 av_log(ctx->avctx, AV_LOG_ERROR, "not enough space for len left\n");
00261 return AVERROR_INVALIDDATA;
00262 }
00263
00264
00265 bytestream2_init(&gb, buf, buf_size);
00266 len = bytestream2_get_be32(&gb);
00267 skip = len * 4 + 4;
00268
00269 if(len == 0)
00270 return 4;
00271
00272 if (len >= INT_MAX/4-1 || len < 0 || skip > buf_size) {
00273 av_log(ctx->avctx, AV_LOG_ERROR, "invalid stream size\n");
00274 return AVERROR_INVALIDDATA;
00275 }
00276
00277 toks = bytestream2_get_be32(&gb);
00278 if(toks & 1) {
00279 len = bytestream2_get_be32(&gb);
00280 if(len == TM2_ESCAPE) {
00281 len = bytestream2_get_be32(&gb);
00282 }
00283 if(len > 0) {
00284 pos = bytestream2_tell(&gb);
00285 if (skip <= pos)
00286 return AVERROR_INVALIDDATA;
00287 init_get_bits(&ctx->gb, buf + pos, (skip - pos) * 8);
00288 if(tm2_read_deltas(ctx, stream_id) == -1)
00289 return AVERROR_INVALIDDATA;
00290 bytestream2_skip(&gb, ((get_bits_count(&ctx->gb) + 31) >> 5) << 2);
00291 }
00292 }
00293
00294 len = bytestream2_get_be32(&gb);
00295 if(len == TM2_ESCAPE) {
00296 bytestream2_skip(&gb, 8);
00297 } else {
00298 bytestream2_skip(&gb, 4);
00299 }
00300
00301 pos = bytestream2_tell(&gb);
00302 if (skip <= pos)
00303 return AVERROR_INVALIDDATA;
00304 init_get_bits(&ctx->gb, buf + pos, (skip - pos) * 8);
00305 if(tm2_build_huff_table(ctx, &codes) == -1)
00306 return AVERROR_INVALIDDATA;
00307 bytestream2_skip(&gb, ((get_bits_count(&ctx->gb) + 31) >> 5) << 2);
00308
00309 toks >>= 1;
00310
00311 if((toks < 0) || (toks > 0xFFFFFF)){
00312 av_log(ctx->avctx, AV_LOG_ERROR, "Incorrect number of tokens: %i\n", toks);
00313 tm2_free_codes(&codes);
00314 return AVERROR_INVALIDDATA;
00315 }
00316 ctx->tokens[stream_id] = av_realloc(ctx->tokens[stream_id], toks * sizeof(int));
00317 ctx->tok_lens[stream_id] = toks;
00318 len = bytestream2_get_be32(&gb);
00319 if(len > 0) {
00320 pos = bytestream2_tell(&gb);
00321 if (skip <= pos)
00322 return AVERROR_INVALIDDATA;
00323 init_get_bits(&ctx->gb, buf + pos, (skip - pos) * 8);
00324 for(i = 0; i < toks; i++) {
00325 if (get_bits_left(&ctx->gb) <= 0) {
00326 av_log(ctx->avctx, AV_LOG_ERROR, "Incorrect number of tokens: %i\n", toks);
00327 return AVERROR_INVALIDDATA;
00328 }
00329 ctx->tokens[stream_id][i] = tm2_get_token(&ctx->gb, &codes);
00330 if (stream_id <= TM2_MOT && ctx->tokens[stream_id][i] >= TM2_DELTAS || ctx->tokens[stream_id][i]<0) {
00331 av_log(ctx->avctx, AV_LOG_ERROR, "Invalid delta token index %d for type %d, n=%d\n",
00332 ctx->tokens[stream_id][i], stream_id, i);
00333 return AVERROR_INVALIDDATA;
00334 }
00335 }
00336 } else {
00337 for(i = 0; i < toks; i++) {
00338 ctx->tokens[stream_id][i] = codes.recode[0];
00339 if (stream_id <= TM2_MOT && ctx->tokens[stream_id][i] >= TM2_DELTAS) {
00340 av_log(ctx->avctx, AV_LOG_ERROR, "Invalid delta token index %d for type %d, n=%d\n",
00341 ctx->tokens[stream_id][i], stream_id, i);
00342 return AVERROR_INVALIDDATA;
00343 }
00344 }
00345 }
00346 tm2_free_codes(&codes);
00347
00348 return skip;
00349 }
00350
00351 static inline int GET_TOK(TM2Context *ctx,int type) {
00352 if(ctx->tok_ptrs[type] >= ctx->tok_lens[type]) {
00353 av_log(ctx->avctx, AV_LOG_ERROR, "Read token from stream %i out of bounds (%i>=%i)\n", type, ctx->tok_ptrs[type], ctx->tok_lens[type]);
00354 return 0;
00355 }
00356 if(type <= TM2_MOT) {
00357 if (ctx->tokens[type][ctx->tok_ptrs[type]] >= TM2_DELTAS) {
00358 av_log(ctx->avctx, AV_LOG_ERROR, "token %d is too large\n", ctx->tokens[type][ctx->tok_ptrs[type]]);
00359 return 0;
00360 }
00361 return ctx->deltas[type][ctx->tokens[type][ctx->tok_ptrs[type]++]];
00362 }
00363 return ctx->tokens[type][ctx->tok_ptrs[type]++];
00364 }
00365
00366
00367
00368
00369 #define TM2_INIT_POINTERS() \
00370 int *last, *clast; \
00371 int *Y, *U, *V;\
00372 int Ystride, Ustride, Vstride;\
00373 \
00374 Ystride = ctx->y_stride;\
00375 Vstride = ctx->uv_stride;\
00376 Ustride = ctx->uv_stride;\
00377 Y = (ctx->cur?ctx->Y2:ctx->Y1) + by * 4 * Ystride + bx * 4;\
00378 V = (ctx->cur?ctx->V2:ctx->V1) + by * 2 * Vstride + bx * 2;\
00379 U = (ctx->cur?ctx->U2:ctx->U1) + by * 2 * Ustride + bx * 2;\
00380 last = ctx->last + bx * 4;\
00381 clast = ctx->clast + bx * 4;
00382
00383 #define TM2_INIT_POINTERS_2() \
00384 int *Yo, *Uo, *Vo;\
00385 int oYstride, oUstride, oVstride;\
00386 \
00387 TM2_INIT_POINTERS();\
00388 oYstride = Ystride;\
00389 oVstride = Vstride;\
00390 oUstride = Ustride;\
00391 Yo = (ctx->cur?ctx->Y1:ctx->Y2) + by * 4 * oYstride + bx * 4;\
00392 Vo = (ctx->cur?ctx->V1:ctx->V2) + by * 2 * oVstride + bx * 2;\
00393 Uo = (ctx->cur?ctx->U1:ctx->U2) + by * 2 * oUstride + bx * 2;
00394
00395
00396 #define TM2_RECALC_BLOCK(CHR, stride, last, CD) {\
00397 CD[0] = CHR[1] - last[1];\
00398 CD[1] = (int)CHR[stride + 1] - (int)CHR[1];\
00399 last[0] = (int)CHR[stride + 0];\
00400 last[1] = (int)CHR[stride + 1];}
00401
00402
00403 static inline void tm2_apply_deltas(TM2Context *ctx, int* Y, int stride, int *deltas, int *last)
00404 {
00405 int ct, d;
00406 int i, j;
00407
00408 for(j = 0; j < 4; j++){
00409 ct = ctx->D[j];
00410 for(i = 0; i < 4; i++){
00411 d = deltas[i + j * 4];
00412 ct += d;
00413 last[i] += ct;
00414 Y[i] = av_clip_uint8(last[i]);
00415 }
00416 Y += stride;
00417 ctx->D[j] = ct;
00418 }
00419 }
00420
00421 static inline void tm2_high_chroma(int *data, int stride, int *last, int *CD, int *deltas)
00422 {
00423 int i, j;
00424 for(j = 0; j < 2; j++){
00425 for(i = 0; i < 2; i++){
00426 CD[j] += deltas[i + j * 2];
00427 last[i] += CD[j];
00428 data[i] = last[i];
00429 }
00430 data += stride;
00431 }
00432 }
00433
00434 static inline void tm2_low_chroma(int *data, int stride, int *clast, int *CD, int *deltas, int bx)
00435 {
00436 int t;
00437 int l;
00438 int prev;
00439
00440 if(bx > 0)
00441 prev = clast[-3];
00442 else
00443 prev = 0;
00444 t = (CD[0] + CD[1]) >> 1;
00445 l = (prev - CD[0] - CD[1] + clast[1]) >> 1;
00446 CD[1] = CD[0] + CD[1] - t;
00447 CD[0] = t;
00448 clast[0] = l;
00449
00450 tm2_high_chroma(data, stride, clast, CD, deltas);
00451 }
00452
00453 static inline void tm2_hi_res_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00454 {
00455 int i;
00456 int deltas[16];
00457 TM2_INIT_POINTERS();
00458
00459
00460 for(i = 0; i < 4; i++) {
00461 deltas[i] = GET_TOK(ctx, TM2_C_HI);
00462 deltas[i + 4] = GET_TOK(ctx, TM2_C_HI);
00463 }
00464 tm2_high_chroma(U, Ustride, clast, ctx->CD, deltas);
00465 tm2_high_chroma(V, Vstride, clast + 2, ctx->CD + 2, deltas + 4);
00466
00467
00468 for(i = 0; i < 16; i++)
00469 deltas[i] = GET_TOK(ctx, TM2_L_HI);
00470
00471 tm2_apply_deltas(ctx, Y, Ystride, deltas, last);
00472 }
00473
00474 static inline void tm2_med_res_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00475 {
00476 int i;
00477 int deltas[16];
00478 TM2_INIT_POINTERS();
00479
00480
00481 deltas[0] = GET_TOK(ctx, TM2_C_LO);
00482 deltas[1] = deltas[2] = deltas[3] = 0;
00483 tm2_low_chroma(U, Ustride, clast, ctx->CD, deltas, bx);
00484
00485 deltas[0] = GET_TOK(ctx, TM2_C_LO);
00486 deltas[1] = deltas[2] = deltas[3] = 0;
00487 tm2_low_chroma(V, Vstride, clast + 2, ctx->CD + 2, deltas, bx);
00488
00489
00490 for(i = 0; i < 16; i++)
00491 deltas[i] = GET_TOK(ctx, TM2_L_HI);
00492
00493 tm2_apply_deltas(ctx, Y, Ystride, deltas, last);
00494 }
00495
00496 static inline void tm2_low_res_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00497 {
00498 int i;
00499 int t1, t2;
00500 int deltas[16];
00501 TM2_INIT_POINTERS();
00502
00503
00504 deltas[0] = GET_TOK(ctx, TM2_C_LO);
00505 deltas[1] = deltas[2] = deltas[3] = 0;
00506 tm2_low_chroma(U, Ustride, clast, ctx->CD, deltas, bx);
00507
00508 deltas[0] = GET_TOK(ctx, TM2_C_LO);
00509 deltas[1] = deltas[2] = deltas[3] = 0;
00510 tm2_low_chroma(V, Vstride, clast + 2, ctx->CD + 2, deltas, bx);
00511
00512
00513 for(i = 0; i < 16; i++)
00514 deltas[i] = 0;
00515
00516 deltas[ 0] = GET_TOK(ctx, TM2_L_LO);
00517 deltas[ 2] = GET_TOK(ctx, TM2_L_LO);
00518 deltas[ 8] = GET_TOK(ctx, TM2_L_LO);
00519 deltas[10] = GET_TOK(ctx, TM2_L_LO);
00520
00521 if(bx > 0)
00522 last[0] = (last[-1] - ctx->D[0] - ctx->D[1] - ctx->D[2] - ctx->D[3] + last[1]) >> 1;
00523 else
00524 last[0] = (last[1] - ctx->D[0] - ctx->D[1] - ctx->D[2] - ctx->D[3])>> 1;
00525 last[2] = (last[1] + last[3]) >> 1;
00526
00527 t1 = ctx->D[0] + ctx->D[1];
00528 ctx->D[0] = t1 >> 1;
00529 ctx->D[1] = t1 - (t1 >> 1);
00530 t2 = ctx->D[2] + ctx->D[3];
00531 ctx->D[2] = t2 >> 1;
00532 ctx->D[3] = t2 - (t2 >> 1);
00533
00534 tm2_apply_deltas(ctx, Y, Ystride, deltas, last);
00535 }
00536
00537 static inline void tm2_null_res_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00538 {
00539 int i;
00540 int ct;
00541 int left, right, diff;
00542 int deltas[16];
00543 TM2_INIT_POINTERS();
00544
00545
00546 deltas[0] = deltas[1] = deltas[2] = deltas[3] = 0;
00547 tm2_low_chroma(U, Ustride, clast, ctx->CD, deltas, bx);
00548
00549 deltas[0] = deltas[1] = deltas[2] = deltas[3] = 0;
00550 tm2_low_chroma(V, Vstride, clast + 2, ctx->CD + 2, deltas, bx);
00551
00552
00553 for(i = 0; i < 16; i++)
00554 deltas[i] = 0;
00555
00556 ct = ctx->D[0] + ctx->D[1] + ctx->D[2] + ctx->D[3];
00557
00558 if(bx > 0)
00559 left = last[-1] - ct;
00560 else
00561 left = 0;
00562
00563 right = last[3];
00564 diff = right - left;
00565 last[0] = left + (diff >> 2);
00566 last[1] = left + (diff >> 1);
00567 last[2] = right - (diff >> 2);
00568 last[3] = right;
00569 {
00570 int tp = left;
00571
00572 ctx->D[0] = (tp + (ct >> 2)) - left;
00573 left += ctx->D[0];
00574 ctx->D[1] = (tp + (ct >> 1)) - left;
00575 left += ctx->D[1];
00576 ctx->D[2] = ((tp + ct) - (ct >> 2)) - left;
00577 left += ctx->D[2];
00578 ctx->D[3] = (tp + ct) - left;
00579 }
00580 tm2_apply_deltas(ctx, Y, Ystride, deltas, last);
00581 }
00582
00583 static inline void tm2_still_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00584 {
00585 int i, j;
00586 TM2_INIT_POINTERS_2();
00587
00588
00589 for(j = 0; j < 2; j++){
00590 for(i = 0; i < 2; i++){
00591 U[i] = Uo[i];
00592 V[i] = Vo[i];
00593 }
00594 U += Ustride; V += Vstride;
00595 Uo += oUstride; Vo += oVstride;
00596 }
00597 U -= Ustride * 2;
00598 V -= Vstride * 2;
00599 TM2_RECALC_BLOCK(U, Ustride, clast, ctx->CD);
00600 TM2_RECALC_BLOCK(V, Vstride, (clast + 2), (ctx->CD + 2));
00601
00602
00603 ctx->D[0] = Yo[3] - last[3];
00604 ctx->D[1] = Yo[3 + oYstride] - Yo[3];
00605 ctx->D[2] = Yo[3 + oYstride * 2] - Yo[3 + oYstride];
00606 ctx->D[3] = Yo[3 + oYstride * 3] - Yo[3 + oYstride * 2];
00607
00608 for(j = 0; j < 4; j++){
00609 for(i = 0; i < 4; i++){
00610 Y[i] = Yo[i];
00611 last[i] = Yo[i];
00612 }
00613 Y += Ystride;
00614 Yo += oYstride;
00615 }
00616 }
00617
00618 static inline void tm2_update_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00619 {
00620 int i, j;
00621 int d;
00622 TM2_INIT_POINTERS_2();
00623
00624
00625 for(j = 0; j < 2; j++){
00626 for(i = 0; i < 2; i++){
00627 U[i] = Uo[i] + GET_TOK(ctx, TM2_UPD);
00628 V[i] = Vo[i] + GET_TOK(ctx, TM2_UPD);
00629 }
00630 U += Ustride; V += Vstride;
00631 Uo += oUstride; Vo += oVstride;
00632 }
00633 U -= Ustride * 2;
00634 V -= Vstride * 2;
00635 TM2_RECALC_BLOCK(U, Ustride, clast, ctx->CD);
00636 TM2_RECALC_BLOCK(V, Vstride, (clast + 2), (ctx->CD + 2));
00637
00638
00639 ctx->D[0] = Yo[3] - last[3];
00640 ctx->D[1] = Yo[3 + oYstride] - Yo[3];
00641 ctx->D[2] = Yo[3 + oYstride * 2] - Yo[3 + oYstride];
00642 ctx->D[3] = Yo[3 + oYstride * 3] - Yo[3 + oYstride * 2];
00643
00644 for(j = 0; j < 4; j++){
00645 d = last[3];
00646 for(i = 0; i < 4; i++){
00647 Y[i] = Yo[i] + GET_TOK(ctx, TM2_UPD);
00648 last[i] = Y[i];
00649 }
00650 ctx->D[j] = last[3] - d;
00651 Y += Ystride;
00652 Yo += oYstride;
00653 }
00654 }
00655
00656 static inline void tm2_motion_block(TM2Context *ctx, AVFrame *pic, int bx, int by)
00657 {
00658 int i, j;
00659 int mx, my;
00660 TM2_INIT_POINTERS_2();
00661
00662 mx = GET_TOK(ctx, TM2_MOT);
00663 my = GET_TOK(ctx, TM2_MOT);
00664 mx = av_clip(mx, -(bx * 4 + 4), ctx->avctx->width - bx * 4);
00665 my = av_clip(my, -(by * 4 + 4), ctx->avctx->height - by * 4);
00666
00667 if (4*bx+mx<0 || 4*by+my<0 || 4*bx+mx+4 > ctx->avctx->width || 4*by+my+4 > ctx->avctx->height) {
00668 av_log(ctx->avctx, AV_LOG_ERROR, "MV out of picture\n");
00669 return;
00670 }
00671
00672 Yo += my * oYstride + mx;
00673 Uo += (my >> 1) * oUstride + (mx >> 1);
00674 Vo += (my >> 1) * oVstride + (mx >> 1);
00675
00676
00677 for(j = 0; j < 2; j++){
00678 for(i = 0; i < 2; i++){
00679 U[i] = Uo[i];
00680 V[i] = Vo[i];
00681 }
00682 U += Ustride; V += Vstride;
00683 Uo += oUstride; Vo += oVstride;
00684 }
00685 U -= Ustride * 2;
00686 V -= Vstride * 2;
00687 TM2_RECALC_BLOCK(U, Ustride, clast, ctx->CD);
00688 TM2_RECALC_BLOCK(V, Vstride, (clast + 2), (ctx->CD + 2));
00689
00690
00691 for(j = 0; j < 4; j++){
00692 for(i = 0; i < 4; i++){
00693 Y[i] = Yo[i];
00694 }
00695 Y += Ystride;
00696 Yo += oYstride;
00697 }
00698
00699 Y -= Ystride * 4;
00700 ctx->D[0] = Y[3] - last[3];
00701 ctx->D[1] = Y[3 + Ystride] - Y[3];
00702 ctx->D[2] = Y[3 + Ystride * 2] - Y[3 + Ystride];
00703 ctx->D[3] = Y[3 + Ystride * 3] - Y[3 + Ystride * 2];
00704 for(i = 0; i < 4; i++)
00705 last[i] = Y[i + Ystride * 3];
00706 }
00707
00708 static int tm2_decode_blocks(TM2Context *ctx, AVFrame *p)
00709 {
00710 int i, j;
00711 int w = ctx->avctx->width, h = ctx->avctx->height, bw = w >> 2, bh = h >> 2, cw = w >> 1;
00712 int type;
00713 int keyframe = 1;
00714 int *Y, *U, *V;
00715 uint8_t *dst;
00716
00717 for(i = 0; i < TM2_NUM_STREAMS; i++)
00718 ctx->tok_ptrs[i] = 0;
00719
00720 if (ctx->tok_lens[TM2_TYPE]<bw*bh){
00721 av_log(ctx->avctx,AV_LOG_ERROR,"Got %i tokens for %i blocks\n",ctx->tok_lens[TM2_TYPE],bw*bh);
00722 return -1;
00723 }
00724
00725 memset(ctx->last, 0, 4 * bw * sizeof(int));
00726 memset(ctx->clast, 0, 4 * bw * sizeof(int));
00727
00728 for(j = 0; j < bh; j++) {
00729 memset(ctx->D, 0, 4 * sizeof(int));
00730 memset(ctx->CD, 0, 4 * sizeof(int));
00731 for(i = 0; i < bw; i++) {
00732 type = GET_TOK(ctx, TM2_TYPE);
00733 switch(type) {
00734 case TM2_HI_RES:
00735 tm2_hi_res_block(ctx, p, i, j);
00736 break;
00737 case TM2_MED_RES:
00738 tm2_med_res_block(ctx, p, i, j);
00739 break;
00740 case TM2_LOW_RES:
00741 tm2_low_res_block(ctx, p, i, j);
00742 break;
00743 case TM2_NULL_RES:
00744 tm2_null_res_block(ctx, p, i, j);
00745 break;
00746 case TM2_UPDATE:
00747 tm2_update_block(ctx, p, i, j);
00748 keyframe = 0;
00749 break;
00750 case TM2_STILL:
00751 tm2_still_block(ctx, p, i, j);
00752 keyframe = 0;
00753 break;
00754 case TM2_MOTION:
00755 tm2_motion_block(ctx, p, i, j);
00756 keyframe = 0;
00757 break;
00758 default:
00759 av_log(ctx->avctx, AV_LOG_ERROR, "Skipping unknown block type %i\n", type);
00760 }
00761 }
00762 }
00763
00764
00765 Y = (ctx->cur?ctx->Y2:ctx->Y1);
00766 U = (ctx->cur?ctx->U2:ctx->U1);
00767 V = (ctx->cur?ctx->V2:ctx->V1);
00768 dst = p->data[0];
00769 for(j = 0; j < h; j++){
00770 for(i = 0; i < w; i++){
00771 int y = Y[i], u = U[i >> 1], v = V[i >> 1];
00772 dst[3*i+0] = av_clip_uint8(y + v);
00773 dst[3*i+1] = av_clip_uint8(y);
00774 dst[3*i+2] = av_clip_uint8(y + u);
00775 }
00776
00777
00778 Y[-4] = Y[-3] = Y[-2] = Y[-1] = Y[0];
00779 Y[w + 3] = Y[w + 2] = Y[w + 1] = Y[w] = Y[w - 1];
00780
00781
00782 if (j == 0) {
00783 memcpy(Y - 4 - 1 * ctx->y_stride, Y - 4, ctx->y_stride);
00784 memcpy(Y - 4 - 2 * ctx->y_stride, Y - 4, ctx->y_stride);
00785 memcpy(Y - 4 - 3 * ctx->y_stride, Y - 4, ctx->y_stride);
00786 memcpy(Y - 4 - 4 * ctx->y_stride, Y - 4, ctx->y_stride);
00787 } else if (j == h - 1) {
00788 memcpy(Y - 4 + 1 * ctx->y_stride, Y - 4, ctx->y_stride);
00789 memcpy(Y - 4 + 2 * ctx->y_stride, Y - 4, ctx->y_stride);
00790 memcpy(Y - 4 + 3 * ctx->y_stride, Y - 4, ctx->y_stride);
00791 memcpy(Y - 4 + 4 * ctx->y_stride, Y - 4, ctx->y_stride);
00792 }
00793
00794 Y += ctx->y_stride;
00795 if (j & 1) {
00796
00797 U[-2] = U[-1] = U[0];
00798 V[-2] = V[-1] = V[0];
00799 U[cw + 1] = U[cw] = U[cw - 1];
00800 V[cw + 1] = V[cw] = V[cw - 1];
00801
00802
00803 if (j == 1) {
00804 memcpy(U - 2 - 1 * ctx->uv_stride, U - 2, ctx->uv_stride);
00805 memcpy(V - 2 - 1 * ctx->uv_stride, V - 2, ctx->uv_stride);
00806 memcpy(U - 2 - 2 * ctx->uv_stride, U - 2, ctx->uv_stride);
00807 memcpy(V - 2 - 2 * ctx->uv_stride, V - 2, ctx->uv_stride);
00808 } else if (j == h - 1) {
00809 memcpy(U - 2 + 1 * ctx->uv_stride, U - 2, ctx->uv_stride);
00810 memcpy(V - 2 + 1 * ctx->uv_stride, V - 2, ctx->uv_stride);
00811 memcpy(U - 2 + 2 * ctx->uv_stride, U - 2, ctx->uv_stride);
00812 memcpy(V - 2 + 2 * ctx->uv_stride, V - 2, ctx->uv_stride);
00813 }
00814
00815 U += ctx->uv_stride;
00816 V += ctx->uv_stride;
00817 }
00818 dst += p->linesize[0];
00819 }
00820
00821 return keyframe;
00822 }
00823
00824 static const int tm2_stream_order[TM2_NUM_STREAMS] = {
00825 TM2_C_HI, TM2_C_LO, TM2_L_HI, TM2_L_LO, TM2_UPD, TM2_MOT, TM2_TYPE
00826 };
00827
00828 static int decode_frame(AVCodecContext *avctx,
00829 void *data, int *got_frame,
00830 AVPacket *avpkt)
00831 {
00832 const uint8_t *buf = avpkt->data;
00833 int buf_size = avpkt->size & ~3;
00834 TM2Context * const l = avctx->priv_data;
00835 AVFrame * const p = &l->pic;
00836 int i, ret, skip, t;
00837
00838 av_fast_padded_malloc(&l->buffer, &l->buffer_size, buf_size);
00839 if(!l->buffer){
00840 av_log(avctx, AV_LOG_ERROR, "Cannot allocate temporary buffer\n");
00841 return AVERROR(ENOMEM);
00842 }
00843 p->reference = 3;
00844 p->buffer_hints = FF_BUFFER_HINTS_VALID | FF_BUFFER_HINTS_PRESERVE | FF_BUFFER_HINTS_REUSABLE;
00845 if((ret = avctx->reget_buffer(avctx, p)) < 0){
00846 av_log(avctx, AV_LOG_ERROR, "get_buffer() failed\n");
00847 return ret;
00848 }
00849
00850 l->dsp.bswap_buf((uint32_t*)l->buffer, (const uint32_t*)buf, buf_size >> 2);
00851 skip = tm2_read_header(l, l->buffer);
00852
00853 if(skip == -1){
00854 return AVERROR_INVALIDDATA;
00855 }
00856
00857 for(i = 0; i < TM2_NUM_STREAMS; i++){
00858 if (skip >= buf_size) {
00859 av_log(avctx, AV_LOG_ERROR, "no space for tm2_read_stream\n");
00860 return AVERROR_INVALIDDATA;
00861 }
00862
00863 t = tm2_read_stream(l, l->buffer + skip, tm2_stream_order[i], buf_size - skip);
00864 if(t < 0){
00865 return t;
00866 }
00867 skip += t;
00868 }
00869 p->key_frame = tm2_decode_blocks(l, p);
00870 if(p->key_frame)
00871 p->pict_type = AV_PICTURE_TYPE_I;
00872 else
00873 p->pict_type = AV_PICTURE_TYPE_P;
00874
00875 l->cur = !l->cur;
00876 *got_frame = 1;
00877 *(AVFrame*)data = l->pic;
00878
00879 return buf_size;
00880 }
00881
00882 static av_cold int decode_init(AVCodecContext *avctx){
00883 TM2Context * const l = avctx->priv_data;
00884 int i, w = avctx->width, h = avctx->height;
00885
00886 if((avctx->width & 3) || (avctx->height & 3)){
00887 av_log(avctx, AV_LOG_ERROR, "Width and height must be multiple of 4\n");
00888 return AVERROR_INVALIDDATA;
00889 }
00890
00891 l->avctx = avctx;
00892 l->pic.data[0]=NULL;
00893 avctx->pix_fmt = AV_PIX_FMT_BGR24;
00894 avcodec_get_frame_defaults(&l->pic);
00895
00896 ff_dsputil_init(&l->dsp, avctx);
00897
00898 l->last = av_malloc(4 * sizeof(*l->last) * (w >> 2));
00899 l->clast = av_malloc(4 * sizeof(*l->clast) * (w >> 2));
00900
00901 for(i = 0; i < TM2_NUM_STREAMS; i++) {
00902 l->tokens[i] = NULL;
00903 l->tok_lens[i] = 0;
00904 }
00905
00906 w += 8;
00907 h += 8;
00908 l->Y1_base = av_malloc(sizeof(*l->Y1_base) * w * h);
00909 l->Y2_base = av_malloc(sizeof(*l->Y2_base) * w * h);
00910 l->y_stride = w;
00911 w = (w + 1) >> 1;
00912 h = (h + 1) >> 1;
00913 l->U1_base = av_malloc(sizeof(*l->U1_base) * w * h);
00914 l->V1_base = av_malloc(sizeof(*l->V1_base) * w * h);
00915 l->U2_base = av_malloc(sizeof(*l->U2_base) * w * h);
00916 l->V2_base = av_malloc(sizeof(*l->V1_base) * w * h);
00917 l->uv_stride = w;
00918 l->cur = 0;
00919 if (!l->Y1_base || !l->Y2_base || !l->U1_base ||
00920 !l->V1_base || !l->U2_base || !l->V2_base ||
00921 !l->last || !l->clast) {
00922 av_freep(l->Y1_base);
00923 av_freep(l->Y2_base);
00924 av_freep(l->U1_base);
00925 av_freep(l->U2_base);
00926 av_freep(l->V1_base);
00927 av_freep(l->V2_base);
00928 av_freep(l->last);
00929 av_freep(l->clast);
00930 return AVERROR(ENOMEM);
00931 }
00932 l->Y1 = l->Y1_base + l->y_stride * 4 + 4;
00933 l->Y2 = l->Y2_base + l->y_stride * 4 + 4;
00934 l->U1 = l->U1_base + l->uv_stride * 2 + 2;
00935 l->U2 = l->U2_base + l->uv_stride * 2 + 2;
00936 l->V1 = l->V1_base + l->uv_stride * 2 + 2;
00937 l->V2 = l->V2_base + l->uv_stride * 2 + 2;
00938
00939 return 0;
00940 }
00941
00942 static av_cold int decode_end(AVCodecContext *avctx){
00943 TM2Context * const l = avctx->priv_data;
00944 AVFrame *pic = &l->pic;
00945 int i;
00946
00947 av_free(l->last);
00948 av_free(l->clast);
00949 for(i = 0; i < TM2_NUM_STREAMS; i++)
00950 av_free(l->tokens[i]);
00951 if(l->Y1){
00952 av_free(l->Y1_base);
00953 av_free(l->U1_base);
00954 av_free(l->V1_base);
00955 av_free(l->Y2_base);
00956 av_free(l->U2_base);
00957 av_free(l->V2_base);
00958 }
00959 av_freep(&l->buffer);
00960 l->buffer_size = 0;
00961
00962 if (pic->data[0])
00963 avctx->release_buffer(avctx, pic);
00964
00965 return 0;
00966 }
00967
00968 AVCodec ff_truemotion2_decoder = {
00969 .name = "truemotion2",
00970 .type = AVMEDIA_TYPE_VIDEO,
00971 .id = AV_CODEC_ID_TRUEMOTION2,
00972 .priv_data_size = sizeof(TM2Context),
00973 .init = decode_init,
00974 .close = decode_end,
00975 .decode = decode_frame,
00976 .capabilities = CODEC_CAP_DR1,
00977 .long_name = NULL_IF_CONFIG_SMALL("Duck TrueMotion 2.0"),
00978 };