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