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