]> git.sesse.net Git - ffmpeg/blob - libavcodec/tiff.c
avcodec/tiff: remove the unnecessary type conversion
[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/error.h"
39 #include "libavutil/intreadwrite.h"
40 #include "libavutil/imgutils.h"
41 #include "libavutil/opt.h"
42 #include "avcodec.h"
43 #include "bytestream.h"
44 #include "faxcompr.h"
45 #include "internal.h"
46 #include "lzw.h"
47 #include "mathops.h"
48 #include "tiff.h"
49 #include "tiff_data.h"
50 #include "mjpegdec.h"
51 #include "thread.h"
52 #include "get_bits.h"
53
54 typedef struct TiffContext {
55     AVClass *class;
56     AVCodecContext *avctx;
57     GetByteContext gb;
58
59     /* JPEG decoding for DNG */
60     AVCodecContext *avctx_mjpeg; // wrapper context for MJPEG
61     AVFrame *jpgframe;           // decoded JPEG tile
62
63     int get_subimage;
64     uint16_t get_page;
65     int get_thumbnail;
66
67     enum TiffType tiff_type;
68     int width, height;
69     unsigned int bpp, bppcount;
70     uint32_t palette[256];
71     int palette_is_set;
72     int le;
73     enum TiffCompr compr;
74     enum TiffPhotometric photometric;
75     int planar;
76     int subsampling[2];
77     int fax_opts;
78     int predictor;
79     int fill_order;
80     uint32_t res[4];
81     int is_thumbnail;
82
83     int is_bayer;
84     uint8_t pattern[4];
85     unsigned black_level;
86     unsigned white_level;
87     uint16_t dng_lut[65536];
88
89     uint32_t sub_ifd;
90     uint16_t cur_page;
91
92     int strips, rps, sstype;
93     int sot;
94     int stripsizesoff, stripsize, stripoff, strippos;
95     LZWState *lzw;
96
97     /* Tile support */
98     int is_tiled;
99     int tile_byte_counts_offset, tile_offsets_offset;
100     int tile_width, tile_length;
101     int tile_count;
102
103     int is_jpeg;
104
105     uint8_t *deinvert_buf;
106     int deinvert_buf_size;
107     uint8_t *yuv_line;
108     unsigned int yuv_line_size;
109     uint8_t *fax_buffer;
110     unsigned int fax_buffer_size;
111
112     int geotag_count;
113     TiffGeoTag *geotags;
114 } TiffContext;
115
116 static void tiff_set_type(TiffContext *s, enum TiffType tiff_type) {
117     if (s->tiff_type < tiff_type) // Prioritize higher-valued entries
118         s->tiff_type = tiff_type;
119 }
120
121 static void free_geotags(TiffContext *const s)
122 {
123     int i;
124     for (i = 0; i < s->geotag_count; i++) {
125         if (s->geotags[i].val)
126             av_freep(&s->geotags[i].val);
127     }
128     av_freep(&s->geotags);
129     s->geotag_count = 0;
130 }
131
132 #define RET_GEOKEY(TYPE, array, element)\
133     if (key >= TIFF_##TYPE##_KEY_ID_OFFSET &&\
134         key - TIFF_##TYPE##_KEY_ID_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_name_type_map))\
135         return ff_tiff_##array##_name_type_map[key - TIFF_##TYPE##_KEY_ID_OFFSET].element;
136
137 static const char *get_geokey_name(int key)
138 {
139     RET_GEOKEY(VERT, vert, name);
140     RET_GEOKEY(PROJ, proj, name);
141     RET_GEOKEY(GEOG, geog, name);
142     RET_GEOKEY(CONF, conf, name);
143
144     return NULL;
145 }
146
147 static int get_geokey_type(int key)
148 {
149     RET_GEOKEY(VERT, vert, type);
150     RET_GEOKEY(PROJ, proj, type);
151     RET_GEOKEY(GEOG, geog, type);
152     RET_GEOKEY(CONF, conf, type);
153
154     return AVERROR_INVALIDDATA;
155 }
156
157 static int cmp_id_key(const void *id, const void *k)
158 {
159     return *(const int*)id - ((const TiffGeoTagKeyName*)k)->key;
160 }
161
162 static const char *search_keyval(const TiffGeoTagKeyName *keys, int n, int id)
163 {
164     TiffGeoTagKeyName *r = bsearch(&id, keys, n, sizeof(keys[0]), cmp_id_key);
165     if(r)
166         return r->name;
167
168     return NULL;
169 }
170
171 static char *get_geokey_val(int key, int val)
172 {
173     char *ap;
174
175     if (val == TIFF_GEO_KEY_UNDEFINED)
176         return av_strdup("undefined");
177     if (val == TIFF_GEO_KEY_USER_DEFINED)
178         return av_strdup("User-Defined");
179
180 #define RET_GEOKEY_VAL(TYPE, array)\
181     if (val >= TIFF_##TYPE##_OFFSET &&\
182         val - TIFF_##TYPE##_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_codes))\
183         return av_strdup(ff_tiff_##array##_codes[val - TIFF_##TYPE##_OFFSET]);
184
185     switch (key) {
186     case TIFF_GT_MODEL_TYPE_GEOKEY:
187         RET_GEOKEY_VAL(GT_MODEL_TYPE, gt_model_type);
188         break;
189     case TIFF_GT_RASTER_TYPE_GEOKEY:
190         RET_GEOKEY_VAL(GT_RASTER_TYPE, gt_raster_type);
191         break;
192     case TIFF_GEOG_LINEAR_UNITS_GEOKEY:
193     case TIFF_PROJ_LINEAR_UNITS_GEOKEY:
194     case TIFF_VERTICAL_UNITS_GEOKEY:
195         RET_GEOKEY_VAL(LINEAR_UNIT, linear_unit);
196         break;
197     case TIFF_GEOG_ANGULAR_UNITS_GEOKEY:
198     case TIFF_GEOG_AZIMUTH_UNITS_GEOKEY:
199         RET_GEOKEY_VAL(ANGULAR_UNIT, angular_unit);
200         break;
201     case TIFF_GEOGRAPHIC_TYPE_GEOKEY:
202         RET_GEOKEY_VAL(GCS_TYPE, gcs_type);
203         RET_GEOKEY_VAL(GCSE_TYPE, gcse_type);
204         break;
205     case TIFF_GEOG_GEODETIC_DATUM_GEOKEY:
206         RET_GEOKEY_VAL(GEODETIC_DATUM, geodetic_datum);
207         RET_GEOKEY_VAL(GEODETIC_DATUM_E, geodetic_datum_e);
208         break;
209     case TIFF_GEOG_ELLIPSOID_GEOKEY:
210         RET_GEOKEY_VAL(ELLIPSOID, ellipsoid);
211         break;
212     case TIFF_GEOG_PRIME_MERIDIAN_GEOKEY:
213         RET_GEOKEY_VAL(PRIME_MERIDIAN, prime_meridian);
214         break;
215     case TIFF_PROJECTED_CS_TYPE_GEOKEY:
216         ap = av_strdup(search_keyval(ff_tiff_proj_cs_type_codes, FF_ARRAY_ELEMS(ff_tiff_proj_cs_type_codes), val));
217         if(ap) return ap;
218         break;
219     case TIFF_PROJECTION_GEOKEY:
220         ap = av_strdup(search_keyval(ff_tiff_projection_codes, FF_ARRAY_ELEMS(ff_tiff_projection_codes), val));
221         if(ap) return ap;
222         break;
223     case TIFF_PROJ_COORD_TRANS_GEOKEY:
224         RET_GEOKEY_VAL(COORD_TRANS, coord_trans);
225         break;
226     case TIFF_VERTICAL_CS_TYPE_GEOKEY:
227         RET_GEOKEY_VAL(VERT_CS, vert_cs);
228         RET_GEOKEY_VAL(ORTHO_VERT_CS, ortho_vert_cs);
229         break;
230
231     }
232
233     ap = av_malloc(14);
234     if (ap)
235         snprintf(ap, 14, "Unknown-%d", val);
236     return ap;
237 }
238
239 static char *doubles2str(double *dp, int count, const char *sep)
240 {
241     int i;
242     char *ap, *ap0;
243     uint64_t component_len;
244     if (!sep) sep = ", ";
245     component_len = 24LL + strlen(sep);
246     if (count >= (INT_MAX - 1)/component_len)
247         return NULL;
248     ap = av_malloc(component_len * count + 1);
249     if (!ap)
250         return NULL;
251     ap0   = ap;
252     ap[0] = '\0';
253     for (i = 0; i < count; i++) {
254         unsigned l = snprintf(ap, component_len, "%.15g%s", dp[i], sep);
255         if(l >= component_len) {
256             av_free(ap0);
257             return NULL;
258         }
259         ap += l;
260     }
261     ap0[strlen(ap0) - strlen(sep)] = '\0';
262     return ap0;
263 }
264
265 static int add_metadata(int count, int type,
266                         const char *name, const char *sep, TiffContext *s, AVFrame *frame)
267 {
268     switch(type) {
269     case TIFF_DOUBLE: return ff_tadd_doubles_metadata(count, name, sep, &s->gb, s->le, &frame->metadata);
270     case TIFF_SHORT : return ff_tadd_shorts_metadata(count, name, sep, &s->gb, s->le, 0, &frame->metadata);
271     case TIFF_STRING: return ff_tadd_string_metadata(count, name, &s->gb, s->le, &frame->metadata);
272     default         : return AVERROR_INVALIDDATA;
273     };
274 }
275
276 static void av_always_inline dng_blit(TiffContext *s, uint8_t *dst, int dst_stride,
277                                       const uint8_t *src, int src_stride, int width, int height,
278                                       int is_single_comp, int is_u16);
279
280 static void av_always_inline horizontal_fill(TiffContext *s,
281                                              unsigned int bpp, uint8_t* dst,
282                                              int usePtr, const uint8_t *src,
283                                              uint8_t c, int width, int offset)
284 {
285     switch (bpp) {
286     case 1:
287         while (--width >= 0) {
288             dst[(width+offset)*8+7] = (usePtr ? src[width] : c)      & 0x1;
289             dst[(width+offset)*8+6] = (usePtr ? src[width] : c) >> 1 & 0x1;
290             dst[(width+offset)*8+5] = (usePtr ? src[width] : c) >> 2 & 0x1;
291             dst[(width+offset)*8+4] = (usePtr ? src[width] : c) >> 3 & 0x1;
292             dst[(width+offset)*8+3] = (usePtr ? src[width] : c) >> 4 & 0x1;
293             dst[(width+offset)*8+2] = (usePtr ? src[width] : c) >> 5 & 0x1;
294             dst[(width+offset)*8+1] = (usePtr ? src[width] : c) >> 6 & 0x1;
295             dst[(width+offset)*8+0] = (usePtr ? src[width] : c) >> 7;
296         }
297         break;
298     case 2:
299         while (--width >= 0) {
300             dst[(width+offset)*4+3] = (usePtr ? src[width] : c) & 0x3;
301             dst[(width+offset)*4+2] = (usePtr ? src[width] : c) >> 2 & 0x3;
302             dst[(width+offset)*4+1] = (usePtr ? src[width] : c) >> 4 & 0x3;
303             dst[(width+offset)*4+0] = (usePtr ? src[width] : c) >> 6;
304         }
305         break;
306     case 4:
307         while (--width >= 0) {
308             dst[(width+offset)*2+1] = (usePtr ? src[width] : c) & 0xF;
309             dst[(width+offset)*2+0] = (usePtr ? src[width] : c) >> 4;
310         }
311         break;
312     case 10:
313     case 12:
314     case 14: {
315             uint16_t *dst16 = (uint16_t *)dst;
316             int is_dng = (s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG);
317             uint8_t shift = is_dng ? 0 : 16 - bpp;
318             GetBitContext gb;
319
320             init_get_bits8(&gb, src, width);
321             for (int i = 0; i < s->width; i++) {
322                 dst16[i] = get_bits(&gb, bpp) << shift;
323             }
324         }
325         break;
326     default:
327         if (usePtr) {
328             memcpy(dst + offset, src, width);
329         } else {
330             memset(dst + offset, c, width);
331         }
332     }
333 }
334
335 static int deinvert_buffer(TiffContext *s, const uint8_t *src, int size)
336 {
337     int i;
338
339     av_fast_padded_malloc(&s->deinvert_buf, &s->deinvert_buf_size, size);
340     if (!s->deinvert_buf)
341         return AVERROR(ENOMEM);
342     for (i = 0; i < size; i++)
343         s->deinvert_buf[i] = ff_reverse[src[i]];
344
345     return 0;
346 }
347
348 static void unpack_gray(TiffContext *s, AVFrame *p,
349                        const uint8_t *src, int lnum, int width, int bpp)
350 {
351     GetBitContext gb;
352     uint16_t *dst = (uint16_t *)(p->data[0] + lnum * p->linesize[0]);
353
354     init_get_bits8(&gb, src, width);
355
356     for (int i = 0; i < s->width; i++) {
357         dst[i] = get_bits(&gb, bpp);
358     }
359 }
360
361 static void unpack_yuv(TiffContext *s, AVFrame *p,
362                        const uint8_t *src, int lnum)
363 {
364     int i, j, k;
365     int w       = (s->width - 1) / s->subsampling[0] + 1;
366     uint8_t *pu = &p->data[1][lnum / s->subsampling[1] * p->linesize[1]];
367     uint8_t *pv = &p->data[2][lnum / s->subsampling[1] * p->linesize[2]];
368     if (s->width % s->subsampling[0] || s->height % s->subsampling[1]) {
369         for (i = 0; i < w; i++) {
370             for (j = 0; j < s->subsampling[1]; j++)
371                 for (k = 0; k < s->subsampling[0]; k++)
372                     p->data[0][FFMIN(lnum + j, s->height-1) * p->linesize[0] +
373                                FFMIN(i * s->subsampling[0] + k, s->width-1)] = *src++;
374             *pu++ = *src++;
375             *pv++ = *src++;
376         }
377     }else{
378         for (i = 0; i < w; i++) {
379             for (j = 0; j < s->subsampling[1]; j++)
380                 for (k = 0; k < s->subsampling[0]; k++)
381                     p->data[0][(lnum + j) * p->linesize[0] +
382                                i * s->subsampling[0] + k] = *src++;
383             *pu++ = *src++;
384             *pv++ = *src++;
385         }
386     }
387 }
388
389 #if CONFIG_ZLIB
390 static int tiff_uncompress(uint8_t *dst, unsigned long *len, const uint8_t *src,
391                            int size)
392 {
393     z_stream zstream = { 0 };
394     int zret;
395
396     zstream.next_in   = src;
397     zstream.avail_in  = size;
398     zstream.next_out  = dst;
399     zstream.avail_out = *len;
400     zret              = inflateInit(&zstream);
401     if (zret != Z_OK) {
402         av_log(NULL, AV_LOG_ERROR, "Inflate init error: %d\n", zret);
403         return zret;
404     }
405     zret = inflate(&zstream, Z_SYNC_FLUSH);
406     inflateEnd(&zstream);
407     *len = zstream.total_out;
408     return zret == Z_STREAM_END ? Z_OK : zret;
409 }
410
411 static int tiff_unpack_zlib(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
412                             const uint8_t *src, int size, int width, int lines,
413                             int strip_start, int is_yuv)
414 {
415     uint8_t *zbuf;
416     unsigned long outlen;
417     int ret, line;
418     outlen = width * lines;
419     zbuf   = av_malloc(outlen);
420     if (!zbuf)
421         return AVERROR(ENOMEM);
422     if (s->fill_order) {
423         if ((ret = deinvert_buffer(s, src, size)) < 0) {
424             av_free(zbuf);
425             return ret;
426         }
427         src = s->deinvert_buf;
428     }
429     ret = tiff_uncompress(zbuf, &outlen, src, size);
430     if (ret != Z_OK) {
431         av_log(s->avctx, AV_LOG_ERROR,
432                "Uncompressing failed (%lu of %lu) with error %d\n", outlen,
433                (unsigned long)width * lines, ret);
434         av_free(zbuf);
435         return AVERROR_UNKNOWN;
436     }
437     src = zbuf;
438     for (line = 0; line < lines; line++) {
439         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
440             horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
441         } else {
442             memcpy(dst, src, width);
443         }
444         if (is_yuv) {
445             unpack_yuv(s, p, dst, strip_start + line);
446             line += s->subsampling[1] - 1;
447         }
448         dst += stride;
449         src += width;
450     }
451     av_free(zbuf);
452     return 0;
453 }
454 #endif
455
456 #if CONFIG_LZMA
457 static int tiff_uncompress_lzma(uint8_t *dst, uint64_t *len, const uint8_t *src,
458                                 int size)
459 {
460     lzma_stream stream = LZMA_STREAM_INIT;
461     lzma_ret ret;
462
463     stream.next_in   = (uint8_t *)src;
464     stream.avail_in  = size;
465     stream.next_out  = dst;
466     stream.avail_out = *len;
467     ret              = lzma_stream_decoder(&stream, UINT64_MAX, 0);
468     if (ret != LZMA_OK) {
469         av_log(NULL, AV_LOG_ERROR, "LZMA init error: %d\n", ret);
470         return ret;
471     }
472     ret = lzma_code(&stream, LZMA_RUN);
473     lzma_end(&stream);
474     *len = stream.total_out;
475     return ret == LZMA_STREAM_END ? LZMA_OK : ret;
476 }
477
478 static int tiff_unpack_lzma(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
479                             const uint8_t *src, int size, int width, int lines,
480                             int strip_start, int is_yuv)
481 {
482     uint64_t outlen = width * (uint64_t)lines;
483     int ret, line;
484     uint8_t *buf = av_malloc(outlen);
485     if (!buf)
486         return AVERROR(ENOMEM);
487     if (s->fill_order) {
488         if ((ret = deinvert_buffer(s, src, size)) < 0) {
489             av_free(buf);
490             return ret;
491         }
492         src = s->deinvert_buf;
493     }
494     ret = tiff_uncompress_lzma(buf, &outlen, src, size);
495     if (ret != LZMA_OK) {
496         av_log(s->avctx, AV_LOG_ERROR,
497                "Uncompressing failed (%"PRIu64" of %"PRIu64") with error %d\n", outlen,
498                (uint64_t)width * lines, ret);
499         av_free(buf);
500         return AVERROR_UNKNOWN;
501     }
502     src = buf;
503     for (line = 0; line < lines; line++) {
504         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
505             horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
506         } else {
507             memcpy(dst, src, width);
508         }
509         if (is_yuv) {
510             unpack_yuv(s, p, dst, strip_start + line);
511             line += s->subsampling[1] - 1;
512         }
513         dst += stride;
514         src += width;
515     }
516     av_free(buf);
517     return 0;
518 }
519 #endif
520
521 static int tiff_unpack_fax(TiffContext *s, uint8_t *dst, int stride,
522                            const uint8_t *src, int size, int width, int lines)
523 {
524     int i, ret = 0;
525     int line;
526     uint8_t *src2;
527
528     av_fast_padded_malloc(&s->fax_buffer, &s->fax_buffer_size, size);
529     src2 = s->fax_buffer;
530
531     if (!src2) {
532         av_log(s->avctx, AV_LOG_ERROR,
533                "Error allocating temporary buffer\n");
534         return AVERROR(ENOMEM);
535     }
536
537     if (!s->fill_order) {
538         memcpy(src2, src, size);
539     } else {
540         for (i = 0; i < size; i++)
541             src2[i] = ff_reverse[src[i]];
542     }
543     memset(src2 + size, 0, AV_INPUT_BUFFER_PADDING_SIZE);
544     ret = ff_ccitt_unpack(s->avctx, src2, size, dst, lines, stride,
545                           s->compr, s->fax_opts);
546     if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
547         for (line = 0; line < lines; line++) {
548             horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
549             dst += stride;
550         }
551     return ret;
552 }
553
554 static int dng_decode_strip(AVCodecContext *avctx, AVFrame *frame);
555
556 static int tiff_unpack_strip(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
557                              const uint8_t *src, int size, int strip_start, int lines)
558 {
559     PutByteContext pb;
560     int c, line, pixels, code, ret;
561     const uint8_t *ssrc = src;
562     int width = ((s->width * s->bpp) + 7) >> 3;
563     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(p->format);
564     int is_yuv = !(desc->flags & AV_PIX_FMT_FLAG_RGB) &&
565                  (desc->flags & AV_PIX_FMT_FLAG_PLANAR) &&
566                  desc->nb_components >= 3;
567     int is_dng;
568
569     if (s->planar)
570         width /= s->bppcount;
571
572     if (size <= 0)
573         return AVERROR_INVALIDDATA;
574
575     if (is_yuv) {
576         int bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
577                             s->subsampling[0] * s->subsampling[1] + 7) >> 3;
578         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
579         if (s->yuv_line == NULL) {
580             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
581             return AVERROR(ENOMEM);
582         }
583         dst = s->yuv_line;
584         stride = 0;
585
586         width = (s->width - 1) / s->subsampling[0] + 1;
587         width = width * s->subsampling[0] * s->subsampling[1] + 2*width;
588         av_assert0(width <= bytes_per_row);
589         av_assert0(s->bpp == 24);
590     }
591     if (s->is_bayer) {
592         av_assert0(width == (s->bpp * s->width + 7) >> 3);
593     }
594     if (p->format == AV_PIX_FMT_GRAY12) {
595         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, width);
596         if (s->yuv_line == NULL) {
597             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
598             return AVERROR(ENOMEM);
599         }
600         dst = s->yuv_line;
601         stride = 0;
602     }
603
604     if (s->compr == TIFF_DEFLATE || s->compr == TIFF_ADOBE_DEFLATE) {
605 #if CONFIG_ZLIB
606         return tiff_unpack_zlib(s, p, dst, stride, src, size, width, lines,
607                                 strip_start, is_yuv);
608 #else
609         av_log(s->avctx, AV_LOG_ERROR,
610                "zlib support not enabled, "
611                "deflate compression not supported\n");
612         return AVERROR(ENOSYS);
613 #endif
614     }
615     if (s->compr == TIFF_LZMA) {
616 #if CONFIG_LZMA
617         return tiff_unpack_lzma(s, p, dst, stride, src, size, width, lines,
618                                 strip_start, is_yuv);
619 #else
620         av_log(s->avctx, AV_LOG_ERROR,
621                "LZMA support not enabled\n");
622         return AVERROR(ENOSYS);
623 #endif
624     }
625     if (s->compr == TIFF_LZW) {
626         if (s->fill_order) {
627             if ((ret = deinvert_buffer(s, src, size)) < 0)
628                 return ret;
629             ssrc = src = s->deinvert_buf;
630         }
631         if (size > 1 && !src[0] && (src[1]&1)) {
632             av_log(s->avctx, AV_LOG_ERROR, "Old style LZW is unsupported\n");
633         }
634         if ((ret = ff_lzw_decode_init(s->lzw, 8, src, size, FF_LZW_TIFF)) < 0) {
635             av_log(s->avctx, AV_LOG_ERROR, "Error initializing LZW decoder\n");
636             return ret;
637         }
638         for (line = 0; line < lines; line++) {
639             pixels = ff_lzw_decode(s->lzw, dst, width);
640             if (pixels < width) {
641                 av_log(s->avctx, AV_LOG_ERROR, "Decoded only %i bytes of %i\n",
642                        pixels, width);
643                 return AVERROR_INVALIDDATA;
644             }
645             if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
646                 horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
647             if (is_yuv) {
648                 unpack_yuv(s, p, dst, strip_start + line);
649                 line += s->subsampling[1] - 1;
650             } else if (p->format == AV_PIX_FMT_GRAY12) {
651                 unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
652             }
653             dst += stride;
654         }
655         return 0;
656     }
657     if (s->compr == TIFF_CCITT_RLE ||
658         s->compr == TIFF_G3        ||
659         s->compr == TIFF_G4) {
660         if (is_yuv || p->format == AV_PIX_FMT_GRAY12)
661             return AVERROR_INVALIDDATA;
662
663         return tiff_unpack_fax(s, dst, stride, src, size, width, lines);
664     }
665
666     bytestream2_init(&s->gb, src, size);
667     bytestream2_init_writer(&pb, dst, is_yuv ? s->yuv_line_size : (stride * lines));
668
669     is_dng = (s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG);
670
671     /* Decode JPEG-encoded DNGs with strips */
672     if (s->compr == TIFF_NEWJPEG && is_dng) {
673         if (s->strips > 1) {
674             av_log(s->avctx, AV_LOG_ERROR, "More than one DNG JPEG strips unsupported\n");
675             return AVERROR_PATCHWELCOME;
676         }
677         if ((ret = dng_decode_strip(s->avctx, p)) < 0)
678             return ret;
679         return 0;
680     }
681
682     for (line = 0; line < lines; line++) {
683         if (src - ssrc > size) {
684             av_log(s->avctx, AV_LOG_ERROR, "Source data overread\n");
685             return AVERROR_INVALIDDATA;
686         }
687
688         if (bytestream2_get_bytes_left(&s->gb) == 0 || bytestream2_get_eof(&pb))
689             break;
690         bytestream2_seek_p(&pb, stride * line, SEEK_SET);
691         switch (s->compr) {
692         case TIFF_RAW:
693             if (ssrc + size - src < width)
694                 return AVERROR_INVALIDDATA;
695
696             if (!s->fill_order) {
697                 horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 || s->is_bayer),
698                                 dst, 1, src, 0, width, 0);
699             } else {
700                 int i;
701                 for (i = 0; i < width; i++)
702                     dst[i] = ff_reverse[src[i]];
703             }
704
705             /* Color processing for DNG images with uncompressed strips (non-tiled) */
706             if (is_dng) {
707                 int is_u16, pixel_size_bytes, pixel_size_bits, elements;
708
709                 is_u16 = (s->bpp > 8);
710                 pixel_size_bits = (is_u16 ? 16 : 8);
711                 pixel_size_bytes = (is_u16 ? sizeof(uint16_t) : sizeof(uint8_t));
712
713                 elements = width / pixel_size_bytes * pixel_size_bits / s->bpp * s->bppcount; // need to account for [1, 16] bpp
714                 av_assert0 (elements * pixel_size_bytes <= FFABS(stride));
715                 dng_blit(s,
716                          dst,
717                          0, // no stride, only 1 line
718                          dst,
719                          0, // no stride, only 1 line
720                          elements,
721                          1,
722                          0, // single-component variation is only preset in JPEG-encoded DNGs
723                          is_u16);
724             }
725
726             src += width;
727             break;
728         case TIFF_PACKBITS:
729             for (pixels = 0; pixels < width;) {
730                 if (ssrc + size - src < 2) {
731                     av_log(s->avctx, AV_LOG_ERROR, "Read went out of bounds\n");
732                     return AVERROR_INVALIDDATA;
733                 }
734                 code = s->fill_order ? (int8_t) ff_reverse[*src++]: (int8_t) *src++;
735                 if (code >= 0) {
736                     code++;
737                     if (pixels + code > width ||
738                         ssrc + size - src < code) {
739                         av_log(s->avctx, AV_LOG_ERROR,
740                                "Copy went out of bounds\n");
741                         return AVERROR_INVALIDDATA;
742                     }
743                     horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
744                                     dst, 1, src, 0, code, pixels);
745                     src    += code;
746                     pixels += code;
747                 } else if (code != -128) { // -127..-1
748                     code = (-code) + 1;
749                     if (pixels + code > width) {
750                         av_log(s->avctx, AV_LOG_ERROR,
751                                "Run went out of bounds\n");
752                         return AVERROR_INVALIDDATA;
753                     }
754                     c = *src++;
755                     horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
756                                     dst, 0, NULL, c, code, pixels);
757                     pixels += code;
758                 }
759             }
760             if (s->fill_order) {
761                 int i;
762                 for (i = 0; i < width; i++)
763                     dst[i] = ff_reverse[dst[i]];
764             }
765             break;
766         }
767         if (is_yuv) {
768             unpack_yuv(s, p, dst, strip_start + line);
769             line += s->subsampling[1] - 1;
770         } else if (p->format == AV_PIX_FMT_GRAY12) {
771             unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
772         }
773         dst += stride;
774     }
775     return 0;
776 }
777
778 /**
779  * Map stored raw sensor values into linear reference values (see: DNG Specification - Chapter 5)
780  */
781 static uint16_t av_always_inline dng_process_color16(uint16_t value,
782                                                      const uint16_t *lut,
783                                                      uint16_t black_level,
784                                                      float scale_factor) {
785     float value_norm;
786
787     // Lookup table lookup
788     if (lut)
789         value = lut[value];
790
791     // Black level subtraction
792     value = av_clip_uint16_c((unsigned)value - black_level);
793
794     // Color scaling
795     value_norm = (float)value * scale_factor;
796
797     value = av_clip_uint16_c(value_norm * 65535);
798
799     return value;
800 }
801
802 static uint16_t av_always_inline dng_process_color8(uint16_t value,
803                                                     const uint16_t *lut,
804                                                     uint16_t black_level,
805                                                     float scale_factor) {
806     return dng_process_color16(value, lut, black_level, scale_factor) >> 8;
807 }
808
809 static void dng_blit(TiffContext *s, uint8_t *dst, int dst_stride,
810                      const uint8_t *src, int src_stride,
811                      int width, int height, int is_single_comp, int is_u16)
812 {
813     int line, col;
814     float scale_factor;
815
816     scale_factor = 1.0f / (s->white_level - s->black_level);
817
818     if (is_single_comp) {
819         if (!is_u16)
820             return; /* <= 8bpp unsupported */
821
822         /* Image is double the width and half the height we need, each row comprises 2 rows of the output
823            (split vertically in the middle). */
824         for (line = 0; line < height / 2; line++) {
825             uint16_t *dst_u16 = (uint16_t *)dst;
826             uint16_t *src_u16 = (uint16_t *)src;
827
828             /* Blit first half of input row row to initial row of output */
829             for (col = 0; col < width; col++)
830                 *dst_u16++ = dng_process_color16(*src_u16++, s->dng_lut, s->black_level, scale_factor);
831
832             /* Advance the destination pointer by a row (source pointer remains in the same place) */
833             dst += dst_stride * sizeof(uint16_t);
834             dst_u16 = (uint16_t *)dst;
835
836             /* Blit second half of input row row to next row of output */
837             for (col = 0; col < width; col++)
838                 *dst_u16++ = dng_process_color16(*src_u16++, s->dng_lut, s->black_level, scale_factor);
839
840             dst += dst_stride * sizeof(uint16_t);
841             src += src_stride * sizeof(uint16_t);
842         }
843     } else {
844         /* Input and output image are the same size and the MJpeg decoder has done per-component
845            deinterleaving, so blitting here is straightforward. */
846         if (is_u16) {
847             for (line = 0; line < height; line++) {
848                 uint16_t *dst_u16 = (uint16_t *)dst;
849                 uint16_t *src_u16 = (uint16_t *)src;
850
851                 for (col = 0; col < width; col++)
852                     *dst_u16++ = dng_process_color16(*src_u16++, s->dng_lut, s->black_level, scale_factor);
853
854                 dst += dst_stride * sizeof(uint16_t);
855                 src += src_stride * sizeof(uint16_t);
856             }
857         } else {
858             for (line = 0; line < height; line++) {
859                 for (col = 0; col < width; col++)
860                     *dst++ = dng_process_color8(*src++, s->dng_lut, s->black_level, scale_factor);
861
862                 dst += dst_stride;
863                 src += src_stride;
864             }
865         }
866     }
867 }
868
869 static int dng_decode_jpeg(AVCodecContext *avctx, AVFrame *frame,
870                            int tile_byte_count, int dst_x, int dst_y, int w, int h)
871 {
872     TiffContext *s = avctx->priv_data;
873     AVPacket jpkt;
874     uint8_t *dst_data, *src_data;
875     uint32_t dst_offset; /* offset from dst buffer in pixels */
876     int is_single_comp, is_u16, pixel_size;
877     int ret;
878
879     /* Prepare a packet and send to the MJPEG decoder */
880     av_init_packet(&jpkt);
881     jpkt.data = (uint8_t*)s->gb.buffer;
882     jpkt.size = tile_byte_count;
883
884     if (s->is_bayer) {
885         MJpegDecodeContext *mjpegdecctx = s->avctx_mjpeg->priv_data;
886         /* We have to set this information here, there is no way to know if a given JPEG is a DNG-embedded
887            image or not from its own data (and we need that information when decoding it). */
888         mjpegdecctx->bayer = 1;
889     }
890
891     ret = avcodec_send_packet(s->avctx_mjpeg, &jpkt);
892     if (ret < 0) {
893         av_log(avctx, AV_LOG_ERROR, "Error submitting a packet for decoding\n");
894         return ret;
895     }
896
897     ret = avcodec_receive_frame(s->avctx_mjpeg, s->jpgframe);
898     if (ret < 0) {
899         av_log(avctx, AV_LOG_ERROR, "JPEG decoding error: %s.\n", av_err2str(ret));
900
901         /* Normally skip, error if explode */
902         if (avctx->err_recognition & AV_EF_EXPLODE)
903             return AVERROR_INVALIDDATA;
904         else
905             return 0;
906     }
907
908     /* Copy the outputted tile's pixels from 'jpgframe' to 'frame' (final buffer) */
909
910     /* See dng_blit for explanation */
911     is_single_comp = (s->avctx_mjpeg->width == w * 2 && s->avctx_mjpeg->height == h / 2);
912
913     is_u16 = (s->bpp > 8);
914     pixel_size = (is_u16 ? sizeof(uint16_t) : sizeof(uint8_t));
915
916     if (is_single_comp && !is_u16) {
917         av_log(s->avctx, AV_LOG_ERROR, "DNGs with bpp <= 8 and 1 component are unsupported\n");
918         av_frame_unref(s->jpgframe);
919         return AVERROR_PATCHWELCOME;
920     }
921
922     dst_offset = dst_x + frame->linesize[0] * dst_y / pixel_size;
923     dst_data = frame->data[0] + dst_offset * pixel_size;
924     src_data = s->jpgframe->data[0];
925
926     dng_blit(s,
927              dst_data,
928              frame->linesize[0] / pixel_size,
929              src_data,
930              s->jpgframe->linesize[0] / pixel_size,
931              w,
932              h,
933              is_single_comp,
934              is_u16);
935
936     av_frame_unref(s->jpgframe);
937
938     return 0;
939 }
940
941 static int dng_decode_tiles(AVCodecContext *avctx, AVFrame *frame, AVPacket *avpkt)
942 {
943     TiffContext *s = avctx->priv_data;
944     int tile_idx;
945     int tile_offset_offset, tile_offset;
946     int tile_byte_count_offset, tile_byte_count;
947     int tile_count_x, tile_count_y;
948     int tile_width, tile_length;
949     int has_width_leftover, has_height_leftover;
950     int tile_x = 0, tile_y = 0;
951     int pos_x = 0, pos_y = 0;
952     int ret;
953
954     s->jpgframe->width  = s->tile_width;
955     s->jpgframe->height = s->tile_length;
956
957     s->avctx_mjpeg->width = s->tile_width;
958     s->avctx_mjpeg->height = s->tile_length;
959
960     has_width_leftover = (s->width % s->tile_width != 0);
961     has_height_leftover = (s->height % s->tile_length != 0);
962
963     /* Calculate tile counts (round up) */
964     tile_count_x = (s->width + s->tile_width - 1) / s->tile_width;
965     tile_count_y = (s->height + s->tile_length - 1) / s->tile_length;
966
967     /* Iterate over the number of tiles */
968     for (tile_idx = 0; tile_idx < s->tile_count; tile_idx++) {
969         tile_x = tile_idx % tile_count_x;
970         tile_y = tile_idx / tile_count_x;
971
972         if (has_width_leftover && tile_x == tile_count_x - 1) // If on the right-most tile
973             tile_width = s->width % s->tile_width;
974         else
975             tile_width = s->tile_width;
976
977         if (has_height_leftover && tile_y == tile_count_y - 1) // If on the bottom-most tile
978             tile_length = s->height % s->tile_length;
979         else
980             tile_length = s->tile_length;
981
982         /* Read tile offset */
983         tile_offset_offset = s->tile_offsets_offset + tile_idx * sizeof(int);
984         bytestream2_seek(&s->gb, tile_offset_offset, SEEK_SET);
985         tile_offset = ff_tget_long(&s->gb, s->le);
986
987         /* Read tile byte size */
988         tile_byte_count_offset = s->tile_byte_counts_offset + tile_idx * sizeof(int);
989         bytestream2_seek(&s->gb, tile_byte_count_offset, SEEK_SET);
990         tile_byte_count = ff_tget_long(&s->gb, s->le);
991
992         /* Seek to tile data */
993         bytestream2_seek(&s->gb, tile_offset, SEEK_SET);
994
995         /* Decode JPEG tile and copy it in the reference frame */
996         ret = dng_decode_jpeg(avctx, frame, tile_byte_count, pos_x, pos_y, tile_width, tile_length);
997
998         if (ret < 0)
999             return ret;
1000
1001         /* Advance current positions */
1002         pos_x += tile_width;
1003         if (tile_x == tile_count_x - 1) { // If on the right edge
1004             pos_x = 0;
1005             pos_y += tile_length;
1006         }
1007     }
1008
1009     /* Frame is ready to be output */
1010     frame->pict_type = AV_PICTURE_TYPE_I;
1011     frame->key_frame = 1;
1012
1013     return avpkt->size;
1014 }
1015
1016 static int dng_decode_strip(AVCodecContext *avctx, AVFrame *frame)
1017 {
1018     TiffContext *s = avctx->priv_data;
1019
1020     s->jpgframe->width  = s->width;
1021     s->jpgframe->height = s->height;
1022
1023     s->avctx_mjpeg->width = s->width;
1024     s->avctx_mjpeg->height = s->height;
1025
1026     return dng_decode_jpeg(avctx, frame, s->stripsize, 0, 0, s->width, s->height);
1027 }
1028
1029 static int init_image(TiffContext *s, ThreadFrame *frame)
1030 {
1031     int ret;
1032     int create_gray_palette = 0;
1033
1034     // make sure there is no aliasing in the following switch
1035     if (s->bpp >= 100 || s->bppcount >= 10) {
1036         av_log(s->avctx, AV_LOG_ERROR,
1037                "Unsupported image parameters: bpp=%d, bppcount=%d\n",
1038                s->bpp, s->bppcount);
1039         return AVERROR_INVALIDDATA;
1040     }
1041
1042     switch (s->planar * 1000 + s->bpp * 10 + s->bppcount + s->is_bayer * 10000) {
1043     case 11:
1044         if (!s->palette_is_set) {
1045             s->avctx->pix_fmt = AV_PIX_FMT_MONOBLACK;
1046             break;
1047         }
1048     case 21:
1049     case 41:
1050         s->avctx->pix_fmt = AV_PIX_FMT_PAL8;
1051         if (!s->palette_is_set) {
1052             create_gray_palette = 1;
1053         }
1054         break;
1055     case 81:
1056         s->avctx->pix_fmt = s->palette_is_set ? AV_PIX_FMT_PAL8 : AV_PIX_FMT_GRAY8;
1057         break;
1058     case 121:
1059         s->avctx->pix_fmt = AV_PIX_FMT_GRAY12;
1060         break;
1061     case 10081:
1062         switch (AV_RL32(s->pattern)) {
1063         case 0x02010100:
1064             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_RGGB8;
1065             break;
1066         case 0x00010102:
1067             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_BGGR8;
1068             break;
1069         case 0x01000201:
1070             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GBRG8;
1071             break;
1072         case 0x01020001:
1073             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GRBG8;
1074             break;
1075         default:
1076             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1077                    AV_RL32(s->pattern));
1078             return AVERROR_PATCHWELCOME;
1079         }
1080         break;
1081     case 10101:
1082     case 10121:
1083     case 10141:
1084     case 10161:
1085         switch (AV_RL32(s->pattern)) {
1086         case 0x02010100:
1087             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_RGGB16;
1088             break;
1089         case 0x00010102:
1090             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_BGGR16;
1091             break;
1092         case 0x01000201:
1093             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GBRG16;
1094             break;
1095         case 0x01020001:
1096             s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GRBG16;
1097             break;
1098         default:
1099             av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1100                    AV_RL32(s->pattern));
1101             return AVERROR_PATCHWELCOME;
1102         }
1103         break;
1104     case 243:
1105         if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1106             if (s->subsampling[0] == 1 && s->subsampling[1] == 1) {
1107                 s->avctx->pix_fmt = AV_PIX_FMT_YUV444P;
1108             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 1) {
1109                 s->avctx->pix_fmt = AV_PIX_FMT_YUV422P;
1110             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 1) {
1111                 s->avctx->pix_fmt = AV_PIX_FMT_YUV411P;
1112             } else if (s->subsampling[0] == 1 && s->subsampling[1] == 2) {
1113                 s->avctx->pix_fmt = AV_PIX_FMT_YUV440P;
1114             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 2) {
1115                 s->avctx->pix_fmt = AV_PIX_FMT_YUV420P;
1116             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 4) {
1117                 s->avctx->pix_fmt = AV_PIX_FMT_YUV410P;
1118             } else {
1119                 av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr subsampling\n");
1120                 return AVERROR_PATCHWELCOME;
1121             }
1122         } else
1123             s->avctx->pix_fmt = AV_PIX_FMT_RGB24;
1124         break;
1125     case 161:
1126         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GRAY16LE : AV_PIX_FMT_GRAY16BE;
1127         break;
1128     case 162:
1129         s->avctx->pix_fmt = AV_PIX_FMT_YA8;
1130         break;
1131     case 322:
1132         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_YA16LE : AV_PIX_FMT_YA16BE;
1133         break;
1134     case 324:
1135         s->avctx->pix_fmt = s->photometric == TIFF_PHOTOMETRIC_SEPARATED ? AV_PIX_FMT_RGB0 : AV_PIX_FMT_RGBA;
1136         break;
1137     case 405:
1138         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED)
1139             s->avctx->pix_fmt = AV_PIX_FMT_RGBA;
1140         else {
1141             av_log(s->avctx, AV_LOG_ERROR,
1142                 "bpp=40 without PHOTOMETRIC_SEPARATED is unsupported\n");
1143             return AVERROR_PATCHWELCOME;
1144         }
1145         break;
1146     case 483:
1147         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGB48LE  : AV_PIX_FMT_RGB48BE;
1148         break;
1149     case 644:
1150         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGBA64LE  : AV_PIX_FMT_RGBA64BE;
1151         break;
1152     case 1243:
1153         s->avctx->pix_fmt = AV_PIX_FMT_GBRP;
1154         break;
1155     case 1324:
1156         s->avctx->pix_fmt = AV_PIX_FMT_GBRAP;
1157         break;
1158     case 1483:
1159         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRP16LE : AV_PIX_FMT_GBRP16BE;
1160         break;
1161     case 1644:
1162         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRAP16LE : AV_PIX_FMT_GBRAP16BE;
1163         break;
1164     default:
1165         av_log(s->avctx, AV_LOG_ERROR,
1166                "This format is not supported (bpp=%d, bppcount=%d)\n",
1167                s->bpp, s->bppcount);
1168         return AVERROR_INVALIDDATA;
1169     }
1170
1171     if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1172         const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(s->avctx->pix_fmt);
1173         if((desc->flags & AV_PIX_FMT_FLAG_RGB) ||
1174            !(desc->flags & AV_PIX_FMT_FLAG_PLANAR) ||
1175            desc->nb_components < 3) {
1176             av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr variant\n");
1177             return AVERROR_INVALIDDATA;
1178         }
1179     }
1180
1181     if (s->width != s->avctx->width || s->height != s->avctx->height) {
1182         ret = ff_set_dimensions(s->avctx, s->width, s->height);
1183         if (ret < 0)
1184             return ret;
1185     }
1186     if ((ret = ff_thread_get_buffer(s->avctx, frame, 0)) < 0)
1187         return ret;
1188     if (s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
1189         if (!create_gray_palette)
1190             memcpy(frame->f->data[1], s->palette, sizeof(s->palette));
1191         else {
1192             /* make default grayscale pal */
1193             int i;
1194             uint32_t *pal = (uint32_t *)frame->f->data[1];
1195             for (i = 0; i < 1<<s->bpp; i++)
1196                 pal[i] = 0xFFU << 24 | i * 255 / ((1<<s->bpp) - 1) * 0x010101;
1197         }
1198     }
1199     return 0;
1200 }
1201
1202 static void set_sar(TiffContext *s, unsigned tag, unsigned num, unsigned den)
1203 {
1204     int offset = tag == TIFF_YRES ? 2 : 0;
1205     s->res[offset++] = num;
1206     s->res[offset]   = den;
1207     if (s->res[0] && s->res[1] && s->res[2] && s->res[3]) {
1208         uint64_t num = s->res[2] * (uint64_t)s->res[1];
1209         uint64_t den = s->res[0] * (uint64_t)s->res[3];
1210         if (num > INT64_MAX || den > INT64_MAX) {
1211             num = num >> 1;
1212             den = den >> 1;
1213         }
1214         av_reduce(&s->avctx->sample_aspect_ratio.num, &s->avctx->sample_aspect_ratio.den,
1215                   num, den, INT32_MAX);
1216         if (!s->avctx->sample_aspect_ratio.den)
1217             s->avctx->sample_aspect_ratio = (AVRational) {0, 1};
1218     }
1219 }
1220
1221 static int tiff_decode_tag(TiffContext *s, AVFrame *frame)
1222 {
1223     AVFrameSideData *sd;
1224     GetByteContext gb_temp;
1225     unsigned tag, type, count, off, value = 0, value2 = 1; // value2 is a denominator so init. to 1
1226     int i, start;
1227     int pos;
1228     int ret;
1229     double *dp;
1230
1231     ret = ff_tread_tag(&s->gb, s->le, &tag, &type, &count, &start);
1232     if (ret < 0) {
1233         goto end;
1234     }
1235
1236     off = bytestream2_tell(&s->gb);
1237     if (count == 1) {
1238         switch (type) {
1239         case TIFF_BYTE:
1240         case TIFF_SHORT:
1241         case TIFF_LONG:
1242             value = ff_tget(&s->gb, type, s->le);
1243             break;
1244         case TIFF_RATIONAL:
1245             value  = ff_tget(&s->gb, TIFF_LONG, s->le);
1246             value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
1247             if (!value2) {
1248                 av_log(s->avctx, AV_LOG_ERROR, "Invalid denominator in rational\n");
1249                 return AVERROR_INVALIDDATA;
1250             }
1251
1252             break;
1253         case TIFF_STRING:
1254             if (count <= 4) {
1255                 break;
1256             }
1257         default:
1258             value = UINT_MAX;
1259         }
1260     }
1261
1262     switch (tag) {
1263     case TIFF_SUBFILE:
1264         s->is_thumbnail = (value != 0);
1265         break;
1266     case TIFF_WIDTH:
1267         s->width = value;
1268         break;
1269     case TIFF_HEIGHT:
1270         s->height = value;
1271         break;
1272     case TIFF_BPP:
1273         if (count > 5U) {
1274             av_log(s->avctx, AV_LOG_ERROR,
1275                    "This format is not supported (bpp=%d, %d components)\n",
1276                    value, count);
1277             return AVERROR_INVALIDDATA;
1278         }
1279         s->bppcount = count;
1280         if (count == 1)
1281             s->bpp = value;
1282         else {
1283             switch (type) {
1284             case TIFF_BYTE:
1285             case TIFF_SHORT:
1286             case TIFF_LONG:
1287                 s->bpp = 0;
1288                 if (bytestream2_get_bytes_left(&s->gb) < type_sizes[type] * count)
1289                     return AVERROR_INVALIDDATA;
1290                 for (i = 0; i < count; i++)
1291                     s->bpp += ff_tget(&s->gb, type, s->le);
1292                 break;
1293             default:
1294                 s->bpp = -1;
1295             }
1296         }
1297         break;
1298     case TIFF_SAMPLES_PER_PIXEL:
1299         if (count != 1) {
1300             av_log(s->avctx, AV_LOG_ERROR,
1301                    "Samples per pixel requires a single value, many provided\n");
1302             return AVERROR_INVALIDDATA;
1303         }
1304         if (value > 5U) {
1305             av_log(s->avctx, AV_LOG_ERROR,
1306                    "Samples per pixel %d is too large\n", value);
1307             return AVERROR_INVALIDDATA;
1308         }
1309         if (s->bppcount == 1)
1310             s->bpp *= value;
1311         s->bppcount = value;
1312         break;
1313     case TIFF_COMPR:
1314         s->compr     = value;
1315         av_log(s->avctx, AV_LOG_DEBUG, "compression: %d\n", s->compr);
1316         s->predictor = 0;
1317         switch (s->compr) {
1318         case TIFF_RAW:
1319         case TIFF_PACKBITS:
1320         case TIFF_LZW:
1321         case TIFF_CCITT_RLE:
1322             break;
1323         case TIFF_G3:
1324         case TIFF_G4:
1325             s->fax_opts = 0;
1326             break;
1327         case TIFF_DEFLATE:
1328         case TIFF_ADOBE_DEFLATE:
1329 #if CONFIG_ZLIB
1330             break;
1331 #else
1332             av_log(s->avctx, AV_LOG_ERROR, "Deflate: ZLib not compiled in\n");
1333             return AVERROR(ENOSYS);
1334 #endif
1335         case TIFF_JPEG:
1336         case TIFF_NEWJPEG:
1337             s->is_jpeg = 1;
1338             break;
1339         case TIFF_LZMA:
1340 #if CONFIG_LZMA
1341             break;
1342 #else
1343             av_log(s->avctx, AV_LOG_ERROR, "LZMA not compiled in\n");
1344             return AVERROR(ENOSYS);
1345 #endif
1346         default:
1347             av_log(s->avctx, AV_LOG_ERROR, "Unknown compression method %i\n",
1348                    s->compr);
1349             return AVERROR_INVALIDDATA;
1350         }
1351         break;
1352     case TIFF_ROWSPERSTRIP:
1353         if (!value || (type == TIFF_LONG && value == UINT_MAX))
1354             value = s->height;
1355         s->rps = FFMIN(value, s->height);
1356         break;
1357     case TIFF_STRIP_OFFS:
1358         if (count == 1) {
1359             if (value > INT_MAX) {
1360                 av_log(s->avctx, AV_LOG_ERROR,
1361                     "strippos %u too large\n", value);
1362                 return AVERROR_INVALIDDATA;
1363             }
1364             s->strippos = 0;
1365             s->stripoff = value;
1366         } else
1367             s->strippos = off;
1368         s->strips = count;
1369         if (s->strips == 1)
1370             s->rps = s->height;
1371         s->sot = type;
1372         break;
1373     case TIFF_STRIP_SIZE:
1374         if (count == 1) {
1375             if (value > INT_MAX) {
1376                 av_log(s->avctx, AV_LOG_ERROR,
1377                     "stripsize %u too large\n", value);
1378                 return AVERROR_INVALIDDATA;
1379             }
1380             s->stripsizesoff = 0;
1381             s->stripsize     = value;
1382             s->strips        = 1;
1383         } else {
1384             s->stripsizesoff = off;
1385         }
1386         s->strips = count;
1387         s->sstype = type;
1388         break;
1389     case TIFF_XRES:
1390     case TIFF_YRES:
1391         set_sar(s, tag, value, value2);
1392         break;
1393     case TIFF_TILE_OFFSETS:
1394         s->tile_offsets_offset = off;
1395         s->tile_count = count;
1396         s->is_tiled = 1;
1397         break;
1398     case TIFF_TILE_BYTE_COUNTS:
1399         s->tile_byte_counts_offset = off;
1400         break;
1401     case TIFF_TILE_LENGTH:
1402         s->tile_length = value;
1403         break;
1404     case TIFF_TILE_WIDTH:
1405         s->tile_width = value;
1406         break;
1407     case TIFF_PREDICTOR:
1408         s->predictor = value;
1409         break;
1410     case TIFF_SUB_IFDS:
1411         if (count == 1)
1412             s->sub_ifd = value;
1413         else if (count > 1)
1414             s->sub_ifd = ff_tget(&s->gb, TIFF_LONG, s->le); /** Only get the first SubIFD */
1415         break;
1416     case DNG_LINEARIZATION_TABLE:
1417         for (int i = 0; i < FFMIN(count, 1 << s->bpp); i++)
1418             s->dng_lut[i] = ff_tget(&s->gb, type, s->le);
1419         break;
1420     case DNG_BLACK_LEVEL:
1421         if (count > 1) {    /* Use the first value in the pattern (assume they're all the same) */
1422             if (type == TIFF_RATIONAL) {
1423                 value  = ff_tget(&s->gb, TIFF_LONG, s->le);
1424                 value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
1425                 if (!value2) {
1426                     av_log(s->avctx, AV_LOG_ERROR, "Invalid black level denominator\n");
1427                     return AVERROR_INVALIDDATA;
1428                 }
1429
1430                 s->black_level = value / value2;
1431             } else
1432                 s->black_level = ff_tget(&s->gb, type, s->le);
1433             av_log(s->avctx, AV_LOG_WARNING, "Assuming black level pattern values are identical\n");
1434         } else {
1435             s->black_level = value / value2;
1436         }
1437         break;
1438     case DNG_WHITE_LEVEL:
1439         s->white_level = value;
1440         break;
1441     case TIFF_CFA_PATTERN_DIM:
1442         if (count != 2 || (ff_tget(&s->gb, type, s->le) != 2 &&
1443                            ff_tget(&s->gb, type, s->le) != 2)) {
1444             av_log(s->avctx, AV_LOG_ERROR, "CFA Pattern dimensions are not 2x2\n");
1445             return AVERROR_INVALIDDATA;
1446         }
1447         break;
1448     case TIFF_CFA_PATTERN:
1449         s->is_bayer = 1;
1450         s->pattern[0] = ff_tget(&s->gb, type, s->le);
1451         s->pattern[1] = ff_tget(&s->gb, type, s->le);
1452         s->pattern[2] = ff_tget(&s->gb, type, s->le);
1453         s->pattern[3] = ff_tget(&s->gb, type, s->le);
1454         break;
1455     case TIFF_PHOTOMETRIC:
1456         switch (value) {
1457         case TIFF_PHOTOMETRIC_WHITE_IS_ZERO:
1458         case TIFF_PHOTOMETRIC_BLACK_IS_ZERO:
1459         case TIFF_PHOTOMETRIC_RGB:
1460         case TIFF_PHOTOMETRIC_PALETTE:
1461         case TIFF_PHOTOMETRIC_SEPARATED:
1462         case TIFF_PHOTOMETRIC_YCBCR:
1463         case TIFF_PHOTOMETRIC_CFA:
1464         case TIFF_PHOTOMETRIC_LINEAR_RAW: // Used by DNG images
1465             s->photometric = value;
1466             break;
1467         case TIFF_PHOTOMETRIC_ALPHA_MASK:
1468         case TIFF_PHOTOMETRIC_CIE_LAB:
1469         case TIFF_PHOTOMETRIC_ICC_LAB:
1470         case TIFF_PHOTOMETRIC_ITU_LAB:
1471         case TIFF_PHOTOMETRIC_LOG_L:
1472         case TIFF_PHOTOMETRIC_LOG_LUV:
1473             avpriv_report_missing_feature(s->avctx,
1474                                           "PhotometricInterpretation 0x%04X",
1475                                           value);
1476             return AVERROR_PATCHWELCOME;
1477         default:
1478             av_log(s->avctx, AV_LOG_ERROR, "PhotometricInterpretation %u is "
1479                    "unknown\n", value);
1480             return AVERROR_INVALIDDATA;
1481         }
1482         break;
1483     case TIFF_FILL_ORDER:
1484         if (value < 1 || value > 2) {
1485             av_log(s->avctx, AV_LOG_ERROR,
1486                    "Unknown FillOrder value %d, trying default one\n", value);
1487             value = 1;
1488         }
1489         s->fill_order = value - 1;
1490         break;
1491     case TIFF_PAL: {
1492         GetByteContext pal_gb[3];
1493         off = type_sizes[type];
1494         if (count / 3 > 256 ||
1495             bytestream2_get_bytes_left(&s->gb) < count / 3 * off * 3)
1496             return AVERROR_INVALIDDATA;
1497
1498         pal_gb[0] = pal_gb[1] = pal_gb[2] = s->gb;
1499         bytestream2_skip(&pal_gb[1], count / 3 * off);
1500         bytestream2_skip(&pal_gb[2], count / 3 * off * 2);
1501
1502         off = (type_sizes[type] - 1) << 3;
1503         if (off > 31U) {
1504             av_log(s->avctx, AV_LOG_ERROR, "palette shift %d is out of range\n", off);
1505             return AVERROR_INVALIDDATA;
1506         }
1507
1508         for (i = 0; i < count / 3; i++) {
1509             uint32_t p = 0xFF000000;
1510             p |= (ff_tget(&pal_gb[0], type, s->le) >> off) << 16;
1511             p |= (ff_tget(&pal_gb[1], type, s->le) >> off) << 8;
1512             p |=  ff_tget(&pal_gb[2], type, s->le) >> off;
1513             s->palette[i] = p;
1514         }
1515         s->palette_is_set = 1;
1516         break;
1517     }
1518     case TIFF_PLANAR:
1519         s->planar = value == 2;
1520         break;
1521     case TIFF_YCBCR_SUBSAMPLING:
1522         if (count != 2) {
1523             av_log(s->avctx, AV_LOG_ERROR, "subsample count invalid\n");
1524             return AVERROR_INVALIDDATA;
1525         }
1526         for (i = 0; i < count; i++) {
1527             s->subsampling[i] = ff_tget(&s->gb, type, s->le);
1528             if (s->subsampling[i] <= 0) {
1529                 av_log(s->avctx, AV_LOG_ERROR, "subsampling %d is invalid\n", s->subsampling[i]);
1530                 s->subsampling[i] = 1;
1531                 return AVERROR_INVALIDDATA;
1532             }
1533         }
1534         break;
1535     case TIFF_T4OPTIONS:
1536         if (s->compr == TIFF_G3)
1537             s->fax_opts = value;
1538         break;
1539     case TIFF_T6OPTIONS:
1540         if (s->compr == TIFF_G4)
1541             s->fax_opts = value;
1542         break;
1543 #define ADD_METADATA(count, name, sep)\
1544     if ((ret = add_metadata(count, type, name, sep, s, frame)) < 0) {\
1545         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");\
1546         goto end;\
1547     }
1548     case TIFF_MODEL_PIXEL_SCALE:
1549         ADD_METADATA(count, "ModelPixelScaleTag", NULL);
1550         break;
1551     case TIFF_MODEL_TRANSFORMATION:
1552         ADD_METADATA(count, "ModelTransformationTag", NULL);
1553         break;
1554     case TIFF_MODEL_TIEPOINT:
1555         ADD_METADATA(count, "ModelTiepointTag", NULL);
1556         break;
1557     case TIFF_GEO_KEY_DIRECTORY:
1558         if (s->geotag_count) {
1559             avpriv_request_sample(s->avctx, "Multiple geo key directories\n");
1560             return AVERROR_INVALIDDATA;
1561         }
1562         ADD_METADATA(1, "GeoTIFF_Version", NULL);
1563         ADD_METADATA(2, "GeoTIFF_Key_Revision", ".");
1564         s->geotag_count   = ff_tget_short(&s->gb, s->le);
1565         if (s->geotag_count > count / 4 - 1) {
1566             s->geotag_count = count / 4 - 1;
1567             av_log(s->avctx, AV_LOG_WARNING, "GeoTIFF key directory buffer shorter than specified\n");
1568         }
1569         if (   bytestream2_get_bytes_left(&s->gb) < s->geotag_count * sizeof(int16_t) * 4
1570             || s->geotag_count == 0) {
1571             s->geotag_count = 0;
1572             return -1;
1573         }
1574         s->geotags = av_mallocz_array(s->geotag_count, sizeof(TiffGeoTag));
1575         if (!s->geotags) {
1576             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1577             s->geotag_count = 0;
1578             goto end;
1579         }
1580         for (i = 0; i < s->geotag_count; i++) {
1581             s->geotags[i].key    = ff_tget_short(&s->gb, s->le);
1582             s->geotags[i].type   = ff_tget_short(&s->gb, s->le);
1583             s->geotags[i].count  = ff_tget_short(&s->gb, s->le);
1584
1585             if (!s->geotags[i].type)
1586                 s->geotags[i].val  = get_geokey_val(s->geotags[i].key, ff_tget_short(&s->gb, s->le));
1587             else
1588                 s->geotags[i].offset = ff_tget_short(&s->gb, s->le);
1589         }
1590         break;
1591     case TIFF_GEO_DOUBLE_PARAMS:
1592         if (count >= INT_MAX / sizeof(int64_t))
1593             return AVERROR_INVALIDDATA;
1594         if (bytestream2_get_bytes_left(&s->gb) < count * sizeof(int64_t))
1595             return AVERROR_INVALIDDATA;
1596         dp = av_malloc_array(count, sizeof(double));
1597         if (!dp) {
1598             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1599             goto end;
1600         }
1601         for (i = 0; i < count; i++)
1602             dp[i] = ff_tget_double(&s->gb, s->le);
1603         for (i = 0; i < s->geotag_count; i++) {
1604             if (s->geotags[i].type == TIFF_GEO_DOUBLE_PARAMS) {
1605                 if (s->geotags[i].count == 0
1606                     || s->geotags[i].offset + s->geotags[i].count > count) {
1607                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1608                 } else if (s->geotags[i].val) {
1609                     av_log(s->avctx, AV_LOG_WARNING, "Duplicate GeoTIFF key %d\n", s->geotags[i].key);
1610                 } else {
1611                     char *ap = doubles2str(&dp[s->geotags[i].offset], s->geotags[i].count, ", ");
1612                     if (!ap) {
1613                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1614                         av_freep(&dp);
1615                         return AVERROR(ENOMEM);
1616                     }
1617                     s->geotags[i].val = ap;
1618                 }
1619             }
1620         }
1621         av_freep(&dp);
1622         break;
1623     case TIFF_GEO_ASCII_PARAMS:
1624         pos = bytestream2_tell(&s->gb);
1625         for (i = 0; i < s->geotag_count; i++) {
1626             if (s->geotags[i].type == TIFF_GEO_ASCII_PARAMS) {
1627                 if (s->geotags[i].count == 0
1628                     || s->geotags[i].offset +  s->geotags[i].count > count) {
1629                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1630                 } else {
1631                     char *ap;
1632
1633                     bytestream2_seek(&s->gb, pos + s->geotags[i].offset, SEEK_SET);
1634                     if (bytestream2_get_bytes_left(&s->gb) < s->geotags[i].count)
1635                         return AVERROR_INVALIDDATA;
1636                     if (s->geotags[i].val)
1637                         return AVERROR_INVALIDDATA;
1638                     ap = av_malloc(s->geotags[i].count);
1639                     if (!ap) {
1640                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1641                         return AVERROR(ENOMEM);
1642                     }
1643                     bytestream2_get_bufferu(&s->gb, ap, s->geotags[i].count);
1644                     ap[s->geotags[i].count - 1] = '\0'; //replace the "|" delimiter with a 0 byte
1645                     s->geotags[i].val = ap;
1646                 }
1647             }
1648         }
1649         break;
1650     case TIFF_ICC_PROFILE:
1651         if (type != TIFF_UNDEFINED)
1652             return AVERROR_INVALIDDATA;
1653
1654         gb_temp = s->gb;
1655         bytestream2_seek(&gb_temp, SEEK_SET, off);
1656
1657         if (bytestream2_get_bytes_left(&gb_temp) < count)
1658             return AVERROR_INVALIDDATA;
1659
1660         sd = av_frame_new_side_data(frame, AV_FRAME_DATA_ICC_PROFILE, count);
1661         if (!sd)
1662             return AVERROR(ENOMEM);
1663
1664         bytestream2_get_bufferu(&gb_temp, sd->data, count);
1665         break;
1666     case TIFF_ARTIST:
1667         ADD_METADATA(count, "artist", NULL);
1668         break;
1669     case TIFF_COPYRIGHT:
1670         ADD_METADATA(count, "copyright", NULL);
1671         break;
1672     case TIFF_DATE:
1673         ADD_METADATA(count, "date", NULL);
1674         break;
1675     case TIFF_DOCUMENT_NAME:
1676         ADD_METADATA(count, "document_name", NULL);
1677         break;
1678     case TIFF_HOST_COMPUTER:
1679         ADD_METADATA(count, "computer", NULL);
1680         break;
1681     case TIFF_IMAGE_DESCRIPTION:
1682         ADD_METADATA(count, "description", NULL);
1683         break;
1684     case TIFF_MAKE:
1685         ADD_METADATA(count, "make", NULL);
1686         break;
1687     case TIFF_MODEL:
1688         ADD_METADATA(count, "model", NULL);
1689         break;
1690     case TIFF_PAGE_NAME:
1691         ADD_METADATA(count, "page_name", NULL);
1692         break;
1693     case TIFF_PAGE_NUMBER:
1694         ADD_METADATA(count, "page_number", " / ");
1695         // need to seek back to re-read the page number
1696         bytestream2_seek(&s->gb, -count * sizeof(uint16_t), SEEK_CUR);
1697         // read the page number
1698         s->cur_page = ff_tget(&s->gb, TIFF_SHORT, s->le);
1699         // get back to where we were before the previous seek
1700         bytestream2_seek(&s->gb, count * sizeof(uint16_t) - sizeof(uint16_t), SEEK_CUR);
1701         break;
1702     case TIFF_SOFTWARE_NAME:
1703         ADD_METADATA(count, "software", NULL);
1704         break;
1705     case DNG_VERSION:
1706         if (count == 4) {
1707             unsigned int ver[4];
1708             ver[0] = ff_tget(&s->gb, type, s->le);
1709             ver[1] = ff_tget(&s->gb, type, s->le);
1710             ver[2] = ff_tget(&s->gb, type, s->le);
1711             ver[3] = ff_tget(&s->gb, type, s->le);
1712
1713             av_log(s->avctx, AV_LOG_DEBUG, "DNG file, version %u.%u.%u.%u\n",
1714                 ver[0], ver[1], ver[2], ver[3]);
1715
1716             tiff_set_type(s, TIFF_TYPE_DNG);
1717         }
1718         break;
1719     case CINEMADNG_TIME_CODES:
1720     case CINEMADNG_FRAME_RATE:
1721     case CINEMADNG_T_STOP:
1722     case CINEMADNG_REEL_NAME:
1723     case CINEMADNG_CAMERA_LABEL:
1724         tiff_set_type(s, TIFF_TYPE_CINEMADNG);
1725         break;
1726     default:
1727         if (s->avctx->err_recognition & AV_EF_EXPLODE) {
1728             av_log(s->avctx, AV_LOG_ERROR,
1729                    "Unknown or unsupported tag %d/0x%0X\n",
1730                    tag, tag);
1731             return AVERROR_INVALIDDATA;
1732         }
1733     }
1734 end:
1735     if (s->bpp > 64U) {
1736         av_log(s->avctx, AV_LOG_ERROR,
1737                 "This format is not supported (bpp=%d, %d components)\n",
1738                 s->bpp, count);
1739         s->bpp = 0;
1740         return AVERROR_INVALIDDATA;
1741     }
1742     bytestream2_seek(&s->gb, start, SEEK_SET);
1743     return 0;
1744 }
1745
1746 static int decode_frame(AVCodecContext *avctx,
1747                         void *data, int *got_frame, AVPacket *avpkt)
1748 {
1749     TiffContext *const s = avctx->priv_data;
1750     AVFrame *const p = data;
1751     ThreadFrame frame = { .f = data };
1752     unsigned off, last_off;
1753     int le, ret, plane, planes;
1754     int i, j, entries, stride;
1755     unsigned soff, ssize;
1756     uint8_t *dst;
1757     GetByteContext stripsizes;
1758     GetByteContext stripdata;
1759     int retry_for_subifd, retry_for_page;
1760     int is_dng;
1761
1762     bytestream2_init(&s->gb, avpkt->data, avpkt->size);
1763
1764     // parse image header
1765     if ((ret = ff_tdecode_header(&s->gb, &le, &off))) {
1766         av_log(avctx, AV_LOG_ERROR, "Invalid TIFF header\n");
1767         return ret;
1768     } else if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1769         av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1770         return AVERROR_INVALIDDATA;
1771     }
1772     s->le          = le;
1773     // TIFF_BPP is not a required tag and defaults to 1
1774
1775     s->tiff_type   = TIFF_TYPE_TIFF;
1776 again:
1777     s->is_thumbnail = 0;
1778     s->bppcount    = s->bpp = 1;
1779     s->photometric = TIFF_PHOTOMETRIC_NONE;
1780     s->compr       = TIFF_RAW;
1781     s->fill_order  = 0;
1782     s->white_level = 0;
1783     s->is_bayer    = 0;
1784     s->is_tiled    = 0;
1785     s->is_jpeg     = 0;
1786     s->cur_page    = 0;
1787
1788     for (i = 0; i < 65536; i++)
1789         s->dng_lut[i] = i;
1790
1791     free_geotags(s);
1792
1793     // Reset these offsets so we can tell if they were set this frame
1794     s->stripsizesoff = s->strippos = 0;
1795     /* parse image file directory */
1796     bytestream2_seek(&s->gb, off, SEEK_SET);
1797     entries = ff_tget_short(&s->gb, le);
1798     if (bytestream2_get_bytes_left(&s->gb) < entries * 12)
1799         return AVERROR_INVALIDDATA;
1800     for (i = 0; i < entries; i++) {
1801         if ((ret = tiff_decode_tag(s, p)) < 0)
1802             return ret;
1803     }
1804
1805     if (s->get_thumbnail && !s->is_thumbnail) {
1806         av_log(avctx, AV_LOG_INFO, "No embedded thumbnail present\n");
1807         return AVERROR_EOF;
1808     }
1809
1810     /** whether we should process this IFD's SubIFD */
1811     retry_for_subifd = s->sub_ifd && (s->get_subimage || (!s->get_thumbnail && s->is_thumbnail));
1812     /** whether we should process this multi-page IFD's next page */
1813     retry_for_page = s->get_page && s->cur_page + 1 < s->get_page;  // get_page is 1-indexed
1814
1815     last_off = off;
1816     if (retry_for_page) {
1817         // set offset to the next IFD
1818         off = ff_tget_long(&s->gb, le);
1819     } else if (retry_for_subifd) {
1820         // set offset to the SubIFD
1821         off = s->sub_ifd;
1822     }
1823
1824     if (retry_for_subifd || retry_for_page) {
1825         if (!off) {
1826             av_log(avctx, AV_LOG_ERROR, "Requested entry not found\n");
1827             return AVERROR_INVALIDDATA;
1828         }
1829         if (off <= last_off) {
1830             avpriv_request_sample(s->avctx, "non increasing IFD offset\n");
1831             return AVERROR_INVALIDDATA;
1832         }
1833         if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1834             av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1835             return AVERROR_INVALIDDATA;
1836         }
1837         s->sub_ifd = 0;
1838         goto again;
1839     }
1840
1841     /* At this point we've decided on which (Sub)IFD to process */
1842
1843     is_dng = (s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG);
1844
1845     for (i = 0; i<s->geotag_count; i++) {
1846         const char *keyname = get_geokey_name(s->geotags[i].key);
1847         if (!keyname) {
1848             av_log(avctx, AV_LOG_WARNING, "Unknown or unsupported GeoTIFF key %d\n", s->geotags[i].key);
1849             continue;
1850         }
1851         if (get_geokey_type(s->geotags[i].key) != s->geotags[i].type) {
1852             av_log(avctx, AV_LOG_WARNING, "Type of GeoTIFF key %d is wrong\n", s->geotags[i].key);
1853             continue;
1854         }
1855         ret = av_dict_set(&p->metadata, keyname, s->geotags[i].val, 0);
1856         if (ret<0) {
1857             av_log(avctx, AV_LOG_ERROR, "Writing metadata with key '%s' failed\n", keyname);
1858             return ret;
1859         }
1860     }
1861
1862     if (is_dng) {
1863         int bps;
1864
1865         if (s->white_level == 0)
1866             s->white_level = (1 << s->bpp) - 1; /* Default value as per the spec */
1867
1868         if (s->white_level <= s->black_level) {
1869             av_log(avctx, AV_LOG_ERROR, "BlackLevel (%"PRId32") must be less than WhiteLevel (%"PRId32")\n",
1870                 s->black_level, s->white_level);
1871             return AVERROR_INVALIDDATA;
1872         }
1873
1874         if (s->bpp % s->bppcount)
1875             return AVERROR_INVALIDDATA;
1876         bps = s->bpp / s->bppcount;
1877         if (bps < 8 || bps > 32)
1878             return AVERROR_INVALIDDATA;
1879         if (s->planar)
1880             return AVERROR_PATCHWELCOME;
1881     }
1882
1883     if (!s->is_tiled && !s->strippos && !s->stripoff) {
1884         av_log(avctx, AV_LOG_ERROR, "Image data is missing\n");
1885         return AVERROR_INVALIDDATA;
1886     }
1887
1888     /* now we have the data and may start decoding */
1889     if ((ret = init_image(s, &frame)) < 0)
1890         return ret;
1891
1892     if (!s->is_tiled) {
1893         if (s->strips == 1 && !s->stripsize) {
1894             av_log(avctx, AV_LOG_WARNING, "Image data size missing\n");
1895             s->stripsize = avpkt->size - s->stripoff;
1896         }
1897
1898         if (s->stripsizesoff) {
1899             if (s->stripsizesoff >= (unsigned)avpkt->size)
1900                 return AVERROR_INVALIDDATA;
1901             bytestream2_init(&stripsizes, avpkt->data + s->stripsizesoff,
1902                             avpkt->size - s->stripsizesoff);
1903         }
1904         if (s->strippos) {
1905             if (s->strippos >= (unsigned)avpkt->size)
1906                 return AVERROR_INVALIDDATA;
1907             bytestream2_init(&stripdata, avpkt->data + s->strippos,
1908                             avpkt->size - s->strippos);
1909         }
1910
1911         if (s->rps <= 0 || s->rps % s->subsampling[1]) {
1912             av_log(avctx, AV_LOG_ERROR, "rps %d invalid\n", s->rps);
1913             return AVERROR_INVALIDDATA;
1914         }
1915     }
1916
1917     if (s->photometric == TIFF_PHOTOMETRIC_LINEAR_RAW ||
1918         s->photometric == TIFF_PHOTOMETRIC_CFA) {
1919         p->color_trc = AVCOL_TRC_LINEAR;
1920     } else if (s->photometric == TIFF_PHOTOMETRIC_BLACK_IS_ZERO) {
1921         p->color_trc = AVCOL_TRC_GAMMA22;
1922     }
1923
1924     /* Handle DNG images with JPEG-compressed tiles */
1925
1926     if (is_dng && s->is_tiled) {
1927         if (!s->is_jpeg) {
1928             avpriv_report_missing_feature(avctx, "DNG uncompressed tiled images");
1929             return AVERROR_PATCHWELCOME;
1930         } else if (!s->is_bayer) {
1931             avpriv_report_missing_feature(avctx, "DNG JPG-compressed tiled non-bayer-encoded images");
1932             return AVERROR_PATCHWELCOME;
1933         } else {
1934             if ((ret = dng_decode_tiles(avctx, (AVFrame*)data, avpkt)) > 0)
1935                 *got_frame = 1;
1936             return ret;
1937         }
1938     }
1939
1940     /* Handle TIFF images and DNG images with uncompressed strips (non-tiled) */
1941
1942     planes = s->planar ? s->bppcount : 1;
1943     for (plane = 0; plane < planes; plane++) {
1944         uint8_t *five_planes = NULL;
1945         int remaining = avpkt->size;
1946         int decoded_height;
1947         stride = p->linesize[plane];
1948         dst = p->data[plane];
1949         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1950             s->avctx->pix_fmt == AV_PIX_FMT_RGBA) {
1951             stride = stride * 5 / 4;
1952             five_planes =
1953             dst = av_malloc(stride * s->height);
1954             if (!dst)
1955                 return AVERROR(ENOMEM);
1956         }
1957         for (i = 0; i < s->height; i += s->rps) {
1958             if (i)
1959                 dst += s->rps * stride;
1960             if (s->stripsizesoff)
1961                 ssize = ff_tget(&stripsizes, s->sstype, le);
1962             else
1963                 ssize = s->stripsize;
1964
1965             if (s->strippos)
1966                 soff = ff_tget(&stripdata, s->sot, le);
1967             else
1968                 soff = s->stripoff;
1969
1970             if (soff > avpkt->size || ssize > avpkt->size - soff || ssize > remaining) {
1971                 av_log(avctx, AV_LOG_ERROR, "Invalid strip size/offset\n");
1972                 av_freep(&five_planes);
1973                 return AVERROR_INVALIDDATA;
1974             }
1975             remaining -= ssize;
1976             if ((ret = tiff_unpack_strip(s, p, dst, stride, avpkt->data + soff, ssize, i,
1977                                          FFMIN(s->rps, s->height - i))) < 0) {
1978                 if (avctx->err_recognition & AV_EF_EXPLODE) {
1979                     av_freep(&five_planes);
1980                     return ret;
1981                 }
1982                 break;
1983             }
1984         }
1985         decoded_height = FFMIN(i, s->height);
1986
1987         if (s->predictor == 2) {
1988             if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1989                 av_log(s->avctx, AV_LOG_ERROR, "predictor == 2 with YUV is unsupported");
1990                 return AVERROR_PATCHWELCOME;
1991             }
1992             dst   = five_planes ? five_planes : p->data[plane];
1993             soff  = s->bpp >> 3;
1994             if (s->planar)
1995                 soff  = FFMAX(soff / s->bppcount, 1);
1996             ssize = s->width * soff;
1997             if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48LE ||
1998                 s->avctx->pix_fmt == AV_PIX_FMT_RGBA64LE ||
1999                 s->avctx->pix_fmt == AV_PIX_FMT_GRAY16LE ||
2000                 s->avctx->pix_fmt == AV_PIX_FMT_YA16LE ||
2001                 s->avctx->pix_fmt == AV_PIX_FMT_GBRP16LE ||
2002                 s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16LE) {
2003                 for (i = 0; i < decoded_height; i++) {
2004                     for (j = soff; j < ssize; j += 2)
2005                         AV_WL16(dst + j, AV_RL16(dst + j) + AV_RL16(dst + j - soff));
2006                     dst += stride;
2007                 }
2008             } else if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48BE ||
2009                        s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE ||
2010                        s->avctx->pix_fmt == AV_PIX_FMT_GRAY16BE ||
2011                        s->avctx->pix_fmt == AV_PIX_FMT_YA16BE ||
2012                        s->avctx->pix_fmt == AV_PIX_FMT_GBRP16BE ||
2013                        s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16BE) {
2014                 for (i = 0; i < decoded_height; i++) {
2015                     for (j = soff; j < ssize; j += 2)
2016                         AV_WB16(dst + j, AV_RB16(dst + j) + AV_RB16(dst + j - soff));
2017                     dst += stride;
2018                 }
2019             } else {
2020                 for (i = 0; i < decoded_height; i++) {
2021                     for (j = soff; j < ssize; j++)
2022                         dst[j] += dst[j - soff];
2023                     dst += stride;
2024                 }
2025             }
2026         }
2027
2028         if (s->photometric == TIFF_PHOTOMETRIC_WHITE_IS_ZERO) {
2029             int c = (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 ? (1<<s->bpp) - 1 : 255);
2030             dst = p->data[plane];
2031             for (i = 0; i < s->height; i++) {
2032                 for (j = 0; j < stride; j++)
2033                     dst[j] = c - dst[j];
2034                 dst += stride;
2035             }
2036         }
2037
2038         if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
2039             (s->avctx->pix_fmt == AV_PIX_FMT_RGB0 || s->avctx->pix_fmt == AV_PIX_FMT_RGBA)) {
2040             int x = s->avctx->pix_fmt == AV_PIX_FMT_RGB0 ? 4 : 5;
2041             uint8_t *src = five_planes ? five_planes : p->data[plane];
2042             dst = p->data[plane];
2043             for (i = 0; i < s->height; i++) {
2044                 for (j = 0; j < s->width; j++) {
2045                     int k =  255 - src[x * j + 3];
2046                     int r = (255 - src[x * j    ]) * k;
2047                     int g = (255 - src[x * j + 1]) * k;
2048                     int b = (255 - src[x * j + 2]) * k;
2049                     dst[4 * j    ] = r * 257 >> 16;
2050                     dst[4 * j + 1] = g * 257 >> 16;
2051                     dst[4 * j + 2] = b * 257 >> 16;
2052                     dst[4 * j + 3] = s->avctx->pix_fmt == AV_PIX_FMT_RGBA ? src[x * j + 4] : 255;
2053                 }
2054                 src += stride;
2055                 dst += p->linesize[plane];
2056             }
2057             av_freep(&five_planes);
2058         } else if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
2059             s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE) {
2060             dst = p->data[plane];
2061             for (i = 0; i < s->height; i++) {
2062                 for (j = 0; j < s->width; j++) {
2063                     uint64_t k =  65535 - AV_RB16(dst + 8 * j + 6);
2064                     uint64_t r = (65535 - AV_RB16(dst + 8 * j    )) * k;
2065                     uint64_t g = (65535 - AV_RB16(dst + 8 * j + 2)) * k;
2066                     uint64_t b = (65535 - AV_RB16(dst + 8 * j + 4)) * k;
2067                     AV_WB16(dst + 8 * j    , r * 65537 >> 32);
2068                     AV_WB16(dst + 8 * j + 2, g * 65537 >> 32);
2069                     AV_WB16(dst + 8 * j + 4, b * 65537 >> 32);
2070                     AV_WB16(dst + 8 * j + 6, 65535);
2071                 }
2072                 dst += p->linesize[plane];
2073             }
2074         }
2075     }
2076
2077     if (s->planar && s->bppcount > 2) {
2078         FFSWAP(uint8_t*, p->data[0],     p->data[2]);
2079         FFSWAP(int,      p->linesize[0], p->linesize[2]);
2080         FFSWAP(uint8_t*, p->data[0],     p->data[1]);
2081         FFSWAP(int,      p->linesize[0], p->linesize[1]);
2082     }
2083
2084     if (s->is_bayer && s->white_level && s->bpp == 16 && !is_dng) {
2085         uint16_t *dst = (uint16_t *)p->data[0];
2086         for (i = 0; i < s->height; i++) {
2087             for (j = 0; j < s->width; j++)
2088                 dst[j] = FFMIN((dst[j] / (float)s->white_level) * 65535, 65535);
2089             dst += stride / 2;
2090         }
2091     }
2092
2093     *got_frame = 1;
2094
2095     return avpkt->size;
2096 }
2097
2098 static av_cold int tiff_init(AVCodecContext *avctx)
2099 {
2100     TiffContext *s = avctx->priv_data;
2101     const AVCodec *codec;
2102     int ret;
2103
2104     s->width  = 0;
2105     s->height = 0;
2106     s->subsampling[0] =
2107     s->subsampling[1] = 1;
2108     s->avctx  = avctx;
2109     ff_lzw_decode_open(&s->lzw);
2110     if (!s->lzw)
2111         return AVERROR(ENOMEM);
2112     ff_ccitt_unpack_init();
2113
2114     /* Allocate JPEG frame */
2115     s->jpgframe = av_frame_alloc();
2116     if (!s->jpgframe)
2117         return AVERROR(ENOMEM);
2118
2119     /* Prepare everything needed for JPEG decoding */
2120     codec = avcodec_find_decoder(AV_CODEC_ID_MJPEG);
2121     if (!codec)
2122         return AVERROR_BUG;
2123     s->avctx_mjpeg = avcodec_alloc_context3(codec);
2124     if (!s->avctx_mjpeg)
2125         return AVERROR(ENOMEM);
2126     s->avctx_mjpeg->flags = avctx->flags;
2127     s->avctx_mjpeg->flags2 = avctx->flags2;
2128     s->avctx_mjpeg->dct_algo = avctx->dct_algo;
2129     s->avctx_mjpeg->idct_algo = avctx->idct_algo;
2130     ret = ff_codec_open2_recursive(s->avctx_mjpeg, codec, NULL);
2131     if (ret < 0) {
2132         return ret;
2133     }
2134
2135     return 0;
2136 }
2137
2138 static av_cold int tiff_end(AVCodecContext *avctx)
2139 {
2140     TiffContext *const s = avctx->priv_data;
2141
2142     free_geotags(s);
2143
2144     ff_lzw_decode_close(&s->lzw);
2145     av_freep(&s->deinvert_buf);
2146     s->deinvert_buf_size = 0;
2147     av_freep(&s->yuv_line);
2148     s->yuv_line_size = 0;
2149     av_freep(&s->fax_buffer);
2150     s->fax_buffer_size = 0;
2151     av_frame_free(&s->jpgframe);
2152     avcodec_free_context(&s->avctx_mjpeg);
2153     return 0;
2154 }
2155
2156 #define OFFSET(x) offsetof(TiffContext, x)
2157 static const AVOption tiff_options[] = {
2158     { "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 },
2159     { "thumbnail", "decode embedded thumbnail subimage instead if available", OFFSET(get_thumbnail), AV_OPT_TYPE_BOOL, {.i64=0},  0, 1, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2160     { "page", "page number of multi-page image to decode (starting from 1)", OFFSET(get_page), AV_OPT_TYPE_INT, {.i64=0}, 0, UINT16_MAX, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2161     { NULL },
2162 };
2163
2164 static const AVClass tiff_decoder_class = {
2165     .class_name = "TIFF decoder",
2166     .item_name  = av_default_item_name,
2167     .option     = tiff_options,
2168     .version    = LIBAVUTIL_VERSION_INT,
2169 };
2170
2171 AVCodec ff_tiff_decoder = {
2172     .name           = "tiff",
2173     .long_name      = NULL_IF_CONFIG_SMALL("TIFF image"),
2174     .type           = AVMEDIA_TYPE_VIDEO,
2175     .id             = AV_CODEC_ID_TIFF,
2176     .priv_data_size = sizeof(TiffContext),
2177     .init           = tiff_init,
2178     .close          = tiff_end,
2179     .decode         = decode_frame,
2180     .capabilities   = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_FRAME_THREADS,
2181     .caps_internal  = FF_CODEC_CAP_INIT_CLEANUP,
2182     .priv_class     = &tiff_decoder_class,
2183 };