[deb_ffmpeg.git] / ffmpeg / libavcodec / tiffenc.c

/*
 * TIFF image encoder
 * Copyright (c) 2007 Bartlomiej Wolowiec
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/**
 * @file
 * TIFF image encoder
 * @author Bartlomiej Wolowiec
 */

#include "config.h"
#if CONFIG_ZLIB
#include <zlib.h>
#endif

#include "libavutil/imgutils.h"
#include "libavutil/log.h"
#include "libavutil/opt.h"
#include "libavutil/pixdesc.h"
#include "avcodec.h"
#include "bytestream.h"
#include "internal.h"
#include "lzw.h"
#include "put_bits.h"
#include "rle.h"
#include "tiff.h"

#define TIFF_MAX_ENTRY 32

/** sizes of various TIFF field types (string size = 1)*/
static const uint8_t type_sizes2[14] = {
    0, 1, 1, 2, 4, 8, 1, 1, 2, 4, 8, 4, 8, 4
};

typedef struct TiffEncoderContext {
    AVClass *class;                         ///< for private options
    AVCodecContext *avctx;

    int width;                              ///< picture width
    int height;                             ///< picture height
    unsigned int bpp;                       ///< bits per pixel
    int compr;                              ///< compression level
    int bpp_tab_size;                       ///< bpp_tab size
    enum TiffPhotometric photometric_interpretation;  ///< photometric interpretation
    int strips;                             ///< number of strips
    uint32_t *strip_sizes;
    unsigned int strip_sizes_size;
    uint32_t *strip_offsets;
    unsigned int strip_offsets_size;
    uint8_t *yuv_line;
    unsigned int yuv_line_size;
    int rps;                                ///< row per strip
    uint8_t entries[TIFF_MAX_ENTRY * 12];   ///< entries in header
    int num_entries;                        ///< number of entries
    uint8_t **buf;                          ///< actual position in buffer
    uint8_t *buf_start;                     ///< pointer to first byte in buffer
    int buf_size;                           ///< buffer size
    uint16_t subsampling[2];                ///< YUV subsampling factors
    struct LZWEncodeState *lzws;            ///< LZW encode state
    uint32_t dpi;                           ///< image resolution in DPI
} TiffEncoderContext;

/**
 * Check free space in buffer.
 *
 * @param s Tiff context
 * @param need Needed bytes
 * @return 0 - ok, 1 - no free space
 */
static inline int check_size(TiffEncoderContext *s, uint64_t need)
{
    if (s->buf_size < *s->buf - s->buf_start + need) {
        *s->buf = s->buf_start + s->buf_size + 1;
        av_log(s->avctx, AV_LOG_ERROR, "Buffer is too small\n");
        return 1;
    }
    return 0;
}

/**
 * Put n values to buffer.
 *
 * @param p pointer to pointer to output buffer
 * @param n number of values
 * @param val pointer to values
 * @param type type of values
 * @param flip = 0 - normal copy, >0 - flip
 */
static void tnput(uint8_t **p, int n, const uint8_t *val, enum TiffTypes type,
                  int flip)
{
    int i;
#if HAVE_BIGENDIAN
    flip ^= ((int[]) { 0, 0, 0, 1, 3, 3 })[type];
#endif
    for (i = 0; i < n * type_sizes2[type]; i++)
        *(*p)++ = val[i ^ flip];
}

/**
 * Add entry to directory in tiff header.
 *
 * @param s Tiff context
 * @param tag tag that identifies the entry
 * @param type entry type
 * @param count the number of values
 * @param ptr_val pointer to values
 */
static void add_entry(TiffEncoderContext *s, enum TiffTags tag,
                      enum TiffTypes type, int count, const void *ptr_val)
{
    uint8_t *entries_ptr = s->entries + 12 * s->num_entries;

    av_assert0(s->num_entries < TIFF_MAX_ENTRY);

    bytestream_put_le16(&entries_ptr, tag);
    bytestream_put_le16(&entries_ptr, type);
    bytestream_put_le32(&entries_ptr, count);

    if (type_sizes[type] * (int64_t)count <= 4) {
        tnput(&entries_ptr, count, ptr_val, type, 0);
    } else {
        bytestream_put_le32(&entries_ptr, *s->buf - s->buf_start);
        check_size(s, count * (int64_t)type_sizes2[type]);
        tnput(s->buf, count, ptr_val, type, 0);
    }

    s->num_entries++;
}

static void add_entry1(TiffEncoderContext *s,
                       enum TiffTags tag, enum TiffTypes type, int val)
{
    uint16_t w  = val;
    uint32_t dw = val;
    add_entry(s, tag, type, 1, type == TIFF_SHORT ? (void *)&w : (void *)&dw);
}

/**
 * Encode one strip in tiff file.
 *
 * @param s Tiff context
 * @param src input buffer
 * @param dst output buffer
 * @param n size of input buffer
 * @param compr compression method
 * @return number of output bytes. If an output error is encountered, -1 is returned
 */
static int encode_strip(TiffEncoderContext *s, const int8_t *src,
                        uint8_t *dst, int n, int compr)
{
    switch (compr) {
#if CONFIG_ZLIB
    case TIFF_DEFLATE:
    case TIFF_ADOBE_DEFLATE:
    {
        unsigned long zlen = s->buf_size - (*s->buf - s->buf_start);
        if (compress(dst, &zlen, src, n) != Z_OK) {
            av_log(s->avctx, AV_LOG_ERROR, "Compressing failed\n");
            return -1;
        }
        return zlen;
    }
#endif
    case TIFF_RAW:
        if (check_size(s, n))
            return -1;
        memcpy(dst, src, n);
        return n;
    case TIFF_PACKBITS:
        return ff_rle_encode(dst, s->buf_size - (*s->buf - s->buf_start),
                             src, 1, n, 2, 0xff, -1, 0);
    case TIFF_LZW:
        return ff_lzw_encode(s->lzws, src, n);
    default:
        return -1;
    }
}

