GCC Code Coverage Report
Directory: ../../../ffmpeg/ Exec Total Coverage
File: src/libavfilter/vf_colorchannelmixer.c Lines: 112 451 24.8 %
Date: 2021-01-27 11:56:10 Branches: 27 164 16.5 %

Line Branch Exec Source
1
/*
2
 * Copyright (c) 2013 Paul B Mahol
3
 *
4
 * This file is part of FFmpeg.
5
 *
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.
10
 *
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.
15
 *
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
19
 */
20
21
#include <float.h>
22
23
#include "libavutil/opt.h"
24
#include "libavutil/pixdesc.h"
25
#include "avfilter.h"
26
#include "drawutils.h"
27
#include "formats.h"
28
#include "internal.h"
29
#include "video.h"
30
31
#define R 0
32
#define G 1
33
#define B 2
34
#define A 3
35
36
typedef struct ThreadData {
37
    AVFrame *in, *out;
38
} ThreadData;
39
40
typedef struct ColorChannelMixerContext {
41
    const AVClass *class;
42
    double rr, rg, rb, ra;
43
    double gr, gg, gb, ga;
44
    double br, bg, bb, ba;
45
    double ar, ag, ab, aa;
46
    double sr, sg, sb;
47
    double preserve_lightness;
48
49
    int *lut[4][4];
50
51
    int *buffer;
52
53
    uint8_t rgba_map[4];
54
55
    int (*filter_slice[2])(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs);
56
} ColorChannelMixerContext;
57
58
#define OFFSET(x) offsetof(ColorChannelMixerContext, x)
59
#define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
60
61
static const AVOption colorchannelmixer_options[] = {
62
    { "rr", "set the red gain for the red channel",     OFFSET(rr), AV_OPT_TYPE_DOUBLE, {.dbl=1}, -2, 2, FLAGS },
63
    { "rg", "set the green gain for the red channel",   OFFSET(rg), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
64
    { "rb", "set the blue gain for the red channel",    OFFSET(rb), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
65
    { "ra", "set the alpha gain for the red channel",   OFFSET(ra), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
66
    { "gr", "set the red gain for the green channel",   OFFSET(gr), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
67
    { "gg", "set the green gain for the green channel", OFFSET(gg), AV_OPT_TYPE_DOUBLE, {.dbl=1}, -2, 2, FLAGS },
68
    { "gb", "set the blue gain for the green channel",  OFFSET(gb), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
69
    { "ga", "set the alpha gain for the green channel", OFFSET(ga), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
70
    { "br", "set the red gain for the blue channel",    OFFSET(br), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
71
    { "bg", "set the green gain for the blue channel",  OFFSET(bg), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
72
    { "bb", "set the blue gain for the blue channel",   OFFSET(bb), AV_OPT_TYPE_DOUBLE, {.dbl=1}, -2, 2, FLAGS },
73
    { "ba", "set the alpha gain for the blue channel",  OFFSET(ba), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
74
    { "ar", "set the red gain for the alpha channel",   OFFSET(ar), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
75
    { "ag", "set the green gain for the alpha channel", OFFSET(ag), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
76
    { "ab", "set the blue gain for the alpha channel",  OFFSET(ab), AV_OPT_TYPE_DOUBLE, {.dbl=0}, -2, 2, FLAGS },
77
    { "aa", "set the alpha gain for the alpha channel", OFFSET(aa), AV_OPT_TYPE_DOUBLE, {.dbl=1}, -2, 2, FLAGS },
78
    { "pl", "preserve lightness",       OFFSET(preserve_lightness), AV_OPT_TYPE_DOUBLE, {.dbl=0},  0, 1, FLAGS },
79
    { NULL }
80
};
81
82
AVFILTER_DEFINE_CLASS(colorchannelmixer);
83
84
1
static int query_formats(AVFilterContext *ctx)
85
{
86
    static const enum AVPixelFormat pix_fmts[] = {
87
        AV_PIX_FMT_RGB24,  AV_PIX_FMT_BGR24,
88
        AV_PIX_FMT_RGBA,   AV_PIX_FMT_BGRA,
89
        AV_PIX_FMT_ARGB,   AV_PIX_FMT_ABGR,
90
        AV_PIX_FMT_0RGB,   AV_PIX_FMT_0BGR,
91
        AV_PIX_FMT_RGB0,   AV_PIX_FMT_BGR0,
92
        AV_PIX_FMT_RGB48,  AV_PIX_FMT_BGR48,
93
        AV_PIX_FMT_RGBA64, AV_PIX_FMT_BGRA64,
94
        AV_PIX_FMT_GBRP,   AV_PIX_FMT_GBRAP,
95
        AV_PIX_FMT_GBRP9,
96
        AV_PIX_FMT_GBRP10, AV_PIX_FMT_GBRAP10,
97
        AV_PIX_FMT_GBRP12, AV_PIX_FMT_GBRAP12,
98
        AV_PIX_FMT_GBRP14,
99
        AV_PIX_FMT_GBRP16, AV_PIX_FMT_GBRAP16,
100
        AV_PIX_FMT_NONE
101
    };
102
103
1
    AVFilterFormats *fmts_list = ff_make_format_list(pix_fmts);
104
1
    if (!fmts_list)
105
        return AVERROR(ENOMEM);
106
1
    return ff_set_common_formats(ctx, fmts_list);
107
}
108
109
static float lerpf(float v0, float v1, float f)
110
{
111
    return v0 + (v1 - v0) * f;
112
}
113
114
static void preservel(float *r, float *g, float *b, float lin, float lout)
115
{
116
    *r *= lout / lin;
117
    *g *= lout / lin;
118
    *b *= lout / lin;
119
}
120
121
static av_always_inline int filter_slice_rgba_planar(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs,
122
                                                     int have_alpha, int pl)
123
{
124
    ColorChannelMixerContext *s = ctx->priv;
125
    ThreadData *td = arg;
126
    AVFrame *in = td->in;
127
    AVFrame *out = td->out;
128
    const float l = s->preserve_lightness;
129
    const float sr = s->sr;
130
    const float sg = s->sg;
131
    const float sb = s->sb;
132
    const int slice_start = (out->height * jobnr) / nb_jobs;
133
    const int slice_end = (out->height * (jobnr+1)) / nb_jobs;
134
    const uint8_t *srcg = in->data[0] + slice_start * in->linesize[0];
135
    const uint8_t *srcb = in->data[1] + slice_start * in->linesize[1];
136
    const uint8_t *srcr = in->data[2] + slice_start * in->linesize[2];
137
    const uint8_t *srca = in->data[3] + slice_start * in->linesize[3];
138
    uint8_t *dstg = out->data[0] + slice_start * out->linesize[0];
139
    uint8_t *dstb = out->data[1] + slice_start * out->linesize[1];
140
    uint8_t *dstr = out->data[2] + slice_start * out->linesize[2];
141
    uint8_t *dsta = out->data[3] + slice_start * out->linesize[3];
142
    int i, j;
143
144
    for (i = slice_start; i < slice_end; i++) {
145
        for (j = 0; j < out->width; j++) {
146
            const uint8_t rin = srcr[j];
147
            const uint8_t gin = srcg[j];
148
            const uint8_t bin = srcb[j];
149
            const uint8_t ain = have_alpha ? srca[j] : 0;
150
            int rout, gout, bout;
151
            float lin;
152
153
            if (pl)
154
                lin = FFMAX3(rin, gin, bin) + FFMIN3(rin, gin, bin);
155
156
            rout = s->lut[R][R][rin] +
157
                   s->lut[R][G][gin] +
158
                   s->lut[R][B][bin] +
159
                   (have_alpha == 1 ? s->lut[R][A][ain] : 0);
160
            gout = s->lut[G][R][rin] +
161
                   s->lut[G][G][gin] +
162
                   s->lut[G][B][bin] +
163
                   (have_alpha == 1 ? s->lut[G][A][ain] : 0);
164
            bout = s->lut[B][R][rin] +
165
                   s->lut[B][G][gin] +
166
                   s->lut[B][B][bin] +
167
                   (have_alpha == 1 ? s->lut[B][A][ain] : 0);
168
169
            if (pl) {
170
                float frout = rout / sr;
171
                float fgout = gout / sg;
172
                float fbout = bout / sb;
173
                float lout = FFMAX3(frout, fgout, fbout) + FFMIN3(frout, fgout, fbout);
174
175
                preservel(&frout, &fgout, &fbout, lin, lout);
176
177
                rout = lrintf(lerpf(rout, frout, l));
178
                gout = lrintf(lerpf(gout, fgout, l));
179
                bout = lrintf(lerpf(bout, fbout, l));
180
            }
181
182
            dstr[j] = av_clip_uint8(rout);
183
            dstg[j] = av_clip_uint8(gout);
184
            dstb[j] = av_clip_uint8(bout);
185
186
            if (have_alpha == 1) {
187
                dsta[j] = av_clip_uint8(s->lut[A][R][rin] +
188
                                        s->lut[A][G][gin] +
189
                                        s->lut[A][B][bin] +
190
                                        s->lut[A][A][ain]);
191
            }
192
        }
193
194
        srcg += in->linesize[0];
195
        srcb += in->linesize[1];
196
        srcr += in->linesize[2];
197
        srca += in->linesize[3];
198
        dstg += out->linesize[0];
199
        dstb += out->linesize[1];
200
        dstr += out->linesize[2];
201
        dsta += out->linesize[3];
202
    }
203
204
    return 0;
205
}
206
207
static av_always_inline int filter_slice_rgba16_planar(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs,
208
                                                       int have_alpha, int depth, int pl)
209
{
210
    ColorChannelMixerContext *s = ctx->priv;
211
    ThreadData *td = arg;
212
    AVFrame *in = td->in;
213
    AVFrame *out = td->out;
214
    const float l = s->preserve_lightness;
215
    const float sr = s->sr;
216
    const float sg = s->sg;
217
    const float sb = s->sb;
218
    const int slice_start = (out->height * jobnr) / nb_jobs;
219
    const int slice_end = (out->height * (jobnr+1)) / nb_jobs;
220
    const uint16_t *srcg = (const uint16_t *)(in->data[0] + slice_start * in->linesize[0]);
221
    const uint16_t *srcb = (const uint16_t *)(in->data[1] + slice_start * in->linesize[1]);
222
    const uint16_t *srcr = (const uint16_t *)(in->data[2] + slice_start * in->linesize[2]);
223
    const uint16_t *srca = (const uint16_t *)(in->data[3] + slice_start * in->linesize[3]);
224
    uint16_t *dstg = (uint16_t *)(out->data[0] + slice_start * out->linesize[0]);
225
    uint16_t *dstb = (uint16_t *)(out->data[1] + slice_start * out->linesize[1]);
226
    uint16_t *dstr = (uint16_t *)(out->data[2] + slice_start * out->linesize[2]);
227
    uint16_t *dsta = (uint16_t *)(out->data[3] + slice_start * out->linesize[3]);
228
    int i, j;
229
230
    for (i = slice_start; i < slice_end; i++) {
231
        for (j = 0; j < out->width; j++) {
232
            const uint16_t rin = srcr[j];
233
            const uint16_t gin = srcg[j];
234
            const uint16_t bin = srcb[j];
235
            const uint16_t ain = have_alpha ? srca[j] : 0;
236
            int rout, gout, bout;
237
            float lin;
238
239
            if (pl)
240
                lin = FFMAX3(rin, gin, bin) + FFMIN3(rin, gin, bin);
241
242
            rout = s->lut[R][R][rin] +
243
                   s->lut[R][G][gin] +
244
                   s->lut[R][B][bin] +
245
                   (have_alpha == 1 ? s->lut[R][A][ain] : 0);
246
            gout = s->lut[G][R][rin] +
247
                   s->lut[G][G][gin] +
248
                   s->lut[G][B][bin] +
249
                   (have_alpha == 1 ? s->lut[G][A][ain] : 0);
250
            bout = s->lut[B][R][rin] +
251
                   s->lut[B][G][gin] +
252
                   s->lut[B][B][bin] +
253
                   (have_alpha == 1 ? s->lut[B][A][ain] : 0);
254
255
            if (pl) {
256
                float frout = rout / sr;
257
                float fgout = gout / sg;
258
                float fbout = bout / sb;
259
                float lout = FFMAX3(frout, fgout, fbout) + FFMIN3(frout, fgout, fbout);
260
261
                preservel(&frout, &fgout, &fbout, lin, lout);
262
263
                rout = lrintf(lerpf(rout, frout, l));
264
                gout = lrintf(lerpf(gout, fgout, l));
265
                bout = lrintf(lerpf(bout, fbout, l));
266
            }
267
268
            dstr[j] = av_clip_uintp2(rout, depth);
269
            dstg[j] = av_clip_uintp2(gout, depth);
270
            dstb[j] = av_clip_uintp2(bout, depth);
271
272
            if (have_alpha == 1) {
273
                dsta[j] = av_clip_uintp2(s->lut[A][R][rin] +
274
                                         s->lut[A][G][gin] +
275
                                         s->lut[A][B][bin] +
276
                                         s->lut[A][A][ain], depth);
277
            }
278
        }
279
280
        srcg += in->linesize[0] / 2;
281
        srcb += in->linesize[1] / 2;
282
        srcr += in->linesize[2] / 2;
283
        srca += in->linesize[3] / 2;
284
        dstg += out->linesize[0] / 2;
285
        dstb += out->linesize[1] / 2;
286
        dstr += out->linesize[2] / 2;
287
        dsta += out->linesize[3] / 2;
288
    }
289
290
    return 0;
291
}
292
293
static int filter_slice_gbrp(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
294
{
295
    return filter_slice_rgba_planar(ctx, arg, jobnr, nb_jobs, 0, 0);
296
}
297
298
static int filter_slice_gbrap(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
299
{
300
    return filter_slice_rgba_planar(ctx, arg, jobnr, nb_jobs, 1, 0);
301
}
302
303
static int filter_slice_gbrp_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
304
{
305
    return filter_slice_rgba_planar(ctx, arg, jobnr, nb_jobs, 0, 1);
306
}
307
308
static int filter_slice_gbrap_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
309
{
310
    return filter_slice_rgba_planar(ctx, arg, jobnr, nb_jobs, 1, 1);
311
}
312
313
static int filter_slice_gbrp9(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
314
{
315
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 9, 0);
316
}
317
318
static int filter_slice_gbrp10(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
319
{
320
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 10, 0);
321
}
322
323
static int filter_slice_gbrap10(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
324
{
325
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 1, 10, 0);
326
}
327
328
static int filter_slice_gbrp12(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
329
{
330
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 12, 0);
331
}
332
333
static int filter_slice_gbrap12(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
334
{
335
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 1, 12, 0);
336
}
337
338
static int filter_slice_gbrp14(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
339
{
340
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 14, 0);
341
}
342
343
static int filter_slice_gbrp16(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
344
{
345
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 16, 0);
346
}
347
348
static int filter_slice_gbrap16(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
349
{
350
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 1, 16, 0);
351
}
352
353
static int filter_slice_gbrp9_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
354
{
355
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 9, 1);
356
}
357
358
static int filter_slice_gbrp10_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
359
{
360
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 10, 1);
361
}
362
363
static int filter_slice_gbrap10_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
364
{
365
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 1, 10, 1);
366
}
367
368
static int filter_slice_gbrp12_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
369
{
370
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 12, 1);
371
}
372
373
static int filter_slice_gbrap12_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
374
{
375
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 1, 12, 1);
376
}
377
378
static int filter_slice_gbrp14_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
379
{
380
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 14, 1);
381
}
382
383
static int filter_slice_gbrp16_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
384
{
385
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 0, 16, 1);
386
}
387
388
static int filter_slice_gbrap16_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
389
{
390
    return filter_slice_rgba16_planar(ctx, arg, jobnr, nb_jobs, 1, 16, 1);
391
}
392
393
450
static av_always_inline int filter_slice_rgba_packed(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs,
394
                                                     int have_alpha, int step, int pl)
395
{
396
450
    ColorChannelMixerContext *s = ctx->priv;
397
450
    ThreadData *td = arg;
398
450
    AVFrame *in = td->in;
399
450
    AVFrame *out = td->out;
400
450
    const float l = s->preserve_lightness;
401
450
    const float sr = s->sr;
402
450
    const float sg = s->sg;
403
450
    const float sb = s->sb;
404
450
    const int slice_start = (out->height * jobnr) / nb_jobs;
405
450
    const int slice_end = (out->height * (jobnr+1)) / nb_jobs;
406
450
    const uint8_t roffset = s->rgba_map[R];
407
450
    const uint8_t goffset = s->rgba_map[G];
408
450
    const uint8_t boffset = s->rgba_map[B];
409
450
    const uint8_t aoffset = s->rgba_map[A];
410
450
    const uint8_t *srcrow = in->data[0] + slice_start * in->linesize[0];
411
450
    uint8_t *dstrow = out->data[0] + slice_start * out->linesize[0];
412
    int i, j;
413
414
14850
    for (i = slice_start; i < slice_end; i++) {
415
14400
        const uint8_t *src = srcrow;
416
14400
        uint8_t *dst = dstrow;
417
418
5083200
        for (j = 0; j < out->width * step; j += step) {
419
5068800
            const uint8_t rin = src[j + roffset];
420
5068800
            const uint8_t gin = src[j + goffset];
421
5068800
            const uint8_t bin = src[j + boffset];
422
5068800
            const uint8_t ain = src[j + aoffset];
423
            int rout, gout, bout;
424
            float lin;
425
426
5068800
            if (pl)
427
                lin = FFMAX3(rin, gin, bin) + FFMIN3(rin, gin, bin);
428
429
10137600
            rout = s->lut[R][R][rin] +
430
5068800
                   s->lut[R][G][gin] +
431
5068800
                   s->lut[R][B][bin] +
432
5068800
                   (have_alpha == 1 ? s->lut[R][A][ain] : 0);
433
10137600
            gout = s->lut[G][R][rin] +
434
5068800
                   s->lut[G][G][gin] +
435
5068800
                   s->lut[G][B][bin] +
436
5068800
                   (have_alpha == 1 ? s->lut[G][A][ain] : 0);
437
10137600
            bout = s->lut[B][R][rin] +
438
5068800
                   s->lut[B][G][gin] +
439
5068800
                   s->lut[B][B][bin] +
440
5068800
                   (have_alpha == 1 ? s->lut[B][A][ain] : 0);
441
442
5068800
            if (pl) {
443
                float frout = rout / sr;
444
                float fgout = gout / sg;
445
                float fbout = bout / sb;
446
                float lout = FFMAX3(frout, fgout, fbout) + FFMIN3(frout, fgout, fbout);
447
448
                preservel(&frout, &fgout, &fbout, lin, lout);
449
450
                rout = lrintf(lerpf(rout, frout, l));
451
                gout = lrintf(lerpf(gout, fgout, l));
452
                bout = lrintf(lerpf(bout, fbout, l));
453
            }
454
455
5068800
            dst[j + roffset] = av_clip_uint8(rout);
456
5068800
            dst[j + goffset] = av_clip_uint8(gout);
457
5068800
            dst[j + boffset] = av_clip_uint8(bout);
458
459
5068800
            if (have_alpha == 1) {
460
                dst[j + aoffset] = av_clip_uint8(s->lut[A][R][rin] +
461
                                                 s->lut[A][G][gin] +
462
                                                 s->lut[A][B][bin] +
463
                                                 s->lut[A][A][ain]);
464

5068800
            } else if (have_alpha == -1 && in != out)
465
                dst[j + aoffset] = 0;
466
        }
467
468
14400
        srcrow += in->linesize[0];
469
14400
        dstrow += out->linesize[0];
470
    }
471
472
450
    return 0;
473
}
474
475
static av_always_inline int filter_slice_rgba16_packed(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs,
476
                                                       int have_alpha, int step, int pl)
477
{
478
    ColorChannelMixerContext *s = ctx->priv;
479
    ThreadData *td = arg;
480
    AVFrame *in = td->in;
481
    AVFrame *out = td->out;
482
    const float l = s->preserve_lightness;
483
    const float sr = s->sr;
484
    const float sg = s->sg;
485
    const float sb = s->sb;
486
    const int slice_start = (out->height * jobnr) / nb_jobs;
487
    const int slice_end = (out->height * (jobnr+1)) / nb_jobs;
488
    const uint8_t roffset = s->rgba_map[R];
489
    const uint8_t goffset = s->rgba_map[G];
490
    const uint8_t boffset = s->rgba_map[B];
491
    const uint8_t aoffset = s->rgba_map[A];
492
    const uint8_t *srcrow = in->data[0] + slice_start * in->linesize[0];
493
    uint8_t *dstrow = out->data[0] + slice_start * out->linesize[0];
494
    int i, j;
495
496
    for (i = slice_start; i < slice_end; i++) {
497
        const uint16_t *src = (const uint16_t *)srcrow;
498
        uint16_t *dst = (uint16_t *)dstrow;
499
500
        for (j = 0; j < out->width * step; j += step) {
501
            const uint16_t rin = src[j + roffset];
502
            const uint16_t gin = src[j + goffset];
503
            const uint16_t bin = src[j + boffset];
504
            const uint16_t ain = src[j + aoffset];
505
            int rout, gout, bout;
506
            float lin;
507
508
            if (pl)
509
                lin = FFMAX3(rin, gin, bin) + FFMIN3(rin, gin, bin);
510
511
            rout = s->lut[R][R][rin] +
512
                   s->lut[R][G][gin] +
513
                   s->lut[R][B][bin] +
514
                   (have_alpha == 1 ? s->lut[R][A][ain] : 0);
515
            gout = s->lut[G][R][rin] +
516
                   s->lut[G][G][gin] +
517
                   s->lut[G][B][bin] +
518
                   (have_alpha == 1 ? s->lut[G][A][ain] : 0);
519
            bout = s->lut[B][R][rin] +
520
                   s->lut[B][G][gin] +
521
                   s->lut[B][B][bin] +
522
                   (have_alpha == 1 ? s->lut[B][A][ain] : 0);
523
524
            if (pl) {
525
                float frout = rout / sr;
526
                float fgout = gout / sg;
527
                float fbout = bout / sb;
528
                float lout = FFMAX3(frout, fgout, fbout) + FFMIN3(frout, fgout, fbout);
529
530
                preservel(&frout, &fgout, &fbout, lin, lout);
531
532
                rout = lrintf(lerpf(rout, frout, l));
533
                gout = lrintf(lerpf(gout, fgout, l));
534
                bout = lrintf(lerpf(bout, fbout, l));
535
            }
536
537
            dst[j + roffset] = av_clip_uint16(rout);
538
            dst[j + goffset] = av_clip_uint16(gout);
539
            dst[j + boffset] = av_clip_uint16(bout);
540
541
            if (have_alpha == 1) {
542
                dst[j + aoffset] = av_clip_uint16(s->lut[A][R][rin] +
543
                                                  s->lut[A][G][gin] +
544
                                                  s->lut[A][B][bin] +
545
                                                  s->lut[A][A][ain]);
546
            }
547
        }
548
549
        srcrow += in->linesize[0];
550
        dstrow += out->linesize[0];
551
    }
552
553
    return 0;
554
}
555
556
static int filter_slice_rgba64(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
557
{
558
    return filter_slice_rgba16_packed(ctx, arg, jobnr, nb_jobs, 1, 4, 0);
559
}
560
561
static int filter_slice_rgb48(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
562
{
563
    return filter_slice_rgba16_packed(ctx, arg, jobnr, nb_jobs, 0, 3, 0);
564
}
565
566
static int filter_slice_rgba64_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
567
{
568
    return filter_slice_rgba16_packed(ctx, arg, jobnr, nb_jobs, 1, 4, 1);
569
}
570
571
static int filter_slice_rgb48_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
572
{
573
    return filter_slice_rgba16_packed(ctx, arg, jobnr, nb_jobs, 0, 3, 1);
574
}
575
576
static int filter_slice_rgba(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
577
{
578
    return filter_slice_rgba_packed(ctx, arg, jobnr, nb_jobs, 1, 4, 0);
579
}
580
581
450
static int filter_slice_rgb24(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
582
{
583
450
    return filter_slice_rgba_packed(ctx, arg, jobnr, nb_jobs, 0, 3, 0);
584
}
585
586
static int filter_slice_rgb0(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
587
{
588
    return filter_slice_rgba_packed(ctx, arg, jobnr, nb_jobs, -1, 4, 0);
589
}
590
591
static int filter_slice_rgba_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
592
{
593
    return filter_slice_rgba_packed(ctx, arg, jobnr, nb_jobs, 1, 4, 1);
594
}
595
596
static int filter_slice_rgb24_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
597
{
598
    return filter_slice_rgba_packed(ctx, arg, jobnr, nb_jobs, 0, 3, 1);
599
}
600
601
static int filter_slice_rgb0_pl(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
602
{
603
    return filter_slice_rgba_packed(ctx, arg, jobnr, nb_jobs, -1, 4, 1);
604
}
605
606
1
static int config_output(AVFilterLink *outlink)
607
{
608
1
    AVFilterContext *ctx = outlink->src;
609
1
    ColorChannelMixerContext *s = ctx->priv;
610
1
    const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(outlink->format);
611
1
    const int depth = desc->comp[0].depth;
612
1
    int i, j, size, *buffer = s->buffer;
613
614
1
    ff_fill_rgba_map(s->rgba_map, outlink->format);
615
616
1
    size = 1 << depth;
617
1
    if (!s->buffer) {
618
1
        s->buffer = buffer = av_malloc(16 * size * sizeof(*s->buffer));
619
1
        if (!s->buffer)
620
            return AVERROR(ENOMEM);
621
622
5
        for (i = 0; i < 4; i++)
623
20
            for (j = 0; j < 4; j++, buffer += size)
624
16
                s->lut[i][j] = buffer;
625
    }
626
627
1
    s->sr = s->rr + s->rg + s->rb + s->ra;
628
1
    s->sg = s->gr + s->gg + s->gb + s->ga;
629
1
    s->sb = s->br + s->bg + s->bb + s->ba;
630
631
1
    if (fabs(s->sr) <= DBL_EPSILON)
632
        s->sr = 1.;
633
634
1
    if (fabs(s->sg) <= DBL_EPSILON)
635
        s->sg = 1.;
636
637
1
    if (fabs(s->sb) <= DBL_EPSILON)
638
        s->sb = 1.;
639
640
257
    for (i = 0; i < size; i++) {
641
256
        s->lut[R][R][i] = lrint(i * s->rr);
642
256
        s->lut[R][G][i] = lrint(i * s->rg);
643
256
        s->lut[R][B][i] = lrint(i * s->rb);
644
256
        s->lut[R][A][i] = lrint(i * s->ra);
645
646
256
        s->lut[G][R][i] = lrint(i * s->gr);
647
256
        s->lut[G][G][i] = lrint(i * s->gg);
648
256
        s->lut[G][B][i] = lrint(i * s->gb);
649
256
        s->lut[G][A][i] = lrint(i * s->ga);
650
651
256
        s->lut[B][R][i] = lrint(i * s->br);
652
256
        s->lut[B][G][i] = lrint(i * s->bg);
653
256
        s->lut[B][B][i] = lrint(i * s->bb);
654
256
        s->lut[B][A][i] = lrint(i * s->ba);
655
656
256
        s->lut[A][R][i] = lrint(i * s->ar);
657
256
        s->lut[A][G][i] = lrint(i * s->ag);
658
256
        s->lut[A][B][i] = lrint(i * s->ab);
659
256
        s->lut[A][A][i] = lrint(i * s->aa);
660
    }
661
662




1
    switch (outlink->format) {
663
1
    case AV_PIX_FMT_BGR24:
664
    case AV_PIX_FMT_RGB24:
665
1
        s->filter_slice[0] = filter_slice_rgb24;
666
1
        s->filter_slice[1] = filter_slice_rgb24_pl;
667
1
        break;
668
    case AV_PIX_FMT_0BGR:
669
    case AV_PIX_FMT_0RGB:
670
    case AV_PIX_FMT_BGR0:
671
    case AV_PIX_FMT_RGB0:
672
        s->filter_slice[0] = filter_slice_rgb0;
673
        s->filter_slice[1] = filter_slice_rgb0_pl;
674
        break;
675
    case AV_PIX_FMT_ABGR:
676
    case AV_PIX_FMT_ARGB:
677
    case AV_PIX_FMT_BGRA:
678
    case AV_PIX_FMT_RGBA:
679
        s->filter_slice[0] = filter_slice_rgba;
680
        s->filter_slice[1] = filter_slice_rgba_pl;
681
        break;
682
    case AV_PIX_FMT_BGR48:
683
    case AV_PIX_FMT_RGB48:
684
        s->filter_slice[0] = filter_slice_rgb48;
685
        s->filter_slice[1] = filter_slice_rgb48_pl;
686
        break;
687
    case AV_PIX_FMT_BGRA64:
688
    case AV_PIX_FMT_RGBA64:
689
        s->filter_slice[0] = filter_slice_rgba64;
690
        s->filter_slice[1] = filter_slice_rgba64_pl;
691
        break;
692
    case AV_PIX_FMT_GBRP:
693
        s->filter_slice[0] = filter_slice_gbrp;
694
        s->filter_slice[1] = filter_slice_gbrp_pl;
695
        break;
696
    case AV_PIX_FMT_GBRAP:
697
        s->filter_slice[0] = filter_slice_gbrap;
698
        s->filter_slice[1] = filter_slice_gbrap_pl;
699
        break;
700
    case AV_PIX_FMT_GBRP9:
701
        s->filter_slice[0] = filter_slice_gbrp9;
702
        s->filter_slice[1] = filter_slice_gbrp9_pl;
703
        break;
704
    case AV_PIX_FMT_GBRP10:
705
        s->filter_slice[0] = filter_slice_gbrp10;
706
        s->filter_slice[1] = filter_slice_gbrp10_pl;
707
        break;
708
    case AV_PIX_FMT_GBRAP10:
709
        s->filter_slice[0] = filter_slice_gbrap10;
710
        s->filter_slice[1] = filter_slice_gbrap10_pl;
711
        break;
712
    case AV_PIX_FMT_GBRP12:
713
        s->filter_slice[0] = filter_slice_gbrp12;
714
        s->filter_slice[1] = filter_slice_gbrp12_pl;
715
        break;
716
    case AV_PIX_FMT_GBRAP12:
717
        s->filter_slice[0] = filter_slice_gbrap12;
718
        s->filter_slice[1] = filter_slice_gbrap12_pl;
719
        break;
720
    case AV_PIX_FMT_GBRP14:
721
        s->filter_slice[0] = filter_slice_gbrp14;
722
        s->filter_slice[1] = filter_slice_gbrp14_pl;
723
        break;
724
    case AV_PIX_FMT_GBRP16:
725
        s->filter_slice[0] = filter_slice_gbrp16;
726
        s->filter_slice[1] = filter_slice_gbrp16_pl;
727
        break;
728
    case AV_PIX_FMT_GBRAP16:
729
        s->filter_slice[0] = filter_slice_gbrap16;
730
        s->filter_slice[1] = filter_slice_gbrap16_pl;
731
        break;
732
    }
733
734
1
    return 0;
735
}
736
737
50
static int filter_frame(AVFilterLink *inlink, AVFrame *in)
738
{
739
50
    AVFilterContext *ctx = inlink->dst;
740
50
    ColorChannelMixerContext *s = ctx->priv;
741
50
    AVFilterLink *outlink = ctx->outputs[0];
742
50
    const int pl = s->preserve_lightness > 0.;
743
    ThreadData td;
744
    AVFrame *out;
745
746
50
    if (av_frame_is_writable(in)) {
747
50
        out = in;
748
    } else {
749
        out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
750
        if (!out) {
751
            av_frame_free(&in);
752
            return AVERROR(ENOMEM);
753
        }
754
        av_frame_copy_props(out, in);
755
    }
756
757
50
    td.in = in;
758
50
    td.out = out;
759
50
    ctx->internal->execute(ctx, s->filter_slice[pl], &td, NULL, FFMIN(outlink->h, ff_filter_get_nb_threads(ctx)));
760
761
50
    if (in != out)
762
        av_frame_free(&in);
763
50
    return ff_filter_frame(outlink, out);
764
}
765
766
static int process_command(AVFilterContext *ctx, const char *cmd, const char *args,
767
                           char *res, int res_len, int flags)
768
{
769
    int ret = ff_filter_process_command(ctx, cmd, args, res, res_len, flags);
770
771
    if (ret < 0)
772
        return ret;
773
774
    return config_output(ctx->outputs[0]);
775
}
776
777
1
static av_cold void uninit(AVFilterContext *ctx)
778
{
779
1
    ColorChannelMixerContext *s = ctx->priv;
780
781
1
    av_freep(&s->buffer);
782
1
}
783
784
static const AVFilterPad colorchannelmixer_inputs[] = {
785
    {
786
        .name         = "default",
787
        .type         = AVMEDIA_TYPE_VIDEO,
788
        .filter_frame = filter_frame,
789
    },
790
    { NULL }
791
};
792
793
static const AVFilterPad colorchannelmixer_outputs[] = {
794
    {
795
        .name         = "default",
796
        .type         = AVMEDIA_TYPE_VIDEO,
797
        .config_props = config_output,
798
    },
799
    { NULL }
800
};
801
802
AVFilter ff_vf_colorchannelmixer = {
803
    .name          = "colorchannelmixer",
804
    .description   = NULL_IF_CONFIG_SMALL("Adjust colors by mixing color channels."),
805
    .priv_size     = sizeof(ColorChannelMixerContext),
806
    .priv_class    = &colorchannelmixer_class,
807
    .uninit        = uninit,
808
    .query_formats = query_formats,
809
    .inputs        = colorchannelmixer_inputs,
810
    .outputs       = colorchannelmixer_outputs,
811
    .flags         = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | AVFILTER_FLAG_SLICE_THREADS,
812
    .process_command = process_command,
813
};