348 lines
11 KiB
C
348 lines
11 KiB
C
/*
|
|
* Copyright (c) 2017 Vittorio Giovara <vittorio.giovara@gmail.com>
|
|
*
|
|
* 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
|
|
* tonemap algorithms
|
|
*/
|
|
|
|
#include <float.h>
|
|
#include <stdio.h>
|
|
#include <string.h>
|
|
|
|
#include "libavutil/imgutils.h"
|
|
#include "libavutil/internal.h"
|
|
#include "libavutil/intreadwrite.h"
|
|
#include "libavutil/opt.h"
|
|
#include "libavutil/pixdesc.h"
|
|
|
|
#include "avfilter.h"
|
|
#include "colorspace.h"
|
|
#include "formats.h"
|
|
#include "internal.h"
|
|
#include "video.h"
|
|
|
|
enum TonemapAlgorithm {
|
|
TONEMAP_NONE,
|
|
TONEMAP_LINEAR,
|
|
TONEMAP_GAMMA,
|
|
TONEMAP_CLIP,
|
|
TONEMAP_REINHARD,
|
|
TONEMAP_HABLE,
|
|
TONEMAP_MOBIUS,
|
|
TONEMAP_MAX,
|
|
};
|
|
|
|
static const struct LumaCoefficients luma_coefficients[AVCOL_SPC_NB] = {
|
|
[AVCOL_SPC_FCC] = { 0.30, 0.59, 0.11 },
|
|
[AVCOL_SPC_BT470BG] = { 0.299, 0.587, 0.114 },
|
|
[AVCOL_SPC_SMPTE170M] = { 0.299, 0.587, 0.114 },
|
|
[AVCOL_SPC_BT709] = { 0.2126, 0.7152, 0.0722 },
|
|
[AVCOL_SPC_SMPTE240M] = { 0.212, 0.701, 0.087 },
|
|
[AVCOL_SPC_BT2020_NCL] = { 0.2627, 0.6780, 0.0593 },
|
|
[AVCOL_SPC_BT2020_CL] = { 0.2627, 0.6780, 0.0593 },
|
|
};
|
|
|
|
typedef struct TonemapContext {
|
|
const AVClass *class;
|
|
|
|
enum TonemapAlgorithm tonemap;
|
|
double param;
|
|
double desat;
|
|
double peak;
|
|
|
|
const struct LumaCoefficients *coeffs;
|
|
} TonemapContext;
|
|
|
|
static const enum AVPixelFormat pix_fmts[] = {
|
|
AV_PIX_FMT_GBRPF32,
|
|
AV_PIX_FMT_GBRAPF32,
|
|
AV_PIX_FMT_NONE,
|
|
};
|
|
|
|
static int query_formats(AVFilterContext *ctx)
|
|
{
|
|
return ff_set_common_formats(ctx, ff_make_format_list(pix_fmts));
|
|
}
|
|
|
|
static av_cold int init(AVFilterContext *ctx)
|
|
{
|
|
TonemapContext *s = ctx->priv;
|
|
|
|
switch(s->tonemap) {
|
|
case TONEMAP_GAMMA:
|
|
if (isnan(s->param))
|
|
s->param = 1.8f;
|
|
break;
|
|
case TONEMAP_REINHARD:
|
|
if (!isnan(s->param))
|
|
s->param = (1.0f - s->param) / s->param;
|
|
break;
|
|
case TONEMAP_MOBIUS:
|
|
if (isnan(s->param))
|
|
s->param = 0.3f;
|
|
break;
|
|
}
|
|
|
|
if (isnan(s->param))
|
|
s->param = 1.0f;
|
|
|
|
return 0;
|
|
}
|
|
|
|
static float hable(float in)
|
|
{
|
|
float a = 0.15f, b = 0.50f, c = 0.10f, d = 0.20f, e = 0.02f, f = 0.30f;
|
|
return (in * (in * a + b * c) + d * e) / (in * (in * a + b) + d * f) - e / f;
|
|
}
|
|
|
|
static float mobius(float in, float j, double peak)
|
|
{
|
|
float a, b;
|
|
|
|
if (in <= j)
|
|
return in;
|
|
|
|
a = -j * j * (peak - 1.0f) / (j * j - 2.0f * j + peak);
|
|
b = (j * j - 2.0f * j * peak + peak) / FFMAX(peak - 1.0f, 1e-6);
|
|
|
|
return (b * b + 2.0f * b * j + j * j) / (b - a) * (in + a) / (in + b);
|
|
}
|
|
|
|
#define MIX(x,y,a) (x) * (1 - (a)) + (y) * (a)
|
|
static void tonemap(TonemapContext *s, AVFrame *out, const AVFrame *in,
|
|
const AVPixFmtDescriptor *desc, int x, int y, double peak)
|
|
{
|
|
const float *r_in = (const float *)(in->data[0] + x * desc->comp[0].step + y * in->linesize[0]);
|
|
const float *b_in = (const float *)(in->data[1] + x * desc->comp[1].step + y * in->linesize[1]);
|
|
const float *g_in = (const float *)(in->data[2] + x * desc->comp[2].step + y * in->linesize[2]);
|
|
float *r_out = (float *)(out->data[0] + x * desc->comp[0].step + y * out->linesize[0]);
|
|
float *b_out = (float *)(out->data[1] + x * desc->comp[1].step + y * out->linesize[1]);
|
|
float *g_out = (float *)(out->data[2] + x * desc->comp[2].step + y * out->linesize[2]);
|
|
float sig, sig_orig;
|
|
|
|
/* load values */
|
|
*r_out = *r_in;
|
|
*b_out = *b_in;
|
|
*g_out = *g_in;
|
|
|
|
/* desaturate to prevent unnatural colors */
|
|
if (s->desat > 0) {
|
|
float luma = s->coeffs->cr * *r_in + s->coeffs->cg * *g_in + s->coeffs->cb * *b_in;
|
|
float overbright = FFMAX(luma - s->desat, 1e-6) / FFMAX(luma, 1e-6);
|
|
*r_out = MIX(*r_in, luma, overbright);
|
|
*g_out = MIX(*g_in, luma, overbright);
|
|
*b_out = MIX(*b_in, luma, overbright);
|
|
}
|
|
|
|
/* pick the brightest component, reducing the value range as necessary
|
|
* to keep the entire signal in range and preventing discoloration due to
|
|
* out-of-bounds clipping */
|
|
sig = FFMAX(FFMAX3(*r_out, *g_out, *b_out), 1e-6);
|
|
sig_orig = sig;
|
|
|
|
switch(s->tonemap) {
|
|
default:
|
|
case TONEMAP_NONE:
|
|
// do nothing
|
|
break;
|
|
case TONEMAP_LINEAR:
|
|
sig = sig * s->param / peak;
|
|
break;
|
|
case TONEMAP_GAMMA:
|
|
sig = sig > 0.05f ? pow(sig / peak, 1.0f / s->param)
|
|
: sig * pow(0.05f / peak, 1.0f / s->param) / 0.05f;
|
|
break;
|
|
case TONEMAP_CLIP:
|
|
sig = av_clipf(sig * s->param, 0, 1.0f);
|
|
break;
|
|
case TONEMAP_HABLE:
|
|
sig = hable(sig) / hable(peak);
|
|
break;
|
|
case TONEMAP_REINHARD:
|
|
sig = sig / (sig + s->param) * (peak + s->param) / peak;
|
|
break;
|
|
case TONEMAP_MOBIUS:
|
|
sig = mobius(sig, s->param, peak);
|
|
break;
|
|
}
|
|
|
|
/* apply the computed scale factor to the color,
|
|
* linearly to prevent discoloration */
|
|
*r_out *= sig / sig_orig;
|
|
*g_out *= sig / sig_orig;
|
|
*b_out *= sig / sig_orig;
|
|
}
|
|
|
|
typedef struct ThreadData {
|
|
AVFrame *in, *out;
|
|
const AVPixFmtDescriptor *desc;
|
|
double peak;
|
|
} ThreadData;
|
|
|
|
static int tonemap_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
|
|
{
|
|
TonemapContext *s = ctx->priv;
|
|
ThreadData *td = arg;
|
|
AVFrame *in = td->in;
|
|
AVFrame *out = td->out;
|
|
const AVPixFmtDescriptor *desc = td->desc;
|
|
const int slice_start = (in->height * jobnr) / nb_jobs;
|
|
const int slice_end = (in->height * (jobnr+1)) / nb_jobs;
|
|
double peak = td->peak;
|
|
|
|
for (int y = slice_start; y < slice_end; y++)
|
|
for (int x = 0; x < out->width; x++)
|
|
tonemap(s, out, in, desc, x, y, peak);
|
|
|
|
return 0;
|
|
}
|
|
|
|
static int filter_frame(AVFilterLink *link, AVFrame *in)
|
|
{
|
|
AVFilterContext *ctx = link->dst;
|
|
TonemapContext *s = ctx->priv;
|
|
AVFilterLink *outlink = ctx->outputs[0];
|
|
ThreadData td;
|
|
AVFrame *out;
|
|
const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(link->format);
|
|
const AVPixFmtDescriptor *odesc = av_pix_fmt_desc_get(outlink->format);
|
|
int ret, x, y;
|
|
double peak = s->peak;
|
|
|
|
if (!desc || !odesc) {
|
|
av_frame_free(&in);
|
|
return AVERROR_BUG;
|
|
}
|
|
|
|
out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
|
|
if (!out) {
|
|
av_frame_free(&in);
|
|
return AVERROR(ENOMEM);
|
|
}
|
|
|
|
ret = av_frame_copy_props(out, in);
|
|
if (ret < 0) {
|
|
av_frame_free(&in);
|
|
av_frame_free(&out);
|
|
return ret;
|
|
}
|
|
|
|
/* input and output transfer will be linear */
|
|
if (in->color_trc == AVCOL_TRC_UNSPECIFIED) {
|
|
av_log(s, AV_LOG_WARNING, "Untagged transfer, assuming linear light\n");
|
|
out->color_trc = AVCOL_TRC_LINEAR;
|
|
} else if (in->color_trc != AVCOL_TRC_LINEAR)
|
|
av_log(s, AV_LOG_WARNING, "Tonemapping works on linear light only\n");
|
|
|
|
/* read peak from side data if not passed in */
|
|
if (!peak) {
|
|
peak = ff_determine_signal_peak(in);
|
|
av_log(s, AV_LOG_DEBUG, "Computed signal peak: %f\n", peak);
|
|
}
|
|
|
|
/* load original color space even if pixel format is RGB to compute overbrights */
|
|
s->coeffs = &luma_coefficients[in->colorspace];
|
|
if (s->desat > 0 && (in->colorspace == AVCOL_SPC_UNSPECIFIED || !s->coeffs)) {
|
|
if (in->colorspace == AVCOL_SPC_UNSPECIFIED)
|
|
av_log(s, AV_LOG_WARNING, "Missing color space information, ");
|
|
else if (!s->coeffs)
|
|
av_log(s, AV_LOG_WARNING, "Unsupported color space '%s', ",
|
|
av_color_space_name(in->colorspace));
|
|
av_log(s, AV_LOG_WARNING, "desaturation is disabled\n");
|
|
s->desat = 0;
|
|
}
|
|
|
|
/* do the tone map */
|
|
td.out = out;
|
|
td.in = in;
|
|
td.desc = desc;
|
|
td.peak = peak;
|
|
ctx->internal->execute(ctx, tonemap_slice, &td, NULL, FFMIN(in->height, ff_filter_get_nb_threads(ctx)));
|
|
|
|
/* copy/generate alpha if needed */
|
|
if (desc->flags & AV_PIX_FMT_FLAG_ALPHA && odesc->flags & AV_PIX_FMT_FLAG_ALPHA) {
|
|
av_image_copy_plane(out->data[3], out->linesize[3],
|
|
in->data[3], in->linesize[3],
|
|
out->linesize[3], outlink->h);
|
|
} else if (odesc->flags & AV_PIX_FMT_FLAG_ALPHA) {
|
|
for (y = 0; y < out->height; y++) {
|
|
for (x = 0; x < out->width; x++) {
|
|
AV_WN32(out->data[3] + x * odesc->comp[3].step + y * out->linesize[3],
|
|
av_float2int(1.0f));
|
|
}
|
|
}
|
|
}
|
|
|
|
av_frame_free(&in);
|
|
|
|
ff_update_hdr_metadata(out, peak);
|
|
|
|
return ff_filter_frame(outlink, out);
|
|
}
|
|
|
|
#define OFFSET(x) offsetof(TonemapContext, x)
|
|
#define FLAGS AV_OPT_FLAG_VIDEO_PARAM | AV_OPT_FLAG_FILTERING_PARAM
|
|
static const AVOption tonemap_options[] = {
|
|
{ "tonemap", "tonemap algorithm selection", OFFSET(tonemap), AV_OPT_TYPE_INT, {.i64 = TONEMAP_NONE}, TONEMAP_NONE, TONEMAP_MAX - 1, FLAGS, "tonemap" },
|
|
{ "none", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_NONE}, 0, 0, FLAGS, "tonemap" },
|
|
{ "linear", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_LINEAR}, 0, 0, FLAGS, "tonemap" },
|
|
{ "gamma", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_GAMMA}, 0, 0, FLAGS, "tonemap" },
|
|
{ "clip", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_CLIP}, 0, 0, FLAGS, "tonemap" },
|
|
{ "reinhard", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_REINHARD}, 0, 0, FLAGS, "tonemap" },
|
|
{ "hable", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_HABLE}, 0, 0, FLAGS, "tonemap" },
|
|
{ "mobius", 0, 0, AV_OPT_TYPE_CONST, {.i64 = TONEMAP_MOBIUS}, 0, 0, FLAGS, "tonemap" },
|
|
{ "param", "tonemap parameter", OFFSET(param), AV_OPT_TYPE_DOUBLE, {.dbl = NAN}, DBL_MIN, DBL_MAX, FLAGS },
|
|
{ "desat", "desaturation strength", OFFSET(desat), AV_OPT_TYPE_DOUBLE, {.dbl = 2}, 0, DBL_MAX, FLAGS },
|
|
{ "peak", "signal peak override", OFFSET(peak), AV_OPT_TYPE_DOUBLE, {.dbl = 0}, 0, DBL_MAX, FLAGS },
|
|
{ NULL }
|
|
};
|
|
|
|
AVFILTER_DEFINE_CLASS(tonemap);
|
|
|
|
static const AVFilterPad tonemap_inputs[] = {
|
|
{
|
|
.name = "default",
|
|
.type = AVMEDIA_TYPE_VIDEO,
|
|
.filter_frame = filter_frame,
|
|
},
|
|
{ NULL }
|
|
};
|
|
|
|
static const AVFilterPad tonemap_outputs[] = {
|
|
{
|
|
.name = "default",
|
|
.type = AVMEDIA_TYPE_VIDEO,
|
|
},
|
|
{ NULL }
|
|
};
|
|
|
|
AVFilter ff_vf_tonemap = {
|
|
.name = "tonemap",
|
|
.description = NULL_IF_CONFIG_SMALL("Conversion to/from different dynamic ranges."),
|
|
.init = init,
|
|
.query_formats = query_formats,
|
|
.priv_size = sizeof(TonemapContext),
|
|
.priv_class = &tonemap_class,
|
|
.inputs = tonemap_inputs,
|
|
.outputs = tonemap_outputs,
|
|
.flags = AVFILTER_FLAG_SLICE_THREADS,
|
|
};
|