static void pack_yuv(TiffEncoderContext *s, const AVFrame *p,
                     uint8_t *dst, int lnum)
{
    int i, j, k;
    int w       = (s->width - 1) / s->subsampling[0] + 1;
    uint8_t *pu = &p->data[1][lnum / s->subsampling[1] * p->linesize[1]];
    uint8_t *pv = &p->data[2][lnum / s->subsampling[1] * p->linesize[2]];
    if (s->width % s->subsampling[0] || s->height % s->subsampling[1]) {
        for (i = 0; i < w; i++) {
            for (j = 0; j < s->subsampling[1]; j++)
                for (k = 0; k < s->subsampling[0]; k++)
                    *dst++ = p->data[0][FFMIN(lnum + j, s->height-1) * p->linesize[0] +
                                        FFMIN(i * s->subsampling[0] + k, s->width-1)];
            *dst++ = *pu++;
            *dst++ = *pv++;
        }
    }else{
        for (i = 0; i < w; i++) {
            for (j = 0; j < s->subsampling[1]; j++)
                for (k = 0; k < s->subsampling[0]; k++)
                    *dst++ = p->data[0][(lnum + j) * p->linesize[0] +
                                        i * s->subsampling[0] + k];
            *dst++ = *pu++;
            *dst++ = *pv++;
        }
    }
}

static int encode_frame(AVCodecContext *avctx, AVPacket *pkt,
                        const AVFrame *pict, int *got_packet)
{
    const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(avctx->pix_fmt);
    TiffEncoderContext *s = avctx->priv_data;
    const AVFrame *const p = pict;
    int i;
    uint8_t *ptr;
    uint8_t *offset;
    uint32_t strips;
    int bytes_per_row;
    uint32_t res[2] = { s->dpi, 1 };    // image resolution (72/1)
    uint16_t bpp_tab[4];
    int ret = 0;
    int is_yuv = 0, alpha = 0;
    int shift_h, shift_v;
    int packet_size;

    s->width          = avctx->width;
    s->height         = avctx->height;
    s->subsampling[0] = 1;
    s->subsampling[1] = 1;

    avctx->bits_per_coded_sample =
    s->bpp          = av_get_bits_per_pixel(desc);
    s->bpp_tab_size = desc->nb_components;

    switch (avctx->pix_fmt) {
    case AV_PIX_FMT_RGBA64LE:
    case AV_PIX_FMT_RGBA:
        alpha = 1;
    case AV_PIX_FMT_RGB48LE:
    case AV_PIX_FMT_RGB24:
        s->photometric_interpretation = TIFF_PHOTOMETRIC_RGB;
        break;
    case AV_PIX_FMT_GRAY8:
        avctx->bits_per_coded_sample = 0x28;
    case AV_PIX_FMT_GRAY8A:
        alpha = avctx->pix_fmt == AV_PIX_FMT_GRAY8A;
    case AV_PIX_FMT_GRAY16LE:
    case AV_PIX_FMT_MONOBLACK:
        s->photometric_interpretation = TIFF_PHOTOMETRIC_BLACK_IS_ZERO;
        break;
    case AV_PIX_FMT_PAL8:
        s->photometric_interpretation = TIFF_PHOTOMETRIC_PALETTE;
        break;
    case AV_PIX_FMT_MONOWHITE:
        s->photometric_interpretation = TIFF_PHOTOMETRIC_WHITE_IS_ZERO;
        break;
    case AV_PIX_FMT_YUV420P:
    case AV_PIX_FMT_YUV422P:
    case AV_PIX_FMT_YUV440P:
    case AV_PIX_FMT_YUV444P:
    case AV_PIX_FMT_YUV410P:
    case AV_PIX_FMT_YUV411P:
        av_pix_fmt_get_chroma_sub_sample(avctx->pix_fmt, &shift_h, &shift_v);
        s->photometric_interpretation = TIFF_PHOTOMETRIC_YCBCR;
        s->subsampling[0]             = 1 << shift_h;
        s->subsampling[1]             = 1 << shift_v;
        is_yuv                        = 1;
        break;
    default:
        av_log(s->avctx, AV_LOG_ERROR,
               "This colors format is not supported\n");
        return -1;
    }

    for (i = 0; i < s->bpp_tab_size; i++)
        bpp_tab[i] = desc->comp[i].depth_minus1 + 1;

    if (s->compr == TIFF_DEFLATE       ||
        s->compr == TIFF_ADOBE_DEFLATE ||
        s->compr == TIFF_LZW)
        // best choice for DEFLATE
        s->rps = s->height;
    else
        // suggest size of strip
        s->rps = FFMAX(8192 / (((s->width * s->bpp) >> 3) + 1), 1);
    // round rps up
    s->rps = ((s->rps - 1) / s->subsampling[1] + 1) * s->subsampling[1];

    strips = (s->height - 1) / s->rps + 1;

    bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
                     s->subsampling[0] * s->subsampling[1] + 7) >> 3;
    packet_size = avctx->height * bytes_per_row * 2 +
                  avctx->height * 4 + FF_MIN_BUFFER_SIZE;

    if ((ret = ff_alloc_packet2(avctx, pkt, packet_size)) < 0)
        return ret;
    ptr          = pkt->data;
    s->buf_start = pkt->data;
    s->buf       = &ptr;
    s->buf_size  = pkt->size;

    if (check_size(s, 8))
        goto fail;

    // write header
    bytestream_put_le16(&ptr, 0x4949);
    bytestream_put_le16(&ptr, 42);

    offset = ptr;
    bytestream_put_le32(&ptr, 0);

    if (strips > INT_MAX / FFMAX(sizeof(s->strip_sizes[0]), sizeof(s->strip_offsets[0]))) {
        ret = AVERROR(ENOMEM);
        goto fail;
    }
    av_fast_padded_mallocz(&s->strip_sizes  , &s->strip_sizes_size  , sizeof(s->strip_sizes  [0]) * strips);
    av_fast_padded_mallocz(&s->strip_offsets, &s->strip_offsets_size, sizeof(s->strip_offsets[0]) * strips);

    if (!s->strip_sizes || !s->strip_offsets) {
        ret = AVERROR(ENOMEM);
        goto fail;
    }

    if (is_yuv) {
        av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
        if (s->yuv_line == NULL) {
            av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
            ret = AVERROR(ENOMEM);
            goto fail;
        }
    }

