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

   1 /*
   2  * Copyright (c) 2006 Konstantin Shishkov
   3  *
   4  * This file is part of FFmpeg.
   5  *
   6  * FFmpeg is free software; you can redistribute it and/or
   7  * modify it under the terms of the GNU Lesser General Public
   8  * License as published by the Free Software Foundation; either
   9  * version 2.1 of the License, or (at your option) any later version.
  10  *
  11  * FFmpeg is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14  * Lesser General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU Lesser General Public
  17  * License along with FFmpeg; if not, write to the Free Software
  18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  19  */
  20
  21 /**
  22  * @file
  23  * TIFF image decoder
  24  * @author Konstantin Shishkov
  25  */
  26
  27 #include "config.h"
  28 #if CONFIG_ZLIB
  29 #include <zlib.h>
  30 #endif
  31 #if CONFIG_LZMA
  32 #define LZMA_API_STATIC
  33 #include <lzma.h>
  34 #endif
  35
  36 #include "libavutil/attributes.h"
  37 #include "libavutil/avstring.h"
  38 #include "libavutil/error.h"
  39 #include "libavutil/intreadwrite.h"
  40 #include "libavutil/imgutils.h"
  41 #include "libavutil/opt.h"
  42 #include "avcodec.h"
  43 #include "bytestream.h"
  44 #include "faxcompr.h"
  45 #include "internal.h"
  46 #include "lzw.h"
  47 #include "mathops.h"
  48 #include "tiff.h"
  49 #include "tiff_data.h"
  50 #include "mjpegdec.h"
  51 #include "thread.h"
  52 #include "get_bits.h"
  53
  54 typedef struct TiffContext {
  55     AVClass *class;
  56     AVCodecContext *avctx;
  57     GetByteContext gb;
  58
  59     /* JPEG decoding for DNG */
  60     AVCodecContext *avctx_mjpeg; // wrapper context for MJPEG
  61     AVFrame *jpgframe;           // decoded JPEG tile
  62
  63     int get_subimage;
  64     uint16_t get_page;
  65     int get_thumbnail;
  66
  67     enum TiffType tiff_type;
  68     int width, height;
  69     unsigned int bpp, bppcount;
  70     uint32_t palette[256];
  71     int palette_is_set;
  72     int le;
  73     enum TiffCompr compr;
  74     enum TiffPhotometric photometric;
  75     int planar;
  76     int subsampling[2];
  77     int fax_opts;
  78     int predictor;
  79     int fill_order;
  80     uint32_t res[4];
  81     int is_thumbnail;
  82
  83     int is_bayer;
  84     uint8_t pattern[4];
  85     unsigned black_level;
  86     unsigned white_level;
  87     const uint16_t *dng_lut; // Pointer to DNG linearization table
  88
  89     uint32_t sub_ifd;
  90     uint16_t cur_page;
  91
  92     int strips, rps, sstype;
  93     int sot;
  94     int stripsizesoff, stripsize, stripoff, strippos;
  95     LZWState *lzw;
  96
  97     /* Tile support */
  98     int is_tiled;
  99     int tile_byte_counts_offset, tile_offsets_offset;
 100     int tile_width, tile_length;
 101     int tile_count;
 102
 103     int is_jpeg;
 104
 105     uint8_t *deinvert_buf;
 106     int deinvert_buf_size;
 107     uint8_t *yuv_line;
 108     unsigned int yuv_line_size;
 109     uint8_t *fax_buffer;
 110     unsigned int fax_buffer_size;
 111
 112     int geotag_count;
 113     TiffGeoTag *geotags;
 114 } TiffContext;
 115
 116 static void tiff_set_type(TiffContext *s, enum TiffType tiff_type) {
 117     if (s->tiff_type < tiff_type) // Prioritize higher-valued entries
 118         s->tiff_type = tiff_type;
 119 }
 120
 121 static void free_geotags(TiffContext *const s)
 122 {
 123     int i;
 124     for (i = 0; i < s->geotag_count; i++) {
 125         if (s->geotags[i].val)
 126             av_freep(&s->geotags[i].val);
 127     }
 128     av_freep(&s->geotags);
 129     s->geotag_count = 0;
 130 }
 131
 132 #define RET_GEOKEY(TYPE, array, element)\
 133     if (key >= TIFF_##TYPE##_KEY_ID_OFFSET &&\
 134         key - TIFF_##TYPE##_KEY_ID_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_name_type_map))\
 135         return ff_tiff_##array##_name_type_map[key - TIFF_##TYPE##_KEY_ID_OFFSET].element;
 136
 137 static const char *get_geokey_name(int key)
 138 {
 139     RET_GEOKEY(VERT, vert, name);
 140     RET_GEOKEY(PROJ, proj, name);
 141     RET_GEOKEY(GEOG, geog, name);
 142     RET_GEOKEY(CONF, conf, name);
 143
 144     return NULL;
 145 }
 146
 147 static int get_geokey_type(int key)
 148 {
 149     RET_GEOKEY(VERT, vert, type);
 150     RET_GEOKEY(PROJ, proj, type);
 151     RET_GEOKEY(GEOG, geog, type);
 152     RET_GEOKEY(CONF, conf, type);
 153
 154     return AVERROR_INVALIDDATA;
 155 }
 156
 157 static int cmp_id_key(const void *id, const void *k)
 158 {
 159     return *(const int*)id - ((const TiffGeoTagKeyName*)k)->key;
 160 }
 161
 162 static const char *search_keyval(const TiffGeoTagKeyName *keys, int n, int id)
 163 {
 164     TiffGeoTagKeyName *r = bsearch(&id, keys, n, sizeof(keys[0]), cmp_id_key);
 165     if(r)
 166         return r->name;
 167
 168     return NULL;
 169 }
 170
 171 static char *get_geokey_val(int key, int val)
 172 {
 173     char *ap;
 174
 175     if (val == TIFF_GEO_KEY_UNDEFINED)
 176         return av_strdup("undefined");
 177     if (val == TIFF_GEO_KEY_USER_DEFINED)
 178         return av_strdup("User-Defined");
 179
 180 #define RET_GEOKEY_VAL(TYPE, array)\
 181     if (val >= TIFF_##TYPE##_OFFSET &&\
 182         val - TIFF_##TYPE##_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_codes))\
 183         return av_strdup(ff_tiff_##array##_codes[val - TIFF_##TYPE##_OFFSET]);
 184
 185     switch (key) {
 186     case TIFF_GT_MODEL_TYPE_GEOKEY:
 187         RET_GEOKEY_VAL(GT_MODEL_TYPE, gt_model_type);
 188         break;
 189     case TIFF_GT_RASTER_TYPE_GEOKEY:
 190         RET_GEOKEY_VAL(GT_RASTER_TYPE, gt_raster_type);
 191         break;
 192     case TIFF_GEOG_LINEAR_UNITS_GEOKEY:
 193     case TIFF_PROJ_LINEAR_UNITS_GEOKEY:
 194     case TIFF_VERTICAL_UNITS_GEOKEY:
 195         RET_GEOKEY_VAL(LINEAR_UNIT, linear_unit);
 196         break;
 197     case TIFF_GEOG_ANGULAR_UNITS_GEOKEY:
 198     case TIFF_GEOG_AZIMUTH_UNITS_GEOKEY:
 199         RET_GEOKEY_VAL(ANGULAR_UNIT, angular_unit);
 200         break;
 201     case TIFF_GEOGRAPHIC_TYPE_GEOKEY:
 202         RET_GEOKEY_VAL(GCS_TYPE, gcs_type);
 203         RET_GEOKEY_VAL(GCSE_TYPE, gcse_type);
 204         break;
 205     case TIFF_GEOG_GEODETIC_DATUM_GEOKEY:
 206         RET_GEOKEY_VAL(GEODETIC_DATUM, geodetic_datum);
 207         RET_GEOKEY_VAL(GEODETIC_DATUM_E, geodetic_datum_e);
 208         break;
 209     case TIFF_GEOG_ELLIPSOID_GEOKEY:
 210         RET_GEOKEY_VAL(ELLIPSOID, ellipsoid);
 211         break;
 212     case TIFF_GEOG_PRIME_MERIDIAN_GEOKEY:
 213         RET_GEOKEY_VAL(PRIME_MERIDIAN, prime_meridian);
 214         break;
 215     case TIFF_PROJECTED_CS_TYPE_GEOKEY:
 216         ap = av_strdup(search_keyval(ff_tiff_proj_cs_type_codes, FF_ARRAY_ELEMS(ff_tiff_proj_cs_type_codes), val));
 217         if(ap) return ap;
 218         break;
 219     case TIFF_PROJECTION_GEOKEY:
 220         ap = av_strdup(search_keyval(ff_tiff_projection_codes, FF_ARRAY_ELEMS(ff_tiff_projection_codes), val));
 221         if(ap) return ap;
 222         break;
 223     case TIFF_PROJ_COORD_TRANS_GEOKEY:
 224         RET_GEOKEY_VAL(COORD_TRANS, coord_trans);
 225         break;
 226     case TIFF_VERTICAL_CS_TYPE_GEOKEY:
 227         RET_GEOKEY_VAL(VERT_CS, vert_cs);
 228         RET_GEOKEY_VAL(ORTHO_VERT_CS, ortho_vert_cs);
 229         break;
 230
 231     }
 232
 233     ap = av_malloc(14);
 234     if (ap)
 235         snprintf(ap, 14, "Unknown-%d", val);
 236     return ap;
 237 }
 238
 239 static char *doubles2str(double *dp, int count, const char *sep)
 240 {
 241     int i;
 242     char *ap, *ap0;
 243     uint64_t component_len;
 244     if (!sep) sep = ", ";
 245     component_len = 24LL + strlen(sep);
 246     if (count >= (INT_MAX - 1)/component_len)
 247         return NULL;
 248     ap = av_malloc(component_len * count + 1);
 249     if (!ap)
 250         return NULL;
 251     ap0   = ap;
 252     ap[0] = '\0';
 253     for (i = 0; i < count; i++) {
 254         unsigned l = snprintf(ap, component_len, "%.15g%s", dp[i], sep);
 255         if(l >= component_len) {
 256             av_free(ap0);
 257             return NULL;
 258         }
 259         ap += l;
 260     }
 261     ap0[strlen(ap0) - strlen(sep)] = '\0';
 262     return ap0;
 263 }
 264
 265 static int add_metadata(int count, int type,
 266                         const char *name, const char *sep, TiffContext *s, AVFrame *frame)
 267 {
 268     switch(type) {
 269     case TIFF_DOUBLE: return ff_tadd_doubles_metadata(count, name, sep, &s->gb, s->le, &frame->metadata);
 270     case TIFF_SHORT : return ff_tadd_shorts_metadata(count, name, sep, &s->gb, s->le, 0, &frame->metadata);
 271     case TIFF_STRING: return ff_tadd_string_metadata(count, name, &s->gb, s->le, &frame->metadata);
 272     default         : return AVERROR_INVALIDDATA;
 273     };
 274 }
 275
 276 static void av_always_inline dng_blit(TiffContext *s, uint8_t *dst, int dst_stride,
 277                                       const uint8_t *src, int src_stride, int width, int height, int is_u16);
 278
 279 static void av_always_inline horizontal_fill(TiffContext *s,
 280                                              unsigned int bpp, uint8_t* dst,
 281                                              int usePtr, const uint8_t *src,
 282                                              uint8_t c, int width, int offset)
 283 {
 284     switch (bpp) {
 285     case 1:
 286         while (--width >= 0) {
 287             dst[(width+offset)*8+7] = (usePtr ? src[width] : c)      & 0x1;
 288             dst[(width+offset)*8+6] = (usePtr ? src[width] : c) >> 1 & 0x1;
 289             dst[(width+offset)*8+5] = (usePtr ? src[width] : c) >> 2 & 0x1;
 290             dst[(width+offset)*8+4] = (usePtr ? src[width] : c) >> 3 & 0x1;
 291             dst[(width+offset)*8+3] = (usePtr ? src[width] : c) >> 4 & 0x1;
 292             dst[(width+offset)*8+2] = (usePtr ? src[width] : c) >> 5 & 0x1;
 293             dst[(width+offset)*8+1] = (usePtr ? src[width] : c) >> 6 & 0x1;
 294             dst[(width+offset)*8+0] = (usePtr ? src[width] : c) >> 7;
 295         }
 296         break;
 297     case 2:
 298         while (--width >= 0) {
 299             dst[(width+offset)*4+3] = (usePtr ? src[width] : c) & 0x3;
 300             dst[(width+offset)*4+2] = (usePtr ? src[width] : c) >> 2 & 0x3;
 301             dst[(width+offset)*4+1] = (usePtr ? src[width] : c) >> 4 & 0x3;
 302             dst[(width+offset)*4+0] = (usePtr ? src[width] : c) >> 6;
 303         }
 304         break;
 305     case 4:
 306         while (--width >= 0) {
 307             dst[(width+offset)*2+1] = (usePtr ? src[width] : c) & 0xF;
 308             dst[(width+offset)*2+0] = (usePtr ? src[width] : c) >> 4;
 309         }
 310         break;
 311     case 12: {
 312                  uint16_t *dst16 = (uint16_t *)dst;
 313                  GetBitContext gb;
 314                  init_get_bits8(&gb, src, width);
 315                  for (int i = 0; i < s->width; i++) {
 316                      dst16[i] = get_bits(&gb, 12) << 4;
 317                  }
 318              }
 319         break;
 320     default:
 321         if (usePtr) {
 322             memcpy(dst + offset, src, width);
 323         } else {
 324             memset(dst + offset, c, width);
 325         }
 326     }
 327 }
 328
 329 static int deinvert_buffer(TiffContext *s, const uint8_t *src, int size)
 330 {
 331     int i;
 332
 333     av_fast_padded_malloc(&s->deinvert_buf, &s->deinvert_buf_size, size);
 334     if (!s->deinvert_buf)
 335         return AVERROR(ENOMEM);
 336     for (i = 0; i < size; i++)
 337         s->deinvert_buf[i] = ff_reverse[src[i]];
 338
 339     return 0;
 340 }
 341
 342 static void unpack_gray(TiffContext *s, AVFrame *p,
 343                        const uint8_t *src, int lnum, int width, int bpp)
 344 {
 345     GetBitContext gb;
 346     uint16_t *dst = (uint16_t *)(p->data[0] + lnum * p->linesize[0]);
 347
 348     init_get_bits8(&gb, src, width);
 349
 350     for (int i = 0; i < s->width; i++) {
 351         dst[i] = get_bits(&gb, bpp);
 352     }
 353 }
 354
 355 static void unpack_yuv(TiffContext *s, AVFrame *p,
 356                        const uint8_t *src, int lnum)
 357 {
 358     int i, j, k;
 359     int w       = (s->width - 1) / s->subsampling[0] + 1;
 360     uint8_t *pu = &p->data[1][lnum / s->subsampling[1] * p->linesize[1]];
 361     uint8_t *pv = &p->data[2][lnum / s->subsampling[1] * p->linesize[2]];
 362     if (s->width % s->subsampling[0] || s->height % s->subsampling[1]) {
 363         for (i = 0; i < w; i++) {
 364             for (j = 0; j < s->subsampling[1]; j++)
 365                 for (k = 0; k < s->subsampling[0]; k++)
 366                     p->data[0][FFMIN(lnum + j, s->height-1) * p->linesize[0] +
 367                                FFMIN(i * s->subsampling[0] + k, s->width-1)] = *src++;
 368             *pu++ = *src++;
 369             *pv++ = *src++;
 370         }
 371     }else{
 372         for (i = 0; i < w; i++) {
 373             for (j = 0; j < s->subsampling[1]; j++)
 374                 for (k = 0; k < s->subsampling[0]; k++)
 375                     p->data[0][(lnum + j) * p->linesize[0] +
 376                                i * s->subsampling[0] + k] = *src++;
 377             *pu++ = *src++;
 378             *pv++ = *src++;
 379         }
 380     }
 381 }
 382
 383 #if CONFIG_ZLIB
 384 static int tiff_uncompress(uint8_t *dst, unsigned long *len, const uint8_t *src,
 385                            int size)
 386 {
 387     z_stream zstream = { 0 };
 388     int zret;
 389
 390     zstream.next_in   = (uint8_t *)src;
 391     zstream.avail_in  = size;
 392     zstream.next_out  = dst;
 393     zstream.avail_out = *len;
 394     zret              = inflateInit(&zstream);
 395     if (zret != Z_OK) {
 396         av_log(NULL, AV_LOG_ERROR, "Inflate init error: %d\n", zret);
 397         return zret;
 398     }
 399     zret = inflate(&zstream, Z_SYNC_FLUSH);
 400     inflateEnd(&zstream);
 401     *len = zstream.total_out;
 402     return zret == Z_STREAM_END ? Z_OK : zret;
 403 }
 404
 405 static int tiff_unpack_zlib(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
 406                             const uint8_t *src, int size, int width, int lines,
 407                             int strip_start, int is_yuv)
 408 {
 409     uint8_t *zbuf;
 410     unsigned long outlen;
 411     int ret, line;
 412     outlen = width * lines;
 413     zbuf   = av_malloc(outlen);
 414     if (!zbuf)
 415         return AVERROR(ENOMEM);
 416     if (s->fill_order) {
 417         if ((ret = deinvert_buffer(s, src, size)) < 0) {
 418             av_free(zbuf);
 419             return ret;
 420         }
 421         src = s->deinvert_buf;
 422     }
 423     ret = tiff_uncompress(zbuf, &outlen, src, size);
 424     if (ret != Z_OK) {
 425         av_log(s->avctx, AV_LOG_ERROR,
 426                "Uncompressing failed (%lu of %lu) with error %d\n", outlen,
 427                (unsigned long)width * lines, ret);
 428         av_free(zbuf);
 429         return AVERROR_UNKNOWN;
 430     }
 431     src = zbuf;
 432     for (line = 0; line < lines; line++) {
 433         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
 434             horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
 435         } else {
 436             memcpy(dst, src, width);
 437         }
 438         if (is_yuv) {
 439             unpack_yuv(s, p, dst, strip_start + line);
 440             line += s->subsampling[1] - 1;
 441         }
 442         dst += stride;
 443         src += width;
 444     }
 445     av_free(zbuf);
 446     return 0;
 447 }
 448 #endif
 449
 450 #if CONFIG_LZMA
 451 static int tiff_uncompress_lzma(uint8_t *dst, uint64_t *len, const uint8_t *src,
 452                                 int size)
 453 {
 454     lzma_stream stream = LZMA_STREAM_INIT;
 455     lzma_ret ret;
 456
 457     stream.next_in   = (uint8_t *)src;
 458     stream.avail_in  = size;
 459     stream.next_out  = dst;
 460     stream.avail_out = *len;
 461     ret              = lzma_stream_decoder(&stream, UINT64_MAX, 0);
 462     if (ret != LZMA_OK) {
 463         av_log(NULL, AV_LOG_ERROR, "LZMA init error: %d\n", ret);
 464         return ret;
 465     }
 466     ret = lzma_code(&stream, LZMA_RUN);
 467     lzma_end(&stream);
 468     *len = stream.total_out;
 469     return ret == LZMA_STREAM_END ? LZMA_OK : ret;
 470 }
 471
 472 static int tiff_unpack_lzma(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
 473                             const uint8_t *src, int size, int width, int lines,
 474                             int strip_start, int is_yuv)
 475 {
 476     uint64_t outlen = width * (uint64_t)lines;
 477     int ret, line;
 478     uint8_t *buf = av_malloc(outlen);
 479     if (!buf)
 480         return AVERROR(ENOMEM);
 481     if (s->fill_order) {
 482         if ((ret = deinvert_buffer(s, src, size)) < 0) {
 483             av_free(buf);
 484             return ret;
 485         }
 486         src = s->deinvert_buf;
 487     }
 488     ret = tiff_uncompress_lzma(buf, &outlen, src, size);
 489     if (ret != LZMA_OK) {
 490         av_log(s->avctx, AV_LOG_ERROR,
 491                "Uncompressing failed (%"PRIu64" of %"PRIu64") with error %d\n", outlen,
 492                (uint64_t)width * lines, ret);
 493         av_free(buf);
 494         return AVERROR_UNKNOWN;
 495     }
 496     src = buf;
 497     for (line = 0; line < lines; line++) {
 498         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
 499             horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
 500         } else {
 501             memcpy(dst, src, width);
 502         }
 503         if (is_yuv) {
 504             unpack_yuv(s, p, dst, strip_start + line);
 505             line += s->subsampling[1] - 1;
 506         }
 507         dst += stride;
 508         src += width;
 509     }
 510     av_free(buf);
 511     return 0;
 512 }
 513 #endif
 514
 515 static int tiff_unpack_fax(TiffContext *s, uint8_t *dst, int stride,
 516                            const uint8_t *src, int size, int width, int lines)
 517 {
 518     int i, ret = 0;
 519     int line;
 520     uint8_t *src2;
 521
 522     av_fast_padded_malloc(&s->fax_buffer, &s->fax_buffer_size, size);
 523     src2 = s->fax_buffer;
 524
 525     if (!src2) {
 526         av_log(s->avctx, AV_LOG_ERROR,
 527                "Error allocating temporary buffer\n");
 528         return AVERROR(ENOMEM);
 529     }
 530
 531     if (!s->fill_order) {
 532         memcpy(src2, src, size);
 533     } else {
 534         for (i = 0; i < size; i++)
 535             src2[i] = ff_reverse[src[i]];
 536     }
 537     memset(src2 + size, 0, AV_INPUT_BUFFER_PADDING_SIZE);
 538     ret = ff_ccitt_unpack(s->avctx, src2, size, dst, lines, stride,
 539                           s->compr, s->fax_opts);
 540     if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
 541         for (line = 0; line < lines; line++) {
 542             horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
 543             dst += stride;
 544         }
 545     return ret;
 546 }
 547
 548 static int tiff_unpack_strip(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
 549                              const uint8_t *src, int size, int strip_start, int lines)
 550 {
 551     PutByteContext pb;
 552     int c, line, pixels, code, ret;
 553     const uint8_t *ssrc = src;
 554     int width = ((s->width * s->bpp) + 7) >> 3;
 555     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(p->format);
 556     int is_yuv = !(desc->flags & AV_PIX_FMT_FLAG_RGB) &&
 557                  (desc->flags & AV_PIX_FMT_FLAG_PLANAR) &&
 558                  desc->nb_components >= 3;
 559     int is_dng;
 560
 561     if (s->planar)
 562         width /= s->bppcount;
 563
 564     if (size <= 0)
 565         return AVERROR_INVALIDDATA;
 566
 567     if (is_yuv) {
 568         int bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
 569                             s->subsampling[0] * s->subsampling[1] + 7) >> 3;
 570         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
 571         if (s->yuv_line == NULL) {
 572             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
 573             return AVERROR(ENOMEM);
 574         }
 575         dst = s->yuv_line;
 576         stride = 0;
 577
 578         width = (s->width - 1) / s->subsampling[0] + 1;
 579         width = width * s->subsampling[0] * s->subsampling[1] + 2*width;
 580         av_assert0(width <= bytes_per_row);
 581         av_assert0(s->bpp == 24);
 582     }
 583     if (s->is_bayer) {
 584         width = (s->bpp * s->width + 7) >> 3;
 585     }
 586     if (p->format == AV_PIX_FMT_GRAY12) {
 587         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, width);
 588         if (s->yuv_line == NULL) {
 589             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
 590             return AVERROR(ENOMEM);
 591         }
 592         dst = s->yuv_line;
 593         stride = 0;
 594     }
 595
 596     if (s->compr == TIFF_DEFLATE || s->compr == TIFF_ADOBE_DEFLATE) {
 597 #if CONFIG_ZLIB
 598         return tiff_unpack_zlib(s, p, dst, stride, src, size, width, lines,
 599                                 strip_start, is_yuv);
 600 #else
 601         av_log(s->avctx, AV_LOG_ERROR,
 602                "zlib support not enabled, "
 603                "deflate compression not supported\n");
 604         return AVERROR(ENOSYS);
 605 #endif
 606     }
 607     if (s->compr == TIFF_LZMA) {
 608 #if CONFIG_LZMA
 609         return tiff_unpack_lzma(s, p, dst, stride, src, size, width, lines,
 610                                 strip_start, is_yuv);
 611 #else
 612         av_log(s->avctx, AV_LOG_ERROR,
 613                "LZMA support not enabled\n");
 614         return AVERROR(ENOSYS);
 615 #endif
 616     }
 617     if (s->compr == TIFF_LZW) {
 618         if (s->fill_order) {
 619             if ((ret = deinvert_buffer(s, src, size)) < 0)
 620                 return ret;
 621             ssrc = src = s->deinvert_buf;
 622         }
 623         if (size > 1 && !src[0] && (src[1]&1)) {
 624             av_log(s->avctx, AV_LOG_ERROR, "Old style LZW is unsupported\n");
 625         }
 626         if ((ret = ff_lzw_decode_init(s->lzw, 8, src, size, FF_LZW_TIFF)) < 0) {
 627             av_log(s->avctx, AV_LOG_ERROR, "Error initializing LZW decoder\n");
 628             return ret;
 629         }
 630         for (line = 0; line < lines; line++) {
 631             pixels = ff_lzw_decode(s->lzw, dst, width);
 632             if (pixels < width) {
 633                 av_log(s->avctx, AV_LOG_ERROR, "Decoded only %i bytes of %i\n",
 634                        pixels, width);
 635                 return AVERROR_INVALIDDATA;
 636             }
 637             if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
 638                 horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
 639             if (is_yuv) {
 640                 unpack_yuv(s, p, dst, strip_start + line);
 641                 line += s->subsampling[1] - 1;
 642             } else if (p->format == AV_PIX_FMT_GRAY12) {
 643                 unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
 644             }
 645             dst += stride;
 646         }
 647         return 0;
 648     }
 649     if (s->compr == TIFF_CCITT_RLE ||
 650         s->compr == TIFF_G3        ||
 651         s->compr == TIFF_G4) {
 652         if (is_yuv || p->format == AV_PIX_FMT_GRAY12)
 653             return AVERROR_INVALIDDATA;
 654
 655         return tiff_unpack_fax(s, dst, stride, src, size, width, lines);
 656     }
 657
 658     bytestream2_init(&s->gb, src, size);
 659     bytestream2_init_writer(&pb, dst, is_yuv ? s->yuv_line_size : (stride * lines));
 660
 661     is_dng = (s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG);
 662
 663     for (line = 0; line < lines; line++) {
 664         if (src - ssrc > size) {
 665             av_log(s->avctx, AV_LOG_ERROR, "Source data overread\n");
 666             return AVERROR_INVALIDDATA;
 667         }
 668
 669         if (bytestream2_get_bytes_left(&s->gb) == 0 || bytestream2_get_eof(&pb))
 670             break;
 671         bytestream2_seek_p(&pb, stride * line, SEEK_SET);
 672         switch (s->compr) {
 673         case TIFF_RAW:
 674             if (ssrc + size - src < width)
 675                 return AVERROR_INVALIDDATA;
 676
 677             if (!s->fill_order) {
 678                 horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 || s->is_bayer),
 679                                 dst, 1, src, 0, width, 0);
 680             } else {
 681                 int i;
 682                 for (i = 0; i < width; i++)
 683                     dst[i] = ff_reverse[src[i]];
 684             }
 685
 686             /* Color processing for DNG images with uncompressed strips (non-tiled) */
 687             if (is_dng) {
 688                 int is_u16, pixel_size_bytes, pixel_size_bits;
 689
 690                 is_u16 = (s->bpp > 8);
 691                 pixel_size_bits = (is_u16 ? 16 : 8);
 692                 pixel_size_bytes = (is_u16 ? sizeof(uint16_t) : sizeof(uint8_t));
 693
 694                 dng_blit(s,
 695                          dst,
 696                          0, // no stride, only 1 line
 697                          dst,
 698                          0, // no stride, only 1 line
 699                          width / pixel_size_bytes * pixel_size_bits / s->bpp * s->bppcount, // need to account for [1, 16] bpp
 700                          1,
 701                          is_u16);
 702             }
 703
 704             src += width;
 705             break;
 706         case TIFF_PACKBITS:
 707             for (pixels = 0; pixels < width;) {
 708                 if (ssrc + size - src < 2) {
 709                     av_log(s->avctx, AV_LOG_ERROR, "Read went out of bounds\n");
 710                     return AVERROR_INVALIDDATA;
 711                 }
 712                 code = s->fill_order ? (int8_t) ff_reverse[*src++]: (int8_t) *src++;
 713                 if (code >= 0) {
 714                     code++;
 715                     if (pixels + code > width ||
 716                         ssrc + size - src < code) {
 717                         av_log(s->avctx, AV_LOG_ERROR,
 718                                "Copy went out of bounds\n");
 719                         return AVERROR_INVALIDDATA;
 720                     }
 721                     horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
 722                                     dst, 1, src, 0, code, pixels);
 723                     src    += code;
 724                     pixels += code;
 725                 } else if (code != -128) { // -127..-1
 726                     code = (-code) + 1;
 727                     if (pixels + code > width) {
 728                         av_log(s->avctx, AV_LOG_ERROR,
 729                                "Run went out of bounds\n");
 730                         return AVERROR_INVALIDDATA;
 731                     }
 732                     c = *src++;
 733                     horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
 734                                     dst, 0, NULL, c, code, pixels);
 735                     pixels += code;
 736                 }
 737             }
 738             if (s->fill_order) {
 739                 int i;
 740                 for (i = 0; i < width; i++)
 741                     dst[i] = ff_reverse[dst[i]];
 742             }
 743             break;
 744         }
 745         if (is_yuv) {
 746             unpack_yuv(s, p, dst, strip_start + line);
 747             line += s->subsampling[1] - 1;
 748         } else if (p->format == AV_PIX_FMT_GRAY12) {
 749             unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
 750         }
 751         dst += stride;
 752     }
 753     return 0;
 754 }
 755
 756 static float av_always_inline linear_to_srgb(float value) {
 757     if (value <= 0.0031308f)
 758         return value * 12.92f;
 759     else
 760         return powf(value * 1.055f, 1.0f / 2.4f) - 0.055f;
 761 }
 762
 763 /**
 764  * Map stored raw sensor values into linear reference values (see: DNG Specification - Chapter 5)
 765  * Then convert to sRGB color space.
 766  */
 767 static uint16_t av_always_inline dng_process_color16(uint16_t value,
 768                                                      const uint16_t *lut,
 769                                                      uint16_t black_level,
 770                                                      float scale_factor) {
 771     float value_norm;
 772
 773     // Lookup table lookup
 774     if (lut)
 775         value = lut[value];
 776
 777     // Black level subtraction
 778     value = av_clip_uint16_c((unsigned)value - black_level);
 779
 780     // Color scaling
 781     value_norm = (float)value * scale_factor;
 782
 783     // Color space conversion (sRGB)
 784     value = av_clip_uint16_c((uint16_t)(linear_to_srgb(value_norm) * 0xFFFF));
 785
 786     return value;
 787 }
 788
 789 static uint16_t av_always_inline dng_process_color8(uint16_t value,
 790                                                     const uint16_t *lut,
 791                                                     uint16_t black_level,
 792                                                     float scale_factor) {
 793     return dng_process_color16(value, lut, black_level, scale_factor) >> 8;
 794 }
 795
 796 static void dng_blit(TiffContext *s, uint8_t *dst, int dst_stride,
 797                      const uint8_t *src, int src_stride,
 798                      int width, int height, int is_u16)
 799 {
 800     int line, col;
 801     float scale_factor;
 802
 803     scale_factor = 1.0f / (s->white_level - s->black_level);
 804
 805     if (is_u16) {
 806         for (line = 0; line < height; line++) {
 807             uint16_t *dst_u16 = (uint16_t *)dst;
 808             uint16_t *src_u16 = (uint16_t *)src;
 809
 810             for (col = 0; col < width; col++)
 811                 *dst_u16++ = dng_process_color16(*src_u16++, s->dng_lut, s->black_level, scale_factor);
 812
 813             dst += dst_stride * sizeof(uint16_t);
 814             src += src_stride * sizeof(uint16_t);
 815         }
 816     } else {
 817         for (line = 0; line < height; line++) {
 818             for (col = 0; col < width; col++)
 819                 *dst++ = dng_process_color8(*src++, s->dng_lut, s->black_level, scale_factor);
 820
 821             dst += dst_stride;
 822             src += src_stride;
 823         }
 824     }
 825 }
 826
 827 static int dng_decode_jpeg_tile(AVCodecContext *avctx, AVFrame *frame,
 828                                 int tile_byte_count, int x, int y, int w, int h)
 829 {
 830     TiffContext *s = avctx->priv_data;
 831     AVPacket jpkt;
 832     uint8_t *dst_data, *src_data;
 833     uint32_t dst_offset; /* offset from dst buffer in pixels */
 834     int is_u16, pixel_size;
 835     int ret;
 836
 837     /* Prepare a packet and send to the MJPEG decoder */
 838     av_init_packet(&jpkt);
 839     jpkt.data = (uint8_t*)s->gb.buffer;
 840     jpkt.size = tile_byte_count;
 841
 842     if (s->is_bayer) {
 843         MJpegDecodeContext *mjpegdecctx = s->avctx_mjpeg->priv_data;
 844         /* We have to set this information here, there is no way to know if a given JPEG is a DNG-embedded
 845            image or not from its own data (and we need that information when decoding it). */
 846         mjpegdecctx->bayer = 1;
 847     }
 848
 849     ret = avcodec_send_packet(s->avctx_mjpeg, &jpkt);
 850     if (ret < 0) {
 851         av_log(avctx, AV_LOG_ERROR, "Error submitting a packet for decoding\n");
 852         return ret;
 853     }
 854
 855     ret = avcodec_receive_frame(s->avctx_mjpeg, s->jpgframe);
 856     if (ret < 0) {
 857         av_log(avctx, AV_LOG_ERROR, "JPEG decoding error: %s.\n", av_err2str(ret));
 858
 859         /* Normally skip, error if explode */
 860         if (avctx->err_recognition & AV_EF_EXPLODE)
 861             return AVERROR_INVALIDDATA;
 862         else
 863             return 0;
 864     }
 865
 866     /* Copy the outputted tile's pixels from 'jpgframe' to 'frame' (final buffer) */
 867
 868     is_u16 = (s->bpp > 8);
 869     pixel_size = (is_u16 ? sizeof(uint16_t) : sizeof(uint8_t));
 870
 871     dst_offset = x + frame->linesize[0] * y / pixel_size;
 872     dst_data = frame->data[0] + dst_offset * pixel_size;
 873     src_data = s->jpgframe->data[0];
 874
 875     dng_blit(s,
 876              dst_data,
 877              frame->linesize[0] / pixel_size,
 878              src_data,
 879              s->jpgframe->linesize[0] / pixel_size,
 880              w,
 881              h,
 882              is_u16);
 883
 884     av_frame_unref(s->jpgframe);
 885
 886     return 0;
 887 }
 888
 889 static int dng_decode_tiles(AVCodecContext *avctx, AVFrame *frame)
 890 {
 891     TiffContext *s = avctx->priv_data;
 892     int tile_idx;
 893     int tile_offset_offset, tile_offset;
 894     int tile_byte_count_offset, tile_byte_count;
 895     int tile_count_x, tile_count_y;
 896     int tile_width, tile_length;
 897     int has_width_leftover, has_height_leftover;
 898     int tile_x = 0, tile_y = 0;
 899     int pos_x = 0, pos_y = 0;
 900     int ret;
 901
 902     has_width_leftover = (s->width % s->tile_width != 0);
 903     has_height_leftover = (s->height % s->tile_length != 0);
 904
 905     /* Calculate tile counts (round up) */
 906     tile_count_x = (s->width + s->tile_width - 1) / s->tile_width;
 907     tile_count_y = (s->height + s->tile_length - 1) / s->tile_length;
 908
 909     /* Iterate over the number of tiles */
 910     for (tile_idx = 0; tile_idx < s->tile_count; tile_idx++) {
 911         tile_x = tile_idx % tile_count_x;
 912         tile_y = tile_idx / tile_count_x;
 913
 914         if (has_width_leftover && tile_x == tile_count_x - 1) // If on the right-most tile
 915             tile_width = s->width % s->tile_width;
 916         else
 917             tile_width = s->tile_width;
 918
 919         if (has_height_leftover && tile_y == tile_count_y - 1) // If on the bottom-most tile
 920             tile_length = s->height % s->tile_length;
 921         else
 922             tile_length = s->tile_length;
 923
 924         /* Read tile offset */
 925         tile_offset_offset = s->tile_offsets_offset + tile_idx * sizeof(int);
 926         bytestream2_seek(&s->gb, tile_offset_offset, SEEK_SET);
 927         tile_offset = ff_tget_long(&s->gb, s->le);
 928
 929         /* Read tile byte size */
 930         tile_byte_count_offset = s->tile_byte_counts_offset + tile_idx * sizeof(int);
 931         bytestream2_seek(&s->gb, tile_byte_count_offset, SEEK_SET);
 932         tile_byte_count = ff_tget_long(&s->gb, s->le);
 933
 934         /* Seek to tile data */
 935         bytestream2_seek(&s->gb, tile_offset, SEEK_SET);
 936
 937         /* Decode JPEG tile and copy it in the reference frame */
 938         ret = dng_decode_jpeg_tile(avctx, frame, tile_byte_count, pos_x, pos_y, tile_width, tile_length);
 939
 940         if (ret < 0)
 941             return ret;
 942
 943         /* Advance current positions */
 944         pos_x += tile_width;
 945         if (tile_x == tile_count_x - 1) { // If on the right edge
 946             pos_x = 0;
 947             pos_y += tile_length;
 948         }
 949     }
 950
 951     return 0;
 952 }
 953
 954 static int dng_decode(AVCodecContext *avctx, AVFrame *frame, AVPacket *avpkt) {
 955     int ret;
 956
 957     TiffContext *s = avctx->priv_data;
 958
 959     s->jpgframe->width  = s->tile_width;
 960     s->jpgframe->height = s->tile_length;
 961
 962     s->avctx_mjpeg->width = s->tile_width;
 963     s->avctx_mjpeg->height = s->tile_length;
 964
 965     /* Decode all tiles in a frame */
 966     ret = dng_decode_tiles(avctx, frame);
 967     if (ret < 0)
 968         return ret;
 969
 970     /* Frame is ready to be output */
 971     frame->pict_type = AV_PICTURE_TYPE_I;
 972     frame->key_frame = 1;
 973
 974     return avpkt->size;
 975 }
 976
 977 static int init_image(TiffContext *s, ThreadFrame *frame)
 978 {
 979     int ret;
 980     int create_gray_palette = 0;
 981
 982     // make sure there is no aliasing in the following switch
 983     if (s->bpp >= 100 || s->bppcount >= 10) {
 984         av_log(s->avctx, AV_LOG_ERROR,
 985                "Unsupported image parameters: bpp=%d, bppcount=%d\n",
 986                s->bpp, s->bppcount);
 987         return AVERROR_INVALIDDATA;
 988     }
 989
 990     switch (s->planar * 1000 + s->bpp * 10 + s->bppcount + s->is_bayer * 10000) {
 991     case 11:
 992         if (!s->palette_is_set) {
 993             s->avctx->pix_fmt = AV_PIX_FMT_MONOBLACK;
 994             break;
 995         }
 996     case 21:
 997     case 41:
 998         s->avctx->pix_fmt = AV_PIX_FMT_PAL8;
 999         if (!s->palette_is_set) {
1000             create_gray_palette = 1;
1001         }
1002         break;
1003     case 81:
1004         s->avctx->pix_fmt = s->palette_is_set ? AV_PIX_FMT_PAL8 : AV_PIX_FMT_GRAY8;
1005         break;
1006     case 121:
1007         s->avctx->pix_fmt = AV_PIX_FMT_GRAY12;
1008         break;
1009     case 10081:
1010         switch (AV_RL32(s->pattern)) {
1011         case 0x02010100:
1012             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_RGGB8;
1013             break;
1014         case 0x00010102:
1015             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_BGGR8;
1016             break;
1017         case 0x01000201:
1018             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GBRG8;
1019             break;
1020         case 0x01020001:
1021             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GRBG8;
1022             break;
1023         default:
1024             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1025                    AV_RL32(s->pattern));
1026             return AVERROR_PATCHWELCOME;
1027         }
1028         break;
1029     case 10121:
1030         switch (AV_RL32(s->pattern)) {
1031         case 0x02010100:
1032             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_RGGB16LE : AV_PIX_FMT_BAYER_RGGB16BE;
1033             break;
1034         case 0x00010102:
1035             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_BGGR16LE : AV_PIX_FMT_BAYER_BGGR16BE;
1036             break;
1037         case 0x01000201:
1038             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GBRG16LE : AV_PIX_FMT_BAYER_GBRG16BE;
1039             break;
1040         case 0x01020001:
1041             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GRBG16LE : AV_PIX_FMT_BAYER_GRBG16BE;
1042             break;
1043         default:
1044             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1045                    AV_RL32(s->pattern));
1046             return AVERROR_PATCHWELCOME;
1047         }
1048         /* Force endianness as mentioned in 'DNG Specification: Chapter 3: BitsPerSample'
1049            NOTE: The spec actually specifies big-endian, not sure why we need little-endian, but
1050                  such images don't work otherwise. Examples are images produced by Zenmuse X7. */
1051         if ((s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG)
1052             && (s->bpp != 8 && s->bpp != 16 && s->bpp != 32)) {
1053             switch (s->avctx->pix_fmt) {
1054             case AV_PIX_FMT_BAYER_RGGB16BE: s->avctx->pix_fmt = AV_PIX_FMT_BAYER_RGGB16LE; break;
1055             case AV_PIX_FMT_BAYER_BGGR16BE: s->avctx->pix_fmt = AV_PIX_FMT_BAYER_BGGR16LE; break;
1056             case AV_PIX_FMT_BAYER_GBRG16BE: s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GBRG16LE; break;
1057             case AV_PIX_FMT_BAYER_GRBG16BE: s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GRBG16LE; break;
1058             }
1059         }
1060         break;
1061     case 10161:
1062         switch (AV_RL32(s->pattern)) {
1063         case 0x02010100:
1064             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_RGGB16LE : AV_PIX_FMT_BAYER_RGGB16BE;
1065             break;
1066         case 0x00010102:
1067             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_BGGR16LE : AV_PIX_FMT_BAYER_BGGR16BE;
1068             break;
1069         case 0x01000201:
1070             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GBRG16LE : AV_PIX_FMT_BAYER_GBRG16BE;
1071             break;
1072         case 0x01020001:
1073             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GRBG16LE : AV_PIX_FMT_BAYER_GRBG16BE;
1074             break;
1075         default:
1076             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1077                    AV_RL32(s->pattern));
1078             return AVERROR_PATCHWELCOME;
1079         }
1080         break;
1081     case 243:
1082         if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1083             if (s->subsampling[0] == 1 && s->subsampling[1] == 1) {
1084                 s->avctx->pix_fmt = AV_PIX_FMT_YUV444P;
1085             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 1) {
1086                 s->avctx->pix_fmt = AV_PIX_FMT_YUV422P;
1087             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 1) {
1088                 s->avctx->pix_fmt = AV_PIX_FMT_YUV411P;
1089             } else if (s->subsampling[0] == 1 && s->subsampling[1] == 2) {
1090                 s->avctx->pix_fmt = AV_PIX_FMT_YUV440P;
1091             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 2) {
1092                 s->avctx->pix_fmt = AV_PIX_FMT_YUV420P;
1093             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 4) {
1094                 s->avctx->pix_fmt = AV_PIX_FMT_YUV410P;
1095             } else {
1096                 av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr subsampling\n");
1097                 return AVERROR_PATCHWELCOME;
1098             }
1099         } else
1100             s->avctx->pix_fmt = AV_PIX_FMT_RGB24;
1101         break;
1102     case 161:
1103         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GRAY16LE : AV_PIX_FMT_GRAY16BE;
1104         break;
1105     case 162:
1106         s->avctx->pix_fmt = AV_PIX_FMT_YA8;
1107         break;
1108     case 322:
1109         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_YA16LE : AV_PIX_FMT_YA16BE;
1110         break;
1111     case 324:
1112         s->avctx->pix_fmt = s->photometric == TIFF_PHOTOMETRIC_SEPARATED ? AV_PIX_FMT_RGB0 : AV_PIX_FMT_RGBA;
1113         break;
1114     case 405:
1115         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED)
1116             s->avctx->pix_fmt = AV_PIX_FMT_RGBA;
1117         else {
1118             av_log(s->avctx, AV_LOG_ERROR,
1119                 "bpp=40 without PHOTOMETRIC_SEPARATED is unsupported\n");
1120             return AVERROR_PATCHWELCOME;
1121         }
1122         break;
1123     case 483:
1124         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGB48LE  : AV_PIX_FMT_RGB48BE;
1125         break;
1126     case 644:
1127         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGBA64LE  : AV_PIX_FMT_RGBA64BE;
1128         break;
1129     case 1243:
1130         s->avctx->pix_fmt = AV_PIX_FMT_GBRP;
1131         break;
1132     case 1324:
1133         s->avctx->pix_fmt = AV_PIX_FMT_GBRAP;
1134         break;
1135     case 1483:
1136         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRP16LE : AV_PIX_FMT_GBRP16BE;
1137         break;
1138     case 1644:
1139         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRAP16LE : AV_PIX_FMT_GBRAP16BE;
1140         break;
1141     default:
1142         av_log(s->avctx, AV_LOG_ERROR,
1143                "This format is not supported (bpp=%d, bppcount=%d)\n",
1144                s->bpp, s->bppcount);
1145         return AVERROR_INVALIDDATA;
1146     }
1147
1148     if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1149         const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(s->avctx->pix_fmt);
1150         if((desc->flags & AV_PIX_FMT_FLAG_RGB) ||
1151            !(desc->flags & AV_PIX_FMT_FLAG_PLANAR) ||
1152            desc->nb_components < 3) {
1153             av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr variant\n");
1154             return AVERROR_INVALIDDATA;
1155         }
1156     }
1157
1158     if (s->width != s->avctx->width || s->height != s->avctx->height) {
1159         ret = ff_set_dimensions(s->avctx, s->width, s->height);
1160         if (ret < 0)
1161             return ret;
1162     }
1163     if ((ret = ff_thread_get_buffer(s->avctx, frame, 0)) < 0)
1164         return ret;
1165     if (s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
1166         if (!create_gray_palette)
1167             memcpy(frame->f->data[1], s->palette, sizeof(s->palette));
1168         else {
1169             /* make default grayscale pal */
1170             int i;
1171             uint32_t *pal = (uint32_t *)frame->f->data[1];
1172             for (i = 0; i < 1<<s->bpp; i++)
1173                 pal[i] = 0xFFU << 24 | i * 255 / ((1<<s->bpp) - 1) * 0x010101;
1174         }
1175     }
1176     return 0;
1177 }
1178
1179 static void set_sar(TiffContext *s, unsigned tag, unsigned num, unsigned den)
1180 {
1181     int offset = tag == TIFF_YRES ? 2 : 0;
1182     s->res[offset++] = num;
1183     s->res[offset]   = den;
1184     if (s->res[0] && s->res[1] && s->res[2] && s->res[3]) {
1185         uint64_t num = s->res[2] * (uint64_t)s->res[1];
1186         uint64_t den = s->res[0] * (uint64_t)s->res[3];
1187         if (num > INT64_MAX || den > INT64_MAX) {
1188             num = num >> 1;
1189             den = den >> 1;
1190         }
1191         av_reduce(&s->avctx->sample_aspect_ratio.num, &s->avctx->sample_aspect_ratio.den,
1192                   num, den, INT32_MAX);
1193         if (!s->avctx->sample_aspect_ratio.den)
1194             s->avctx->sample_aspect_ratio = (AVRational) {0, 1};
1195     }
1196 }
1197
1198 static int tiff_decode_tag(TiffContext *s, AVFrame *frame)
1199 {
1200     unsigned tag, type, count, off, value = 0, value2 = 1; // value2 is a denominator so init. to 1
1201     int i, start;
1202     int pos;
1203     int ret;
1204     double *dp;
1205
1206     ret = ff_tread_tag(&s->gb, s->le, &tag, &type, &count, &start);
1207     if (ret < 0) {
1208         goto end;
1209     }
1210
1211     off = bytestream2_tell(&s->gb);
1212     if (count == 1) {
1213         switch (type) {
1214         case TIFF_BYTE:
1215         case TIFF_SHORT:
1216         case TIFF_LONG:
1217             value = ff_tget(&s->gb, type, s->le);
1218             break;
1219         case TIFF_RATIONAL:
1220             value  = ff_tget(&s->gb, TIFF_LONG, s->le);
1221             value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
1222             break;
1223         case TIFF_STRING:
1224             if (count <= 4) {
1225                 break;
1226             }
1227         default:
1228             value = UINT_MAX;
1229         }
1230     }
1231
1232     switch (tag) {
1233     case TIFF_SUBFILE:
1234         s->is_thumbnail = (value != 0);
1235         break;
1236     case TIFF_WIDTH:
1237         s->width = value;
1238         break;
1239     case TIFF_HEIGHT:
1240         s->height = value;
1241         break;
1242     case TIFF_BPP:
1243         if (count > 5U) {
1244             av_log(s->avctx, AV_LOG_ERROR,
1245                    "This format is not supported (bpp=%d, %d components)\n",
1246                    value, count);
1247             return AVERROR_INVALIDDATA;
1248         }
1249         s->bppcount = count;
1250         if (count == 1)
1251             s->bpp = value;
1252         else {
1253             switch (type) {
1254             case TIFF_BYTE:
1255             case TIFF_SHORT:
1256             case TIFF_LONG:
1257                 s->bpp = 0;
1258                 if (bytestream2_get_bytes_left(&s->gb) < type_sizes[type] * count)
1259                     return AVERROR_INVALIDDATA;
1260                 for (i = 0; i < count; i++)
1261                     s->bpp += ff_tget(&s->gb, type, s->le);
1262                 break;
1263             default:
1264                 s->bpp = -1;
1265             }
1266         }
1267         break;
1268     case TIFF_SAMPLES_PER_PIXEL:
1269         if (count != 1) {
1270             av_log(s->avctx, AV_LOG_ERROR,
1271                    "Samples per pixel requires a single value, many provided\n");
1272             return AVERROR_INVALIDDATA;
1273         }
1274         if (value > 5U) {
1275             av_log(s->avctx, AV_LOG_ERROR,
1276                    "Samples per pixel %d is too large\n", value);
1277             return AVERROR_INVALIDDATA;
1278         }
1279         if (s->bppcount == 1)
1280             s->bpp *= value;
1281         s->bppcount = value;
1282         break;
1283     case TIFF_COMPR:
1284         s->compr     = value;
1285         av_log(s->avctx, AV_LOG_DEBUG, "compression: %d\n", s->compr);
1286         s->predictor = 0;
1287         switch (s->compr) {
1288         case TIFF_RAW:
1289         case TIFF_PACKBITS:
1290         case TIFF_LZW:
1291         case TIFF_CCITT_RLE:
1292             break;
1293         case TIFF_G3:
1294         case TIFF_G4:
1295             s->fax_opts = 0;
1296             break;
1297         case TIFF_DEFLATE:
1298         case TIFF_ADOBE_DEFLATE:
1299 #if CONFIG_ZLIB
1300             break;
1301 #else
1302             av_log(s->avctx, AV_LOG_ERROR, "Deflate: ZLib not compiled in\n");
1303             return AVERROR(ENOSYS);
1304 #endif
1305         case TIFF_JPEG:
1306         case TIFF_NEWJPEG:
1307             s->is_jpeg = 1;
1308             break;
1309         case TIFF_LZMA:
1310 #if CONFIG_LZMA
1311             break;
1312 #else
1313             av_log(s->avctx, AV_LOG_ERROR, "LZMA not compiled in\n");
1314             return AVERROR(ENOSYS);
1315 #endif
1316         default:
1317             av_log(s->avctx, AV_LOG_ERROR, "Unknown compression method %i\n",
1318                    s->compr);
1319             return AVERROR_INVALIDDATA;
1320         }
1321         break;
1322     case TIFF_ROWSPERSTRIP:
1323         if (!value || (type == TIFF_LONG && value == UINT_MAX))
1324             value = s->height;
1325         s->rps = FFMIN(value, s->height);
1326         break;
1327     case TIFF_STRIP_OFFS:
1328         if (count == 1) {
1329             if (value > INT_MAX) {
1330                 av_log(s->avctx, AV_LOG_ERROR,
1331                     "strippos %u too large\n", value);
1332                 return AVERROR_INVALIDDATA;
1333             }
1334             s->strippos = 0;
1335             s->stripoff = value;
1336         } else
1337             s->strippos = off;
1338         s->strips = count;
1339         if (s->strips == 1)
1340             s->rps = s->height;
1341         s->sot = type;
1342         break;
1343     case TIFF_STRIP_SIZE:
1344         if (count == 1) {
1345             if (value > INT_MAX) {
1346                 av_log(s->avctx, AV_LOG_ERROR,
1347                     "stripsize %u too large\n", value);
1348                 return AVERROR_INVALIDDATA;
1349             }
1350             s->stripsizesoff = 0;
1351             s->stripsize     = value;
1352             s->strips        = 1;
1353         } else {
1354             s->stripsizesoff = off;
1355         }
1356         s->strips = count;
1357         s->sstype = type;
1358         break;
1359     case TIFF_XRES:
1360     case TIFF_YRES:
1361         set_sar(s, tag, value, value2);
1362         break;
1363     case TIFF_TILE_OFFSETS:
1364         s->tile_offsets_offset = off;
1365         s->tile_count = count;
1366         s->is_tiled = 1;
1367         break;
1368     case TIFF_TILE_BYTE_COUNTS:
1369         s->tile_byte_counts_offset = off;
1370         break;
1371     case TIFF_TILE_LENGTH:
1372         s->tile_length = value;
1373         break;
1374     case TIFF_TILE_WIDTH:
1375         s->tile_width = value;
1376         break;
1377     case TIFF_PREDICTOR:
1378         s->predictor = value;
1379         break;
1380     case TIFF_SUB_IFDS:
1381         if (count == 1)
1382             s->sub_ifd = value;
1383         else if (count > 1)
1384             s->sub_ifd = ff_tget(&s->gb, TIFF_LONG, s->le); /** Only get the first SubIFD */
1385         break;
1386     case DNG_LINEARIZATION_TABLE: {
1387         uint32_t lut_offset = value;
1388         uint32_t lut_size = count;
1389         uint32_t lut_wanted_size = 1 << s->bpp;
1390         if (lut_wanted_size != lut_size)
1391             av_log(s->avctx, AV_LOG_WARNING, "DNG contains LUT with invalid size (%"PRIu32"), disabling LUT\n", lut_size);
1392         else if (lut_offset >= bytestream2_size(&s->gb))
1393             av_log(s->avctx, AV_LOG_WARNING, "DNG contains LUT with invalid offset (%"PRIu32"), disabling LUT\n", lut_offset);
1394         else
1395             s->dng_lut = (uint16_t*)(s->gb.buffer + lut_offset);
1396         break;
1397     }
1398     case DNG_BLACK_LEVEL:
1399         if (count > 1) {    /* Use the first value in the pattern (assume they're all the same) */
1400             if (type == TIFF_RATIONAL) {
1401                 value  = ff_tget(&s->gb, TIFF_LONG, s->le);
1402                 value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
1403
1404                 s->black_level = value / value2;
1405             } else
1406                 s->black_level = ff_tget(&s->gb, type, s->le);
1407             av_log(s->avctx, AV_LOG_WARNING, "Assuming black level pattern values are identical\n");
1408         } else {
1409             s->black_level = value / value2;
1410         }
1411         break;
1412     case DNG_WHITE_LEVEL:
1413         s->white_level = value;
1414         break;
1415     case TIFF_CFA_PATTERN_DIM:
1416         if (count != 2 || (ff_tget(&s->gb, type, s->le) != 2 &&
1417                            ff_tget(&s->gb, type, s->le) != 2)) {
1418             av_log(s->avctx, AV_LOG_ERROR, "CFA Pattern dimensions are not 2x2\n");
1419             return AVERROR_INVALIDDATA;
1420         }
1421         break;
1422     case TIFF_CFA_PATTERN:
1423         s->is_bayer = 1;
1424         s->pattern[0] = ff_tget(&s->gb, type, s->le);
1425         s->pattern[1] = ff_tget(&s->gb, type, s->le);
1426         s->pattern[2] = ff_tget(&s->gb, type, s->le);
1427         s->pattern[3] = ff_tget(&s->gb, type, s->le);
1428         break;
1429     case TIFF_PHOTOMETRIC:
1430         switch (value) {
1431         case TIFF_PHOTOMETRIC_WHITE_IS_ZERO:
1432         case TIFF_PHOTOMETRIC_BLACK_IS_ZERO:
1433         case TIFF_PHOTOMETRIC_RGB:
1434         case TIFF_PHOTOMETRIC_PALETTE:
1435         case TIFF_PHOTOMETRIC_SEPARATED:
1436         case TIFF_PHOTOMETRIC_YCBCR:
1437         case TIFF_PHOTOMETRIC_CFA:
1438             s->photometric = value;
1439             break;
1440         case TIFF_PHOTOMETRIC_ALPHA_MASK:
1441         case TIFF_PHOTOMETRIC_CIE_LAB:
1442         case TIFF_PHOTOMETRIC_ICC_LAB:
1443         case TIFF_PHOTOMETRIC_ITU_LAB:
1444         case TIFF_PHOTOMETRIC_LOG_L:
1445         case TIFF_PHOTOMETRIC_LOG_LUV:
1446         case TIFF_PHOTOMETRIC_LINEAR_RAW:
1447             avpriv_report_missing_feature(s->avctx,
1448                                           "PhotometricInterpretation 0x%04X",
1449                                           value);
1450             return AVERROR_PATCHWELCOME;
1451         default:
1452             av_log(s->avctx, AV_LOG_ERROR, "PhotometricInterpretation %u is "
1453                    "unknown\n", value);
1454             return AVERROR_INVALIDDATA;
1455         }
1456         break;
1457     case TIFF_FILL_ORDER:
1458         if (value < 1 || value > 2) {
1459             av_log(s->avctx, AV_LOG_ERROR,
1460                    "Unknown FillOrder value %d, trying default one\n", value);
1461             value = 1;
1462         }
1463         s->fill_order = value - 1;
1464         break;
1465     case TIFF_PAL: {
1466         GetByteContext pal_gb[3];
1467         off = type_sizes[type];
1468         if (count / 3 > 256 ||
1469             bytestream2_get_bytes_left(&s->gb) < count / 3 * off * 3)
1470             return AVERROR_INVALIDDATA;
1471
1472         pal_gb[0] = pal_gb[1] = pal_gb[2] = s->gb;
1473         bytestream2_skip(&pal_gb[1], count / 3 * off);
1474         bytestream2_skip(&pal_gb[2], count / 3 * off * 2);
1475
1476         off = (type_sizes[type] - 1) << 3;
1477         if (off > 31U) {
1478             av_log(s->avctx, AV_LOG_ERROR, "palette shift %d is out of range\n", off);
1479             return AVERROR_INVALIDDATA;
1480         }
1481
1482         for (i = 0; i < count / 3; i++) {
1483             uint32_t p = 0xFF000000;
1484             p |= (ff_tget(&pal_gb[0], type, s->le) >> off) << 16;
1485             p |= (ff_tget(&pal_gb[1], type, s->le) >> off) << 8;
1486             p |=  ff_tget(&pal_gb[2], type, s->le) >> off;
1487             s->palette[i] = p;
1488         }
1489         s->palette_is_set = 1;
1490         break;
1491     }
1492     case TIFF_PLANAR:
1493         s->planar = value == 2;
1494         break;
1495     case TIFF_YCBCR_SUBSAMPLING:
1496         if (count != 2) {
1497             av_log(s->avctx, AV_LOG_ERROR, "subsample count invalid\n");
1498             return AVERROR_INVALIDDATA;
1499         }
1500         for (i = 0; i < count; i++) {
1501             s->subsampling[i] = ff_tget(&s->gb, type, s->le);
1502             if (s->subsampling[i] <= 0) {
1503                 av_log(s->avctx, AV_LOG_ERROR, "subsampling %d is invalid\n", s->subsampling[i]);
1504                 s->subsampling[i] = 1;
1505                 return AVERROR_INVALIDDATA;
1506             }
1507         }
1508         break;
1509     case TIFF_T4OPTIONS:
1510         if (s->compr == TIFF_G3)
1511             s->fax_opts = value;
1512         break;
1513     case TIFF_T6OPTIONS:
1514         if (s->compr == TIFF_G4)
1515             s->fax_opts = value;
1516         break;
1517 #define ADD_METADATA(count, name, sep)\
1518     if ((ret = add_metadata(count, type, name, sep, s, frame)) < 0) {\
1519         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");\
1520         goto end;\
1521     }
1522     case TIFF_MODEL_PIXEL_SCALE:
1523         ADD_METADATA(count, "ModelPixelScaleTag", NULL);
1524         break;
1525     case TIFF_MODEL_TRANSFORMATION:
1526         ADD_METADATA(count, "ModelTransformationTag", NULL);
1527         break;
1528     case TIFF_MODEL_TIEPOINT:
1529         ADD_METADATA(count, "ModelTiepointTag", NULL);
1530         break;
1531     case TIFF_GEO_KEY_DIRECTORY:
1532         if (s->geotag_count) {
1533             avpriv_request_sample(s->avctx, "Multiple geo key directories\n");
1534             return AVERROR_INVALIDDATA;
1535         }
1536         ADD_METADATA(1, "GeoTIFF_Version", NULL);
1537         ADD_METADATA(2, "GeoTIFF_Key_Revision", ".");
1538         s->geotag_count   = ff_tget_short(&s->gb, s->le);
1539         if (s->geotag_count > count / 4 - 1) {
1540             s->geotag_count = count / 4 - 1;
1541             av_log(s->avctx, AV_LOG_WARNING, "GeoTIFF key directory buffer shorter than specified\n");
1542         }
1543         if (   bytestream2_get_bytes_left(&s->gb) < s->geotag_count * sizeof(int16_t) * 4
1544             || s->geotag_count == 0) {
1545             s->geotag_count = 0;
1546             return -1;
1547         }
1548         s->geotags = av_mallocz_array(s->geotag_count, sizeof(TiffGeoTag));
1549         if (!s->geotags) {
1550             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1551             s->geotag_count = 0;
1552             goto end;
1553         }
1554         for (i = 0; i < s->geotag_count; i++) {
1555             s->geotags[i].key    = ff_tget_short(&s->gb, s->le);
1556             s->geotags[i].type   = ff_tget_short(&s->gb, s->le);
1557             s->geotags[i].count  = ff_tget_short(&s->gb, s->le);
1558
1559             if (!s->geotags[i].type)
1560                 s->geotags[i].val  = get_geokey_val(s->geotags[i].key, ff_tget_short(&s->gb, s->le));
1561             else
1562                 s->geotags[i].offset = ff_tget_short(&s->gb, s->le);
1563         }
1564         break;
1565     case TIFF_GEO_DOUBLE_PARAMS:
1566         if (count >= INT_MAX / sizeof(int64_t))
1567             return AVERROR_INVALIDDATA;
1568         if (bytestream2_get_bytes_left(&s->gb) < count * sizeof(int64_t))
1569             return AVERROR_INVALIDDATA;
1570         dp = av_malloc_array(count, sizeof(double));
1571         if (!dp) {
1572             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1573             goto end;
1574         }
1575         for (i = 0; i < count; i++)
1576             dp[i] = ff_tget_double(&s->gb, s->le);
1577         for (i = 0; i < s->geotag_count; i++) {
1578             if (s->geotags[i].type == TIFF_GEO_DOUBLE_PARAMS) {
1579                 if (s->geotags[i].count == 0
1580                     || s->geotags[i].offset + s->geotags[i].count > count) {
1581                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1582                 } else if (s->geotags[i].val) {
1583                     av_log(s->avctx, AV_LOG_WARNING, "Duplicate GeoTIFF key %d\n", s->geotags[i].key);
1584                 } else {
1585                     char *ap = doubles2str(&dp[s->geotags[i].offset], s->geotags[i].count, ", ");
1586                     if (!ap) {
1587                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1588                         av_freep(&dp);
1589                         return AVERROR(ENOMEM);
1590                     }
1591                     s->geotags[i].val = ap;
1592                 }
1593             }
1594         }
1595         av_freep(&dp);
1596         break;
1597     case TIFF_GEO_ASCII_PARAMS:
1598         pos = bytestream2_tell(&s->gb);
1599         for (i = 0; i < s->geotag_count; i++) {
1600             if (s->geotags[i].type == TIFF_GEO_ASCII_PARAMS) {
1601                 if (s->geotags[i].count == 0
1602                     || s->geotags[i].offset +  s->geotags[i].count > count) {
1603                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1604                 } else {
1605                     char *ap;
1606
1607                     bytestream2_seek(&s->gb, pos + s->geotags[i].offset, SEEK_SET);
1608                     if (bytestream2_get_bytes_left(&s->gb) < s->geotags[i].count)
1609                         return AVERROR_INVALIDDATA;
1610                     if (s->geotags[i].val)
1611                         return AVERROR_INVALIDDATA;
1612                     ap = av_malloc(s->geotags[i].count);
1613                     if (!ap) {
1614                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1615                         return AVERROR(ENOMEM);
1616                     }
1617                     bytestream2_get_bufferu(&s->gb, ap, s->geotags[i].count);
1618                     ap[s->geotags[i].count - 1] = '\0'; //replace the "|" delimiter with a 0 byte
1619                     s->geotags[i].val = ap;
1620                 }
1621             }
1622         }
1623         break;
1624     case TIFF_ARTIST:
1625         ADD_METADATA(count, "artist", NULL);
1626         break;
1627     case TIFF_COPYRIGHT:
1628         ADD_METADATA(count, "copyright", NULL);
1629         break;
1630     case TIFF_DATE:
1631         ADD_METADATA(count, "date", NULL);
1632         break;
1633     case TIFF_DOCUMENT_NAME:
1634         ADD_METADATA(count, "document_name", NULL);
1635         break;
1636     case TIFF_HOST_COMPUTER:
1637         ADD_METADATA(count, "computer", NULL);
1638         break;
1639     case TIFF_IMAGE_DESCRIPTION:
1640         ADD_METADATA(count, "description", NULL);
1641         break;
1642     case TIFF_MAKE:
1643         ADD_METADATA(count, "make", NULL);
1644         break;
1645     case TIFF_MODEL:
1646         ADD_METADATA(count, "model", NULL);
1647         break;
1648     case TIFF_PAGE_NAME:
1649         ADD_METADATA(count, "page_name", NULL);
1650         break;
1651     case TIFF_PAGE_NUMBER:
1652         ADD_METADATA(count, "page_number", " / ");
1653         // need to seek back to re-read the page number
1654         bytestream2_seek(&s->gb, -count * sizeof(uint16_t), SEEK_CUR);
1655         // read the page number
1656         s->cur_page = ff_tget(&s->gb, TIFF_SHORT, s->le);
1657         // get back to where we were before the previous seek
1658         bytestream2_seek(&s->gb, count * sizeof(uint16_t) - sizeof(uint16_t), SEEK_CUR);
1659         break;
1660     case TIFF_SOFTWARE_NAME:
1661         ADD_METADATA(count, "software", NULL);
1662         break;
1663     case DNG_VERSION:
1664         if (count == 4) {
1665             unsigned int ver[4];
1666             ver[0] = ff_tget(&s->gb, type, s->le);
1667             ver[1] = ff_tget(&s->gb, type, s->le);
1668             ver[2] = ff_tget(&s->gb, type, s->le);
1669             ver[3] = ff_tget(&s->gb, type, s->le);
1670
1671             av_log(s->avctx, AV_LOG_DEBUG, "DNG file, version %u.%u.%u.%u\n",
1672                 ver[0], ver[1], ver[2], ver[3]);
1673
1674             tiff_set_type(s, TIFF_TYPE_DNG);
1675         }
1676         break;
1677     case CINEMADNG_TIME_CODES:
1678     case CINEMADNG_FRAME_RATE:
1679     case CINEMADNG_T_STOP:
1680     case CINEMADNG_REEL_NAME:
1681     case CINEMADNG_CAMERA_LABEL:
1682         tiff_set_type(s, TIFF_TYPE_CINEMADNG);
1683         break;
1684     default:
1685         if (s->avctx->err_recognition & AV_EF_EXPLODE) {
1686             av_log(s->avctx, AV_LOG_ERROR,
1687                    "Unknown or unsupported tag %d/0x%0X\n",
1688                    tag, tag);
1689             return AVERROR_INVALIDDATA;
1690         }
1691     }
1692 end:
1693     if (s->bpp > 64U) {
1694         av_log(s->avctx, AV_LOG_ERROR,
1695                 "This format is not supported (bpp=%d, %d components)\n",
1696                 s->bpp, count);
1697         s->bpp = 0;
1698         return AVERROR_INVALIDDATA;
1699     }
1700     bytestream2_seek(&s->gb, start, SEEK_SET);
1701     return 0;
1702 }
1703
1704 static int decode_frame(AVCodecContext *avctx,
1705                         void *data, int *got_frame, AVPacket *avpkt)
1706 {
1707     TiffContext *const s = avctx->priv_data;
1708     AVFrame *const p = data;
1709     ThreadFrame frame = { .f = data };
1710     unsigned off, last_off;
1711     int le, ret, plane, planes;
1712     int i, j, entries, stride;
1713     unsigned soff, ssize;
1714     uint8_t *dst;
1715     GetByteContext stripsizes;
1716     GetByteContext stripdata;
1717     int retry_for_subifd, retry_for_page;
1718
1719     bytestream2_init(&s->gb, avpkt->data, avpkt->size);
1720
1721     // parse image header
1722     if ((ret = ff_tdecode_header(&s->gb, &le, &off))) {
1723         av_log(avctx, AV_LOG_ERROR, "Invalid TIFF header\n");
1724         return ret;
1725     } else if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1726         av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1727         return AVERROR_INVALIDDATA;
1728     }
1729     s->le          = le;
1730     // TIFF_BPP is not a required tag and defaults to 1
1731
1732     s->tiff_type   = TIFF_TYPE_TIFF;
1733 again:
1734     s->is_thumbnail = 0;
1735     s->bppcount    = s->bpp = 1;
1736     s->photometric = TIFF_PHOTOMETRIC_NONE;
1737     s->compr       = TIFF_RAW;
1738     s->fill_order  = 0;
1739     s->white_level = 0;
1740     s->is_bayer    = 0;
1741     s->is_tiled    = 0;
1742     s->is_jpeg     = 0;
1743     s->cur_page    = 0;
1744     s->dng_lut     = NULL;
1745     free_geotags(s);
1746
1747     // Reset these offsets so we can tell if they were set this frame
1748     s->stripsizesoff = s->strippos = 0;
1749     /* parse image file directory */
1750     bytestream2_seek(&s->gb, off, SEEK_SET);
1751     entries = ff_tget_short(&s->gb, le);
1752     if (bytestream2_get_bytes_left(&s->gb) < entries * 12)
1753         return AVERROR_INVALIDDATA;
1754     for (i = 0; i < entries; i++) {
1755         if ((ret = tiff_decode_tag(s, p)) < 0)
1756             return ret;
1757     }
1758
1759     if (s->get_thumbnail && !s->is_thumbnail) {
1760         av_log(avctx, AV_LOG_INFO, "No embedded thumbnail present\n");
1761         return AVERROR_EOF;
1762     }
1763
1764     /** whether we should process this IFD's SubIFD */
1765     retry_for_subifd = s->sub_ifd && (s->get_subimage || (!s->get_thumbnail && s->is_thumbnail));
1766     /** whether we should process this multi-page IFD's next page */
1767     retry_for_page = s->get_page && s->cur_page + 1 < s->get_page;  // get_page is 1-indexed
1768
1769     last_off = off;
1770     if (retry_for_page) {
1771         // set offset to the next IFD
1772         off = ff_tget_long(&s->gb, le);
1773     } else if (retry_for_subifd) {
1774         // set offset to the SubIFD
1775         off = s->sub_ifd;
1776     }
1777
1778     if (retry_for_subifd || retry_for_page) {
1779         if (!off) {
1780             av_log(avctx, AV_LOG_ERROR, "Requested entry not found\n");
1781             return AVERROR_INVALIDDATA;
1782         }
1783         if (off <= last_off) {
1784             avpriv_request_sample(s->avctx, "non increasing IFD offset\n");
1785             return AVERROR_INVALIDDATA;
1786         }
1787         if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1788             av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1789             return AVERROR_INVALIDDATA;
1790         }
1791         s->sub_ifd = 0;
1792         goto again;
1793     }
1794
1795     for (i = 0; i<s->geotag_count; i++) {
1796         const char *keyname = get_geokey_name(s->geotags[i].key);
1797         if (!keyname) {
1798             av_log(avctx, AV_LOG_WARNING, "Unknown or unsupported GeoTIFF key %d\n", s->geotags[i].key);
1799             continue;
1800         }
1801         if (get_geokey_type(s->geotags[i].key) != s->geotags[i].type) {
1802             av_log(avctx, AV_LOG_WARNING, "Type of GeoTIFF key %d is wrong\n", s->geotags[i].key);
1803             continue;
1804         }
1805         ret = av_dict_set(&p->metadata, keyname, s->geotags[i].val, 0);
1806         if (ret<0) {
1807             av_log(avctx, AV_LOG_ERROR, "Writing metadata with key '%s' failed\n", keyname);
1808             return ret;
1809         }
1810     }
1811
1812     if (!s->is_tiled && !s->strippos && !s->stripoff) {
1813         av_log(avctx, AV_LOG_ERROR, "Image data is missing\n");
1814         return AVERROR_INVALIDDATA;
1815     }
1816     /* now we have the data and may start decoding */
1817     if ((ret = init_image(s, &frame)) < 0)
1818         return ret;
1819
1820     if (!s->is_tiled) {
1821         if (s->strips == 1 && !s->stripsize) {
1822             av_log(avctx, AV_LOG_WARNING, "Image data size missing\n");
1823             s->stripsize = avpkt->size - s->stripoff;
1824         }
1825
1826         if (s->stripsizesoff) {
1827             if (s->stripsizesoff >= (unsigned)avpkt->size)
1828                 return AVERROR_INVALIDDATA;
1829             bytestream2_init(&stripsizes, avpkt->data + s->stripsizesoff,
1830                             avpkt->size - s->stripsizesoff);
1831         }
1832         if (s->strippos) {
1833             if (s->strippos >= (unsigned)avpkt->size)
1834                 return AVERROR_INVALIDDATA;
1835             bytestream2_init(&stripdata, avpkt->data + s->strippos,
1836                             avpkt->size - s->strippos);
1837         }
1838
1839         if (s->rps <= 0 || s->rps % s->subsampling[1]) {
1840             av_log(avctx, AV_LOG_ERROR, "rps %d invalid\n", s->rps);
1841             return AVERROR_INVALIDDATA;
1842         }
1843     }
1844
1845     /* Handle DNG images with JPEG-compressed tiles */
1846
1847     if ((s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG) && s->is_tiled) {
1848         if (!s->is_jpeg) {
1849             avpriv_report_missing_feature(avctx, "DNG uncompressed tiled images");
1850             return AVERROR_PATCHWELCOME;
1851         } else if (!s->is_bayer) {
1852             avpriv_report_missing_feature(avctx, "DNG JPG-compressed tiled non-bayer-encoded images");
1853             return AVERROR_PATCHWELCOME;
1854         } else {
1855             if ((ret = dng_decode(avctx, (AVFrame*)data, avpkt)) > 0)
1856                 *got_frame = 1;
1857             return ret;
1858         }
1859     }
1860
1861     /* Handle TIFF images and DNG images with uncompressed strips (non-tiled) */
1862
1863     planes = s->planar ? s->bppcount : 1;
1864     for (plane = 0; plane < planes; plane++) {
1865         uint8_t *five_planes = NULL;
1866         int remaining = avpkt->size;
1867         int decoded_height;
1868         stride = p->linesize[plane];
1869         dst = p->data[plane];
1870         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1871             s->avctx->pix_fmt == AV_PIX_FMT_RGBA) {
1872             stride = stride * 5 / 4;
1873             five_planes =
1874             dst = av_malloc(stride * s->height);
1875             if (!dst)
1876                 return AVERROR(ENOMEM);
1877         }
1878         for (i = 0; i < s->height; i += s->rps) {
1879             if (i)
1880                 dst += s->rps * stride;
1881             if (s->stripsizesoff)
1882                 ssize = ff_tget(&stripsizes, s->sstype, le);
1883             else
1884                 ssize = s->stripsize;
1885
1886             if (s->strippos)
1887                 soff = ff_tget(&stripdata, s->sot, le);
1888             else
1889                 soff = s->stripoff;
1890
1891             if (soff > avpkt->size || ssize > avpkt->size - soff || ssize > remaining) {
1892                 av_log(avctx, AV_LOG_ERROR, "Invalid strip size/offset\n");
1893                 av_freep(&five_planes);
1894                 return AVERROR_INVALIDDATA;
1895             }
1896             remaining -= ssize;
1897             if ((ret = tiff_unpack_strip(s, p, dst, stride, avpkt->data + soff, ssize, i,
1898                                          FFMIN(s->rps, s->height - i))) < 0) {
1899                 if (avctx->err_recognition & AV_EF_EXPLODE) {
1900                     av_freep(&five_planes);
1901                     return ret;
1902                 }
1903                 break;
1904             }
1905         }
1906         decoded_height = FFMIN(i, s->height);
1907
1908         if (s->predictor == 2) {
1909             if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1910                 av_log(s->avctx, AV_LOG_ERROR, "predictor == 2 with YUV is unsupported");
1911                 return AVERROR_PATCHWELCOME;
1912             }
1913             dst   = five_planes ? five_planes : p->data[plane];
1914             soff  = s->bpp >> 3;
1915             if (s->planar)
1916                 soff  = FFMAX(soff / s->bppcount, 1);
1917             ssize = s->width * soff;
1918             if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48LE ||
1919                 s->avctx->pix_fmt == AV_PIX_FMT_RGBA64LE ||
1920                 s->avctx->pix_fmt == AV_PIX_FMT_GRAY16LE ||
1921                 s->avctx->pix_fmt == AV_PIX_FMT_YA16LE ||
1922                 s->avctx->pix_fmt == AV_PIX_FMT_GBRP16LE ||
1923                 s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16LE) {
1924                 for (i = 0; i < decoded_height; i++) {
1925                     for (j = soff; j < ssize; j += 2)
1926                         AV_WL16(dst + j, AV_RL16(dst + j) + AV_RL16(dst + j - soff));
1927                     dst += stride;
1928                 }
1929             } else if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48BE ||
1930                        s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE ||
1931                        s->avctx->pix_fmt == AV_PIX_FMT_GRAY16BE ||
1932                        s->avctx->pix_fmt == AV_PIX_FMT_YA16BE ||
1933                        s->avctx->pix_fmt == AV_PIX_FMT_GBRP16BE ||
1934                        s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16BE) {
1935                 for (i = 0; i < decoded_height; i++) {
1936                     for (j = soff; j < ssize; j += 2)
1937                         AV_WB16(dst + j, AV_RB16(dst + j) + AV_RB16(dst + j - soff));
1938                     dst += stride;
1939                 }
1940             } else {
1941                 for (i = 0; i < decoded_height; i++) {
1942                     for (j = soff; j < ssize; j++)
1943                         dst[j] += dst[j - soff];
1944                     dst += stride;
1945                 }
1946             }
1947         }
1948
1949         if (s->photometric == TIFF_PHOTOMETRIC_WHITE_IS_ZERO) {
1950             int c = (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 ? (1<<s->bpp) - 1 : 255);
1951             dst = p->data[plane];
1952             for (i = 0; i < s->height; i++) {
1953                 for (j = 0; j < stride; j++)
1954                     dst[j] = c - dst[j];
1955                 dst += stride;
1956             }
1957         }
1958
1959         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1960             (s->avctx->pix_fmt == AV_PIX_FMT_RGB0 || s->avctx->pix_fmt == AV_PIX_FMT_RGBA)) {
1961             int x = s->avctx->pix_fmt == AV_PIX_FMT_RGB0 ? 4 : 5;
1962             uint8_t *src = five_planes ? five_planes : p->data[plane];
1963             dst = p->data[plane];
1964             for (i = 0; i < s->height; i++) {
1965                 for (j = 0; j < s->width; j++) {
1966                     int k =  255 - src[x * j + 3];
1967                     int r = (255 - src[x * j    ]) * k;
1968                     int g = (255 - src[x * j + 1]) * k;
1969                     int b = (255 - src[x * j + 2]) * k;
1970                     dst[4 * j    ] = r * 257 >> 16;
1971                     dst[4 * j + 1] = g * 257 >> 16;
1972                     dst[4 * j + 2] = b * 257 >> 16;
1973                     dst[4 * j + 3] = s->avctx->pix_fmt == AV_PIX_FMT_RGBA ? src[x * j + 4] : 255;
1974                 }
1975                 src += stride;
1976                 dst += p->linesize[plane];
1977             }
1978             av_freep(&five_planes);
1979         } else if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1980             s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE) {
1981             dst = p->data[plane];
1982             for (i = 0; i < s->height; i++) {
1983                 for (j = 0; j < s->width; j++) {
1984                     uint64_t k =  65535 - AV_RB16(dst + 8 * j + 6);
1985                     uint64_t r = (65535 - AV_RB16(dst + 8 * j    )) * k;
1986                     uint64_t g = (65535 - AV_RB16(dst + 8 * j + 2)) * k;
1987                     uint64_t b = (65535 - AV_RB16(dst + 8 * j + 4)) * k;
1988                     AV_WB16(dst + 8 * j    , r * 65537 >> 32);
1989                     AV_WB16(dst + 8 * j + 2, g * 65537 >> 32);
1990                     AV_WB16(dst + 8 * j + 4, b * 65537 >> 32);
1991                     AV_WB16(dst + 8 * j + 6, 65535);
1992                 }
1993                 dst += p->linesize[plane];
1994             }
1995         }
1996     }
1997
1998     if (s->planar && s->bppcount > 2) {
1999         FFSWAP(uint8_t*, p->data[0],     p->data[2]);
2000         FFSWAP(int,      p->linesize[0], p->linesize[2]);
2001         FFSWAP(uint8_t*, p->data[0],     p->data[1]);
2002         FFSWAP(int,      p->linesize[0], p->linesize[1]);
2003     }
2004
2005     if (s->is_bayer && s->white_level && s->bpp == 16 &&
2006         !(s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG)) {
2007         uint16_t *dst = (uint16_t *)p->data[0];
2008         for (i = 0; i < s->height; i++) {
2009             for (j = 0; j < s->width; j++)
2010                 dst[j] = FFMIN((dst[j] / (float)s->white_level) * 65535, 65535);
2011             dst += stride / 2;
2012         }
2013     }
2014
2015     *got_frame = 1;
2016
2017     return avpkt->size;
2018 }
2019
2020 static av_cold int tiff_init(AVCodecContext *avctx)
2021 {
2022     TiffContext *s = avctx->priv_data;
2023     const AVCodec *codec;
2024     int ret;
2025
2026     s->width  = 0;
2027     s->height = 0;
2028     s->subsampling[0] =
2029     s->subsampling[1] = 1;
2030     s->avctx  = avctx;
2031     ff_lzw_decode_open(&s->lzw);
2032     if (!s->lzw)
2033         return AVERROR(ENOMEM);
2034     ff_ccitt_unpack_init();
2035
2036     /* Allocate JPEG frame */
2037     s->jpgframe = av_frame_alloc();
2038     if (!s->jpgframe)
2039         return AVERROR(ENOMEM);
2040
2041     /* Prepare everything needed for JPEG decoding */
2042     codec = avcodec_find_decoder(AV_CODEC_ID_MJPEG);
2043     if (!codec)
2044         return AVERROR_BUG;
2045     s->avctx_mjpeg = avcodec_alloc_context3(codec);
2046     if (!s->avctx_mjpeg)
2047         return AVERROR(ENOMEM);
2048     s->avctx_mjpeg->flags = avctx->flags;
2049     s->avctx_mjpeg->flags2 = avctx->flags2;
2050     s->avctx_mjpeg->dct_algo = avctx->dct_algo;
2051     s->avctx_mjpeg->idct_algo = avctx->idct_algo;
2052     ret = ff_codec_open2_recursive(s->avctx_mjpeg, codec, NULL);
2053     if (ret < 0) {
2054         av_frame_free(&s->jpgframe);
2055         avcodec_free_context(&s->avctx_mjpeg);
2056         return ret;
2057     }
2058
2059     return 0;
2060 }
2061
2062 static av_cold int tiff_end(AVCodecContext *avctx)
2063 {
2064     TiffContext *const s = avctx->priv_data;
2065
2066     free_geotags(s);
2067
2068     ff_lzw_decode_close(&s->lzw);
2069     av_freep(&s->deinvert_buf);
2070     s->deinvert_buf_size = 0;
2071     av_freep(&s->yuv_line);
2072     s->yuv_line_size = 0;
2073     av_freep(&s->fax_buffer);
2074     s->fax_buffer_size = 0;
2075     av_frame_free(&s->jpgframe);
2076     avcodec_free_context(&s->avctx_mjpeg);
2077     return 0;
2078 }
2079
2080 #define OFFSET(x) offsetof(TiffContext, x)
2081 static const AVOption tiff_options[] = {
2082     { "subimage", "decode subimage instead if available", OFFSET(get_subimage), AV_OPT_TYPE_BOOL, {.i64=0},  0, 1, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2083     { "thumbnail", "decode embedded thumbnail subimage instead if available", OFFSET(get_thumbnail), AV_OPT_TYPE_BOOL, {.i64=0},  0, 1, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2084     { "page", "page number of multi-page image to decode (starting from 1)", OFFSET(get_page), AV_OPT_TYPE_INT, {.i64=0}, 0, UINT16_MAX, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2085     { NULL },
2086 };
2087
2088 static const AVClass tiff_decoder_class = {
2089     .class_name = "TIFF decoder",
2090     .item_name  = av_default_item_name,
2091     .option     = tiff_options,
2092     .version    = LIBAVUTIL_VERSION_INT,
2093 };
2094
2095 AVCodec ff_tiff_decoder = {
2096     .name           = "tiff",
2097     .long_name      = NULL_IF_CONFIG_SMALL("TIFF image"),
2098     .type           = AVMEDIA_TYPE_VIDEO,
2099     .id             = AV_CODEC_ID_TIFF,
2100     .priv_data_size = sizeof(TiffContext),
2101     .init           = tiff_init,
2102     .close          = tiff_end,
2103     .decode         = decode_frame,
2104     .init_thread_copy = ONLY_IF_THREADS_ENABLED(tiff_init),
2105     .capabilities   = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_FRAME_THREADS,
2106     .priv_class     = &tiff_decoder_class,
2107 };