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