#if CONFIG_ZLIB
    if (s->compr == TIFF_DEFLATE || s->compr == TIFF_ADOBE_DEFLATE) {
        uint8_t *zbuf;
        int zlen, zn;
        int j;

        zlen = bytes_per_row * s->rps;
        zbuf = av_malloc(zlen);
        if (!zbuf) {
            ret = AVERROR(ENOMEM);
            goto fail;
        }
        s->strip_offsets[0] = ptr - pkt->data;
        zn               = 0;
        for (j = 0; j < s->rps; j++) {
            if (is_yuv) {
                pack_yuv(s, p, s->yuv_line, j);
                memcpy(zbuf + zn, s->yuv_line, bytes_per_row);
                j += s->subsampling[1] - 1;
            } else
                memcpy(zbuf + j * bytes_per_row,
                       p->data[0] + j * p->linesize[0], bytes_per_row);
            zn += bytes_per_row;
        }
        ret = encode_strip(s, zbuf, ptr, zn, s->compr);
        av_free(zbuf);
        if (ret < 0) {
            av_log(s->avctx, AV_LOG_ERROR, "Encode strip failed\n");
            goto fail;
        }
        ptr           += ret;
        s->strip_sizes[0] = ptr - pkt->data - s->strip_offsets[0];
    } else
#endif
    {
    if (s->compr == TIFF_LZW) {
        s->lzws = av_malloc(ff_lzw_encode_state_size);
        if (!s->lzws) {
            ret = AVERROR(ENOMEM);
            goto fail;
        }
    }
    for (i = 0; i < s->height; i++) {
        if (s->strip_sizes[i / s->rps] == 0) {
            if (s->compr == TIFF_LZW) {
                ff_lzw_encode_init(s->lzws, ptr,
                                   s->buf_size - (*s->buf - s->buf_start),
                                   12, FF_LZW_TIFF, put_bits);
            }
            s->strip_offsets[i / s->rps] = ptr - pkt->data;
        }
        if (is_yuv) {
            pack_yuv(s, p, s->yuv_line, i);
            ret = encode_strip(s, s->yuv_line, ptr, bytes_per_row, s->compr);
            i  += s->subsampling[1] - 1;
        } else
            ret = encode_strip(s, p->data[0] + i * p->linesize[0],
                               ptr, bytes_per_row, s->compr);
        if (ret < 0) {
            av_log(s->avctx, AV_LOG_ERROR, "Encode strip failed\n");
            goto fail;
        }
        s->strip_sizes[i / s->rps] += ret;
        ptr                     += ret;
        if (s->compr == TIFF_LZW &&
            (i == s->height - 1 || i % s->rps == s->rps - 1)) {
            ret = ff_lzw_encode_flush(s->lzws, flush_put_bits);
            s->strip_sizes[(i / s->rps)] += ret;
            ptr                          += ret;
        }
    }
    if (s->compr == TIFF_LZW)
        av_freep(&s->lzws);
    }

    s->num_entries = 0;

    add_entry1(s, TIFF_SUBFILE, TIFF_LONG, 0);
    add_entry1(s, TIFF_WIDTH,   TIFF_LONG, s->width);
    add_entry1(s, TIFF_HEIGHT,  TIFF_LONG, s->height);

    if (s->bpp_tab_size)
        add_entry(s, TIFF_BPP, TIFF_SHORT, s->bpp_tab_size, bpp_tab);

    add_entry1(s, TIFF_COMPR,       TIFF_SHORT, s->compr);
    add_entry1(s, TIFF_PHOTOMETRIC, TIFF_SHORT, s->photometric_interpretation);
    add_entry(s,  TIFF_STRIP_OFFS,  TIFF_LONG,  strips, s->strip_offsets);

    if (s->bpp_tab_size)
        add_entry1(s, TIFF_SAMPLES_PER_PIXEL, TIFF_SHORT, s->bpp_tab_size);

    add_entry1(s, TIFF_ROWSPERSTRIP, TIFF_LONG,     s->rps);
    add_entry(s,  TIFF_STRIP_SIZE,   TIFF_LONG,     strips, s->strip_sizes);
    add_entry(s,  TIFF_XRES,         TIFF_RATIONAL, 1,      res);
    if (avctx->sample_aspect_ratio.num > 0 &&
        avctx->sample_aspect_ratio.den > 0) {
        AVRational y = av_mul_q(av_make_q(s->dpi, 1),
                                avctx->sample_aspect_ratio);
        res[0] = y.num;
        res[1] = y.den;
    }
    add_entry(s,  TIFF_YRES,         TIFF_RATIONAL, 1,      res);
    add_entry1(s, TIFF_RES_UNIT,     TIFF_SHORT,    2);

    if (!(avctx->flags & CODEC_FLAG_BITEXACT))
        add_entry(s, TIFF_SOFTWARE_NAME, TIFF_STRING,
                  strlen(LIBAVCODEC_IDENT) + 1, LIBAVCODEC_IDENT);

    if (avctx->pix_fmt == AV_PIX_FMT_PAL8) {
        uint16_t pal[256 * 3];
        for (i = 0; i < 256; i++) {
            uint32_t rgb = *(uint32_t *) (p->data[1] + i * 4);
            pal[i]       = ((rgb >> 16) & 0xff) * 257;
            pal[i + 256] = ((rgb >>  8) & 0xff) * 257;
            pal[i + 512] =  (rgb        & 0xff) * 257;
        }
        add_entry(s, TIFF_PAL, TIFF_SHORT, 256 * 3, pal);
    }
    if (alpha)
        add_entry1(s,TIFF_EXTRASAMPLES,      TIFF_SHORT,            2);
    if (is_yuv) {
        /** according to CCIR Recommendation 601.1 */
        uint32_t refbw[12] = { 15, 1, 235, 1, 128, 1, 240, 1, 128, 1, 240, 1 };
        add_entry(s, TIFF_YCBCR_SUBSAMPLING, TIFF_SHORT,    2, s->subsampling);
        if (avctx->chroma_sample_location == AVCHROMA_LOC_TOPLEFT)
            add_entry1(s, TIFF_YCBCR_POSITIONING, TIFF_SHORT, 2);
        add_entry(s, TIFF_REFERENCE_BW,      TIFF_RATIONAL, 6, refbw);
    }
    // write offset to dir
    bytestream_put_le32(&offset, ptr - pkt->data);

    if (check_size(s, 6 + s->num_entries * 12)) {
        ret = AVERROR(EINVAL);
        goto fail;
    }
    bytestream_put_le16(&ptr, s->num_entries);  // write tag count
    bytestream_put_buffer(&ptr, s->entries, s->num_entries * 12);
    bytestream_put_le32(&ptr, 0);

    pkt->size   = ptr - pkt->data;
    pkt->flags |= AV_PKT_FLAG_KEY;
    *got_packet = 1;

