diff options
Diffstat (limited to 'libavfilter/vf_vectorscope.c')
-rw-r--r-- | libavfilter/vf_vectorscope.c | 1359 |
1 files changed, 1359 insertions, 0 deletions
diff --git a/libavfilter/vf_vectorscope.c b/libavfilter/vf_vectorscope.c new file mode 100644 index 0000000000..987bc66bd4 --- /dev/null +++ b/libavfilter/vf_vectorscope.c @@ -0,0 +1,1359 @@ +/* + * Copyright (c) 2015 Paul B Mahol + * + * 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 + */ + +#include "libavutil/avassert.h" +#include "libavutil/intreadwrite.h" +#include "libavutil/opt.h" +#include "libavutil/parseutils.h" +#include "libavutil/pixdesc.h" +#include "libavutil/xga_font_data.h" +#include "avfilter.h" +#include "formats.h" +#include "internal.h" +#include "video.h" + +enum VectorscopeMode { + GRAY, + COLOR, + COLOR2, + COLOR3, + COLOR4, + COLOR5, + MODE_NB +}; + +typedef struct VectorscopeContext { + const AVClass *class; + int mode; + int intensity; + float fintensity; + uint16_t bg_color[4]; + int planewidth[4]; + int planeheight[4]; + int hsub, vsub; + int x, y, pd; + int is_yuv; + int size; + int depth; + int mult; + int envelope; + int graticule; + float opacity; + float bgopacity; + float lthreshold; + float hthreshold; + int tmin; + int tmax; + int flags; + int colorspace; + int cs; + uint8_t *peak_memory; + uint8_t **peak; + + void (*vectorscope)(struct VectorscopeContext *s, + AVFrame *in, AVFrame *out, int pd); + void (*graticulef)(struct VectorscopeContext *s, AVFrame *out, + int X, int Y, int D, int P); +} VectorscopeContext; + +#define OFFSET(x) offsetof(VectorscopeContext, x) +#define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM + +static const AVOption vectorscope_options[] = { + { "mode", "set vectorscope mode", OFFSET(mode), AV_OPT_TYPE_INT, {.i64=0}, 0, MODE_NB-1, FLAGS, "mode"}, + { "m", "set vectorscope mode", OFFSET(mode), AV_OPT_TYPE_INT, {.i64=0}, 0, MODE_NB-1, FLAGS, "mode"}, + { "gray", 0, 0, AV_OPT_TYPE_CONST, {.i64=GRAY}, 0, 0, FLAGS, "mode" }, + { "color", 0, 0, AV_OPT_TYPE_CONST, {.i64=COLOR}, 0, 0, FLAGS, "mode" }, + { "color2", 0, 0, AV_OPT_TYPE_CONST, {.i64=COLOR2}, 0, 0, FLAGS, "mode" }, + { "color3", 0, 0, AV_OPT_TYPE_CONST, {.i64=COLOR3}, 0, 0, FLAGS, "mode" }, + { "color4", 0, 0, AV_OPT_TYPE_CONST, {.i64=COLOR4}, 0, 0, FLAGS, "mode" }, + { "color5", 0, 0, AV_OPT_TYPE_CONST, {.i64=COLOR5}, 0, 0, FLAGS, "mode" }, + { "x", "set color component on X axis", OFFSET(x), AV_OPT_TYPE_INT, {.i64=1}, 0, 2, FLAGS}, + { "y", "set color component on Y axis", OFFSET(y), AV_OPT_TYPE_INT, {.i64=2}, 0, 2, FLAGS}, + { "intensity", "set intensity", OFFSET(fintensity), AV_OPT_TYPE_FLOAT, {.dbl=0.004}, 0, 1, FLAGS}, + { "i", "set intensity", OFFSET(fintensity), AV_OPT_TYPE_FLOAT, {.dbl=0.004}, 0, 1, FLAGS}, + { "envelope", "set envelope", OFFSET(envelope), AV_OPT_TYPE_INT, {.i64=0}, 0, 3, FLAGS, "envelope"}, + { "e", "set envelope", OFFSET(envelope), AV_OPT_TYPE_INT, {.i64=0}, 0, 3, FLAGS, "envelope"}, + { "none", 0, 0, AV_OPT_TYPE_CONST, {.i64=0}, 0, 0, FLAGS, "envelope" }, + { "instant", 0, 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, FLAGS, "envelope" }, + { "peak", 0, 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, FLAGS, "envelope" }, + { "peak+instant", 0, 0, AV_OPT_TYPE_CONST, {.i64=3}, 0, 0, FLAGS, "envelope" }, + { "graticule", "set graticule", OFFSET(graticule), AV_OPT_TYPE_INT, {.i64=0}, 0, 2, FLAGS, "graticule"}, + { "g", "set graticule", OFFSET(graticule), AV_OPT_TYPE_INT, {.i64=0}, 0, 2, FLAGS, "graticule"}, + { "none", 0, 0, AV_OPT_TYPE_CONST, {.i64=0}, 0, 0, FLAGS, "graticule" }, + { "green", 0, 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, FLAGS, "graticule" }, + { "color", 0, 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, FLAGS, "graticule" }, + { "opacity", "set graticule opacity", OFFSET(opacity), AV_OPT_TYPE_FLOAT, {.dbl=0.75}, 0, 1, FLAGS}, + { "o", "set graticule opacity", OFFSET(opacity), AV_OPT_TYPE_FLOAT, {.dbl=0.75}, 0, 1, FLAGS}, + { "flags", "set graticule flags", OFFSET(flags), AV_OPT_TYPE_FLAGS, {.i64=4}, 0, 7, FLAGS, "flags"}, + { "f", "set graticule flags", OFFSET(flags), AV_OPT_TYPE_FLAGS, {.i64=4}, 0, 7, FLAGS, "flags"}, + { "white", "draw white point", 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, FLAGS, "flags" }, + { "black", "draw black point", 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, FLAGS, "flags" }, + { "name", "draw point name", 0, AV_OPT_TYPE_CONST, {.i64=4}, 0, 0, FLAGS, "flags" }, + { "bgopacity", "set background opacity", OFFSET(bgopacity), AV_OPT_TYPE_FLOAT, {.dbl=0.3}, 0, 1, FLAGS}, + { "b", "set background opacity", OFFSET(bgopacity), AV_OPT_TYPE_FLOAT, {.dbl=0.3}, 0, 1, FLAGS}, + { "lthreshold", "set low threshold", OFFSET(lthreshold), AV_OPT_TYPE_FLOAT, {.dbl=0}, 0, 1, FLAGS}, + { "l", "set low threshold", OFFSET(lthreshold), AV_OPT_TYPE_FLOAT, {.dbl=0}, 0, 1, FLAGS}, + { "hthreshold", "set high threshold", OFFSET(hthreshold), AV_OPT_TYPE_FLOAT, {.dbl=1}, 0, 1, FLAGS}, + { "h", "set high threshold", OFFSET(hthreshold), AV_OPT_TYPE_FLOAT, {.dbl=1}, 0, 1, FLAGS}, + { "colorspace", "set colorspace", OFFSET(colorspace), AV_OPT_TYPE_INT, {.i64=0}, 0, 2, FLAGS, "colorspace"}, + { "c", "set colorspace", OFFSET(colorspace), AV_OPT_TYPE_INT, {.i64=0}, 0, 2, FLAGS, "colorspace"}, + { "auto", 0, 0, AV_OPT_TYPE_CONST, {.i64=0}, 0, 0, FLAGS, "colorspace" }, + { "601", 0, 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, FLAGS, "colorspace" }, + { "709", 0, 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, FLAGS, "colorspace" }, + { NULL } +}; + +AVFILTER_DEFINE_CLASS(vectorscope); + +static const enum AVPixelFormat out_yuv8_pix_fmts[] = { + AV_PIX_FMT_YUVA444P, AV_PIX_FMT_YUV444P, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_yuv9_pix_fmts[] = { + AV_PIX_FMT_YUVA444P9, AV_PIX_FMT_YUV444P9, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_yuv10_pix_fmts[] = { + AV_PIX_FMT_YUVA444P10, AV_PIX_FMT_YUV444P10, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_yuv12_pix_fmts[] = { + AV_PIX_FMT_YUV444P12, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_rgb8_pix_fmts[] = { + AV_PIX_FMT_GBRAP, AV_PIX_FMT_GBRP, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_rgb9_pix_fmts[] = { + AV_PIX_FMT_GBRP9, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_rgb10_pix_fmts[] = { + AV_PIX_FMT_GBRP10, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat out_rgb12_pix_fmts[] = { + AV_PIX_FMT_GBRP12, AV_PIX_FMT_GBRAP12, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat in1_pix_fmts[] = { + AV_PIX_FMT_YUVA444P, AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUVJ444P, + AV_PIX_FMT_YUV444P9, AV_PIX_FMT_YUV444P10, + AV_PIX_FMT_YUVA444P9, AV_PIX_FMT_YUVA444P10, + AV_PIX_FMT_YUV444P12, + AV_PIX_FMT_GBRAP, AV_PIX_FMT_GBRP, + AV_PIX_FMT_GBRP9, AV_PIX_FMT_GBRP10, + AV_PIX_FMT_GBRP12, AV_PIX_FMT_GBRAP12, + AV_PIX_FMT_NONE +}; + +static const enum AVPixelFormat in2_pix_fmts[] = { + AV_PIX_FMT_YUVA420P, AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUVJ420P, + AV_PIX_FMT_YUVA422P, AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUVJ422P, + AV_PIX_FMT_YUVA444P, AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUVJ444P, + AV_PIX_FMT_YUV411P, AV_PIX_FMT_YUVJ411P, + AV_PIX_FMT_YUV440P, AV_PIX_FMT_YUV410P, + AV_PIX_FMT_GBRAP, AV_PIX_FMT_GBRP, + AV_PIX_FMT_GBRP9, AV_PIX_FMT_GBRP10, + AV_PIX_FMT_GBRP12, AV_PIX_FMT_GBRAP12, + AV_PIX_FMT_YUV420P9, AV_PIX_FMT_YUV422P9, AV_PIX_FMT_YUV444P9, + AV_PIX_FMT_YUV420P10, AV_PIX_FMT_YUV422P10, AV_PIX_FMT_YUV444P10, + AV_PIX_FMT_YUVA420P9, AV_PIX_FMT_YUVA422P9, AV_PIX_FMT_YUVA444P9, + AV_PIX_FMT_YUVA420P10, AV_PIX_FMT_YUVA422P10, AV_PIX_FMT_YUVA444P10, + AV_PIX_FMT_YUV420P12, AV_PIX_FMT_YUV422P12, AV_PIX_FMT_YUV444P12, AV_PIX_FMT_YUV440P12, + AV_PIX_FMT_NONE +}; + +static int query_formats(AVFilterContext *ctx) +{ + VectorscopeContext *s = ctx->priv; + const enum AVPixelFormat *out_pix_fmts; + const AVPixFmtDescriptor *desc; + AVFilterFormats *avff; + int depth, rgb, i, ret; + + if (!ctx->inputs[0]->in_formats || + !ctx->inputs[0]->in_formats->nb_formats) { + return AVERROR(EAGAIN); + } + + if (!ctx->inputs[0]->out_formats) { + const enum AVPixelFormat *in_pix_fmts; + + if ((s->x == 1 && s->y == 2) || (s->x == 2 && s->y == 1)) + in_pix_fmts = in2_pix_fmts; + else + in_pix_fmts = in1_pix_fmts; + if ((ret = ff_formats_ref(ff_make_format_list(in_pix_fmts), &ctx->inputs[0]->out_formats)) < 0) + return ret; + } + + avff = ctx->inputs[0]->in_formats; + desc = av_pix_fmt_desc_get(avff->formats[0]); + rgb = desc->flags & AV_PIX_FMT_FLAG_RGB; + depth = desc->comp[0].depth; + for (i = 1; i < avff->nb_formats; i++) { + desc = av_pix_fmt_desc_get(avff->formats[i]); + if (rgb != (desc->flags & AV_PIX_FMT_FLAG_RGB) || + depth != desc->comp[0].depth) + return AVERROR(EAGAIN); + } + + if (rgb && depth == 8) + out_pix_fmts = out_rgb8_pix_fmts; + else if (rgb && depth == 9) + out_pix_fmts = out_rgb9_pix_fmts; + else if (rgb && depth == 10) + out_pix_fmts = out_rgb10_pix_fmts; + else if (rgb && depth == 12) + out_pix_fmts = out_rgb12_pix_fmts; + else if (depth == 8) + out_pix_fmts = out_yuv8_pix_fmts; + else if (depth == 9) + out_pix_fmts = out_yuv9_pix_fmts; + else if (depth == 10) + out_pix_fmts = out_yuv10_pix_fmts; + else if (depth == 12) + out_pix_fmts = out_yuv12_pix_fmts; + else + return AVERROR(EAGAIN); + if ((ret = ff_formats_ref(ff_make_format_list(out_pix_fmts), &ctx->outputs[0]->in_formats)) < 0) + return ret; + + return 0; +} + +static int config_output(AVFilterLink *outlink) +{ + VectorscopeContext *s = outlink->src->priv; + int i; + + s->intensity = s->fintensity * (s->size - 1); + outlink->h = outlink->w = s->size; + outlink->sample_aspect_ratio = (AVRational){1,1}; + + s->peak_memory = av_calloc(s->size, s->size); + if (!s->peak_memory) + return AVERROR(ENOMEM); + + s->peak = av_calloc(s->size, sizeof(*s->peak)); + if (!s->peak) + return AVERROR(ENOMEM); + + for (i = 0; i < s->size; i++) + s->peak[i] = s->peak_memory + s->size * i; + + return 0; +} + +static void envelope_instant16(VectorscopeContext *s, AVFrame *out) +{ + const int dlinesize = out->linesize[0] / 2; + uint16_t *dpd = s->mode == COLOR || !s->is_yuv ? (uint16_t *)out->data[s->pd] : (uint16_t *)out->data[0]; + const int max = s->size - 1; + int i, j; + + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + const int pos = i * dlinesize + j; + const int poa = (i - 1) * dlinesize + j; + const int pob = (i + 1) * dlinesize + j; + + if (dpd[pos] && (((!j || !dpd[pos - 1]) || ((j == (out->width - 1)) || !dpd[pos + 1])) + || ((!i || !dpd[poa]) || ((i == (out->height - 1)) || !dpd[pob])))) { + dpd[pos] = max; + } + } + } +} + +static void envelope_peak16(VectorscopeContext *s, AVFrame *out) +{ + const int dlinesize = out->linesize[0] / 2; + uint16_t *dpd = s->mode == COLOR || !s->is_yuv ? (uint16_t *)out->data[s->pd] : (uint16_t *)out->data[0]; + const int max = s->size - 1; + int i, j; + + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + const int pos = i * dlinesize + j; + + if (dpd[pos]) + s->peak[i][j] = 1; + } + } + + if (s->envelope == 3) + envelope_instant16(s, out); + + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + const int pos = i * dlinesize + j; + + if (s->peak[i][j] && (((!j || !s->peak[i][j-1]) || ((j == (out->width - 1)) || !s->peak[i][j + 1])) + || ((!i || !s->peak[i-1][j]) || ((i == (out->height - 1)) || !s->peak[i + 1][j])))) { + dpd[pos] = max; + } + } + } +} + +static void envelope_instant(VectorscopeContext *s, AVFrame *out) +{ + const int dlinesize = out->linesize[0]; + uint8_t *dpd = s->mode == COLOR || !s->is_yuv ? out->data[s->pd] : out->data[0]; + int i, j; + + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + const int pos = i * dlinesize + j; + const int poa = (i - 1) * dlinesize + j; + const int pob = (i + 1) * dlinesize + j; + + if (dpd[pos] && (((!j || !dpd[pos - 1]) || ((j == (out->width - 1)) || !dpd[pos + 1])) + || ((!i || !dpd[poa]) || ((i == (out->height - 1)) || !dpd[pob])))) { + dpd[pos] = 255; + } + } + } +} + +static void envelope_peak(VectorscopeContext *s, AVFrame *out) +{ + const int dlinesize = out->linesize[0]; + uint8_t *dpd = s->mode == COLOR || !s->is_yuv ? out->data[s->pd] : out->data[0]; + int i, j; + + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + const int pos = i * dlinesize + j; + + if (dpd[pos]) + s->peak[i][j] = 1; + } + } + + if (s->envelope == 3) + envelope_instant(s, out); + + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + const int pos = i * dlinesize + j; + + if (s->peak[i][j] && (((!j || !s->peak[i][j-1]) || ((j == (out->width - 1)) || !s->peak[i][j + 1])) + || ((!i || !s->peak[i-1][j]) || ((i == (out->height - 1)) || !s->peak[i + 1][j])))) { + dpd[pos] = 255; + } + } + } +} + +static void envelope16(VectorscopeContext *s, AVFrame *out) +{ + if (!s->envelope) { + return; + } else if (s->envelope == 1) { + envelope_instant16(s, out); + } else { + envelope_peak16(s, out); + } +} + +static void envelope(VectorscopeContext *s, AVFrame *out) +{ + if (!s->envelope) { + return; + } else if (s->envelope == 1) { + envelope_instant(s, out); + } else { + envelope_peak(s, out); + } +} + +static void vectorscope16(VectorscopeContext *s, AVFrame *in, AVFrame *out, int pd) +{ + const uint16_t * const *src = (const uint16_t * const *)in->data; + const int slinesizex = in->linesize[s->x] / 2; + const int slinesizey = in->linesize[s->y] / 2; + const int slinesized = in->linesize[pd] / 2; + const int dlinesize = out->linesize[0] / 2; + const int intensity = s->intensity; + const int px = s->x, py = s->y; + const int h = s->planeheight[py]; + const int w = s->planewidth[px]; + const uint16_t *spx = src[px]; + const uint16_t *spy = src[py]; + const uint16_t *spd = src[pd]; + const int hsub = s->hsub; + const int vsub = s->vsub; + uint16_t **dst = (uint16_t **)out->data; + uint16_t *dpx = dst[px]; + uint16_t *dpy = dst[py]; + uint16_t *dpd = dst[pd]; + const int max = s->size - 1; + const int mid = s->size / 2; + const int tmin = s->tmin; + const int tmax = s->tmax; + int i, j, k; + + for (k = 0; k < 4 && dst[k]; k++) { + for (i = 0; i < out->height ; i++) + for (j = 0; j < out->width; j++) + AV_WN16(out->data[k] + i * out->linesize[k] + j * 2, + (s->mode == COLOR || s->mode == COLOR5) && k == s->pd ? 0 : s->bg_color[k]); + } + + switch (s->mode) { + case COLOR: + case COLOR5: + case GRAY: + if (s->is_yuv) { + for (i = 0; i < h; i++) { + const int iwx = i * slinesizex; + const int iwy = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = FFMIN(spx[iwx + j], max); + const int y = FFMIN(spy[iwy + j], max); + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dpd[pos] = FFMIN(dpd[pos] + intensity, max); + } + } + } else { + for (i = 0; i < h; i++) { + const int iwx = i * slinesizex; + const int iwy = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = FFMIN(spx[iwx + j], max); + const int y = FFMIN(spy[iwy + j], max); + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dst[0][pos] = FFMIN(dst[0][pos] + intensity, max); + dst[1][pos] = FFMIN(dst[1][pos] + intensity, max); + dst[2][pos] = FFMIN(dst[2][pos] + intensity, max); + } + } + } + break; + case COLOR2: + if (s->is_yuv) { + for (i = 0; i < h; i++) { + const int iw1 = i * slinesizex; + const int iw2 = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = FFMIN(spx[iw1 + j], max); + const int y = FFMIN(spy[iw2 + j], max); + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + if (!dpd[pos]) + dpd[pos] = FFABS(mid - x) + FFABS(mid - y); + dpx[pos] = x; + dpy[pos] = y; + } + } + } else { + for (i = 0; i < h; i++) { + const int iw1 = i * slinesizex; + const int iw2 = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = FFMIN(spx[iw1 + j], max); + const int y = FFMIN(spy[iw2 + j], max); + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + if (!dpd[pos]) + dpd[pos] = FFMIN(x + y, max); + dpx[pos] = x; + dpy[pos] = y; + } + } + } + break; + case COLOR3: + for (i = 0; i < h; i++) { + const int iw1 = i * slinesizex; + const int iw2 = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = FFMIN(spx[iw1 + j], max); + const int y = FFMIN(spy[iw2 + j], max); + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dpd[pos] = FFMIN(max, dpd[pos] + intensity); + dpx[pos] = x; + dpy[pos] = y; + } + } + break; + case COLOR4: + for (i = 0; i < in->height; i++) { + const int iwx = (i >> vsub) * slinesizex; + const int iwy = (i >> vsub) * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < in->width; j++) { + const int x = FFMIN(spx[iwx + (j >> hsub)], max); + const int y = FFMIN(spy[iwy + (j >> hsub)], max); + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dpd[pos] = FFMAX(z, dpd[pos]); + dpx[pos] = x; + dpy[pos] = y; + } + } + break; + default: + av_assert0(0); + } + + envelope16(s, out); + + if (dst[3]) { + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + int pos = i * dlinesize + j; + + if (dpd[pos]) + dst[3][pos] = max; + } + } + } + + if (s->mode == COLOR) { + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + if (!dpd[i * dlinesize + j]) { + dpx[i * dlinesize + j] = j; + dpy[i * dlinesize + j] = i; + dpd[i * dlinesize + j] = mid; + } + } + } + } else if (s->mode == COLOR5) { + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + if (!dpd[i * dlinesize + j]) { + dpx[i * dlinesize + j] = j; + dpy[i * dlinesize + j] = i; + dpd[i * dlinesize + j] = mid * M_SQRT2 - hypot(i - mid, j - mid); + } + } + } + } +} + +static void vectorscope8(VectorscopeContext *s, AVFrame *in, AVFrame *out, int pd) +{ + const uint8_t * const *src = (const uint8_t * const *)in->data; + const int slinesizex = in->linesize[s->x]; + const int slinesizey = in->linesize[s->y]; + const int slinesized = in->linesize[pd]; + const int dlinesize = out->linesize[0]; + const int intensity = s->intensity; + const int px = s->x, py = s->y; + const int h = s->planeheight[py]; + const int w = s->planewidth[px]; + const uint8_t *spx = src[px]; + const uint8_t *spy = src[py]; + const uint8_t *spd = src[pd]; + const int hsub = s->hsub; + const int vsub = s->vsub; + uint8_t **dst = out->data; + uint8_t *dpx = dst[px]; + uint8_t *dpy = dst[py]; + uint8_t *dpd = dst[pd]; + const int tmin = s->tmin; + const int tmax = s->tmax; + int i, j, k; + + for (k = 0; k < 4 && dst[k]; k++) + for (i = 0; i < out->height ; i++) + memset(dst[k] + i * out->linesize[k], + (s->mode == COLOR || s->mode == COLOR5) && k == s->pd ? 0 : s->bg_color[k], out->width); + + switch (s->mode) { + case COLOR5: + case COLOR: + case GRAY: + if (s->is_yuv) { + for (i = 0; i < h; i++) { + const int iwx = i * slinesizex; + const int iwy = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = spx[iwx + j]; + const int y = spy[iwy + j]; + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dpd[pos] = FFMIN(dpd[pos] + intensity, 255); + } + } + } else { + for (i = 0; i < h; i++) { + const int iwx = i * slinesizex; + const int iwy = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = spx[iwx + j]; + const int y = spy[iwy + j]; + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dst[0][pos] = FFMIN(dst[0][pos] + intensity, 255); + dst[1][pos] = FFMIN(dst[1][pos] + intensity, 255); + dst[2][pos] = FFMIN(dst[2][pos] + intensity, 255); + } + } + } + break; + case COLOR2: + if (s->is_yuv) { + for (i = 0; i < h; i++) { + const int iw1 = i * slinesizex; + const int iw2 = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = spx[iw1 + j]; + const int y = spy[iw2 + j]; + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + if (!dpd[pos]) + dpd[pos] = FFABS(128 - x) + FFABS(128 - y); + dpx[pos] = x; + dpy[pos] = y; + } + } + } else { + for (i = 0; i < h; i++) { + const int iw1 = i * slinesizex; + const int iw2 = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = spx[iw1 + j]; + const int y = spy[iw2 + j]; + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + if (!dpd[pos]) + dpd[pos] = FFMIN(x + y, 255); + dpx[pos] = x; + dpy[pos] = y; + } + } + } + break; + case COLOR3: + for (i = 0; i < h; i++) { + const int iw1 = i * slinesizex; + const int iw2 = i * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < w; j++) { + const int x = spx[iw1 + j]; + const int y = spy[iw2 + j]; + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dpd[pos] = FFMIN(255, dpd[pos] + intensity); + dpx[pos] = x; + dpy[pos] = y; + } + } + break; + case COLOR4: + for (i = 0; i < in->height; i++) { + const int iwx = (i >> vsub) * slinesizex; + const int iwy = (i >> vsub) * slinesizey; + const int iwd = i * slinesized; + for (j = 0; j < in->width; j++) { + const int x = spx[iwx + (j >> hsub)]; + const int y = spy[iwy + (j >> hsub)]; + const int z = spd[iwd + j]; + const int pos = y * dlinesize + x; + + if (z < tmin || z > tmax) + continue; + + dpd[pos] = FFMAX(z, dpd[pos]); + dpx[pos] = x; + dpy[pos] = y; + } + } + break; + default: + av_assert0(0); + } + + envelope(s, out); + + if (dst[3]) { + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + int pos = i * dlinesize + j; + + if (dpd[pos]) + dst[3][pos] = 255; + } + } + } + + if (s->mode == COLOR) { + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + if (!dpd[i * out->linesize[pd] + j]) { + dpx[i * out->linesize[px] + j] = j; + dpy[i * out->linesize[py] + j] = i; + dpd[i * out->linesize[pd] + j] = 128; + } + } + } + } else if (s->mode == COLOR5) { + for (i = 0; i < out->height; i++) { + for (j = 0; j < out->width; j++) { + if (!dpd[i * out->linesize[pd] + j]) { + dpx[i * out->linesize[px] + j] = j; + dpy[i * out->linesize[py] + j] = i; + dpd[i * out->linesize[pd] + j] = 128 * M_SQRT2 - hypot(i - 128, j - 128); + } + } + } + } +} + +const static char *positions_name[] = { + "R", "B", "Cy", "Yl", "G", "Mg", +}; + +const static uint16_t positions[][14][3] = { + { + { 81, 90, 240 }, { 41, 240, 110 }, { 170, 166, 16 }, + { 210, 16, 146 }, { 145, 54, 34 }, { 106, 202, 222 }, + { 162, 44, 142 }, { 131, 156, 44 }, { 112, 72, 58 }, + { 84, 184, 198 }, { 65, 100, 212 }, { 35, 212, 114 }, + { 235, 128, 128 }, { 16, 128, 128 } }, + { { 63, 102, 240 }, { 32, 240, 118 }, { 188, 154, 16 }, + { 219, 16, 138 }, { 173, 42, 26 }, { 78, 214, 230 }, + { 28, 212, 120 }, { 51, 109, 212 }, { 63, 193, 204 }, + { 133, 63, 52 }, { 145, 147, 44 }, { 168, 44, 136 }, + { 235, 128, 128 }, { 16, 128, 128 } }, + { { 81*2, 90*2, 240*2 }, { 41*2, 240*2, 110*2 }, { 170*2, 166*2, 16*2 }, + { 210*2, 16*2, 146*2 }, { 145*2, 54*2, 34*2 }, { 106*2, 202*2, 222*2 }, + { 162*2, 44*2, 142*2 }, { 131*2, 156*2, 44*2 }, { 112*2, 72*2, 58*2 }, + { 84*2, 184*2, 198*2 }, { 65*2, 100*2, 212*2 }, { 35*2, 212*2, 114*2 }, + { 470, 256, 256 }, { 32, 256, 256 } }, + { { 63*2, 102*2, 240*2 }, { 32*2, 240*2, 118*2 }, { 188*2, 154*2, 16*2 }, + { 219*2, 16*2, 138*2 }, { 173*2, 42*2, 26*2 }, { 78*2, 214*2, 230*2 }, + { 28*2, 212*2, 120*2 }, { 51*2, 109*2, 212*2 }, { 63*2, 193*2, 204*2 }, + { 133*2, 63*2, 52*2 }, { 145*2, 147*2, 44*2 }, { 168*2, 44*2, 136*2 }, + { 470, 256, 256 }, { 32, 256, 256 } }, + { { 81*4, 90*4, 240*4 }, { 41*4, 240*4, 110*4 }, { 170*4, 166*4, 16*4 }, + { 210*4, 16*4, 146*4 }, { 145*4, 54*4, 34*4 }, { 106*4, 202*4, 222*4 }, + { 162*4, 44*4, 142*4 }, { 131*4, 156*4, 44*4 }, { 112*4, 72*4, 58*4 }, + { 84*4, 184*4, 198*4 }, { 65*4, 100*4, 212*4 }, { 35*4, 212*4, 114*4 }, + { 940, 512, 512 }, { 64, 512, 512 } }, + { { 63*4, 102*4, 240*4 }, { 32*4, 240*4, 118*4 }, { 188*4, 154*4, 16*4 }, + { 219*4, 16*4, 138*4 }, { 173*4, 42*4, 26*4 }, { 78*4, 214*4, 230*4 }, + { 28*4, 212*4, 120*4 }, { 51*4, 109*4, 212*4 }, { 63*4, 193*4, 204*4 }, + { 133*4, 63*4, 52*4 }, { 145*4, 147*4, 44*4 }, { 168*4, 44*4, 136*4 }, + { 940, 512, 512 }, { 64, 512, 512 } }, + { { 81*8, 90*4, 240*8 }, { 41*8, 240*8, 110*8 }, { 170*8, 166*8, 16*8 }, + { 210*8, 16*4, 146*8 }, { 145*8, 54*8, 34*8 }, { 106*8, 202*8, 222*8 }, + { 162*8, 44*4, 142*8 }, { 131*8, 156*8, 44*8 }, { 112*8, 72*8, 58*8 }, + { 84*8, 184*4, 198*8 }, { 65*8, 100*8, 212*8 }, { 35*8, 212*8, 114*8 }, + { 1880, 1024, 1024 }, { 128, 1024, 1024 } }, + { { 63*8, 102*8, 240*8 }, { 32*8, 240*8, 118*8 }, { 188*8, 154*8, 16*8 }, + { 219*8, 16*8, 138*8 }, { 173*8, 42*8, 26*8 }, { 78*8, 214*8, 230*8 }, + { 28*8, 212*8, 120*8 }, { 51*8, 109*8, 212*8 }, { 63*8, 193*8, 204*8 }, + { 133*8, 63*8, 52*8 }, { 145*8, 147*8, 44*8 }, { 168*8, 44*8, 136*8 }, + { 1880, 1024, 1024 }, { 128, 1024, 1024 } }, + { { 81*16, 90*16, 240*16 }, { 41*16, 240*16, 110*16 }, { 170*16, 166*16, 16*16 }, + { 210*16, 16*16, 146*16 }, { 145*16, 54*16, 34*16 }, { 106*16, 202*16, 222*16 }, + { 162*16, 44*16, 142*16 }, { 131*16, 156*16, 44*16 }, { 112*16, 72*16, 58*16 }, + { 84*16, 184*16, 198*16 }, { 65*16, 100*16, 212*16 }, { 35*16, 212*16, 114*16 }, + { 3760, 2048, 2048 }, { 256, 2048, 2048 } }, + { { 63*16, 102*16, 240*16 }, { 32*16, 240*16, 118*16 }, { 188*16, 154*16, 16*16 }, + { 219*16, 16*16, 138*16 }, { 173*16, 42*16, 26*16 }, { 78*16, 214*16, 230*16 }, + { 28*16, 212*16, 120*16 }, { 51*16, 109*16, 212*16 }, { 63*16, 193*16, 204*16 }, + { 133*16, 63*16, 52*16 }, { 145*16, 147*16, 44*16 }, { 168*16, 44*16, 136*16 }, + { 3760, 2048, 2048 }, { 256, 2048, 2048 } }, +}; + +static void draw_dots(uint8_t *dst, int L, int v, float o) +{ + const float f = 1. - o; + const float V = o * v; + int l = L * 2; + + dst[ l - 3] = dst[ l - 3] * f + V; + dst[ l + 3] = dst[ l + 3] * f + V; + dst[-l - 3] = dst[-l - 3] * f + V; + dst[-l + 3] = dst[-l + 3] * f + V; + + l += L; + + dst[ l - 3] = dst[ l - 3] * f + V; + dst[ l + 3] = dst[ l + 3] * f + V; + dst[ l - 2] = dst[ l - 2] * f + V; + dst[ l + 2] = dst[ l + 2] * f + V; + dst[-l - 3] = dst[-l - 3] * f + V; + dst[-l + 3] = dst[-l + 3] * f + V; + dst[-l - 2] = dst[-l - 2] * f + V; + dst[-l + 2] = dst[-l + 2] * f + V; +} + +static void draw_dots16(uint16_t *dst, int L, int v, float o) +{ + const float f = 1. - o; + const float V = o * v; + int l = L * 2; + + dst[ l - 3] = dst[ l - 3] * f + V; + dst[ l + 3] = dst[ l + 3] * f + V; + dst[-l - 3] = dst[-l - 3] * f + V; + dst[-l + 3] = dst[-l + 3] * f + V; + + l += L; + + dst[ l - 3] = dst[ l - 3] * f + V; + dst[ l + 3] = dst[ l + 3] * f + V; + dst[ l - 2] = dst[ l - 2] * f + V; + dst[ l + 2] = dst[ l + 2] * f + V; + dst[-l - 3] = dst[-l - 3] * f + V; + dst[-l + 3] = dst[-l + 3] * f + V; + dst[-l - 2] = dst[-l - 2] * f + V; + dst[-l + 2] = dst[-l + 2] * f + V; +} + +static void none_graticule(VectorscopeContext *s, AVFrame *out, int X, int Y, int D, int P) +{ +} + +static void draw_htext(AVFrame *out, int x, int y, float o1, float o2, const char *txt, const uint8_t color[4]) +{ + const uint8_t *font; + int font_height; + int i, plane; + + font = avpriv_cga_font, font_height = 8; + + for (plane = 0; plane < 4 && out->data[plane]; plane++) { + for (i = 0; txt[i]; i++) { + int char_y, mask; + int v = color[plane]; + + uint8_t *p = out->data[plane] + y * out->linesize[plane] + (x + i * 8); + for (char_y = font_height - 1; char_y >= 0; char_y--) { + for (mask = 0x80; mask; mask >>= 1) { + if (font[txt[i] * font_height + char_y] & mask) + p[0] = p[0] * o2 + v * o1; + p++; + } + p += out->linesize[plane] - 8; + } + } + } +} + +static void draw_htext16(AVFrame *out, int x, int y, float o1, float o2, const char *txt, const uint16_t color[4]) +{ + const uint8_t *font; + int font_height; + int i, plane; + + font = avpriv_cga_font, font_height = 8; + + for (plane = 0; plane < 4 && out->data[plane]; plane++) { + for (i = 0; txt[i]; i++) { + int char_y, mask; + int v = color[plane]; + + uint16_t *p = (uint16_t *)(out->data[plane] + y * out->linesize[plane]) + (x + i * 8); + for (char_y = font_height - 1; char_y >= 0; char_y--) { + for (mask = 0x80; mask; mask >>= 1) { + if (font[txt[i] * font_height + char_y] & mask) + p[0] = p[0] * o2 + v * o1; + p++; + } + p += out->linesize[plane] / 2 - 8; + } + } + } +} + +static void color_graticule16(VectorscopeContext *s, AVFrame *out, int X, int Y, int D, int P) +{ + const int max = s->size - 1; + const float o = s->opacity; + int i; + + for (i = 0; i < 12; i++) { + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + int d = positions[P][i][D]; + + draw_dots16((uint16_t *)(out->data[D] + y * out->linesize[D] + x * 2), out->linesize[D] / 2, d, o); + draw_dots16((uint16_t *)(out->data[X] + y * out->linesize[X] + x * 2), out->linesize[X] / 2, x, o); + draw_dots16((uint16_t *)(out->data[Y] + y * out->linesize[Y] + x * 2), out->linesize[Y] / 2, y, o); + if (out->data[3]) + draw_dots16((uint16_t *)(out->data[3] + y * out->linesize[3] + x * 2), out->linesize[3] / 2, max, o); + } + + if (s->flags & 1) { + int x = positions[P][12][X]; + int y = positions[P][12][Y]; + int d = positions[P][12][D]; + + draw_dots16((uint16_t *)(out->data[D] + y * out->linesize[D] + x * 2), out->linesize[D] / 2, d, o); + draw_dots16((uint16_t *)(out->data[X] + y * out->linesize[X] + x * 2), out->linesize[X] / 2, x, o); + draw_dots16((uint16_t *)(out->data[Y] + y * out->linesize[Y] + x * 2), out->linesize[Y] / 2, y, o); + if (out->data[3]) + draw_dots16((uint16_t *)(out->data[3] + y * out->linesize[3] + x * 2), out->linesize[3] / 2, max, o); + } + + if (s->flags & 2) { + int x = positions[P][13][X]; + int y = positions[P][13][Y]; + int d = positions[P][13][D]; + + draw_dots16((uint16_t *)(out->data[D] + y * out->linesize[D] + x * 2), out->linesize[D] / 2, d, o); + draw_dots16((uint16_t *)(out->data[X] + y * out->linesize[X] + x * 2), out->linesize[X] / 2, x, o); + draw_dots16((uint16_t *)(out->data[Y] + y * out->linesize[Y] + x * 2), out->linesize[Y] / 2, y, o); + if (out->data[3]) + draw_dots16((uint16_t *)(out->data[3] + y * out->linesize[3] + x * 2), out->linesize[3] / 2, max, o); + } + + for (i = 0; i < 6 && s->flags & 4; i++) { + uint16_t color[4] = { 0, 0, 0, 0 }; + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + int d = positions[P][i][D]; + + color[D] = d; + color[X] = x; + color[Y] = y; + color[3] = max; + + if (x > max / 2) + x += 8; + else + x -= 14; + if (y > max / 2) + y += 8; + else + y -= 14; + + x = av_clip(x, 0, out->width - 9); + y = av_clip(y, 0, out->height - 9); + draw_htext16(out, x, y, o, 1. - o, positions_name[i], color); + } +} + +static void color_graticule(VectorscopeContext *s, AVFrame *out, int X, int Y, int D, int P) +{ + const float o = s->opacity; + int i; + + for (i = 0; i < 12; i++) { + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + int d = positions[P][i][D]; + + draw_dots(out->data[D] + y * out->linesize[D] + x, out->linesize[D], d, o); + draw_dots(out->data[X] + y * out->linesize[X] + x, out->linesize[X], x, o); + draw_dots(out->data[Y] + y * out->linesize[Y] + x, out->linesize[Y], y, o); + if (out->data[3]) + draw_dots(out->data[3] + y * out->linesize[3] + x, out->linesize[3], 255, o); + } + + if (s->flags & 1) { + int x = positions[P][12][X]; + int y = positions[P][12][Y]; + int d = positions[P][12][D]; + + draw_dots(out->data[D] + y * out->linesize[D] + x, out->linesize[D], d, o); + draw_dots(out->data[X] + y * out->linesize[X] + x, out->linesize[X], x, o); + draw_dots(out->data[Y] + y * out->linesize[Y] + x, out->linesize[Y], y, o); + if (out->data[3]) + draw_dots(out->data[3] + y * out->linesize[3] + x, out->linesize[3], 255, o); + } + + if (s->flags & 2) { + int x = positions[P][13][X]; + int y = positions[P][13][Y]; + int d = positions[P][12][D]; + + draw_dots(out->data[D] + y * out->linesize[D] + x, out->linesize[D], d, o); + draw_dots(out->data[X] + y * out->linesize[X] + x, out->linesize[X], x, o); + draw_dots(out->data[Y] + y * out->linesize[Y] + x, out->linesize[Y], y, o); + if (out->data[3]) + draw_dots(out->data[3] + y * out->linesize[3] + x, out->linesize[3], 255, o); + } + + for (i = 0; i < 6 && s->flags & 4; i++) { + uint8_t color[4] = { 0, 0, 0, 255 }; + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + int d = positions[P][i][D]; + + color[D] = d; + color[X] = x; + color[Y] = y; + + if (x > 128) + x += 8; + else + x -= 14; + if (y > 128) + y += 8; + else + y -= 14; + + x = av_clip(x, 0, out->width - 9); + y = av_clip(y, 0, out->height - 9); + draw_htext(out, x, y, o, 1. - o, positions_name[i], color); + } +} + +static void green_graticule16(VectorscopeContext *s, AVFrame *out, int X, int Y, int D, int P) +{ + const int max = s->size - 1; + const float o = s->opacity; + const int m = s->mult; + int i; + + for (i = 0; i < 12; i++) { + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + + draw_dots16((uint16_t *)(out->data[0] + y * out->linesize[0] + x * 2), out->linesize[0] / 2, 128 * m, o); + draw_dots16((uint16_t *)(out->data[1] + y * out->linesize[1] + x * 2), out->linesize[1] / 2, 0, o); + draw_dots16((uint16_t *)(out->data[2] + y * out->linesize[2] + x * 2), out->linesize[2] / 2, 0, o); + if (out->data[3]) + draw_dots16((uint16_t *)(out->data[3] + y * out->linesize[3] + x * 2), out->linesize[3] / 2, max, o); + } + + if (s->flags & 1) { + int x = positions[P][12][X]; + int y = positions[P][12][Y]; + + draw_dots16((uint16_t *)(out->data[0] + y * out->linesize[0] + x * 2), out->linesize[0] / 2, 128 * m, o); + draw_dots16((uint16_t *)(out->data[1] + y * out->linesize[1] + x * 2), out->linesize[1] / 2, 0, o); + draw_dots16((uint16_t *)(out->data[2] + y * out->linesize[2] + x * 2), out->linesize[2] / 2, 0, o); + if (out->data[3]) + draw_dots16((uint16_t *)(out->data[3] + y * out->linesize[3] + x * 2), out->linesize[3] / 2, max, o); + } + + if (s->flags & 2) { + int x = positions[P][13][X]; + int y = positions[P][13][Y]; + + draw_dots16((uint16_t *)(out->data[0] + y * out->linesize[0] + x * 2), out->linesize[0] / 2, 128 * m, o); + draw_dots16((uint16_t *)(out->data[1] + y * out->linesize[1] + x * 2), out->linesize[1] / 2, 0, o); + draw_dots16((uint16_t *)(out->data[2] + y * out->linesize[2] + x * 2), out->linesize[2] / 2, 0, o); + if (out->data[3]) + draw_dots16((uint16_t *)(out->data[3] + y * out->linesize[3] + x * 2), out->linesize[3] / 2, max, o); + } + + for (i = 0; i < 6 && s->flags & 4; i++) { + const uint16_t color[4] = { 128 * m, 0, 0, max }; + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + + if (x > max / 2) + x += 8; + else + x -= 14; + if (y > max / 2) + y += 8; + else + y -= 14; + + x = av_clip(x, 0, out->width - 9); + y = av_clip(y, 0, out->height - 9); + draw_htext16(out, x, y, o, 1. - o, positions_name[i], color); + } +} + +static void green_graticule(VectorscopeContext *s, AVFrame *out, int X, int Y, int D, int P) +{ + const float o = s->opacity; + int i; + + for (i = 0; i < 12; i++) { + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + + draw_dots(out->data[0] + y * out->linesize[0] + x, out->linesize[0], 128, o); + draw_dots(out->data[1] + y * out->linesize[1] + x, out->linesize[1], 0, o); + draw_dots(out->data[2] + y * out->linesize[2] + x, out->linesize[2], 0, o); + if (out->data[3]) + draw_dots(out->data[3] + y * out->linesize[3] + x, out->linesize[3], 255, o); + } + + if (s->flags & 1) { + int x = positions[P][12][X]; + int y = positions[P][12][Y]; + + draw_dots(out->data[0] + y * out->linesize[0] + x, out->linesize[0], 128, o); + draw_dots(out->data[1] + y * out->linesize[1] + x, out->linesize[1], 0, o); + draw_dots(out->data[2] + y * out->linesize[2] + x, out->linesize[2], 0, o); + if (out->data[3]) + draw_dots(out->data[3] + y * out->linesize[3] + x, out->linesize[3], 255, o); + } + + if (s->flags & 2) { + int x = positions[P][13][X]; + int y = positions[P][13][Y]; + + draw_dots(out->data[0] + y * out->linesize[0] + x, out->linesize[0], 128, o); + draw_dots(out->data[1] + y * out->linesize[1] + x, out->linesize[1], 0, o); + draw_dots(out->data[2] + y * out->linesize[2] + x, out->linesize[2], 0, o); + if (out->data[3]) + draw_dots(out->data[3] + y * out->linesize[3] + x, out->linesize[3], 255, o); + } + + for (i = 0; i < 6 && s->flags & 4; i++) { + const uint8_t color[4] = { 128, 0, 0, 255 }; + int x = positions[P][i][X]; + int y = positions[P][i][Y]; + + if (x > 128) + x += 8; + else + x -= 14; + if (y > 128) + y += 8; + else + y -= 14; + + x = av_clip(x, 0, out->width - 9); + y = av_clip(y, 0, out->height - 9); + draw_htext(out, x, y, o, 1. - o, positions_name[i], color); + } +} + +static int filter_frame(AVFilterLink *inlink, AVFrame *in) +{ + AVFilterContext *ctx = inlink->dst; + VectorscopeContext *s = ctx->priv; + AVFilterLink *outlink = ctx->outputs[0]; + AVFrame *out; + int plane; + + if (s->colorspace) { + s->cs = (s->depth - 8) * 2 + s->colorspace - 1; + } else { + switch (av_frame_get_colorspace(in)) { + case AVCOL_SPC_SMPTE170M: + case AVCOL_SPC_BT470BG: + s->cs = (s->depth - 8) * 2 + 0; + break; + case AVCOL_SPC_BT709: + default: + s->cs = (s->depth - 8) * 2 + 1; + } + } + + out = ff_get_video_buffer(outlink, outlink->w, outlink->h); + if (!out) { + av_frame_free(&in); + return AVERROR(ENOMEM); + } + av_frame_copy_props(out, in); + + s->vectorscope(s, in, out, s->pd); + s->graticulef(s, out, s->x, s->y, s->pd, s->cs); + + for (plane = 0; plane < 4; plane++) { + if (out->data[plane]) { + out->data[plane] += (s->size - 1) * out->linesize[plane]; + out->linesize[plane] = -out->linesize[plane]; + } + } + + av_frame_free(&in); + return ff_filter_frame(outlink, out); +} + +static int config_input(AVFilterLink *inlink) +{ + const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format); + AVFilterContext *ctx = inlink->dst; + VectorscopeContext *s = ctx->priv; + + s->is_yuv = !(desc->flags & AV_PIX_FMT_FLAG_RGB); + s->size = 1 << desc->comp[0].depth; + s->mult = s->size / 256; + s->depth = desc->comp[0].depth; + s->tmin = s->lthreshold * (s->size - 1); + s->tmax = s->hthreshold * (s->size - 1); + + if (s->tmin > s->tmax) { + av_log(ctx, AV_LOG_ERROR, "low threshold should be less than high threshold\n"); + return AVERROR(EINVAL); + } + + if (s->mode == GRAY && s->is_yuv) + s->pd = 0; + else { + if ((s->x == 1 && s->y == 2) || (s->x == 2 && s->y == 1)) + s->pd = 0; + else if ((s->x == 0 && s->y == 2) || (s->x == 2 && s->y == 0)) + s->pd = 1; + else if ((s->x == 0 && s->y == 1) || (s->x == 1 && s->y == 0)) + s->pd = 2; + } + + if (s->size == 256) + s->vectorscope = vectorscope8; + else + s->vectorscope = vectorscope16; + + s->graticulef = none_graticule; + + if (s->is_yuv && s->size == 256) { + if (s->graticule == 1) + s->graticulef = green_graticule; + else if (s->graticule == 2) + s->graticulef = color_graticule; + } else if (s->is_yuv) { + if (s->graticule == 1) + s->graticulef = green_graticule16; + else if (s->graticule == 2) + s->graticulef = color_graticule16; + } + + s->bg_color[3] = s->bgopacity * (s->size - 1); + + switch (inlink->format) { + case AV_PIX_FMT_GBRP12: + case AV_PIX_FMT_GBRP10: + case AV_PIX_FMT_GBRP9: + case AV_PIX_FMT_GBRAP: + case AV_PIX_FMT_GBRP: + s->bg_color[0] = 0; + s->bg_color[1] = 0; + s->bg_color[2] = 0; + break; + default: + s->bg_color[0] = 0; + s->bg_color[1] = s->size / 2 - 1; + s->bg_color[2] = s->size / 2 - 1; + } + + s->hsub = desc->log2_chroma_w; + s->vsub = desc->log2_chroma_h; + s->planeheight[1] = s->planeheight[2] = AV_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h); + s->planeheight[0] = s->planeheight[3] = inlink->h; + s->planewidth[1] = s->planewidth[2] = AV_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w); + s->planewidth[0] = s->planewidth[3] = inlink->w; + + return 0; +} + +static av_cold void uninit(AVFilterContext *ctx) +{ + VectorscopeContext *s = ctx->priv; + + av_freep(&s->peak); + av_freep(&s->peak_memory); +} + +static const AVFilterPad inputs[] = { + { + .name = "default", + .type = AVMEDIA_TYPE_VIDEO, + .filter_frame = filter_frame, + .config_props = config_input, + }, + { NULL } +}; + +static const AVFilterPad outputs[] = { + { + .name = "default", + .type = AVMEDIA_TYPE_VIDEO, + .config_props = config_output, + }, + { NULL } +}; + +AVFilter ff_vf_vectorscope = { + .name = "vectorscope", + .description = NULL_IF_CONFIG_SMALL("Video vectorscope."), + .priv_size = sizeof(VectorscopeContext), + .priv_class = &vectorscope_class, + .query_formats = query_formats, + .uninit = uninit, + .inputs = inputs, + .outputs = outputs, +}; |