1 /*
2 * Apple Pixlet decoder
3 * Copyright (c) 2016 Paul B Mahol
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 #include <stdint.h>
23
27
34
36
37 #define PIXLET_MAGIC 0xDEADBEEF
38
41
47
50
53
57
63
65 {
68 return 0;
69 }
70
72 {
74
78 }
79
81 {
86 return 0;
87 }
88
90 {
93
99
100 for (plane = 0; plane < 3; plane++) {
101 unsigned shift = plane > 0;
104
108
109 for (i = 0; i < NB_LEVELS * 3; i++) {
110 unsigned scale = ctx->
levels - (i / 3);
111
115
116 ctx->
band[
plane][i + 1].
x = (w >> scale) * (((i + 1) % 3) != 2);
117 ctx->
band[
plane][i + 1].
y = (h >> scale) * (((i + 1) % 3) != 1);
118 }
119 }
120
121 return 0;
122 }
123
126 {
129 unsigned cnt1, nbits, k, j = 0, i = 0;
131 int rlen, escape,
flag = 0;
132
133 while (i < size) {
134 nbits =
FFMIN(
ff_clz((state >> 8) + 3) ^ 0x1F, 14);
135
137 if (cnt1 < 8) {
139 if (value <= 1) {
141 escape = ((1 << nbits) - 1) * cnt1;
142 } else {
144 escape = value + ((1 << nbits) - 1) * cnt1 - 1;
145 }
146 } else {
148 }
149
150 value = -((escape +
flag) & 1) | 1;
151 dst[j++] = value * ((escape + flag + 1) >> 1);
152 i++;
153 if (j == width) {
154 j = 0;
156 }
157 state = 120 * (escape +
flag) + state - (120 * state >> 8);
158 flag = 0;
159
160 if (state * 4ULL > 0xFF || i >= size)
161 continue;
162
163 nbits = ((state + 8) >> 5) + (state ?
ff_clz(state) : 32) - 24;
164 escape = av_mod_uintp2(16383, nbits);
166 if (cnt1 > 7) {
168 } else {
170 if (value > 1) {
172 rlen = value + escape * cnt1 - 1;
173 } else {
175 rlen = escape * cnt1;
176 }
177 }
178
179 if (rlen > size - i)
181 i += rlen;
182
183 for (k = 0; k < rlen; k++) {
184 dst[j++] = 0;
185 if (j == width) {
186 j = 0;
188 }
189 }
190
191 state = 0;
192 flag = rlen < 0xFFFF ? 1 : 0;
193 }
194
197 }
198
200 int size,
int c,
int a,
int d,
202 {
205 unsigned cnt1, shbits, rlen, nbits,
length, i = 0, j = 0, k;
206 int ret, escape, pfx,
value, yflag, xflag,
flag = 0;
208
210 if (ret < 0)
211 return ret;
212
213 if (a ^ (a >> 31)) {
214 nbits = 33 -
ff_clz(a ^ (a >> 31));
215 if (nbits > 16)
217 } else {
218 nbits = 1;
219 }
220
221 length = 25 - nbits;
222
223 while (i < size) {
224 if (state >> 8 != -3)
225 value =
ff_clz((state >> 8) + 3) ^ 0x1F;
226 else
227 value = -1;
228
230 if (cnt1 >= length) {
232 } else {
233 pfx = 14 + ((((uint64_t)(value - 14)) >> 32) & (value - 14));
234 if (pfx < 1 || pfx > 25)
236 cnt1 *= (1 << pfx) - 1;
238 if (shbits <= 1) {
240 } else {
242 cnt1 += shbits - 1;
243 }
244 }
245
246 xflag = flag + cnt1;
247 yflag = xflag;
248
249 if (flag + cnt1 == 0) {
250 value = 0;
251 } else {
253 tmp = (int64_t)c * ((yflag + 1) >> 1) + (c >> 1);
254 value = xflag + (
tmp ^ -xflag);
255 }
256
257 i++;
259 if (j == width) {
260 j = 0;
262 }
263 state += (int64_t)d * (uint64_t)yflag - ((int64_t)(d * (uint64_t)
state) >> 8);
264
265 flag = 0;
266
267 if ((uint64_t)state > 0xFF / 4 || i >= size)
268 continue;
269
270 pfx = ((state + 8) >> 5) + (state ?
ff_clz(state) : 32) - 24;
271 escape = av_mod_uintp2(16383, pfx);
273 if (cnt1 < 8) {
274 if (pfx < 1 || pfx > 25)
276
278 if (value > 1) {
280 rlen = value + escape * cnt1 - 1;
281 } else {
283 rlen = escape * cnt1;
284 }
285 } else {
288 else
290
291 rlen = value + 8 * escape;
292 }
293
294 if (rlen > 0xFFFF || i + rlen > size)
296 i += rlen;
297
298 for (k = 0; k < rlen; k++) {
299 dst[j++] = 0;
300 if (j == width) {
301 j = 0;
303 }
304 }
305
306 state = 0;
307 flag = rlen < 0xFFFF ? 1 : 0;
308 }
309
312 }
313
316 {
319 int i, ret;
320
321 for (i = 0; i < ctx->
levels * 3; i++) {
325 int32_t d = bytestream2_get_be32(&ctx->
gb);
326 int16_t *dest = (int16_t *)frame->
data[plane] +
327 ctx->
band[plane][i + 1].
x +
328 ctx->
band[plane][i + 1].
y * stride;
330 uint32_t magic = bytestream2_get_be32(&ctx->
gb);
331
334 "wrong magic number: 0x%08"PRIX32" for plane %d, band %d\n",
335 magic, plane, i);
337 }
338
339 if (a == INT32_MIN)
341
343 c, (b >=
FFABS(a)) ? b : a, d,
345 if (ret < 0) {
347 "error in highpass coefficients for plane %d, band %d\n",
348 plane, i);
349 return ret;
350 }
352 }
353
354 return 0;
355 }
356
359 {
361 int i, j;
362
363 memset(pred, 0, width * sizeof(*pred));
364
365 for (i = 0; i <
height; i++) {
366 val = pred[0] + dst[0];
367 dst[0] = pred[0] =
val;
368 for (j = 1; j <
width; j++) {
369 val = pred[j] + dst[j];
370 dst[j] = pred[j] =
val;
371 dst[j] += dst[j-1];
372 }
374 }
375 }
376
378 {
379 int16_t *low, *high, *ll, *lh, *hl, *hh;
380 int hsize, i, j;
382
383 hsize = size >> 1;
384 low = tmp + 4;
385 high = &low[hsize + 8];
386
387 memcpy(low, dest, size);
388 memcpy(high, dest + hsize, size);
389
390 ll = &low[hsize];
391 lh = &low[hsize];
392 hl = &high[hsize];
393 hh = hl;
394 for (i = 4, j = 2; i; i--, j++, ll--, hh++, lh++, hl--) {
395 low[i - 5] = low[j - 1];
396 lh[0] = ll[-1];
397 high[i - 5] = high[j - 2];
398 hh[0] = hl[-2];
399 }
400
401 for (i = 0; i < hsize; i++) {
402 value = (int64_t) low [i + 1] * -INT64_C(325392907) +
403 (int64_t) low [i + 0] * INT64_C(3687786320) +
404 (int64_t) low [i - 1] * -INT64_C(325392907) +
405 (int64_t) high[i + 0] * INT64_C(1518500249) +
406 (int64_t) high[i - 1] * INT64_C(1518500249);
407 dest[i * 2] = av_clip_int16(((value >> 32) * scale) >> 32);
408 }
409
410 for (i = 0; i < hsize; i++) {
411 value = (int64_t) low [i + 2] * -INT64_C(65078576) +
412 (int64_t) low [i + 1] * INT64_C(1583578880) +
413 (int64_t) low [i + 0] * INT64_C(1583578880) +
414 (int64_t) low [i - 1] * -INT64_C(65078576) +
415 (int64_t) high[i + 1] * INT64_C(303700064) +
416 (int64_t) high[i + 0] * -INT64_C(3644400640) +
417 (int64_t) high[i - 1] * INT64_C(303700064);
418 dest[i * 2 + 1] = av_clip_int16(((value >> 32) * scale) >> 32);
419 }
420 }
421
424 int64_t *scaling_h, int64_t *scaling_v)
425 {
427 unsigned scaled_width, scaled_height;
429 int i, j, k;
430
434
436 int64_t scale_v = scaling_v[i];
437 int64_t scale_h = scaling_h[i];
438 scaled_width <<= 1;
439 scaled_height <<= 1;
440
441 ptr = dest;
442 for (j = 0; j < scaled_height; j++) {
445 }
446
447 for (j = 0; j < scaled_width; j++) {
448 ptr = dest + j;
449 for (k = 0; k < scaled_height; k++) {
450 tmp[k] = *ptr;
452 }
453
455
456 ptr = dest + j;
457 for (k = 0; k < scaled_height; k++) {
458 *ptr = tmp[k];
460 }
461 }
462 }
463 }
464
466 {
467 uint16_t *dsty = (uint16_t *)frame->
data[0];
468 int16_t *srcy = (int16_t *)frame->
data[0];
469 ptrdiff_t stridey = frame->
linesize[0] / 2;
470 int i, j;
471
472 for (j = 0; j <
h; j++) {
473 for (i = 0; i <
w; i++) {
474 if (srcy[i] <= 0)
475 dsty[i] = 0;
476 else if (srcy[i] > ((1 << depth) - 1))
477 dsty[i] = 65535;
478 else
479 dsty[i] = ((int64_t) srcy[i] * srcy[i] * 65535) /
480 ((1 << depth) - 1) / ((1 << depth) - 1);
481 }
482 dsty += stridey;
483 srcy += stridey;
484 }
485 }
486
488 {
489 uint16_t *dstu = (uint16_t *)frame->
data[1];
490 uint16_t *dstv = (uint16_t *)frame->
data[2];
491 int16_t *srcu = (int16_t *)frame->
data[1];
492 int16_t *srcv = (int16_t *)frame->
data[2];
493 ptrdiff_t strideu = frame->
linesize[1] / 2;
494 ptrdiff_t stridev = frame->
linesize[2] / 2;
495 const unsigned add = 1 << (depth - 1);
496 const unsigned shift = 16 - depth;
497 int i, j;
498
499 for (j = 0; j <
h; j++) {
500 for (i = 0; i <
w; i++) {
503 }
504 dstu += strideu;
505 dstv += stridev;
506 srcu += strideu;
507 srcv += stridev;
508 }
509 }
510
513 {
516 unsigned shift = plane > 0;
517 int16_t *dst;
518 int i, ret;
519
520 for (i = ctx->
levels - 1; i >= 0; i--) {
523
524 if (!h || !v)
526
529 }
530
532
533 dst = (int16_t *)frame->
data[plane];
535
538 if (ret < 0)
539 return ret;
540
543 if (ret < 0) {
545 "error in lowpass coefficients for plane %d, top row\n", plane);
546 return ret;
547 }
548
551 if (ret < 0) {
553 "error in lowpass coefficients for plane %d, left column\n",
554 plane);
555 return ret;
556 }
557
561 if (ret < 0) {
563 "error in lowpass coefficients for plane %d, rest\n", plane);
564 return ret;
565 }
566
571 }
572
574 if (ret < 0)
575 return ret;
576
579
581 ctx->
h >> shift, stride, ctx->
scaling[plane][
H],
583
584 return 0;
585 }
586
589 {
594 uint32_t pktsize;
595
597
598 pktsize = bytestream2_get_be32(&ctx->
gb);
602 }
603
604 version = bytestream2_get_le32(&ctx->
gb);
605 if (version != 1)
607
609 if (bytestream2_get_be32(&ctx->
gb) != 1)
612
613 width = bytestream2_get_be32(&ctx->
gb);
614 height = bytestream2_get_be32(&ctx->
gb);
615
616 if ( width > INT_MAX - (1
U << (
NB_LEVELS + 1))
617 || height > INT_MAX - (1
U << (
NB_LEVELS + 1)))
619
622
623 ctx->
levels = bytestream2_get_be32(&ctx->
gb);
626 ctx->
depth = bytestream2_get_be32(&ctx->
gb);
630 }
631
633 if (ret < 0)
634 return ret;
637
638 if (ctx->
w != w || ctx->
h != h) {
642
644 if (ret < 0) {
648 return ret;
649 }
650 }
651
653
657
659 if (ret < 0)
660 return ret;
661
662 for (i = 0; i < 3; i++) {
664 if (ret < 0)
665 return ret;
667 break;
668 }
669
672
673 *got_frame = 1;
674
675 return pktsize;
676 }
677
678 #if HAVE_THREADS
680 {
682
688
689 return 0;
690 }
691 #endif /* HAVE_THREADS */
692
707 };
#define FF_CODEC_CAP_INIT_CLEANUP
The codec allows calling the close function for deallocation even if the init function returned a fai...
SubBand band[4][NB_LEVELS *3+1]
const char const char void * val
#define AVERROR_INVALIDDATA
Invalid data found when processing input.
static int shift(int a, int b)
This structure describes decoded (raw) audio or video data.
ptrdiff_t const GLvoid * data
static unsigned int get_bits(GetBitContext *s, int n)
Read 1-25 bits.
int ff_set_dimensions(AVCodecContext *s, int width, int height)
Check that the provided frame dimensions are valid and set them on the codec context.
enum AVColorRange color_range
MPEG vs JPEG YUV range.
enum AVPixelFormat pix_fmt
Pixel format, see AV_PIX_FMT_xxx.
static av_always_inline void bytestream2_init(GetByteContext *g, const uint8_t *buf, int buf_size)
static int read_low_coeffs(AVCodecContext *avctx, int16_t *dst, int size, int width, ptrdiff_t stride)
void void avpriv_request_sample(void *avc, const char *msg,...) av_printf_format(2
Log a generic warning message about a missing feature.
#define FF_CODEC_CAP_INIT_THREADSAFE
The codec does not modify any global variables in the init function, allowing to call the init functi...
Multithreading support functions.
#define u(width, name, range_min, range_max)
static void free_buffers(AVCodecContext *avctx)
static int get_bits_count(const GetBitContext *s)
bitstream reader API header.
static void filterfn(int16_t *dest, int16_t *tmp, unsigned size, int64_t scale)
#define AV_CODEC_FLAG_GRAY
Only decode/encode grayscale.
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
static int decode_plane(AVCodecContext *avctx, int plane, AVPacket *avpkt, AVFrame *frame)
static av_always_inline void bytestream2_skip(GetByteContext *g, unsigned int size)
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
enum AVColorRange color_range
MPEG vs JPEG YUV range.
static av_always_inline unsigned int bytestream2_get_bytes_left(GetByteContext *g)
int flags
AV_CODEC_FLAG_*.
static void postprocess_luma(AVFrame *frame, int w, int h, int depth)
const char * name
Name of the codec implementation.
#define AV_CODEC_CAP_FRAME_THREADS
Codec supports frame-level multithreading.
#define ONLY_IF_THREADS_ENABLED(x)
Define a function with only the non-default version specified.
static int read_highpass(AVCodecContext *avctx, uint8_t *ptr, int plane, AVFrame *frame)
enum AVPictureType pict_type
Picture type of the frame.
int width
picture width / height.
GLsizei GLboolean const GLfloat * value
static unsigned int show_bits(GetBitContext *s, int n)
Show 1-25 bits.
#define FFABS(a)
Absolute value, Note, INT_MIN / INT64_MIN result in undefined behavior as they are not representable ...
static int init_decoder(AVCodecContext *avctx)
the normal 2^n-1 "JPEG" YUV ranges
static const float pred[4]
#define AV_PIX_FMT_YUV420P16
static av_always_inline int bytestream2_tell(GetByteContext *g)
int64_t scaling[4][2][NB_LEVELS]
static av_cold int pixlet_init(AVCodecContext *avctx)
Libavcodec external API header.
int linesize[AV_NUM_DATA_POINTERS]
For video, size in bytes of each picture line.
static int init_get_bits8(GetBitContext *s, const uint8_t *buffer, int byte_size)
Initialize GetBitContext.
int ff_thread_get_buffer(AVCodecContext *avctx, ThreadFrame *f, int flags)
Wrapper around get_buffer() for frame-multithreaded codecs.
main external API structure.
static unsigned int get_bits1(GetBitContext *s)
Describe the class of an AVClass context structure.
static void skip_bits(GetBitContext *s, int n)
static void postprocess_chroma(AVFrame *frame, int w, int h, int depth)
static av_const int sign_extend(int val, unsigned bits)
static int pixlet_decode_frame(AVCodecContext *avctx, void *data, int *got_frame, AVPacket *avpkt)
static int read_high_coeffs(AVCodecContext *avctx, uint8_t *src, int16_t *dst, int size, int c, int a, int d, int width, ptrdiff_t stride)
uint8_t * data[AV_NUM_DATA_POINTERS]
pointer to the picture/channel planes.
GLint GLenum GLboolean GLsizei stride
common internal api header.
static int get_unary(GetBitContext *gb, int stop, int len)
Get unary code of limited length.
int key_frame
1 -> keyframe, 0-> not
static const uint8_t * align_get_bits(GetBitContext *s)
static void reconstruction(AVCodecContext *avctx, int16_t *dest, unsigned width, unsigned height, ptrdiff_t stride, int64_t *scaling_h, int64_t *scaling_v)
#define av_malloc_array(a, b)
static void lowpass_prediction(int16_t *dst, int16_t *pred, int width, int height, ptrdiff_t stride)
static av_cold int pixlet_close(AVCodecContext *avctx)
This structure stores compressed data.
AVCodec ff_pixlet_decoder
#define AV_CODEC_CAP_DR1
Codec uses get_buffer() for allocating buffers and supports custom allocators.
static av_always_inline av_const unsigned av_clip_uintp2_c(int a, int p)
Clip a signed integer to an unsigned power of two range.