2 * Copyright (c) 2017 Paul B Mahol
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
21 #include "libavutil/attributes.h"
22 #include "libavutil/common.h"
23 #include "libavutil/eval.h"
24 #include "libavutil/imgutils.h"
25 #include "libavutil/opt.h"
26 #include "libavutil/pixdesc.h"
32 static const char *const var_names[] = {
33 "w", ///< width of the input video
34 "h", ///< height of the input video
35 "val", ///< input value for the pixel
86 typedef struct Curve {
94 typedef struct Range {
98 typedef struct Preset {
105 static const Range full_range = {0, 256};
106 static const Range spec1_range[] = {{0, 16}, {16, 236}, {236, 256}};
107 static const Range spec2_range[] = {{0, 16}, {16, 22}, {22, 226}, {226, 236}, {236, 256}};
108 static const Range shadows_range[] = {{0, 32}, {32, 256}};
109 static const Range highlights_range[] = {{0, 214}, {214, 224}, {224, 256}};
111 static const Fill spec1_fills[] = {{{0.5f, 0.f, .5f, 1.f}}, {{-1.f, -1.f, -1.f, 1.f}}, {{1.f, 0.f, 0.f, 1.f}}};
112 static const Fill spec2_fills[] = {{{0.5f, 0.f, .5f, 1.f}}, {{0.f, 1.f, 1.f, 1.f}}, {{-1.f, -1.f, -1.f, 1.f}}, {{1.f, 1.f, 0.f, 1.f}}, {{1.f, 0.f, 0.f, 1.f}}};
113 static const Fill shadows_fills[] = {{{0.8f, 0.4f, .8f, 1.f}}, {{-1.f, -1.f, -1.f, 1.f}}};
114 static const Fill highlights_fills[] = {{{-1.f, -1.f, -1.f, 1.f}}, {{1.f, 0.3f, 0.6f, 1.f}}, {{1.f, 0.2f, .5f, 1.f}}};
116 static const Curve curves[] =
119 {-7.5631093e-16, 7.4289183e-13, -2.8525484e-10, 5.4446085e-08, -5.5596238e-06, 3.0569325e-04, -2.3137421e-03, 1.2152095e-02 },
120 { 1.3217636e-15, -1.2214648e-12, 4.4319712e-10, -8.0197993e-08, 7.6598370e-06, -3.6523704e-04, 8.4836670e-03, -2.5536888e-02 },
121 {-1.1446568e-15, 1.0013446e-12, -3.5651575e-10, 6.6775016e-08, -6.7120346e-06, 2.7346619e-04, 4.7969657e-03, 1.1971441e-02 },
124 {-3.9848859e-18, 9.4821649e-14, -6.7371977e-11, 1.8469937e-08, -2.5359307e-06, 1.7959053e-04, 3.9782564e-04, 2.8845935e-04 },
125 { 6.8408539e-16, -6.5499979e-13, 2.4562526e-10, -4.5989298e-08, 4.5723324e-06, -2.2111913e-04, 5.2023164e-03, -1.1226064e-02 },
126 {-2.9921470e-15, 2.5864165e-12, -8.7403799e-10, 1.4713388e-07, -1.2701505e-05, 4.5159935e-04, 3.1087989e-03, 1.9122831e-02 },
129 { 3.6196089e-16, -3.3623041e-13, 1.2324010e-10, -2.2769060e-08, 2.2297792e-06, -1.2567829e-04, 9.9791629e-03, 5.7247918e-02 },
130 { 5.0262888e-16, -5.3193896e-13, 2.2451715e-10, -4.7529623e-08, 5.1374873e-06, -2.3260136e-04, 3.1502825e-03, 1.5362491e-02 },
131 {-1.7782261e-16, 2.2487839e-13, -1.0610236e-10, 2.4112644e-08, -2.6331623e-06, 8.9499751e-05, 2.1386328e-03, 5.3824268e-01 },
134 { 9.4850045e-16, -8.6629383e-13, 3.0310944e-10, -5.1340396e-08, 4.6024275e-06, -2.2744239e-04, 4.5559993e-03, 2.5662350e-01 },
135 { 9.6461041e-17, -6.9209477e-14, 1.7625397e-11, -2.0229773e-09, 1.4900110e-07, -1.9315187e-05, 5.8967339e-03, 3.9544827e-03 },
136 { 5.1785449e-16, -3.6663004e-13, 1.0249990e-10, -1.5431998e-08, 1.5007941e-06, -1.2001502e-04, 7.6951526e-03, 3.2292815e-01 },
139 {-4.3683890e-15, 3.7020347e-12, -1.1712592e-09, 1.6401790e-07, -8.6842919e-06, -1.8542465e-06, 8.4485325e-03, 1.6267077e-01 },
140 {-4.0011069e-16, 2.7861423e-13, -6.3388921e-11, 5.8872238e-09, -5.4466522e-07, 1.8037114e-05, 1.0599869e-02, 7.6914696e-02 },
141 {-2.8242609e-15, 2.9234108e-12, -1.1726546e-09, 2.2552115e-07, -2.0059387e-05, 5.0595552e-04, 1.7714932e-02, 2.7271836e-01 },
144 {-9.5484131e-16, 9.6988184e-13, -4.0058766e-10, 8.5743924e-08, -9.9644797e-06, 5.9197908e-04, -1.0361579e-02, 3.3164429e-02 },
145 { 1.2731941e-17, -9.4238449e-15, 2.2808841e-12, -1.1548296e-10, -2.3888913e-08, 3.8986680e-06, 2.5879330e-03, 1.2769733e-01 },
146 { 4.6004608e-16, -5.0686849e-13, 2.2753449e-10, -5.3074099e-08, 6.7196096e-06, -4.4120020e-04, 1.3435551e-02, 2.8293355e-01 },
150 static const Preset presets[] =
152 [PRESET_MAGMA] = { 1, &full_range, &curves[MAGMA], NULL },
153 [PRESET_INFERNO] = { 1, &full_range, &curves[INFERNO], NULL },
154 [PRESET_PLASMA] = { 1, &full_range, &curves[PLASMA], NULL },
155 [PRESET_VIRIDIS] = { 1, &full_range, &curves[VIRIDIS], NULL },
156 [PRESET_TURBO] = { 1, &full_range, &curves[TURBO], NULL },
157 [PRESET_CIVIDIS] = { 1, &full_range, &curves[CIVIDIS], NULL },
158 [PRESET_RANGE1] = { 3, spec1_range, NULL, spec1_fills },
159 [PRESET_RANGE2] = { 5, spec2_range, NULL, spec2_fills },
160 [PRESET_SHADOWS] = { 2, shadows_range, NULL, shadows_fills },
161 [PRESET_HIGHLIGHTS] = { 3, highlights_range, NULL, highlights_fills },
164 typedef struct PseudoColorContext {
165 const AVClass *class;
173 int width[4], height[4];
174 double var_values[VAR_VARS_NB];
175 char *comp_expr_str[4];
176 AVExpr *comp_expr[4];
177 float lut[4][256*256];
179 void (*filter[4])(int max, int width, int height,
180 const uint8_t *index, const uint8_t *src,
187 } PseudoColorContext;
189 #define OFFSET(x) offsetof(PseudoColorContext, x)
190 #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
192 static const AVOption pseudocolor_options[] = {
193 { "c0", "set component #0 expression", OFFSET(comp_expr_str[0]), AV_OPT_TYPE_STRING, {.str="val"}, .flags = FLAGS },
194 { "c1", "set component #1 expression", OFFSET(comp_expr_str[1]), AV_OPT_TYPE_STRING, {.str="val"}, .flags = FLAGS },
195 { "c2", "set component #2 expression", OFFSET(comp_expr_str[2]), AV_OPT_TYPE_STRING, {.str="val"}, .flags = FLAGS },
196 { "c3", "set component #3 expression", OFFSET(comp_expr_str[3]), AV_OPT_TYPE_STRING, {.str="val"}, .flags = FLAGS },
197 { "index", "set component as base", OFFSET(index), AV_OPT_TYPE_INT, {.i64=0}, 0, 3, .flags = FLAGS },
198 { "i", "set component as base", OFFSET(index), AV_OPT_TYPE_INT, {.i64=0}, 0, 3, .flags = FLAGS },
199 { "preset", "set preset", OFFSET(preset), AV_OPT_TYPE_INT, {.i64=-1},-1, NB_PRESETS-1, .flags = FLAGS, "preset" },
200 { "p", "set preset", OFFSET(preset), AV_OPT_TYPE_INT, {.i64=-1},-1, NB_PRESETS-1, .flags = FLAGS, "preset" },
201 { "none", NULL, 0, AV_OPT_TYPE_CONST, {.i64=-1}, .flags = FLAGS, "preset" },
202 { "magma", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_MAGMA}, .flags = FLAGS, "preset" },
203 { "inferno", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_INFERNO}, .flags = FLAGS, "preset" },
204 { "plasma", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_PLASMA}, .flags = FLAGS, "preset" },
205 { "viridis", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_VIRIDIS}, .flags = FLAGS, "preset" },
206 { "turbo", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_TURBO}, .flags = FLAGS, "preset" },
207 { "cividis", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_CIVIDIS}, .flags = FLAGS, "preset" },
208 { "range1", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_RANGE1}, .flags = FLAGS, "preset" },
209 { "range2", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_RANGE2}, .flags = FLAGS, "preset" },
210 { "shadows", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_SHADOWS}, .flags = FLAGS, "preset" },
211 { "highlights", NULL, 0, AV_OPT_TYPE_CONST, {.i64=PRESET_HIGHLIGHTS},.flags=FLAGS, "preset" },
212 { "opacity", "set pseudocolor opacity",OFFSET(opacity), AV_OPT_TYPE_FLOAT, {.dbl=1}, 0, 1, .flags = FLAGS },
216 static const enum AVPixelFormat pix_fmts[] = {
217 AV_PIX_FMT_GRAY8, AV_PIX_FMT_GRAY9, AV_PIX_FMT_GRAY10, AV_PIX_FMT_GRAY12, AV_PIX_FMT_GRAY14, AV_PIX_FMT_GRAY16,
218 AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUVA420P,
219 AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUVA422P,
220 AV_PIX_FMT_YUV444P, AV_PIX_FMT_GBRP,
221 AV_PIX_FMT_YUVA444P, AV_PIX_FMT_GBRAP,
222 AV_PIX_FMT_YUV422P9, AV_PIX_FMT_YUVA422P9,
223 AV_PIX_FMT_YUV420P9, AV_PIX_FMT_YUVA420P9,
224 AV_PIX_FMT_YUV444P9, AV_PIX_FMT_YUVA444P9,
225 AV_PIX_FMT_YUV420P10, AV_PIX_FMT_YUVA420P10,
226 AV_PIX_FMT_YUV422P10, AV_PIX_FMT_YUVA422P10,
227 AV_PIX_FMT_YUV444P10, AV_PIX_FMT_YUVA444P10,
228 AV_PIX_FMT_YUV420P12,
229 AV_PIX_FMT_YUV422P12,
230 AV_PIX_FMT_YUV444P12,
231 AV_PIX_FMT_YUV420P14,
232 AV_PIX_FMT_YUV422P14,
233 AV_PIX_FMT_YUV444P14,
234 AV_PIX_FMT_YUV420P16, AV_PIX_FMT_YUVA420P16,
235 AV_PIX_FMT_YUV422P16, AV_PIX_FMT_YUVA422P16,
236 AV_PIX_FMT_YUV444P16, AV_PIX_FMT_YUVA444P16,
238 AV_PIX_FMT_GBRP10, AV_PIX_FMT_GBRAP10,
239 AV_PIX_FMT_GBRP12, AV_PIX_FMT_GBRAP12,
241 AV_PIX_FMT_GBRP16, AV_PIX_FMT_GBRAP16,
245 static int query_formats(AVFilterContext *ctx)
247 AVFilterFormats *fmts_list = ff_make_format_list(pix_fmts);
249 return AVERROR(ENOMEM);
250 return ff_set_common_formats(ctx, fmts_list);
253 static inline float lerpf(float v0, float v1, float f)
255 return v0 + (v1 - v0) * f;
258 #define PCLIP(v, max, dst, src, x) \
259 if (v >= 0 && v <= max) { \
260 dst[x] = lerpf(src[x], v, opacity);\
265 static void pseudocolor_filter(int max, int width, int height,
266 const uint8_t *index,
277 for (y = 0; y < height; y++) {
278 for (x = 0; x < width; x++) {
279 int v = lut[index[x]];
281 PCLIP(v, max, dst, src, x);
289 static void pseudocolor_filter_11(int max, int width, int height,
290 const uint8_t *index,
301 for (y = 0; y < height; y++) {
302 for (x = 0; x < width; x++) {
303 int v = lut[index[(y << 1) * ilinesize + (x << 1)]];
305 PCLIP(v, max, dst, src, x);
312 static void pseudocolor_filter_11d(int max, int width, int height,
313 const uint8_t *index,
324 for (y = 0; y < height; y++) {
325 for (x = 0; x < width; x++) {
326 int v = lut[index[(y >> 1) * ilinesize + (x >> 1)]];
328 PCLIP(v, max, dst, src, x);
335 static void pseudocolor_filter_10(int max, int width, int height,
336 const uint8_t *index,
347 for (y = 0; y < height; y++) {
348 for (x = 0; x < width; x++) {
349 int v = lut[index[x << 1]];
351 PCLIP(v, max, dst, src, x);
359 static void pseudocolor_filter_10d(int max, int width, int height,
360 const uint8_t *index,
371 for (y = 0; y < height; y++) {
372 for (x = 0; x < width; x++) {
373 int v = lut[index[x >> 1]];
375 PCLIP(v, max, dst, src, x);
383 static void pseudocolor_filter_16(int max, int width, int height,
384 const uint8_t *iindex,
393 const uint16_t *index = (const uint16_t *)iindex;
394 const uint16_t *src = (const uint16_t *)ssrc;
395 uint16_t *dst = (uint16_t *)ddst;
398 for (y = 0; y < height; y++) {
399 for (x = 0; x < width; x++) {
400 int v = lut[index[x]];
402 PCLIP(v, max, dst, src, x);
404 index += ilinesize / 2;
405 src += slinesize / 2;
406 dst += dlinesize / 2;
410 static void pseudocolor_filter_16_10(int max, int width, int height,
411 const uint8_t *iindex,
420 const uint16_t *index = (const uint16_t *)iindex;
421 const uint16_t *src = (const uint16_t *)ssrc;
422 uint16_t *dst = (uint16_t *)ddst;
425 for (y = 0; y < height; y++) {
426 for (x = 0; x < width; x++) {
427 int v = lut[index[x << 1]];
429 PCLIP(v, max, dst, src, x);
431 index += ilinesize / 2;
432 src += slinesize / 2;
433 dst += dlinesize / 2;
437 static void pseudocolor_filter_16_10d(int max, int width, int height,
438 const uint8_t *iindex,
447 const uint16_t *index = (const uint16_t *)iindex;
448 const uint16_t *src = (const uint16_t *)ssrc;
449 uint16_t *dst = (uint16_t *)ddst;
452 for (y = 0; y < height; y++) {
453 for (x = 0; x < width; x++) {
454 int v = lut[index[x >> 1]];
456 PCLIP(v, max, dst, src, x);
458 index += ilinesize / 2;
459 src += slinesize / 2;
460 dst += dlinesize / 2;
464 static void pseudocolor_filter_16_11(int max, int width, int height,
465 const uint8_t *iindex,
474 const uint16_t *index = (const uint16_t *)iindex;
475 const uint16_t *src = (const uint16_t *)ssrc;
476 uint16_t *dst = (uint16_t *)ddst;
483 for (y = 0; y < height; y++) {
484 for (x = 0; x < width; x++) {
485 int v = lut[index[(y << 1) * ilinesize + (x << 1)]];
487 PCLIP(v, max, dst, src, x);
494 static void pseudocolor_filter_16_11d(int max, int width, int height,
495 const uint8_t *iindex,
504 const uint16_t *index = (const uint16_t *)iindex;
505 const uint16_t *src = (const uint16_t *)ssrc;
506 uint16_t *dst = (uint16_t *)ddst;
513 for (y = 0; y < height; y++) {
514 for (x = 0; x < width; x++) {
515 int v = lut[index[(y >> 1) * ilinesize + (x >> 1)]];
517 PCLIP(v, max, dst, src, x);
524 #define RGB_TO_Y_BT709(r, g, b) \
525 ((0.21260*219.0/255.0) * (r) + (0.71520*219.0/255.0) * (g) + \
526 (0.07220*219.0/255.0) * (b))
528 #define RGB_TO_U_BT709(r1, g1, b1, max) \
529 (-(0.11457*224.0/255.0) * r1 - (0.38543*224.0/255.0) * g1 + \
530 (0.50000*224.0/255.0) * b1 + max * 0.5)
532 #define RGB_TO_V_BT709(r1, g1, b1, max) \
533 ((0.50000*224.0/255.0) * r1 - (0.45415*224.0/255.0) * g1 - \
534 (0.04585*224.0/255.0) * b1 + max * 0.5)
536 static double poly_eval(const double *const poly, double x)
540 for (int i = 0; i < 8; i++) {
541 res += pow(x, i) * poly[7-i];
544 return av_clipd(res, 0., 1.);
547 static int config_input(AVFilterLink *inlink)
549 AVFilterContext *ctx = inlink->dst;
550 PseudoColorContext *s = ctx->priv;
551 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
552 int depth, ret, hsub, vsub, color, factor, rgb;
554 rgb = desc->flags & AV_PIX_FMT_FLAG_RGB;
555 depth = desc->comp[0].depth;
556 factor = 1 << (depth - 8);
557 s->max = (1 << depth) - 1;
558 s->nb_planes = av_pix_fmt_count_planes(inlink->format);
560 if (s->index >= s->nb_planes) {
561 av_log(ctx, AV_LOG_ERROR, "index out of allowed range\n");
562 return AVERROR(EINVAL);
565 if ((ret = av_image_fill_linesizes(s->linesize, inlink->format, inlink->w)) < 0)
568 hsub = desc->log2_chroma_w;
569 vsub = desc->log2_chroma_h;
570 s->height[1] = s->height[2] = AV_CEIL_RSHIFT(inlink->h, vsub);
571 s->height[0] = s->height[3] = inlink->h;
572 s->width[1] = s->width[2] = AV_CEIL_RSHIFT(inlink->w, hsub);
573 s->width[0] = s->width[3] = inlink->w;
575 s->var_values[VAR_W] = inlink->w;
576 s->var_values[VAR_H] = inlink->h;
578 s->var_values[VAR_YMIN] = 16 * (1 << (depth - 8));
579 s->var_values[VAR_UMIN] = 16 * (1 << (depth - 8));
580 s->var_values[VAR_VMIN] = 16 * (1 << (depth - 8));
581 s->var_values[VAR_AMIN] = 0;
582 s->var_values[VAR_YMAX] = 235 * (1 << (depth - 8));
583 s->var_values[VAR_UMAX] = 240 * (1 << (depth - 8));
584 s->var_values[VAR_VMAX] = 240 * (1 << (depth - 8));
585 s->var_values[VAR_AMAX] = s->max;
587 for (color = 0; color < s->nb_planes && s->preset < 0; color++) {
591 /* create the parsed expression */
592 av_expr_free(s->comp_expr[color]);
593 s->comp_expr[color] = NULL;
594 ret = av_expr_parse(&s->comp_expr[color], s->comp_expr_str[color],
595 var_names, NULL, NULL, NULL, NULL, 0, ctx);
597 av_log(ctx, AV_LOG_ERROR,
598 "Error when parsing the expression '%s' for the component %d and color %d.\n",
599 s->comp_expr_str[color], color, color);
600 return AVERROR(EINVAL);
603 /* compute the lut */
604 for (val = 0; val < FF_ARRAY_ELEMS(s->lut[color]); val++) {
605 s->var_values[VAR_VAL] = val;
607 res = av_expr_eval(s->comp_expr[color], s->var_values, s);
609 av_log(ctx, AV_LOG_ERROR,
610 "Error when evaluating the expression '%s' for the value %d for the component %d.\n",
611 s->comp_expr_str[color], val, color);
612 return AVERROR(EINVAL);
614 s->lut[color][val] = res;
618 if (s->preset >= 0) {
619 int nb_segments = presets[s->preset].nb_segments;
621 for (int seg = 0; seg < nb_segments; seg++) {
622 int start = presets[s->preset].ranges[seg].start;
623 int end = presets[s->preset].ranges[seg].end;
625 for (int i = start; i < end; i++) {
626 if (!presets[s->preset].curves) {
627 const Fill fill = presets[s->preset].fills[seg];
629 for (int j = 0; j < factor; j++) {
637 if (g >= 0.f && b >= 0.f && r >= 0.f) {
643 double y = RGB_TO_Y_BT709(r, g, b);
644 double u = RGB_TO_U_BT709(r, g, b, s->max);
645 double v = RGB_TO_V_BT709(r, g, b, s->max);
653 s->lut[0][i*factor+j] = g;
654 s->lut[1][i*factor+j] = b;
655 s->lut[2][i*factor+j] = r;
656 s->lut[3][i*factor+j] = a * s->max;
659 const Curve curve = presets[s->preset].curves[seg];
661 for (int j = 0; j < factor; j++) {
662 const double lf = j / (double)factor;
665 g = poly_eval(curve.coef[1], i + lf) * s->max;
666 b = poly_eval(curve.coef[2], i + lf) * s->max;
667 r = poly_eval(curve.coef[0], i + lf) * s->max;
670 double y = RGB_TO_Y_BT709(r, g, b);
671 double u = RGB_TO_U_BT709(r, g, b, s->max);
672 double v = RGB_TO_V_BT709(r, g, b, s->max);
679 s->lut[0][i*factor+j] = g;
680 s->lut[1][i*factor+j] = b;
681 s->lut[2][i*factor+j] = r;
682 s->lut[3][i*factor+j] = 1.f * s->max;
689 switch (inlink->format) {
690 case AV_PIX_FMT_YUV444P:
691 case AV_PIX_FMT_YUVA444P:
692 case AV_PIX_FMT_GBRP:
693 case AV_PIX_FMT_GBRAP:
694 case AV_PIX_FMT_GRAY8:
695 s->filter[0] = s->filter[1] = s->filter[2] = s->filter[3] = pseudocolor_filter;
697 case AV_PIX_FMT_YUV420P:
698 case AV_PIX_FMT_YUVA420P:
702 s->filter[0] = s->filter[3] = pseudocolor_filter;
703 s->filter[1] = s->filter[2] = pseudocolor_filter_11;
707 s->filter[0] = s->filter[3] = pseudocolor_filter_11d;
708 s->filter[1] = s->filter[2] = pseudocolor_filter;
712 case AV_PIX_FMT_YUV422P:
713 case AV_PIX_FMT_YUVA422P:
717 s->filter[0] = s->filter[3] = pseudocolor_filter;
718 s->filter[1] = s->filter[2] = pseudocolor_filter_10;
722 s->filter[0] = s->filter[3] = pseudocolor_filter_10d;
723 s->filter[1] = s->filter[2] = pseudocolor_filter;
727 case AV_PIX_FMT_YUV444P9:
728 case AV_PIX_FMT_YUVA444P9:
729 case AV_PIX_FMT_YUV444P10:
730 case AV_PIX_FMT_YUVA444P10:
731 case AV_PIX_FMT_YUV444P12:
732 case AV_PIX_FMT_YUV444P14:
733 case AV_PIX_FMT_YUV444P16:
734 case AV_PIX_FMT_YUVA444P16:
735 case AV_PIX_FMT_GBRP9:
736 case AV_PIX_FMT_GBRP10:
737 case AV_PIX_FMT_GBRP12:
738 case AV_PIX_FMT_GBRP14:
739 case AV_PIX_FMT_GBRP16:
740 case AV_PIX_FMT_GBRAP10:
741 case AV_PIX_FMT_GBRAP12:
742 case AV_PIX_FMT_GBRAP16:
743 case AV_PIX_FMT_GRAY9:
744 case AV_PIX_FMT_GRAY10:
745 case AV_PIX_FMT_GRAY12:
746 case AV_PIX_FMT_GRAY14:
747 case AV_PIX_FMT_GRAY16:
748 s->filter[0] = s->filter[1] = s->filter[2] = s->filter[3] = pseudocolor_filter_16;
750 case AV_PIX_FMT_YUV422P9:
751 case AV_PIX_FMT_YUVA422P9:
752 case AV_PIX_FMT_YUV422P10:
753 case AV_PIX_FMT_YUVA422P10:
754 case AV_PIX_FMT_YUV422P12:
755 case AV_PIX_FMT_YUV422P14:
756 case AV_PIX_FMT_YUV422P16:
757 case AV_PIX_FMT_YUVA422P16:
761 s->filter[0] = s->filter[3] = pseudocolor_filter_16;
762 s->filter[1] = s->filter[2] = pseudocolor_filter_16_10;
766 s->filter[0] = s->filter[3] = pseudocolor_filter_16_10d;
767 s->filter[1] = s->filter[2] = pseudocolor_filter_16;
771 case AV_PIX_FMT_YUV420P9:
772 case AV_PIX_FMT_YUVA420P9:
773 case AV_PIX_FMT_YUV420P10:
774 case AV_PIX_FMT_YUVA420P10:
775 case AV_PIX_FMT_YUV420P12:
776 case AV_PIX_FMT_YUV420P14:
777 case AV_PIX_FMT_YUV420P16:
778 case AV_PIX_FMT_YUVA420P16:
782 s->filter[0] = s->filter[3] = pseudocolor_filter_16;
783 s->filter[1] = s->filter[2] = pseudocolor_filter_16_11;
787 s->filter[0] = s->filter[3] = pseudocolor_filter_16_11d;
788 s->filter[1] = s->filter[2] = pseudocolor_filter_16;
797 typedef struct ThreadData {
801 static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
803 PseudoColorContext *s = ctx->priv;
804 ThreadData *td = arg;
805 AVFrame *in = td->in;
806 AVFrame *out = td->out;
808 for (int plane = 0; plane < s->nb_planes; plane++) {
809 const int slice_start = (s->height[plane] * jobnr) / nb_jobs;
810 const int slice_end = (s->height[plane] * (jobnr+1)) / nb_jobs;
811 const int islice_start = (s->height[s->index] * jobnr) / nb_jobs;
812 ptrdiff_t ilinesize = in->linesize[s->index];
813 ptrdiff_t slinesize = in->linesize[plane];
814 ptrdiff_t dlinesize = out->linesize[plane];
815 const uint8_t *index = in->data[s->index] + islice_start * ilinesize;
816 const uint8_t *src = in->data[plane] + slice_start * slinesize;
817 uint8_t *dst = out->data[plane] + slice_start * dlinesize;
819 s->filter[plane](s->max, s->width[plane], slice_end - slice_start,
820 index, src, dst, ilinesize, slinesize,
821 dlinesize, s->lut[plane], s->opacity);
827 static int filter_frame(AVFilterLink *inlink, AVFrame *in)
829 AVFilterContext *ctx = inlink->dst;
830 PseudoColorContext *s = ctx->priv;
831 AVFilterLink *outlink = ctx->outputs[0];
835 out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
838 return AVERROR(ENOMEM);
840 av_frame_copy_props(out, in);
842 td.out = out, td.in = in;
843 ctx->internal->execute(ctx, filter_slice, &td, NULL, FFMIN(s->height[1], ff_filter_get_nb_threads(ctx)));
846 return ff_filter_frame(outlink, out);
849 static int process_command(AVFilterContext *ctx, const char *cmd, const char *args,
850 char *res, int res_len, int flags)
852 int ret = ff_filter_process_command(ctx, cmd, args, res, res_len, flags);
857 return config_input(ctx->inputs[0]);
860 static const AVFilterPad inputs[] = {
863 .type = AVMEDIA_TYPE_VIDEO,
864 .filter_frame = filter_frame,
865 .config_props = config_input,
870 static const AVFilterPad outputs[] = {
873 .type = AVMEDIA_TYPE_VIDEO,
878 static av_cold void uninit(AVFilterContext *ctx)
880 PseudoColorContext *s = ctx->priv;
883 for (i = 0; i < 4; i++) {
884 av_expr_free(s->comp_expr[i]);
885 s->comp_expr[i] = NULL;
889 AVFILTER_DEFINE_CLASS(pseudocolor);
891 AVFilter ff_vf_pseudocolor = {
892 .name = "pseudocolor",
893 .description = NULL_IF_CONFIG_SMALL("Make pseudocolored video frames."),
894 .priv_size = sizeof(PseudoColorContext),
895 .priv_class = &pseudocolor_class,
897 .query_formats = query_formats,
900 .flags = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | AVFILTER_FLAG_SLICE_THREADS,
901 .process_command = process_command,