| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | /* | ||
| 2 | * Microsoft Screen 3 (aka Microsoft ATC Screen) decoder | ||
| 3 | * Copyright (c) 2012 Konstantin Shishkov | ||
| 4 | * | ||
| 5 | * This file is part of FFmpeg. | ||
| 6 | * | ||
| 7 | * FFmpeg is free software; you can redistribute it and/or | ||
| 8 | * modify it under the terms of the GNU Lesser General Public | ||
| 9 | * License as published by the Free Software Foundation; either | ||
| 10 | * version 2.1 of the License, or (at your option) any later version. | ||
| 11 | * | ||
| 12 | * FFmpeg 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 GNU | ||
| 15 | * Lesser General Public License for more details. | ||
| 16 | * | ||
| 17 | * You should have received a copy of the GNU Lesser General Public | ||
| 18 | * License along with FFmpeg; if not, write to the Free Software | ||
| 19 | * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA | ||
| 20 | */ | ||
| 21 | |||
| 22 | /** | ||
| 23 | * @file | ||
| 24 | * Microsoft Screen 3 (aka Microsoft ATC Screen) decoder | ||
| 25 | */ | ||
| 26 | |||
| 27 | #include "libavutil/mem.h" | ||
| 28 | #include "avcodec.h" | ||
| 29 | #include "bytestream.h" | ||
| 30 | #include "codec_internal.h" | ||
| 31 | #include "decode.h" | ||
| 32 | #include "mathops.h" | ||
| 33 | #include "mss34dsp.h" | ||
| 34 | |||
| 35 | #define HEADER_SIZE 27 | ||
| 36 | |||
| 37 | #define MODEL2_SCALE 13 | ||
| 38 | #define MODEL_SCALE 15 | ||
| 39 | #define MODEL256_SEC_SCALE 9 | ||
| 40 | |||
| 41 | typedef struct Model2 { | ||
| 42 | int upd_val, till_rescale; | ||
| 43 | unsigned zero_freq, zero_weight; | ||
| 44 | unsigned total_freq, total_weight; | ||
| 45 | } Model2; | ||
| 46 | |||
| 47 | typedef struct Model { | ||
| 48 | int weights[16], freqs[16]; | ||
| 49 | int num_syms; | ||
| 50 | int tot_weight; | ||
| 51 | int upd_val, max_upd_val, till_rescale; | ||
| 52 | } Model; | ||
| 53 | |||
| 54 | typedef struct Model256 { | ||
| 55 | int weights[256], freqs[256]; | ||
| 56 | int tot_weight; | ||
| 57 | int secondary[68]; | ||
| 58 | int sec_size; | ||
| 59 | int upd_val, max_upd_val, till_rescale; | ||
| 60 | } Model256; | ||
| 61 | |||
| 62 | #define RAC_BOTTOM 0x01000000 | ||
| 63 | typedef struct RangeCoder { | ||
| 64 | const uint8_t *src, *src_end; | ||
| 65 | |||
| 66 | uint32_t range, low; | ||
| 67 | int got_error; | ||
| 68 | } RangeCoder; | ||
| 69 | |||
| 70 | enum BlockType { | ||
| 71 | FILL_BLOCK = 0, | ||
| 72 | IMAGE_BLOCK, | ||
| 73 | DCT_BLOCK, | ||
| 74 | HAAR_BLOCK, | ||
| 75 | SKIP_BLOCK | ||
| 76 | }; | ||
| 77 | |||
| 78 | typedef struct BlockTypeContext { | ||
| 79 | int last_type; | ||
| 80 | Model bt_model[5]; | ||
| 81 | } BlockTypeContext; | ||
| 82 | |||
| 83 | typedef struct FillBlockCoder { | ||
| 84 | int fill_val; | ||
| 85 | Model coef_model; | ||
| 86 | } FillBlockCoder; | ||
| 87 | |||
| 88 | typedef struct ImageBlockCoder { | ||
| 89 | Model256 esc_model, vec_entry_model; | ||
| 90 | Model vec_size_model; | ||
| 91 | Model vq_model[125]; | ||
| 92 | } ImageBlockCoder; | ||
| 93 | |||
| 94 | typedef struct DCTBlockCoder { | ||
| 95 | int *prev_dc; | ||
| 96 | ptrdiff_t prev_dc_stride; | ||
| 97 | int prev_dc_height; | ||
| 98 | int quality; | ||
| 99 | uint16_t qmat[64]; | ||
| 100 | Model dc_model; | ||
| 101 | Model2 sign_model; | ||
| 102 | Model256 ac_model; | ||
| 103 | } DCTBlockCoder; | ||
| 104 | |||
| 105 | typedef struct HaarBlockCoder { | ||
| 106 | int quality, scale; | ||
| 107 | Model256 coef_model; | ||
| 108 | Model coef_hi_model; | ||
| 109 | } HaarBlockCoder; | ||
| 110 | |||
| 111 | typedef struct MSS3Context { | ||
| 112 | AVCodecContext *avctx; | ||
| 113 | AVFrame *pic; | ||
| 114 | |||
| 115 | int got_error; | ||
| 116 | RangeCoder coder; | ||
| 117 | BlockTypeContext btype[3]; | ||
| 118 | FillBlockCoder fill_coder[3]; | ||
| 119 | ImageBlockCoder image_coder[3]; | ||
| 120 | DCTBlockCoder dct_coder[3]; | ||
| 121 | HaarBlockCoder haar_coder[3]; | ||
| 122 | |||
| 123 | int dctblock[64]; | ||
| 124 | int hblock[16 * 16]; | ||
| 125 | } MSS3Context; | ||
| 126 | |||
| 127 | |||
| 128 | ✗ | static void model2_reset(Model2 *m) | |
| 129 | { | ||
| 130 | ✗ | m->zero_weight = 1; | |
| 131 | ✗ | m->total_weight = 2; | |
| 132 | ✗ | m->zero_freq = 0x1000; | |
| 133 | ✗ | m->total_freq = 0x2000; | |
| 134 | ✗ | m->upd_val = 4; | |
| 135 | ✗ | m->till_rescale = 4; | |
| 136 | ✗ | } | |
| 137 | |||
| 138 | ✗ | static void model2_update(Model2 *m, int bit) | |
| 139 | { | ||
| 140 | unsigned scale; | ||
| 141 | |||
| 142 | ✗ | if (!bit) | |
| 143 | ✗ | m->zero_weight++; | |
| 144 | ✗ | m->till_rescale--; | |
| 145 | ✗ | if (m->till_rescale) | |
| 146 | ✗ | return; | |
| 147 | |||
| 148 | ✗ | m->total_weight += m->upd_val; | |
| 149 | ✗ | if (m->total_weight > 0x2000) { | |
| 150 | ✗ | m->total_weight = (m->total_weight + 1) >> 1; | |
| 151 | ✗ | m->zero_weight = (m->zero_weight + 1) >> 1; | |
| 152 | ✗ | if (m->total_weight == m->zero_weight) | |
| 153 | ✗ | m->total_weight = m->zero_weight + 1; | |
| 154 | } | ||
| 155 | ✗ | m->upd_val = m->upd_val * 5 >> 2; | |
| 156 | ✗ | if (m->upd_val > 64) | |
| 157 | ✗ | m->upd_val = 64; | |
| 158 | ✗ | scale = 0x80000000u / m->total_weight; | |
| 159 | ✗ | m->zero_freq = m->zero_weight * scale >> 18; | |
| 160 | ✗ | m->total_freq = m->total_weight * scale >> 18; | |
| 161 | ✗ | m->till_rescale = m->upd_val; | |
| 162 | } | ||
| 163 | |||
| 164 | ✗ | static void model_update(Model *m, int val) | |
| 165 | { | ||
| 166 | ✗ | int i, sum = 0; | |
| 167 | unsigned scale; | ||
| 168 | |||
| 169 | ✗ | m->weights[val]++; | |
| 170 | ✗ | m->till_rescale--; | |
| 171 | ✗ | if (m->till_rescale) | |
| 172 | ✗ | return; | |
| 173 | ✗ | m->tot_weight += m->upd_val; | |
| 174 | |||
| 175 | ✗ | if (m->tot_weight > 0x8000) { | |
| 176 | ✗ | m->tot_weight = 0; | |
| 177 | ✗ | for (i = 0; i < m->num_syms; i++) { | |
| 178 | ✗ | m->weights[i] = (m->weights[i] + 1) >> 1; | |
| 179 | ✗ | m->tot_weight += m->weights[i]; | |
| 180 | } | ||
| 181 | } | ||
| 182 | ✗ | scale = 0x80000000u / m->tot_weight; | |
| 183 | ✗ | for (i = 0; i < m->num_syms; i++) { | |
| 184 | ✗ | m->freqs[i] = sum * scale >> 16; | |
| 185 | ✗ | sum += m->weights[i]; | |
| 186 | } | ||
| 187 | |||
| 188 | ✗ | m->upd_val = m->upd_val * 5 >> 2; | |
| 189 | ✗ | if (m->upd_val > m->max_upd_val) | |
| 190 | ✗ | m->upd_val = m->max_upd_val; | |
| 191 | ✗ | m->till_rescale = m->upd_val; | |
| 192 | } | ||
| 193 | |||
| 194 | ✗ | static void model_reset(Model *m) | |
| 195 | { | ||
| 196 | int i; | ||
| 197 | |||
| 198 | ✗ | m->tot_weight = 0; | |
| 199 | ✗ | for (i = 0; i < m->num_syms - 1; i++) | |
| 200 | ✗ | m->weights[i] = 1; | |
| 201 | ✗ | m->weights[m->num_syms - 1] = 0; | |
| 202 | |||
| 203 | ✗ | m->upd_val = m->num_syms; | |
| 204 | ✗ | m->till_rescale = 1; | |
| 205 | ✗ | model_update(m, m->num_syms - 1); | |
| 206 | ✗ | m->till_rescale = | |
| 207 | ✗ | m->upd_val = (m->num_syms + 6) >> 1; | |
| 208 | ✗ | } | |
| 209 | |||
| 210 | ✗ | static av_cold void model_init(Model *m, int num_syms) | |
| 211 | { | ||
| 212 | ✗ | m->num_syms = num_syms; | |
| 213 | ✗ | m->max_upd_val = 8 * num_syms + 48; | |
| 214 | |||
| 215 | ✗ | model_reset(m); | |
| 216 | ✗ | } | |
| 217 | |||
| 218 | ✗ | static void model256_update(Model256 *m, int val) | |
| 219 | { | ||
| 220 | ✗ | int i, sum = 0; | |
| 221 | unsigned scale; | ||
| 222 | ✗ | int send, sidx = 1; | |
| 223 | |||
| 224 | ✗ | m->weights[val]++; | |
| 225 | ✗ | m->till_rescale--; | |
| 226 | ✗ | if (m->till_rescale) | |
| 227 | ✗ | return; | |
| 228 | ✗ | m->tot_weight += m->upd_val; | |
| 229 | |||
| 230 | ✗ | if (m->tot_weight > 0x8000) { | |
| 231 | ✗ | m->tot_weight = 0; | |
| 232 | ✗ | for (i = 0; i < 256; i++) { | |
| 233 | ✗ | m->weights[i] = (m->weights[i] + 1) >> 1; | |
| 234 | ✗ | m->tot_weight += m->weights[i]; | |
| 235 | } | ||
| 236 | } | ||
| 237 | ✗ | scale = 0x80000000u / m->tot_weight; | |
| 238 | ✗ | m->secondary[0] = 0; | |
| 239 | ✗ | for (i = 0; i < 256; i++) { | |
| 240 | ✗ | m->freqs[i] = sum * scale >> 16; | |
| 241 | ✗ | sum += m->weights[i]; | |
| 242 | ✗ | send = m->freqs[i] >> MODEL256_SEC_SCALE; | |
| 243 | ✗ | while (sidx <= send) | |
| 244 | ✗ | m->secondary[sidx++] = i - 1; | |
| 245 | } | ||
| 246 | ✗ | while (sidx < m->sec_size) | |
| 247 | ✗ | m->secondary[sidx++] = 255; | |
| 248 | |||
| 249 | ✗ | m->upd_val = m->upd_val * 5 >> 2; | |
| 250 | ✗ | if (m->upd_val > m->max_upd_val) | |
| 251 | ✗ | m->upd_val = m->max_upd_val; | |
| 252 | ✗ | m->till_rescale = m->upd_val; | |
| 253 | } | ||
| 254 | |||
| 255 | ✗ | static void model256_reset(Model256 *m) | |
| 256 | { | ||
| 257 | int i; | ||
| 258 | |||
| 259 | ✗ | for (i = 0; i < 255; i++) | |
| 260 | ✗ | m->weights[i] = 1; | |
| 261 | ✗ | m->weights[255] = 0; | |
| 262 | |||
| 263 | ✗ | m->tot_weight = 0; | |
| 264 | ✗ | m->upd_val = 256; | |
| 265 | ✗ | m->till_rescale = 1; | |
| 266 | ✗ | model256_update(m, 255); | |
| 267 | ✗ | m->till_rescale = | |
| 268 | ✗ | m->upd_val = (256 + 6) >> 1; | |
| 269 | ✗ | } | |
| 270 | |||
| 271 | ✗ | static av_cold void model256_init(Model256 *m) | |
| 272 | { | ||
| 273 | ✗ | m->max_upd_val = 8 * 256 + 48; | |
| 274 | ✗ | m->sec_size = (1 << 6) + 2; | |
| 275 | |||
| 276 | ✗ | model256_reset(m); | |
| 277 | ✗ | } | |
| 278 | |||
| 279 | ✗ | static void rac_init(RangeCoder *c, const uint8_t *src, int size) | |
| 280 | { | ||
| 281 | int i; | ||
| 282 | |||
| 283 | ✗ | c->src = src; | |
| 284 | ✗ | c->src_end = src + size; | |
| 285 | ✗ | c->low = 0; | |
| 286 | ✗ | for (i = 0; i < FFMIN(size, 4); i++) | |
| 287 | ✗ | c->low = (c->low << 8) | *c->src++; | |
| 288 | ✗ | c->range = 0xFFFFFFFF; | |
| 289 | ✗ | c->got_error = 0; | |
| 290 | ✗ | } | |
| 291 | |||
| 292 | ✗ | static void rac_normalise(RangeCoder *c) | |
| 293 | { | ||
| 294 | for (;;) { | ||
| 295 | ✗ | c->range <<= 8; | |
| 296 | ✗ | c->low <<= 8; | |
| 297 | ✗ | if (c->src < c->src_end) { | |
| 298 | ✗ | c->low |= *c->src++; | |
| 299 | ✗ | } else if (!c->low) { | |
| 300 | ✗ | c->got_error = 1; | |
| 301 | ✗ | c->low = 1; | |
| 302 | } | ||
| 303 | ✗ | if (c->low > c->range) { | |
| 304 | ✗ | c->got_error = 1; | |
| 305 | ✗ | c->low = 1; | |
| 306 | } | ||
| 307 | ✗ | if (c->range >= RAC_BOTTOM) | |
| 308 | ✗ | return; | |
| 309 | } | ||
| 310 | } | ||
| 311 | |||
| 312 | ✗ | static int rac_get_bit(RangeCoder *c) | |
| 313 | { | ||
| 314 | int bit; | ||
| 315 | |||
| 316 | ✗ | c->range >>= 1; | |
| 317 | |||
| 318 | ✗ | bit = (c->range <= c->low); | |
| 319 | ✗ | if (bit) | |
| 320 | ✗ | c->low -= c->range; | |
| 321 | |||
| 322 | ✗ | if (c->range < RAC_BOTTOM) | |
| 323 | ✗ | rac_normalise(c); | |
| 324 | |||
| 325 | ✗ | return bit; | |
| 326 | } | ||
| 327 | |||
| 328 | ✗ | static int rac_get_bits(RangeCoder *c, int nbits) | |
| 329 | { | ||
| 330 | int val; | ||
| 331 | |||
| 332 | ✗ | c->range >>= nbits; | |
| 333 | ✗ | val = c->low / c->range; | |
| 334 | ✗ | c->low -= c->range * val; | |
| 335 | |||
| 336 | ✗ | if (c->range < RAC_BOTTOM) | |
| 337 | ✗ | rac_normalise(c); | |
| 338 | |||
| 339 | ✗ | return val; | |
| 340 | } | ||
| 341 | |||
| 342 | ✗ | static int rac_get_model2_sym(RangeCoder *c, Model2 *m) | |
| 343 | { | ||
| 344 | int bit, helper; | ||
| 345 | |||
| 346 | ✗ | helper = m->zero_freq * (c->range >> MODEL2_SCALE); | |
| 347 | ✗ | bit = (c->low >= helper); | |
| 348 | ✗ | if (bit) { | |
| 349 | ✗ | c->low -= helper; | |
| 350 | ✗ | c->range -= helper; | |
| 351 | } else { | ||
| 352 | ✗ | c->range = helper; | |
| 353 | } | ||
| 354 | |||
| 355 | ✗ | if (c->range < RAC_BOTTOM) | |
| 356 | ✗ | rac_normalise(c); | |
| 357 | |||
| 358 | ✗ | model2_update(m, bit); | |
| 359 | |||
| 360 | ✗ | return bit; | |
| 361 | } | ||
| 362 | |||
| 363 | ✗ | static int rac_get_model_sym(RangeCoder *c, Model *m) | |
| 364 | { | ||
| 365 | int val; | ||
| 366 | int end, end2; | ||
| 367 | unsigned prob, prob2, helper; | ||
| 368 | |||
| 369 | ✗ | prob = 0; | |
| 370 | ✗ | prob2 = c->range; | |
| 371 | ✗ | c->range >>= MODEL_SCALE; | |
| 372 | ✗ | val = 0; | |
| 373 | ✗ | end = m->num_syms >> 1; | |
| 374 | ✗ | end2 = m->num_syms; | |
| 375 | do { | ||
| 376 | ✗ | helper = m->freqs[end] * c->range; | |
| 377 | ✗ | if (helper <= c->low) { | |
| 378 | ✗ | val = end; | |
| 379 | ✗ | prob = helper; | |
| 380 | } else { | ||
| 381 | ✗ | end2 = end; | |
| 382 | ✗ | prob2 = helper; | |
| 383 | } | ||
| 384 | ✗ | end = (end2 + val) >> 1; | |
| 385 | ✗ | } while (end != val); | |
| 386 | ✗ | c->low -= prob; | |
| 387 | ✗ | c->range = prob2 - prob; | |
| 388 | ✗ | if (c->range < RAC_BOTTOM) | |
| 389 | ✗ | rac_normalise(c); | |
| 390 | |||
| 391 | ✗ | model_update(m, val); | |
| 392 | |||
| 393 | ✗ | return val; | |
| 394 | } | ||
| 395 | |||
| 396 | ✗ | static int rac_get_model256_sym(RangeCoder *c, Model256 *m) | |
| 397 | { | ||
| 398 | int val; | ||
| 399 | int start, end; | ||
| 400 | int ssym; | ||
| 401 | unsigned prob, prob2, helper; | ||
| 402 | |||
| 403 | ✗ | prob2 = c->range; | |
| 404 | ✗ | c->range >>= MODEL_SCALE; | |
| 405 | |||
| 406 | ✗ | helper = c->low / c->range; | |
| 407 | ✗ | ssym = helper >> MODEL256_SEC_SCALE; | |
| 408 | ✗ | val = m->secondary[ssym]; | |
| 409 | |||
| 410 | ✗ | end = start = m->secondary[ssym + 1] + 1; | |
| 411 | ✗ | while (end > val + 1) { | |
| 412 | ✗ | ssym = (end + val) >> 1; | |
| 413 | ✗ | if (m->freqs[ssym] <= helper) { | |
| 414 | ✗ | end = start; | |
| 415 | ✗ | val = ssym; | |
| 416 | } else { | ||
| 417 | ✗ | end = (end + val) >> 1; | |
| 418 | ✗ | start = ssym; | |
| 419 | } | ||
| 420 | } | ||
| 421 | ✗ | prob = m->freqs[val] * c->range; | |
| 422 | ✗ | if (val != 255) | |
| 423 | ✗ | prob2 = m->freqs[val + 1] * c->range; | |
| 424 | |||
| 425 | ✗ | c->low -= prob; | |
| 426 | ✗ | c->range = prob2 - prob; | |
| 427 | ✗ | if (c->range < RAC_BOTTOM) | |
| 428 | ✗ | rac_normalise(c); | |
| 429 | |||
| 430 | ✗ | model256_update(m, val); | |
| 431 | |||
| 432 | ✗ | return val; | |
| 433 | } | ||
| 434 | |||
| 435 | ✗ | static int decode_block_type(RangeCoder *c, BlockTypeContext *bt) | |
| 436 | { | ||
| 437 | ✗ | bt->last_type = rac_get_model_sym(c, &bt->bt_model[bt->last_type]); | |
| 438 | |||
| 439 | ✗ | return bt->last_type; | |
| 440 | } | ||
| 441 | |||
| 442 | ✗ | static int decode_coeff(RangeCoder *c, Model *m) | |
| 443 | { | ||
| 444 | int val, sign; | ||
| 445 | |||
| 446 | ✗ | val = rac_get_model_sym(c, m); | |
| 447 | ✗ | if (val) { | |
| 448 | ✗ | sign = rac_get_bit(c); | |
| 449 | ✗ | if (val > 1) { | |
| 450 | ✗ | val--; | |
| 451 | ✗ | val = (1 << val) + rac_get_bits(c, val); | |
| 452 | } | ||
| 453 | ✗ | if (!sign) | |
| 454 | ✗ | val = -val; | |
| 455 | } | ||
| 456 | |||
| 457 | ✗ | return val; | |
| 458 | } | ||
| 459 | |||
| 460 | ✗ | static void decode_fill_block(RangeCoder *c, FillBlockCoder *fc, | |
| 461 | uint8_t *dst, ptrdiff_t stride, int block_size) | ||
| 462 | { | ||
| 463 | int i; | ||
| 464 | |||
| 465 | ✗ | fc->fill_val += decode_coeff(c, &fc->coef_model); | |
| 466 | |||
| 467 | ✗ | for (i = 0; i < block_size; i++, dst += stride) | |
| 468 | ✗ | memset(dst, fc->fill_val, block_size); | |
| 469 | ✗ | } | |
| 470 | |||
| 471 | ✗ | static void decode_image_block(RangeCoder *c, ImageBlockCoder *ic, | |
| 472 | uint8_t *dst, ptrdiff_t stride, int block_size) | ||
| 473 | { | ||
| 474 | int i, j; | ||
| 475 | int vec_size; | ||
| 476 | int vec[4]; | ||
| 477 | int prev_line[16]; | ||
| 478 | int A, B, C; | ||
| 479 | |||
| 480 | ✗ | vec_size = rac_get_model_sym(c, &ic->vec_size_model) + 2; | |
| 481 | ✗ | for (i = 0; i < vec_size; i++) | |
| 482 | ✗ | vec[i] = rac_get_model256_sym(c, &ic->vec_entry_model); | |
| 483 | ✗ | for (; i < 4; i++) | |
| 484 | ✗ | vec[i] = 0; | |
| 485 | ✗ | memset(prev_line, 0, sizeof(prev_line)); | |
| 486 | |||
| 487 | ✗ | for (j = 0; j < block_size; j++) { | |
| 488 | ✗ | A = 0; | |
| 489 | ✗ | B = 0; | |
| 490 | ✗ | for (i = 0; i < block_size; i++) { | |
| 491 | ✗ | C = B; | |
| 492 | ✗ | B = prev_line[i]; | |
| 493 | ✗ | A = rac_get_model_sym(c, &ic->vq_model[A + B * 5 + C * 25]); | |
| 494 | |||
| 495 | ✗ | prev_line[i] = A; | |
| 496 | ✗ | if (A < 4) | |
| 497 | ✗ | dst[i] = vec[A]; | |
| 498 | else | ||
| 499 | ✗ | dst[i] = rac_get_model256_sym(c, &ic->esc_model); | |
| 500 | } | ||
| 501 | ✗ | dst += stride; | |
| 502 | } | ||
| 503 | ✗ | } | |
| 504 | |||
| 505 | ✗ | static int decode_dct(RangeCoder *c, DCTBlockCoder *bc, int *block, | |
| 506 | int bx, int by) | ||
| 507 | { | ||
| 508 | ✗ | int skip, val, sign, pos = 1, zz_pos, dc; | |
| 509 | ✗ | int blk_pos = bx + by * bc->prev_dc_stride; | |
| 510 | |||
| 511 | ✗ | memset(block, 0, sizeof(*block) * 64); | |
| 512 | |||
| 513 | ✗ | dc = decode_coeff(c, &bc->dc_model); | |
| 514 | ✗ | if (by) { | |
| 515 | ✗ | if (bx) { | |
| 516 | int l, tl, t; | ||
| 517 | |||
| 518 | ✗ | l = bc->prev_dc[blk_pos - 1]; | |
| 519 | ✗ | tl = bc->prev_dc[blk_pos - 1 - bc->prev_dc_stride]; | |
| 520 | ✗ | t = bc->prev_dc[blk_pos - bc->prev_dc_stride]; | |
| 521 | |||
| 522 | ✗ | if (FFABS(t - tl) <= FFABS(l - tl)) | |
| 523 | ✗ | dc += l; | |
| 524 | else | ||
| 525 | ✗ | dc += t; | |
| 526 | } else { | ||
| 527 | ✗ | dc += bc->prev_dc[blk_pos - bc->prev_dc_stride]; | |
| 528 | } | ||
| 529 | ✗ | } else if (bx) { | |
| 530 | ✗ | dc += bc->prev_dc[bx - 1]; | |
| 531 | } | ||
| 532 | ✗ | bc->prev_dc[blk_pos] = dc; | |
| 533 | ✗ | block[0] = dc * bc->qmat[0]; | |
| 534 | |||
| 535 | ✗ | while (pos < 64) { | |
| 536 | ✗ | val = rac_get_model256_sym(c, &bc->ac_model); | |
| 537 | ✗ | if (!val) | |
| 538 | ✗ | return 0; | |
| 539 | ✗ | if (val == 0xF0) { | |
| 540 | ✗ | pos += 16; | |
| 541 | ✗ | continue; | |
| 542 | } | ||
| 543 | ✗ | skip = val >> 4; | |
| 544 | ✗ | val = val & 0xF; | |
| 545 | ✗ | if (!val) | |
| 546 | ✗ | return -1; | |
| 547 | ✗ | pos += skip; | |
| 548 | ✗ | if (pos >= 64) | |
| 549 | ✗ | return -1; | |
| 550 | |||
| 551 | ✗ | sign = rac_get_model2_sym(c, &bc->sign_model); | |
| 552 | ✗ | if (val > 1) { | |
| 553 | ✗ | val--; | |
| 554 | ✗ | val = (1 << val) + rac_get_bits(c, val); | |
| 555 | } | ||
| 556 | ✗ | if (!sign) | |
| 557 | ✗ | val = -val; | |
| 558 | |||
| 559 | ✗ | zz_pos = ff_zigzag_direct[pos]; | |
| 560 | ✗ | block[zz_pos] = val * bc->qmat[zz_pos]; | |
| 561 | ✗ | pos++; | |
| 562 | } | ||
| 563 | |||
| 564 | ✗ | return pos == 64 ? 0 : -1; | |
| 565 | } | ||
| 566 | |||
| 567 | ✗ | static void decode_dct_block(RangeCoder *c, DCTBlockCoder *bc, | |
| 568 | uint8_t *dst, ptrdiff_t stride, int block_size, | ||
| 569 | int *block, int mb_x, int mb_y) | ||
| 570 | { | ||
| 571 | int i, j; | ||
| 572 | int bx, by; | ||
| 573 | ✗ | int nblocks = block_size >> 3; | |
| 574 | |||
| 575 | ✗ | bx = mb_x * nblocks; | |
| 576 | ✗ | by = mb_y * nblocks; | |
| 577 | |||
| 578 | ✗ | for (j = 0; j < nblocks; j++) { | |
| 579 | ✗ | for (i = 0; i < nblocks; i++) { | |
| 580 | ✗ | if (decode_dct(c, bc, block, bx + i, by + j)) { | |
| 581 | ✗ | c->got_error = 1; | |
| 582 | ✗ | return; | |
| 583 | } | ||
| 584 | ✗ | ff_mss34_dct_put(dst + i * 8, stride, block); | |
| 585 | } | ||
| 586 | ✗ | dst += 8 * stride; | |
| 587 | } | ||
| 588 | } | ||
| 589 | |||
| 590 | ✗ | static void decode_haar_block(RangeCoder *c, HaarBlockCoder *hc, | |
| 591 | uint8_t *dst, ptrdiff_t stride, | ||
| 592 | int block_size, int *block) | ||
| 593 | { | ||
| 594 | ✗ | const int hsize = block_size >> 1; | |
| 595 | int A, B, C, D, t1, t2, t3, t4; | ||
| 596 | int i, j; | ||
| 597 | |||
| 598 | ✗ | for (j = 0; j < block_size; j++) { | |
| 599 | ✗ | for (i = 0; i < block_size; i++) { | |
| 600 | ✗ | if (i < hsize && j < hsize) | |
| 601 | ✗ | block[i] = rac_get_model256_sym(c, &hc->coef_model); | |
| 602 | else | ||
| 603 | ✗ | block[i] = decode_coeff(c, &hc->coef_hi_model); | |
| 604 | ✗ | block[i] *= hc->scale; | |
| 605 | } | ||
| 606 | ✗ | block += block_size; | |
| 607 | } | ||
| 608 | ✗ | block -= block_size * block_size; | |
| 609 | |||
| 610 | ✗ | for (j = 0; j < hsize; j++) { | |
| 611 | ✗ | for (i = 0; i < hsize; i++) { | |
| 612 | ✗ | A = block[i]; | |
| 613 | ✗ | B = block[i + hsize]; | |
| 614 | ✗ | C = block[i + hsize * block_size]; | |
| 615 | ✗ | D = block[i + hsize * block_size + hsize]; | |
| 616 | |||
| 617 | ✗ | t1 = A - B; | |
| 618 | ✗ | t2 = C - D; | |
| 619 | ✗ | t3 = A + B; | |
| 620 | ✗ | t4 = C + D; | |
| 621 | ✗ | dst[i * 2] = av_clip_uint8(t1 - t2); | |
| 622 | ✗ | dst[i * 2 + stride] = av_clip_uint8(t1 + t2); | |
| 623 | ✗ | dst[i * 2 + 1] = av_clip_uint8(t3 - t4); | |
| 624 | ✗ | dst[i * 2 + 1 + stride] = av_clip_uint8(t3 + t4); | |
| 625 | } | ||
| 626 | ✗ | block += block_size; | |
| 627 | ✗ | dst += stride * 2; | |
| 628 | } | ||
| 629 | ✗ | } | |
| 630 | |||
| 631 | ✗ | static void reset_coders(MSS3Context *ctx, int quality) | |
| 632 | { | ||
| 633 | int i, j; | ||
| 634 | |||
| 635 | ✗ | for (i = 0; i < 3; i++) { | |
| 636 | ✗ | ctx->btype[i].last_type = SKIP_BLOCK; | |
| 637 | ✗ | for (j = 0; j < 5; j++) | |
| 638 | ✗ | model_reset(&ctx->btype[i].bt_model[j]); | |
| 639 | ✗ | ctx->fill_coder[i].fill_val = 0; | |
| 640 | ✗ | model_reset(&ctx->fill_coder[i].coef_model); | |
| 641 | ✗ | model256_reset(&ctx->image_coder[i].esc_model); | |
| 642 | ✗ | model256_reset(&ctx->image_coder[i].vec_entry_model); | |
| 643 | ✗ | model_reset(&ctx->image_coder[i].vec_size_model); | |
| 644 | ✗ | for (j = 0; j < 125; j++) | |
| 645 | ✗ | model_reset(&ctx->image_coder[i].vq_model[j]); | |
| 646 | ✗ | if (ctx->dct_coder[i].quality != quality) { | |
| 647 | ✗ | ctx->dct_coder[i].quality = quality; | |
| 648 | ✗ | ff_mss34_gen_quant_mat(ctx->dct_coder[i].qmat, quality, !i); | |
| 649 | } | ||
| 650 | ✗ | memset(ctx->dct_coder[i].prev_dc, 0, | |
| 651 | sizeof(*ctx->dct_coder[i].prev_dc) * | ||
| 652 | ✗ | ctx->dct_coder[i].prev_dc_stride * | |
| 653 | ✗ | ctx->dct_coder[i].prev_dc_height); | |
| 654 | ✗ | model_reset(&ctx->dct_coder[i].dc_model); | |
| 655 | ✗ | model2_reset(&ctx->dct_coder[i].sign_model); | |
| 656 | ✗ | model256_reset(&ctx->dct_coder[i].ac_model); | |
| 657 | ✗ | if (ctx->haar_coder[i].quality != quality) { | |
| 658 | ✗ | ctx->haar_coder[i].quality = quality; | |
| 659 | ✗ | ctx->haar_coder[i].scale = 17 - 7 * quality / 50; | |
| 660 | } | ||
| 661 | ✗ | model_reset(&ctx->haar_coder[i].coef_hi_model); | |
| 662 | ✗ | model256_reset(&ctx->haar_coder[i].coef_model); | |
| 663 | } | ||
| 664 | ✗ | } | |
| 665 | |||
| 666 | ✗ | static av_cold void init_coders(MSS3Context *ctx) | |
| 667 | { | ||
| 668 | int i, j; | ||
| 669 | |||
| 670 | ✗ | for (i = 0; i < 3; i++) { | |
| 671 | ✗ | for (j = 0; j < 5; j++) | |
| 672 | ✗ | model_init(&ctx->btype[i].bt_model[j], 5); | |
| 673 | ✗ | model_init(&ctx->fill_coder[i].coef_model, 12); | |
| 674 | ✗ | model256_init(&ctx->image_coder[i].esc_model); | |
| 675 | ✗ | model256_init(&ctx->image_coder[i].vec_entry_model); | |
| 676 | ✗ | model_init(&ctx->image_coder[i].vec_size_model, 3); | |
| 677 | ✗ | for (j = 0; j < 125; j++) | |
| 678 | ✗ | model_init(&ctx->image_coder[i].vq_model[j], 5); | |
| 679 | ✗ | model_init(&ctx->dct_coder[i].dc_model, 12); | |
| 680 | ✗ | model256_init(&ctx->dct_coder[i].ac_model); | |
| 681 | ✗ | model_init(&ctx->haar_coder[i].coef_hi_model, 12); | |
| 682 | ✗ | model256_init(&ctx->haar_coder[i].coef_model); | |
| 683 | } | ||
| 684 | ✗ | } | |
| 685 | |||
| 686 | ✗ | static int mss3_decode_frame(AVCodecContext *avctx, AVFrame *rframe, | |
| 687 | int *got_frame, AVPacket *avpkt) | ||
| 688 | { | ||
| 689 | ✗ | const uint8_t *buf = avpkt->data; | |
| 690 | ✗ | int buf_size = avpkt->size; | |
| 691 | ✗ | MSS3Context *c = avctx->priv_data; | |
| 692 | ✗ | RangeCoder *acoder = &c->coder; | |
| 693 | GetByteContext gb; | ||
| 694 | uint8_t *dst[3]; | ||
| 695 | int dec_width, dec_height, dec_x, dec_y, quality, keyframe; | ||
| 696 | int x, y, i, mb_width, mb_height, blk_size, btype; | ||
| 697 | int ret; | ||
| 698 | |||
| 699 | ✗ | if (buf_size < HEADER_SIZE) { | |
| 700 | ✗ | av_log(avctx, AV_LOG_ERROR, | |
| 701 | "Frame should have at least %d bytes, got %d instead\n", | ||
| 702 | HEADER_SIZE, buf_size); | ||
| 703 | ✗ | return AVERROR_INVALIDDATA; | |
| 704 | } | ||
| 705 | |||
| 706 | ✗ | bytestream2_init(&gb, buf, buf_size); | |
| 707 | ✗ | keyframe = bytestream2_get_be32(&gb); | |
| 708 | ✗ | if (keyframe & ~0x301) { | |
| 709 | ✗ | av_log(avctx, AV_LOG_ERROR, "Invalid frame type %X\n", keyframe); | |
| 710 | ✗ | return AVERROR_INVALIDDATA; | |
| 711 | } | ||
| 712 | ✗ | keyframe = !(keyframe & 1); | |
| 713 | ✗ | bytestream2_skip(&gb, 6); | |
| 714 | ✗ | dec_x = bytestream2_get_be16(&gb); | |
| 715 | ✗ | dec_y = bytestream2_get_be16(&gb); | |
| 716 | ✗ | dec_width = bytestream2_get_be16(&gb); | |
| 717 | ✗ | dec_height = bytestream2_get_be16(&gb); | |
| 718 | |||
| 719 | ✗ | if (dec_x + dec_width > avctx->width || | |
| 720 | ✗ | dec_y + dec_height > avctx->height || | |
| 721 | ✗ | (dec_width | dec_height) & 0xF) { | |
| 722 | ✗ | av_log(avctx, AV_LOG_ERROR, "Invalid frame dimensions %dx%d +%d,%d\n", | |
| 723 | dec_width, dec_height, dec_x, dec_y); | ||
| 724 | ✗ | return AVERROR_INVALIDDATA; | |
| 725 | } | ||
| 726 | ✗ | bytestream2_skip(&gb, 4); | |
| 727 | ✗ | quality = bytestream2_get_byte(&gb); | |
| 728 | ✗ | if (quality < 1 || quality > 100) { | |
| 729 | ✗ | av_log(avctx, AV_LOG_ERROR, "Invalid quality setting %d\n", quality); | |
| 730 | ✗ | return AVERROR_INVALIDDATA; | |
| 731 | } | ||
| 732 | ✗ | bytestream2_skip(&gb, 4); | |
| 733 | |||
| 734 | ✗ | if (keyframe && !bytestream2_get_bytes_left(&gb)) { | |
| 735 | ✗ | av_log(avctx, AV_LOG_ERROR, "Keyframe without data found\n"); | |
| 736 | ✗ | return AVERROR_INVALIDDATA; | |
| 737 | } | ||
| 738 | ✗ | if (!keyframe && c->got_error) | |
| 739 | ✗ | return buf_size; | |
| 740 | ✗ | c->got_error = 0; | |
| 741 | |||
| 742 | ✗ | if ((ret = ff_reget_buffer(avctx, c->pic, 0)) < 0) | |
| 743 | ✗ | return ret; | |
| 744 | ✗ | if (keyframe) | |
| 745 | ✗ | c->pic->flags |= AV_FRAME_FLAG_KEY; | |
| 746 | else | ||
| 747 | ✗ | c->pic->flags &= ~AV_FRAME_FLAG_KEY; | |
| 748 | ✗ | c->pic->pict_type = keyframe ? AV_PICTURE_TYPE_I : AV_PICTURE_TYPE_P; | |
| 749 | ✗ | if (!bytestream2_get_bytes_left(&gb)) { | |
| 750 | ✗ | if ((ret = av_frame_ref(rframe, c->pic)) < 0) | |
| 751 | ✗ | return ret; | |
| 752 | ✗ | *got_frame = 1; | |
| 753 | |||
| 754 | ✗ | return buf_size; | |
| 755 | } | ||
| 756 | |||
| 757 | ✗ | reset_coders(c, quality); | |
| 758 | |||
| 759 | ✗ | rac_init(acoder, buf + HEADER_SIZE, buf_size - HEADER_SIZE); | |
| 760 | |||
| 761 | ✗ | mb_width = dec_width >> 4; | |
| 762 | ✗ | mb_height = dec_height >> 4; | |
| 763 | ✗ | dst[0] = c->pic->data[0] + dec_x + dec_y * c->pic->linesize[0]; | |
| 764 | ✗ | dst[1] = c->pic->data[1] + dec_x / 2 + (dec_y / 2) * c->pic->linesize[1]; | |
| 765 | ✗ | dst[2] = c->pic->data[2] + dec_x / 2 + (dec_y / 2) * c->pic->linesize[2]; | |
| 766 | ✗ | for (y = 0; y < mb_height; y++) { | |
| 767 | ✗ | for (x = 0; x < mb_width; x++) { | |
| 768 | ✗ | for (i = 0; i < 3; i++) { | |
| 769 | ✗ | blk_size = 8 << !i; | |
| 770 | |||
| 771 | ✗ | btype = decode_block_type(acoder, c->btype + i); | |
| 772 | ✗ | switch (btype) { | |
| 773 | ✗ | case FILL_BLOCK: | |
| 774 | ✗ | decode_fill_block(acoder, c->fill_coder + i, | |
| 775 | ✗ | dst[i] + x * blk_size, | |
| 776 | ✗ | c->pic->linesize[i], blk_size); | |
| 777 | ✗ | break; | |
| 778 | ✗ | case IMAGE_BLOCK: | |
| 779 | ✗ | decode_image_block(acoder, c->image_coder + i, | |
| 780 | ✗ | dst[i] + x * blk_size, | |
| 781 | ✗ | c->pic->linesize[i], blk_size); | |
| 782 | ✗ | break; | |
| 783 | ✗ | case DCT_BLOCK: | |
| 784 | ✗ | decode_dct_block(acoder, c->dct_coder + i, | |
| 785 | ✗ | dst[i] + x * blk_size, | |
| 786 | ✗ | c->pic->linesize[i], blk_size, | |
| 787 | ✗ | c->dctblock, x, y); | |
| 788 | ✗ | break; | |
| 789 | ✗ | case HAAR_BLOCK: | |
| 790 | ✗ | decode_haar_block(acoder, c->haar_coder + i, | |
| 791 | ✗ | dst[i] + x * blk_size, | |
| 792 | ✗ | c->pic->linesize[i], blk_size, | |
| 793 | ✗ | c->hblock); | |
| 794 | ✗ | break; | |
| 795 | } | ||
| 796 | ✗ | if (c->got_error || acoder->got_error) { | |
| 797 | ✗ | av_log(avctx, AV_LOG_ERROR, "Error decoding block %d,%d\n", | |
| 798 | x, y); | ||
| 799 | ✗ | c->got_error = 1; | |
| 800 | ✗ | return AVERROR_INVALIDDATA; | |
| 801 | } | ||
| 802 | } | ||
| 803 | } | ||
| 804 | ✗ | dst[0] += c->pic->linesize[0] * 16; | |
| 805 | ✗ | dst[1] += c->pic->linesize[1] * 8; | |
| 806 | ✗ | dst[2] += c->pic->linesize[2] * 8; | |
| 807 | } | ||
| 808 | |||
| 809 | ✗ | if ((ret = av_frame_ref(rframe, c->pic)) < 0) | |
| 810 | ✗ | return ret; | |
| 811 | |||
| 812 | ✗ | *got_frame = 1; | |
| 813 | |||
| 814 | ✗ | return buf_size; | |
| 815 | } | ||
| 816 | |||
| 817 | ✗ | static av_cold int mss3_decode_end(AVCodecContext *avctx) | |
| 818 | { | ||
| 819 | ✗ | MSS3Context * const c = avctx->priv_data; | |
| 820 | int i; | ||
| 821 | |||
| 822 | ✗ | av_frame_free(&c->pic); | |
| 823 | ✗ | for (i = 0; i < 3; i++) | |
| 824 | ✗ | av_freep(&c->dct_coder[i].prev_dc); | |
| 825 | |||
| 826 | ✗ | return 0; | |
| 827 | } | ||
| 828 | |||
| 829 | ✗ | static av_cold int mss3_decode_init(AVCodecContext *avctx) | |
| 830 | { | ||
| 831 | ✗ | MSS3Context * const c = avctx->priv_data; | |
| 832 | int i; | ||
| 833 | |||
| 834 | ✗ | c->avctx = avctx; | |
| 835 | |||
| 836 | ✗ | if ((avctx->width & 0xF) || (avctx->height & 0xF)) { | |
| 837 | ✗ | av_log(avctx, AV_LOG_ERROR, | |
| 838 | "Image dimensions should be a multiple of 16.\n"); | ||
| 839 | ✗ | return AVERROR_INVALIDDATA; | |
| 840 | } | ||
| 841 | |||
| 842 | ✗ | c->got_error = 0; | |
| 843 | ✗ | for (i = 0; i < 3; i++) { | |
| 844 | ✗ | int b_width = avctx->width >> (2 + !!i); | |
| 845 | ✗ | int b_height = avctx->height >> (2 + !!i); | |
| 846 | ✗ | c->dct_coder[i].prev_dc_stride = b_width; | |
| 847 | ✗ | c->dct_coder[i].prev_dc_height = b_height; | |
| 848 | ✗ | c->dct_coder[i].prev_dc = av_malloc(sizeof(*c->dct_coder[i].prev_dc) * | |
| 849 | ✗ | b_width * b_height); | |
| 850 | ✗ | if (!c->dct_coder[i].prev_dc) { | |
| 851 | ✗ | av_log(avctx, AV_LOG_ERROR, "Cannot allocate buffer\n"); | |
| 852 | ✗ | return AVERROR(ENOMEM); | |
| 853 | } | ||
| 854 | } | ||
| 855 | |||
| 856 | ✗ | c->pic = av_frame_alloc(); | |
| 857 | ✗ | if (!c->pic) | |
| 858 | ✗ | return AVERROR(ENOMEM); | |
| 859 | |||
| 860 | ✗ | avctx->pix_fmt = AV_PIX_FMT_YUV420P; | |
| 861 | |||
| 862 | ✗ | init_coders(c); | |
| 863 | |||
| 864 | ✗ | return 0; | |
| 865 | } | ||
| 866 | |||
| 867 | const FFCodec ff_msa1_decoder = { | ||
| 868 | .p.name = "msa1", | ||
| 869 | CODEC_LONG_NAME("MS ATC Screen"), | ||
| 870 | .p.type = AVMEDIA_TYPE_VIDEO, | ||
| 871 | .p.id = AV_CODEC_ID_MSA1, | ||
| 872 | .priv_data_size = sizeof(MSS3Context), | ||
| 873 | .init = mss3_decode_init, | ||
| 874 | .close = mss3_decode_end, | ||
| 875 | FF_CODEC_DECODE_CB(mss3_decode_frame), | ||
| 876 | .p.capabilities = AV_CODEC_CAP_DR1, | ||
| 877 | .caps_internal = FF_CODEC_CAP_INIT_CLEANUP, | ||
| 878 | }; | ||
| 879 |