2 * Copyright (c) 2006 Konstantin Shishkov
4 * This file is part of FFmpeg.
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.
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.
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
24 * @author Konstantin Shishkov
32 #define LZMA_API_STATIC
36 #include "libavutil/attributes.h"
37 #include "libavutil/avstring.h"
38 #include "libavutil/error.h"
39 #include "libavutil/intreadwrite.h"
40 #include "libavutil/imgutils.h"
41 #include "libavutil/opt.h"
43 #include "bytestream.h"
49 #include "tiff_data.h"
54 typedef struct TiffContext {
56 AVCodecContext *avctx;
59 /* JPEG decoding for DNG */
60 AVCodecContext *avctx_mjpeg; // wrapper context for MJPEG
61 AVFrame *jpgframe; // decoded JPEG tile
67 enum TiffType tiff_type;
69 unsigned int bpp, bppcount;
70 uint32_t palette[256];
74 enum TiffPhotometric photometric;
87 const uint16_t *dng_lut; // Pointer to DNG linearization table
92 int strips, rps, sstype;
94 int stripsizesoff, stripsize, stripoff, strippos;
99 int tile_byte_counts_offset, tile_offsets_offset;
100 int tile_width, tile_length;
105 uint8_t *deinvert_buf;
106 int deinvert_buf_size;
108 unsigned int yuv_line_size;
110 unsigned int fax_buffer_size;
116 static void tiff_set_type(TiffContext *s, enum TiffType tiff_type) {
117 if (s->tiff_type < tiff_type) // Prioritize higher-valued entries
118 s->tiff_type = tiff_type;
121 static void free_geotags(TiffContext *const s)
124 for (i = 0; i < s->geotag_count; i++) {
125 if (s->geotags[i].val)
126 av_freep(&s->geotags[i].val);
128 av_freep(&s->geotags);
132 #define RET_GEOKEY(TYPE, array, element)\
133 if (key >= TIFF_##TYPE##_KEY_ID_OFFSET &&\
134 key - TIFF_##TYPE##_KEY_ID_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_name_type_map))\
135 return ff_tiff_##array##_name_type_map[key - TIFF_##TYPE##_KEY_ID_OFFSET].element;
137 static const char *get_geokey_name(int key)
139 RET_GEOKEY(VERT, vert, name);
140 RET_GEOKEY(PROJ, proj, name);
141 RET_GEOKEY(GEOG, geog, name);
142 RET_GEOKEY(CONF, conf, name);
147 static int get_geokey_type(int key)
149 RET_GEOKEY(VERT, vert, type);
150 RET_GEOKEY(PROJ, proj, type);
151 RET_GEOKEY(GEOG, geog, type);
152 RET_GEOKEY(CONF, conf, type);
154 return AVERROR_INVALIDDATA;
157 static int cmp_id_key(const void *id, const void *k)
159 return *(const int*)id - ((const TiffGeoTagKeyName*)k)->key;
162 static const char *search_keyval(const TiffGeoTagKeyName *keys, int n, int id)
164 TiffGeoTagKeyName *r = bsearch(&id, keys, n, sizeof(keys[0]), cmp_id_key);
171 static char *get_geokey_val(int key, int val)
175 if (val == TIFF_GEO_KEY_UNDEFINED)
176 return av_strdup("undefined");
177 if (val == TIFF_GEO_KEY_USER_DEFINED)
178 return av_strdup("User-Defined");
180 #define RET_GEOKEY_VAL(TYPE, array)\
181 if (val >= TIFF_##TYPE##_OFFSET &&\
182 val - TIFF_##TYPE##_OFFSET < FF_ARRAY_ELEMS(ff_tiff_##array##_codes))\
183 return av_strdup(ff_tiff_##array##_codes[val - TIFF_##TYPE##_OFFSET]);
186 case TIFF_GT_MODEL_TYPE_GEOKEY:
187 RET_GEOKEY_VAL(GT_MODEL_TYPE, gt_model_type);
189 case TIFF_GT_RASTER_TYPE_GEOKEY:
190 RET_GEOKEY_VAL(GT_RASTER_TYPE, gt_raster_type);
192 case TIFF_GEOG_LINEAR_UNITS_GEOKEY:
193 case TIFF_PROJ_LINEAR_UNITS_GEOKEY:
194 case TIFF_VERTICAL_UNITS_GEOKEY:
195 RET_GEOKEY_VAL(LINEAR_UNIT, linear_unit);
197 case TIFF_GEOG_ANGULAR_UNITS_GEOKEY:
198 case TIFF_GEOG_AZIMUTH_UNITS_GEOKEY:
199 RET_GEOKEY_VAL(ANGULAR_UNIT, angular_unit);
201 case TIFF_GEOGRAPHIC_TYPE_GEOKEY:
202 RET_GEOKEY_VAL(GCS_TYPE, gcs_type);
203 RET_GEOKEY_VAL(GCSE_TYPE, gcse_type);
205 case TIFF_GEOG_GEODETIC_DATUM_GEOKEY:
206 RET_GEOKEY_VAL(GEODETIC_DATUM, geodetic_datum);
207 RET_GEOKEY_VAL(GEODETIC_DATUM_E, geodetic_datum_e);
209 case TIFF_GEOG_ELLIPSOID_GEOKEY:
210 RET_GEOKEY_VAL(ELLIPSOID, ellipsoid);
212 case TIFF_GEOG_PRIME_MERIDIAN_GEOKEY:
213 RET_GEOKEY_VAL(PRIME_MERIDIAN, prime_meridian);
215 case TIFF_PROJECTED_CS_TYPE_GEOKEY:
216 ap = av_strdup(search_keyval(ff_tiff_proj_cs_type_codes, FF_ARRAY_ELEMS(ff_tiff_proj_cs_type_codes), val));
219 case TIFF_PROJECTION_GEOKEY:
220 ap = av_strdup(search_keyval(ff_tiff_projection_codes, FF_ARRAY_ELEMS(ff_tiff_projection_codes), val));
223 case TIFF_PROJ_COORD_TRANS_GEOKEY:
224 RET_GEOKEY_VAL(COORD_TRANS, coord_trans);
226 case TIFF_VERTICAL_CS_TYPE_GEOKEY:
227 RET_GEOKEY_VAL(VERT_CS, vert_cs);
228 RET_GEOKEY_VAL(ORTHO_VERT_CS, ortho_vert_cs);
235 snprintf(ap, 14, "Unknown-%d", val);
239 static char *doubles2str(double *dp, int count, const char *sep)
243 uint64_t component_len;
244 if (!sep) sep = ", ";
245 component_len = 24LL + strlen(sep);
246 if (count >= (INT_MAX - 1)/component_len)
248 ap = av_malloc(component_len * count + 1);
253 for (i = 0; i < count; i++) {
254 unsigned l = snprintf(ap, component_len, "%.15g%s", dp[i], sep);
255 if(l >= component_len) {
261 ap0[strlen(ap0) - strlen(sep)] = '\0';
265 static int add_metadata(int count, int type,
266 const char *name, const char *sep, TiffContext *s, AVFrame *frame)
269 case TIFF_DOUBLE: return ff_tadd_doubles_metadata(count, name, sep, &s->gb, s->le, &frame->metadata);
270 case TIFF_SHORT : return ff_tadd_shorts_metadata(count, name, sep, &s->gb, s->le, 0, &frame->metadata);
271 case TIFF_STRING: return ff_tadd_string_metadata(count, name, &s->gb, s->le, &frame->metadata);
272 default : return AVERROR_INVALIDDATA;
276 static void av_always_inline dng_blit(TiffContext *s, uint8_t *dst, int dst_stride,
277 const uint8_t *src, int src_stride, int width, int height, int is_u16);
279 static void av_always_inline horizontal_fill(TiffContext *s,
280 unsigned int bpp, uint8_t* dst,
281 int usePtr, const uint8_t *src,
282 uint8_t c, int width, int offset)
286 while (--width >= 0) {
287 dst[(width+offset)*8+7] = (usePtr ? src[width] : c) & 0x1;
288 dst[(width+offset)*8+6] = (usePtr ? src[width] : c) >> 1 & 0x1;
289 dst[(width+offset)*8+5] = (usePtr ? src[width] : c) >> 2 & 0x1;
290 dst[(width+offset)*8+4] = (usePtr ? src[width] : c) >> 3 & 0x1;
291 dst[(width+offset)*8+3] = (usePtr ? src[width] : c) >> 4 & 0x1;
292 dst[(width+offset)*8+2] = (usePtr ? src[width] : c) >> 5 & 0x1;
293 dst[(width+offset)*8+1] = (usePtr ? src[width] : c) >> 6 & 0x1;
294 dst[(width+offset)*8+0] = (usePtr ? src[width] : c) >> 7;
298 while (--width >= 0) {
299 dst[(width+offset)*4+3] = (usePtr ? src[width] : c) & 0x3;
300 dst[(width+offset)*4+2] = (usePtr ? src[width] : c) >> 2 & 0x3;
301 dst[(width+offset)*4+1] = (usePtr ? src[width] : c) >> 4 & 0x3;
302 dst[(width+offset)*4+0] = (usePtr ? src[width] : c) >> 6;
306 while (--width >= 0) {
307 dst[(width+offset)*2+1] = (usePtr ? src[width] : c) & 0xF;
308 dst[(width+offset)*2+0] = (usePtr ? src[width] : c) >> 4;
312 uint16_t *dst16 = (uint16_t *)dst;
314 init_get_bits8(&gb, src, width);
315 for (int i = 0; i < s->width; i++) {
316 dst16[i] = get_bits(&gb, 12) << 4;
322 memcpy(dst + offset, src, width);
324 memset(dst + offset, c, width);
329 static int deinvert_buffer(TiffContext *s, const uint8_t *src, int size)
333 av_fast_padded_malloc(&s->deinvert_buf, &s->deinvert_buf_size, size);
334 if (!s->deinvert_buf)
335 return AVERROR(ENOMEM);
336 for (i = 0; i < size; i++)
337 s->deinvert_buf[i] = ff_reverse[src[i]];
342 static void unpack_gray(TiffContext *s, AVFrame *p,
343 const uint8_t *src, int lnum, int width, int bpp)
346 uint16_t *dst = (uint16_t *)(p->data[0] + lnum * p->linesize[0]);
348 init_get_bits8(&gb, src, width);
350 for (int i = 0; i < s->width; i++) {
351 dst[i] = get_bits(&gb, bpp);
355 static void unpack_yuv(TiffContext *s, AVFrame *p,
356 const uint8_t *src, int lnum)
359 int w = (s->width - 1) / s->subsampling[0] + 1;
360 uint8_t *pu = &p->data[1][lnum / s->subsampling[1] * p->linesize[1]];
361 uint8_t *pv = &p->data[2][lnum / s->subsampling[1] * p->linesize[2]];
362 if (s->width % s->subsampling[0] || s->height % s->subsampling[1]) {
363 for (i = 0; i < w; i++) {
364 for (j = 0; j < s->subsampling[1]; j++)
365 for (k = 0; k < s->subsampling[0]; k++)
366 p->data[0][FFMIN(lnum + j, s->height-1) * p->linesize[0] +
367 FFMIN(i * s->subsampling[0] + k, s->width-1)] = *src++;
372 for (i = 0; i < w; i++) {
373 for (j = 0; j < s->subsampling[1]; j++)
374 for (k = 0; k < s->subsampling[0]; k++)
375 p->data[0][(lnum + j) * p->linesize[0] +
376 i * s->subsampling[0] + k] = *src++;
384 static int tiff_uncompress(uint8_t *dst, unsigned long *len, const uint8_t *src,
387 z_stream zstream = { 0 };
390 zstream.next_in = (uint8_t *)src;
391 zstream.avail_in = size;
392 zstream.next_out = dst;
393 zstream.avail_out = *len;
394 zret = inflateInit(&zstream);
396 av_log(NULL, AV_LOG_ERROR, "Inflate init error: %d\n", zret);
399 zret = inflate(&zstream, Z_SYNC_FLUSH);
400 inflateEnd(&zstream);
401 *len = zstream.total_out;
402 return zret == Z_STREAM_END ? Z_OK : zret;
405 static int tiff_unpack_zlib(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
406 const uint8_t *src, int size, int width, int lines,
407 int strip_start, int is_yuv)
410 unsigned long outlen;
412 outlen = width * lines;
413 zbuf = av_malloc(outlen);
415 return AVERROR(ENOMEM);
417 if ((ret = deinvert_buffer(s, src, size)) < 0) {
421 src = s->deinvert_buf;
423 ret = tiff_uncompress(zbuf, &outlen, src, size);
425 av_log(s->avctx, AV_LOG_ERROR,
426 "Uncompressing failed (%lu of %lu) with error %d\n", outlen,
427 (unsigned long)width * lines, ret);
429 return AVERROR_UNKNOWN;
432 for (line = 0; line < lines; line++) {
433 if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
434 horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
436 memcpy(dst, src, width);
439 unpack_yuv(s, p, dst, strip_start + line);
440 line += s->subsampling[1] - 1;
451 static int tiff_uncompress_lzma(uint8_t *dst, uint64_t *len, const uint8_t *src,
454 lzma_stream stream = LZMA_STREAM_INIT;
457 stream.next_in = (uint8_t *)src;
458 stream.avail_in = size;
459 stream.next_out = dst;
460 stream.avail_out = *len;
461 ret = lzma_stream_decoder(&stream, UINT64_MAX, 0);
462 if (ret != LZMA_OK) {
463 av_log(NULL, AV_LOG_ERROR, "LZMA init error: %d\n", ret);
466 ret = lzma_code(&stream, LZMA_RUN);
468 *len = stream.total_out;
469 return ret == LZMA_STREAM_END ? LZMA_OK : ret;
472 static int tiff_unpack_lzma(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
473 const uint8_t *src, int size, int width, int lines,
474 int strip_start, int is_yuv)
476 uint64_t outlen = width * (uint64_t)lines;
478 uint8_t *buf = av_malloc(outlen);
480 return AVERROR(ENOMEM);
482 if ((ret = deinvert_buffer(s, src, size)) < 0) {
486 src = s->deinvert_buf;
488 ret = tiff_uncompress_lzma(buf, &outlen, src, size);
489 if (ret != LZMA_OK) {
490 av_log(s->avctx, AV_LOG_ERROR,
491 "Uncompressing failed (%"PRIu64" of %"PRIu64") with error %d\n", outlen,
492 (uint64_t)width * lines, ret);
494 return AVERROR_UNKNOWN;
497 for (line = 0; line < lines; line++) {
498 if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
499 horizontal_fill(s, s->bpp, dst, 1, src, 0, width, 0);
501 memcpy(dst, src, width);
504 unpack_yuv(s, p, dst, strip_start + line);
505 line += s->subsampling[1] - 1;
515 static int tiff_unpack_fax(TiffContext *s, uint8_t *dst, int stride,
516 const uint8_t *src, int size, int width, int lines)
522 av_fast_padded_malloc(&s->fax_buffer, &s->fax_buffer_size, size);
523 src2 = s->fax_buffer;
526 av_log(s->avctx, AV_LOG_ERROR,
527 "Error allocating temporary buffer\n");
528 return AVERROR(ENOMEM);
531 if (!s->fill_order) {
532 memcpy(src2, src, size);
534 for (i = 0; i < size; i++)
535 src2[i] = ff_reverse[src[i]];
537 memset(src2 + size, 0, AV_INPUT_BUFFER_PADDING_SIZE);
538 ret = ff_ccitt_unpack(s->avctx, src2, size, dst, lines, stride,
539 s->compr, s->fax_opts);
540 if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
541 for (line = 0; line < lines; line++) {
542 horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
548 static int tiff_unpack_strip(TiffContext *s, AVFrame *p, uint8_t *dst, int stride,
549 const uint8_t *src, int size, int strip_start, int lines)
552 int c, line, pixels, code, ret;
553 const uint8_t *ssrc = src;
554 int width = ((s->width * s->bpp) + 7) >> 3;
555 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(p->format);
556 int is_yuv = !(desc->flags & AV_PIX_FMT_FLAG_RGB) &&
557 (desc->flags & AV_PIX_FMT_FLAG_PLANAR) &&
558 desc->nb_components >= 3;
562 width /= s->bppcount;
565 return AVERROR_INVALIDDATA;
568 int bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
569 s->subsampling[0] * s->subsampling[1] + 7) >> 3;
570 av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
571 if (s->yuv_line == NULL) {
572 av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
573 return AVERROR(ENOMEM);
578 width = (s->width - 1) / s->subsampling[0] + 1;
579 width = width * s->subsampling[0] * s->subsampling[1] + 2*width;
580 av_assert0(width <= bytes_per_row);
581 av_assert0(s->bpp == 24);
584 width = (s->bpp * s->width + 7) >> 3;
586 if (p->format == AV_PIX_FMT_GRAY12) {
587 av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, width);
588 if (s->yuv_line == NULL) {
589 av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
590 return AVERROR(ENOMEM);
596 if (s->compr == TIFF_DEFLATE || s->compr == TIFF_ADOBE_DEFLATE) {
598 return tiff_unpack_zlib(s, p, dst, stride, src, size, width, lines,
599 strip_start, is_yuv);
601 av_log(s->avctx, AV_LOG_ERROR,
602 "zlib support not enabled, "
603 "deflate compression not supported\n");
604 return AVERROR(ENOSYS);
607 if (s->compr == TIFF_LZMA) {
609 return tiff_unpack_lzma(s, p, dst, stride, src, size, width, lines,
610 strip_start, is_yuv);
612 av_log(s->avctx, AV_LOG_ERROR,
613 "LZMA support not enabled\n");
614 return AVERROR(ENOSYS);
617 if (s->compr == TIFF_LZW) {
619 if ((ret = deinvert_buffer(s, src, size)) < 0)
621 ssrc = src = s->deinvert_buf;
623 if (size > 1 && !src[0] && (src[1]&1)) {
624 av_log(s->avctx, AV_LOG_ERROR, "Old style LZW is unsupported\n");
626 if ((ret = ff_lzw_decode_init(s->lzw, 8, src, size, FF_LZW_TIFF)) < 0) {
627 av_log(s->avctx, AV_LOG_ERROR, "Error initializing LZW decoder\n");
630 for (line = 0; line < lines; line++) {
631 pixels = ff_lzw_decode(s->lzw, dst, width);
632 if (pixels < width) {
633 av_log(s->avctx, AV_LOG_ERROR, "Decoded only %i bytes of %i\n",
635 return AVERROR_INVALIDDATA;
637 if (s->bpp < 8 && s->avctx->pix_fmt == AV_PIX_FMT_PAL8)
638 horizontal_fill(s, s->bpp, dst, 1, dst, 0, width, 0);
640 unpack_yuv(s, p, dst, strip_start + line);
641 line += s->subsampling[1] - 1;
642 } else if (p->format == AV_PIX_FMT_GRAY12) {
643 unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
649 if (s->compr == TIFF_CCITT_RLE ||
650 s->compr == TIFF_G3 ||
651 s->compr == TIFF_G4) {
652 if (is_yuv || p->format == AV_PIX_FMT_GRAY12)
653 return AVERROR_INVALIDDATA;
655 return tiff_unpack_fax(s, dst, stride, src, size, width, lines);
658 bytestream2_init(&s->gb, src, size);
659 bytestream2_init_writer(&pb, dst, is_yuv ? s->yuv_line_size : (stride * lines));
661 is_dng = (s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG);
663 for (line = 0; line < lines; line++) {
664 if (src - ssrc > size) {
665 av_log(s->avctx, AV_LOG_ERROR, "Source data overread\n");
666 return AVERROR_INVALIDDATA;
669 if (bytestream2_get_bytes_left(&s->gb) == 0 || bytestream2_get_eof(&pb))
671 bytestream2_seek_p(&pb, stride * line, SEEK_SET);
674 if (ssrc + size - src < width)
675 return AVERROR_INVALIDDATA;
677 if (!s->fill_order) {
678 horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 || s->is_bayer),
679 dst, 1, src, 0, width, 0);
682 for (i = 0; i < width; i++)
683 dst[i] = ff_reverse[src[i]];
686 /* Color processing for DNG images with uncompressed strips (non-tiled) */
688 int is_u16, pixel_size_bytes, pixel_size_bits;
690 is_u16 = (s->bpp > 8);
691 pixel_size_bits = (is_u16 ? 16 : 8);
692 pixel_size_bytes = (is_u16 ? sizeof(uint16_t) : sizeof(uint8_t));
696 0, // no stride, only 1 line
698 0, // no stride, only 1 line
699 width / pixel_size_bytes * pixel_size_bits / s->bpp * s->bppcount, // need to account for [1, 16] bpp
707 for (pixels = 0; pixels < width;) {
708 if (ssrc + size - src < 2) {
709 av_log(s->avctx, AV_LOG_ERROR, "Read went out of bounds\n");
710 return AVERROR_INVALIDDATA;
712 code = s->fill_order ? (int8_t) ff_reverse[*src++]: (int8_t) *src++;
715 if (pixels + code > width ||
716 ssrc + size - src < code) {
717 av_log(s->avctx, AV_LOG_ERROR,
718 "Copy went out of bounds\n");
719 return AVERROR_INVALIDDATA;
721 horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
722 dst, 1, src, 0, code, pixels);
725 } else if (code != -128) { // -127..-1
727 if (pixels + code > width) {
728 av_log(s->avctx, AV_LOG_ERROR,
729 "Run went out of bounds\n");
730 return AVERROR_INVALIDDATA;
733 horizontal_fill(s, s->bpp * (s->avctx->pix_fmt == AV_PIX_FMT_PAL8),
734 dst, 0, NULL, c, code, pixels);
740 for (i = 0; i < width; i++)
741 dst[i] = ff_reverse[dst[i]];
746 unpack_yuv(s, p, dst, strip_start + line);
747 line += s->subsampling[1] - 1;
748 } else if (p->format == AV_PIX_FMT_GRAY12) {
749 unpack_gray(s, p, dst, strip_start + line, width, s->bpp);
756 static float av_always_inline linear_to_srgb(float value) {
757 if (value <= 0.0031308f)
758 return value * 12.92f;
760 return powf(value * 1.055f, 1.0f / 2.4f) - 0.055f;
764 * Map stored raw sensor values into linear reference values (see: DNG Specification - Chapter 5)
765 * Then convert to sRGB color space.
767 static uint16_t av_always_inline dng_process_color16(uint16_t value,
769 uint16_t black_level,
770 float scale_factor) {
773 // Lookup table lookup
777 // Black level subtraction
778 value = av_clip_uint16_c((unsigned)value - black_level);
781 value_norm = (float)value * scale_factor;
783 // Color space conversion (sRGB)
784 value = av_clip_uint16_c((uint16_t)(linear_to_srgb(value_norm) * 0xFFFF));
789 static uint16_t av_always_inline dng_process_color8(uint16_t value,
791 uint16_t black_level,
792 float scale_factor) {
793 return dng_process_color16(value, lut, black_level, scale_factor) >> 8;
796 static void dng_blit(TiffContext *s, uint8_t *dst, int dst_stride,
797 const uint8_t *src, int src_stride,
798 int width, int height, int is_u16)
803 scale_factor = 1.0f / (s->white_level - s->black_level);
806 for (line = 0; line < height; line++) {
807 uint16_t *dst_u16 = (uint16_t *)dst;
808 uint16_t *src_u16 = (uint16_t *)src;
810 for (col = 0; col < width; col++)
811 *dst_u16++ = dng_process_color16(*src_u16++, s->dng_lut, s->black_level, scale_factor);
813 dst += dst_stride * sizeof(uint16_t);
814 src += src_stride * sizeof(uint16_t);
817 for (line = 0; line < height; line++) {
818 for (col = 0; col < width; col++)
819 *dst++ = dng_process_color8(*src++, s->dng_lut, s->black_level, scale_factor);
827 static int dng_decode_jpeg_tile(AVCodecContext *avctx, AVFrame *frame,
828 int tile_byte_count, int x, int y, int w, int h)
830 TiffContext *s = avctx->priv_data;
832 uint8_t *dst_data, *src_data;
833 uint32_t dst_offset; /* offset from dst buffer in pixels */
834 int is_u16, pixel_size;
837 /* Prepare a packet and send to the MJPEG decoder */
838 av_init_packet(&jpkt);
839 jpkt.data = (uint8_t*)s->gb.buffer;
840 jpkt.size = tile_byte_count;
842 ret = avcodec_send_packet(s->avctx_mjpeg, &jpkt);
844 av_log(avctx, AV_LOG_ERROR, "Error submitting a packet for decoding\n");
848 ret = avcodec_receive_frame(s->avctx_mjpeg, s->jpgframe);
850 av_log(avctx, AV_LOG_ERROR, "JPEG decoding error: %s.\n", av_err2str(ret));
852 /* Normally skip, error if explode */
853 if (avctx->err_recognition & AV_EF_EXPLODE)
854 return AVERROR_INVALIDDATA;
859 /* Copy the outputted tile's pixels from 'jpgframe' to 'frame' (final buffer) */
861 is_u16 = (s->bpp > 8);
862 pixel_size = (is_u16 ? sizeof(uint16_t) : sizeof(uint8_t));
864 dst_offset = x + frame->linesize[0] * y / pixel_size;
865 dst_data = frame->data[0] + dst_offset * pixel_size;
866 src_data = s->jpgframe->data[0];
870 frame->linesize[0] / pixel_size,
872 s->jpgframe->linesize[0] / pixel_size,
877 av_frame_unref(s->jpgframe);
882 static int dng_decode_tiles(AVCodecContext *avctx, AVFrame *frame)
884 TiffContext *s = avctx->priv_data;
886 int tile_offset_offset, tile_offset;
887 int tile_byte_count_offset, tile_byte_count;
888 int tile_count_x, tile_count_y;
889 int tile_width, tile_length;
890 int has_width_leftover, has_height_leftover;
891 int tile_x = 0, tile_y = 0;
892 int pos_x = 0, pos_y = 0;
895 has_width_leftover = (s->width % s->tile_width != 0);
896 has_height_leftover = (s->height % s->tile_length != 0);
898 /* Calculate tile counts (round up) */
899 tile_count_x = (s->width + s->tile_width - 1) / s->tile_width;
900 tile_count_y = (s->height + s->tile_length - 1) / s->tile_length;
902 /* Iterate over the number of tiles */
903 for (tile_idx = 0; tile_idx < s->tile_count; tile_idx++) {
904 tile_x = tile_idx % tile_count_x;
905 tile_y = tile_idx / tile_count_x;
907 if (has_width_leftover && tile_x == tile_count_x - 1) // If on the right-most tile
908 tile_width = s->width % s->tile_width;
910 tile_width = s->tile_width;
912 if (has_height_leftover && tile_y == tile_count_y - 1) // If on the bottom-most tile
913 tile_length = s->height % s->tile_length;
915 tile_length = s->tile_length;
917 /* Read tile offset */
918 tile_offset_offset = s->tile_offsets_offset + tile_idx * sizeof(int);
919 bytestream2_seek(&s->gb, tile_offset_offset, SEEK_SET);
920 tile_offset = ff_tget_long(&s->gb, s->le);
922 /* Read tile byte size */
923 tile_byte_count_offset = s->tile_byte_counts_offset + tile_idx * sizeof(int);
924 bytestream2_seek(&s->gb, tile_byte_count_offset, SEEK_SET);
925 tile_byte_count = ff_tget_long(&s->gb, s->le);
927 /* Seek to tile data */
928 bytestream2_seek(&s->gb, tile_offset, SEEK_SET);
930 /* Decode JPEG tile and copy it in the reference frame */
931 ret = dng_decode_jpeg_tile(avctx, frame, tile_byte_count, pos_x, pos_y, tile_width, tile_length);
936 /* Advance current positions */
938 if (tile_x == tile_count_x - 1) { // If on the right edge
940 pos_y += tile_length;
947 static int dng_decode(AVCodecContext *avctx, AVFrame *frame, AVPacket *avpkt) {
950 TiffContext *s = avctx->priv_data;
952 s->jpgframe->width = s->tile_width;
953 s->jpgframe->height = s->tile_length;
955 s->avctx_mjpeg->width = s->tile_width;
956 s->avctx_mjpeg->height = s->tile_length;
958 /* Decode all tiles in a frame */
959 ret = dng_decode_tiles(avctx, frame);
963 /* Frame is ready to be output */
964 frame->pict_type = AV_PICTURE_TYPE_I;
965 frame->key_frame = 1;
970 static int init_image(TiffContext *s, ThreadFrame *frame)
973 int create_gray_palette = 0;
975 // make sure there is no aliasing in the following switch
976 if (s->bpp >= 100 || s->bppcount >= 10) {
977 av_log(s->avctx, AV_LOG_ERROR,
978 "Unsupported image parameters: bpp=%d, bppcount=%d\n",
979 s->bpp, s->bppcount);
980 return AVERROR_INVALIDDATA;
983 switch (s->planar * 1000 + s->bpp * 10 + s->bppcount + s->is_bayer * 10000) {
985 if (!s->palette_is_set) {
986 s->avctx->pix_fmt = AV_PIX_FMT_MONOBLACK;
991 s->avctx->pix_fmt = AV_PIX_FMT_PAL8;
992 if (!s->palette_is_set) {
993 create_gray_palette = 1;
997 s->avctx->pix_fmt = s->palette_is_set ? AV_PIX_FMT_PAL8 : AV_PIX_FMT_GRAY8;
1000 s->avctx->pix_fmt = AV_PIX_FMT_GRAY12;
1003 switch (AV_RL32(s->pattern)) {
1005 s->avctx->pix_fmt = AV_PIX_FMT_BAYER_RGGB8;
1008 s->avctx->pix_fmt = AV_PIX_FMT_BAYER_BGGR8;
1011 s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GBRG8;
1014 s->avctx->pix_fmt = AV_PIX_FMT_BAYER_GRBG8;
1017 av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1018 AV_RL32(s->pattern));
1019 return AVERROR_PATCHWELCOME;
1023 switch (AV_RL32(s->pattern)) {
1025 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_RGGB16LE : AV_PIX_FMT_BAYER_RGGB16BE;
1028 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_BGGR16LE : AV_PIX_FMT_BAYER_BGGR16BE;
1031 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GBRG16LE : AV_PIX_FMT_BAYER_GBRG16BE;
1034 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GRBG16LE : AV_PIX_FMT_BAYER_GRBG16BE;
1037 av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1038 AV_RL32(s->pattern));
1039 return AVERROR_PATCHWELCOME;
1043 switch (AV_RL32(s->pattern)) {
1045 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_RGGB16LE : AV_PIX_FMT_BAYER_RGGB16BE;
1048 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_BGGR16LE : AV_PIX_FMT_BAYER_BGGR16BE;
1051 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GBRG16LE : AV_PIX_FMT_BAYER_GBRG16BE;
1054 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_BAYER_GRBG16LE : AV_PIX_FMT_BAYER_GRBG16BE;
1057 av_log(s->avctx, AV_LOG_ERROR, "Unsupported Bayer pattern: 0x%X\n",
1058 AV_RL32(s->pattern));
1059 return AVERROR_PATCHWELCOME;
1063 if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1064 if (s->subsampling[0] == 1 && s->subsampling[1] == 1) {
1065 s->avctx->pix_fmt = AV_PIX_FMT_YUV444P;
1066 } else if (s->subsampling[0] == 2 && s->subsampling[1] == 1) {
1067 s->avctx->pix_fmt = AV_PIX_FMT_YUV422P;
1068 } else if (s->subsampling[0] == 4 && s->subsampling[1] == 1) {
1069 s->avctx->pix_fmt = AV_PIX_FMT_YUV411P;
1070 } else if (s->subsampling[0] == 1 && s->subsampling[1] == 2) {
1071 s->avctx->pix_fmt = AV_PIX_FMT_YUV440P;
1072 } else if (s->subsampling[0] == 2 && s->subsampling[1] == 2) {
1073 s->avctx->pix_fmt = AV_PIX_FMT_YUV420P;
1074 } else if (s->subsampling[0] == 4 && s->subsampling[1] == 4) {
1075 s->avctx->pix_fmt = AV_PIX_FMT_YUV410P;
1077 av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr subsampling\n");
1078 return AVERROR_PATCHWELCOME;
1081 s->avctx->pix_fmt = AV_PIX_FMT_RGB24;
1084 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GRAY16LE : AV_PIX_FMT_GRAY16BE;
1087 s->avctx->pix_fmt = AV_PIX_FMT_YA8;
1090 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_YA16LE : AV_PIX_FMT_YA16BE;
1093 s->avctx->pix_fmt = s->photometric == TIFF_PHOTOMETRIC_SEPARATED ? AV_PIX_FMT_RGB0 : AV_PIX_FMT_RGBA;
1096 if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED)
1097 s->avctx->pix_fmt = AV_PIX_FMT_RGBA;
1099 av_log(s->avctx, AV_LOG_ERROR,
1100 "bpp=40 without PHOTOMETRIC_SEPARATED is unsupported\n");
1101 return AVERROR_PATCHWELCOME;
1105 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGB48LE : AV_PIX_FMT_RGB48BE;
1108 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_RGBA64LE : AV_PIX_FMT_RGBA64BE;
1111 s->avctx->pix_fmt = AV_PIX_FMT_GBRP;
1114 s->avctx->pix_fmt = AV_PIX_FMT_GBRAP;
1117 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRP16LE : AV_PIX_FMT_GBRP16BE;
1120 s->avctx->pix_fmt = s->le ? AV_PIX_FMT_GBRAP16LE : AV_PIX_FMT_GBRAP16BE;
1123 av_log(s->avctx, AV_LOG_ERROR,
1124 "This format is not supported (bpp=%d, bppcount=%d)\n",
1125 s->bpp, s->bppcount);
1126 return AVERROR_INVALIDDATA;
1129 if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1130 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(s->avctx->pix_fmt);
1131 if((desc->flags & AV_PIX_FMT_FLAG_RGB) ||
1132 !(desc->flags & AV_PIX_FMT_FLAG_PLANAR) ||
1133 desc->nb_components < 3) {
1134 av_log(s->avctx, AV_LOG_ERROR, "Unsupported YCbCr variant\n");
1135 return AVERROR_INVALIDDATA;
1139 if (s->width != s->avctx->width || s->height != s->avctx->height) {
1140 ret = ff_set_dimensions(s->avctx, s->width, s->height);
1144 if ((ret = ff_thread_get_buffer(s->avctx, frame, 0)) < 0)
1146 if (s->avctx->pix_fmt == AV_PIX_FMT_PAL8) {
1147 if (!create_gray_palette)
1148 memcpy(frame->f->data[1], s->palette, sizeof(s->palette));
1150 /* make default grayscale pal */
1152 uint32_t *pal = (uint32_t *)frame->f->data[1];
1153 for (i = 0; i < 1<<s->bpp; i++)
1154 pal[i] = 0xFFU << 24 | i * 255 / ((1<<s->bpp) - 1) * 0x010101;
1160 static void set_sar(TiffContext *s, unsigned tag, unsigned num, unsigned den)
1162 int offset = tag == TIFF_YRES ? 2 : 0;
1163 s->res[offset++] = num;
1164 s->res[offset] = den;
1165 if (s->res[0] && s->res[1] && s->res[2] && s->res[3]) {
1166 uint64_t num = s->res[2] * (uint64_t)s->res[1];
1167 uint64_t den = s->res[0] * (uint64_t)s->res[3];
1168 if (num > INT64_MAX || den > INT64_MAX) {
1172 av_reduce(&s->avctx->sample_aspect_ratio.num, &s->avctx->sample_aspect_ratio.den,
1173 num, den, INT32_MAX);
1174 if (!s->avctx->sample_aspect_ratio.den)
1175 s->avctx->sample_aspect_ratio = (AVRational) {0, 1};
1179 static int tiff_decode_tag(TiffContext *s, AVFrame *frame)
1181 unsigned tag, type, count, off, value = 0, value2 = 1; // value2 is a denominator so init. to 1
1187 ret = ff_tread_tag(&s->gb, s->le, &tag, &type, &count, &start);
1192 off = bytestream2_tell(&s->gb);
1198 value = ff_tget(&s->gb, type, s->le);
1201 value = ff_tget(&s->gb, TIFF_LONG, s->le);
1202 value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
1215 s->is_thumbnail = (value != 0);
1225 av_log(s->avctx, AV_LOG_ERROR,
1226 "This format is not supported (bpp=%d, %d components)\n",
1228 return AVERROR_INVALIDDATA;
1230 s->bppcount = count;
1239 if (bytestream2_get_bytes_left(&s->gb) < type_sizes[type] * count)
1240 return AVERROR_INVALIDDATA;
1241 for (i = 0; i < count; i++)
1242 s->bpp += ff_tget(&s->gb, type, s->le);
1249 case TIFF_SAMPLES_PER_PIXEL:
1251 av_log(s->avctx, AV_LOG_ERROR,
1252 "Samples per pixel requires a single value, many provided\n");
1253 return AVERROR_INVALIDDATA;
1256 av_log(s->avctx, AV_LOG_ERROR,
1257 "Samples per pixel %d is too large\n", value);
1258 return AVERROR_INVALIDDATA;
1260 if (s->bppcount == 1)
1262 s->bppcount = value;
1266 av_log(s->avctx, AV_LOG_DEBUG, "compression: %d\n", s->compr);
1272 case TIFF_CCITT_RLE:
1279 case TIFF_ADOBE_DEFLATE:
1283 av_log(s->avctx, AV_LOG_ERROR, "Deflate: ZLib not compiled in\n");
1284 return AVERROR(ENOSYS);
1294 av_log(s->avctx, AV_LOG_ERROR, "LZMA not compiled in\n");
1295 return AVERROR(ENOSYS);
1298 av_log(s->avctx, AV_LOG_ERROR, "Unknown compression method %i\n",
1300 return AVERROR_INVALIDDATA;
1303 case TIFF_ROWSPERSTRIP:
1304 if (!value || (type == TIFF_LONG && value == UINT_MAX))
1306 s->rps = FFMIN(value, s->height);
1308 case TIFF_STRIP_OFFS:
1310 if (value > INT_MAX) {
1311 av_log(s->avctx, AV_LOG_ERROR,
1312 "strippos %u too large\n", value);
1313 return AVERROR_INVALIDDATA;
1316 s->stripoff = value;
1324 case TIFF_STRIP_SIZE:
1326 if (value > INT_MAX) {
1327 av_log(s->avctx, AV_LOG_ERROR,
1328 "stripsize %u too large\n", value);
1329 return AVERROR_INVALIDDATA;
1331 s->stripsizesoff = 0;
1332 s->stripsize = value;
1335 s->stripsizesoff = off;
1342 set_sar(s, tag, value, value2);
1344 case TIFF_TILE_OFFSETS:
1345 s->tile_offsets_offset = off;
1346 s->tile_count = count;
1349 case TIFF_TILE_BYTE_COUNTS:
1350 s->tile_byte_counts_offset = off;
1352 case TIFF_TILE_LENGTH:
1353 s->tile_length = value;
1355 case TIFF_TILE_WIDTH:
1356 s->tile_width = value;
1358 case TIFF_PREDICTOR:
1359 s->predictor = value;
1365 s->sub_ifd = ff_tget(&s->gb, TIFF_LONG, s->le); /** Only get the first SubIFD */
1367 case DNG_LINEARIZATION_TABLE: {
1368 uint32_t lut_offset = value;
1369 uint32_t lut_size = count;
1370 uint32_t lut_wanted_size = 1 << s->bpp;
1371 if (lut_wanted_size != lut_size)
1372 av_log(s->avctx, AV_LOG_WARNING, "DNG contains LUT with invalid size (%"PRIu32"), disabling LUT\n", lut_size);
1373 else if (lut_offset >= bytestream2_size(&s->gb))
1374 av_log(s->avctx, AV_LOG_WARNING, "DNG contains LUT with invalid offset (%"PRIu32"), disabling LUT\n", lut_offset);
1376 s->dng_lut = (uint16_t*)(s->gb.buffer + lut_offset);
1379 case DNG_BLACK_LEVEL:
1380 if (count > 1) { /* Use the first value in the pattern (assume they're all the same) */
1381 if (type == TIFF_RATIONAL) {
1382 value = ff_tget(&s->gb, TIFF_LONG, s->le);
1383 value2 = ff_tget(&s->gb, TIFF_LONG, s->le);
1385 s->black_level = value / value2;
1387 s->black_level = ff_tget(&s->gb, type, s->le);
1388 av_log(s->avctx, AV_LOG_WARNING, "Assuming black level pattern values are identical\n");
1390 s->black_level = value / value2;
1393 case DNG_WHITE_LEVEL:
1394 s->white_level = value;
1396 case TIFF_CFA_PATTERN_DIM:
1397 if (count != 2 || (ff_tget(&s->gb, type, s->le) != 2 &&
1398 ff_tget(&s->gb, type, s->le) != 2)) {
1399 av_log(s->avctx, AV_LOG_ERROR, "CFA Pattern dimensions are not 2x2\n");
1400 return AVERROR_INVALIDDATA;
1403 case TIFF_CFA_PATTERN:
1405 s->pattern[0] = ff_tget(&s->gb, type, s->le);
1406 s->pattern[1] = ff_tget(&s->gb, type, s->le);
1407 s->pattern[2] = ff_tget(&s->gb, type, s->le);
1408 s->pattern[3] = ff_tget(&s->gb, type, s->le);
1410 case TIFF_PHOTOMETRIC:
1412 case TIFF_PHOTOMETRIC_WHITE_IS_ZERO:
1413 case TIFF_PHOTOMETRIC_BLACK_IS_ZERO:
1414 case TIFF_PHOTOMETRIC_RGB:
1415 case TIFF_PHOTOMETRIC_PALETTE:
1416 case TIFF_PHOTOMETRIC_SEPARATED:
1417 case TIFF_PHOTOMETRIC_YCBCR:
1418 case TIFF_PHOTOMETRIC_CFA:
1419 s->photometric = value;
1421 case TIFF_PHOTOMETRIC_ALPHA_MASK:
1422 case TIFF_PHOTOMETRIC_CIE_LAB:
1423 case TIFF_PHOTOMETRIC_ICC_LAB:
1424 case TIFF_PHOTOMETRIC_ITU_LAB:
1425 case TIFF_PHOTOMETRIC_LOG_L:
1426 case TIFF_PHOTOMETRIC_LOG_LUV:
1427 case TIFF_PHOTOMETRIC_LINEAR_RAW:
1428 avpriv_report_missing_feature(s->avctx,
1429 "PhotometricInterpretation 0x%04X",
1431 return AVERROR_PATCHWELCOME;
1433 av_log(s->avctx, AV_LOG_ERROR, "PhotometricInterpretation %u is "
1434 "unknown\n", value);
1435 return AVERROR_INVALIDDATA;
1438 case TIFF_FILL_ORDER:
1439 if (value < 1 || value > 2) {
1440 av_log(s->avctx, AV_LOG_ERROR,
1441 "Unknown FillOrder value %d, trying default one\n", value);
1444 s->fill_order = value - 1;
1447 GetByteContext pal_gb[3];
1448 off = type_sizes[type];
1449 if (count / 3 > 256 ||
1450 bytestream2_get_bytes_left(&s->gb) < count / 3 * off * 3)
1451 return AVERROR_INVALIDDATA;
1453 pal_gb[0] = pal_gb[1] = pal_gb[2] = s->gb;
1454 bytestream2_skip(&pal_gb[1], count / 3 * off);
1455 bytestream2_skip(&pal_gb[2], count / 3 * off * 2);
1457 off = (type_sizes[type] - 1) << 3;
1459 av_log(s->avctx, AV_LOG_ERROR, "palette shift %d is out of range\n", off);
1460 return AVERROR_INVALIDDATA;
1463 for (i = 0; i < count / 3; i++) {
1464 uint32_t p = 0xFF000000;
1465 p |= (ff_tget(&pal_gb[0], type, s->le) >> off) << 16;
1466 p |= (ff_tget(&pal_gb[1], type, s->le) >> off) << 8;
1467 p |= ff_tget(&pal_gb[2], type, s->le) >> off;
1470 s->palette_is_set = 1;
1474 s->planar = value == 2;
1476 case TIFF_YCBCR_SUBSAMPLING:
1478 av_log(s->avctx, AV_LOG_ERROR, "subsample count invalid\n");
1479 return AVERROR_INVALIDDATA;
1481 for (i = 0; i < count; i++) {
1482 s->subsampling[i] = ff_tget(&s->gb, type, s->le);
1483 if (s->subsampling[i] <= 0) {
1484 av_log(s->avctx, AV_LOG_ERROR, "subsampling %d is invalid\n", s->subsampling[i]);
1485 s->subsampling[i] = 1;
1486 return AVERROR_INVALIDDATA;
1490 case TIFF_T4OPTIONS:
1491 if (s->compr == TIFF_G3)
1492 s->fax_opts = value;
1494 case TIFF_T6OPTIONS:
1495 if (s->compr == TIFF_G4)
1496 s->fax_opts = value;
1498 #define ADD_METADATA(count, name, sep)\
1499 if ((ret = add_metadata(count, type, name, sep, s, frame)) < 0) {\
1500 av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");\
1503 case TIFF_MODEL_PIXEL_SCALE:
1504 ADD_METADATA(count, "ModelPixelScaleTag", NULL);
1506 case TIFF_MODEL_TRANSFORMATION:
1507 ADD_METADATA(count, "ModelTransformationTag", NULL);
1509 case TIFF_MODEL_TIEPOINT:
1510 ADD_METADATA(count, "ModelTiepointTag", NULL);
1512 case TIFF_GEO_KEY_DIRECTORY:
1513 if (s->geotag_count) {
1514 avpriv_request_sample(s->avctx, "Multiple geo key directories\n");
1515 return AVERROR_INVALIDDATA;
1517 ADD_METADATA(1, "GeoTIFF_Version", NULL);
1518 ADD_METADATA(2, "GeoTIFF_Key_Revision", ".");
1519 s->geotag_count = ff_tget_short(&s->gb, s->le);
1520 if (s->geotag_count > count / 4 - 1) {
1521 s->geotag_count = count / 4 - 1;
1522 av_log(s->avctx, AV_LOG_WARNING, "GeoTIFF key directory buffer shorter than specified\n");
1524 if ( bytestream2_get_bytes_left(&s->gb) < s->geotag_count * sizeof(int16_t) * 4
1525 || s->geotag_count == 0) {
1526 s->geotag_count = 0;
1529 s->geotags = av_mallocz_array(s->geotag_count, sizeof(TiffGeoTag));
1531 av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1532 s->geotag_count = 0;
1535 for (i = 0; i < s->geotag_count; i++) {
1536 s->geotags[i].key = ff_tget_short(&s->gb, s->le);
1537 s->geotags[i].type = ff_tget_short(&s->gb, s->le);
1538 s->geotags[i].count = ff_tget_short(&s->gb, s->le);
1540 if (!s->geotags[i].type)
1541 s->geotags[i].val = get_geokey_val(s->geotags[i].key, ff_tget_short(&s->gb, s->le));
1543 s->geotags[i].offset = ff_tget_short(&s->gb, s->le);
1546 case TIFF_GEO_DOUBLE_PARAMS:
1547 if (count >= INT_MAX / sizeof(int64_t))
1548 return AVERROR_INVALIDDATA;
1549 if (bytestream2_get_bytes_left(&s->gb) < count * sizeof(int64_t))
1550 return AVERROR_INVALIDDATA;
1551 dp = av_malloc_array(count, sizeof(double));
1553 av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1556 for (i = 0; i < count; i++)
1557 dp[i] = ff_tget_double(&s->gb, s->le);
1558 for (i = 0; i < s->geotag_count; i++) {
1559 if (s->geotags[i].type == TIFF_GEO_DOUBLE_PARAMS) {
1560 if (s->geotags[i].count == 0
1561 || s->geotags[i].offset + s->geotags[i].count > count) {
1562 av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1563 } else if (s->geotags[i].val) {
1564 av_log(s->avctx, AV_LOG_WARNING, "Duplicate GeoTIFF key %d\n", s->geotags[i].key);
1566 char *ap = doubles2str(&dp[s->geotags[i].offset], s->geotags[i].count, ", ");
1568 av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1570 return AVERROR(ENOMEM);
1572 s->geotags[i].val = ap;
1578 case TIFF_GEO_ASCII_PARAMS:
1579 pos = bytestream2_tell(&s->gb);
1580 for (i = 0; i < s->geotag_count; i++) {
1581 if (s->geotags[i].type == TIFF_GEO_ASCII_PARAMS) {
1582 if (s->geotags[i].count == 0
1583 || s->geotags[i].offset + s->geotags[i].count > count) {
1584 av_log(s->avctx, AV_LOG_WARNING, "Invalid GeoTIFF key %d\n", s->geotags[i].key);
1588 bytestream2_seek(&s->gb, pos + s->geotags[i].offset, SEEK_SET);
1589 if (bytestream2_get_bytes_left(&s->gb) < s->geotags[i].count)
1590 return AVERROR_INVALIDDATA;
1591 if (s->geotags[i].val)
1592 return AVERROR_INVALIDDATA;
1593 ap = av_malloc(s->geotags[i].count);
1595 av_log(s->avctx, AV_LOG_ERROR, "Error allocating temporary buffer\n");
1596 return AVERROR(ENOMEM);
1598 bytestream2_get_bufferu(&s->gb, ap, s->geotags[i].count);
1599 ap[s->geotags[i].count - 1] = '\0'; //replace the "|" delimiter with a 0 byte
1600 s->geotags[i].val = ap;
1606 ADD_METADATA(count, "artist", NULL);
1608 case TIFF_COPYRIGHT:
1609 ADD_METADATA(count, "copyright", NULL);
1612 ADD_METADATA(count, "date", NULL);
1614 case TIFF_DOCUMENT_NAME:
1615 ADD_METADATA(count, "document_name", NULL);
1617 case TIFF_HOST_COMPUTER:
1618 ADD_METADATA(count, "computer", NULL);
1620 case TIFF_IMAGE_DESCRIPTION:
1621 ADD_METADATA(count, "description", NULL);
1624 ADD_METADATA(count, "make", NULL);
1627 ADD_METADATA(count, "model", NULL);
1629 case TIFF_PAGE_NAME:
1630 ADD_METADATA(count, "page_name", NULL);
1632 case TIFF_PAGE_NUMBER:
1633 ADD_METADATA(count, "page_number", " / ");
1634 // need to seek back to re-read the page number
1635 bytestream2_seek(&s->gb, -count * sizeof(uint16_t), SEEK_CUR);
1636 // read the page number
1637 s->cur_page = ff_tget(&s->gb, TIFF_SHORT, s->le);
1638 // get back to where we were before the previous seek
1639 bytestream2_seek(&s->gb, count * sizeof(uint16_t) - sizeof(uint16_t), SEEK_CUR);
1641 case TIFF_SOFTWARE_NAME:
1642 ADD_METADATA(count, "software", NULL);
1646 unsigned int ver[4];
1647 ver[0] = ff_tget(&s->gb, type, s->le);
1648 ver[1] = ff_tget(&s->gb, type, s->le);
1649 ver[2] = ff_tget(&s->gb, type, s->le);
1650 ver[3] = ff_tget(&s->gb, type, s->le);
1652 av_log(s->avctx, AV_LOG_DEBUG, "DNG file, version %u.%u.%u.%u\n",
1653 ver[0], ver[1], ver[2], ver[3]);
1655 tiff_set_type(s, TIFF_TYPE_DNG);
1658 case CINEMADNG_TIME_CODES:
1659 case CINEMADNG_FRAME_RATE:
1660 case CINEMADNG_T_STOP:
1661 case CINEMADNG_REEL_NAME:
1662 case CINEMADNG_CAMERA_LABEL:
1663 tiff_set_type(s, TIFF_TYPE_CINEMADNG);
1666 if (s->avctx->err_recognition & AV_EF_EXPLODE) {
1667 av_log(s->avctx, AV_LOG_ERROR,
1668 "Unknown or unsupported tag %d/0x%0X\n",
1670 return AVERROR_INVALIDDATA;
1675 av_log(s->avctx, AV_LOG_ERROR,
1676 "This format is not supported (bpp=%d, %d components)\n",
1679 return AVERROR_INVALIDDATA;
1681 bytestream2_seek(&s->gb, start, SEEK_SET);
1685 static int decode_frame(AVCodecContext *avctx,
1686 void *data, int *got_frame, AVPacket *avpkt)
1688 TiffContext *const s = avctx->priv_data;
1689 AVFrame *const p = data;
1690 ThreadFrame frame = { .f = data };
1691 unsigned off, last_off;
1692 int le, ret, plane, planes;
1693 int i, j, entries, stride;
1694 unsigned soff, ssize;
1696 GetByteContext stripsizes;
1697 GetByteContext stripdata;
1698 int retry_for_subifd, retry_for_page;
1700 bytestream2_init(&s->gb, avpkt->data, avpkt->size);
1702 // parse image header
1703 if ((ret = ff_tdecode_header(&s->gb, &le, &off))) {
1704 av_log(avctx, AV_LOG_ERROR, "Invalid TIFF header\n");
1706 } else if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1707 av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1708 return AVERROR_INVALIDDATA;
1711 // TIFF_BPP is not a required tag and defaults to 1
1713 s->tiff_type = TIFF_TYPE_TIFF;
1715 s->is_thumbnail = 0;
1716 s->bppcount = s->bpp = 1;
1717 s->photometric = TIFF_PHOTOMETRIC_NONE;
1718 s->compr = TIFF_RAW;
1728 // Reset these offsets so we can tell if they were set this frame
1729 s->stripsizesoff = s->strippos = 0;
1730 /* parse image file directory */
1731 bytestream2_seek(&s->gb, off, SEEK_SET);
1732 entries = ff_tget_short(&s->gb, le);
1733 if (bytestream2_get_bytes_left(&s->gb) < entries * 12)
1734 return AVERROR_INVALIDDATA;
1735 for (i = 0; i < entries; i++) {
1736 if ((ret = tiff_decode_tag(s, p)) < 0)
1740 if (s->get_thumbnail && !s->is_thumbnail) {
1741 av_log(avctx, AV_LOG_INFO, "No embedded thumbnail present\n");
1745 /** whether we should process this IFD's SubIFD */
1746 retry_for_subifd = s->sub_ifd && (s->get_subimage || (!s->get_thumbnail && s->is_thumbnail));
1747 /** whether we should process this multi-page IFD's next page */
1748 retry_for_page = s->get_page && s->cur_page + 1 < s->get_page; // get_page is 1-indexed
1751 if (retry_for_page) {
1752 // set offset to the next IFD
1753 off = ff_tget_long(&s->gb, le);
1754 } else if (retry_for_subifd) {
1755 // set offset to the SubIFD
1759 if (retry_for_subifd || retry_for_page) {
1761 av_log(avctx, AV_LOG_ERROR, "Requested entry not found\n");
1762 return AVERROR_INVALIDDATA;
1764 if (off <= last_off) {
1765 avpriv_request_sample(s->avctx, "non increasing IFD offset\n");
1766 return AVERROR_INVALIDDATA;
1768 if (off >= UINT_MAX - 14 || avpkt->size < off + 14) {
1769 av_log(avctx, AV_LOG_ERROR, "IFD offset is greater than image size\n");
1770 return AVERROR_INVALIDDATA;
1776 for (i = 0; i<s->geotag_count; i++) {
1777 const char *keyname = get_geokey_name(s->geotags[i].key);
1779 av_log(avctx, AV_LOG_WARNING, "Unknown or unsupported GeoTIFF key %d\n", s->geotags[i].key);
1782 if (get_geokey_type(s->geotags[i].key) != s->geotags[i].type) {
1783 av_log(avctx, AV_LOG_WARNING, "Type of GeoTIFF key %d is wrong\n", s->geotags[i].key);
1786 ret = av_dict_set(&p->metadata, keyname, s->geotags[i].val, 0);
1788 av_log(avctx, AV_LOG_ERROR, "Writing metadata with key '%s' failed\n", keyname);
1793 if (!s->is_tiled && !s->strippos && !s->stripoff) {
1794 av_log(avctx, AV_LOG_ERROR, "Image data is missing\n");
1795 return AVERROR_INVALIDDATA;
1797 /* now we have the data and may start decoding */
1798 if ((ret = init_image(s, &frame)) < 0)
1802 if (s->strips == 1 && !s->stripsize) {
1803 av_log(avctx, AV_LOG_WARNING, "Image data size missing\n");
1804 s->stripsize = avpkt->size - s->stripoff;
1807 if (s->stripsizesoff) {
1808 if (s->stripsizesoff >= (unsigned)avpkt->size)
1809 return AVERROR_INVALIDDATA;
1810 bytestream2_init(&stripsizes, avpkt->data + s->stripsizesoff,
1811 avpkt->size - s->stripsizesoff);
1814 if (s->strippos >= (unsigned)avpkt->size)
1815 return AVERROR_INVALIDDATA;
1816 bytestream2_init(&stripdata, avpkt->data + s->strippos,
1817 avpkt->size - s->strippos);
1820 if (s->rps <= 0 || s->rps % s->subsampling[1]) {
1821 av_log(avctx, AV_LOG_ERROR, "rps %d invalid\n", s->rps);
1822 return AVERROR_INVALIDDATA;
1826 /* Handle DNG images with JPEG-compressed tiles */
1828 if ((s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG) && s->is_tiled) {
1830 avpriv_report_missing_feature(avctx, "DNG uncompressed tiled images");
1831 return AVERROR_PATCHWELCOME;
1832 } else if (!s->is_bayer) {
1833 avpriv_report_missing_feature(avctx, "DNG JPG-compressed tiled non-bayer-encoded images");
1834 return AVERROR_PATCHWELCOME;
1836 if ((ret = dng_decode(avctx, (AVFrame*)data, avpkt)) > 0)
1842 /* Handle TIFF images and DNG images with uncompressed strips (non-tiled) */
1844 planes = s->planar ? s->bppcount : 1;
1845 for (plane = 0; plane < planes; plane++) {
1846 uint8_t *five_planes = NULL;
1847 int remaining = avpkt->size;
1849 stride = p->linesize[plane];
1850 dst = p->data[plane];
1851 if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1852 s->avctx->pix_fmt == AV_PIX_FMT_RGBA) {
1853 stride = stride * 5 / 4;
1855 dst = av_malloc(stride * s->height);
1857 return AVERROR(ENOMEM);
1859 for (i = 0; i < s->height; i += s->rps) {
1861 dst += s->rps * stride;
1862 if (s->stripsizesoff)
1863 ssize = ff_tget(&stripsizes, s->sstype, le);
1865 ssize = s->stripsize;
1868 soff = ff_tget(&stripdata, s->sot, le);
1872 if (soff > avpkt->size || ssize > avpkt->size - soff || ssize > remaining) {
1873 av_log(avctx, AV_LOG_ERROR, "Invalid strip size/offset\n");
1874 av_freep(&five_planes);
1875 return AVERROR_INVALIDDATA;
1878 if ((ret = tiff_unpack_strip(s, p, dst, stride, avpkt->data + soff, ssize, i,
1879 FFMIN(s->rps, s->height - i))) < 0) {
1880 if (avctx->err_recognition & AV_EF_EXPLODE) {
1881 av_freep(&five_planes);
1887 decoded_height = FFMIN(i, s->height);
1889 if (s->predictor == 2) {
1890 if (s->photometric == TIFF_PHOTOMETRIC_YCBCR) {
1891 av_log(s->avctx, AV_LOG_ERROR, "predictor == 2 with YUV is unsupported");
1892 return AVERROR_PATCHWELCOME;
1894 dst = five_planes ? five_planes : p->data[plane];
1897 soff = FFMAX(soff / s->bppcount, 1);
1898 ssize = s->width * soff;
1899 if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48LE ||
1900 s->avctx->pix_fmt == AV_PIX_FMT_RGBA64LE ||
1901 s->avctx->pix_fmt == AV_PIX_FMT_GRAY16LE ||
1902 s->avctx->pix_fmt == AV_PIX_FMT_YA16LE ||
1903 s->avctx->pix_fmt == AV_PIX_FMT_GBRP16LE ||
1904 s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16LE) {
1905 for (i = 0; i < decoded_height; i++) {
1906 for (j = soff; j < ssize; j += 2)
1907 AV_WL16(dst + j, AV_RL16(dst + j) + AV_RL16(dst + j - soff));
1910 } else if (s->avctx->pix_fmt == AV_PIX_FMT_RGB48BE ||
1911 s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE ||
1912 s->avctx->pix_fmt == AV_PIX_FMT_GRAY16BE ||
1913 s->avctx->pix_fmt == AV_PIX_FMT_YA16BE ||
1914 s->avctx->pix_fmt == AV_PIX_FMT_GBRP16BE ||
1915 s->avctx->pix_fmt == AV_PIX_FMT_GBRAP16BE) {
1916 for (i = 0; i < decoded_height; i++) {
1917 for (j = soff; j < ssize; j += 2)
1918 AV_WB16(dst + j, AV_RB16(dst + j) + AV_RB16(dst + j - soff));
1922 for (i = 0; i < decoded_height; i++) {
1923 for (j = soff; j < ssize; j++)
1924 dst[j] += dst[j - soff];
1930 if (s->photometric == TIFF_PHOTOMETRIC_WHITE_IS_ZERO) {
1931 int c = (s->avctx->pix_fmt == AV_PIX_FMT_PAL8 ? (1<<s->bpp) - 1 : 255);
1932 dst = p->data[plane];
1933 for (i = 0; i < s->height; i++) {
1934 for (j = 0; j < stride; j++)
1935 dst[j] = c - dst[j];
1940 if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1941 (s->avctx->pix_fmt == AV_PIX_FMT_RGB0 || s->avctx->pix_fmt == AV_PIX_FMT_RGBA)) {
1942 int x = s->avctx->pix_fmt == AV_PIX_FMT_RGB0 ? 4 : 5;
1943 uint8_t *src = five_planes ? five_planes : p->data[plane];
1944 dst = p->data[plane];
1945 for (i = 0; i < s->height; i++) {
1946 for (j = 0; j < s->width; j++) {
1947 int k = 255 - src[x * j + 3];
1948 int r = (255 - src[x * j ]) * k;
1949 int g = (255 - src[x * j + 1]) * k;
1950 int b = (255 - src[x * j + 2]) * k;
1951 dst[4 * j ] = r * 257 >> 16;
1952 dst[4 * j + 1] = g * 257 >> 16;
1953 dst[4 * j + 2] = b * 257 >> 16;
1954 dst[4 * j + 3] = s->avctx->pix_fmt == AV_PIX_FMT_RGBA ? src[x * j + 4] : 255;
1957 dst += p->linesize[plane];
1959 av_freep(&five_planes);
1960 } else if (s->photometric == TIFF_PHOTOMETRIC_SEPARATED &&
1961 s->avctx->pix_fmt == AV_PIX_FMT_RGBA64BE) {
1962 dst = p->data[plane];
1963 for (i = 0; i < s->height; i++) {
1964 for (j = 0; j < s->width; j++) {
1965 uint64_t k = 65535 - AV_RB16(dst + 8 * j + 6);
1966 uint64_t r = (65535 - AV_RB16(dst + 8 * j )) * k;
1967 uint64_t g = (65535 - AV_RB16(dst + 8 * j + 2)) * k;
1968 uint64_t b = (65535 - AV_RB16(dst + 8 * j + 4)) * k;
1969 AV_WB16(dst + 8 * j , r * 65537 >> 32);
1970 AV_WB16(dst + 8 * j + 2, g * 65537 >> 32);
1971 AV_WB16(dst + 8 * j + 4, b * 65537 >> 32);
1972 AV_WB16(dst + 8 * j + 6, 65535);
1974 dst += p->linesize[plane];
1979 if (s->planar && s->bppcount > 2) {
1980 FFSWAP(uint8_t*, p->data[0], p->data[2]);
1981 FFSWAP(int, p->linesize[0], p->linesize[2]);
1982 FFSWAP(uint8_t*, p->data[0], p->data[1]);
1983 FFSWAP(int, p->linesize[0], p->linesize[1]);
1986 if (s->is_bayer && s->white_level && s->bpp == 16 &&
1987 !(s->tiff_type == TIFF_TYPE_DNG || s->tiff_type == TIFF_TYPE_CINEMADNG)) {
1988 uint16_t *dst = (uint16_t *)p->data[0];
1989 for (i = 0; i < s->height; i++) {
1990 for (j = 0; j < s->width; j++)
1991 dst[j] = FFMIN((dst[j] / (float)s->white_level) * 65535, 65535);
2001 static av_cold int tiff_init(AVCodecContext *avctx)
2003 TiffContext *s = avctx->priv_data;
2004 const AVCodec *codec;
2010 s->subsampling[1] = 1;
2012 ff_lzw_decode_open(&s->lzw);
2014 return AVERROR(ENOMEM);
2015 ff_ccitt_unpack_init();
2017 /* Allocate JPEG frame */
2018 s->jpgframe = av_frame_alloc();
2020 return AVERROR(ENOMEM);
2022 /* Prepare everything needed for JPEG decoding */
2023 codec = avcodec_find_decoder(AV_CODEC_ID_MJPEG);
2026 s->avctx_mjpeg = avcodec_alloc_context3(codec);
2027 if (!s->avctx_mjpeg)
2028 return AVERROR(ENOMEM);
2029 s->avctx_mjpeg->flags = avctx->flags;
2030 s->avctx_mjpeg->flags2 = avctx->flags2;
2031 s->avctx_mjpeg->dct_algo = avctx->dct_algo;
2032 s->avctx_mjpeg->idct_algo = avctx->idct_algo;
2033 ret = ff_codec_open2_recursive(s->avctx_mjpeg, codec, NULL);
2035 av_frame_free(&s->jpgframe);
2036 avcodec_free_context(&s->avctx_mjpeg);
2043 static av_cold int tiff_end(AVCodecContext *avctx)
2045 TiffContext *const s = avctx->priv_data;
2049 ff_lzw_decode_close(&s->lzw);
2050 av_freep(&s->deinvert_buf);
2051 s->deinvert_buf_size = 0;
2052 av_freep(&s->yuv_line);
2053 s->yuv_line_size = 0;
2054 av_freep(&s->fax_buffer);
2055 s->fax_buffer_size = 0;
2056 av_frame_free(&s->jpgframe);
2057 avcodec_free_context(&s->avctx_mjpeg);
2061 #define OFFSET(x) offsetof(TiffContext, x)
2062 static const AVOption tiff_options[] = {
2063 { "subimage", "decode subimage instead if available", OFFSET(get_subimage), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2064 { "thumbnail", "decode embedded thumbnail subimage instead if available", OFFSET(get_thumbnail), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2065 { "page", "page number of multi-page image to decode (starting from 1)", OFFSET(get_page), AV_OPT_TYPE_INT, {.i64=0}, 0, UINT16_MAX, AV_OPT_FLAG_DECODING_PARAM | AV_OPT_FLAG_VIDEO_PARAM },
2069 static const AVClass tiff_decoder_class = {
2070 .class_name = "TIFF decoder",
2071 .item_name = av_default_item_name,
2072 .option = tiff_options,
2073 .version = LIBAVUTIL_VERSION_INT,
2076 AVCodec ff_tiff_decoder = {
2078 .long_name = NULL_IF_CONFIG_SMALL("TIFF image"),
2079 .type = AVMEDIA_TYPE_VIDEO,
2080 .id = AV_CODEC_ID_TIFF,
2081 .priv_data_size = sizeof(TiffContext),
2084 .decode = decode_frame,
2085 .init_thread_copy = ONLY_IF_THREADS_ENABLED(tiff_init),
2086 .capabilities = AV_CODEC_CAP_DR1 | AV_CODEC_CAP_FRAME_THREADS,
2087 .priv_class = &tiff_decoder_class,