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