fail:
    return ret < 0 ? ret : 0;
}

static av_cold int encode_init(AVCodecContext *avctx)
{
    TiffEncoderContext *s = avctx->priv_data;

    avctx->coded_frame = av_frame_alloc();
    if (!avctx->coded_frame)
        return AVERROR(ENOMEM);

    avctx->coded_frame->pict_type = AV_PICTURE_TYPE_I;
    avctx->coded_frame->key_frame = 1;
    s->avctx = avctx;

    return 0;
}

static av_cold int encode_close(AVCodecContext *avctx)
{
    TiffEncoderContext *s = avctx->priv_data;

    av_frame_free(&avctx->coded_frame);
    av_freep(&s->strip_sizes);
    av_freep(&s->strip_offsets);
    av_freep(&s->yuv_line);

    return 0;
}

#define OFFSET(x) offsetof(TiffEncoderContext, x)
#define VE AV_OPT_FLAG_VIDEO_PARAM | AV_OPT_FLAG_ENCODING_PARAM
static const AVOption options[] = {
    {"dpi", "set the image resolution (in dpi)", OFFSET(dpi), AV_OPT_TYPE_INT, {.i64 = 72}, 1, 0x10000, AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_ENCODING_PARAM},
    { "compression_algo", NULL, OFFSET(compr), AV_OPT_TYPE_INT,   { .i64 = TIFF_PACKBITS }, TIFF_RAW, TIFF_DEFLATE, VE, "compression_algo" },
    { "packbits",         NULL, 0,             AV_OPT_TYPE_CONST, { .i64 = TIFF_PACKBITS }, 0,        0,            VE, "compression_algo" },
    { "raw",              NULL, 0,             AV_OPT_TYPE_CONST, { .i64 = TIFF_RAW      }, 0,        0,            VE, "compression_algo" },
    { "lzw",              NULL, 0,             AV_OPT_TYPE_CONST, { .i64 = TIFF_LZW      }, 0,        0,            VE, "compression_algo" },
#if CONFIG_ZLIB
    { "deflate",          NULL, 0,             AV_OPT_TYPE_CONST, { .i64 = TIFF_DEFLATE  }, 0,        0,            VE, "compression_algo" },
#endif
    { NULL },
};

static const AVClass tiffenc_class = {
    .class_name = "TIFF encoder",
    .item_name  = av_default_item_name,
    .option     = options,
    .version    = LIBAVUTIL_VERSION_INT,
};

