lavc/tiff: add support for LZMA compression
[ffmpeg.git] / libavcodec / tiff.c
1 /*
2  * Copyright (c) 2006 Konstantin Shishkov
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20
21 /**
22  * @file
23  * TIFF image decoder
24  * @author Konstantin Shishkov
25  */
26
27 #include "config.h"
28 #if CONFIG_ZLIB
29 #include <zlib.h>
30 #endif
31 #if CONFIG_LZMA
32 #include <lzma.h>
33 #endif
34
35 #include "libavutil/attributes.h"
36 #include "libavutil/avstring.h"
37 #include "libavutil/intreadwrite.h"
38 #include "libavutil/imgutils.h"
39 #include "avcodec.h"
40 #include "bytestream.h"
41 #include "faxcompr.h"
42 #include "internal.h"
43 #include "lzw.h"
44 #include "mathops.h"
45 #include "tiff.h"
46 #include "tiff_data.h"
47 #include "thread.h"
48
49 typedef struct TiffContext {
50     AVCodecContext *avctx;
51     GetByteContext gb;
52
53     int width, height;
54     unsigned int bpp, bppcount;
55     uint32_t palette[256];
56     int palette_is_set;
57     int le;
58     enum TiffCompr compr;
59     enum TiffPhotometric photometric;
60     int planar;
61     int subsampling[2];
62     int fax_opts;
63     int predictor;
64     int fill_order;
65     uint32_t res[4];
66
67     int strips, rps, sstype;
68     int sot;
69     int stripsizesoff, stripsize, stripoff, strippos;
70     LZWState *lzw;
71
72     uint8_t *deinvert_buf;
73     int deinvert_buf_size;
74     uint8_t *yuv_line;
75     unsigned int yuv_line_size;
76
77     int geotag_count;
78     TiffGeoTag *geotags;
79 } TiffContext;
80
81 static void free_geotags(TiffContext *const s)
82 {
83     int i;
84     for (i = 0; i < s->geotag_count; i++) {
85         if (s->geotags[i].val)
86             av_freep(&s->geotags[i].val);
87     }
88     av_freep(&s->geotags);
89     s->geotag_count = 0;
90 }
91
92 #define RET_GEOKEY(TYPE, array, element)\
93     if (key >= TIFF_##TYPE##_KEY_ID_OFFSET &&\
94         key - TIFF_##TYPE##_KEY_ID_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_name_type_map))\
95         return ff_tiff_##array##_name_type_map[key - TIFF_##TYPE##_KEY_ID_OFFSET].element;
96
97 static const char *get_geokey_name(int key)
98 {
99     RET_GEOKEY(VERT, vert, name);
100     RET_GEOKEY(PROJ, proj, name);
101     RET_GEOKEY(GEOG, geog, name);
102     RET_GEOKEY(CONF, conf, name);
103
104     return NULL;
105 }
106
107 static int get_geokey_type(int key)
108 {
109     RET_GEOKEY(VERT, vert, type);
110     RET_GEOKEY(PROJ, proj, type);
111     RET_GEOKEY(GEOG, geog, type);
112     RET_GEOKEY(CONF, conf, type);
113
114     return AVERROR_INVALIDDATA;
115 }
116
117 static int cmp_id_key(const void *id, const void *k)
118 {
119     return *(const int*)id - ((const TiffGeoTagKeyName*)k)->key;
120 }
121
122 static const char *search_keyval(const TiffGeoTagKeyName *keys, int n, int id)
123 {
124     TiffGeoTagKeyName *r = bsearch(&id, keys, n, sizeof(keys[0]), cmp_id_key);
125     if(r)
126         return r->name;
127
128     return NULL;
129 }
130
131 static char *get_geokey_val(int key, int val)
132 {
133     char *ap;
134
135     if (val == TIFF_GEO_KEY_UNDEFINED)
136         return av_strdup("undefined");
137     if (val == TIFF_GEO_KEY_USER_DEFINED)
138         return av_strdup("User-Defined");
139
140 #define RET_GEOKEY_VAL(TYPE, array)\
141     if (val >= TIFF_##TYPE##_OFFSET &&\
142         val - TIFF_##TYPE##_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_codes))\
143         return av_strdup(ff_tiff_##array##_codes[val - TIFF_##TYPE##_OFFSET]);
144
145     switch (key) {
146     case TIFF_GT_MODEL_TYPE_GEOKEY:
147         RET_GEOKEY_VAL(GT_MODEL_TYPE, gt_model_type);
148         break;
149     case TIFF_GT_RASTER_TYPE_GEOKEY:
150         RET_GEOKEY_VAL(GT_RASTER_TYPE, gt_raster_type);
151         break;
152     case TIFF_GEOG_LINEAR_UNITS_GEOKEY:
153     case TIFF_PROJ_LINEAR_UNITS_GEOKEY:
154     case TIFF_VERTICAL_UNITS_GEOKEY:
155         RET_GEOKEY_VAL(LINEAR_UNIT, linear_unit);
156         break;
157     case TIFF_GEOG_ANGULAR_UNITS_GEOKEY:
158     case TIFF_GEOG_AZIMUTH_UNITS_GEOKEY:
159         RET_GEOKEY_VAL(ANGULAR_UNIT, angular_unit);
160         break;
161     case TIFF_GEOGRAPHIC_TYPE_GEOKEY:
162         RET_GEOKEY_VAL(GCS_TYPE, gcs_type);
163         RET_GEOKEY_VAL(GCSE_TYPE, gcse_type);
164         break;
165     case TIFF_GEOG_GEODETIC_DATUM_GEOKEY:
166         RET_GEOKEY_VAL(GEODETIC_DATUM, geodetic_datum);
167         RET_GEOKEY_VAL(GEODETIC_DATUM_E, geodetic_datum_e);
168         break;
169     case TIFF_GEOG_ELLIPSOID_GEOKEY:
170         RET_GEOKEY_VAL(ELLIPSOID, ellipsoid);
171         break;
172     case TIFF_GEOG_PRIME_MERIDIAN_GEOKEY:
173         RET_GEOKEY_VAL(PRIME_MERIDIAN, prime_meridian);
174         break;
175     case TIFF_PROJECTED_CS_TYPE_GEOKEY:
176         ap = av_strdup(search_keyval(ff_tiff_proj_cs_type_codes, FF_ARRAY_ELEMS(ff_tiff_proj_cs_type_codes), val));
177         if(ap) return ap;
178         break;
179     case TIFF_PROJECTION_GEOKEY:
180         ap = av_strdup(search_keyval(ff_tiff_projection_codes, FF_ARRAY_ELEMS(ff_tiff_projection_codes), val));
181         if(ap) return ap;
182         break;
183     case TIFF_PROJ_COORD_TRANS_GEOKEY:
184         RET_GEOKEY_VAL(COORD_TRANS, coord_trans);
185         break;
186     case TIFF_VERTICAL_CS_TYPE_GEOKEY:
187         RET_GEOKEY_VAL(VERT_CS, vert_cs);
188         RET_GEOKEY_VAL(ORTHO_VERT_CS, ortho_vert_cs);
189         break;
190
191     }
192
193     ap = av_malloc(14);
194     if (ap)
195         snprintf(ap, 14, "Unknown-%d", val);
196     return ap;
197 }
198
199 static char *doubles2str(double *dp, int count, const char *sep)
200 {
201     int i;
202     char *ap, *ap0;
203     uint64_t component_len;
204     if (!sep) sep = ", ";
205     component_len = 24LL + strlen(sep);
206     if (count >= (INT_MAX - 1)/component_len)
207         return NULL;
208     ap = av_malloc(component_len * count + 1);
209     if (!ap)
210         return NULL;
211     ap0   = ap;
212     ap[0] = '\0';
213     for (i = 0; i < count; i++) {
214         unsigned l = snprintf(ap, component_len, "%.15g%s", dp[i], sep);
215         if(l >= component_len) {
216             av_free(ap0);
217             return NULL;
218         }
219         ap += l;
220     }
221     ap0[strlen(ap0) - strlen(sep)] = '\0';
222     return ap0;
223 }
224
225 static int add_metadata(int count, int type,
226                         const char *name, const char *sep, TiffContext *s, AVFrame *frame)
227 {
228     switch(type) {
229     case TIFF_DOUBLE: return ff_tadd_doubles_metadata(count, name, sep, &s->gb, s->le, avpriv_frame_get_metadatap(frame));
230     case TIFF_SHORT : return ff_tadd_shorts_metadata(count, name, sep, &s->gb, s->le, 0, avpriv_frame_get_metadatap(frame));
231     case TIFF_STRING: return ff_tadd_string_metadata(count, name, &s->gb, s->le, avpriv_frame_get_metadatap(frame));
232     default         : return AVERROR_INVALIDDATA;
233     };
234 }
235
236 static void av_always_inline horizontal_fill(unsigned int bpp, uint8_t* dst,
237                                              int usePtr, const uint8_t *src,
238                                              uint8_t c, int width, int offset)
239 {
240     switch (bpp) {
241     case 1:
242         while (--width >= 0) {
243             dst[(width+offset)*8+7] = (usePtr ? src[width] : c)      & 0x1;
244             dst[(width+offset)*8+6] = (usePtr ? src[width] : c) >> 1 & 0x1;
245             dst[(width+offset)*8+5] = (usePtr ? src[width] : c) >> 2 & 0x1;
246             dst[(width+offset)*8+4] = (usePtr ? src[width] : c) >> 3 & 0x1;
247             dst[(width+offset)*8+3] = (usePtr ? src[width] : c) >> 4 & 0x1;
248             dst[(width+offset)*8+2] = (usePtr ? src[width] : c) >> 5 & 0x1;
249             dst[(width+offset)*8+1] = (usePtr ? src[width] : c) >> 6 & 0x1;
250             dst[(width+offset)*8+0] = (usePtr ? src[width] : c) >> 7;
251         }
252         break;
253     case 2:
254         while (--width >= 0) {
255             dst[(width+offset)*4+3] = (usePtr ? src[width] : c) & 0x3;
256             dst[(width+offset)*4+2] = (usePtr ? src[width] : c) >> 2 & 0x3;
257             dst[(width+offset)*4+1] = (usePtr ? src[width] : c) >> 4 & 0x3;
258             dst[(width+offset)*4+0] = (usePtr ? src[width] : c) >> 6;
259         }
260         break;
261     case 4:
262         while (--width >= 0) {
263             dst[(width+offset)*2+1] = (usePtr ? src[width] : c) & 0xF;
264             dst[(width+offset)*2+0] = (usePtr ? src[width] : c) >> 4;
265         }
266         break;
267     default:
268         if (usePtr) {
269             memcpy(dst + offset, src, width);
270         } else {
271             memset(dst + offset, c, width);
272         }
273     }
274 }
275
276 static int deinvert_buffer(TiffContext *s, const uint8_t *src, int size)
277 {
278     int i;
279
280     av_fast_padded_malloc(&s->deinvert_buf, &s->deinvert_buf_size, size);
281     if (!s->deinvert_buf)
282         return AVERROR(ENOMEM);
283     for (i = 0; i < size; i++)
284         s->deinvert_buf[i] = ff_reverse[src[i]];
285
286     return 0;
287 }
288
289 static void unpack_yuv(TiffContext *s, AVFrame *p,
290                        const uint8_t *src, int lnum)
291 {
292     int i, j, k;
293     int w       = (s->width - 1) / s->subsampling[0] + 1;
294     uint8_t *pu = &p->data[1][lnum / s->subsampling[1] * p->linesize[1]];
295     uint8_t *pv = &p->data[2][lnum / s->subsampling[1] * p->linesize[2]];
296     if (s->width % s->subsampling[0] || s->height % s->subsampling[1]) {
297         for (i = 0; i < w; i++) {
298             for (j = 0; j < s->subsampling[1]; j++)
299                 for (k = 0; k < s->subsampling[0]; k++)
300                     p->data[0][FFMIN(lnum + j, s->height-1) * p->linesize[0] +
301                                FFMIN(i * s->subsampling[0] + k, s->width-1)] = *src++;
302             *pu++ = *src++;
303             *pv++ = *src++;
304         }
305     }else{
306         for (i = 0; i < w; i++) {
307             for (j = 0; j < s->subsampling[1]; j++)
308                 for (k = 0; k < s->subsampling[0]; k++)
309                     p->data[0][(lnum + j) * p->linesize[0] +
310                                i * s->subsampling[0] + k] = *src++;
311             *pu++ = *src++;
312             *pv++ = *src++;
313         }
314     }
315 }
316
317 #if CONFIG_ZLIB
318 static int tiff_uncompress(uint8_t *dst, unsigned long *len, const uint8_t *src,
319                            int size)
320 {
321     z_stream zstream = { 0 };
322     int zret;
323
324     zstream.next_in   = (uint8_t *)src;
325     zstream.avail_in  = size;
326     zstream.next_out  = dst;
327     zstream.avail_out = *len;
328     zret              = inflateInit(&zstream);
329     if (zret != Z_OK) {
330         av_log(NULL, AV_LOG_ERROR, "Inflate init error: %d\n", zret);
331         return zret;
332     }
333     zret = inflate(&zstream, Z_SYNC_FLUSH);
334     inflateEnd(&zstream);
335     *len = zstream.total_out;
336     return zret == Z_STREAM_END ? Z_OK : zret;
337 }
338
339 static int tiff_unpack_zlib(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
340                             const uint8_t *src, int size, int width, int lines,
341                             int strip_start, int is_yuv)
342 {
343     uint8_t *zbuf;
344     unsigned long outlen;
345     int ret, line;
346     outlen = width * lines;
347     zbuf   = av_malloc(outlen);
348     if (!zbuf)
349         return AVERROR(ENOMEM);
350     if (s->fill_order) {
351         if ((ret = deinvert_buffer(s, src, size)) < 0) {
352             av_free(zbuf);
353             return ret;
354         }
355         src = s->deinvert_buf;
356     }
357     ret = tiff_uncompress(zbuf, &outlen, src, size);
358     if (ret != Z_OK) {
359         av_log(s->avctx, AV_LOG_ERROR,
360                "Uncompressing failed (%lu of %lu) with error %d\n", outlen,
361                (unsigned long)width * lines, ret);
362         av_free(zbuf);
363         return AVERROR_UNKNOWN;
364     }
365     src = zbuf;
366     for (line = 0; line < lines; line++) {
367         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
368             horizontal_fill(s->bpp, dst, 1, src, 0, width, 0);
369         } else {
370             memcpy(dst, src, width);
371         }
372         if (is_yuv) {
373             unpack_yuv(s, p, dst, strip_start + line);
374             line += s->subsampling[1] - 1;
375         }
376         dst += stride;
377         src += width;
378     }
379     av_free(zbuf);
380     return 0;
381 }
382 #endif
383
384 #if CONFIG_LZMA
385 static int tiff_uncompress_lzma(uint8_t *dst, uint64_t *len, const uint8_t *src,
386                                 int size)
387 {
388     lzma_stream stream = LZMA_STREAM_INIT;
389     lzma_ret ret;
390
391     stream.next_in   = (uint8_t *)src;
392     stream.avail_in  = size;
393     stream.next_out  = dst;
394     stream.avail_out = *len;
395     ret              = lzma_stream_decoder(&stream, UINT64_MAX, 0);
396     if (ret != LZMA_OK) {
397         av_log(NULL, AV_LOG_ERROR, "LZMA init error: %d\n", ret);
398         return ret;
399     }
400     ret = lzma_code(&stream, LZMA_RUN);
401     lzma_end(&stream);
402     *len = stream.total_out;
403     return ret == LZMA_STREAM_END ? LZMA_OK : ret;
404 }
405
406 static int tiff_unpack_lzma(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
407                             const uint8_t *src, int size, int width, int lines,
408                             int strip_start, int is_yuv)
409 {
410     uint64_t outlen = width * lines;
411     int ret, line;
412     uint8_t *buf = av_malloc(outlen);
413     if (!buf)
414         return AVERROR(ENOMEM);
415     if (s->fill_order) {
416         if ((ret = deinvert_buffer(s, src, size)) < 0) {
417             av_free(buf);
418             return ret;
419         }
420         src = s->deinvert_buf;
421     }
422     ret = tiff_uncompress_lzma(buf, &outlen, src, size);
423     if (ret != LZMA_OK) {
424         av_log(s->avctx, AV_LOG_ERROR,
425                "Uncompressing failed (%"PRIu64" of %"PRIu64") with error %d\n", outlen,
426                (uint64_t)width * lines, ret);
427         av_free(buf);
428         return AVERROR_UNKNOWN;
429     }
430     src = buf;
431     for (line = 0; line < lines; line++) {
432         if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
433             horizontal_fill(s->bpp, dst, 1, src, 0, width, 0);
434         } else {
435             memcpy(dst, src, width);
436         }
437         if (is_yuv) {
438             unpack_yuv(s, p, dst, strip_start + line);
439             line += s->subsampling[1] - 1;
440         }
441         dst += stride;
442         src += width;
443     }
444     av_free(buf);
445     return 0;
446 }
447 #endif
448
449 static int tiff_unpack_fax(TiffContext *s, uint8_t *dst, int stride,
450                            const uint8_t *src, int size, int width, int lines)
451 {
452     int i, ret = 0;
453     int line;
454     uint8_t *src2 = av_malloc((unsigned)size +
455                               FF_INPUT_BUFFER_PADDING_SIZE);
456
457     if (!src2) {
458         av_log(s->avctx, AV_LOG_ERROR,
459                "Error allocating temporary buffer\n");
460         return AVERROR(ENOMEM);
461     }
462     if (s->fax_opts & 2) {
463         avpriv_request_sample(s->avctx, "Uncompressed fax mode");
464         av_free(src2);
465         return AVERROR_PATCHWELCOME;
466     }
467     if (!s->fill_order) {
468         memcpy(src2, src, size);
469     } else {
470         for (i = 0; i < size; i++)
471             src2[i] = ff_reverse[src[i]];
472     }
473     memset(src2 + size, 0, FF_INPUT_BUFFER_PADDING_SIZE);
474     ret = ff_ccitt_unpack(s->avctx, src2, size, dst, lines, stride,
475                           s->compr, s->fax_opts);
476     if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
477         for (line = 0; line < lines; line++) {
478             horizontal_fill(s->bpp, dst, 1, dst, 0, width, 0);
479             dst += stride;
480         }
481     av_free(src2);
482     return ret;
483 }
484
485 static int tiff_unpack_strip(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
486                              const uint8_t *src, int size, int strip_start, int lines)
487 {
488     PutByteContext pb;
489     int c, line, pixels, code, ret;
490     const uint8_t *ssrc = src;
491     int width = ((s->width * s->bpp) + 7) >> 3;
492     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(p->format);
493     int is_yuv = !(desc->flags & AV_PIX_FMT_FLAG_RGB) && desc->nb_components >= 2;
494
495     if (s->planar)
496         width /= s->bppcount;
497
498     if (size <= 0)
499         return AVERROR_INVALIDDATA;
500
501     if (is_yuv) {
502         int bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
503                             s->subsampling[0] * s->subsampling[1] + 7) >> 3;
504         av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
505         if (s->yuv_line == NULL) {
506             av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
507             return AVERROR(ENOMEM);
508         }
509         dst = s->yuv_line;
510         stride = 0;
511         width = s->width * s->subsampling[1] + 2*(s->width / s->subsampling[0]);
512         av_assert0(width <= bytes_per_row);
513         av_assert0(s->bpp == 24);
514     }
515
516     if (s->compr == TIFF_DEFLATE || s->compr == TIFF_ADOBE_DEFLATE) {
517 #if CONFIG_ZLIB
518         return tiff_unpack_zlib(s, p, dst, stride, src, size, width, lines,
519                                 strip_start, is_yuv);
520 #else
521         av_log(s->avctx, AV_LOG_ERROR,
522                "zlib support not enabled, "
523                "deflate compression not supported\n");
524         return AVERROR(ENOSYS);
525 #endif
526     }
527     if (s->compr == TIFF_LZMA) {
528 #if CONFIG_LZMA
529         return tiff_unpack_lzma(s, p, dst, stride, src, size, width, lines,
530                                 strip_start, is_yuv);
531 #else
532         av_log(s->avctx, AV_LOG_ERROR,
533                "LZMA support not enabled\n");
534         return AVERROR(ENOSYS);
535 #endif
536     }
537     if (s->compr == TIFF_LZW) {
538         if (s->fill_order) {
539             if ((ret = deinvert_buffer(s, src, size)) < 0)
540                 return ret;
541             ssrc = src = s->deinvert_buf;
542         }
543         if (size > 1 && !src[0] && (src[1]&1)) {
544             av_log(s->avctx, AV_LOG_ERROR, "Old style LZW is unsupported\n");
545         }
546         if ((ret = ff_lzw_decode_init(s->lzw, 8, src, size, FF_LZW_TIFF)) < 0) {
547             av_log(s->avctx, AV_LOG_ERROR, "Error initializing LZW decoder\n");
548             return ret;
549         }
550         for (line = 0; line < lines; line++) {
551             pixels = ff_lzw_decode(s->lzw, dst, width);
552             if (pixels < width) {
553                 av_log(s->avctx, AV_LOG_ERROR, "Decoded only %i bytes of %i\n",
554                        pixels, width);
555                 return AVERROR_INVALIDDATA;
556             }
557             if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
558                 horizontal_fill(s->bpp, dst, 1, dst, 0, width, 0);
559             if (is_yuv) {
560                 unpack_yuv(s, p, dst, strip_start + line);
561                 line += s->subsampling[1] - 1;
562             }
563             dst += stride;
564         }
565         return 0;
566     }
567     if (s->compr == TIFF_CCITT_RLE ||
568         s->compr == TIFF_G3        ||
569         s->compr == TIFF_G4) {
570         if (is_yuv)
571             return AVERROR_INVALIDDATA;
572
573         return tiff_unpack_fax(s, dst, stride, src, size, width, lines);
574     }
575
576     bytestream2_init(&s->gb, src, size);
577     bytestream2_init_writer(&pb, dst, is_yuv ? s->yuv_line_size : (stride * lines));
578
579     for (line = 0; line < lines; line++) {
580         if (src - ssrc > size) {
581             av_log(s->avctx, AV_LOG_ERROR, "Source data overread\n");
582             return AVERROR_INVALIDDATA;
583         }
584
585         if (bytestream2_get_bytes_left(&s->gb) == 0 || bytestream2_get_eof(&pb))
586             break;
587         bytestream2_seek_p(&pb, stride * line, SEEK_SET);
588         switch (s->compr) {
589         case TIFF_RAW:
590             if (ssrc + size - src < width)
591                 return AVERROR_INVALIDDATA;
592
593             if (!s->fill_order) {
594                 horizontal_fill(s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
595                                 dst, 1, src, 0, width, 0);
596             } else {
597                 int i;
598                 for (i = 0; i < width; i++)
599                     dst[i] = ff_reverse[src[i]];
600             }
601             src += width;
602             break;
603         case TIFF_PACKBITS:
604             for (pixels = 0; pixels < width;) {
605                 if (ssrc + size - src < 2) {
606                     av_log(s->avctx, AV_LOG_ERROR, "Read went out of bounds\n");
607                     return AVERROR_INVALIDDATA;
608                 }
609                 code = s->fill_order ? (int8_t) ff_reverse[*src++]: (int8_t) *src++;
610                 if (code >= 0) {
611                     code++;
612                     if (pixels + code > width ||
613                         ssrc + size - src < code) {
614                         av_log(s->avctx, AV_LOG_ERROR,
615                                "Copy went out of bounds\n");
616                         return AVERROR_INVALIDDATA;
617                     }
618                     horizontal_fill(s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
619                                     dst, 1, src, 0, code, pixels);
620                     src    += code;
621                     pixels += code;
622                 } else if (code != -128) { // -127..-1
623                     code = (-code) + 1;
624                     if (pixels + code > width) {
625                         av_log(s->avctx, AV_LOG_ERROR,
626                                "Run went out of bounds\n");
627                         return AVERROR_INVALIDDATA;
628                     }
629                     c = *src++;
630                     horizontal_fill(s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
631                                     dst, 0, NULL, c, code, pixels);
632                     pixels += code;
633                 }
634             }
635             if (s->fill_order) {
636                 int i;
637                 for (i = 0; i < width; i++)
638                     dst[i] = ff_reverse[dst[i]];
639             }
640             break;
641         }
642         if (is_yuv) {
643             unpack_yuv(s, p, dst, strip_start + line);
644             line += s->subsampling[1] - 1;
645         }
646         dst += stride;
647     }
648     return 0;
649 }
650
651 static int init_image(TiffContext *s, ThreadFrame *frame)
652 {
653     int ret;
654
655     switch (s->planar * 1000 + s->bpp * 10 + s->bppcount) {
656     case 11:
657         if (!s->palette_is_set) {
658             s->avctx->pix_fmt = AV_PIX_FMT_MONOBLACK;
659             break;
660         }
661     case 21:
662     case 41:
663     case 81:
664         s->avctx->pix_fmt = s->palette_is_set ? AV_PIX_FMT_PAL8 : AV_PIX_FMT_GRAY8;
665         break;
666     case 243:
667         if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
668             if (s->subsampling[0] == 1 && s->subsampling[1] == 1) {
669                 s->avctx->pix_fmt = AV_PIX_FMT_YUV444P;
670             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 1) {
671                 s->avctx->pix_fmt = AV_PIX_FMT_YUV422P;
672             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 1) {
673                 s->avctx->pix_fmt = AV_PIX_FMT_YUV411P;
674             } else if (s->subsampling[0] == 1 && s->subsampling[1] == 2) {
675                 s->avctx->pix_fmt = AV_PIX_FMT_YUV440P;
676             } else if (s->subsampling[0] == 2 && s->subsampling[1] == 2) {
677                 s->avctx->pix_fmt = AV_PIX_FMT_YUV420P;
678             } else if (s->subsampling[0] == 4 && s->subsampling[1] == 4) {
679                 s->avctx->pix_fmt = AV_PIX_FMT_YUV410P;
680             } else {
681                 av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr subsampling\n");
682                 return AVERROR_PATCHWELCOME;
683             }
684         } else
685             s->avctx->pix_fmt = AV_PIX_FMT_RGB24;
686         break;
687     case 161:
688         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GRAY16LE : AV_PIX_FMT_GRAY16BE;
689         break;
690     case 162:
691         s->avctx->pix_fmt = AV_PIX_FMT_YA8;
692         break;
693     case 322:
694         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_YA16LE : AV_PIX_FMT_YA16BE;
695         break;
696     case 324:
697         s->avctx->pix_fmt = AV_PIX_FMT_RGBA;
698         break;
699     case 483:
700         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGB48LE  : AV_PIX_FMT_RGB48BE;
701         break;
702     case 644:
703         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGBA64LE  : AV_PIX_FMT_RGBA64BE;
704         break;
705     case 1243:
706         s->avctx->pix_fmt = AV_PIX_FMT_GBRP;
707         break;
708     case 1324:
709         s->avctx->pix_fmt = AV_PIX_FMT_GBRAP;
710         break;
711     case 1483:
712         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRP16LE : AV_PIX_FMT_GBRP16BE;
713         break;
714     case 1644:
715         s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRAP16LE : AV_PIX_FMT_GBRAP16BE;
716         break;
717     default:
718         av_log(s->avctx, AV_LOG_ERROR,
719                "This format is not supported (bpp=%d, bppcount=%d)\n",
720                s->bpp, s->bppcount);
721         return AVERROR_INVALIDDATA;
722     }
723
724     if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
725         const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(s->avctx->pix_fmt);
726         if((desc->flags & AV_PIX_FMT_FLAG_RGB) || desc->nb_components < 3) {
727             av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr variant\n");
728             return AVERROR_INVALIDDATA;
729         }
730     }
731
732     if (s->width != s->avctx->width || s->height != s->avctx->height) {
733         ret = ff_set_dimensions(s->avctx, s->width, s->height);
734         if (ret < 0)
735             return ret;
736     }
737     if ((ret = ff_thread_get_buffer(s->avctx, frame, 0)) < 0)
738         return ret;
739     if (s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
740         memcpy(frame->f->data[1], s->palette, sizeof(s->palette));
741     }
742     return 0;
743 }
744
745 static void set_sar(TiffContext *s, unsigned tag, unsigned num, unsigned den)
746 {
747     int offset = tag == TIFF_YRES ? 2 : 0;
748     s->res[offset++] = num;
749     s->res[offset]   = den;
750     if (s->res[0] && s->res[1] && s->res[2] && s->res[3])
751         av_reduce(&s->avctx->sample_aspect_ratio.num, &s->avctx->sample_aspect_ratio.den,
752                   s->res[2] * (uint64_t)s->res[1], s->res[0] * (uint64_t)s->res[3], INT32_MAX);
753 }
754
755 static int tiff_decode_tag(TiffContext *s, AVFrame *frame)
756 {
757     unsigned tag, type, count, off, value = 0, value2 = 0;
758     int i, start;
759     int pos;
760     int ret;
761     double *dp;
762
763     ret = ff_tread_tag(&s->gb, s->le, &tag, &type, &count, &start);
764     if (ret < 0) {
765         goto end;
766     }
767
768     off = bytestream2_tell(&s->gb);
769     if (count == 1) {
770         switch (type) {
771         case TIFF_BYTE:
772         case TIFF_SHORT:
773         case TIFF_LONG:
774             value = ff_tget(&s->gb, type, s->le);
775             break;
776         case TIFF_RATIONAL:
777             value  = ff_tget(&s->gb, TIFF_LONG, s->le);
778             value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
779             break;
780         case TIFF_STRING:
781             if (count <= 4) {
782                 break;
783             }
784         default:
785             value = UINT_MAX;
786         }
787     }
788
789     switch (tag) {
790     case TIFF_WIDTH:
791         s->width = value;
792         break;
793     case TIFF_HEIGHT:
794         s->height = value;
795         break;
796     case TIFF_BPP:
797         s->bppcount = count;
798         if (count > 4) {
799             av_log(s->avctx, AV_LOG_ERROR,
800                    "This format is not supported (bpp=%d, %d components)\n",
801                    s->bpp, count);
802             return AVERROR_INVALIDDATA;
803         }
804         if (count == 1)
805             s->bpp = value;
806         else {
807             switch (type) {
808             case TIFF_BYTE:
809             case TIFF_SHORT:
810             case TIFF_LONG:
811                 s->bpp = 0;
812                 if (bytestream2_get_bytes_left(&s->gb) < type_sizes[type] * count)
813                     return AVERROR_INVALIDDATA;
814                 for (i = 0; i < count; i++)
815                     s->bpp += ff_tget(&s->gb, type, s->le);
816                 break;
817             default:
818                 s->bpp = -1;
819             }
820         }
821         break;
822     case TIFF_SAMPLES_PER_PIXEL:
823         if (count != 1) {
824             av_log(s->avctx, AV_LOG_ERROR,
825                    "Samples per pixel requires a single value, many provided\n");
826             return AVERROR_INVALIDDATA;
827         }
828         if (value > 4U) {
829             av_log(s->avctx, AV_LOG_ERROR,
830                    "Samples per pixel %d is too large\n", value);
831             return AVERROR_INVALIDDATA;
832         }
833         if (s->bppcount == 1)
834             s->bpp *= value;
835         s->bppcount = value;
836         break;
837     case TIFF_COMPR:
838         s->compr     = value;
839         s->predictor = 0;
840         switch (s->compr) {
841         case TIFF_RAW:
842         case TIFF_PACKBITS:
843         case TIFF_LZW:
844         case TIFF_CCITT_RLE:
845             break;
846         case TIFF_G3:
847         case TIFF_G4:
848             s->fax_opts = 0;
849             break;
850         case TIFF_DEFLATE:
851         case TIFF_ADOBE_DEFLATE:
852 #if CONFIG_ZLIB
853             break;
854 #else
855             av_log(s->avctx, AV_LOG_ERROR, "Deflate: ZLib not compiled in\n");
856             return AVERROR(ENOSYS);
857 #endif
858         case TIFF_JPEG:
859         case TIFF_NEWJPEG:
860             avpriv_report_missing_feature(s->avctx, "JPEG compression");
861             return AVERROR_PATCHWELCOME;
862         case TIFF_LZMA:
863 #if CONFIG_LZMA
864             break;
865 #else
866             av_log(s->avctx, AV_LOG_ERROR, "LZMA not compiled in\n");
867             return AVERROR(ENOSYS);
868 #endif
869         default:
870             av_log(s->avctx, AV_LOG_ERROR, "Unknown compression method %i\n",
871                    s->compr);
872             return AVERROR_INVALIDDATA;
873         }
874         break;
875     case TIFF_ROWSPERSTRIP:
876         if (!value || (type == TIFF_LONG && value == UINT_MAX))
877             value = s->height;
878         s->rps = FFMIN(value, s->height);
879         break;
880     case TIFF_STRIP_OFFS:
881         if (count == 1) {
882             s->strippos = 0;
883             s->stripoff = value;
884         } else
885             s->strippos = off;
886         s->strips = count;
887         if (s->strips == 1)
888             s->rps = s->height;
889         s->sot = type;
890         break;
891     case TIFF_STRIP_SIZE:
892         if (count == 1) {
893             s->stripsizesoff = 0;
894             s->stripsize     = value;
895             s->strips        = 1;
896         } else {
897             s->stripsizesoff = off;
898         }
899         s->strips = count;
900         s->sstype = type;
901         break;
902     case TIFF_XRES:
903     case TIFF_YRES:
904         set_sar(s, tag, value, value2);
905         break;
906     case TIFF_TILE_BYTE_COUNTS:
907     case TIFF_TILE_LENGTH:
908     case TIFF_TILE_OFFSETS:
909     case TIFF_TILE_WIDTH:
910         av_log(s->avctx, AV_LOG_ERROR, "Tiled images are not supported\n");
911         return AVERROR_PATCHWELCOME;
912         break;
913     case TIFF_PREDICTOR:
914         s->predictor = value;
915         break;
916     case TIFF_PHOTOMETRIC:
917         switch (value) {
918         case TIFF_PHOTOMETRIC_WHITE_IS_ZERO:
919         case TIFF_PHOTOMETRIC_BLACK_IS_ZERO:
920         case TIFF_PHOTOMETRIC_RGB:
921         case TIFF_PHOTOMETRIC_PALETTE:
922         case TIFF_PHOTOMETRIC_YCBCR:
923             s->photometric = value;
924             break;
925         case TIFF_PHOTOMETRIC_ALPHA_MASK:
926         case TIFF_PHOTOMETRIC_SEPARATED:
927         case TIFF_PHOTOMETRIC_CIE_LAB:
928         case TIFF_PHOTOMETRIC_ICC_LAB:
929         case TIFF_PHOTOMETRIC_ITU_LAB:
930         case TIFF_PHOTOMETRIC_CFA:
931         case TIFF_PHOTOMETRIC_LOG_L:
932         case TIFF_PHOTOMETRIC_LOG_LUV:
933         case TIFF_PHOTOMETRIC_LINEAR_RAW:
934             avpriv_report_missing_feature(s->avctx,
935                                           "PhotometricInterpretation 0x%04X",
936                                           value);
937             return AVERROR_PATCHWELCOME;
938         default:
939             av_log(s->avctx, AV_LOG_ERROR, "PhotometricInterpretation %u is "
940                    "unknown\n", value);
941             return AVERROR_INVALIDDATA;
942         }
943         break;
944     case TIFF_FILL_ORDER:
945         if (value < 1 || value > 2) {
946             av_log(s->avctx, AV_LOG_ERROR,
947                    "Unknown FillOrder value %d, trying default one\n", value);
948             value = 1;
949         }
950         s->fill_order = value - 1;
951         break;
952     case TIFF_PAL: {
953         GetByteContext pal_gb[3];
954         off = type_sizes[type];
955         if (count / 3 > 256 ||
956             bytestream2_get_bytes_left(&s->gb) < count / 3 * off * 3)
957             return AVERROR_INVALIDDATA;
958
959         pal_gb[0] = pal_gb[1] = pal_gb[2] = s->gb;
960         bytestream2_skip(&pal_gb[1], count / 3 * off);
961         bytestream2_skip(&pal_gb[2], count / 3 * off * 2);
962
963         off = (type_sizes[type] - 1) << 3;
964         for (i = 0; i < count / 3; i++) {
965             uint32_t p = 0xFF000000;
966             p |= (ff_tget(&pal_gb[0], type, s->le) >> off) << 16;
967             p |= (ff_tget(&pal_gb[1], type, s->le) >> off) << 8;
968             p |=  ff_tget(&pal_gb[2], type, s->le) >> off;
969             s->palette[i] = p;
970         }
971         s->palette_is_set = 1;
972         break;
973     }
974     case TIFF_PLANAR:
975         s->planar = value == 2;
976         break;
977     case TIFF_YCBCR_SUBSAMPLING:
978         if (count != 2) {
979             av_log(s->avctx, AV_LOG_ERROR, "subsample count invalid\n");
980             return AVERROR_INVALIDDATA;
981         }
982         for (i = 0; i < count; i++)
983             s->subsampling[i] = ff_tget(&s->gb, type, s->le);
984         break;
985     case TIFF_T4OPTIONS:
986         if (s->compr == TIFF_G3)
987             s->fax_opts = value;
988         break;
989     case TIFF_T6OPTIONS:
990         if (s->compr == TIFF_G4)
991             s->fax_opts = value;
992         break;
993 #define ADD_METADATA(count, name, sep)\
994     if ((ret = add_metadata(count, type, name, sep, s, frame)) < 0) {\
995         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");\
996         goto end;\
997     }
998     case TIFF_MODEL_PIXEL_SCALE:
999         ADD_METADATA(count, "ModelPixelScaleTag", NULL);
1000         break;
1001     case TIFF_MODEL_TRANSFORMATION:
1002         ADD_METADATA(count, "ModelTransformationTag", NULL);
1003         break;
1004     case TIFF_MODEL_TIEPOINT:
1005         ADD_METADATA(count, "ModelTiepointTag", NULL);
1006         break;
1007     case TIFF_GEO_KEY_DIRECTORY:
1008         ADD_METADATA(1, "GeoTIFF_Version", NULL);
1009         ADD_METADATA(2, "GeoTIFF_Key_Revision", ".");
1010         s->geotag_count   = ff_tget_short(&s->gb, s->le);
1011         if (s->geotag_count > count / 4 - 1) {
1012             s->geotag_count = count / 4 - 1;
1013             av_log(s->avctx, AV_LOG_WARNING, "GeoTIFF key directory buffer shorter than specified\n");
1014         }
1015         if (bytestream2_get_bytes_left(&s->gb) < s->geotag_count * sizeof(int16_t) * 4) {
1016             s->geotag_count = 0;
1017             return -1;
1018         }
1019         s->geotags = av_mallocz_array(s->geotag_count, sizeof(TiffGeoTag));
1020         if (!s->geotags) {
1021             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1022             s->geotag_count = 0;
1023             goto end;
1024         }
1025         for (i = 0; i < s->geotag_count; i++) {
1026             s->geotags[i].key    = ff_tget_short(&s->gb, s->le);
1027             s->geotags[i].type   = ff_tget_short(&s->gb, s->le);
1028             s->geotags[i].count  = ff_tget_short(&s->gb, s->le);
1029
1030             if (!s->geotags[i].type)
1031                 s->geotags[i].val  = get_geokey_val(s->geotags[i].key, ff_tget_short(&s->gb, s->le));
1032             else
1033                 s->geotags[i].offset = ff_tget_short(&s->gb, s->le);
1034         }
1035         break;
1036     case TIFF_GEO_DOUBLE_PARAMS:
1037         if (count >= INT_MAX / sizeof(int64_t))
1038             return AVERROR_INVALIDDATA;
1039         if (bytestream2_get_bytes_left(&s->gb) < count * sizeof(int64_t))
1040             return AVERROR_INVALIDDATA;
1041         dp = av_malloc_array(count, sizeof(double));
1042         if (!dp) {
1043             av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1044             goto end;
1045         }
1046         for (i = 0; i < count; i++)
1047             dp[i] = ff_tget_double(&s->gb, s->le);
1048         for (i = 0; i < s->geotag_count; i++) {
1049             if (s->geotags[i].type == TIFF_GEO_DOUBLE_PARAMS) {
1050                 if (s->geotags[i].count == 0
1051                     || s->geotags[i].offset + s->geotags[i].count > count) {
1052                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1053                 } else {
1054                     char *ap = doubles2str(&dp[s->geotags[i].offset], s->geotags[i].count, ", ");
1055                     if (!ap) {
1056                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1057                         av_freep(&dp);
1058                         return AVERROR(ENOMEM);
1059                     }
1060                     s->geotags[i].val = ap;
1061                 }
1062             }
1063         }
1064         av_freep(&dp);
1065         break;
1066     case TIFF_GEO_ASCII_PARAMS:
1067         pos = bytestream2_tell(&s->gb);
1068         for (i = 0; i < s->geotag_count; i++) {
1069             if (s->geotags[i].type == TIFF_GEO_ASCII_PARAMS) {
1070                 if (s->geotags[i].count == 0
1071                     || s->geotags[i].offset +  s->geotags[i].count > count) {
1072                     av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1073                 } else {
1074                     char *ap;
1075
1076                     bytestream2_seek(&s->gb, pos + s->geotags[i].offset, SEEK_SET);
1077                     if (bytestream2_get_bytes_left(&s->gb) < s->geotags[i].count)
1078                         return AVERROR_INVALIDDATA;
1079                     ap = av_malloc(s->geotags[i].count);
1080                     if (!ap) {
1081                         av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1082                         return AVERROR(ENOMEM);
1083                     }
1084                     bytestream2_get_bufferu(&s->gb, ap, s->geotags[i].count);
1085                     ap[s->geotags[i].count - 1] = '\0'; //replace the "|" delimiter with a 0 byte
1086                     s->geotags[i].val = ap;
1087                 }
1088             }
1089         }
1090         break;
1091     case TIFF_ARTIST:
1092         ADD_METADATA(count, "artist", NULL);
1093         break;
1094     case TIFF_COPYRIGHT:
1095         ADD_METADATA(count, "copyright", NULL);
1096         break;
1097     case TIFF_DATE:
1098         ADD_METADATA(count, "date", NULL);
1099         break;
1100     case TIFF_DOCUMENT_NAME:
1101         ADD_METADATA(count, "document_name", NULL);
1102         break;
1103     case TIFF_HOST_COMPUTER:
1104         ADD_METADATA(count, "computer", NULL);
1105         break;
1106     case TIFF_IMAGE_DESCRIPTION:
1107         ADD_METADATA(count, "description", NULL);
1108         break;
1109     case TIFF_MAKE:
1110         ADD_METADATA(count, "make", NULL);
1111         break;
1112     case TIFF_MODEL:
1113         ADD_METADATA(count, "model", NULL);
1114         break;
1115     case TIFF_PAGE_NAME:
1116         ADD_METADATA(count, "page_name", NULL);
1117         break;
1118     case TIFF_PAGE_NUMBER:
1119         ADD_METADATA(count, "page_number", " / ");
1120         break;
1121     case TIFF_SOFTWARE_NAME:
1122         ADD_METADATA(count, "software", NULL);
1123         break;
1124     default:
1125         if (s->avctx->err_recognition & AV_EF_EXPLODE) {
1126             av_log(s->avctx, AV_LOG_ERROR,
1127                    "Unknown or unsupported tag %d/0X%0X\n",
1128                    tag, tag);
1129             return AVERROR_INVALIDDATA;
1130         }
1131     }
1132 end:
1133     bytestream2_seek(&s->gb, start, SEEK_SET);
1134     return 0;
1135 }
1136
1137 static int decode_frame(AVCodecContext *avctx,
1138                         void *data, int *got_frame, AVPacket *avpkt)
1139 {
1140     TiffContext *const s = avctx->priv_data;
1141     AVFrame *const p = data;
1142     ThreadFrame frame = { .f = data };
1143     unsigned off;
1144     int le, ret, plane, planes;
1145     int i, j, entries, stride;
1146     unsigned soff, ssize;
1147     uint8_t *dst;
1148     GetByteContext stripsizes;
1149     GetByteContext stripdata;
1150
1151     bytestream2_init(&s->gb, avpkt->data, avpkt->size);
1152
1153     // parse image header
1154     if ((ret = ff_tdecode_header(&s->gb, &le, &off))) {
1155         av_log(avctx, AV_LOG_ERROR, "Invalid TIFF header\n");
1156         return ret;
1157     } else if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1158         av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1159         return AVERROR_INVALIDDATA;
1160     }
1161     s->le          = le;
1162     // TIFF_BPP is not a required tag and defaults to 1
1163     s->bppcount    = s->bpp = 1;
1164     s->photometric = TIFF_PHOTOMETRIC_NONE;
1165     s->compr       = TIFF_RAW;
1166     s->fill_order  = 0;
1167     free_geotags(s);
1168
1169     // Reset these offsets so we can tell if they were set this frame
1170     s->stripsizesoff = s->strippos = 0;
1171     /* parse image file directory */
1172     bytestream2_seek(&s->gb, off, SEEK_SET);
1173     entries = ff_tget_short(&s->gb, le);
1174     if (bytestream2_get_bytes_left(&s->gb) < entries * 12)
1175         return AVERROR_INVALIDDATA;
1176     for (i = 0; i < entries; i++) {
1177         if ((ret = tiff_decode_tag(s, p)) < 0)
1178             return ret;
1179     }
1180
1181     for (i = 0; i<s->geotag_count; i++) {
1182         const char *keyname = get_geokey_name(s->geotags[i].key);
1183         if (!keyname) {
1184             av_log(avctx, AV_LOG_WARNING, "Unknown or unsupported GeoTIFF key %d\n", s->geotags[i].key);
1185             continue;
1186         }
1187         if (get_geokey_type(s->geotags[i].key) != s->geotags[i].type) {
1188             av_log(avctx, AV_LOG_WARNING, "Type of GeoTIFF key %d is wrong\n", s->geotags[i].key);
1189             continue;
1190         }
1191         ret = av_dict_set(avpriv_frame_get_metadatap(p), keyname, s->geotags[i].val, 0);
1192         if (ret<0) {
1193             av_log(avctx, AV_LOG_ERROR, "Writing metadata with key '%s' failed\n", keyname);
1194             return ret;
1195         }
1196     }
1197
1198     if (!s->strippos && !s->stripoff) {
1199         av_log(avctx, AV_LOG_ERROR, "Image data is missing\n");
1200         return AVERROR_INVALIDDATA;
1201     }
1202     /* now we have the data and may start decoding */
1203     if ((ret = init_image(s, &frame)) < 0)
1204         return ret;
1205
1206     if (s->strips == 1 && !s->stripsize) {
1207         av_log(avctx, AV_LOG_WARNING, "Image data size missing\n");
1208         s->stripsize = avpkt->size - s->stripoff;
1209     }
1210
1211     if (s->stripsizesoff) {
1212         if (s->stripsizesoff >= (unsigned)avpkt->size)
1213             return AVERROR_INVALIDDATA;
1214         bytestream2_init(&stripsizes, avpkt->data + s->stripsizesoff,
1215                          avpkt->size - s->stripsizesoff);
1216     }
1217     if (s->strippos) {
1218         if (s->strippos >= (unsigned)avpkt->size)
1219             return AVERROR_INVALIDDATA;
1220         bytestream2_init(&stripdata, avpkt->data + s->strippos,
1221                          avpkt->size - s->strippos);
1222     }
1223
1224     if (s->rps <= 0) {
1225         av_log(avctx, AV_LOG_ERROR, "rps %d invalid\n", s->rps);
1226         return AVERROR_INVALIDDATA;
1227     }
1228
1229     planes = s->planar ? s->bppcount : 1;
1230     for (plane = 0; plane < planes; plane++) {
1231         stride = p->linesize[plane];
1232         dst    = p->data[plane];
1233     for (i = 0; i < s->height; i += s->rps) {
1234         if (s->stripsizesoff)
1235             ssize = ff_tget(&stripsizes, s->sstype, le);
1236         else
1237             ssize = s->stripsize;
1238
1239         if (s->strippos)
1240             soff = ff_tget(&stripdata, s->sot, le);
1241         else
1242             soff = s->stripoff;
1243
1244         if (soff > avpkt->size || ssize > avpkt->size - soff) {
1245             av_log(avctx, AV_LOG_ERROR, "Invalid strip size/offset\n");
1246             return AVERROR_INVALIDDATA;
1247         }
1248         if ((ret = tiff_unpack_strip(s, p, dst, stride, avpkt->data + soff, ssize, i,
1249                                      FFMIN(s->rps, s->height - i))) < 0) {
1250             if (avctx->err_recognition & AV_EF_EXPLODE)
1251                 return ret;
1252             break;
1253         }
1254         dst += s->rps * stride;
1255     }
1256     if (s->predictor == 2) {
1257         if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1258             av_log(s->avctx, AV_LOG_ERROR, "predictor == 2 with YUV is unsupported");
1259             return AVERROR_PATCHWELCOME;
1260         }
1261         dst   = p->data[plane];
1262         soff  = s->bpp >> 3;
1263         if (s->planar)
1264             soff  = FFMAX(soff / s->bppcount, 1);
1265         ssize = s->width * soff;
1266         if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48LE ||
1267             s->avctx->pix_fmt == AV_PIX_FMT_RGBA64LE ||
1268             s->avctx->pix_fmt == AV_PIX_FMT_GBRP16LE ||
1269             s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16LE) {
1270             for (i = 0; i < s->height; i++) {
1271                 for (j = soff; j < ssize; j += 2)
1272                     AV_WL16(dst + j, AV_RL16(dst + j) + AV_RL16(dst + j - soff));
1273                 dst += stride;
1274             }
1275         } else if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48BE ||
1276                    s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE ||
1277                    s->avctx->pix_fmt == AV_PIX_FMT_GBRP16BE ||
1278                    s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16BE) {
1279             for (i = 0; i < s->height; i++) {
1280                 for (j = soff; j < ssize; j += 2)
1281                     AV_WB16(dst + j, AV_RB16(dst + j) + AV_RB16(dst + j - soff));
1282                 dst += stride;
1283             }
1284         } else {
1285             for (i = 0; i < s->height; i++) {
1286                 for (j = soff; j < ssize; j++)
1287                     dst[j] += dst[j - soff];
1288                 dst += stride;
1289             }
1290         }
1291     }
1292
1293     if (s->photometric == TIFF_PHOTOMETRIC_WHITE_IS_ZERO) {
1294         dst = p->data[plane];
1295         for (i = 0; i < s->height; i++) {
1296             for (j = 0; j < p->linesize[plane]; j++)
1297                 dst[j] = (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 ? (1<<s->bpp) - 1 : 255) - dst[j];
1298             dst += stride;
1299         }
1300     }
1301     }
1302
1303     if (s->planar && s->bppcount > 2) {
1304         FFSWAP(uint8_t*, p->data[0],     p->data[2]);
1305         FFSWAP(int,      p->linesize[0], p->linesize[2]);
1306         FFSWAP(uint8_t*, p->data[0],     p->data[1]);
1307         FFSWAP(int,      p->linesize[0], p->linesize[1]);
1308     }
1309
1310     *got_frame = 1;
1311
1312     return avpkt->size;
1313 }
1314
1315 static av_cold int tiff_init(AVCodecContext *avctx)
1316 {
1317     TiffContext *s = avctx->priv_data;
1318
1319     s->width  = 0;
1320     s->height = 0;
1321     s->subsampling[0] =
1322     s->subsampling[1] = 1;
1323     s->avctx  = avctx;
1324     ff_lzw_decode_open(&s->lzw);
1325     ff_ccitt_unpack_init();
1326
1327     return 0;
1328 }
1329
1330 static av_cold int tiff_end(AVCodecContext *avctx)
1331 {
1332     TiffContext *const s = avctx->priv_data;
1333
1334     free_geotags(s);
1335
1336     ff_lzw_decode_close(&s->lzw);
1337     av_freep(&s->deinvert_buf);
1338     return 0;
1339 }
1340
1341 AVCodec ff_tiff_decoder = {
1342     .name           = "tiff",
1343     .long_name      = NULL_IF_CONFIG_SMALL("TIFF image"),
1344     .type           = AVMEDIA_TYPE_VIDEO,
1345     .id             = AV_CODEC_ID_TIFF,
1346     .priv_data_size = sizeof(TiffContext),
1347     .init           = tiff_init,
1348     .close          = tiff_end,
1349     .decode         = decode_frame,
1350     .init_thread_copy = ONLY_IF_THREADS_ENABLED(tiff_init),
1351     .capabilities   = CODEC_CAP_DR1 | CODEC_CAP_FRAME_THREADS,
1352 };