FFmpeg coverage


Directory: ../../../ffmpeg/
File: src/libavfilter/vsrc_gradients.c
Date: 2024-02-16 17:37:06
Exec Total Coverage
Lines: 0 228 0.0%
Functions: 0 12 0.0%
Branches: 0 97 0.0%

Line Branch Exec Source
1 /*
2 * Copyright (c) 2020 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 "avfilter.h"
22 #include "filters.h"
23 #include "video.h"
24 #include "internal.h"
25 #include "libavutil/imgutils.h"
26 #include "libavutil/opt.h"
27 #include "libavutil/lfg.h"
28 #include "libavutil/random_seed.h"
29 #include <float.h>
30 #include <math.h>
31
32 typedef struct GradientsContext {
33 const AVClass *class;
34 int w, h;
35 int type;
36 AVRational frame_rate;
37 int64_t pts;
38 int64_t duration; ///< duration expressed in microseconds
39 float speed;
40 float angle;
41
42 uint8_t color_rgba[8][4];
43 float color_rgbaf[8][4];
44 int nb_colors;
45 int x0, y0, x1, y1;
46 float fx0, fy0, fx1, fy1;
47
48 int64_t seed;
49
50 AVLFG lfg;
51 int (*draw_slice)(AVFilterContext *ctx, void *arg, int job, int nb_jobs);
52 } GradientsContext;
53
54 #define OFFSET(x) offsetof(GradientsContext, x)
55 #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
56 #define VFT AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM
57
58 static const AVOption gradients_options[] = {
59 {"size", "set frame size", OFFSET(w), AV_OPT_TYPE_IMAGE_SIZE, {.str="640x480"}, 0, 0, FLAGS },
60 {"s", "set frame size", OFFSET(w), AV_OPT_TYPE_IMAGE_SIZE, {.str="640x480"}, 0, 0, FLAGS },
61 {"rate", "set frame rate", OFFSET(frame_rate), AV_OPT_TYPE_VIDEO_RATE, {.str="25"}, 0, INT_MAX, FLAGS },
62 {"r", "set frame rate", OFFSET(frame_rate), AV_OPT_TYPE_VIDEO_RATE, {.str="25"}, 0, INT_MAX, FLAGS },
63 {"c0", "set 1st color", OFFSET(color_rgba[0]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
64 {"c1", "set 2nd color", OFFSET(color_rgba[1]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
65 {"c2", "set 3rd color", OFFSET(color_rgba[2]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
66 {"c3", "set 4th color", OFFSET(color_rgba[3]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
67 {"c4", "set 5th color", OFFSET(color_rgba[4]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
68 {"c5", "set 6th color", OFFSET(color_rgba[5]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
69 {"c6", "set 7th color", OFFSET(color_rgba[6]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
70 {"c7", "set 8th color", OFFSET(color_rgba[7]), AV_OPT_TYPE_COLOR, {.str = "random"}, 0, 0, FLAGS },
71 {"x0", "set gradient line source x0", OFFSET(x0), AV_OPT_TYPE_INT, {.i64=-1}, -1, INT_MAX, FLAGS },
72 {"y0", "set gradient line source y0", OFFSET(y0), AV_OPT_TYPE_INT, {.i64=-1}, -1, INT_MAX, FLAGS },
73 {"x1", "set gradient line destination x1", OFFSET(x1), AV_OPT_TYPE_INT, {.i64=-1}, -1, INT_MAX, FLAGS },
74 {"y1", "set gradient line destination y1", OFFSET(y1), AV_OPT_TYPE_INT, {.i64=-1}, -1, INT_MAX, FLAGS },
75 {"nb_colors", "set the number of colors", OFFSET(nb_colors), AV_OPT_TYPE_INT, {.i64=2}, 2, 8, FLAGS },
76 {"n", "set the number of colors", OFFSET(nb_colors), AV_OPT_TYPE_INT, {.i64=2}, 2, 8, FLAGS },
77 {"seed", "set the seed", OFFSET(seed), AV_OPT_TYPE_INT64, {.i64=-1}, -1, UINT32_MAX, FLAGS },
78 {"duration", "set video duration", OFFSET(duration), AV_OPT_TYPE_DURATION, {.i64=-1}, -1, INT64_MAX, FLAGS },
79 {"d", "set video duration", OFFSET(duration), AV_OPT_TYPE_DURATION, {.i64=-1}, -1, INT64_MAX, FLAGS },
80 {"speed", "set gradients rotation speed", OFFSET(speed), AV_OPT_TYPE_FLOAT,{.dbl=0.01}, 0, 1, VFT },
81 {"type", "set gradient type", OFFSET(type), AV_OPT_TYPE_INT, {.i64=0}, 0, 4, VFT, .unit = "type" },
82 {"t", "set gradient type", OFFSET(type), AV_OPT_TYPE_INT, {.i64=0}, 0, 4, VFT, .unit = "type" },
83 { "linear", "set linear gradient", 0, AV_OPT_TYPE_CONST, {.i64=0}, 0, 0, VFT, .unit = "type" },
84 { "radial", "set radial gradient", 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, VFT, .unit = "type" },
85 { "circular", "set circular gradient", 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, VFT, .unit = "type" },
86 { "spiral", "set spiral gradient", 0, AV_OPT_TYPE_CONST, {.i64=3}, 0, 0, VFT, .unit = "type" },
87 { "square", "set square gradient", 0, AV_OPT_TYPE_CONST, {.i64=4}, 0, 0, VFT, .unit = "type" },
88 {NULL},
89 };
90
91 AVFILTER_DEFINE_CLASS(gradients);
92
93 static float lerpf(float a, float b, float x)
94 {
95 const float y = 1.f - x;
96
97 return a * y + b * x;
98 }
99
100 static uint32_t lerp_color(uint8_t c0[4], uint8_t c1[4], float x)
101 {
102 const float y = 1.f - x;
103
104 return (lrintf(c0[0] * y + c1[0] * x)) << 0 |
105 (lrintf(c0[1] * y + c1[1] * x)) << 8 |
106 (lrintf(c0[2] * y + c1[2] * x)) << 16 |
107 (lrintf(c0[3] * y + c1[3] * x)) << 24;
108 }
109
110 static uint64_t lerp_color16(uint8_t c0[4], uint8_t c1[4], float x)
111 {
112 const float y = 1.f - x;
113
114 return ((uint64_t)llrintf((c0[0] * y + c1[0] * x) * 256)) << 0 |
115 ((uint64_t)llrintf((c0[1] * y + c1[1] * x) * 256)) << 16 |
116 ((uint64_t)llrintf((c0[2] * y + c1[2] * x) * 256)) << 32 |
117 ((uint64_t)llrintf((c0[3] * y + c1[3] * x) * 256)) << 48;
118 }
119
120 static uint32_t lerp_colors(uint8_t arr[8][4], int nb_colors, int nb_wrap_colors, float step)
121 {
122 float scl;
123 int i, j;
124
125 if (nb_colors == 1 || step <= 0.0) {
126 return arr[0][0] | (arr[0][1] << 8) | (arr[0][2] << 16) | (arr[0][3] << 24);
127 } else if (step >= 1.0) {
128 i = nb_colors - 1;
129 return arr[i][0] | (arr[i][1] << 8) | (arr[i][2] << 16) | (arr[i][3] << 24);
130 }
131
132 scl = step * (nb_wrap_colors - 1);
133 i = floorf(scl);
134 j = i + 1;
135 if (i >= nb_colors - 1) {
136 i = nb_colors - 1;
137 j = 0;
138 }
139
140 return lerp_color(arr[i], arr[j], scl - i);
141 }
142
143 static uint64_t lerp_colors16(uint8_t arr[8][4], int nb_colors, int nb_wrap_colors, float step)
144 {
145 float scl;
146 int i, j;
147
148 if (nb_colors == 1 || step <= 0.0) {
149 return ((uint64_t)arr[0][0] << 8) | ((uint64_t)arr[0][1] << 24) | ((uint64_t)arr[0][2] << 40) | ((uint64_t)arr[0][3] << 56);
150 } else if (step >= 1.0) {
151 i = nb_colors - 1;
152 return ((uint64_t)arr[i][0] << 8) | ((uint64_t)arr[i][1] << 24) | ((uint64_t)arr[i][2] << 40) | ((uint64_t)arr[i][3] << 56);
153 }
154
155 scl = step * (nb_wrap_colors - 1);
156 i = floorf(scl);
157 j = i + 1;
158 if (i >= nb_colors - 1) {
159 i = nb_colors - 1;
160 j = 0;
161 }
162
163 return lerp_color16(arr[i], arr[j], scl - i);
164 }
165
166 static void lerp_colors32(float arr[8][4], int nb_colors,
167 int nb_wrap_colors, float step,
168 float *r, float *g, float *b, float *a)
169 {
170 float scl, x;
171 int i, j;
172
173 if (nb_colors == 1 || step <= 0.0) {
174 *r = arr[0][0];
175 *g = arr[0][1];
176 *b = arr[0][2];
177 *a = arr[0][3];
178 return;
179 } else if (step >= 1.0) {
180 i = nb_colors - 1;
181 *r = arr[i][0];
182 *g = arr[i][1];
183 *b = arr[i][2];
184 *a = arr[i][3];
185 return;
186 }
187
188 scl = step * (nb_wrap_colors - 1);
189 i = floorf(scl);
190 j = i + 1;
191 if (i >= nb_colors - 1) {
192 i = nb_colors - 1;
193 j = 0;
194 }
195 x = scl - i;
196
197 *r = lerpf(arr[i][0], arr[j][0], x);
198 *g = lerpf(arr[i][1], arr[j][1], x);
199 *b = lerpf(arr[i][2], arr[j][2], x);
200 *a = lerpf(arr[i][3], arr[j][3], x);
201 }
202
203 static float project(float origin_x, float origin_y,
204 float dest_x, float dest_y,
205 float point_x, float point_y, int type)
206 {
207 float op_x = point_x - origin_x;
208 float op_y = point_y - origin_y;
209 float od_x = dest_x - origin_x;
210 float od_y = dest_y - origin_y;
211 float op_x_od;
212 float od_s_q;
213
214 switch (type) {
215 case 0:
216 od_s_q = od_x * od_x + od_y * od_y;
217 break;
218 case 1:
219 od_s_q = sqrtf(od_x * od_x + od_y * od_y);
220 break;
221 case 2:
222 case 3:
223 od_s_q = M_PI * 2.f;
224 break;
225 case 4:
226 od_s_q = fmaxf(fabsf(od_x), fabsf(od_y));
227 break;
228 }
229
230 switch (type) {
231 case 0:
232 op_x_od = op_x * od_x + op_y * od_y;
233 break;
234 case 1:
235 op_x_od = sqrtf(op_x * op_x + op_y * op_y);
236 break;
237 case 2:
238 op_x_od = atan2f(op_x, op_y) + M_PI;
239 break;
240 case 3:
241 op_x_od = fmodf(atan2f(op_x, op_y) + M_PI + point_x / fmaxf(origin_x, dest_x), 2.f * M_PI);
242 break;
243 case 4:
244 op_x_od = fmaxf(fabsf(op_x), fabsf(op_y));
245 break;
246 }
247
248 // Normalize and clamp range.
249 return av_clipf(op_x_od / od_s_q, 0.f, 1.f);
250 }
251
252 static int draw_gradients_slice(AVFilterContext *ctx, void *arg, int job, int nb_jobs)
253 {
254 GradientsContext *s = ctx->priv;
255 AVFrame *frame = arg;
256 const int width = frame->width;
257 const int height = frame->height;
258 const int start = (height * job ) / nb_jobs;
259 const int end = (height * (job+1)) / nb_jobs;
260 const ptrdiff_t linesize = frame->linesize[0] / 4;
261 uint32_t *dst = (uint32_t *)frame->data[0] + start * linesize;
262 const int type = s->type;
263
264 for (int y = start; y < end; y++) {
265 for (int x = 0; x < width; x++) {
266 float factor = project(s->fx0, s->fy0, s->fx1, s->fy1, x, y, type);
267 dst[x] = lerp_colors(s->color_rgba, s->nb_colors, s->nb_colors + (type >= 2 && type <= 3), factor);
268 }
269
270 dst += linesize;
271 }
272
273 return 0;
274 }
275
276 static int draw_gradients_slice16(AVFilterContext *ctx, void *arg, int job, int nb_jobs)
277 {
278 GradientsContext *s = ctx->priv;
279 AVFrame *frame = arg;
280 const int width = frame->width;
281 const int height = frame->height;
282 const int start = (height * job ) / nb_jobs;
283 const int end = (height * (job+1)) / nb_jobs;
284 const ptrdiff_t linesize = frame->linesize[0] / 8;
285 uint64_t *dst = (uint64_t *)frame->data[0] + start * linesize;
286 const int type = s->type;
287
288 for (int y = start; y < end; y++) {
289 for (int x = 0; x < width; x++) {
290 float factor = project(s->fx0, s->fy0, s->fx1, s->fy1, x, y, type);
291 dst[x] = lerp_colors16(s->color_rgba, s->nb_colors, s->nb_colors + (type >= 2 && type <= 3), factor);
292 }
293
294 dst += linesize;
295 }
296
297 return 0;
298 }
299
300 static int draw_gradients_slice32_planar(AVFilterContext *ctx, void *arg, int job, int nb_jobs)
301 {
302 GradientsContext *s = ctx->priv;
303 AVFrame *frame = arg;
304 const int width = frame->width;
305 const int height = frame->height;
306 const int start = (height * job ) / nb_jobs;
307 const int end = (height * (job+1)) / nb_jobs;
308 const ptrdiff_t linesize_g = frame->linesize[0] / 4;
309 const ptrdiff_t linesize_b = frame->linesize[1] / 4;
310 const ptrdiff_t linesize_r = frame->linesize[2] / 4;
311 const ptrdiff_t linesize_a = frame->linesize[3] / 4;
312 float *dst_g = (float *)frame->data[0] + start * linesize_g;
313 float *dst_b = (float *)frame->data[1] + start * linesize_b;
314 float *dst_r = (float *)frame->data[2] + start * linesize_r;
315 float *dst_a = (float *)frame->data[3] + start * linesize_a;
316 const int type = s->type;
317
318 for (int y = start; y < end; y++) {
319 for (int x = 0; x < width; x++) {
320 float factor = project(s->fx0, s->fy0, s->fx1, s->fy1, x, y, type);
321 lerp_colors32(s->color_rgbaf, s->nb_colors, s->nb_colors + (type >= 2 && type <= 3), factor,
322 &dst_r[x], &dst_g[x], &dst_b[x], &dst_a[x]);
323 }
324
325 dst_g += linesize_g;
326 dst_b += linesize_b;
327 dst_r += linesize_r;
328 dst_a += linesize_a;
329 }
330
331 return 0;
332 }
333
334 static int config_output(AVFilterLink *outlink)
335 {
336 AVFilterContext *ctx = outlink->src;
337 GradientsContext *s = ctx->priv;
338 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(outlink->format);
339
340 if (av_image_check_size(s->w, s->h, 0, ctx) < 0)
341 return AVERROR(EINVAL);
342
343 outlink->w = s->w;
344 outlink->h = s->h;
345 outlink->time_base = av_inv_q(s->frame_rate);
346 outlink->sample_aspect_ratio = (AVRational) {1, 1};
347 outlink->frame_rate = s->frame_rate;
348 if (s->seed == -1)
349 s->seed = av_get_random_seed();
350 av_lfg_init(&s->lfg, s->seed);
351
352 switch (desc->comp[0].depth) {
353 case 8:
354 s->draw_slice = draw_gradients_slice;
355 break;
356 case 16:
357 s->draw_slice = draw_gradients_slice16;
358 break;
359 case 32:
360 s->draw_slice = draw_gradients_slice32_planar;
361 break;
362 default:
363 return AVERROR_BUG;
364 }
365
366 if (s->x0 < 0 || s->x0 >= s->w)
367 s->x0 = av_lfg_get(&s->lfg) % s->w;
368 if (s->y0 < 0 || s->y0 >= s->h)
369 s->y0 = av_lfg_get(&s->lfg) % s->h;
370 if (s->x1 < 0 || s->x1 >= s->w)
371 s->x1 = av_lfg_get(&s->lfg) % s->w;
372 if (s->y1 < 0 || s->y1 >= s->h)
373 s->y1 = av_lfg_get(&s->lfg) % s->h;
374
375 for (int n = 0; n < 8; n++) {
376 for (int c = 0; c < 4; c++)
377 s->color_rgbaf[n][c] = s->color_rgba[n][c] / 255.f;
378 }
379
380 return 0;
381 }
382
383 static int activate(AVFilterContext *ctx)
384 {
385 GradientsContext *s = ctx->priv;
386 AVFilterLink *outlink = ctx->outputs[0];
387
388 if (s->duration >= 0 &&
389 av_rescale_q(s->pts, outlink->time_base, AV_TIME_BASE_Q) >= s->duration) {
390 ff_outlink_set_status(outlink, AVERROR_EOF, s->pts);
391 return 0;
392 }
393
394 if (ff_outlink_frame_wanted(outlink)) {
395 AVFrame *frame = ff_get_video_buffer(outlink, s->w, s->h);
396 float angle = fmodf(s->angle, 2.f * M_PI);
397 const float w2 = s->w / 2.f;
398 const float h2 = s->h / 2.f;
399
400 s->angle = angle + s->speed;
401
402 s->fx0 = (s->x0 - w2) * cosf(angle) - (s->y0 - h2) * sinf(angle) + w2;
403 s->fy0 = (s->x0 - w2) * sinf(angle) + (s->y0 - h2) * cosf(angle) + h2;
404
405 s->fx1 = (s->x1 - w2) * cosf(angle) - (s->y1 - h2) * sinf(angle) + w2;
406 s->fy1 = (s->x1 - w2) * sinf(angle) + (s->y1 - h2) * cosf(angle) + h2;
407
408 if (!frame)
409 return AVERROR(ENOMEM);
410
411 #if FF_API_FRAME_KEY
412 FF_DISABLE_DEPRECATION_WARNINGS
413 frame->key_frame = 1;
414 FF_ENABLE_DEPRECATION_WARNINGS
415 #endif
416
417 frame->flags |= AV_FRAME_FLAG_KEY;
418 #if FF_API_INTERLACED_FRAME
419 FF_DISABLE_DEPRECATION_WARNINGS
420 frame->interlaced_frame = 0;
421 FF_ENABLE_DEPRECATION_WARNINGS
422 #endif
423 frame->flags &= ~AV_FRAME_FLAG_INTERLACED;
424 frame->pict_type = AV_PICTURE_TYPE_I;
425 frame->sample_aspect_ratio = (AVRational) {1, 1};
426 frame->pts = s->pts++;
427 frame->duration = 1;
428
429 ff_filter_execute(ctx, s->draw_slice, frame, NULL,
430 FFMIN(outlink->h, ff_filter_get_nb_threads(ctx)));
431
432 return ff_filter_frame(outlink, frame);
433 }
434
435 return FFERROR_NOT_READY;
436 }
437
438 static const AVFilterPad gradients_outputs[] = {
439 {
440 .name = "default",
441 .type = AVMEDIA_TYPE_VIDEO,
442 .config_props = config_output,
443 },
444 };
445
446 const AVFilter ff_vsrc_gradients = {
447 .name = "gradients",
448 .description = NULL_IF_CONFIG_SMALL("Draw a gradients."),
449 .priv_size = sizeof(GradientsContext),
450 .priv_class = &gradients_class,
451 .inputs = NULL,
452 FILTER_OUTPUTS(gradients_outputs),
453 FILTER_PIXFMTS(AV_PIX_FMT_RGBA, AV_PIX_FMT_RGBA64, AV_PIX_FMT_GBRAPF32),
454 .activate = activate,
455 .flags = AVFILTER_FLAG_SLICE_THREADS,
456 .process_command = ff_filter_process_command,
457 };
458