FFmpeg
vf_fftfilt.c
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2015 Arwa Arif <arwaarif1994@gmail.com>
3  * Copyright (c) 2017 Paul B Mahol
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg is free software; you can redistribute it and/or modify it
8  * under the terms of the GNU Lesser General Public License as published
9  * by the Free Software Foundation; either version 2.1 of the License,
10  * or (at your option) any later version.
11  *
12  * FFmpeg is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21 
22 /**
23  * @file
24  * FFT domain filtering.
25  */
26 
27 #include "libavfilter/internal.h"
28 #include "libavutil/common.h"
29 #include "libavutil/cpu.h"
30 #include "libavutil/imgutils.h"
31 #include "libavutil/opt.h"
32 #include "libavutil/pixdesc.h"
33 #include "libavutil/tx.h"
34 #include "libavutil/eval.h"
35 
36 #define MAX_THREADS 32
37 #define MAX_PLANES 4
38 
39 enum EvalMode {
40  EVAL_MODE_INIT,
41  EVAL_MODE_FRAME,
42  EVAL_MODE_NB
43 };
44 
45 typedef struct FFTFILTContext {
46  const AVClass *class;
47 
48  int eval_mode;
49  int depth;
50  int nb_planes;
51  int nb_threads;
52  int planewidth[MAX_PLANES];
53  int planeheight[MAX_PLANES];
54 
55  AVTXContext *hrdft[MAX_THREADS][MAX_PLANES];
56  AVTXContext *vrdft[MAX_THREADS][MAX_PLANES];
57  AVTXContext *ihrdft[MAX_THREADS][MAX_PLANES];
58  AVTXContext *ivrdft[MAX_THREADS][MAX_PLANES];
59 
60  av_tx_fn htx_fn, ihtx_fn;
61  av_tx_fn vtx_fn, ivtx_fn;
62 
63  int rdft_hbits[MAX_PLANES];
64  int rdft_vbits[MAX_PLANES];
65  size_t rdft_hstride[MAX_PLANES];
66  size_t rdft_vstride[MAX_PLANES];
67  size_t rdft_hlen[MAX_PLANES];
68  size_t rdft_vlen[MAX_PLANES];
69  float *rdft_hdata_in[MAX_PLANES];
70  float *rdft_vdata_in[MAX_PLANES];
71  float *rdft_hdata_out[MAX_PLANES];
72  float *rdft_vdata_out[MAX_PLANES];
73 
74  int dc[MAX_PLANES];
75  char *weight_str[MAX_PLANES];
76  AVExpr *weight_expr[MAX_PLANES];
77  double *weight[MAX_PLANES];
78 
79  int (*rdft_horizontal)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs);
80  int (*irdft_horizontal)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs);
81 } FFTFILTContext;
82 
83 static const char *const var_names[] = { "X", "Y", "W", "H", "N", "WS", "HS", NULL };
84 enum { VAR_X, VAR_Y, VAR_W, VAR_H, VAR_N, VAR_WS, VAR_HS, VAR_VARS_NB };
85 
86 enum { Y = 0, U, V };
87 
88 #define OFFSET(x) offsetof(FFTFILTContext, x)
89 #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
90 
91 static const AVOption fftfilt_options[] = {
92  { "dc_Y", "adjust gain in Y plane", OFFSET(dc[Y]), AV_OPT_TYPE_INT, {.i64 = 0}, 0, 1000, FLAGS },
93  { "dc_U", "adjust gain in U plane", OFFSET(dc[U]), AV_OPT_TYPE_INT, {.i64 = 0}, 0, 1000, FLAGS },
94  { "dc_V", "adjust gain in V plane", OFFSET(dc[V]), AV_OPT_TYPE_INT, {.i64 = 0}, 0, 1000, FLAGS },
95  { "weight_Y", "set luminance expression in Y plane", OFFSET(weight_str[Y]), AV_OPT_TYPE_STRING, {.str = "1"}, 0, 0, FLAGS },
96  { "weight_U", "set chrominance expression in U plane", OFFSET(weight_str[U]), AV_OPT_TYPE_STRING, {.str = NULL}, 0, 0, FLAGS },
97  { "weight_V", "set chrominance expression in V plane", OFFSET(weight_str[V]), AV_OPT_TYPE_STRING, {.str = NULL}, 0, 0, FLAGS },
98  { "eval", "specify when to evaluate expressions", OFFSET(eval_mode), AV_OPT_TYPE_INT, {.i64 = EVAL_MODE_INIT}, 0, EVAL_MODE_NB-1, FLAGS, "eval" },
99  { "init", "eval expressions once during initialization", 0, AV_OPT_TYPE_CONST, {.i64=EVAL_MODE_INIT}, .flags = FLAGS, .unit = "eval" },
100  { "frame", "eval expressions per-frame", 0, AV_OPT_TYPE_CONST, {.i64=EVAL_MODE_FRAME}, .flags = FLAGS, .unit = "eval" },
101  {NULL},
102 };
103 
104 AVFILTER_DEFINE_CLASS(fftfilt);
105 
106 static inline double lum(void *priv, double x, double y, int plane)
107 {
108  FFTFILTContext *s = priv;
109  return s->rdft_vdata_out[plane][(int)x * s->rdft_vstride[plane] + (int)y];
110 }
111 
112 static double weight_Y(void *priv, double x, double y) { return lum(priv, x, y, Y); }
113 static double weight_U(void *priv, double x, double y) { return lum(priv, x, y, U); }
114 static double weight_V(void *priv, double x, double y) { return lum(priv, x, y, V); }
115 
116 static void copy_rev(float *dest, int w, int w2)
117 {
118  int i;
119 
120  for (i = w; i < w + (w2-w)/2; i++)
121  dest[i] = dest[2*w - i - 1];
122 
123  for (; i < w2; i++)
124  dest[i] = dest[w2 - i];
125 }
126 
127 static int rdft_horizontal8(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
128 {
129  FFTFILTContext *s = ctx->priv;
130  AVFrame *in = arg;
131 
132  for (int plane = 0; plane < s->nb_planes; plane++) {
133  const int w = s->planewidth[plane];
134  const int h = s->planeheight[plane];
135  const int slice_start = (h * jobnr) / nb_jobs;
136  const int slice_end = (h * (jobnr+1)) / nb_jobs;
137 
138  for (int i = slice_start; i < slice_end; i++) {
139  const uint8_t *src = in->data[plane] + i * in->linesize[plane];
140  float *hdata_in = s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane];
141 
142  for (int j = 0; j < w; j++)
143  hdata_in[j] = src[j];
144 
145  copy_rev(s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane], w, s->rdft_hlen[plane]);
146  }
147 
148  for (int i = slice_start; i < slice_end; i++)
149  s->htx_fn(s->hrdft[jobnr][plane],
150  s->rdft_hdata_out[plane] + i * s->rdft_hstride[plane],
151  s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane],
152  sizeof(float));
153  }
154 
155  return 0;
156 }
157 
158 static int rdft_horizontal16(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
159 {
160  FFTFILTContext *s = ctx->priv;
161  AVFrame *in = arg;
162 
163  for (int plane = 0; plane < s->nb_planes; plane++) {
164  const int w = s->planewidth[plane];
165  const int h = s->planeheight[plane];
166  const int slice_start = (h * jobnr) / nb_jobs;
167  const int slice_end = (h * (jobnr+1)) / nb_jobs;
168 
169  for (int i = slice_start; i < slice_end; i++) {
170  const uint16_t *src = (const uint16_t *)(in->data[plane] + i * in->linesize[plane]);
171  float *hdata_in = s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane];
172 
173  for (int j = 0; j < w; j++)
174  hdata_in[j] = src[j];
175 
176  copy_rev(s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane], w, s->rdft_hlen[plane]);
177  }
178 
179  for (int i = slice_start; i < slice_end; i++)
180  s->htx_fn(s->hrdft[jobnr][plane],
181  s->rdft_hdata_out[plane] + i * s->rdft_hstride[plane],
182  s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane],
183  sizeof(float));
184  }
185 
186  return 0;
187 }
188 
189 static int irdft_horizontal8(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
190 {
191  FFTFILTContext *s = ctx->priv;
192  AVFrame *out = arg;
193 
194  for (int plane = 0; plane < s->nb_planes; plane++) {
195  const int w = s->planewidth[plane];
196  const int h = s->planeheight[plane];
197  const int slice_start = (h * jobnr) / nb_jobs;
198  const int slice_end = (h * (jobnr+1)) / nb_jobs;
199 
200  for (int i = slice_start; i < slice_end; i++)
201  s->ihtx_fn(s->ihrdft[jobnr][plane],
202  s->rdft_hdata_out[plane] + i * s->rdft_hstride[plane],
203  s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane],
204  sizeof(float));
205 
206  for (int i = slice_start; i < slice_end; i++) {
207  const float scale = 1.f / (s->rdft_hlen[plane] * s->rdft_vlen[plane]);
208  const float *src = s->rdft_hdata_out[plane] + i * s->rdft_hstride[plane];
209  uint8_t *dst = out->data[plane] + i * out->linesize[plane];
210 
211  for (int j = 0; j < w; j++)
212  dst[j] = av_clip_uint8(lrintf(src[j] * scale));
213  }
214  }
215 
216  return 0;
217 }
218 
219 static int irdft_horizontal16(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
220 {
221  FFTFILTContext *s = ctx->priv;
222  AVFrame *out = arg;
223 
224  for (int plane = 0; plane < s->nb_planes; plane++) {
225  int max = (1 << s->depth) - 1;
226  const int w = s->planewidth[plane];
227  const int h = s->planeheight[plane];
228  const int slice_start = (h * jobnr) / nb_jobs;
229  const int slice_end = (h * (jobnr+1)) / nb_jobs;
230 
231  for (int i = slice_start; i < slice_end; i++)
232  s->ihtx_fn(s->ihrdft[jobnr][plane],
233  s->rdft_hdata_out[plane] + i * s->rdft_hstride[plane],
234  s->rdft_hdata_in[plane] + i * s->rdft_hstride[plane],
235  sizeof(float));
236 
237  for (int i = slice_start; i < slice_end; i++) {
238  const float scale = 1.f / (s->rdft_hlen[plane] * s->rdft_vlen[plane]);
239  const float *src = s->rdft_hdata_out[plane] + i * s->rdft_hstride[plane];
240  uint16_t *dst = (uint16_t *)(out->data[plane] + i * out->linesize[plane]);
241 
242  for (int j = 0; j < w; j++)
243  dst[j] = av_clip(lrintf(src[j] * scale), 0, max);
244  }
245  }
246 
247  return 0;
248 }
249 
250 static av_cold int initialize(AVFilterContext *ctx)
251 {
252  FFTFILTContext *s = ctx->priv;
253  int ret = 0, plane;
254 
255  if (!s->dc[U] && !s->dc[V]) {
256  s->dc[U] = s->dc[Y];
257  s->dc[V] = s->dc[Y];
258  } else {
259  if (!s->dc[U]) s->dc[U] = s->dc[V];
260  if (!s->dc[V]) s->dc[V] = s->dc[U];
261  }
262 
263  if (!s->weight_str[U] && !s->weight_str[V]) {
264  s->weight_str[U] = av_strdup(s->weight_str[Y]);
265  s->weight_str[V] = av_strdup(s->weight_str[Y]);
266  } else {
267  if (!s->weight_str[U]) s->weight_str[U] = av_strdup(s->weight_str[V]);
268  if (!s->weight_str[V]) s->weight_str[V] = av_strdup(s->weight_str[U]);
269  }
270 
271  for (plane = 0; plane < 3; plane++) {
272  static double (*p[])(void *, double, double) = { weight_Y, weight_U, weight_V };
273  const char *const func2_names[] = {"weight_Y", "weight_U", "weight_V", NULL };
274  double (*func2[])(void *, double, double) = { weight_Y, weight_U, weight_V, p[plane], NULL };
275 
276  ret = av_expr_parse(&s->weight_expr[plane], s->weight_str[plane], var_names,
277  NULL, NULL, func2_names, func2, 0, ctx);
278  if (ret < 0)
279  break;
280  }
281  return ret;
282 }
283 
284 static void do_eval(FFTFILTContext *s, AVFilterLink *inlink, int plane)
285 {
286  double values[VAR_VARS_NB];
287  int i, j;
288 
289  values[VAR_N] = inlink->frame_count_out;
290  values[VAR_W] = s->planewidth[plane];
291  values[VAR_H] = s->planeheight[plane];
292  values[VAR_WS] = s->rdft_hlen[plane];
293  values[VAR_HS] = s->rdft_vlen[plane];
294 
295  for (i = 0; i < s->rdft_hlen[plane]; i++) {
296  values[VAR_X] = i;
297  for (j = 0; j < s->rdft_vlen[plane]; j++) {
298  values[VAR_Y] = j;
299  s->weight[plane][i * s->rdft_vlen[plane] + j] =
300  av_expr_eval(s->weight_expr[plane], values, s);
301  }
302  }
303 }
304 
305 static int config_props(AVFilterLink *inlink)
306 {
307  FFTFILTContext *s = inlink->dst->priv;
308  const AVPixFmtDescriptor *desc;
309  int i, plane;
310 
311  desc = av_pix_fmt_desc_get(inlink->format);
312  s->depth = desc->comp[0].depth;
313  s->planewidth[1] = s->planewidth[2] = AV_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w);
314  s->planewidth[0] = s->planewidth[3] = inlink->w;
315  s->planeheight[1] = s->planeheight[2] = AV_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h);
316  s->planeheight[0] = s->planeheight[3] = inlink->h;
317 
318  s->nb_planes = av_pix_fmt_count_planes(inlink->format);
319  s->nb_threads = FFMIN(32, ff_filter_get_nb_threads(inlink->dst));
320 
321  for (i = 0; i < desc->nb_components; i++) {
322  int w = s->planewidth[i];
323  int h = s->planeheight[i];
324 
325  /* RDFT - Array initialization for Horizontal pass*/
326  s->rdft_hlen[i] = 1 << (32 - ff_clz(w));
327  s->rdft_hstride[i] = FFALIGN(s->rdft_hlen[i] + 2, av_cpu_max_align());
328  s->rdft_hbits[i] = av_log2(s->rdft_hlen[i]);
329  if (!(s->rdft_hdata_in[i] = av_calloc(h, s->rdft_hstride[i] * sizeof(float))))
330  return AVERROR(ENOMEM);
331 
332  if (!(s->rdft_hdata_out[i] = av_calloc(h, s->rdft_hstride[i] * sizeof(float))))
333  return AVERROR(ENOMEM);
334 
335  for (int j = 0; j < s->nb_threads; j++) {
336  float scale = 1.f, iscale = 1.f;
337 
338  av_tx_init(&s->hrdft[j][i], &s->htx_fn, AV_TX_FLOAT_RDFT, 0, 1 << s->rdft_hbits[i], &scale, 0);
339  if (!s->hrdft[j][i])
340  return AVERROR(ENOMEM);
341  av_tx_init(&s->ihrdft[j][i], &s->ihtx_fn, AV_TX_FLOAT_RDFT, 1, 1 << s->rdft_hbits[i], &iscale, 0);
342  if (!s->ihrdft[j][i])
343  return AVERROR(ENOMEM);
344  }
345 
346  /* RDFT - Array initialization for Vertical pass*/
347  s->rdft_vlen[i] = 1 << (32 - ff_clz(h));
348  s->rdft_vstride[i] = FFALIGN(s->rdft_vlen[i] + 2, av_cpu_max_align());
349  s->rdft_vbits[i] = av_log2(s->rdft_vlen[i]);
350  if (!(s->rdft_vdata_in[i] = av_calloc(s->rdft_hstride[i], s->rdft_vstride[i] * sizeof(float))))
351  return AVERROR(ENOMEM);
352 
353  if (!(s->rdft_vdata_out[i] = av_calloc(s->rdft_hstride[i], s->rdft_vstride[i] * sizeof(float))))
354  return AVERROR(ENOMEM);
355 
356  for (int j = 0; j < s->nb_threads; j++) {
357  float scale = 1.f, iscale = 1.f;
358 
359  av_tx_init(&s->vrdft[j][i], &s->vtx_fn, AV_TX_FLOAT_RDFT, 0, 1 << s->rdft_vbits[i], &scale, 0);
360  if (!s->vrdft[j][i])
361  return AVERROR(ENOMEM);
362  av_tx_init(&s->ivrdft[j][i], &s->ivtx_fn, AV_TX_FLOAT_RDFT, 1, 1 << s->rdft_vbits[i], &iscale, 0);
363  if (!s->ivrdft[j][i])
364  return AVERROR(ENOMEM);
365  }
366  }
367 
368  /*Luminance value - Array initialization*/
369  for (plane = 0; plane < 3; plane++) {
370  if(!(s->weight[plane] = av_calloc(s->rdft_hlen[plane], s->rdft_vlen[plane] * sizeof(double))))
371  return AVERROR(ENOMEM);
372 
373  if (s->eval_mode == EVAL_MODE_INIT)
374  do_eval(s, inlink, plane);
375  }
376 
377  if (s->depth <= 8) {
378  s->rdft_horizontal = rdft_horizontal8;
379  s->irdft_horizontal = irdft_horizontal8;
380  } else if (s->depth > 8) {
381  s->rdft_horizontal = rdft_horizontal16;
382  s->irdft_horizontal = irdft_horizontal16;
383  } else {
384  return AVERROR_BUG;
385  }
386  return 0;
387 }
388 
389 static int multiply_data(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
390 {
391  FFTFILTContext *s = ctx->priv;
392 
393  for (int plane = 0; plane < s->nb_planes; plane++) {
394  const int height = s->rdft_hlen[plane];
395  const int slice_start = (height * jobnr) / nb_jobs;
396  const int slice_end = (height * (jobnr+1)) / nb_jobs;
397  /*Change user defined parameters*/
398  for (int i = slice_start; i < slice_end; i++) {
399  const double *weight = s->weight[plane] + i * s->rdft_vlen[plane];
400  float *vdata = s->rdft_vdata_out[plane] + i * s->rdft_vstride[plane];
401 
402  for (int j = 0; j < s->rdft_vlen[plane]; j++)
403  vdata[j] *= weight[j];
404  }
405  }
406 
407  return 0;
408 }
409 
410 static int copy_vertical(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
411 {
412  FFTFILTContext *s = ctx->priv;
413 
414  for (int plane = 0; plane < s->nb_planes; plane++) {
415  const int hlen = s->rdft_hlen[plane];
416  const int vlen = s->rdft_vlen[plane];
417  const int hstride = s->rdft_hstride[plane];
418  const int vstride = s->rdft_vstride[plane];
419  const int slice_start = (hlen * jobnr) / nb_jobs;
420  const int slice_end = (hlen * (jobnr+1)) / nb_jobs;
421  const int h = s->planeheight[plane];
422  float *hdata = s->rdft_hdata_out[plane];
423  float *vdata = s->rdft_vdata_in[plane];
424 
425  for (int i = slice_start; i < slice_end; i++) {
426  for (int j = 0; j < h; j++)
427  vdata[i * vstride + j] = hdata[j * hstride + i];
428  copy_rev(vdata + i * vstride, h, vlen);
429  }
430  }
431 
432  return 0;
433 }
434 
435 static int rdft_vertical(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
436 {
437  FFTFILTContext *s = ctx->priv;
438 
439  for (int plane = 0; plane < s->nb_planes; plane++) {
440  const int height = s->rdft_hlen[plane];
441  const int slice_start = (height * jobnr) / nb_jobs;
442  const int slice_end = (height * (jobnr+1)) / nb_jobs;
443 
444  for (int i = slice_start; i < slice_end; i++)
445  s->vtx_fn(s->vrdft[jobnr][plane],
446  s->rdft_vdata_out[plane] + i * s->rdft_vstride[plane],
447  s->rdft_vdata_in[plane] + i * s->rdft_vstride[plane],
448  sizeof(float));
449  }
450 
451  return 0;
452 }
453 
454 static int irdft_vertical(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
455 {
456  FFTFILTContext *s = ctx->priv;
457 
458  for (int plane = 0; plane < s->nb_planes; plane++) {
459  const int height = s->rdft_hlen[plane];
460  const int slice_start = (height * jobnr) / nb_jobs;
461  const int slice_end = (height * (jobnr+1)) / nb_jobs;
462 
463  for (int i = slice_start; i < slice_end; i++)
464  s->ivtx_fn(s->ivrdft[jobnr][plane],
465  s->rdft_vdata_in[plane] + i * s->rdft_vstride[plane],
466  s->rdft_vdata_out[plane] + i * s->rdft_vstride[plane],
467  sizeof(float));
468  }
469 
470  return 0;
471 }
472 
473 static int copy_horizontal(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
474 {
475  FFTFILTContext *s = ctx->priv;
476 
477  for (int plane = 0; plane < s->nb_planes; plane++) {
478  const int hlen = s->rdft_hlen[plane];
479  const int hstride = s->rdft_hstride[plane];
480  const int vstride = s->rdft_vstride[plane];
481  const int slice_start = (hlen * jobnr) / nb_jobs;
482  const int slice_end = (hlen * (jobnr+1)) / nb_jobs;
483  const int h = s->planeheight[plane];
484  float *hdata = s->rdft_hdata_in[plane];
485  float *vdata = s->rdft_vdata_in[plane];
486 
487  for (int i = slice_start; i < slice_end; i++)
488  for (int j = 0; j < h; j++)
489  hdata[j * hstride + i] = vdata[i * vstride + j];
490  }
491 
492  return 0;
493 }
494 
495 static int filter_frame(AVFilterLink *inlink, AVFrame *in)
496 {
497  AVFilterContext *ctx = inlink->dst;
498  AVFilterLink *outlink = inlink->dst->outputs[0];
499  FFTFILTContext *s = ctx->priv;
500  AVFrame *out;
501 
502  out = ff_get_video_buffer(outlink, inlink->w, inlink->h);
503  if (!out) {
504  av_frame_free(&in);
505  return AVERROR(ENOMEM);
506  }
507 
508  av_frame_copy_props(out, in);
509 
510  ff_filter_execute(ctx, s->rdft_horizontal, in, NULL,
511  FFMIN(s->planeheight[1], s->nb_threads));
512 
513  ff_filter_execute(ctx, copy_vertical, NULL, NULL,
514  FFMIN(s->planeheight[1], s->nb_threads));
515 
516  ff_filter_execute(ctx, rdft_vertical, NULL, NULL,
517  FFMIN(s->planeheight[1], s->nb_threads));
518 
519  for (int plane = 0; plane < s->nb_planes; plane++) {
520  if (s->eval_mode == EVAL_MODE_FRAME)
521  do_eval(s, inlink, plane);
522  }
523 
524  ff_filter_execute(ctx, multiply_data, NULL, NULL,
525  FFMIN(s->planeheight[1], s->nb_threads));
526 
527  for (int plane = 0; plane < s->nb_planes; plane++)
528  s->rdft_vdata_out[plane][0] += s->rdft_hlen[plane] * s->rdft_vlen[plane] * s->dc[plane] * (1 << (s->depth - 8));
529 
530  ff_filter_execute(ctx, irdft_vertical, NULL, NULL,
531  FFMIN(s->planeheight[1], s->nb_threads));
532 
533  ff_filter_execute(ctx, copy_horizontal, NULL, NULL,
534  FFMIN(s->planeheight[1], s->nb_threads));
535 
536  ff_filter_execute(ctx, s->irdft_horizontal, out, NULL,
537  FFMIN(s->planeheight[1], s->nb_threads));
538 
539  av_frame_free(&in);
540  return ff_filter_frame(outlink, out);
541 }
542 
543 static av_cold void uninit(AVFilterContext *ctx)
544 {
545  FFTFILTContext *s = ctx->priv;
546 
547  for (int i = 0; i < MAX_PLANES; i++) {
548  av_freep(&s->rdft_hdata_in[i]);
549  av_freep(&s->rdft_vdata_in[i]);
550  av_freep(&s->rdft_hdata_out[i]);
551  av_freep(&s->rdft_vdata_out[i]);
552  av_expr_free(s->weight_expr[i]);
553  av_freep(&s->weight[i]);
554  for (int j = 0; j < s->nb_threads; j++) {
555  av_tx_uninit(&s->hrdft[j][i]);
556  av_tx_uninit(&s->ihrdft[j][i]);
557  av_tx_uninit(&s->vrdft[j][i]);
558  av_tx_uninit(&s->ivrdft[j][i]);
559  }
560  }
561 }
562 
563 static const enum AVPixelFormat pixel_fmts_fftfilt[] = {
564  AV_PIX_FMT_GRAY8,
565  AV_PIX_FMT_GRAY9, AV_PIX_FMT_GRAY10, AV_PIX_FMT_GRAY12,
566  AV_PIX_FMT_GRAY14, AV_PIX_FMT_GRAY16,
567  AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUVJ444P,
568  AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUVJ420P,
569  AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUVJ422P,
570  AV_PIX_FMT_YUV420P9, AV_PIX_FMT_YUV420P10,
571  AV_PIX_FMT_YUV420P12, AV_PIX_FMT_YUV420P14,
572  AV_PIX_FMT_YUV420P16,
573  AV_PIX_FMT_YUV422P9, AV_PIX_FMT_YUV422P10,
574  AV_PIX_FMT_YUV422P12, AV_PIX_FMT_YUV422P14,
575  AV_PIX_FMT_YUV422P16,
576  AV_PIX_FMT_YUV444P9, AV_PIX_FMT_YUV444P10,
577  AV_PIX_FMT_YUV444P12, AV_PIX_FMT_YUV444P14,
578  AV_PIX_FMT_YUV444P16,
579  AV_PIX_FMT_NONE
580 };
581 
582 static const AVFilterPad fftfilt_inputs[] = {
583  {
584  .name = "default",
585  .type = AVMEDIA_TYPE_VIDEO,
586  .config_props = config_props,
587  .filter_frame = filter_frame,
588  },
589 };
590 
591 static const AVFilterPad fftfilt_outputs[] = {
592  {
593  .name = "default",
594  .type = AVMEDIA_TYPE_VIDEO,
595  },
596 };
597 
598 const AVFilter ff_vf_fftfilt = {
599  .name = "fftfilt",
600  .description = NULL_IF_CONFIG_SMALL("Apply arbitrary expressions to pixels in frequency domain."),
601  .priv_size = sizeof(FFTFILTContext),
602  .priv_class = &fftfilt_class,
603  FILTER_INPUTS(fftfilt_inputs),
604  FILTER_OUTPUTS(fftfilt_outputs),
605  FILTER_PIXFMTS_ARRAY(pixel_fmts_fftfilt),
606  .init = initialize,
607  .uninit = uninit,
608  .flags = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | AVFILTER_FLAG_SLICE_THREADS,
609 };
ff_get_video_buffer
AVFrame * ff_get_video_buffer(AVFilterLink *link, int w, int h)
Request a picture buffer with a specific set of permissions.
Definition: video.c:101
FFTFILTContext::hrdft
AVTXContext * hrdft[MAX_THREADS][MAX_PLANES]
Definition: vf_fftfilt.c:55
AVPixelFormat
AVPixelFormat
Pixel format.
Definition: pixfmt.h:64
av_clip
#define av_clip
Definition: common.h:95
AVERROR
Filter the word “frame” indicates either a video frame or a group of audio as stored in an AVFrame structure Format for each input and each output the list of supported formats For video that means pixel format For audio that means channel sample they are references to shared objects When the negotiation mechanism computes the intersection of the formats supported at each end of a all references to both lists are replaced with a reference to the intersection And when a single format is eventually chosen for a link amongst the remaining all references to the list are updated That means that if a filter requires that its input and output have the same format amongst a supported all it has to do is use a reference to the same list of formats query_formats can leave some formats unset and return AVERROR(EAGAIN) to cause the negotiation mechanism toagain later. That can be used by filters with complex requirements to use the format negotiated on one link to set the formats supported on another. Frame references ownership and permissions
opt.h
rdft_horizontal8
static int rdft_horizontal8(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:127
fftfilt_outputs
static const AVFilterPad fftfilt_outputs[]
Definition: vf_fftfilt.c:591
rdft_vertical
static int rdft_vertical(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:435
U
@ U
Definition: vf_fftfilt.c:86
out
FILE * out
Definition: movenc.c:54
V
@ V
Definition: vf_fftfilt.c:86
ff_filter_frame
int ff_filter_frame(AVFilterLink *link, AVFrame *frame)
Send a frame of data to the next filter.
Definition: avfilter.c:999
av_pix_fmt_desc_get
const AVPixFmtDescriptor * av_pix_fmt_desc_get(enum AVPixelFormat pix_fmt)
Definition: pixdesc.c:2662
config_props
static int config_props(AVFilterLink *inlink)
Definition: vf_fftfilt.c:305
FFTFILTContext::rdft_vdata_in
float * rdft_vdata_in[MAX_PLANES]
Definition: vf_fftfilt.c:70
irdft_vertical
static int irdft_vertical(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:454
FILTER_PIXFMTS_ARRAY
#define FILTER_PIXFMTS_ARRAY(array)
Definition: internal.h:170
AVTXContext
Definition: tx_priv.h:201
inlink
The exact code depends on how similar the blocks are and how related they are to the and needs to apply these operations to the correct inlink or outlink if there are several Macros are available to factor that when no extra processing is inlink
Definition: filter_design.txt:212
FFTFILTContext::ivtx_fn
av_tx_fn ivtx_fn
Definition: vf_fftfilt.c:61
ff_clz
#define ff_clz
Definition: intmath.h:142
av_frame_free
void av_frame_free(AVFrame **frame)
Free the frame and any dynamically allocated objects in it, e.g.
Definition: frame.c:111
AVFrame
This structure describes decoded (raw) audio or video data.
Definition: frame.h:325
pixdesc.h
w
uint8_t w
Definition: llviddspenc.c:38
FFTFILTContext::rdft_vbits
int rdft_vbits[MAX_PLANES]
Definition: vf_fftfilt.c:64
AVOption
AVOption.
Definition: opt.h:251
AV_PIX_FMT_YUV420P10
#define AV_PIX_FMT_YUV420P10
Definition: pixfmt.h:406
irdft_horizontal16
static int irdft_horizontal16(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:219
max
#define max(a, b)
Definition: cuda_runtime.h:33
AVFilter::name
const char * name
Filter name.
Definition: avfilter.h:175
VAR_HS
@ VAR_HS
Definition: vf_fftfilt.c:84
av_tx_init
av_cold int av_tx_init(AVTXContext **ctx, av_tx_fn *tx, enum AVTXType type, int inv, int len, const void *scale, uint64_t flags)
Initialize a transform context with the given configuration (i)MDCTs with an odd length are currently...
Definition: tx.c:649
func2_names
static const char *const func2_names[]
Definition: af_afftfilt.c:98
AV_PIX_FMT_GRAY9
#define AV_PIX_FMT_GRAY9
Definition: pixfmt.h:386
AVFrame::data
uint8_t * data[AV_NUM_DATA_POINTERS]
pointer to the picture/channel planes.
Definition: frame.h:346
uninit
static av_cold void uninit(AVFilterContext *ctx)
Definition: vf_fftfilt.c:543
FFTFILTContext::irdft_horizontal
int(* irdft_horizontal)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:80
av_expr_parse
int av_expr_parse(AVExpr **expr, const char *s, const char *const *const_names, const char *const *func1_names, double(*const *funcs1)(void *, double), const char *const *func2_names, double(*const *funcs2)(void *, double, double), int log_offset, void *log_ctx)
Parse an expression.
Definition: eval.c:685
av_pix_fmt_count_planes
int av_pix_fmt_count_planes(enum AVPixelFormat pix_fmt)
Definition: pixdesc.c:2702
FFTFILTContext::nb_threads
int nb_threads
Definition: vf_fftfilt.c:51
EVAL_MODE_FRAME
@ EVAL_MODE_FRAME
Definition: vf_fftfilt.c:41
do_eval
static void do_eval(FFTFILTContext *s, AVFilterLink *inlink, int plane)
Definition: vf_fftfilt.c:284
AV_PIX_FMT_YUV422P9
#define AV_PIX_FMT_YUV422P9
Definition: pixfmt.h:404
MAX_THREADS
#define MAX_THREADS
Definition: vf_fftfilt.c:36
FFTFILTContext::rdft_vstride
size_t rdft_vstride[MAX_PLANES]
Definition: vf_fftfilt.c:66
scale
static av_always_inline float scale(float x, float s)
Definition: vf_v360.c:1389
AV_PIX_FMT_GRAY16
#define AV_PIX_FMT_GRAY16
Definition: pixfmt.h:390
av_expr_free
void av_expr_free(AVExpr *e)
Free a parsed expression previously created with av_expr_parse().
Definition: eval.c:336
copy_rev
static void copy_rev(float *dest, int w, int w2)
Definition: vf_fftfilt.c:116
AVFilterPad
A filter pad used for either input or output.
Definition: internal.h:49
FFTFILTContext::planewidth
int planewidth[MAX_PLANES]
Definition: vf_fftfilt.c:52
FFTFILTContext::eval_mode
int eval_mode
Definition: vf_fftfilt.c:48
AV_PIX_FMT_YUV444P10
#define AV_PIX_FMT_YUV444P10
Definition: pixfmt.h:409
func2
static double(*const func2[])(void *, double, double)
Definition: af_afftfilt.c:99
FFTFILTContext::ivrdft
AVTXContext * ivrdft[MAX_THREADS][MAX_PLANES]
Definition: vf_fftfilt.c:58
av_cold
#define av_cold
Definition: attributes.h:90
AV_PIX_FMT_YUV422P16
#define AV_PIX_FMT_YUV422P16
Definition: pixfmt.h:418
av_tx_fn
void(* av_tx_fn)(AVTXContext *s, void *out, void *in, ptrdiff_t stride)
Function pointer to a function to perform the transform.
Definition: tx.h:111
FFTFILTContext::rdft_horizontal
int(* rdft_horizontal)(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:79
AV_PIX_FMT_YUVJ422P
@ AV_PIX_FMT_YUVJ422P
planar YUV 4:2:2, 16bpp, full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV422P and setting col...
Definition: pixfmt.h:79
s
#define s(width, name)
Definition: cbs_vp9.c:256
AV_PIX_FMT_YUV444P16
#define AV_PIX_FMT_YUV444P16
Definition: pixfmt.h:419
AV_CEIL_RSHIFT
#define AV_CEIL_RSHIFT(a, b)
Definition: common.h:50
FFTFILTContext::rdft_vlen
size_t rdft_vlen[MAX_PLANES]
Definition: vf_fftfilt.c:68
FFTFILTContext::dc
int dc[MAX_PLANES]
Definition: vf_fftfilt.c:74
slice_end
static int slice_end(AVCodecContext *avctx, AVFrame *pict)
Handle slice ends.
Definition: mpeg12dec.c:2013
AV_PIX_FMT_YUV420P9
#define AV_PIX_FMT_YUV420P9
Definition: pixfmt.h:403
AV_PIX_FMT_YUV420P16
#define AV_PIX_FMT_YUV420P16
Definition: pixfmt.h:417
ctx
AVFormatContext * ctx
Definition: movenc.c:48
av_expr_eval
double av_expr_eval(AVExpr *e, const double *const_values, void *opaque)
Evaluate a previously parsed expression.
Definition: eval.c:766
MAX_PLANES
#define MAX_PLANES
Definition: vf_fftfilt.c:37
OFFSET
#define OFFSET(x)
Definition: vf_fftfilt.c:88
AV_PIX_FMT_GRAY14
#define AV_PIX_FMT_GRAY14
Definition: pixfmt.h:389
FFTFILTContext::depth
int depth
Definition: vf_fftfilt.c:49
AVExpr
Definition: eval.c:157
AV_PIX_FMT_YUV420P
@ AV_PIX_FMT_YUV420P
planar YUV 4:2:0, 12bpp, (1 Cr & Cb sample per 2x2 Y samples)
Definition: pixfmt.h:66
ff_vf_fftfilt
const AVFilter ff_vf_fftfilt
Definition: vf_fftfilt.c:598
FILTER_INPUTS
#define FILTER_INPUTS(array)
Definition: internal.h:190
AV_PIX_FMT_YUVJ444P
@ AV_PIX_FMT_YUVJ444P
planar YUV 4:4:4, 24bpp, full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV444P and setting col...
Definition: pixfmt.h:80
arg
const char * arg
Definition: jacosubdec.c:67
AV_PIX_FMT_GRAY10
#define AV_PIX_FMT_GRAY10
Definition: pixfmt.h:387
FFTFILTContext::ihtx_fn
av_tx_fn ihtx_fn
Definition: vf_fftfilt.c:60
VAR_VARS_NB
@ VAR_VARS_NB
Definition: vf_fftfilt.c:84
irdft_horizontal8
static int irdft_horizontal8(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:189
AVClass
Describe the class of an AVClass context structure.
Definition: log.h:66
NULL
#define NULL
Definition: coverity.c:32
FFTFILTContext::planeheight
int planeheight[MAX_PLANES]
Definition: vf_fftfilt.c:53
av_frame_copy_props
int av_frame_copy_props(AVFrame *dst, const AVFrame *src)
Copy only "metadata" fields from src to dst.
Definition: frame.c:596
FFTFILTContext::rdft_hlen
size_t rdft_hlen[MAX_PLANES]
Definition: vf_fftfilt.c:67
AV_PIX_FMT_YUVJ420P
@ AV_PIX_FMT_YUVJ420P
planar YUV 4:2:0, 12bpp, full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV420P and setting col...
Definition: pixfmt.h:78
copy_horizontal
static int copy_horizontal(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:473
VAR_WS
@ VAR_WS
Definition: vf_fftfilt.c:84
double
double
Definition: af_crystalizer.c:132
av_cpu_max_align
size_t av_cpu_max_align(void)
Get the maximum data alignment that may be required by FFmpeg.
Definition: cpu.c:254
FFTFILTContext::htx_fn
av_tx_fn htx_fn
Definition: vf_fftfilt.c:60
AV_PIX_FMT_YUV422P10
#define AV_PIX_FMT_YUV422P10
Definition: pixfmt.h:407
AV_PIX_FMT_GRAY8
@ AV_PIX_FMT_GRAY8
Y , 8bpp.
Definition: pixfmt.h:74
initialize
static av_cold int initialize(AVFilterContext *ctx)
Definition: vf_fftfilt.c:250
fftfilt_options
static const AVOption fftfilt_options[]
Definition: vf_fftfilt.c:91
weight
static int weight(int i, int blen, int offset)
Definition: diracdec.c:1562
eval.h
FFTFILTContext::rdft_vdata_out
float * rdft_vdata_out[MAX_PLANES]
Definition: vf_fftfilt.c:72
FFTFILTContext::weight_str
char * weight_str[MAX_PLANES]
Definition: vf_fftfilt.c:75
dc
Tag MUST be and< 10hcoeff half pel interpolation filter coefficients, hcoeff[0] are the 2 middle coefficients[1] are the next outer ones and so on, resulting in a filter like:...eff[2], hcoeff[1], hcoeff[0], hcoeff[0], hcoeff[1], hcoeff[2] ... the sign of the coefficients is not explicitly stored but alternates after each coeff and coeff[0] is positive, so ...,+,-,+,-,+,+,-,+,-,+,... hcoeff[0] is not explicitly stored but found by subtracting the sum of all stored coefficients with signs from 32 hcoeff[0]=32 - hcoeff[1] - hcoeff[2] - ... a good choice for hcoeff and htaps is htaps=6 hcoeff={40,-10, 2} an alternative which requires more computations at both encoder and decoder side and may or may not be better is htaps=8 hcoeff={42,-14, 6,-2}ref_frames minimum of the number of available reference frames and max_ref_frames for example the first frame after a key frame always has ref_frames=1spatial_decomposition_type wavelet type 0 is a 9/7 symmetric compact integer wavelet 1 is a 5/3 symmetric compact integer wavelet others are reserved stored as delta from last, last is reset to 0 if always_reset||keyframeqlog quality(logarithmic quantizer scale) stored as delta from last, last is reset to 0 if always_reset||keyframemv_scale stored as delta from last, last is reset to 0 if always_reset||keyframe FIXME check that everything works fine if this changes between framesqbias dequantization bias stored as delta from last, last is reset to 0 if always_reset||keyframeblock_max_depth maximum depth of the block tree stored as delta from last, last is reset to 0 if always_reset||keyframequant_table quantization tableHighlevel bitstream structure:==============================--------------------------------------------|Header|--------------------------------------------|------------------------------------|||Block0||||split?||||yes no||||......... intra?||||:Block01 :yes no||||:Block02 :....... ..........||||:Block03 ::y DC ::ref index:||||:Block04 ::cb DC ::motion x :||||......... :cr DC ::motion y :||||....... ..........|||------------------------------------||------------------------------------|||Block1|||...|--------------------------------------------|------------ ------------ ------------|||Y subbands||Cb subbands||Cr subbands||||--- ---||--- ---||--- ---|||||LL0||HL0||||LL0||HL0||||LL0||HL0|||||--- ---||--- ---||--- ---||||--- ---||--- ---||--- ---|||||LH0||HH0||||LH0||HH0||||LH0||HH0|||||--- ---||--- ---||--- ---||||--- ---||--- ---||--- ---|||||HL1||LH1||||HL1||LH1||||HL1||LH1|||||--- ---||--- ---||--- ---||||--- ---||--- ---||--- ---|||||HH1||HL2||||HH1||HL2||||HH1||HL2|||||...||...||...|||------------ ------------ ------------|--------------------------------------------Decoding process:=================------------|||Subbands|------------||||------------|Intra DC||||LL0 subband prediction ------------|\ Dequantization ------------------- \||Reference frames|\ IDWT|------- -------|Motion \|||Frame 0||Frame 1||Compensation . OBMC v -------|------- -------|--------------. \------> Frame n output Frame Frame<----------------------------------/|...|------------------- Range Coder:============Binary Range Coder:------------------- The implemented range coder is an adapted version based upon "Range encoding: an algorithm for removing redundancy from a digitised message." by G. N. N. Martin. The symbols encoded by the Snow range coder are bits(0|1). The associated probabilities are not fix but change depending on the symbol mix seen so far. bit seen|new state ---------+----------------------------------------------- 0|256 - state_transition_table[256 - old_state];1|state_transition_table[old_state];state_transition_table={ 0, 0, 0, 0, 0, 0, 0, 0, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 190, 191, 192, 194, 194, 195, 196, 197, 198, 199, 200, 201, 202, 202, 204, 205, 206, 207, 208, 209, 209, 210, 211, 212, 213, 215, 215, 216, 217, 218, 219, 220, 220, 222, 223, 224, 225, 226, 227, 227, 229, 229, 230, 231, 232, 234, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 248, 0, 0, 0, 0, 0, 0, 0};FIXME Range Coding of integers:------------------------- FIXME Neighboring Blocks:===================left and top are set to the respective blocks unless they are outside of the image in which case they are set to the Null block top-left is set to the top left block unless it is outside of the image in which case it is set to the left block if this block has no larger parent block or it is at the left side of its parent block and the top right block is not outside of the image then the top right block is used for top-right else the top-left block is used Null block y, cb, cr are 128 level, ref, mx and my are 0 Motion Vector Prediction:=========================1. the motion vectors of all the neighboring blocks are scaled to compensate for the difference of reference frames scaled_mv=(mv *(256 *(current_reference+1)/(mv.reference+1))+128)> the median of the scaled top and top right vectors is used as motion vector prediction the used motion vector is the sum of the predictor and(mvx_diff, mvy_diff) *mv_scale Intra DC Prediction block[y][x] dc[1]
Definition: snow.txt:400
NULL_IF_CONFIG_SMALL
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification.
Definition: internal.h:117
cpu.h
FFTFILTContext
Definition: vf_fftfilt.c:45
FFTFILTContext::vtx_fn
av_tx_fn vtx_fn
Definition: vf_fftfilt.c:61
AV_PIX_FMT_YUV422P12
#define AV_PIX_FMT_YUV422P12
Definition: pixfmt.h:411
AV_PIX_FMT_YUV444P12
#define AV_PIX_FMT_YUV444P12
Definition: pixfmt.h:413
height
#define height
copy_vertical
static int copy_vertical(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:410
av_tx_uninit
av_cold void av_tx_uninit(AVTXContext **ctx)
Frees a context and sets *ctx to NULL, does nothing when *ctx == NULL.
Definition: tx.c:251
internal.h
AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC
#define AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC
Some filters support a generic "enable" expression option that can be used to enable or disable a fil...
Definition: avfilter.h:152
VAR_H
@ VAR_H
Definition: vf_fftfilt.c:84
lrintf
#define lrintf(x)
Definition: libm_mips.h:72
i
#define i(width, name, range_min, range_max)
Definition: cbs_h2645.c:269
FFTFILTContext::rdft_hbits
int rdft_hbits[MAX_PLANES]
Definition: vf_fftfilt.c:63
FFTFILTContext::rdft_hdata_out
float * rdft_hdata_out[MAX_PLANES]
Definition: vf_fftfilt.c:71
common.h
ff_filter_get_nb_threads
int ff_filter_get_nb_threads(AVFilterContext *ctx)
Get number of threads for current filter instance.
Definition: avfilter.c:783
EvalMode
EvalMode
Definition: af_volume.h:39
FFMIN
#define FFMIN(a, b)
Definition: macros.h:49
AVFILTER_DEFINE_CLASS
AVFILTER_DEFINE_CLASS(fftfilt)
AVFilterPad::name
const char * name
Pad name.
Definition: internal.h:55
var_names
static const char *const var_names[]
Definition: vf_fftfilt.c:83
FFTFILTContext::weight_expr
AVExpr * weight_expr[MAX_PLANES]
Definition: vf_fftfilt.c:76
av_calloc
void * av_calloc(size_t nmemb, size_t size)
Definition: mem.c:272
AV_PIX_FMT_YUV444P9
#define AV_PIX_FMT_YUV444P9
Definition: pixfmt.h:405
weight_U
static double weight_U(void *priv, double x, double y)
Definition: vf_fftfilt.c:113
AVFilter
Filter definition.
Definition: avfilter.h:171
FFTFILTContext::vrdft
AVTXContext * vrdft[MAX_THREADS][MAX_PLANES]
Definition: vf_fftfilt.c:56
FFTFILTContext::weight
double * weight[MAX_PLANES]
Definition: vf_fftfilt.c:77
ret
ret
Definition: filter_design.txt:187
fftfilt_inputs
static const AVFilterPad fftfilt_inputs[]
Definition: vf_fftfilt.c:582
AV_PIX_FMT_YUV420P12
#define AV_PIX_FMT_YUV420P12
Definition: pixfmt.h:410
AV_TX_FLOAT_RDFT
@ AV_TX_FLOAT_RDFT
Real to complex and complex to real DFTs.
Definition: tx.h:88
AV_PIX_FMT_YUV422P14
#define AV_PIX_FMT_YUV422P14
Definition: pixfmt.h:415
filter_frame
static int filter_frame(AVFilterLink *inlink, AVFrame *in)
Definition: vf_fftfilt.c:495
FFTFILTContext::rdft_hdata_in
float * rdft_hdata_in[MAX_PLANES]
Definition: vf_fftfilt.c:69
VAR_X
@ VAR_X
Definition: vf_fftfilt.c:84
AV_PIX_FMT_NONE
@ AV_PIX_FMT_NONE
Definition: pixfmt.h:65
AV_OPT_TYPE_INT
@ AV_OPT_TYPE_INT
Definition: opt.h:225
values
these buffered frames must be flushed immediately if a new input produces new the filter must not call request_frame to get more It must just process the frame or queue it The task of requesting more frames is left to the filter s request_frame method or the application If a filter has several the filter must be ready for frames arriving randomly on any input any filter with several inputs will most likely require some kind of queuing mechanism It is perfectly acceptable to have a limited queue and to drop frames when the inputs are too unbalanced request_frame For filters that do not use the this method is called when a frame is wanted on an output For a it should directly call filter_frame on the corresponding output For a if there are queued frames already one of these frames should be pushed If the filter should request a frame on one of its repeatedly until at least one frame has been pushed Return values
Definition: filter_design.txt:263
rdft_horizontal16
static int rdft_horizontal16(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:158
pixel_fmts_fftfilt
static enum AVPixelFormat pixel_fmts_fftfilt[]
Definition: vf_fftfilt.c:563
av_clip_uint8
#define av_clip_uint8
Definition: common.h:101
AV_PIX_FMT_YUV444P
@ AV_PIX_FMT_YUV444P
planar YUV 4:4:4, 24bpp, (1 Cr & Cb sample per 1x1 Y samples)
Definition: pixfmt.h:71
AVFilterContext
An instance of a filter.
Definition: avfilter.h:408
weight_V
static double weight_V(void *priv, double x, double y)
Definition: vf_fftfilt.c:114
VAR_Y
@ VAR_Y
Definition: vf_fftfilt.c:84
AVFILTER_FLAG_SLICE_THREADS
#define AVFILTER_FLAG_SLICE_THREADS
The filter supports multithreading by splitting frames into multiple parts and processing them concur...
Definition: avfilter.h:127
av_strdup
char * av_strdup(const char *s)
Duplicate a string.
Definition: mem.c:280
desc
const char * desc
Definition: libsvtav1.c:83
AVMEDIA_TYPE_VIDEO
@ AVMEDIA_TYPE_VIDEO
Definition: avutil.h:201
AV_PIX_FMT_YUV422P
@ AV_PIX_FMT_YUV422P
planar YUV 4:2:2, 16bpp, (1 Cr & Cb sample per 2x1 Y samples)
Definition: pixfmt.h:70
VAR_N
@ VAR_N
Definition: vf_fftfilt.c:84
AVPixFmtDescriptor
Descriptor that unambiguously describes how the bits of a pixel are stored in the up to 4 data planes...
Definition: pixdesc.h:69
EVAL_MODE_INIT
@ EVAL_MODE_INIT
Definition: vf_fftfilt.c:40
FFALIGN
#define FFALIGN(x, a)
Definition: macros.h:78
lum
static double lum(void *priv, double x, double y, int plane)
Definition: vf_fftfilt.c:106
FILTER_OUTPUTS
#define FILTER_OUTPUTS(array)
Definition: internal.h:191
av_freep
#define av_freep(p)
Definition: tableprint_vlc.h:34
src
INIT_CLIP pixel * src
Definition: h264pred_template.c:418
VAR_W
@ VAR_W
Definition: vf_fftfilt.c:84
imgutils.h
AVERROR_BUG
#define AVERROR_BUG
Internal bug, also see AVERROR_BUG2.
Definition: error.h:52
AVFrame::linesize
int linesize[AV_NUM_DATA_POINTERS]
For video, a positive or negative value, which is typically indicating the size in bytes of each pict...
Definition: frame.h:370
FFTFILTContext::rdft_hstride
size_t rdft_hstride[MAX_PLANES]
Definition: vf_fftfilt.c:65
weight_Y
static double weight_Y(void *priv, double x, double y)
Definition: vf_fftfilt.c:112
h
h
Definition: vp9dsp_template.c:2038
AV_PIX_FMT_YUV444P14
#define AV_PIX_FMT_YUV444P14
Definition: pixfmt.h:416
AV_OPT_TYPE_STRING
@ AV_OPT_TYPE_STRING
Definition: opt.h:229
FLAGS
#define FLAGS
Definition: vf_fftfilt.c:89
EVAL_MODE_NB
@ EVAL_MODE_NB
Definition: vf_fftfilt.c:42
AV_PIX_FMT_GRAY12
#define AV_PIX_FMT_GRAY12
Definition: pixfmt.h:388
Y
@ Y
Definition: vf_fftfilt.c:86
ff_filter_execute
static av_always_inline int ff_filter_execute(AVFilterContext *ctx, avfilter_action_func *func, void *arg, int *ret, int nb_jobs)
Definition: internal.h:142
FFTFILTContext::nb_planes
int nb_planes
Definition: vf_fftfilt.c:50
int
int
Definition: ffmpeg_filter.c:153
FFTFILTContext::ihrdft
AVTXContext * ihrdft[MAX_THREADS][MAX_PLANES]
Definition: vf_fftfilt.c:57
AV_OPT_TYPE_CONST
@ AV_OPT_TYPE_CONST
Definition: opt.h:234
av_log2
int av_log2(unsigned v)
Definition: intmath.c:26
AV_PIX_FMT_YUV420P14
#define AV_PIX_FMT_YUV420P14
Definition: pixfmt.h:414
tx.h
multiply_data
static int multiply_data(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
Definition: vf_fftfilt.c:389
Generated on Wed Aug 24 2022 21:42:05 for FFmpeg by   doxygen 1.8.17