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