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