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