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/intreadwrite.h"
  39 #include "libavutil/imgutils.h"
  40 #include "libavutil/opt.h"
  41 #include "avcodec.h"
  42 #include "bytestream.h"
  43 #include "faxcompr.h"
  44 #include "internal.h"
  45 #include "lzw.h"
  46 #include "mathops.h"
  47 #include "tiff.h"
  48 #include "tiff_data.h"
  49 #include "thread.h"
  50 #include "get_bits.h"
  51
  52 typedef struct TiffContext {
  53     AVClass *class;
  54     AVCodecContext *avctx;
  55     GetByteContext gb;
  56
  57     int get_subimage;
  58     uint16_t get_page;
  59     int get_thumbnail;
  60
  61     enum TiffType tiff_type;
  62     int width, height;
  63     unsigned int bpp, bppcount;
  64     uint32_t palette[256];
  65     int palette_is_set;
  66     int le;
  67     enum TiffCompr compr;
  68     enum TiffPhotometric photometric;
  69     int planar;
  70     int subsampling[2];
  71     int fax_opts;
  72     int predictor;
  73     int fill_order;
  74     uint32_t res[4];
  75     int is_thumbnail;
  76
  77     int is_bayer;
  78     uint8_t pattern[4];
  79     unsigned white_level;
  80
  81     uint32_t sub_ifd;
  82     uint16_t cur_page;
  83
  84     int strips, rps, sstype;
  85     int sot;
  86     int stripsizesoff, stripsize, stripoff, strippos;
  87     LZWState *lzw;
  88
  89     uint8_t *deinvert_buf;
  90     int deinvert_buf_size;
  91     uint8_t *yuv_line;
  92     unsigned int yuv_line_size;
  93     uint8_t *fax_buffer;
  94     unsigned int fax_buffer_size;
  95
  96     int geotag_count;
  97     TiffGeoTag *geotags;
  98 } TiffContext;
  99
 100 static void tiff_set_type(TiffContext *s, enum TiffType tiff_type) {
 101     if (s->tiff_type < tiff_type) // Prioritize higher-valued entries
 102         s->tiff_type = tiff_type;
 103 }
 104
 105 static void free_geotags(TiffContext *const s)
 106 {
 107     int i;
 108     for (i = 0; i < s->geotag_count; i++) {
 109         if (s->geotags[i].val)
 110             av_freep(&s->geotags[i].val);
 111     }
 112     av_freep(&s->geotags);
 113     s->geotag_count = 0;
 114 }
 115
 116 #define RET_GEOKEY(TYPE, array, element)\
 117     if (key >= TIFF_##TYPE##_KEY_ID_OFFSET &&\
 118         key - TIFF_##TYPE##_KEY_ID_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_name_type_map))\
 119         return ff_tiff_##array##_name_type_map[key - TIFF_##TYPE##_KEY_ID_OFFSET].element;
 120
 121 static const char *get_geokey_name(int key)
 122 {
 123     RET_GEOKEY(VERT, vert, name);
 124     RET_GEOKEY(PROJ, proj, name);
 125     RET_GEOKEY(GEOG, geog, name);
 126     RET_GEOKEY(CONF, conf, name);
 127
 128     return NULL;
 129 }
 130
 131 static int get_geokey_type(int key)
 132 {
 133     RET_GEOKEY(VERT, vert, type);
 134     RET_GEOKEY(PROJ, proj, type);
 135     RET_GEOKEY(GEOG, geog, type);
 136     RET_GEOKEY(CONF, conf, type);
 137
 138     return AVERROR_INVALIDDATA;
 139 }
 140
 141 static int cmp_id_key(const void *id, const void *k)
 142 {
 143     return *(const int*)id - ((const TiffGeoTagKeyName*)k)->key;
 144 }
 145
 146 static const char *search_keyval(const TiffGeoTagKeyName *keys, int n, int id)
 147 {
 148     TiffGeoTagKeyName *r = bsearch(&id, keys, n, sizeof(keys[0]), cmp_id_key);
 149     if(r)
 150         return r->name;
 151
 152     return NULL;
 153 }
 154
 155 static char *get_geokey_val(int key, int val)
 156 {
 157     char *ap;
 158
 159     if (val == TIFF_GEO_KEY_UNDEFINED)
 160         return av_strdup("undefined");
 161     if (val == TIFF_GEO_KEY_USER_DEFINED)
 162         return av_strdup("User-Defined");
 163
 164 #define RET_GEOKEY_VAL(TYPE, array)\
 165     if (val >= TIFF_##TYPE##_OFFSET &&\
 166         val - TIFF_##TYPE##_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_codes))\
 167         return av_strdup(ff_tiff_##array##_codes[val - TIFF_##TYPE##_OFFSET]);
 168
 169     switch (key) {
 170     case TIFF_GT_MODEL_TYPE_GEOKEY:
 171         RET_GEOKEY_VAL(GT_MODEL_TYPE, gt_model_type);
 172         break;
 173     case TIFF_GT_RASTER_TYPE_GEOKEY:
 174         RET_GEOKEY_VAL(GT_RASTER_TYPE, gt_raster_type);
 175         break;
 176     case TIFF_GEOG_LINEAR_UNITS_GEOKEY:
 177     case TIFF_PROJ_LINEAR_UNITS_GEOKEY:
 178     case TIFF_VERTICAL_UNITS_GEOKEY:
 179         RET_GEOKEY_VAL(LINEAR_UNIT, linear_unit);
 180         break;
 181     case TIFF_GEOG_ANGULAR_UNITS_GEOKEY:
 182     case TIFF_GEOG_AZIMUTH_UNITS_GEOKEY:
 183         RET_GEOKEY_VAL(ANGULAR_UNIT, angular_unit);
 184         break;
 185     case TIFF_GEOGRAPHIC_TYPE_GEOKEY:
 186         RET_GEOKEY_VAL(GCS_TYPE, gcs_type);
 187         RET_GEOKEY_VAL(GCSE_TYPE, gcse_type);
 188         break;
 189     case TIFF_GEOG_GEODETIC_DATUM_GEOKEY:
 190         RET_GEOKEY_VAL(GEODETIC_DATUM, geodetic_datum);
 191         RET_GEOKEY_VAL(GEODETIC_DATUM_E, geodetic_datum_e);
 192         break;
 193     case TIFF_GEOG_ELLIPSOID_GEOKEY:
 194         RET_GEOKEY_VAL(ELLIPSOID, ellipsoid);
 195         break;
 196     case TIFF_GEOG_PRIME_MERIDIAN_GEOKEY:
 197         RET_GEOKEY_VAL(PRIME_MERIDIAN, prime_meridian);
 198         break;
 199     case TIFF_PROJECTED_CS_TYPE_GEOKEY:
 200         ap = av_strdup(search_keyval(ff_tiff_proj_cs_type_codes, FF_ARRAY_ELEMS(ff_tiff_proj_cs_type_codes), val));
 201         if(ap) return ap;
 202         break;
 203     case TIFF_PROJECTION_GEOKEY:
 204         ap = av_strdup(search_keyval(ff_tiff_projection_codes, FF_ARRAY_ELEMS(ff_tiff_projection_codes), val));
 205         if(ap) return ap;
 206         break;
 207     case TIFF_PROJ_COORD_TRANS_GEOKEY:
 208         RET_GEOKEY_VAL(COORD_TRANS, coord_trans);
 209         break;
 210     case TIFF_VERTICAL_CS_TYPE_GEOKEY:
 211         RET_GEOKEY_VAL(VERT_CS, vert_cs);
 212         RET_GEOKEY_VAL(ORTHO_VERT_CS, ortho_vert_cs);
 213         break;
 214
 215     }
 216
 217     ap = av_malloc(14);
 218     if (ap)
 219         snprintf(ap, 14, "Unknown-%d", val);
 220     return ap;
 221 }
 222
 223 static char *doubles2str(double *dp, int count, const char *sep)
 224 {
 225     int i;
 226     char *ap, *ap0;
 227     uint64_t component_len;
 228     if (!sep) sep = ", ";
 229     component_len = 24LL + strlen(sep);
 230     if (count >= (INT_MAX - 1)/component_len)
 231         return NULL;
 232     ap = av_malloc(component_len * count + 1);
 233     if (!ap)
 234         return NULL;
 235     ap0   = ap;
 236     ap[0] = '\0';
 237     for (i = 0; i < count; i++) {
 238         unsigned l = snprintf(ap, component_len, "%.15g%s", dp[i], sep);
 239         if(l >= component_len) {
 240             av_free(ap0);
 241             return NULL;
 242         }
 243         ap += l;
 244     }
 245     ap0[strlen(ap0) - strlen(sep)] = '\0';
 246     return ap0;
 247 }
 248
 249 static int add_metadata(int count, int type,
 250                         const char *name, const char *sep, TiffContext *s, AVFrame *frame)
 251 {
 252     switch(type) {
 253     case TIFF_DOUBLE: return ff_tadd_doubles_metadata(count, name, sep, &s->gb, s->le, &frame->metadata);
 254     case TIFF_SHORT : return ff_tadd_shorts_metadata(count, name, sep, &s->gb, s->le, 0, &frame->metadata);
 255     case TIFF_STRING: return ff_tadd_string_metadata(count, name, &s->gb, s->le, &frame->metadata);
 256     default         : return AVERROR_INVALIDDATA;
 257     };
 258 }
 259
 260 static void av_always_inline horizontal_fill(TiffContext *s,
 261                                              unsigned int bpp, uint8_t* dst,
 262                                              int usePtr, const uint8_t *src,
 263                                              uint8_t c, int width, int offset)
 264 {
 265     switch (bpp) {
 266     case 1:
 267         while (--width >= 0) {
 268             dst[(width+offset)*8+7] = (usePtr ? src[width] : c)      & 0x1;
 269             dst[(width+offset)*8+6] = (usePtr ? src[width] : c) >> 1 & 0x1;
 270             dst[(width+offset)*8+5] = (usePtr ? src[width] : c) >> 2 & 0x1;
 271             dst[(width+offset)*8+4] = (usePtr ? src[width] : c) >> 3 & 0x1;
 272             dst[(width+offset)*8+3] = (usePtr ? src[width] : c) >> 4 & 0x1;
 273             dst[(width+offset)*8+2] = (usePtr ? src[width] : c) >> 5 & 0x1;
 274             dst[(width+offset)*8+1] = (usePtr ? src[width] : c) >> 6 & 0x1;
 275             dst[(width+offset)*8+0] = (usePtr ? src[width] : c) >> 7;
 276         }
 277         break;
 278     case 2:
 279         while (--width >= 0) {
 280             dst[(width+offset)*4+3] = (usePtr ? src[width] : c) & 0x3;
 281             dst[(width+offset)*4+2] = (usePtr ? src[width] : c) >> 2 & 0x3;
 282             dst[(width+offset)*4+1] = (usePtr ? src[width] : c) >> 4 & 0x3;
 283             dst[(width+offset)*4+0] = (usePtr ? src[width] : c) >> 6;
 284         }
 285         break;
 286     case 4:
 287         while (--width >= 0) {
 288             dst[(width+offset)*2+1] = (usePtr ? src[width] : c) & 0xF;
 289             dst[(width+offset)*2+0] = (usePtr ? src[width] : c) >> 4;
 290         }
 291         break;
 292     case 12: {
 293                  uint16_t *dst16 = (uint16_t *)dst;
 294                  GetBitContext gb;
 295                  init_get_bits8(&gb, src, width);
 296                  for (int i = 0; i < s->width; i++) {
 297                      dst16[i] = get_bits(&gb, 12) << 4;
 298                  }
 299              }
 300         break;
 301     default:
 302         if (usePtr) {
 303             memcpy(dst + offset, src, width);
 304         } else {
 305             memset(dst + offset, c, width);
 306         }
 307     }
 308 }
 309
 310 static int deinvert_buffer(TiffContext *s, const uint8_t *src, int size)
 311 {
 312     int i;
 313
 314     av_fast_padded_malloc(&s->deinvert_buf, &s->deinvert_buf_size, size);
 315     if (!s->deinvert_buf)
 316         return AVERROR(ENOMEM);
 317     for (i = 0; i < size; i++)
 318         s->deinvert_buf[i] = ff_reverse[src[i]];
 319
 320     return 0;
 321 }
 322
 323 static void unpack_gray(TiffContext *s, AVFrame *p,
 324                        const uint8_t *src, int lnum, int width, int bpp)
 325 {
 326     GetBitContext gb;
 327     uint16_t *dst = (uint16_t *)(p->data[0] + lnum * p->linesize[0]);
 328
 329     init_get_bits8(&gb, src, width);
 330
 331     for (int i = 0; i < s->width; i++) {
 332         dst[i] = get_bits(&gb, bpp);
 333     }
 334 }
 335
 336 static void unpack_yuv(TiffContext *s, AVFrame *p,
 337                        const uint8_t *src, int lnum)
 338 {
 339     int i, j, k;
 340     int w       = (s->width - 1) / s->subsampling[0] + 1;
 341     uint8_t *pu = &p->data[1][lnum / s->subsampling[1] * p->linesize[1]];
 342     uint8_t *pv = &p->data[2][lnum / s->subsampling[1] * p->linesize[2]];
 343     if (s->width % s->subsampling[0] || s->height % s->subsampling[1]) {
 344         for (i = 0; i < w; i++) {
 345             for (j = 0; j < s->subsampling[1]; j++)
 346                 for (k = 0; k < s->subsampling[0]; k++)
 347                     p->data[0][FFMIN(lnum + j, s->height-1) * p->linesize[0] +
 348                                FFMIN(i * s->subsampling[0] + k, s->width-1)] = *src++;
 349             *pu++ = *src++;
 350             *pv++ = *src++;
 351         }
 352     }else{
 353         for (i = 0; i < w; i++) {
 354             for (j = 0; j < s->subsampling[1]; j++)
 355                 for (k = 0; k < s->subsampling[0]; k++)
 356                     p->data[0][(lnum + j) * p->linesize[0] +
 357                                i * s->subsampling[0] + k] = *src++;
 358             *pu++ = *src++;
 359             *pv++ = *src++;
 360         }
 361     }
 362 }
 363
 364 #if CONFIG_ZLIB
 365 static int tiff_uncompress(uint8_t *dst, unsigned long *len, const uint8_t *src,
 366                            int size)
 367 {
 368     z_stream zstream = { 0 };
 369     int zret;
 370
 371     zstream.next_in   = (uint8_t *)src;
 372     zstream.avail_in  = size;
 373     zstream.next_out  = dst;
 374     zstream.avail_out = *len;
 375     zret              = inflateInit(&zstream);
 376     if (zret != Z_OK) {
 377         av_log(NULL, AV_LOG_ERROR, "Inflate init error: %d\n", zret);
 378         return zret;
 379     }
 380     zret = inflate(&zstream, Z_SYNC_FLUSH);
 381     inflateEnd(&zstream);
 382     *len = zstream.total_out;
 383     return zret == Z_STREAM_END ? Z_OK : zret;
 384 }
 385
 386 static int tiff_unpack_zlib(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
 387                             const uint8_t *src, int size, int width, int lines,
 388                             int strip_start, int is_yuv)
 389 {
 390     uint8_t *zbuf;
 391     unsigned long outlen;
 392     int ret, line;
 393     outlen = width * lines;
 394     zbuf   = av_malloc(outlen);
 395     if (!zbuf)
 396         return AVERROR(ENOMEM);
 397     if (s->fill_order) {
 398         if ((ret = deinvert_buffer(s, src, size)) < 0) {
 399             av_free(zbuf);
 400             return ret;
 401         }
 402         src = s->deinvert_buf;
 403     }
 404     ret = tiff_uncompress(zbuf, &outlen, src, size);
 405     if (ret != Z_OK) {
 406         av_log(s->avctx, AV_LOG_ERROR,
 407                "Uncompressing failed (%lu of %lu) with error %d\n", outlen,
 408                (unsigned long)width * lines, ret);
 409         av_free(zbuf);
 410         return AVERROR_UNKNOWN;
 411     }
 412     src = zbuf;
 413     for (line = 0; line < lines; line++) {
 414         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
 415             horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
 416         } else {
 417             memcpy(dst, src, width);
 418         }
 419         if (is_yuv) {
 420             unpack_yuv(s, p, dst, strip_start + line);
 421             line += s->subsampling[1] - 1;
 422         }
 423         dst += stride;
 424         src += width;
 425     }
 426     av_free(zbuf);
 427     return 0;
 428 }
 429 #endif
 430
 431 #if CONFIG_LZMA
 432 static int tiff_uncompress_lzma(uint8_t *dst, uint64_t *len, const uint8_t *src,
 433                                 int size)
 434 {
 435     lzma_stream stream = LZMA_STREAM_INIT;
 436     lzma_ret ret;
 437
 438     stream.next_in   = (uint8_t *)src;
 439     stream.avail_in  = size;
 440     stream.next_out  = dst;
 441     stream.avail_out = *len;
 442     ret              = lzma_stream_decoder(&stream, UINT64_MAX, 0);
 443     if (ret != LZMA_OK) {
 444         av_log(NULL, AV_LOG_ERROR, "LZMA init error: %d\n", ret);
 445         return ret;
 446     }
 447     ret = lzma_code(&stream, LZMA_RUN);
 448     lzma_end(&stream);
 449     *len = stream.total_out;
 450     return ret == LZMA_STREAM_END ? LZMA_OK : ret;
 451 }
 452
 453 static int tiff_unpack_lzma(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
 454                             const uint8_t *src, int size, int width, int lines,
 455                             int strip_start, int is_yuv)
 456 {
 457     uint64_t outlen = width * (uint64_t)lines;
 458     int ret, line;
 459     uint8_t *buf = av_malloc(outlen);
 460     if (!buf)
 461         return AVERROR(ENOMEM);
 462     if (s->fill_order) {
 463         if ((ret = deinvert_buffer(s, src, size)) < 0) {
 464             av_free(buf);
 465             return ret;
 466         }
 467         src = s->deinvert_buf;
 468     }
 469     ret = tiff_uncompress_lzma(buf, &outlen, src, size);
 470     if (ret != LZMA_OK) {
 471         av_log(s->avctx, AV_LOG_ERROR,
 472                "Uncompressing failed (%"PRIu64" of %"PRIu64") with error %d\n", outlen,
 473                (uint64_t)width * lines, ret);
 474         av_free(buf);
 475         return AVERROR_UNKNOWN;
 476     }
 477     src = buf;
 478     for (line = 0; line < lines; line++) {
 479         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
 480             horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
 481         } else {
 482             memcpy(dst, src, width);
 483         }
 484         if (is_yuv) {
 485             unpack_yuv(s, p, dst, strip_start + line);
 486             line += s->subsampling[1] - 1;
 487         }
 488         dst += stride;
 489         src += width;
 490     }
 491     av_free(buf);
 492     return 0;
 493 }
 494 #endif
 495
 496 static int tiff_unpack_fax(TiffContext *s, uint8_t *dst, int stride,
 497                            const uint8_t *src, int size, int width, int lines)
 498 {
 499     int i, ret = 0;
 500     int line;
 501     uint8_t *src2;
 502
 503     av_fast_padded_malloc(&s->fax_buffer, &s->fax_buffer_size, size);
 504     src2 = s->fax_buffer;
 505
 506     if (!src2) {
 507         av_log(s->avctx, AV_LOG_ERROR,
 508                "Error allocating temporary buffer\n");
 509         return AVERROR(ENOMEM);
 510     }
 511
 512     if (!s->fill_order) {
 513         memcpy(src2, src, size);
 514     } else {
 515         for (i = 0; i < size; i++)
 516             src2[i] = ff_reverse[src[i]];
 517     }
 518     memset(src2 + size, 0, AV_INPUT_BUFFER_PADDING_SIZE);
 519     ret = ff_ccitt_unpack(s->avctx, src2, size, dst, lines, stride,
 520                           s->compr, s->fax_opts);
 521     if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
 522         for (line = 0; line < lines; line++) {
 523             horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
 524             dst += stride;
 525         }
 526     return ret;
 527 }
 528
 529 static int tiff_unpack_strip(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
 530                              const uint8_t *src, int size, int strip_start, int lines)
 531 {
 532     PutByteContext pb;
 533     int c, line, pixels, code, ret;
 534     const uint8_t *ssrc = src;
 535     int width = ((s->width * s->bpp) + 7) >> 3;
 536     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(p->format);
 537     int is_yuv = !(desc->flags & AV_PIX_FMT_FLAG_RGB) &&
 538                  (desc->flags & AV_PIX_FMT_FLAG_PLANAR) &&
 539                  desc->nb_components >= 3;
 540
 541     if (s->planar)
 542         width /= s->bppcount;
 543
 544     if (size <= 0)
 545         return AVERROR_INVALIDDATA;
 546
 547     if (is_yuv) {
 548         int bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
 549                             s->subsampling[0] * s->subsampling[1] + 7) >> 3;
 550         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
 551         if (s->yuv_line == NULL) {
 552             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
 553             return AVERROR(ENOMEM);
 554         }
 555         dst = s->yuv_line;
 556         stride = 0;
 557
 558         width = (s->width - 1) / s->subsampling[0] + 1;
 559         width = width * s->subsampling[0] * s->subsampling[1] + 2*width;
 560         av_assert0(width <= bytes_per_row);
 561         av_assert0(s->bpp == 24);
 562     }
 563     if (s->is_bayer) {
 564         width = (s->bpp * s->width + 7) >> 3;
 565     }
 566     if (p->format == AV_PIX_FMT_GRAY12) {
 567         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, width);
 568         if (s->yuv_line == NULL) {
 569             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
 570             return AVERROR(ENOMEM);
 571         }
 572         dst = s->yuv_line;
 573         stride = 0;
 574     }
 575
 576     if (s->compr == TIFF_DEFLATE || s->compr == TIFF_ADOBE_DEFLATE) {
 577 #if CONFIG_ZLIB
 578         return tiff_unpack_zlib(s, p, dst, stride, src, size, width, lines,
 579                                 strip_start, is_yuv);
 580 #else
 581         av_log(s->avctx, AV_LOG_ERROR,
 582                "zlib support not enabled, "
 583                "deflate compression not supported\n");
 584         return AVERROR(ENOSYS);
 585 #endif
 586     }
 587     if (s->compr == TIFF_LZMA) {
 588 #if CONFIG_LZMA
 589         return tiff_unpack_lzma(s, p, dst, stride, src, size, width, lines,
 590                                 strip_start, is_yuv);
 591 #else
 592         av_log(s->avctx, AV_LOG_ERROR,
 593                "LZMA support not enabled\n");
 594         return AVERROR(ENOSYS);
 595 #endif
 596     }
 597     if (s->compr == TIFF_LZW) {
 598         if (s->fill_order) {
 599             if ((ret = deinvert_buffer(s, src, size)) < 0)
 600                 return ret;
 601             ssrc = src = s->deinvert_buf;
 602         }
 603         if (size > 1 && !src[0] && (src[1]&1)) {
 604             av_log(s->avctx, AV_LOG_ERROR, "Old style LZW is unsupported\n");
 605         }
 606         if ((ret = ff_lzw_decode_init(s->lzw, 8, src, size, FF_LZW_TIFF)) < 0) {
 607             av_log(s->avctx, AV_LOG_ERROR, "Error initializing LZW decoder\n");
 608             return ret;
 609         }
 610         for (line = 0; line < lines; line++) {
 611             pixels = ff_lzw_decode(s->lzw, dst, width);
 612             if (pixels < width) {
 613                 av_log(s->avctx, AV_LOG_ERROR, "Decoded only %i bytes of %i\n",
 614                        pixels, width);
 615                 return AVERROR_INVALIDDATA;
 616             }
 617             if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
 618                 horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
 619             if (is_yuv) {
 620                 unpack_yuv(s, p, dst, strip_start + line);
 621                 line += s->subsampling[1] - 1;
 622             } else if (p->format == AV_PIX_FMT_GRAY12) {
 623                 unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
 624             }
 625             dst += stride;
 626         }
 627         return 0;
 628     }
 629     if (s->compr == TIFF_CCITT_RLE ||
 630         s->compr == TIFF_G3        ||
 631         s->compr == TIFF_G4) {
 632         if (is_yuv || p->format == AV_PIX_FMT_GRAY12)
 633             return AVERROR_INVALIDDATA;
 634
 635         return tiff_unpack_fax(s, dst, stride, src, size, width, lines);
 636     }
 637
 638     bytestream2_init(&s->gb, src, size);
 639     bytestream2_init_writer(&pb, dst, is_yuv ? s->yuv_line_size : (stride * lines));
 640
 641     for (line = 0; line < lines; line++) {
 642         if (src - ssrc > size) {
 643             av_log(s->avctx, AV_LOG_ERROR, "Source data overread\n");
 644             return AVERROR_INVALIDDATA;
 645         }
 646
 647         if (bytestream2_get_bytes_left(&s->gb) == 0 || bytestream2_get_eof(&pb))
 648             break;
 649         bytestream2_seek_p(&pb, stride * line, SEEK_SET);
 650         switch (s->compr) {
 651         case TIFF_RAW:
 652             if (ssrc + size - src < width)
 653                 return AVERROR_INVALIDDATA;
 654
 655             if (!s->fill_order) {
 656                 horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 || s->is_bayer),
 657                                 dst, 1, src, 0, width, 0);
 658             } else {
 659                 int i;
 660                 for (i = 0; i < width; i++)
 661                     dst[i] = ff_reverse[src[i]];
 662             }
 663             src += width;
 664             break;
 665         case TIFF_PACKBITS:
 666             for (pixels = 0; pixels < width;) {
 667                 if (ssrc + size - src < 2) {
 668                     av_log(s->avctx, AV_LOG_ERROR, "Read went out of bounds\n");
 669                     return AVERROR_INVALIDDATA;
 670                 }
 671                 code = s->fill_order ? (int8_t) ff_reverse[*src++]: (int8_t) *src++;
 672                 if (code >= 0) {
 673                     code++;
 674                     if (pixels + code > width ||
 675                         ssrc + size - src < code) {
 676                         av_log(s->avctx, AV_LOG_ERROR,
 677                                "Copy went out of bounds\n");
 678                         return AVERROR_INVALIDDATA;
 679                     }
 680                     horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
 681                                     dst, 1, src, 0, code, pixels);
 682                     src    += code;
 683                     pixels += code;
 684                 } else if (code != -128) { // -127..-1
 685                     code = (-code) + 1;
 686                     if (pixels + code > width) {
 687                         av_log(s->avctx, AV_LOG_ERROR,
 688                                "Run went out of bounds\n");
 689                         return AVERROR_INVALIDDATA;
 690                     }
 691                     c = *src++;
 692                     horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
 693                                     dst, 0, NULL, c, code, pixels);
 694                     pixels += code;
 695                 }
 696             }
 697             if (s->fill_order) {
 698                 int i;
 699                 for (i = 0; i < width; i++)
 700                     dst[i] = ff_reverse[dst[i]];
 701             }
 702             break;
 703         }
 704         if (is_yuv) {
 705             unpack_yuv(s, p, dst, strip_start + line);
 706             line += s->subsampling[1] - 1;
 707         } else if (p->format == AV_PIX_FMT_GRAY12) {
 708             unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
 709         }
 710         dst += stride;
 711     }
 712     return 0;
 713 }
 714
 715 static int init_image(TiffContext *s, ThreadFrame *frame)
 716 {
 717     int ret;
 718     int create_gray_palette = 0;
 719
 720     // make sure there is no aliasing in the following switch
 721     if (s->bpp >= 100 || s->bppcount >= 10) {
 722         av_log(s->avctx, AV_LOG_ERROR,
 723                "Unsupported image parameters: bpp=%d, bppcount=%d\n",
 724                s->bpp, s->bppcount);
 725         return AVERROR_INVALIDDATA;
 726     }
 727
 728     switch (s->planar * 1000 + s->bpp * 10 + s->bppcount + s->is_bayer * 10000) {
 729     case 11:
 730         if (!s->palette_is_set) {
 731             s->avctx->pix_fmt = AV_PIX_FMT_MONOBLACK;
 732             break;
 733         }
 734     case 21:
 735     case 41:
 736         s->avctx->pix_fmt = AV_PIX_FMT_PAL8;
 737         if (!s->palette_is_set) {
 738             create_gray_palette = 1;
 739         }
 740         break;
 741     case 81:
 742         s->avctx->pix_fmt = s->palette_is_set ? AV_PIX_FMT_PAL8 : AV_PIX_FMT_GRAY8;
 743         break;
 744     case 121:
 745         s->avctx->pix_fmt = AV_PIX_FMT_GRAY12;
 746         break;
 747     case 10081:
 748         switch (AV_RL32(s->pattern)) {
 749         case 0x02010100:
 750             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_RGGB8;
 751             break;
 752         case 0x00010102:
 753             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_BGGR8;
 754             break;
 755         case 0x01000201:
 756             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GBRG8;
 757             break;
 758         case 0x01020001:
 759             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GRBG8;
 760             break;
 761         default:
 762             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
 763                    AV_RL32(s->pattern));
 764             return AVERROR_PATCHWELCOME;
 765         }
 766         break;
 767     case 10121:
 768         switch (AV_RL32(s->pattern)) {
 769         case 0x02010100:
 770             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_RGGB16LE : AV_PIX_FMT_BAYER_RGGB16BE;
 771             break;
 772         case 0x00010102:
 773             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_BGGR16LE : AV_PIX_FMT_BAYER_BGGR16BE;
 774             break;
 775         case 0x01000201:
 776             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GBRG16LE : AV_PIX_FMT_BAYER_GBRG16BE;
 777             break;
 778         case 0x01020001:
 779             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GRBG16LE : AV_PIX_FMT_BAYER_GRBG16BE;
 780             break;
 781         default:
 782             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
 783                    AV_RL32(s->pattern));
 784             return AVERROR_PATCHWELCOME;
 785         }
 786         break;
 787     case 10161:
 788         switch (AV_RL32(s->pattern)) {
 789         case 0x02010100:
 790             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_RGGB16LE : AV_PIX_FMT_BAYER_RGGB16BE;
 791             break;
 792         case 0x00010102:
 793             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_BGGR16LE : AV_PIX_FMT_BAYER_BGGR16BE;
 794             break;
 795         case 0x01000201:
 796             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GBRG16LE : AV_PIX_FMT_BAYER_GBRG16BE;
 797             break;
 798         case 0x01020001:
 799             s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GRBG16LE : AV_PIX_FMT_BAYER_GRBG16BE;
 800             break;
 801         default:
 802             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
 803                    AV_RL32(s->pattern));
 804             return AVERROR_PATCHWELCOME;
 805         }
 806         break;
 807     case 243:
 808         if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
 809             if (s->subsampling[0] == 1 && s->subsampling[1] == 1) {
 810                 s->avctx->pix_fmt = AV_PIX_FMT_YUV444P;
 811             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 1) {
 812                 s->avctx->pix_fmt = AV_PIX_FMT_YUV422P;
 813             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 1) {
 814                 s->avctx->pix_fmt = AV_PIX_FMT_YUV411P;
 815             } else if (s->subsampling[0] == 1 && s->subsampling[1] == 2) {
 816                 s->avctx->pix_fmt = AV_PIX_FMT_YUV440P;
 817             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 2) {
 818                 s->avctx->pix_fmt = AV_PIX_FMT_YUV420P;
 819             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 4) {
 820                 s->avctx->pix_fmt = AV_PIX_FMT_YUV410P;
 821             } else {
 822                 av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr subsampling\n");
 823                 return AVERROR_PATCHWELCOME;
 824             }
 825         } else
 826             s->avctx->pix_fmt = AV_PIX_FMT_RGB24;
 827         break;
 828     case 161:
 829         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GRAY16LE : AV_PIX_FMT_GRAY16BE;
 830         break;
 831     case 162:
 832         s->avctx->pix_fmt = AV_PIX_FMT_YA8;
 833         break;
 834     case 322:
 835         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_YA16LE : AV_PIX_FMT_YA16BE;
 836         break;
 837     case 324:
 838         s->avctx->pix_fmt = s->photometric == TIFF_PHOTOMETRIC_SEPARATED ? AV_PIX_FMT_RGB0 : AV_PIX_FMT_RGBA;
 839         break;
 840     case 405:
 841         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED)
 842             s->avctx->pix_fmt = AV_PIX_FMT_RGBA;
 843         else {
 844             av_log(s->avctx, AV_LOG_ERROR,
 845                 "bpp=40 without PHOTOMETRIC_SEPARATED is unsupported\n");
 846             return AVERROR_PATCHWELCOME;
 847         }
 848         break;
 849     case 483:
 850         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGB48LE  : AV_PIX_FMT_RGB48BE;
 851         break;
 852     case 644:
 853         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGBA64LE  : AV_PIX_FMT_RGBA64BE;
 854         break;
 855     case 1243:
 856         s->avctx->pix_fmt = AV_PIX_FMT_GBRP;
 857         break;
 858     case 1324:
 859         s->avctx->pix_fmt = AV_PIX_FMT_GBRAP;
 860         break;
 861     case 1483:
 862         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRP16LE : AV_PIX_FMT_GBRP16BE;
 863         break;
 864     case 1644:
 865         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRAP16LE : AV_PIX_FMT_GBRAP16BE;
 866         break;
 867     default:
 868         av_log(s->avctx, AV_LOG_ERROR,
 869                "This format is not supported (bpp=%d, bppcount=%d)\n",
 870                s->bpp, s->bppcount);
 871         return AVERROR_INVALIDDATA;
 872     }
 873
 874     if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
 875         const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(s->avctx->pix_fmt);
 876         if((desc->flags & AV_PIX_FMT_FLAG_RGB) ||
 877            !(desc->flags & AV_PIX_FMT_FLAG_PLANAR) ||
 878            desc->nb_components < 3) {
 879             av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr variant\n");
 880             return AVERROR_INVALIDDATA;
 881         }
 882     }
 883
 884     if (s->width != s->avctx->width || s->height != s->avctx->height) {
 885         ret = ff_set_dimensions(s->avctx, s->width, s->height);
 886         if (ret < 0)
 887             return ret;
 888     }
 889     if ((ret = ff_thread_get_buffer(s->avctx, frame, 0)) < 0)
 890         return ret;
 891     if (s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
 892         if (!create_gray_palette)
 893             memcpy(frame->f->data[1], s->palette, sizeof(s->palette));
 894         else {
 895             /* make default grayscale pal */
 896             int i;
 897             uint32_t *pal = (uint32_t *)frame->f->data[1];
 898             for (i = 0; i < 1<<s->bpp; i++)
 899                 pal[i] = 0xFFU << 24 | i * 255 / ((1<<s->bpp) - 1) * 0x010101;
 900         }
 901     }
 902     return 0;
 903 }
 904
 905 static void set_sar(TiffContext *s, unsigned tag, unsigned num, unsigned den)
 906 {
 907     int offset = tag == TIFF_YRES ? 2 : 0;
 908     s->res[offset++] = num;
 909     s->res[offset]   = den;
 910     if (s->res[0] && s->res[1] && s->res[2] && s->res[3]) {
 911         uint64_t num = s->res[2] * (uint64_t)s->res[1];
 912         uint64_t den = s->res[0] * (uint64_t)s->res[3];
 913         if (num > INT64_MAX || den > INT64_MAX) {
 914             num = num >> 1;
 915             den = den >> 1;
 916         }
 917         av_reduce(&s->avctx->sample_aspect_ratio.num, &s->avctx->sample_aspect_ratio.den,
 918                   num, den, INT32_MAX);
 919         if (!s->avctx->sample_aspect_ratio.den)
 920             s->avctx->sample_aspect_ratio = (AVRational) {0, 1};
 921     }
 922 }
 923
 924 static int tiff_decode_tag(TiffContext *s, AVFrame *frame)
 925 {
 926     unsigned tag, type, count, off, value = 0, value2 = 0;
 927     int i, start;
 928     int pos;
 929     int ret;
 930     double *dp;
 931
 932     ret = ff_tread_tag(&s->gb, s->le, &tag, &type, &count, &start);
 933     if (ret < 0) {
 934         goto end;
 935     }
 936
 937     off = bytestream2_tell(&s->gb);
 938     if (count == 1) {
 939         switch (type) {
 940         case TIFF_BYTE:
 941         case TIFF_SHORT:
 942         case TIFF_LONG:
 943             value = ff_tget(&s->gb, type, s->le);
 944             break;
 945         case TIFF_RATIONAL:
 946             value  = ff_tget(&s->gb, TIFF_LONG, s->le);
 947             value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
 948             break;
 949         case TIFF_STRING:
 950             if (count <= 4) {
 951                 break;
 952             }
 953         default:
 954             value = UINT_MAX;
 955         }
 956     }
 957
 958     switch (tag) {
 959     case TIFF_SUBFILE:
 960         s->is_thumbnail = (value != 0);
 961         break;
 962     case TIFF_WIDTH:
 963         s->width = value;
 964         break;
 965     case TIFF_HEIGHT:
 966         s->height = value;
 967         break;
 968     case TIFF_BPP:
 969         if (count > 5U) {
 970             av_log(s->avctx, AV_LOG_ERROR,
 971                    "This format is not supported (bpp=%d, %d components)\n",
 972                    value, count);
 973             return AVERROR_INVALIDDATA;
 974         }
 975         s->bppcount = count;
 976         if (count == 1)
 977             s->bpp = value;
 978         else {
 979             switch (type) {
 980             case TIFF_BYTE:
 981             case TIFF_SHORT:
 982             case TIFF_LONG:
 983                 s->bpp = 0;
 984                 if (bytestream2_get_bytes_left(&s->gb) < type_sizes[type] * count)
 985                     return AVERROR_INVALIDDATA;
 986                 for (i = 0; i < count; i++)
 987                     s->bpp += ff_tget(&s->gb, type, s->le);
 988                 break;
 989             default:
 990                 s->bpp = -1;
 991             }
 992         }
 993         break;
 994     case TIFF_SAMPLES_PER_PIXEL:
 995         if (count != 1) {
 996             av_log(s->avctx, AV_LOG_ERROR,
 997                    "Samples per pixel requires a single value, many provided\n");
 998             return AVERROR_INVALIDDATA;
 999         }
1000         if (value > 5U) {
1001             av_log(s->avctx, AV_LOG_ERROR,
1002                    "Samples per pixel %d is too large\n", value);
1003             return AVERROR_INVALIDDATA;
1004         }
1005         if (s->bppcount == 1)
1006             s->bpp *= value;
1007         s->bppcount = value;
1008         break;
1009     case TIFF_COMPR:
1010         s->compr     = value;
1011         av_log(s->avctx, AV_LOG_DEBUG, "compression: %d\n", s->compr);
1012         s->predictor = 0;
1013         switch (s->compr) {
1014         case TIFF_RAW:
1015         case TIFF_PACKBITS:
1016         case TIFF_LZW:
1017         case TIFF_CCITT_RLE:
1018             break;
1019         case TIFF_G3:
1020         case TIFF_G4:
1021             s->fax_opts = 0;
1022             break;
1023         case TIFF_DEFLATE:
1024         case TIFF_ADOBE_DEFLATE:
1025 #if CONFIG_ZLIB
1026             break;
1027 #else
1028             av_log(s->avctx, AV_LOG_ERROR, "Deflate: ZLib not compiled in\n");
1029             return AVERROR(ENOSYS);
1030 #endif
1031         case TIFF_JPEG:
1032         case TIFF_NEWJPEG:
1033             avpriv_report_missing_feature(s->avctx, "JPEG compression");
1034             return AVERROR_PATCHWELCOME;
1035         case TIFF_LZMA:
1036 #if CONFIG_LZMA
1037             break;
1038 #else
1039             av_log(s->avctx, AV_LOG_ERROR, "LZMA not compiled in\n");
1040             return AVERROR(ENOSYS);
1041 #endif
1042         default:
1043             av_log(s->avctx, AV_LOG_ERROR, "Unknown compression method %i\n",
1044                    s->compr);
1045             return AVERROR_INVALIDDATA;
1046         }
1047         break;
1048     case TIFF_ROWSPERSTRIP:
1049         if (!value || (type == TIFF_LONG && value == UINT_MAX))
1050             value = s->height;
1051         s->rps = FFMIN(value, s->height);
1052         break;
1053     case TIFF_STRIP_OFFS:
1054         if (count == 1) {
1055             if (value > INT_MAX) {
1056                 av_log(s->avctx, AV_LOG_ERROR,
1057                     "strippos %u too large\n", value);
1058                 return AVERROR_INVALIDDATA;
1059             }
1060             s->strippos = 0;
1061             s->stripoff = value;
1062         } else
1063             s->strippos = off;
1064         s->strips = count;
1065         if (s->strips == 1)
1066             s->rps = s->height;
1067         s->sot = type;
1068         break;
1069     case TIFF_STRIP_SIZE:
1070         if (count == 1) {
1071             if (value > INT_MAX) {
1072                 av_log(s->avctx, AV_LOG_ERROR,
1073                     "stripsize %u too large\n", value);
1074                 return AVERROR_INVALIDDATA;
1075             }
1076             s->stripsizesoff = 0;
1077             s->stripsize     = value;
1078             s->strips        = 1;
1079         } else {
1080             s->stripsizesoff = off;
1081         }
1082         s->strips = count;
1083         s->sstype = type;
1084         break;
1085     case TIFF_XRES:
1086     case TIFF_YRES:
1087         set_sar(s, tag, value, value2);
1088         break;
1089     case TIFF_TILE_BYTE_COUNTS:
1090     case TIFF_TILE_LENGTH:
1091     case TIFF_TILE_OFFSETS:
1092     case TIFF_TILE_WIDTH:
1093         av_log(s->avctx, AV_LOG_ERROR, "Tiled images are not supported\n");
1094         return AVERROR_PATCHWELCOME;
1095         break;
1096     case TIFF_PREDICTOR:
1097         s->predictor = value;
1098         break;
1099     case TIFF_SUB_IFDS:
1100         if (count == 1)
1101             s->sub_ifd = value;
1102         else if (count > 1)
1103             s->sub_ifd = ff_tget(&s->gb, TIFF_LONG, s->le); /** Only get the first SubIFD */
1104         break;
1105     case DNG_WHITE_LEVEL:
1106         s->white_level = value;
1107         break;
1108     case TIFF_CFA_PATTERN_DIM:
1109         if (count != 2 || (ff_tget(&s->gb, type, s->le) != 2 &&
1110                            ff_tget(&s->gb, type, s->le) != 2)) {
1111             av_log(s->avctx, AV_LOG_ERROR, "CFA Pattern dimensions are not 2x2\n");
1112             return AVERROR_INVALIDDATA;
1113         }
1114         break;
1115     case TIFF_CFA_PATTERN:
1116         s->is_bayer = 1;
1117         s->pattern[0] = ff_tget(&s->gb, type, s->le);
1118         s->pattern[1] = ff_tget(&s->gb, type, s->le);
1119         s->pattern[2] = ff_tget(&s->gb, type, s->le);
1120         s->pattern[3] = ff_tget(&s->gb, type, s->le);
1121         break;
1122     case TIFF_PHOTOMETRIC:
1123         switch (value) {
1124         case TIFF_PHOTOMETRIC_WHITE_IS_ZERO:
1125         case TIFF_PHOTOMETRIC_BLACK_IS_ZERO:
1126         case TIFF_PHOTOMETRIC_RGB:
1127         case TIFF_PHOTOMETRIC_PALETTE:
1128         case TIFF_PHOTOMETRIC_SEPARATED:
1129         case TIFF_PHOTOMETRIC_YCBCR:
1130         case TIFF_PHOTOMETRIC_CFA:
1131             s->photometric = value;
1132             break;
1133         case TIFF_PHOTOMETRIC_ALPHA_MASK:
1134         case TIFF_PHOTOMETRIC_CIE_LAB:
1135         case TIFF_PHOTOMETRIC_ICC_LAB:
1136         case TIFF_PHOTOMETRIC_ITU_LAB:
1137         case TIFF_PHOTOMETRIC_LOG_L:
1138         case TIFF_PHOTOMETRIC_LOG_LUV:
1139         case TIFF_PHOTOMETRIC_LINEAR_RAW:
1140             avpriv_report_missing_feature(s->avctx,
1141                                           "PhotometricInterpretation 0x%04X",
1142                                           value);
1143             return AVERROR_PATCHWELCOME;
1144         default:
1145             av_log(s->avctx, AV_LOG_ERROR, "PhotometricInterpretation %u is "
1146                    "unknown\n", value);
1147             return AVERROR_INVALIDDATA;
1148         }
1149         break;
1150     case TIFF_FILL_ORDER:
1151         if (value < 1 || value > 2) {
1152             av_log(s->avctx, AV_LOG_ERROR,
1153                    "Unknown FillOrder value %d, trying default one\n", value);
1154             value = 1;
1155         }
1156         s->fill_order = value - 1;
1157         break;
1158     case TIFF_PAL: {
1159         GetByteContext pal_gb[3];
1160         off = type_sizes[type];
1161         if (count / 3 > 256 ||
1162             bytestream2_get_bytes_left(&s->gb) < count / 3 * off * 3)
1163             return AVERROR_INVALIDDATA;
1164
1165         pal_gb[0] = pal_gb[1] = pal_gb[2] = s->gb;
1166         bytestream2_skip(&pal_gb[1], count / 3 * off);
1167         bytestream2_skip(&pal_gb[2], count / 3 * off * 2);
1168
1169         off = (type_sizes[type] - 1) << 3;
1170         if (off > 31U) {
1171             av_log(s->avctx, AV_LOG_ERROR, "palette shift %d is out of range\n", off);
1172             return AVERROR_INVALIDDATA;
1173         }
1174
1175         for (i = 0; i < count / 3; i++) {
1176             uint32_t p = 0xFF000000;
1177             p |= (ff_tget(&pal_gb[0], type, s->le) >> off) << 16;
1178             p |= (ff_tget(&pal_gb[1], type, s->le) >> off) << 8;
1179             p |=  ff_tget(&pal_gb[2], type, s->le) >> off;
1180             s->palette[i] = p;
1181         }
1182         s->palette_is_set = 1;
1183         break;
1184     }
1185     case TIFF_PLANAR:
1186         s->planar = value == 2;
1187         break;
1188     case TIFF_YCBCR_SUBSAMPLING:
1189         if (count != 2) {
1190             av_log(s->avctx, AV_LOG_ERROR, "subsample count invalid\n");
1191             return AVERROR_INVALIDDATA;
1192         }
1193         for (i = 0; i < count; i++) {
1194             s->subsampling[i] = ff_tget(&s->gb, type, s->le);
1195             if (s->subsampling[i] <= 0) {
1196                 av_log(s->avctx, AV_LOG_ERROR, "subsampling %d is invalid\n", s->subsampling[i]);
1197                 s->subsampling[i] = 1;
1198                 return AVERROR_INVALIDDATA;
1199             }
1200         }
1201         break;
1202     case TIFF_T4OPTIONS:
1203         if (s->compr == TIFF_G3)
1204             s->fax_opts = value;
1205         break;
1206     case TIFF_T6OPTIONS:
1207         if (s->compr == TIFF_G4)
1208             s->fax_opts = value;
1209         break;
1210 #define ADD_METADATA(count, name, sep)\
1211     if ((ret = add_metadata(count, type, name, sep, s, frame)) < 0) {\
1212         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");\
1213         goto end;\
1214     }
1215     case TIFF_MODEL_PIXEL_SCALE:
1216         ADD_METADATA(count, "ModelPixelScaleTag", NULL);
1217         break;
1218     case TIFF_MODEL_TRANSFORMATION:
1219         ADD_METADATA(count, "ModelTransformationTag", NULL);
1220         break;
1221     case TIFF_MODEL_TIEPOINT:
1222         ADD_METADATA(count, "ModelTiepointTag", NULL);
1223         break;
1224     case TIFF_GEO_KEY_DIRECTORY:
1225         if (s->geotag_count) {
1226             avpriv_request_sample(s->avctx, "Multiple geo key directories\n");
1227             return AVERROR_INVALIDDATA;
1228         }
1229         ADD_METADATA(1, "GeoTIFF_Version", NULL);
1230         ADD_METADATA(2, "GeoTIFF_Key_Revision", ".");
1231         s->geotag_count   = ff_tget_short(&s->gb, s->le);
1232         if (s->geotag_count > count / 4 - 1) {
1233             s->geotag_count = count / 4 - 1;
1234             av_log(s->avctx, AV_LOG_WARNING, "GeoTIFF key directory buffer shorter than specified\n");
1235         }
1236         if (   bytestream2_get_bytes_left(&s->gb) < s->geotag_count * sizeof(int16_t) * 4
1237             || s->geotag_count == 0) {
1238             s->geotag_count = 0;
1239             return -1;
1240         }
1241         s->geotags = av_mallocz_array(s->geotag_count, sizeof(TiffGeoTag));
1242         if (!s->geotags) {
1243             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1244             s->geotag_count = 0;
1245             goto end;
1246         }
1247         for (i = 0; i < s->geotag_count; i++) {
1248             s->geotags[i].key    = ff_tget_short(&s->gb, s->le);
1249             s->geotags[i].type   = ff_tget_short(&s->gb, s->le);
1250             s->geotags[i].count  = ff_tget_short(&s->gb, s->le);
1251
1252             if (!s->geotags[i].type)
1253                 s->geotags[i].val  = get_geokey_val(s->geotags[i].key, ff_tget_short(&s->gb, s->le));
1254             else
1255                 s->geotags[i].offset = ff_tget_short(&s->gb, s->le);
1256         }
1257         break;
1258     case TIFF_GEO_DOUBLE_PARAMS:
1259         if (count >= INT_MAX / sizeof(int64_t))
1260             return AVERROR_INVALIDDATA;
1261         if (bytestream2_get_bytes_left(&s->gb) < count * sizeof(int64_t))
1262             return AVERROR_INVALIDDATA;
1263         dp = av_malloc_array(count, sizeof(double));
1264         if (!dp) {
1265             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1266             goto end;
1267         }
1268         for (i = 0; i < count; i++)
1269             dp[i] = ff_tget_double(&s->gb, s->le);
1270         for (i = 0; i < s->geotag_count; i++) {
1271             if (s->geotags[i].type == TIFF_GEO_DOUBLE_PARAMS) {
1272                 if (s->geotags[i].count == 0
1273                     || s->geotags[i].offset + s->geotags[i].count > count) {
1274                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1275                 } else if (s->geotags[i].val) {
1276                     av_log(s->avctx, AV_LOG_WARNING, "Duplicate GeoTIFF key %d\n", s->geotags[i].key);
1277                 } else {
1278                     char *ap = doubles2str(&dp[s->geotags[i].offset], s->geotags[i].count, ", ");
1279                     if (!ap) {
1280                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1281                         av_freep(&dp);
1282                         return AVERROR(ENOMEM);
1283                     }
1284                     s->geotags[i].val = ap;
1285                 }
1286             }
1287         }
1288         av_freep(&dp);
1289         break;
1290     case TIFF_GEO_ASCII_PARAMS:
1291         pos = bytestream2_tell(&s->gb);
1292         for (i = 0; i < s->geotag_count; i++) {
1293             if (s->geotags[i].type == TIFF_GEO_ASCII_PARAMS) {
1294                 if (s->geotags[i].count == 0
1295                     || s->geotags[i].offset +  s->geotags[i].count > count) {
1296                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1297                 } else {
1298                     char *ap;
1299
1300                     bytestream2_seek(&s->gb, pos + s->geotags[i].offset, SEEK_SET);
1301                     if (bytestream2_get_bytes_left(&s->gb) < s->geotags[i].count)
1302                         return AVERROR_INVALIDDATA;
1303                     if (s->geotags[i].val)
1304                         return AVERROR_INVALIDDATA;
1305                     ap = av_malloc(s->geotags[i].count);
1306                     if (!ap) {
1307                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1308                         return AVERROR(ENOMEM);
1309                     }
1310                     bytestream2_get_bufferu(&s->gb, ap, s->geotags[i].count);
1311                     ap[s->geotags[i].count - 1] = '\0'; //replace the "|" delimiter with a 0 byte
1312                     s->geotags[i].val = ap;
1313                 }
1314             }
1315         }
1316         break;
1317     case TIFF_ARTIST:
1318         ADD_METADATA(count, "artist", NULL);
1319         break;
1320     case TIFF_COPYRIGHT:
1321         ADD_METADATA(count, "copyright", NULL);
1322         break;
1323     case TIFF_DATE:
1324         ADD_METADATA(count, "date", NULL);
1325         break;
1326     case TIFF_DOCUMENT_NAME:
1327         ADD_METADATA(count, "document_name", NULL);
1328         break;
1329     case TIFF_HOST_COMPUTER:
1330         ADD_METADATA(count, "computer", NULL);
1331         break;
1332     case TIFF_IMAGE_DESCRIPTION:
1333         ADD_METADATA(count, "description", NULL);
1334         break;
1335     case TIFF_MAKE:
1336         ADD_METADATA(count, "make", NULL);
1337         break;
1338     case TIFF_MODEL:
1339         ADD_METADATA(count, "model", NULL);
1340         break;
1341     case TIFF_PAGE_NAME:
1342         ADD_METADATA(count, "page_name", NULL);
1343         break;
1344     case TIFF_PAGE_NUMBER:
1345         ADD_METADATA(count, "page_number", " / ");
1346         // need to seek back to re-read the page number
1347         bytestream2_seek(&s->gb, -count * sizeof(uint16_t), SEEK_CUR);
1348         // read the page number
1349         s->cur_page = ff_tget(&s->gb, TIFF_SHORT, s->le);
1350         // get back to where we were before the previous seek
1351         bytestream2_seek(&s->gb, count * sizeof(uint16_t) - sizeof(uint16_t), SEEK_CUR);
1352         break;
1353     case TIFF_SOFTWARE_NAME:
1354         ADD_METADATA(count, "software", NULL);
1355         break;
1356     case DNG_VERSION:
1357         if (count == 4) {
1358             unsigned int ver[4];
1359             ver[0] = ff_tget(&s->gb, type, s->le);
1360             ver[1] = ff_tget(&s->gb, type, s->le);
1361             ver[2] = ff_tget(&s->gb, type, s->le);
1362             ver[3] = ff_tget(&s->gb, type, s->le);
1363
1364             av_log(s->avctx, AV_LOG_DEBUG, "DNG file, version %u.%u.%u.%u\n",
1365                 ver[0], ver[1], ver[2], ver[3]);
1366
1367             tiff_set_type(s, TIFF_TYPE_DNG);
1368         }
1369         break;
1370     case CINEMADNG_TIME_CODES:
1371     case CINEMADNG_FRAME_RATE:
1372     case CINEMADNG_T_STOP:
1373     case CINEMADNG_REEL_NAME:
1374     case CINEMADNG_CAMERA_LABEL:
1375         tiff_set_type(s, TIFF_TYPE_CINEMADNG);
1376         break;
1377     default:
1378         if (s->avctx->err_recognition & AV_EF_EXPLODE) {
1379             av_log(s->avctx, AV_LOG_ERROR,
1380                    "Unknown or unsupported tag %d/0x%0X\n",
1381                    tag, tag);
1382             return AVERROR_INVALIDDATA;
1383         }
1384     }
1385 end:
1386     if (s->bpp > 64U) {
1387         av_log(s->avctx, AV_LOG_ERROR,
1388                 "This format is not supported (bpp=%d, %d components)\n",
1389                 s->bpp, count);
1390         s->bpp = 0;
1391         return AVERROR_INVALIDDATA;
1392     }
1393     bytestream2_seek(&s->gb, start, SEEK_SET);
1394     return 0;
1395 }
1396
1397 static int decode_frame(AVCodecContext *avctx,
1398                         void *data, int *got_frame, AVPacket *avpkt)
1399 {
1400     TiffContext *const s = avctx->priv_data;
1401     AVFrame *const p = data;
1402     ThreadFrame frame = { .f = data };
1403     unsigned off, last_off;
1404     int le, ret, plane, planes;
1405     int i, j, entries, stride;
1406     unsigned soff, ssize;
1407     uint8_t *dst;
1408     GetByteContext stripsizes;
1409     GetByteContext stripdata;
1410     int retry_for_subifd, retry_for_page;
1411
1412     bytestream2_init(&s->gb, avpkt->data, avpkt->size);
1413
1414     // parse image header
1415     if ((ret = ff_tdecode_header(&s->gb, &le, &off))) {
1416         av_log(avctx, AV_LOG_ERROR, "Invalid TIFF header\n");
1417         return ret;
1418     } else if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1419         av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1420         return AVERROR_INVALIDDATA;
1421     }
1422     s->le          = le;
1423     // TIFF_BPP is not a required tag and defaults to 1
1424 again:
1425     s->is_thumbnail = 0;
1426     s->bppcount    = s->bpp = 1;
1427     s->photometric = TIFF_PHOTOMETRIC_NONE;
1428     s->compr       = TIFF_RAW;
1429     s->fill_order  = 0;
1430     s->white_level = 0;
1431     s->is_bayer    = 0;
1432     s->cur_page    = 0;
1433     s->tiff_type   = TIFF_TYPE_TIFF;
1434     free_geotags(s);
1435
1436     // Reset these offsets so we can tell if they were set this frame
1437     s->stripsizesoff = s->strippos = 0;
1438     /* parse image file directory */
1439     bytestream2_seek(&s->gb, off, SEEK_SET);
1440     entries = ff_tget_short(&s->gb, le);
1441     if (bytestream2_get_bytes_left(&s->gb) < entries * 12)
1442         return AVERROR_INVALIDDATA;
1443     for (i = 0; i < entries; i++) {
1444         if ((ret = tiff_decode_tag(s, p)) < 0)
1445             return ret;
1446     }
1447
1448     if (s->get_thumbnail && !s->is_thumbnail) {
1449         av_log(avctx, AV_LOG_INFO, "No embedded thumbnail present\n");
1450         return AVERROR_EOF;
1451     }
1452
1453     /** whether we should process this IFD's SubIFD */
1454     retry_for_subifd = s->sub_ifd && (s->get_subimage || (!s->get_thumbnail && s->is_thumbnail));
1455     /** whether we should process this multi-page IFD's next page */
1456     retry_for_page = s->get_page && s->cur_page + 1 < s->get_page;  // get_page is 1-indexed
1457
1458     last_off = off;
1459     if (retry_for_page) {
1460         // set offset to the next IFD
1461         off = ff_tget_long(&s->gb, le);
1462     } else if (retry_for_subifd) {
1463         // set offset to the SubIFD
1464         off = s->sub_ifd;
1465     }
1466
1467     if (retry_for_subifd || retry_for_page) {
1468         if (!off) {
1469             av_log(avctx, AV_LOG_ERROR, "Requested entry not found\n");
1470             return AVERROR_INVALIDDATA;
1471         }
1472         if (off <= last_off) {
1473             avpriv_request_sample(s->avctx, "non increasing IFD offset\n");
1474             return AVERROR_INVALIDDATA;
1475         }
1476         if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1477             av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1478             return AVERROR_INVALIDDATA;
1479         }
1480         s->sub_ifd = 0;
1481         goto again;
1482     }
1483
1484     for (i = 0; i<s->geotag_count; i++) {
1485         const char *keyname = get_geokey_name(s->geotags[i].key);
1486         if (!keyname) {
1487             av_log(avctx, AV_LOG_WARNING, "Unknown or unsupported GeoTIFF key %d\n", s->geotags[i].key);
1488             continue;
1489         }
1490         if (get_geokey_type(s->geotags[i].key) != s->geotags[i].type) {
1491             av_log(avctx, AV_LOG_WARNING, "Type of GeoTIFF key %d is wrong\n", s->geotags[i].key);
1492             continue;
1493         }
1494         ret = av_dict_set(&p->metadata, keyname, s->geotags[i].val, 0);
1495         if (ret<0) {
1496             av_log(avctx, AV_LOG_ERROR, "Writing metadata with key '%s' failed\n", keyname);
1497             return ret;
1498         }
1499     }
1500
1501     if (!s->strippos && !s->stripoff) {
1502         av_log(avctx, AV_LOG_ERROR, "Image data is missing\n");
1503         return AVERROR_INVALIDDATA;
1504     }
1505     /* now we have the data and may start decoding */
1506     if ((ret = init_image(s, &frame)) < 0)
1507         return ret;
1508
1509     if (s->strips == 1 && !s->stripsize) {
1510         av_log(avctx, AV_LOG_WARNING, "Image data size missing\n");
1511         s->stripsize = avpkt->size - s->stripoff;
1512     }
1513
1514     if (s->stripsizesoff) {
1515         if (s->stripsizesoff >= (unsigned)avpkt->size)
1516             return AVERROR_INVALIDDATA;
1517         bytestream2_init(&stripsizes, avpkt->data + s->stripsizesoff,
1518                          avpkt->size - s->stripsizesoff);
1519     }
1520     if (s->strippos) {
1521         if (s->strippos >= (unsigned)avpkt->size)
1522             return AVERROR_INVALIDDATA;
1523         bytestream2_init(&stripdata, avpkt->data + s->strippos,
1524                          avpkt->size - s->strippos);
1525     }
1526
1527     if (s->rps <= 0 || s->rps % s->subsampling[1]) {
1528         av_log(avctx, AV_LOG_ERROR, "rps %d invalid\n", s->rps);
1529         return AVERROR_INVALIDDATA;
1530     }
1531
1532     planes = s->planar ? s->bppcount : 1;
1533     for (plane = 0; plane < planes; plane++) {
1534         uint8_t *five_planes = NULL;
1535         int remaining = avpkt->size;
1536         int decoded_height;
1537         stride = p->linesize[plane];
1538         dst = p->data[plane];
1539         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1540             s->avctx->pix_fmt == AV_PIX_FMT_RGBA) {
1541             stride = stride * 5 / 4;
1542             five_planes =
1543             dst = av_malloc(stride * s->height);
1544             if (!dst)
1545                 return AVERROR(ENOMEM);
1546         }
1547         for (i = 0; i < s->height; i += s->rps) {
1548             if (i)
1549                 dst += s->rps * stride;
1550             if (s->stripsizesoff)
1551                 ssize = ff_tget(&stripsizes, s->sstype, le);
1552             else
1553                 ssize = s->stripsize;
1554
1555             if (s->strippos)
1556                 soff = ff_tget(&stripdata, s->sot, le);
1557             else
1558                 soff = s->stripoff;
1559
1560             if (soff > avpkt->size || ssize > avpkt->size - soff || ssize > remaining) {
1561                 av_log(avctx, AV_LOG_ERROR, "Invalid strip size/offset\n");
1562                 av_freep(&five_planes);
1563                 return AVERROR_INVALIDDATA;
1564             }
1565             remaining -= ssize;
1566             if ((ret = tiff_unpack_strip(s, p, dst, stride, avpkt->data + soff, ssize, i,
1567                                          FFMIN(s->rps, s->height - i))) < 0) {
1568                 if (avctx->err_recognition & AV_EF_EXPLODE) {
1569                     av_freep(&five_planes);
1570                     return ret;
1571                 }
1572                 break;
1573             }
1574         }
1575         decoded_height = FFMIN(i, s->height);
1576
1577         if (s->predictor == 2) {
1578             if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1579                 av_log(s->avctx, AV_LOG_ERROR, "predictor == 2 with YUV is unsupported");
1580                 return AVERROR_PATCHWELCOME;
1581             }
1582             dst   = five_planes ? five_planes : p->data[plane];
1583             soff  = s->bpp >> 3;
1584             if (s->planar)
1585                 soff  = FFMAX(soff / s->bppcount, 1);
1586             ssize = s->width * soff;
1587             if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48LE ||
1588                 s->avctx->pix_fmt == AV_PIX_FMT_RGBA64LE ||
1589                 s->avctx->pix_fmt == AV_PIX_FMT_GRAY16LE ||
1590                 s->avctx->pix_fmt == AV_PIX_FMT_YA16LE ||
1591                 s->avctx->pix_fmt == AV_PIX_FMT_GBRP16LE ||
1592                 s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16LE) {
1593                 for (i = 0; i < decoded_height; i++) {
1594                     for (j = soff; j < ssize; j += 2)
1595                         AV_WL16(dst + j, AV_RL16(dst + j) + AV_RL16(dst + j - soff));
1596                     dst += stride;
1597                 }
1598             } else if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48BE ||
1599                        s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE ||
1600                        s->avctx->pix_fmt == AV_PIX_FMT_GRAY16BE ||
1601                        s->avctx->pix_fmt == AV_PIX_FMT_YA16BE ||
1602                        s->avctx->pix_fmt == AV_PIX_FMT_GBRP16BE ||
1603                        s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16BE) {
1604                 for (i = 0; i < decoded_height; i++) {
1605                     for (j = soff; j < ssize; j += 2)
1606                         AV_WB16(dst + j, AV_RB16(dst + j) + AV_RB16(dst + j - soff));
1607                     dst += stride;
1608                 }
1609             } else {
1610                 for (i = 0; i < decoded_height; i++) {
1611                     for (j = soff; j < ssize; j++)
1612                         dst[j] += dst[j - soff];
1613                     dst += stride;
1614                 }
1615             }
1616         }
1617
1618         if (s->photometric == TIFF_PHOTOMETRIC_WHITE_IS_ZERO) {
1619             int c = (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 ? (1<<s->bpp) - 1 : 255);
1620             dst = p->data[plane];
1621             for (i = 0; i < s->height; i++) {
1622                 for (j = 0; j < stride; j++)
1623                     dst[j] = c - dst[j];
1624                 dst += stride;
1625             }
1626         }
1627
1628         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1629             (s->avctx->pix_fmt == AV_PIX_FMT_RGB0 || s->avctx->pix_fmt == AV_PIX_FMT_RGBA)) {
1630             int x = s->avctx->pix_fmt == AV_PIX_FMT_RGB0 ? 4 : 5;
1631             uint8_t *src = five_planes ? five_planes : p->data[plane];
1632             dst = p->data[plane];
1633             for (i = 0; i < s->height; i++) {
1634                 for (j = 0; j < s->width; j++) {
1635                     int k =  255 - src[x * j + 3];
1636                     int r = (255 - src[x * j    ]) * k;
1637                     int g = (255 - src[x * j + 1]) * k;
1638                     int b = (255 - src[x * j + 2]) * k;
1639                     dst[4 * j    ] = r * 257 >> 16;
1640                     dst[4 * j + 1] = g * 257 >> 16;
1641                     dst[4 * j + 2] = b * 257 >> 16;
1642                     dst[4 * j + 3] = s->avctx->pix_fmt == AV_PIX_FMT_RGBA ? src[x * j + 4] : 255;
1643                 }
1644                 src += stride;
1645                 dst += p->linesize[plane];
1646             }
1647             av_freep(&five_planes);
1648         } else if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1649             s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE) {
1650             dst = p->data[plane];
1651             for (i = 0; i < s->height; i++) {
1652                 for (j = 0; j < s->width; j++) {
1653                     uint64_t k =  65535 - AV_RB16(dst + 8 * j + 6);
1654                     uint64_t r = (65535 - AV_RB16(dst + 8 * j    )) * k;
1655                     uint64_t g = (65535 - AV_RB16(dst + 8 * j + 2)) * k;
1656                     uint64_t b = (65535 - AV_RB16(dst + 8 * j + 4)) * k;
1657                     AV_WB16(dst + 8 * j    , r * 65537 >> 32);
1658                     AV_WB16(dst + 8 * j + 2, g * 65537 >> 32);
1659                     AV_WB16(dst + 8 * j + 4, b * 65537 >> 32);
1660                     AV_WB16(dst + 8 * j + 6, 65535);
1661                 }
1662                 dst += p->linesize[plane];
1663             }
1664         }
1665     }
1666
1667     if (s->planar && s->bppcount > 2) {
1668         FFSWAP(uint8_t*, p->data[0],     p->data[2]);
1669         FFSWAP(int,      p->linesize[0], p->linesize[2]);
1670         FFSWAP(uint8_t*, p->data[0],     p->data[1]);
1671         FFSWAP(int,      p->linesize[0], p->linesize[1]);
1672     }
1673
1674     if (s->is_bayer && s->white_level && s->bpp == 16) {
1675         uint16_t *dst = (uint16_t *)p->data[0];
1676         for (i = 0; i < s->height; i++) {
1677             for (j = 0; j < s->width; j++)
1678                 dst[j] = FFMIN((dst[j] / (float)s->white_level) * 65535, 65535);
1679             dst += stride / 2;
1680         }
1681     }
1682
1683     *got_frame = 1;
1684
1685     return avpkt->size;
1686 }
1687
1688 static av_cold int tiff_init(AVCodecContext *avctx)
1689 {
1690     TiffContext *s = avctx->priv_data;
1691
1692     s->width  = 0;
1693     s->height = 0;
1694     s->subsampling[0] =
1695     s->subsampling[1] = 1;
1696     s->avctx  = avctx;
1697     ff_lzw_decode_open(&s->lzw);
1698     if (!s->lzw)
1699         return AVERROR(ENOMEM);
1700     ff_ccitt_unpack_init();
1701
1702     return 0;
1703 }
1704
1705 static av_cold int tiff_end(AVCodecContext *avctx)
1706 {
1707     TiffContext *const s = avctx->priv_data;
1708
1709     free_geotags(s);
1710
1711     ff_lzw_decode_close(&s->lzw);
1712     av_freep(&s->deinvert_buf);
1713     s->deinvert_buf_size = 0;
1714     av_freep(&s->yuv_line);
1715     s->yuv_line_size = 0;
1716     av_freep(&s->fax_buffer);
1717     s->fax_buffer_size = 0;
1718     return 0;
1719 }
1720
1721 #define OFFSET(x) offsetof(TiffContext, x)
1722 static const AVOption tiff_options[] = {
1723     { "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 },
1724     { "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 },
1725     { "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 },
1726     { NULL },
1727 };
1728
1729 static const AVClass tiff_decoder_class = {
1730     .class_name = "TIFF decoder",
1731     .item_name  = av_default_item_name,
1732     .option     = tiff_options,
1733     .version    = LIBAVUTIL_VERSION_INT,
1734 };
1735
1736 AVCodec ff_tiff_decoder = {
1737     .name           = "tiff",
1738     .long_name      = NULL_IF_CONFIG_SMALL("TIFF image"),
1739     .type           = AVMEDIA_TYPE_VIDEO,
1740     .id             = AV_CODEC_ID_TIFF,
1741     .priv_data_size = sizeof(TiffContext),
1742     .init           = tiff_init,
1743     .close          = tiff_end,
1744     .decode         = decode_frame,
1745     .init_thread_copy = ONLY_IF_THREADS_ENABLED(tiff_init),
1746     .capabilities   = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_FRAME_THREADS,
1747     .priv_class     = &tiff_decoder_class,
1748 };