AVCodec ff_tiff_encoder = {
    .name           = "tiff",
    .long_name      = NULL_IF_CONFIG_SMALL("TIFF image"),
    .type           = AVMEDIA_TYPE_VIDEO,
    .id             = AV_CODEC_ID_TIFF,
    .priv_data_size = sizeof(TiffEncoderContext),
    .init           = encode_init,
    .close          = encode_close,
    .capabilities   = CODEC_CAP_FRAME_THREADS | CODEC_CAP_INTRA_ONLY,
    .encode2        = encode_frame,
    .pix_fmts       = (const enum AVPixelFormat[]) {
        AV_PIX_FMT_RGB24, AV_PIX_FMT_PAL8, AV_PIX_FMT_GRAY8,
        AV_PIX_FMT_GRAY8A, AV_PIX_FMT_GRAY16LE,
        AV_PIX_FMT_MONOBLACK, AV_PIX_FMT_MONOWHITE,
        AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUV440P, AV_PIX_FMT_YUV444P,
        AV_PIX_FMT_YUV410P, AV_PIX_FMT_YUV411P, AV_PIX_FMT_RGB48LE,
        AV_PIX_FMT_RGBA, AV_PIX_FMT_RGBA64LE,
        AV_PIX_FMT_NONE
    },
    .priv_class     = &tiffenc_class,
};
Commit	Line	Data
2ba45a60 DM	1	/*
	2	* TIFF image encoder
	3	* Copyright (c) 2007 Bartlomiej Wolowiec
	4	*
	5	* This file is part of FFmpeg.
	6	*
	7	* FFmpeg is free software; you can redistribute it and/or
	8	* modify it under the terms of the GNU Lesser General Public
	9	* License as published by the Free Software Foundation; either
	10	* version 2.1 of the License, or (at your option) any later version.
	11	*
	12	* FFmpeg is distributed in the hope that it will be useful,
	13	* but WITHOUT ANY WARRANTY; without even the implied warranty of
	14	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
	15	* Lesser General Public License for more details.
	16	*
	17	* You should have received a copy of the GNU Lesser General Public
	18	* License along with FFmpeg; if not, write to the Free Software
	19	* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
	20	*/
	21
	22	/**
	23	* @file
	24	* TIFF image encoder
	25	* @author Bartlomiej Wolowiec
	26	*/
	27
	28	#include "config.h"
	29	#if CONFIG_ZLIB
	30	#include <zlib.h>
	31	#endif
	32
	33	#include "libavutil/imgutils.h"
	34	#include "libavutil/log.h"
	35	#include "libavutil/opt.h"
	36	#include "libavutil/pixdesc.h"
	37	#include "avcodec.h"
	38	#include "bytestream.h"
	39	#include "internal.h"
	40	#include "lzw.h"
	41	#include "put_bits.h"
	42	#include "rle.h"
	43	#include "tiff.h"
	44
	45	#define TIFF_MAX_ENTRY 32
	46
	47	/** sizes of various TIFF field types (string size = 1)*/
	48	static const uint8_t type_sizes2[14] = {
	49	0, 1, 1, 2, 4, 8, 1, 1, 2, 4, 8, 4, 8, 4
	50	};
	51
	52	typedef struct TiffEncoderContext {
	53	AVClass *class; ///< for private options
	54	AVCodecContext *avctx;
	55
	56	int width; ///< picture width
	57	int height; ///< picture height
	58	unsigned int bpp; ///< bits per pixel
	59	int compr; ///< compression level
	60	int bpp_tab_size; ///< bpp_tab size
	61	enum TiffPhotometric photometric_interpretation; ///< photometric interpretation
	62	int strips; ///< number of strips
	63	uint32_t *strip_sizes;
	64	unsigned int strip_sizes_size;
65	uint32_t *strip_offsets;
66	unsigned int strip_offsets_size;
67	uint8_t *yuv_line;
68	unsigned int yuv_line_size;
69	int rps; ///< row per strip
70	uint8_t entries[TIFF_MAX_ENTRY * 12]; ///< entries in header
71	int num_entries; ///< number of entries
72	uint8_t **buf; ///< actual position in buffer
73	uint8_t *buf_start; ///< pointer to first byte in buffer
74	int buf_size; ///< buffer size
75	uint16_t subsampling[2]; ///< YUV subsampling factors
76	struct LZWEncodeState *lzws; ///< LZW encode state
77	uint32_t dpi; ///< image resolution in DPI
78	} TiffEncoderContext;
79
80	/**
81	* Check free space in buffer.
82	*
83	* @param s Tiff context
84	* @param need Needed bytes
85	* @return 0 - ok, 1 - no free space
86	*/
87	static inline int check_size(TiffEncoderContext *s, uint64_t need)
88	{
89	if (s->buf_size < *s->buf - s->buf_start + need) {
90	*s->buf = s->buf_start + s->buf_size + 1;
91	av_log(s->avctx, AV_LOG_ERROR, "Buffer is too small\n");
92	return 1;
93	}
94	return 0;
95	}
96
97	/**
98	* Put n values to buffer.
99	*
100	* @param p pointer to pointer to output buffer
101	* @param n number of values
102	* @param val pointer to values
103	* @param type type of values
104	* @param flip = 0 - normal copy, >0 - flip
105	*/
106	static void tnput(uint8_t *p, int n, const uint8_t val, enum TiffTypes type,
107	int flip)
108	{
109	int i;
110	#if HAVE_BIGENDIAN
111	flip ^= ((int[]) { 0, 0, 0, 1, 3, 3 })[type];
112	#endif
113	for (i = 0; i < n * type_sizes2[type]; i++)
114	(p)++ = val[i ^ flip];
115	}
116
117	/**
118	* Add entry to directory in tiff header.
119	*
120	* @param s Tiff context
121	* @param tag tag that identifies the entry
122	* @param type entry type
123	* @param count the number of values
124	* @param ptr_val pointer to values
125	*/
126	static void add_entry(TiffEncoderContext *s, enum TiffTags tag,
127	enum TiffTypes type, int count, const void *ptr_val)
128	{
129	uint8_t entries_ptr = s->entries + 12 s->num_entries;
130
131	av_assert0(s->num_entries < TIFF_MAX_ENTRY);
132
133	bytestream_put_le16(&entries_ptr, tag);
134	bytestream_put_le16(&entries_ptr, type);
135	bytestream_put_le32(&entries_ptr, count);
136
137	if (type_sizes[type] * (int64_t)count <= 4) {
138	tnput(&entries_ptr, count, ptr_val, type, 0);
139	} else {
140	bytestream_put_le32(&entries_ptr, *s->buf - s->buf_start);
141	check_size(s, count * (int64_t)type_sizes2[type]);
142	tnput(s->buf, count, ptr_val, type, 0);
143	}
144
145	s->num_entries++;
146	}
147
148	static void add_entry1(TiffEncoderContext *s,
149	enum TiffTags tag, enum TiffTypes type, int val)
150	{
151	uint16_t w = val;
152	uint32_t dw = val;
153	add_entry(s, tag, type, 1, type == TIFF_SHORT ? (void )&w : (void )&dw);
154	}
155
156	/**
157	* Encode one strip in tiff file.
158	*
159	* @param s Tiff context
160	* @param src input buffer
161	* @param dst output buffer
162	* @param n size of input buffer
163	* @param compr compression method
164	* @return number of output bytes. If an output error is encountered, -1 is returned
165	*/
166	static int encode_strip(TiffEncoderContext s, const int8_t src,
167	uint8_t *dst, int n, int compr)
168	{
169	switch (compr) {
170	#if CONFIG_ZLIB
171	case TIFF_DEFLATE:
172	case TIFF_ADOBE_DEFLATE:
173	{
174	unsigned long zlen = s->buf_size - (*s->buf - s->buf_start);
175	if (compress(dst, &zlen, src, n) != Z_OK) {
176	av_log(s->avctx, AV_LOG_ERROR, "Compressing failed\n");
177	return -1;
178	}
179	return zlen;
180	}
181	#endif
182	case TIFF_RAW:
183	if (check_size(s, n))
184	return -1;
185	memcpy(dst, src, n);
186	return n;
187	case TIFF_PACKBITS:
188	return ff_rle_encode(dst, s->buf_size - (*s->buf - s->buf_start),
189	src, 1, n, 2, 0xff, -1, 0);
190	case TIFF_LZW:
191	return ff_lzw_encode(s->lzws, src, n);
192	default:
193	return -1;
194	}
195	}
196
197	static void pack_yuv(TiffEncoderContext s, const AVFrame p,
198	uint8_t *dst, int lnum)
199	{
200	int i, j, k;
201	int w = (s->width - 1) / s->subsampling[0] + 1;
202	uint8_t pu = &p->data[1][lnum / s->subsampling[1] p->linesize[1]];
203	uint8_t pv = &p->data[2][lnum / s->subsampling[1] p->linesize[2]];
204	if (s->width % s->subsampling[0] \|\| s->height % s->subsampling[1]) {
205	for (i = 0; i < w; i++) {
206	for (j = 0; j < s->subsampling[1]; j++)
207	for (k = 0; k < s->subsampling[0]; k++)
208	dst++ = p->data[0][FFMIN(lnum + j, s->height-1) p->linesize[0] +
209	FFMIN(i * s->subsampling[0] + k, s->width-1)];
210	dst++ = pu++;
211	dst++ = pv++;
212	}
213	}else{
214	for (i = 0; i < w; i++) {
215	for (j = 0; j < s->subsampling[1]; j++)
216	for (k = 0; k < s->subsampling[0]; k++)
217	dst++ = p->data[0][(lnum + j) p->linesize[0] +
218	i * s->subsampling[0] + k];
219	dst++ = pu++;
220	dst++ = pv++;
221	}
222	}
223	}
224
225	static int encode_frame(AVCodecContext avctx, AVPacket pkt,
226	const AVFrame pict, int got_packet)
227	{
228	const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(avctx->pix_fmt);
229	TiffEncoderContext *s = avctx->priv_data;
230	const AVFrame *const p = pict;
231	int i;
232	uint8_t *ptr;
233	uint8_t *offset;
234	uint32_t strips;
235	int bytes_per_row;
236	uint32_t res[2] = { s->dpi, 1 }; // image resolution (72/1)
237	uint16_t bpp_tab[4];
f6fa7814	238	int ret = 0;
2ba45a60 DM	239	int is_yuv = 0, alpha = 0;
	240	int shift_h, shift_v;
	241	int packet_size;
	242
	243	s->width = avctx->width;
	244	s->height = avctx->height;
	245	s->subsampling[0] = 1;
	246	s->subsampling[1] = 1;
	247
	248	avctx->bits_per_coded_sample =
	249	s->bpp = av_get_bits_per_pixel(desc);
	250	s->bpp_tab_size = desc->nb_components;
	251
	252	switch (avctx->pix_fmt) {
	253	case AV_PIX_FMT_RGBA64LE:
	254	case AV_PIX_FMT_RGBA:
	255	alpha = 1;
	256	case AV_PIX_FMT_RGB48LE:
	257	case AV_PIX_FMT_RGB24:
	258	s->photometric_interpretation = TIFF_PHOTOMETRIC_RGB;
	259	break;
	260	case AV_PIX_FMT_GRAY8:
	261	avctx->bits_per_coded_sample = 0x28;
	262	case AV_PIX_FMT_GRAY8A:
	263	alpha = avctx->pix_fmt == AV_PIX_FMT_GRAY8A;
	264	case AV_PIX_FMT_GRAY16LE:
	265	case AV_PIX_FMT_MONOBLACK:
	266	s->photometric_interpretation = TIFF_PHOTOMETRIC_BLACK_IS_ZERO;
	267	break;
	268	case AV_PIX_FMT_PAL8:
	269	s->photometric_interpretation = TIFF_PHOTOMETRIC_PALETTE;
	270	break;
	271	case AV_PIX_FMT_MONOWHITE:
	272	s->photometric_interpretation = TIFF_PHOTOMETRIC_WHITE_IS_ZERO;
	273	break;
	274	case AV_PIX_FMT_YUV420P:
	275	case AV_PIX_FMT_YUV422P:
	276	case AV_PIX_FMT_YUV440P:
	277	case AV_PIX_FMT_YUV444P:
	278	case AV_PIX_FMT_YUV410P:
	279	case AV_PIX_FMT_YUV411P:
	280	av_pix_fmt_get_chroma_sub_sample(avctx->pix_fmt, &shift_h, &shift_v);
	281	s->photometric_interpretation = TIFF_PHOTOMETRIC_YCBCR;
	282	s->subsampling[0] = 1 << shift_h;
	283	s->subsampling[1] = 1 << shift_v;
	284	is_yuv = 1;
	285	break;
	286	default:
	287	av_log(s->avctx, AV_LOG_ERROR,
	288	"This colors format is not supported\n");
	289	return -1;
	290	}
	291
	292	for (i = 0; i < s->bpp_tab_size; i++)
	293	bpp_tab[i] = desc->comp[i].depth_minus1 + 1;
	294
	295	if (s->compr == TIFF_DEFLATE \|\|
	296	s->compr == TIFF_ADOBE_DEFLATE \|\|
	297	s->compr == TIFF_LZW)
	298	// best choice for DEFLATE
	299	s->rps = s->height;
	300	else
	301	// suggest size of strip
	302	s->rps = FFMAX(8192 / (((s->width * s->bpp) >> 3) + 1), 1);
303	// round rps up
304	s->rps = ((s->rps - 1) / s->subsampling[1] + 1) * s->subsampling[1];
305
306	strips = (s->height - 1) / s->rps + 1;
307
308	bytes_per_row = (((s->width - 1) / s->subsampling[0] + 1) * s->bpp *
309	s->subsampling[0] * s->subsampling[1] + 7) >> 3;
310	packet_size = avctx->height * bytes_per_row * 2 +
311	avctx->height * 4 + FF_MIN_BUFFER_SIZE;
312
313	if ((ret = ff_alloc_packet2(avctx, pkt, packet_size)) < 0)
314	return ret;
315	ptr = pkt->data;
316	s->buf_start = pkt->data;
317	s->buf = &ptr;
318	s->buf_size = pkt->size;
319
320	if (check_size(s, 8))
321	goto fail;
322
323	// write header
324	bytestream_put_le16(&ptr, 0x4949);
325	bytestream_put_le16(&ptr, 42);
326
327	offset = ptr;
328	bytestream_put_le32(&ptr, 0);
329
f6fa7814 DM	330	if (strips > INT_MAX / FFMAX(sizeof(s->strip_sizes[0]), sizeof(s->strip_offsets[0]))) {
	331	ret = AVERROR(ENOMEM);
	332	goto fail;
	333	}
2ba45a60 DM	334	av_fast_padded_mallocz(&s->strip_sizes , &s->strip_sizes_size , sizeof(s->strip_sizes [0]) * strips);
	335	av_fast_padded_mallocz(&s->strip_offsets, &s->strip_offsets_size, sizeof(s->strip_offsets[0]) * strips);
	336
	337	if (!s->strip_sizes \|\| !s->strip_offsets) {
	338	ret = AVERROR(ENOMEM);
	339	goto fail;
	340	}
	341
	342	if (is_yuv) {
	343	av_fast_padded_malloc(&s->yuv_line, &s->yuv_line_size, bytes_per_row);
	344	if (s->yuv_line == NULL) {
	345	av_log(s->avctx, AV_LOG_ERROR, "Not enough memory\n");
	346	ret = AVERROR(ENOMEM);
	347	goto fail;
	348	}
	349	}
	350
	351	#if CONFIG_ZLIB
	352	if (s->compr == TIFF_DEFLATE \|\| s->compr == TIFF_ADOBE_DEFLATE) {
	353	uint8_t *zbuf;
	354	int zlen, zn;
	355	int j;
	356
	357	zlen = bytes_per_row * s->rps;
	358	zbuf = av_malloc(zlen);
	359	if (!zbuf) {
	360	ret = AVERROR(ENOMEM);
	361	goto fail;
	362	}
	363	s->strip_offsets[0] = ptr - pkt->data;
	364	zn = 0;
	365	for (j = 0; j < s->rps; j++) {
	366	if (is_yuv) {
	367	pack_yuv(s, p, s->yuv_line, j);
	368	memcpy(zbuf + zn, s->yuv_line, bytes_per_row);
	369	j += s->subsampling[1] - 1;
	370	} else
	371	memcpy(zbuf + j * bytes_per_row,
	372	p->data[0] + j * p->linesize[0], bytes_per_row);
	373	zn += bytes_per_row;
	374	}
	375	ret = encode_strip(s, zbuf, ptr, zn, s->compr);
	376	av_free(zbuf);
	377	if (ret < 0) {
	378	av_log(s->avctx, AV_LOG_ERROR, "Encode strip failed\n");
	379	goto fail;
	380	}
	381	ptr += ret;
	382	s->strip_sizes[0] = ptr - pkt->data - s->strip_offsets[0];
	383	} else
	384	#endif
	385	{
	386	if (s->compr == TIFF_LZW) {
	387	s->lzws = av_malloc(ff_lzw_encode_state_size);
	388	if (!s->lzws) {
	389	ret = AVERROR(ENOMEM);
	390	goto fail;
	391	}
	392	}
	393	for (i = 0; i < s->height; i++) {
	394	if (s->strip_sizes[i / s->rps] == 0) {
	395	if (s->compr == TIFF_LZW) {
	396	ff_lzw_encode_init(s->lzws, ptr,
	397	s->buf_size - (*s->buf - s->buf_start),
398	12, FF_LZW_TIFF, put_bits);
399	}
400	s->strip_offsets[i / s->rps] = ptr - pkt->data;
401	}
402	if (is_yuv) {
403	pack_yuv(s, p, s->yuv_line, i);
404	ret = encode_strip(s, s->yuv_line, ptr, bytes_per_row, s->compr);
405	i += s->subsampling[1] - 1;
406	} else
407	ret = encode_strip(s, p->data[0] + i * p->linesize[0],
408	ptr, bytes_per_row, s->compr);
409	if (ret < 0) {
410	av_log(s->avctx, AV_LOG_ERROR, "Encode strip failed\n");
411	goto fail;
412	}
413	s->strip_sizes[i / s->rps] += ret;
414	ptr += ret;
415	if (s->compr == TIFF_LZW &&
416	(i == s->height - 1 \|\| i % s->rps == s->rps - 1)) {
417	ret = ff_lzw_encode_flush(s->lzws, flush_put_bits);
418	s->strip_sizes[(i / s->rps)] += ret;
419	ptr += ret;
420	}
421	}
422	if (s->compr == TIFF_LZW)
f6fa7814	423	av_freep(&s->lzws);
2ba45a60 DM	424	}
	425
	426	s->num_entries = 0;
	427
	428	add_entry1(s, TIFF_SUBFILE, TIFF_LONG, 0);
	429	add_entry1(s, TIFF_WIDTH, TIFF_LONG, s->width);
	430	add_entry1(s, TIFF_HEIGHT, TIFF_LONG, s->height);
	431
	432	if (s->bpp_tab_size)
	433	add_entry(s, TIFF_BPP, TIFF_SHORT, s->bpp_tab_size, bpp_tab);
	434
	435	add_entry1(s, TIFF_COMPR, TIFF_SHORT, s->compr);
	436	add_entry1(s, TIFF_PHOTOMETRIC, TIFF_SHORT, s->photometric_interpretation);
	437	add_entry(s, TIFF_STRIP_OFFS, TIFF_LONG, strips, s->strip_offsets);
	438
	439	if (s->bpp_tab_size)
	440	add_entry1(s, TIFF_SAMPLES_PER_PIXEL, TIFF_SHORT, s->bpp_tab_size);
	441
	442	add_entry1(s, TIFF_ROWSPERSTRIP, TIFF_LONG, s->rps);
	443	add_entry(s, TIFF_STRIP_SIZE, TIFF_LONG, strips, s->strip_sizes);
	444	add_entry(s, TIFF_XRES, TIFF_RATIONAL, 1, res);
	445	if (avctx->sample_aspect_ratio.num > 0 &&
	446	avctx->sample_aspect_ratio.den > 0) {
	447	AVRational y = av_mul_q(av_make_q(s->dpi, 1),
	448	avctx->sample_aspect_ratio);
	449	res[0] = y.num;
	450	res[1] = y.den;
	451	}
	452	add_entry(s, TIFF_YRES, TIFF_RATIONAL, 1, res);
	453	add_entry1(s, TIFF_RES_UNIT, TIFF_SHORT, 2);
	454
	455	if (!(avctx->flags & CODEC_FLAG_BITEXACT))
	456	add_entry(s, TIFF_SOFTWARE_NAME, TIFF_STRING,
	457	strlen(LIBAVCODEC_IDENT) + 1, LIBAVCODEC_IDENT);
	458
	459	if (avctx->pix_fmt == AV_PIX_FMT_PAL8) {
	460	uint16_t pal[256 * 3];
	461	for (i = 0; i < 256; i++) {
	462	uint32_t rgb = (uint32_t ) (p->data[1] + i * 4);
	463	pal[i] = ((rgb >> 16) & 0xff) * 257;
	464	pal[i + 256] = ((rgb >> 8) & 0xff) * 257;
	465	pal[i + 512] = (rgb & 0xff) * 257;
	466	}
	467	add_entry(s, TIFF_PAL, TIFF_SHORT, 256 * 3, pal);
	468	}
	469	if (alpha)
	470	add_entry1(s,TIFF_EXTRASAMPLES, TIFF_SHORT, 2);
	471	if (is_yuv) {
	472	/** according to CCIR Recommendation 601.1 */
	473	uint32_t refbw[12] = { 15, 1, 235, 1, 128, 1, 240, 1, 128, 1, 240, 1 };
	474	add_entry(s, TIFF_YCBCR_SUBSAMPLING, TIFF_SHORT, 2, s->subsampling);
	475	if (avctx->chroma_sample_location == AVCHROMA_LOC_TOPLEFT)
	476	add_entry1(s, TIFF_YCBCR_POSITIONING, TIFF_SHORT, 2);
	477	add_entry(s, TIFF_REFERENCE_BW, TIFF_RATIONAL, 6, refbw);
	478	}
	479	// write offset to dir
	480	bytestream_put_le32(&offset, ptr - pkt->data);
	481
	482	if (check_size(s, 6 + s->num_entries * 12)) {
	483	ret = AVERROR(EINVAL);
	484	goto fail;
	485	}
	486	bytestream_put_le16(&ptr, s->num_entries); // write tag count
	487	bytestream_put_buffer(&ptr, s->entries, s->num_entries * 12);
488	bytestream_put_le32(&ptr, 0);
489
490	pkt->size = ptr - pkt->data;
491	pkt->flags \|= AV_PKT_FLAG_KEY;
492	*got_packet = 1;
493
494	fail:
495	return ret < 0 ? ret : 0;
496	}
497
498	static av_cold int encode_init(AVCodecContext *avctx)
499	{
500	TiffEncoderContext *s = avctx->priv_data;
501
502	avctx->coded_frame = av_frame_alloc();
503	if (!avctx->coded_frame)
504	return AVERROR(ENOMEM);
505
506	avctx->coded_frame->pict_type = AV_PICTURE_TYPE_I;
507	avctx->coded_frame->key_frame = 1;
508	s->avctx = avctx;
509
510	return 0;
511	}
512
513	static av_cold int encode_close(AVCodecContext *avctx)
514	{
515	TiffEncoderContext *s = avctx->priv_data;
516
517	av_frame_free(&avctx->coded_frame);
518	av_freep(&s->strip_sizes);
519	av_freep(&s->strip_offsets);
520	av_freep(&s->yuv_line);
521
522	return 0;
523	}
524
525	#define OFFSET(x) offsetof(TiffEncoderContext, x)
526	#define VE AV_OPT_FLAG_VIDEO_PARAM \| AV_OPT_FLAG_ENCODING_PARAM
527	static const AVOption options[] = {
528	{"dpi", "set the image resolution (in dpi)", OFFSET(dpi), AV_OPT_TYPE_INT, {.i64 = 72}, 1, 0x10000, AV_OPT_FLAG_VIDEO_PARAM\|AV_OPT_FLAG_ENCODING_PARAM},
529	{ "compression_algo", NULL, OFFSET(compr), AV_OPT_TYPE_INT, { .i64 = TIFF_PACKBITS }, TIFF_RAW, TIFF_DEFLATE, VE, "compression_algo" },
530	{ "packbits", NULL, 0, AV_OPT_TYPE_CONST, { .i64 = TIFF_PACKBITS }, 0, 0, VE, "compression_algo" },
531	{ "raw", NULL, 0, AV_OPT_TYPE_CONST, { .i64 = TIFF_RAW }, 0, 0, VE, "compression_algo" },
532	{ "lzw", NULL, 0, AV_OPT_TYPE_CONST, { .i64 = TIFF_LZW }, 0, 0, VE, "compression_algo" },
533	#if CONFIG_ZLIB
534	{ "deflate", NULL, 0, AV_OPT_TYPE_CONST, { .i64 = TIFF_DEFLATE }, 0, 0, VE, "compression_algo" },
535	#endif
536	{ NULL },
537	};
538
539	static const AVClass tiffenc_class = {
540	.class_name = "TIFF encoder",
541	.item_name = av_default_item_name,
542	.option = options,
543	.version = LIBAVUTIL_VERSION_INT,
544	};
545
546	AVCodec ff_tiff_encoder = {
547	.name = "tiff",
548	.long_name = NULL_IF_CONFIG_SMALL("TIFF image"),
549	.type = AVMEDIA_TYPE_VIDEO,
550	.id = AV_CODEC_ID_TIFF,
551	.priv_data_size = sizeof(TiffEncoderContext),
552	.init = encode_init,
553	.close = encode_close,
554	.capabilities = CODEC_CAP_FRAME_THREADS \| CODEC_CAP_INTRA_ONLY,
555	.encode2 = encode_frame,
556	.pix_fmts = (const enum AVPixelFormat[]) {
557	AV_PIX_FMT_RGB24, AV_PIX_FMT_PAL8, AV_PIX_FMT_GRAY8,
558	AV_PIX_FMT_GRAY8A, AV_PIX_FMT_GRAY16LE,
559	AV_PIX_FMT_MONOBLACK, AV_PIX_FMT_MONOWHITE,
560	AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUV440P, AV_PIX_FMT_YUV444P,
561	AV_PIX_FMT_YUV410P, AV_PIX_FMT_YUV411P, AV_PIX_FMT_RGB48LE,
562	AV_PIX_FMT_RGBA, AV_PIX_FMT_RGBA64LE,
563	AV_PIX_FMT_NONE
564	},
565	.priv_class = &tiffenc_class,
566	};