git.sesse.net Git - ffmpeg/blob - libavcodec/pngdec.c

   1 /*
   2  * PNG image format
   3  * Copyright (c) 2003 Fabrice Bellard
   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 //#define DEBUG
  23
  24 #include "libavutil/avassert.h"
  25 #include "libavutil/bprint.h"
  26 #include "libavutil/imgutils.h"
  27 #include "libavutil/stereo3d.h"
  28
  29 #include "avcodec.h"
  30 #include "bytestream.h"
  31 #include "internal.h"
  32 #include "apng.h"
  33 #include "png.h"
  34 #include "pngdsp.h"
  35 #include "thread.h"
  36
  37 #include <zlib.h>
  38
  39 enum PNGHeaderState {
  40     PNG_IHDR = 1 << 0,
  41     PNG_PLTE = 1 << 1,
  42 };
  43
  44 enum PNGImageState {
  45     PNG_IDAT     = 1 << 0,
  46     PNG_ALLIMAGE = 1 << 1,
  47 };
  48
  49 typedef struct PNGDecContext {
  50     PNGDSPContext dsp;
  51     AVCodecContext *avctx;
  52
  53     GetByteContext gb;
  54     ThreadFrame previous_picture;
  55     ThreadFrame last_picture;
  56     ThreadFrame picture;
  57
  58     enum PNGHeaderState hdr_state;
  59     enum PNGImageState pic_state;
  60     int width, height;
  61     int cur_w, cur_h;
  62     int last_w, last_h;
  63     int x_offset, y_offset;
  64     int last_x_offset, last_y_offset;
  65     uint8_t dispose_op, blend_op;
  66     uint8_t last_dispose_op;
  67     int bit_depth;
  68     int color_type;
  69     int compression_type;
  70     int interlace_type;
  71     int filter_type;
  72     int channels;
  73     int bits_per_pixel;
  74     int bpp;
  75     int has_trns;
  76     uint8_t transparent_color_be[6];
  77
  78     uint8_t *image_buf;
  79     int image_linesize;
  80     uint32_t palette[256];
  81     uint8_t *crow_buf;
  82     uint8_t *last_row;
  83     unsigned int last_row_size;
  84     uint8_t *tmp_row;
  85     unsigned int tmp_row_size;
  86     uint8_t *buffer;
  87     int buffer_size;
  88     int pass;
  89     int crow_size; /* compressed row size (include filter type) */
  90     int row_size; /* decompressed row size */
  91     int pass_row_size; /* decompress row size of the current pass */
  92     int y;
  93     z_stream zstream;
  94 } PNGDecContext;
  95
  96 /* Mask to determine which pixels are valid in a pass */
  97 static const uint8_t png_pass_mask[NB_PASSES] = {
  98     0x01, 0x01, 0x11, 0x11, 0x55, 0x55, 0xff,
  99 };
 100
 101 /* Mask to determine which y pixels can be written in a pass */
 102 static const uint8_t png_pass_dsp_ymask[NB_PASSES] = {
 103     0xff, 0xff, 0x0f, 0xff, 0x33, 0xff, 0x55,
 104 };
 105
 106 /* Mask to determine which pixels to overwrite while displaying */
 107 static const uint8_t png_pass_dsp_mask[NB_PASSES] = {
 108     0xff, 0x0f, 0xff, 0x33, 0xff, 0x55, 0xff
 109 };
 110
 111 /* NOTE: we try to construct a good looking image at each pass. width
 112  * is the original image width. We also do pixel format conversion at
 113  * this stage */
 114 static void png_put_interlaced_row(uint8_t *dst, int width,
 115                                    int bits_per_pixel, int pass,
 116                                    int color_type, const uint8_t *src)
 117 {
 118     int x, mask, dsp_mask, j, src_x, b, bpp;
 119     uint8_t *d;
 120     const uint8_t *s;
 121
 122     mask     = png_pass_mask[pass];
 123     dsp_mask = png_pass_dsp_mask[pass];
 124
 125     switch (bits_per_pixel) {
 126     case 1:
 127         src_x = 0;
 128         for (x = 0; x < width; x++) {
 129             j = (x & 7);
 130             if ((dsp_mask << j) & 0x80) {
 131                 b = (src[src_x >> 3] >> (7 - (src_x & 7))) & 1;
 132                 dst[x >> 3] &= 0xFF7F>>j;
 133                 dst[x >> 3] |= b << (7 - j);
 134             }
 135             if ((mask << j) & 0x80)
 136                 src_x++;
 137         }
 138         break;
 139     case 2:
 140         src_x = 0;
 141         for (x = 0; x < width; x++) {
 142             int j2 = 2 * (x & 3);
 143             j = (x & 7);
 144             if ((dsp_mask << j) & 0x80) {
 145                 b = (src[src_x >> 2] >> (6 - 2*(src_x & 3))) & 3;
 146                 dst[x >> 2] &= 0xFF3F>>j2;
 147                 dst[x >> 2] |= b << (6 - j2);
 148             }
 149             if ((mask << j) & 0x80)
 150                 src_x++;
 151         }
 152         break;
 153     case 4:
 154         src_x = 0;
 155         for (x = 0; x < width; x++) {
 156             int j2 = 4*(x&1);
 157             j = (x & 7);
 158             if ((dsp_mask << j) & 0x80) {
 159                 b = (src[src_x >> 1] >> (4 - 4*(src_x & 1))) & 15;
 160                 dst[x >> 1] &= 0xFF0F>>j2;
 161                 dst[x >> 1] |= b << (4 - j2);
 162             }
 163             if ((mask << j) & 0x80)
 164                 src_x++;
 165         }
 166         break;
 167     default:
 168         bpp = bits_per_pixel >> 3;
 169         d   = dst;
 170         s   = src;
 171             for (x = 0; x < width; x++) {
 172                 j = x & 7;
 173                 if ((dsp_mask << j) & 0x80) {
 174                     memcpy(d, s, bpp);
 175                 }
 176                 d += bpp;
 177                 if ((mask << j) & 0x80)
 178                     s += bpp;
 179             }
 180         break;
 181     }
 182 }
 183
 184 void ff_add_png_paeth_prediction(uint8_t *dst, uint8_t *src, uint8_t *top,
 185                                  int w, int bpp)
 186 {
 187     int i;
 188     for (i = 0; i < w; i++) {
 189         int a, b, c, p, pa, pb, pc;
 190
 191         a = dst[i - bpp];
 192         b = top[i];
 193         c = top[i - bpp];
 194
 195         p  = b - c;
 196         pc = a - c;
 197
 198         pa = abs(p);
 199         pb = abs(pc);
 200         pc = abs(p + pc);
 201
 202         if (pa <= pb && pa <= pc)
 203             p = a;
 204         else if (pb <= pc)
 205             p = b;
 206         else
 207             p = c;
 208         dst[i] = p + src[i];
 209     }
 210 }
 211
 212 #define UNROLL1(bpp, op)                                                      \
 213     {                                                                         \
 214         r = dst[0];                                                           \
 215         if (bpp >= 2)                                                         \
 216             g = dst[1];                                                       \
 217         if (bpp >= 3)                                                         \
 218             b = dst[2];                                                       \
 219         if (bpp >= 4)                                                         \
 220             a = dst[3];                                                       \
 221         for (; i <= size - bpp; i += bpp) {                                   \
 222             dst[i + 0] = r = op(r, src[i + 0], last[i + 0]);                  \
 223             if (bpp == 1)                                                     \
 224                 continue;                                                     \
 225             dst[i + 1] = g = op(g, src[i + 1], last[i + 1]);                  \
 226             if (bpp == 2)                                                     \
 227                 continue;                                                     \
 228             dst[i + 2] = b = op(b, src[i + 2], last[i + 2]);                  \
 229             if (bpp == 3)                                                     \
 230                 continue;                                                     \
 231             dst[i + 3] = a = op(a, src[i + 3], last[i + 3]);                  \
 232         }                                                                     \
 233     }
 234
 235 #define UNROLL_FILTER(op)                                                     \
 236     if (bpp == 1) {                                                           \
 237         UNROLL1(1, op)                                                        \
 238     } else if (bpp == 2) {                                                    \
 239         UNROLL1(2, op)                                                        \
 240     } else if (bpp == 3) {                                                    \
 241         UNROLL1(3, op)                                                        \
 242     } else if (bpp == 4) {                                                    \
 243         UNROLL1(4, op)                                                        \
 244     }                                                                         \
 245     for (; i < size; i++) {                                                   \
 246         dst[i] = op(dst[i - bpp], src[i], last[i]);                           \
 247     }
 248
 249 /* NOTE: 'dst' can be equal to 'last' */
 250 static void png_filter_row(PNGDSPContext *dsp, uint8_t *dst, int filter_type,
 251                            uint8_t *src, uint8_t *last, int size, int bpp)
 252 {
 253     int i, p, r, g, b, a;
 254
 255     switch (filter_type) {
 256     case PNG_FILTER_VALUE_NONE:
 257         memcpy(dst, src, size);
 258         break;
 259     case PNG_FILTER_VALUE_SUB:
 260         for (i = 0; i < bpp; i++)
 261             dst[i] = src[i];
 262         if (bpp == 4) {
 263             p = *(int *)dst;
 264             for (; i < size; i += bpp) {
 265                 unsigned s = *(int *)(src + i);
 266                 p = ((s & 0x7f7f7f7f) + (p & 0x7f7f7f7f)) ^ ((s ^ p) & 0x80808080);
 267                 *(int *)(dst + i) = p;
 268             }
 269         } else {
 270 #define OP_SUB(x, s, l) ((x) + (s))
 271             UNROLL_FILTER(OP_SUB);
 272         }
 273         break;
 274     case PNG_FILTER_VALUE_UP:
 275         dsp->add_bytes_l2(dst, src, last, size);
 276         break;
 277     case PNG_FILTER_VALUE_AVG:
 278         for (i = 0; i < bpp; i++) {
 279             p      = (last[i] >> 1);
 280             dst[i] = p + src[i];
 281         }
 282 #define OP_AVG(x, s, l) (((((x) + (l)) >> 1) + (s)) & 0xff)
 283         UNROLL_FILTER(OP_AVG);
 284         break;
 285     case PNG_FILTER_VALUE_PAETH:
 286         for (i = 0; i < bpp; i++) {
 287             p      = last[i];
 288             dst[i] = p + src[i];
 289         }
 290         if (bpp > 2 && size > 4) {
 291             /* would write off the end of the array if we let it process
 292              * the last pixel with bpp=3 */
 293             int w = (bpp & 3) ? size - 3 : size;
 294
 295             if (w > i) {
 296                 dsp->add_paeth_prediction(dst + i, src + i, last + i, size - i, bpp);
 297                 i = w;
 298             }
 299         }
 300         ff_add_png_paeth_prediction(dst + i, src + i, last + i, size - i, bpp);
 301         break;
 302     }
 303 }
 304
 305 /* This used to be called "deloco" in FFmpeg
 306  * and is actually an inverse reversible colorspace transformation */
 307 #define YUV2RGB(NAME, TYPE) \
 308 static void deloco_ ## NAME(TYPE *dst, int size, int alpha) \
 309 { \
 310     int i; \
 311     for (i = 0; i < size; i += 3 + alpha) { \
 312         int g = dst [i + 1]; \
 313         dst[i + 0] += g; \
 314         dst[i + 2] += g; \
 315     } \
 316 }
 317
 318 YUV2RGB(rgb8, uint8_t)
 319 YUV2RGB(rgb16, uint16_t)
 320
 321 /* process exactly one decompressed row */
 322 static void png_handle_row(PNGDecContext *s)
 323 {
 324     uint8_t *ptr, *last_row;
 325     int got_line;
 326
 327     if (!s->interlace_type) {
 328         ptr = s->image_buf + s->image_linesize * (s->y + s->y_offset) + s->x_offset * s->bpp;
 329         if (s->y == 0)
 330             last_row = s->last_row;
 331         else
 332             last_row = ptr - s->image_linesize;
 333
 334         png_filter_row(&s->dsp, ptr, s->crow_buf[0], s->crow_buf + 1,
 335                        last_row, s->row_size, s->bpp);
 336         /* loco lags by 1 row so that it doesn't interfere with top prediction */
 337         if (s->filter_type == PNG_FILTER_TYPE_LOCO && s->y > 0) {
 338             if (s->bit_depth == 16) {
 339                 deloco_rgb16((uint16_t *)(ptr - s->image_linesize), s->row_size / 2,
 340                              s->color_type == PNG_COLOR_TYPE_RGB_ALPHA);
 341             } else {
 342                 deloco_rgb8(ptr - s->image_linesize, s->row_size,
 343                             s->color_type == PNG_COLOR_TYPE_RGB_ALPHA);
 344             }
 345         }
 346         s->y++;
 347         if (s->y == s->cur_h) {
 348             s->pic_state |= PNG_ALLIMAGE;
 349             if (s->filter_type == PNG_FILTER_TYPE_LOCO) {
 350                 if (s->bit_depth == 16) {
 351                     deloco_rgb16((uint16_t *)ptr, s->row_size / 2,
 352                                  s->color_type == PNG_COLOR_TYPE_RGB_ALPHA);
 353                 } else {
 354                     deloco_rgb8(ptr, s->row_size,
 355                                 s->color_type == PNG_COLOR_TYPE_RGB_ALPHA);
 356                 }
 357             }
 358         }
 359     } else {
 360         got_line = 0;
 361         for (;;) {
 362             ptr = s->image_buf + s->image_linesize * (s->y + s->y_offset) + s->x_offset * s->bpp;
 363             if ((ff_png_pass_ymask[s->pass] << (s->y & 7)) & 0x80) {
 364                 /* if we already read one row, it is time to stop to
 365                  * wait for the next one */
 366                 if (got_line)
 367                     break;
 368                 png_filter_row(&s->dsp, s->tmp_row, s->crow_buf[0], s->crow_buf + 1,
 369                                s->last_row, s->pass_row_size, s->bpp);
 370                 FFSWAP(uint8_t *, s->last_row, s->tmp_row);
 371                 FFSWAP(unsigned int, s->last_row_size, s->tmp_row_size);
 372                 got_line = 1;
 373             }
 374             if ((png_pass_dsp_ymask[s->pass] << (s->y & 7)) & 0x80) {
 375                 png_put_interlaced_row(ptr, s->cur_w, s->bits_per_pixel, s->pass,
 376                                        s->color_type, s->last_row);
 377             }
 378             s->y++;
 379             if (s->y == s->cur_h) {
 380                 memset(s->last_row, 0, s->row_size);
 381                 for (;;) {
 382                     if (s->pass == NB_PASSES - 1) {
 383                         s->pic_state |= PNG_ALLIMAGE;
 384                         goto the_end;
 385                     } else {
 386                         s->pass++;
 387                         s->y = 0;
 388                         s->pass_row_size = ff_png_pass_row_size(s->pass,
 389                                                                 s->bits_per_pixel,
 390                                                                 s->cur_w);
 391                         s->crow_size = s->pass_row_size + 1;
 392                         if (s->pass_row_size != 0)
 393                             break;
 394                         /* skip pass if empty row */
 395                     }
 396                 }
 397             }
 398         }
 399 the_end:;
 400     }
 401 }
 402
 403 static int png_decode_idat(PNGDecContext *s, int length)
 404 {
 405     int ret;
 406     s->zstream.avail_in = FFMIN(length, bytestream2_get_bytes_left(&s->gb));
 407     s->zstream.next_in  = (unsigned char *)s->gb.buffer;
 408     bytestream2_skip(&s->gb, length);
 409
 410     /* decode one line if possible */
 411     while (s->zstream.avail_in > 0) {
 412         ret = inflate(&s->zstream, Z_PARTIAL_FLUSH);
 413         if (ret != Z_OK && ret != Z_STREAM_END) {
 414             av_log(s->avctx, AV_LOG_ERROR, "inflate returned error %d\n", ret);
 415             return AVERROR_EXTERNAL;
 416         }
 417         if (s->zstream.avail_out == 0) {
 418             if (!(s->pic_state & PNG_ALLIMAGE)) {
 419                 png_handle_row(s);
 420             }
 421             s->zstream.avail_out = s->crow_size;
 422             s->zstream.next_out  = s->crow_buf;
 423         }
 424         if (ret == Z_STREAM_END && s->zstream.avail_in > 0) {
 425             av_log(NULL, AV_LOG_WARNING,
 426                    "%d undecompressed bytes left in buffer\n", s->zstream.avail_in);
 427             return 0;
 428         }
 429     }
 430     return 0;
 431 }
 432
 433 static int decode_zbuf(AVBPrint *bp, const uint8_t *data,
 434                        const uint8_t *data_end)
 435 {
 436     z_stream zstream;
 437     unsigned char *buf;
 438     unsigned buf_size;
 439     int ret;
 440
 441     zstream.zalloc = ff_png_zalloc;
 442     zstream.zfree  = ff_png_zfree;
 443     zstream.opaque = NULL;
 444     if (inflateInit(&zstream) != Z_OK)
 445         return AVERROR_EXTERNAL;
 446     zstream.next_in  = (unsigned char *)data;
 447     zstream.avail_in = data_end - data;
 448     av_bprint_init(bp, 0, -1);
 449
 450     while (zstream.avail_in > 0) {
 451         av_bprint_get_buffer(bp, 2, &buf, &buf_size);
 452         if (buf_size < 2) {
 453             ret = AVERROR(ENOMEM);
 454             goto fail;
 455         }
 456         zstream.next_out  = buf;
 457         zstream.avail_out = buf_size - 1;
 458         ret = inflate(&zstream, Z_PARTIAL_FLUSH);
 459         if (ret != Z_OK && ret != Z_STREAM_END) {
 460             ret = AVERROR_EXTERNAL;
 461             goto fail;
 462         }
 463         bp->len += zstream.next_out - buf;
 464         if (ret == Z_STREAM_END)
 465             break;
 466     }
 467     inflateEnd(&zstream);
 468     bp->str[bp->len] = 0;
 469     return 0;
 470
 471 fail:
 472     inflateEnd(&zstream);
 473     av_bprint_finalize(bp, NULL);
 474     return ret;
 475 }
 476
 477 static uint8_t *iso88591_to_utf8(const uint8_t *in, size_t size_in)
 478 {
 479     size_t extra = 0, i;
 480     uint8_t *out, *q;
 481
 482     for (i = 0; i < size_in; i++)
 483         extra += in[i] >= 0x80;
 484     if (size_in == SIZE_MAX || extra > SIZE_MAX - size_in - 1)
 485         return NULL;
 486     q = out = av_malloc(size_in + extra + 1);
 487     if (!out)
 488         return NULL;
 489     for (i = 0; i < size_in; i++) {
 490         if (in[i] >= 0x80) {
 491             *(q++) = 0xC0 | (in[i] >> 6);
 492             *(q++) = 0x80 | (in[i] & 0x3F);
 493         } else {
 494             *(q++) = in[i];
 495         }
 496     }
 497     *(q++) = 0;
 498     return out;
 499 }
 500
 501 static int decode_text_chunk(PNGDecContext *s, uint32_t length, int compressed,
 502                              AVDictionary **dict)
 503 {
 504     int ret, method;
 505     const uint8_t *data        = s->gb.buffer;
 506     const uint8_t *data_end    = data + length;
 507     const uint8_t *keyword     = data;
 508     const uint8_t *keyword_end = memchr(keyword, 0, data_end - keyword);
 509     uint8_t *kw_utf8 = NULL, *text, *txt_utf8 = NULL;
 510     unsigned text_len;
 511     AVBPrint bp;
 512
 513     if (!keyword_end)
 514         return AVERROR_INVALIDDATA;
 515     data = keyword_end + 1;
 516
 517     if (compressed) {
 518         if (data == data_end)
 519             return AVERROR_INVALIDDATA;
 520         method = *(data++);
 521         if (method)
 522             return AVERROR_INVALIDDATA;
 523         if ((ret = decode_zbuf(&bp, data, data_end)) < 0)
 524             return ret;
 525         text_len = bp.len;
 526         av_bprint_finalize(&bp, (char **)&text);
 527         if (!text)
 528             return AVERROR(ENOMEM);
 529     } else {
 530         text = (uint8_t *)data;
 531         text_len = data_end - text;
 532     }
 533
 534     kw_utf8  = iso88591_to_utf8(keyword, keyword_end - keyword);
 535     txt_utf8 = iso88591_to_utf8(text, text_len);
 536     if (text != data)
 537         av_free(text);
 538     if (!(kw_utf8 && txt_utf8)) {
 539         av_free(kw_utf8);
 540         av_free(txt_utf8);
 541         return AVERROR(ENOMEM);
 542     }
 543
 544     av_dict_set(dict, kw_utf8, txt_utf8,
 545                 AV_DICT_DONT_STRDUP_KEY | AV_DICT_DONT_STRDUP_VAL);
 546     return 0;
 547 }
 548
 549 static int decode_ihdr_chunk(AVCodecContext *avctx, PNGDecContext *s,
 550                              uint32_t length)
 551 {
 552     if (length != 13)
 553         return AVERROR_INVALIDDATA;
 554
 555     if (s->pic_state & PNG_IDAT) {
 556         av_log(avctx, AV_LOG_ERROR, "IHDR after IDAT\n");
 557         return AVERROR_INVALIDDATA;
 558     }
 559
 560     if (s->hdr_state & PNG_IHDR) {
 561         av_log(avctx, AV_LOG_ERROR, "Multiple IHDR\n");
 562         return AVERROR_INVALIDDATA;
 563     }
 564
 565     s->width  = s->cur_w = bytestream2_get_be32(&s->gb);
 566     s->height = s->cur_h = bytestream2_get_be32(&s->gb);
 567     if (av_image_check_size(s->width, s->height, 0, avctx)) {
 568         s->cur_w = s->cur_h = s->width = s->height = 0;
 569         av_log(avctx, AV_LOG_ERROR, "Invalid image size\n");
 570         return AVERROR_INVALIDDATA;
 571     }
 572     s->bit_depth        = bytestream2_get_byte(&s->gb);
 573     if (s->bit_depth != 1 && s->bit_depth != 2 && s->bit_depth != 4 &&
 574         s->bit_depth != 8 && s->bit_depth != 16) {
 575         av_log(avctx, AV_LOG_ERROR, "Invalid bit depth\n");
 576         goto error;
 577     }
 578     s->color_type       = bytestream2_get_byte(&s->gb);
 579     s->compression_type = bytestream2_get_byte(&s->gb);
 580     s->filter_type      = bytestream2_get_byte(&s->gb);
 581     s->interlace_type   = bytestream2_get_byte(&s->gb);
 582     bytestream2_skip(&s->gb, 4); /* crc */
 583     s->hdr_state |= PNG_IHDR;
 584     if (avctx->debug & FF_DEBUG_PICT_INFO)
 585         av_log(avctx, AV_LOG_DEBUG, "width=%d height=%d depth=%d color_type=%d "
 586                 "compression_type=%d filter_type=%d interlace_type=%d\n",
 587                 s->width, s->height, s->bit_depth, s->color_type,
 588                 s->compression_type, s->filter_type, s->interlace_type);
 589
 590     return 0;
 591 error:
 592     s->cur_w = s->cur_h = s->width = s->height = 0;
 593     s->bit_depth = 8;
 594     return AVERROR_INVALIDDATA;
 595 }
 596
 597 static int decode_phys_chunk(AVCodecContext *avctx, PNGDecContext *s)
 598 {
 599     if (s->pic_state & PNG_IDAT) {
 600         av_log(avctx, AV_LOG_ERROR, "pHYs after IDAT\n");
 601         return AVERROR_INVALIDDATA;
 602     }
 603     avctx->sample_aspect_ratio.num = bytestream2_get_be32(&s->gb);
 604     avctx->sample_aspect_ratio.den = bytestream2_get_be32(&s->gb);
 605     if (avctx->sample_aspect_ratio.num < 0 || avctx->sample_aspect_ratio.den < 0)
 606         avctx->sample_aspect_ratio = (AVRational){ 0, 1 };
 607     bytestream2_skip(&s->gb, 1); /* unit specifier */
 608     bytestream2_skip(&s->gb, 4); /* crc */
 609
 610     return 0;
 611 }
 612
 613 static int decode_idat_chunk(AVCodecContext *avctx, PNGDecContext *s,
 614                              uint32_t length, AVFrame *p)
 615 {
 616     int ret;
 617     size_t byte_depth = s->bit_depth > 8 ? 2 : 1;
 618
 619     if (!(s->hdr_state & PNG_IHDR)) {
 620         av_log(avctx, AV_LOG_ERROR, "IDAT without IHDR\n");
 621         return AVERROR_INVALIDDATA;
 622     }
 623     if (!(s->pic_state & PNG_IDAT)) {
 624         /* init image info */
 625         ret = ff_set_dimensions(avctx, s->width, s->height);
 626         if (ret < 0)
 627             return ret;
 628
 629         s->channels       = ff_png_get_nb_channels(s->color_type);
 630         s->bits_per_pixel = s->bit_depth * s->channels;
 631         s->bpp            = (s->bits_per_pixel + 7) >> 3;
 632         s->row_size       = (s->cur_w * s->bits_per_pixel + 7) >> 3;
 633
 634         if ((s->bit_depth == 2 || s->bit_depth == 4 || s->bit_depth == 8) &&
 635                 s->color_type == PNG_COLOR_TYPE_RGB) {
 636             avctx->pix_fmt = AV_PIX_FMT_RGB24;
 637         } else if ((s->bit_depth == 2 || s->bit_depth == 4 || s->bit_depth == 8) &&
 638                 s->color_type == PNG_COLOR_TYPE_RGB_ALPHA) {
 639             avctx->pix_fmt = AV_PIX_FMT_RGBA;
 640         } else if ((s->bit_depth == 2 || s->bit_depth == 4 || s->bit_depth == 8) &&
 641                 s->color_type == PNG_COLOR_TYPE_GRAY) {
 642             avctx->pix_fmt = AV_PIX_FMT_GRAY8;
 643         } else if (s->bit_depth == 16 &&
 644                 s->color_type == PNG_COLOR_TYPE_GRAY) {
 645             avctx->pix_fmt = AV_PIX_FMT_GRAY16BE;
 646         } else if (s->bit_depth == 16 &&
 647                 s->color_type == PNG_COLOR_TYPE_RGB) {
 648             avctx->pix_fmt = AV_PIX_FMT_RGB48BE;
 649         } else if (s->bit_depth == 16 &&
 650                 s->color_type == PNG_COLOR_TYPE_RGB_ALPHA) {
 651             avctx->pix_fmt = AV_PIX_FMT_RGBA64BE;
 652         } else if ((s->bits_per_pixel == 1 || s->bits_per_pixel == 2 || s->bits_per_pixel == 4 || s->bits_per_pixel == 8) &&
 653                 s->color_type == PNG_COLOR_TYPE_PALETTE) {
 654             avctx->pix_fmt = AV_PIX_FMT_PAL8;
 655         } else if (s->bit_depth == 1 && s->bits_per_pixel == 1 && avctx->codec_id != AV_CODEC_ID_APNG) {
 656             avctx->pix_fmt = AV_PIX_FMT_MONOBLACK;
 657         } else if (s->bit_depth == 8 &&
 658                 s->color_type == PNG_COLOR_TYPE_GRAY_ALPHA) {
 659             avctx->pix_fmt = AV_PIX_FMT_YA8;
 660         } else if (s->bit_depth == 16 &&
 661                 s->color_type == PNG_COLOR_TYPE_GRAY_ALPHA) {
 662             avctx->pix_fmt = AV_PIX_FMT_YA16BE;
 663         } else {
 664             av_log(avctx, AV_LOG_ERROR, "unsupported bit depth %d "
 665                     "and color type %d\n",
 666                     s->bit_depth, s->color_type);
 667             return AVERROR_INVALIDDATA;
 668         }
 669
 670         if (s->has_trns && s->color_type != PNG_COLOR_TYPE_PALETTE) {
 671             switch (avctx->pix_fmt) {
 672             case AV_PIX_FMT_RGB24:
 673                 avctx->pix_fmt = AV_PIX_FMT_RGBA;
 674                 break;
 675
 676             case AV_PIX_FMT_RGB48BE:
 677                 avctx->pix_fmt = AV_PIX_FMT_RGBA64BE;
 678                 break;
 679
 680             case AV_PIX_FMT_GRAY8:
 681                 avctx->pix_fmt = AV_PIX_FMT_YA8;
 682                 break;
 683
 684             case AV_PIX_FMT_GRAY16BE:
 685                 avctx->pix_fmt = AV_PIX_FMT_YA16BE;
 686                 break;
 687
 688             default:
 689                 avpriv_request_sample(avctx, "bit depth %d "
 690                         "and color type %d with TRNS",
 691                         s->bit_depth, s->color_type);
 692                 return AVERROR_INVALIDDATA;
 693             }
 694
 695             s->bpp += byte_depth;
 696         }
 697
 698         if ((ret = ff_thread_get_buffer(avctx, &s->picture, AV_GET_BUFFER_FLAG_REF)) < 0)
 699             return ret;
 700         if (avctx->codec_id == AV_CODEC_ID_APNG && s->last_dispose_op != APNG_DISPOSE_OP_PREVIOUS) {
 701             ff_thread_release_buffer(avctx, &s->previous_picture);
 702             if ((ret = ff_thread_get_buffer(avctx, &s->previous_picture, AV_GET_BUFFER_FLAG_REF)) < 0)
 703                 return ret;
 704         }
 705         p->pict_type        = AV_PICTURE_TYPE_I;
 706         p->key_frame        = 1;
 707         p->interlaced_frame = !!s->interlace_type;
 708
 709         ff_thread_finish_setup(avctx);
 710
 711         /* compute the compressed row size */
 712         if (!s->interlace_type) {
 713             s->crow_size = s->row_size + 1;
 714         } else {
 715             s->pass          = 0;
 716             s->pass_row_size = ff_png_pass_row_size(s->pass,
 717                     s->bits_per_pixel,
 718                     s->cur_w);
 719             s->crow_size = s->pass_row_size + 1;
 720         }
 721         ff_dlog(avctx, "row_size=%d crow_size =%d\n",
 722                 s->row_size, s->crow_size);
 723         s->image_buf      = p->data[0];
 724         s->image_linesize = p->linesize[0];
 725         /* copy the palette if needed */
 726         if (avctx->pix_fmt == AV_PIX_FMT_PAL8)
 727             memcpy(p->data[1], s->palette, 256 * sizeof(uint32_t));
 728         /* empty row is used if differencing to the first row */
 729         av_fast_padded_mallocz(&s->last_row, &s->last_row_size, s->row_size);
 730         if (!s->last_row)
 731             return AVERROR_INVALIDDATA;
 732         if (s->interlace_type ||
 733                 s->color_type == PNG_COLOR_TYPE_RGB_ALPHA) {
 734             av_fast_padded_malloc(&s->tmp_row, &s->tmp_row_size, s->row_size);
 735             if (!s->tmp_row)
 736                 return AVERROR_INVALIDDATA;
 737         }
 738         /* compressed row */
 739         av_fast_padded_malloc(&s->buffer, &s->buffer_size, s->row_size + 16);
 740         if (!s->buffer)
 741             return AVERROR(ENOMEM);
 742
 743         /* we want crow_buf+1 to be 16-byte aligned */
 744         s->crow_buf          = s->buffer + 15;
 745         s->zstream.avail_out = s->crow_size;
 746         s->zstream.next_out  = s->crow_buf;
 747     }
 748
 749     s->pic_state |= PNG_IDAT;
 750
 751     /* set image to non-transparent bpp while decompressing */
 752     if (s->has_trns && s->color_type != PNG_COLOR_TYPE_PALETTE)
 753         s->bpp -= byte_depth;
 754
 755     ret = png_decode_idat(s, length);
 756
 757     if (s->has_trns && s->color_type != PNG_COLOR_TYPE_PALETTE)
 758         s->bpp += byte_depth;
 759
 760     if (ret < 0)
 761         return ret;
 762
 763     bytestream2_skip(&s->gb, 4); /* crc */
 764
 765     return 0;
 766 }
 767
 768 static int decode_plte_chunk(AVCodecContext *avctx, PNGDecContext *s,
 769                              uint32_t length)
 770 {
 771     int n, i, r, g, b;
 772
 773     if ((length % 3) != 0 || length > 256 * 3)
 774         return AVERROR_INVALIDDATA;
 775     /* read the palette */
 776     n = length / 3;
 777     for (i = 0; i < n; i++) {
 778         r = bytestream2_get_byte(&s->gb);
 779         g = bytestream2_get_byte(&s->gb);
 780         b = bytestream2_get_byte(&s->gb);
 781         s->palette[i] = (0xFFU << 24) | (r << 16) | (g << 8) | b;
 782     }
 783     for (; i < 256; i++)
 784         s->palette[i] = (0xFFU << 24);
 785     s->hdr_state |= PNG_PLTE;
 786     bytestream2_skip(&s->gb, 4);     /* crc */
 787
 788     return 0;
 789 }
 790
 791 static int decode_trns_chunk(AVCodecContext *avctx, PNGDecContext *s,
 792                              uint32_t length)
 793 {
 794     int v, i;
 795
 796     if (!(s->hdr_state & PNG_IHDR)) {
 797         av_log(avctx, AV_LOG_ERROR, "trns before IHDR\n");
 798         return AVERROR_INVALIDDATA;
 799     }
 800
 801     if (s->pic_state & PNG_IDAT) {
 802         av_log(avctx, AV_LOG_ERROR, "trns after IDAT\n");
 803         return AVERROR_INVALIDDATA;
 804     }
 805
 806     if (s->color_type == PNG_COLOR_TYPE_PALETTE) {
 807         if (length > 256 || !(s->hdr_state & PNG_PLTE))
 808             return AVERROR_INVALIDDATA;
 809
 810         for (i = 0; i < length; i++) {
 811             unsigned v = bytestream2_get_byte(&s->gb);
 812             s->palette[i] = (s->palette[i] & 0x00ffffff) | (v << 24);
 813         }
 814     } else if (s->color_type == PNG_COLOR_TYPE_GRAY || s->color_type == PNG_COLOR_TYPE_RGB) {
 815         if ((s->color_type == PNG_COLOR_TYPE_GRAY && length != 2) ||
 816             (s->color_type == PNG_COLOR_TYPE_RGB && length != 6) ||
 817             s->bit_depth == 1)
 818             return AVERROR_INVALIDDATA;
 819
 820         for (i = 0; i < length / 2; i++) {
 821             /* only use the least significant bits */
 822             v = av_mod_uintp2(bytestream2_get_be16(&s->gb), s->bit_depth);
 823
 824             if (s->bit_depth > 8)
 825                 AV_WB16(&s->transparent_color_be[2 * i], v);
 826             else
 827                 s->transparent_color_be[i] = v;
 828         }
 829     } else {
 830         return AVERROR_INVALIDDATA;
 831     }
 832
 833     bytestream2_skip(&s->gb, 4); /* crc */
 834     s->has_trns = 1;
 835
 836     return 0;
 837 }
 838
 839 static int decode_iccp_chunk(PNGDecContext *s, int length, AVFrame *f)
 840 {
 841     int ret, cnt = 0;
 842     uint8_t *data, profile_name[82];
 843     AVBPrint bp;
 844     AVFrameSideData *sd;
 845
 846     while ((profile_name[cnt++] = bytestream2_get_byte(&s->gb)) && cnt < 81);
 847     if (cnt > 80) {
 848         av_log(s->avctx, AV_LOG_ERROR, "iCCP with invalid name!\n");
 849         return AVERROR_INVALIDDATA;
 850     }
 851
 852     length = FFMAX(length - cnt, 0);
 853
 854     if (bytestream2_get_byte(&s->gb) != 0) {
 855         av_log(s->avctx, AV_LOG_ERROR, "iCCP with invalid compression!\n");
 856         return AVERROR_INVALIDDATA;
 857     }
 858
 859     length = FFMAX(length - 1, 0);
 860
 861     if ((ret = decode_zbuf(&bp, s->gb.buffer, s->gb.buffer + length)) < 0)
 862         return ret;
 863
 864     av_bprint_finalize(&bp, (char **)&data);
 865
 866     sd = av_frame_new_side_data(f, AV_FRAME_DATA_ICC_PROFILE, bp.len);
 867     if (!sd) {
 868         av_free(data);
 869         return AVERROR(ENOMEM);
 870     }
 871
 872     av_dict_set(&sd->metadata, "name", profile_name, 0);
 873     memcpy(sd->data, data, bp.len);
 874     av_free(data);
 875
 876     /* ICC compressed data and CRC */
 877     bytestream2_skip(&s->gb, length + 4);
 878
 879     return 0;
 880 }
 881
 882 static void handle_small_bpp(PNGDecContext *s, AVFrame *p)
 883 {
 884     if (s->bits_per_pixel == 1 && s->color_type == PNG_COLOR_TYPE_PALETTE) {
 885         int i, j, k;
 886         uint8_t *pd = p->data[0];
 887         for (j = 0; j < s->height; j++) {
 888             i = s->width / 8;
 889             for (k = 7; k >= 1; k--)
 890                 if ((s->width&7) >= k)
 891                     pd[8*i + k - 1] = (pd[i]>>8-k) & 1;
 892             for (i--; i >= 0; i--) {
 893                 pd[8*i + 7]=  pd[i]     & 1;
 894                 pd[8*i + 6]= (pd[i]>>1) & 1;
 895                 pd[8*i + 5]= (pd[i]>>2) & 1;
 896                 pd[8*i + 4]= (pd[i]>>3) & 1;
 897                 pd[8*i + 3]= (pd[i]>>4) & 1;
 898                 pd[8*i + 2]= (pd[i]>>5) & 1;
 899                 pd[8*i + 1]= (pd[i]>>6) & 1;
 900                 pd[8*i + 0]=  pd[i]>>7;
 901             }
 902             pd += s->image_linesize;
 903         }
 904     } else if (s->bits_per_pixel == 2) {
 905         int i, j;
 906         uint8_t *pd = p->data[0];
 907         for (j = 0; j < s->height; j++) {
 908             i = s->width / 4;
 909             if (s->color_type == PNG_COLOR_TYPE_PALETTE) {
 910                 if ((s->width&3) >= 3) pd[4*i + 2]= (pd[i] >> 2) & 3;
 911                 if ((s->width&3) >= 2) pd[4*i + 1]= (pd[i] >> 4) & 3;
 912                 if ((s->width&3) >= 1) pd[4*i + 0]=  pd[i] >> 6;
 913                 for (i--; i >= 0; i--) {
 914                     pd[4*i + 3]=  pd[i]     & 3;
 915                     pd[4*i + 2]= (pd[i]>>2) & 3;
 916                     pd[4*i + 1]= (pd[i]>>4) & 3;
 917                     pd[4*i + 0]=  pd[i]>>6;
 918                 }
 919             } else {
 920                 if ((s->width&3) >= 3) pd[4*i + 2]= ((pd[i]>>2) & 3)*0x55;
 921                 if ((s->width&3) >= 2) pd[4*i + 1]= ((pd[i]>>4) & 3)*0x55;
 922                 if ((s->width&3) >= 1) pd[4*i + 0]= ( pd[i]>>6     )*0x55;
 923                 for (i--; i >= 0; i--) {
 924                     pd[4*i + 3]= ( pd[i]     & 3)*0x55;
 925                     pd[4*i + 2]= ((pd[i]>>2) & 3)*0x55;
 926                     pd[4*i + 1]= ((pd[i]>>4) & 3)*0x55;
 927                     pd[4*i + 0]= ( pd[i]>>6     )*0x55;
 928                 }
 929             }
 930             pd += s->image_linesize;
 931         }
 932     } else if (s->bits_per_pixel == 4) {
 933         int i, j;
 934         uint8_t *pd = p->data[0];
 935         for (j = 0; j < s->height; j++) {
 936             i = s->width/2;
 937             if (s->color_type == PNG_COLOR_TYPE_PALETTE) {
 938                 if (s->width&1) pd[2*i+0]= pd[i]>>4;
 939                 for (i--; i >= 0; i--) {
 940                     pd[2*i + 1] = pd[i] & 15;
 941                     pd[2*i + 0] = pd[i] >> 4;
 942                 }
 943             } else {
 944                 if (s->width & 1) pd[2*i + 0]= (pd[i] >> 4) * 0x11;
 945                 for (i--; i >= 0; i--) {
 946                     pd[2*i + 1] = (pd[i] & 15) * 0x11;
 947                     pd[2*i + 0] = (pd[i] >> 4) * 0x11;
 948                 }
 949             }
 950             pd += s->image_linesize;
 951         }
 952     }
 953 }
 954
 955 static int decode_fctl_chunk(AVCodecContext *avctx, PNGDecContext *s,
 956                              uint32_t length)
 957 {
 958     uint32_t sequence_number;
 959     int cur_w, cur_h, x_offset, y_offset, dispose_op, blend_op;
 960
 961     if (length != 26)
 962         return AVERROR_INVALIDDATA;
 963
 964     if (!(s->hdr_state & PNG_IHDR)) {
 965         av_log(avctx, AV_LOG_ERROR, "fctl before IHDR\n");
 966         return AVERROR_INVALIDDATA;
 967     }
 968
 969     s->last_w = s->cur_w;
 970     s->last_h = s->cur_h;
 971     s->last_x_offset = s->x_offset;
 972     s->last_y_offset = s->y_offset;
 973     s->last_dispose_op = s->dispose_op;
 974
 975     sequence_number = bytestream2_get_be32(&s->gb);
 976     cur_w           = bytestream2_get_be32(&s->gb);
 977     cur_h           = bytestream2_get_be32(&s->gb);
 978     x_offset        = bytestream2_get_be32(&s->gb);
 979     y_offset        = bytestream2_get_be32(&s->gb);
 980     bytestream2_skip(&s->gb, 4); /* delay_num (2), delay_den (2) */
 981     dispose_op      = bytestream2_get_byte(&s->gb);
 982     blend_op        = bytestream2_get_byte(&s->gb);
 983     bytestream2_skip(&s->gb, 4); /* crc */
 984
 985     if (sequence_number == 0 &&
 986         (cur_w != s->width ||
 987          cur_h != s->height ||
 988          x_offset != 0 ||
 989          y_offset != 0) ||
 990         cur_w <= 0 || cur_h <= 0 ||
 991         x_offset < 0 || y_offset < 0 ||
 992         cur_w > s->width - x_offset|| cur_h > s->height - y_offset)
 993             return AVERROR_INVALIDDATA;
 994
 995     if (blend_op != APNG_BLEND_OP_OVER && blend_op != APNG_BLEND_OP_SOURCE) {
 996         av_log(avctx, AV_LOG_ERROR, "Invalid blend_op %d\n", blend_op);
 997         return AVERROR_INVALIDDATA;
 998     }
 999
1000     if ((sequence_number == 0 || !s->previous_picture.f->data[0]) &&
1001         dispose_op == APNG_DISPOSE_OP_PREVIOUS) {
1002         // No previous frame to revert to for the first frame
1003         // Spec says to just treat it as a APNG_DISPOSE_OP_BACKGROUND
1004         dispose_op = APNG_DISPOSE_OP_BACKGROUND;
1005     }
1006
1007     if (blend_op == APNG_BLEND_OP_OVER && !s->has_trns && (
1008             avctx->pix_fmt == AV_PIX_FMT_RGB24 ||
1009             avctx->pix_fmt == AV_PIX_FMT_RGB48BE ||
1010             avctx->pix_fmt == AV_PIX_FMT_PAL8 ||
1011             avctx->pix_fmt == AV_PIX_FMT_GRAY8 ||
1012             avctx->pix_fmt == AV_PIX_FMT_GRAY16BE ||
1013             avctx->pix_fmt == AV_PIX_FMT_MONOBLACK
1014         )) {
1015         // APNG_BLEND_OP_OVER is the same as APNG_BLEND_OP_SOURCE when there is no alpha channel
1016         blend_op = APNG_BLEND_OP_SOURCE;
1017     }
1018
1019     s->cur_w      = cur_w;
1020     s->cur_h      = cur_h;
1021     s->x_offset   = x_offset;
1022     s->y_offset   = y_offset;
1023     s->dispose_op = dispose_op;
1024     s->blend_op   = blend_op;
1025
1026     return 0;
1027 }
1028
1029 static void handle_p_frame_png(PNGDecContext *s, AVFrame *p)
1030 {
1031     int i, j;
1032     uint8_t *pd      = p->data[0];
1033     uint8_t *pd_last = s->last_picture.f->data[0];
1034     int ls = FFMIN(av_image_get_linesize(p->format, s->width, 0), s->width * s->bpp);
1035
1036     ff_thread_await_progress(&s->last_picture, INT_MAX, 0);
1037     for (j = 0; j < s->height; j++) {
1038         for (i = 0; i < ls; i++)
1039             pd[i] += pd_last[i];
1040         pd      += s->image_linesize;
1041         pd_last += s->image_linesize;
1042     }
1043 }
1044
1045 // divide by 255 and round to nearest
1046 // apply a fast variant: (X+127)/255 = ((X+127)*257+257)>>16 = ((X+128)*257)>>16
1047 #define FAST_DIV255(x) ((((x) + 128) * 257) >> 16)
1048
1049 static int handle_p_frame_apng(AVCodecContext *avctx, PNGDecContext *s,
1050                                AVFrame *p)
1051 {
1052     size_t x, y;
1053     uint8_t *buffer;
1054
1055     if (s->blend_op == APNG_BLEND_OP_OVER &&
1056         avctx->pix_fmt != AV_PIX_FMT_RGBA &&
1057         avctx->pix_fmt != AV_PIX_FMT_GRAY8A &&
1058         avctx->pix_fmt != AV_PIX_FMT_PAL8) {
1059         avpriv_request_sample(avctx, "Blending with pixel format %s",
1060                               av_get_pix_fmt_name(avctx->pix_fmt));
1061         return AVERROR_PATCHWELCOME;
1062     }
1063
1064     buffer = av_malloc_array(s->image_linesize, s->height);
1065     if (!buffer)
1066         return AVERROR(ENOMEM);
1067
1068
1069     // Do the disposal operation specified by the last frame on the frame
1070     if (s->last_dispose_op != APNG_DISPOSE_OP_PREVIOUS) {
1071         ff_thread_await_progress(&s->last_picture, INT_MAX, 0);
1072         memcpy(buffer, s->last_picture.f->data[0], s->image_linesize * s->height);
1073
1074         if (s->last_dispose_op == APNG_DISPOSE_OP_BACKGROUND)
1075             for (y = s->last_y_offset; y < s->last_y_offset + s->last_h; ++y)
1076                 memset(buffer + s->image_linesize * y + s->bpp * s->last_x_offset, 0, s->bpp * s->last_w);
1077
1078         memcpy(s->previous_picture.f->data[0], buffer, s->image_linesize * s->height);
1079         ff_thread_report_progress(&s->previous_picture, INT_MAX, 0);
1080     } else {
1081         ff_thread_await_progress(&s->previous_picture, INT_MAX, 0);
1082         memcpy(buffer, s->previous_picture.f->data[0], s->image_linesize * s->height);
1083     }
1084
1085     // Perform blending
1086     if (s->blend_op == APNG_BLEND_OP_SOURCE) {
1087         for (y = s->y_offset; y < s->y_offset + s->cur_h; ++y) {
1088             size_t row_start = s->image_linesize * y + s->bpp * s->x_offset;
1089             memcpy(buffer + row_start, p->data[0] + row_start, s->bpp * s->cur_w);
1090         }
1091     } else { // APNG_BLEND_OP_OVER
1092         for (y = s->y_offset; y < s->y_offset + s->cur_h; ++y) {
1093             uint8_t *foreground = p->data[0] + s->image_linesize * y + s->bpp * s->x_offset;
1094             uint8_t *background = buffer + s->image_linesize * y + s->bpp * s->x_offset;
1095             for (x = s->x_offset; x < s->x_offset + s->cur_w; ++x, foreground += s->bpp, background += s->bpp) {
1096                 size_t b;
1097                 uint8_t foreground_alpha, background_alpha, output_alpha;
1098                 uint8_t output[10];
1099
1100                 // Since we might be blending alpha onto alpha, we use the following equations:
1101                 // output_alpha = foreground_alpha + (1 - foreground_alpha) * background_alpha
1102                 // output = (foreground_alpha * foreground + (1 - foreground_alpha) * background_alpha * background) / output_alpha
1103
1104                 switch (avctx->pix_fmt) {
1105                 case AV_PIX_FMT_RGBA:
1106                     foreground_alpha = foreground[3];
1107                     background_alpha = background[3];
1108                     break;
1109
1110                 case AV_PIX_FMT_GRAY8A:
1111                     foreground_alpha = foreground[1];
1112                     background_alpha = background[1];
1113                     break;
1114
1115                 case AV_PIX_FMT_PAL8:
1116                     foreground_alpha = s->palette[foreground[0]] >> 24;
1117                     background_alpha = s->palette[background[0]] >> 24;
1118                     break;
1119                 }
1120
1121                 if (foreground_alpha == 0)
1122                     continue;
1123
1124                 if (foreground_alpha == 255) {
1125                     memcpy(background, foreground, s->bpp);
1126                     continue;
1127                 }
1128
1129                 if (avctx->pix_fmt == AV_PIX_FMT_PAL8) {
1130                     // TODO: Alpha blending with PAL8 will likely need the entire image converted over to RGBA first
1131                     avpriv_request_sample(avctx, "Alpha blending palette samples");
1132                     background[0] = foreground[0];
1133                     continue;
1134                 }
1135
1136                 output_alpha = foreground_alpha + FAST_DIV255((255 - foreground_alpha) * background_alpha);
1137
1138                 av_assert0(s->bpp <= 10);
1139
1140                 for (b = 0; b < s->bpp - 1; ++b) {
1141                     if (output_alpha == 0) {
1142                         output[b] = 0;
1143                     } else if (background_alpha == 255) {
1144                         output[b] = FAST_DIV255(foreground_alpha * foreground[b] + (255 - foreground_alpha) * background[b]);
1145                     } else {
1146                         output[b] = (255 * foreground_alpha * foreground[b] + (255 - foreground_alpha) * background_alpha * background[b]) / (255 * output_alpha);
1147                     }
1148                 }
1149                 output[b] = output_alpha;
1150                 memcpy(background, output, s->bpp);
1151             }
1152         }
1153     }
1154
1155     // Copy blended buffer into the frame and free
1156     memcpy(p->data[0], buffer, s->image_linesize * s->height);
1157     av_free(buffer);
1158
1159     return 0;
1160 }
1161
1162 static int decode_frame_common(AVCodecContext *avctx, PNGDecContext *s,
1163                                AVFrame *p, AVPacket *avpkt)
1164 {
1165     AVDictionary **metadatap = NULL;
1166     uint32_t tag, length;
1167     int decode_next_dat = 0;
1168     int ret;
1169
1170     for (;;) {
1171         length = bytestream2_get_bytes_left(&s->gb);
1172         if (length <= 0) {
1173
1174             if (avctx->codec_id == AV_CODEC_ID_PNG &&
1175                 avctx->skip_frame == AVDISCARD_ALL) {
1176                 return 0;
1177             }
1178
1179             if (CONFIG_APNG_DECODER && avctx->codec_id == AV_CODEC_ID_APNG && length == 0) {
1180                 if (!(s->pic_state & PNG_IDAT))
1181                     return 0;
1182                 else
1183                     goto exit_loop;
1184             }
1185             av_log(avctx, AV_LOG_ERROR, "%d bytes left\n", length);
1186             if (   s->pic_state & PNG_ALLIMAGE
1187                 && avctx->strict_std_compliance <= FF_COMPLIANCE_NORMAL)
1188                 goto exit_loop;
1189             ret = AVERROR_INVALIDDATA;
1190             goto fail;
1191         }
1192
1193         length = bytestream2_get_be32(&s->gb);
1194         if (length > 0x7fffffff || length > bytestream2_get_bytes_left(&s->gb)) {
1195             av_log(avctx, AV_LOG_ERROR, "chunk too big\n");
1196             ret = AVERROR_INVALIDDATA;
1197             goto fail;
1198         }
1199         tag = bytestream2_get_le32(&s->gb);
1200         if (avctx->debug & FF_DEBUG_STARTCODE)
1201             av_log(avctx, AV_LOG_DEBUG, "png: tag=%s length=%u\n",
1202                    av_fourcc2str(tag), length);
1203
1204         if (avctx->codec_id == AV_CODEC_ID_PNG &&
1205             avctx->skip_frame == AVDISCARD_ALL) {
1206             switch(tag) {
1207             case MKTAG('I', 'H', 'D', 'R'):
1208             case MKTAG('p', 'H', 'Y', 's'):
1209             case MKTAG('t', 'E', 'X', 't'):
1210             case MKTAG('I', 'D', 'A', 'T'):
1211             case MKTAG('t', 'R', 'N', 'S'):
1212                 break;
1213             default:
1214                 goto skip_tag;
1215             }
1216         }
1217
1218         metadatap = &p->metadata;
1219         switch (tag) {
1220         case MKTAG('I', 'H', 'D', 'R'):
1221             if ((ret = decode_ihdr_chunk(avctx, s, length)) < 0)
1222                 goto fail;
1223             break;
1224         case MKTAG('p', 'H', 'Y', 's'):
1225             if ((ret = decode_phys_chunk(avctx, s)) < 0)
1226                 goto fail;
1227             break;
1228         case MKTAG('f', 'c', 'T', 'L'):
1229             if (!CONFIG_APNG_DECODER || avctx->codec_id != AV_CODEC_ID_APNG)
1230                 goto skip_tag;
1231             if ((ret = decode_fctl_chunk(avctx, s, length)) < 0)
1232                 goto fail;
1233             decode_next_dat = 1;
1234             break;
1235         case MKTAG('f', 'd', 'A', 'T'):
1236             if (!CONFIG_APNG_DECODER || avctx->codec_id != AV_CODEC_ID_APNG)
1237                 goto skip_tag;
1238             if (!decode_next_dat) {
1239                 ret = AVERROR_INVALIDDATA;
1240                 goto fail;
1241             }
1242             bytestream2_get_be32(&s->gb);
1243             length -= 4;
1244             /* fallthrough */
1245         case MKTAG('I', 'D', 'A', 'T'):
1246             if (CONFIG_APNG_DECODER && avctx->codec_id == AV_CODEC_ID_APNG && !decode_next_dat)
1247                 goto skip_tag;
1248             if ((ret = decode_idat_chunk(avctx, s, length, p)) < 0)
1249                 goto fail;
1250             break;
1251         case MKTAG('P', 'L', 'T', 'E'):
1252             if (decode_plte_chunk(avctx, s, length) < 0)
1253                 goto skip_tag;
1254             break;
1255         case MKTAG('t', 'R', 'N', 'S'):
1256             if (decode_trns_chunk(avctx, s, length) < 0)
1257                 goto skip_tag;
1258             break;
1259         case MKTAG('t', 'E', 'X', 't'):
1260             if (decode_text_chunk(s, length, 0, metadatap) < 0)
1261                 av_log(avctx, AV_LOG_WARNING, "Broken tEXt chunk\n");
1262             bytestream2_skip(&s->gb, length + 4);
1263             break;
1264         case MKTAG('z', 'T', 'X', 't'):
1265             if (decode_text_chunk(s, length, 1, metadatap) < 0)
1266                 av_log(avctx, AV_LOG_WARNING, "Broken zTXt chunk\n");
1267             bytestream2_skip(&s->gb, length + 4);
1268             break;
1269         case MKTAG('s', 'T', 'E', 'R'): {
1270             int mode = bytestream2_get_byte(&s->gb);
1271             AVStereo3D *stereo3d = av_stereo3d_create_side_data(p);
1272             if (!stereo3d)
1273                 goto fail;
1274
1275             if (mode == 0 || mode == 1) {
1276                 stereo3d->type  = AV_STEREO3D_SIDEBYSIDE;
1277                 stereo3d->flags = mode ? 0 : AV_STEREO3D_FLAG_INVERT;
1278             } else {
1279                  av_log(avctx, AV_LOG_WARNING,
1280                         "Unknown value in sTER chunk (%d)\n", mode);
1281             }
1282             bytestream2_skip(&s->gb, 4); /* crc */
1283             break;
1284         }
1285         case MKTAG('i', 'C', 'C', 'P'): {
1286             if (decode_iccp_chunk(s, length, p) < 0)
1287                 goto fail;
1288             break;
1289         }
1290         case MKTAG('I', 'E', 'N', 'D'):
1291             if (!(s->pic_state & PNG_ALLIMAGE))
1292                 av_log(avctx, AV_LOG_ERROR, "IEND without all image\n");
1293             if (!(s->pic_state & (PNG_ALLIMAGE|PNG_IDAT))) {
1294                 ret = AVERROR_INVALIDDATA;
1295                 goto fail;
1296             }
1297             bytestream2_skip(&s->gb, 4); /* crc */
1298             goto exit_loop;
1299         default:
1300             /* skip tag */
1301 skip_tag:
1302             bytestream2_skip(&s->gb, length + 4);
1303             break;
1304         }
1305     }
1306 exit_loop:
1307
1308     if (avctx->codec_id == AV_CODEC_ID_PNG &&
1309         avctx->skip_frame == AVDISCARD_ALL) {
1310         return 0;
1311     }
1312
1313     if (s->bits_per_pixel <= 4)
1314         handle_small_bpp(s, p);
1315
1316     /* apply transparency if needed */
1317     if (s->has_trns && s->color_type != PNG_COLOR_TYPE_PALETTE) {
1318         size_t byte_depth = s->bit_depth > 8 ? 2 : 1;
1319         size_t raw_bpp = s->bpp - byte_depth;
1320         unsigned x, y;
1321
1322         av_assert0(s->bit_depth > 1);
1323
1324         for (y = 0; y < s->height; ++y) {
1325             uint8_t *row = &s->image_buf[s->image_linesize * y];
1326
1327             /* since we're updating in-place, we have to go from right to left */
1328             for (x = s->width; x > 0; --x) {
1329                 uint8_t *pixel = &row[s->bpp * (x - 1)];
1330                 memmove(pixel, &row[raw_bpp * (x - 1)], raw_bpp);
1331
1332                 if (!memcmp(pixel, s->transparent_color_be, raw_bpp)) {
1333                     memset(&pixel[raw_bpp], 0, byte_depth);
1334                 } else {
1335                     memset(&pixel[raw_bpp], 0xff, byte_depth);
1336                 }
1337             }
1338         }
1339     }
1340
1341     /* handle P-frames only if a predecessor frame is available */
1342     if (s->last_picture.f->data[0]) {
1343         if (   !(avpkt->flags & AV_PKT_FLAG_KEY) && avctx->codec_tag != AV_RL32("MPNG")
1344             && s->last_picture.f->width == p->width
1345             && s->last_picture.f->height== p->height
1346             && s->last_picture.f->format== p->format
1347          ) {
1348             if (CONFIG_PNG_DECODER && avctx->codec_id != AV_CODEC_ID_APNG)
1349                 handle_p_frame_png(s, p);
1350             else if (CONFIG_APNG_DECODER &&
1351                      avctx->codec_id == AV_CODEC_ID_APNG &&
1352                      (ret = handle_p_frame_apng(avctx, s, p)) < 0)
1353                 goto fail;
1354         }
1355     }
1356     ff_thread_report_progress(&s->picture, INT_MAX, 0);
1357     ff_thread_report_progress(&s->previous_picture, INT_MAX, 0);
1358
1359     return 0;
1360
1361 fail:
1362     ff_thread_report_progress(&s->picture, INT_MAX, 0);
1363     ff_thread_report_progress(&s->previous_picture, INT_MAX, 0);
1364     return ret;
1365 }
1366
1367 #if CONFIG_PNG_DECODER
1368 static int decode_frame_png(AVCodecContext *avctx,
1369                         void *data, int *got_frame,
1370                         AVPacket *avpkt)
1371 {
1372     PNGDecContext *const s = avctx->priv_data;
1373     const uint8_t *buf     = avpkt->data;
1374     int buf_size           = avpkt->size;
1375     AVFrame *p;
1376     int64_t sig;
1377     int ret;
1378
1379     ff_thread_release_buffer(avctx, &s->last_picture);
1380     FFSWAP(ThreadFrame, s->picture, s->last_picture);
1381     p = s->picture.f;
1382
1383     bytestream2_init(&s->gb, buf, buf_size);
1384
1385     /* check signature */
1386     sig = bytestream2_get_be64(&s->gb);
1387     if (sig != PNGSIG &&
1388         sig != MNGSIG) {
1389         av_log(avctx, AV_LOG_ERROR, "Invalid PNG signature 0x%08"PRIX64".\n", sig);
1390         return AVERROR_INVALIDDATA;
1391     }
1392
1393     s->y = s->has_trns = 0;
1394     s->hdr_state = 0;
1395     s->pic_state = 0;
1396
1397     /* init the zlib */
1398     s->zstream.zalloc = ff_png_zalloc;
1399     s->zstream.zfree  = ff_png_zfree;
1400     s->zstream.opaque = NULL;
1401     ret = inflateInit(&s->zstream);
1402     if (ret != Z_OK) {
1403         av_log(avctx, AV_LOG_ERROR, "inflateInit returned error %d\n", ret);
1404         return AVERROR_EXTERNAL;
1405     }
1406
1407     if ((ret = decode_frame_common(avctx, s, p, avpkt)) < 0)
1408         goto the_end;
1409
1410     if (avctx->skip_frame == AVDISCARD_ALL) {
1411         *got_frame = 0;
1412         ret = bytestream2_tell(&s->gb);
1413         goto the_end;
1414     }
1415
1416     if ((ret = av_frame_ref(data, s->picture.f)) < 0)
1417         goto the_end;
1418
1419     *got_frame = 1;
1420
1421     ret = bytestream2_tell(&s->gb);
1422 the_end:
1423     inflateEnd(&s->zstream);
1424     s->crow_buf = NULL;
1425     return ret;
1426 }
1427 #endif
1428
1429 #if CONFIG_APNG_DECODER
1430 static int decode_frame_apng(AVCodecContext *avctx,
1431                         void *data, int *got_frame,
1432                         AVPacket *avpkt)
1433 {
1434     PNGDecContext *const s = avctx->priv_data;
1435     int ret;
1436     AVFrame *p;
1437
1438     ff_thread_release_buffer(avctx, &s->last_picture);
1439     FFSWAP(ThreadFrame, s->picture, s->last_picture);
1440     p = s->picture.f;
1441
1442     if (!(s->hdr_state & PNG_IHDR)) {
1443         if (!avctx->extradata_size)
1444             return AVERROR_INVALIDDATA;
1445
1446         /* only init fields, there is no zlib use in extradata */
1447         s->zstream.zalloc = ff_png_zalloc;
1448         s->zstream.zfree  = ff_png_zfree;
1449
1450         bytestream2_init(&s->gb, avctx->extradata, avctx->extradata_size);
1451         if ((ret = decode_frame_common(avctx, s, p, avpkt)) < 0)
1452             goto end;
1453     }
1454
1455     /* reset state for a new frame */
1456     if ((ret = inflateInit(&s->zstream)) != Z_OK) {
1457         av_log(avctx, AV_LOG_ERROR, "inflateInit returned error %d\n", ret);
1458         ret = AVERROR_EXTERNAL;
1459         goto end;
1460     }
1461     s->y = 0;
1462     s->pic_state = 0;
1463     bytestream2_init(&s->gb, avpkt->data, avpkt->size);
1464     if ((ret = decode_frame_common(avctx, s, p, avpkt)) < 0)
1465         goto end;
1466
1467     if (!(s->pic_state & PNG_ALLIMAGE))
1468         av_log(avctx, AV_LOG_WARNING, "Frame did not contain a complete image\n");
1469     if (!(s->pic_state & (PNG_ALLIMAGE|PNG_IDAT))) {
1470         ret = AVERROR_INVALIDDATA;
1471         goto end;
1472     }
1473     if ((ret = av_frame_ref(data, s->picture.f)) < 0)
1474         goto end;
1475
1476     *got_frame = 1;
1477     ret = bytestream2_tell(&s->gb);
1478
1479 end:
1480     inflateEnd(&s->zstream);
1481     return ret;
1482 }
1483 #endif
1484
1485 #if HAVE_THREADS
1486 static int update_thread_context(AVCodecContext *dst, const AVCodecContext *src)
1487 {
1488     PNGDecContext *psrc = src->priv_data;
1489     PNGDecContext *pdst = dst->priv_data;
1490     int ret;
1491
1492     if (dst == src)
1493         return 0;
1494
1495     ff_thread_release_buffer(dst, &pdst->picture);
1496     if (psrc->picture.f->data[0] &&
1497         (ret = ff_thread_ref_frame(&pdst->picture, &psrc->picture)) < 0)
1498         return ret;
1499     if (CONFIG_APNG_DECODER && dst->codec_id == AV_CODEC_ID_APNG) {
1500         pdst->width             = psrc->width;
1501         pdst->height            = psrc->height;
1502         pdst->bit_depth         = psrc->bit_depth;
1503         pdst->color_type        = psrc->color_type;
1504         pdst->compression_type  = psrc->compression_type;
1505         pdst->interlace_type    = psrc->interlace_type;
1506         pdst->filter_type       = psrc->filter_type;
1507         pdst->cur_w = psrc->cur_w;
1508         pdst->cur_h = psrc->cur_h;
1509         pdst->x_offset = psrc->x_offset;
1510         pdst->y_offset = psrc->y_offset;
1511         pdst->has_trns = psrc->has_trns;
1512         memcpy(pdst->transparent_color_be, psrc->transparent_color_be, sizeof(pdst->transparent_color_be));
1513
1514         pdst->dispose_op = psrc->dispose_op;
1515
1516         memcpy(pdst->palette, psrc->palette, sizeof(pdst->palette));
1517
1518         pdst->hdr_state |= psrc->hdr_state;
1519
1520         ff_thread_release_buffer(dst, &pdst->last_picture);
1521         if (psrc->last_picture.f->data[0] &&
1522             (ret = ff_thread_ref_frame(&pdst->last_picture, &psrc->last_picture)) < 0)
1523             return ret;
1524
1525         ff_thread_release_buffer(dst, &pdst->previous_picture);
1526         if (psrc->previous_picture.f->data[0] &&
1527             (ret = ff_thread_ref_frame(&pdst->previous_picture, &psrc->previous_picture)) < 0)
1528             return ret;
1529     }
1530
1531     return 0;
1532 }
1533 #endif
1534
1535 static av_cold int png_dec_init(AVCodecContext *avctx)
1536 {
1537     PNGDecContext *s = avctx->priv_data;
1538
1539     avctx->color_range = AVCOL_RANGE_JPEG;
1540
1541     s->avctx = avctx;
1542     s->previous_picture.f = av_frame_alloc();
1543     s->last_picture.f = av_frame_alloc();
1544     s->picture.f = av_frame_alloc();
1545     if (!s->previous_picture.f || !s->last_picture.f || !s->picture.f) {
1546         av_frame_free(&s->previous_picture.f);
1547         av_frame_free(&s->last_picture.f);
1548         av_frame_free(&s->picture.f);
1549         return AVERROR(ENOMEM);
1550     }
1551
1552     if (!avctx->internal->is_copy) {
1553         avctx->internal->allocate_progress = 1;
1554         ff_pngdsp_init(&s->dsp);
1555     }
1556
1557     return 0;
1558 }
1559
1560 static av_cold int png_dec_end(AVCodecContext *avctx)
1561 {
1562     PNGDecContext *s = avctx->priv_data;
1563
1564     ff_thread_release_buffer(avctx, &s->previous_picture);
1565     av_frame_free(&s->previous_picture.f);
1566     ff_thread_release_buffer(avctx, &s->last_picture);
1567     av_frame_free(&s->last_picture.f);
1568     ff_thread_release_buffer(avctx, &s->picture);
1569     av_frame_free(&s->picture.f);
1570     av_freep(&s->buffer);
1571     s->buffer_size = 0;
1572     av_freep(&s->last_row);
1573     s->last_row_size = 0;
1574     av_freep(&s->tmp_row);
1575     s->tmp_row_size = 0;
1576
1577     return 0;
1578 }
1579
1580 #if CONFIG_APNG_DECODER
1581 AVCodec ff_apng_decoder = {
1582     .name           = "apng",
1583     .long_name      = NULL_IF_CONFIG_SMALL("APNG (Animated Portable Network Graphics) image"),
1584     .type           = AVMEDIA_TYPE_VIDEO,
1585     .id             = AV_CODEC_ID_APNG,
1586     .priv_data_size = sizeof(PNGDecContext),
1587     .init           = png_dec_init,
1588     .close          = png_dec_end,
1589     .decode         = decode_frame_apng,
1590     .init_thread_copy = ONLY_IF_THREADS_ENABLED(png_dec_init),
1591     .update_thread_context = ONLY_IF_THREADS_ENABLED(update_thread_context),
1592     .capabilities   = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_FRAME_THREADS /*| AV_CODEC_CAP_DRAW_HORIZ_BAND*/,
1593     .caps_internal  = FF_CODEC_CAP_INIT_THREADSAFE,
1594 };
1595 #endif
1596
1597 #if CONFIG_PNG_DECODER
1598 AVCodec ff_png_decoder = {
1599     .name           = "png",
1600     .long_name      = NULL_IF_CONFIG_SMALL("PNG (Portable Network Graphics) image"),
1601     .type           = AVMEDIA_TYPE_VIDEO,
1602     .id             = AV_CODEC_ID_PNG,
1603     .priv_data_size = sizeof(PNGDecContext),
1604     .init           = png_dec_init,
1605     .close          = png_dec_end,
1606     .decode         = decode_frame_png,
1607     .init_thread_copy = ONLY_IF_THREADS_ENABLED(png_dec_init),
1608     .update_thread_context = ONLY_IF_THREADS_ENABLED(update_thread_context),
1609     .capabilities   = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_FRAME_THREADS /*| AV_CODEC_CAP_DRAW_HORIZ_BAND*/,
1610     .caps_internal  = FF_CODEC_CAP_SKIP_FRAME_FILL_PARAM | FF_CODEC_CAP_INIT_THREADSAFE,
1611 };
1612 #endif