FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
vf_removegrain.c
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2012 Laurent de Soras
3  * Copyright (c) 2013 Fredrik Mellbin
4  * Copyright (c) 2015 Paul B Mahol
5  * Copyright (c) 2015 James Darnley
6  *
7  * This file is part of FFmpeg.
8  *
9  * FFmpeg is free software; you can redistribute it and/or
10  * modify it under the terms of the GNU Lesser General Public
11  * License as published by the Free Software Foundation; either
12  * version 2.1 of the License, or (at your option) any later version.
13  *
14  * FFmpeg is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17  * Lesser General Public License for more details.
18  *
19  * You should have received a copy of the GNU Lesser General Public
20  * License along with FFmpeg; if not, write to the Free Software
21  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22  */
23 
24 #include "libavutil/imgutils.h"
25 #include "libavutil/opt.h"
26 #include "libavutil/pixdesc.h"
27 #include "avfilter.h"
28 #include "formats.h"
29 #include "internal.h"
30 #include "removegrain.h"
31 #include "video.h"
32 
33 #define OFFSET(x) offsetof(RemoveGrainContext, x)
34 #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM
35 
36 static const AVOption removegrain_options[] = {
37  { "m0", "set mode for 1st plane", OFFSET(mode[0]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
38  { "m1", "set mode for 2nd plane", OFFSET(mode[1]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
39  { "m2", "set mode for 3rd plane", OFFSET(mode[2]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
40  { "m3", "set mode for 4th plane", OFFSET(mode[3]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
41  {NULL}
42 };
43 
44 AVFILTER_DEFINE_CLASS(removegrain);
45 
47 {
48  static const enum AVPixelFormat pix_fmts[] = {
57  };
58 
59  AVFilterFormats *fmts_list = ff_make_format_list(pix_fmts);
60  if (!fmts_list)
61  return AVERROR(ENOMEM);
62  return ff_set_common_formats(ctx, fmts_list);
63 }
64 
65 #define REMOVE_GRAIN_SORT_AXIS \
66  const int ma1 = FFMAX(a1, a8); \
67  const int mi1 = FFMIN(a1, a8); \
68  const int ma2 = FFMAX(a2, a7); \
69  const int mi2 = FFMIN(a2, a7); \
70  const int ma3 = FFMAX(a3, a6); \
71  const int mi3 = FFMIN(a3, a6); \
72  const int ma4 = FFMAX(a4, a5); \
73  const int mi4 = FFMIN(a4, a5);
74 
75 static int mode01(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
76 {
77  const int mi = FFMIN(FFMIN(FFMIN(a1, a2), FFMIN(a3, a4)), FFMIN(FFMIN(a5, a6), FFMIN(a7, a8)));
78  const int ma = FFMAX(FFMAX(FFMAX(a1, a2), FFMAX(a3, a4)), FFMAX(FFMAX(a5, a6), FFMAX(a7, a8)));
79 
80  return av_clip(c, mi, ma);
81 }
82 
83 static int cmp_int(const void *p1, const void *p2)
84 {
85  int left = *(const int *)p1;
86  int right = *(const int *)p2;
87 
88  return ((left > right) - (left < right));
89 }
90 
91 static int mode02(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
92 {
93  int a[8] = { a1, a2, a3, a4, a5, a6, a7, a8 };
94 
95  qsort(&a, 8, sizeof(a[0]), cmp_int);
96 
97  return av_clip(c, a[2 - 1 ], a[7 - 1]);
98 }
99 
100 static int mode03(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
101 {
102  int a[8] = { a1, a2, a3, a4, a5, a6, a7, a8 };
103 
104  qsort(&a, 8, sizeof(a[0]), cmp_int);
105 
106  return av_clip(c, a[3 - 1 ], a[6 - 1]);
107 }
108 
109 static int mode04(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
110 {
111  int a[8] = { a1, a2, a3, a4, a5, a6, a7, a8 };
112 
113  qsort(&a, 8, sizeof(a[0]), cmp_int);
114 
115  return av_clip(c, a[4 - 1 ], a[5 - 1]);
116 }
117 
118 static int mode05(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
119 {
121 
122  const int c1 = FFABS(c - av_clip(c, mi1, ma1));
123  const int c2 = FFABS(c - av_clip(c, mi2, ma2));
124  const int c3 = FFABS(c - av_clip(c, mi3, ma3));
125  const int c4 = FFABS(c - av_clip(c, mi4, ma4));
126 
127  const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
128 
129  /* When adding SIMD notice the return order here: 4, 2, 3, 1. */
130  if (mindiff == c4) {
131  return av_clip(c, mi4, ma4);
132  } else if (mindiff == c2) {
133  return av_clip(c, mi2, ma2);
134  } else if (mindiff == c3) {
135  return av_clip(c, mi3, ma3);
136  }
137 
138  return av_clip(c, mi1, ma1);
139 }
140 
141 static int mode06(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
142 {
144 
145  const int d1 = ma1 - mi1;
146  const int d2 = ma2 - mi2;
147  const int d3 = ma3 - mi3;
148  const int d4 = ma4 - mi4;
149 
150  const int cli1 = av_clip(c, mi1, ma1);
151  const int cli2 = av_clip(c, mi2, ma2);
152  const int cli3 = av_clip(c, mi3, ma3);
153  const int cli4 = av_clip(c, mi4, ma4);
154 
155  const int c1 = av_clip_uint16((FFABS(c - cli1) << 1) + d1);
156  const int c2 = av_clip_uint16((FFABS(c - cli2) << 1) + d2);
157  const int c3 = av_clip_uint16((FFABS(c - cli3) << 1) + d3);
158  const int c4 = av_clip_uint16((FFABS(c - cli4) << 1) + d4);
159 
160  const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
161 
162  if (mindiff == c4) {
163  return cli4;
164  } else if (mindiff == c2) {
165  return cli2;
166  } else if (mindiff == c3) {
167  return cli3;
168  }
169 
170  return cli1;
171 }
172 
173 static int mode07(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
174 {
176 
177  const int d1 = ma1 - mi1;
178  const int d2 = ma2 - mi2;
179  const int d3 = ma3 - mi3;
180  const int d4 = ma4 - mi4;
181 
182  const int cli1 = av_clip(c, mi1, ma1);
183  const int cli2 = av_clip(c, mi2, ma2);
184  const int cli3 = av_clip(c, mi3, ma3);
185  const int cli4 = av_clip(c, mi4, ma4);
186 
187  const int c1 = FFABS(c - cli1) + d1;
188  const int c2 = FFABS(c - cli2) + d2;
189  const int c3 = FFABS(c - cli3) + d3;
190  const int c4 = FFABS(c - cli4) + d4;
191 
192  const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
193 
194  if (mindiff == c4) {
195  return cli4;
196  } else if (mindiff == c2) {
197  return cli2;
198  } else if (mindiff == c3) {
199  return cli3;
200  }
201 
202  return cli1;
203 }
204 
205 static int mode08(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
206 {
208 
209  const int d1 = ma1 - mi1;
210  const int d2 = ma2 - mi2;
211  const int d3 = ma3 - mi3;
212  const int d4 = ma4 - mi4;
213 
214  const int cli1 = av_clip(c, mi1, ma1);
215  const int cli2 = av_clip(c, mi2, ma2);
216  const int cli3 = av_clip(c, mi3, ma3);
217  const int cli4 = av_clip(c, mi4, ma4);
218 
219  const int c1 = av_clip_uint16(FFABS(c - cli1) + (d1 << 1));
220  const int c2 = av_clip_uint16(FFABS(c - cli2) + (d2 << 1));
221  const int c3 = av_clip_uint16(FFABS(c - cli3) + (d3 << 1));
222  const int c4 = av_clip_uint16(FFABS(c - cli4) + (d4 << 1));
223 
224  const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
225 
226  if (mindiff == c4) {
227  return cli4;
228  } else if (mindiff == c2) {
229  return cli2;
230  } else if (mindiff == c3) {
231  return cli3;
232  }
233 
234  return cli1;
235 }
236 
237 static int mode09(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
238 {
240 
241  const int d1 = ma1 - mi1;
242  const int d2 = ma2 - mi2;
243  const int d3 = ma3 - mi3;
244  const int d4 = ma4 - mi4;
245 
246  const int mindiff = FFMIN(FFMIN(d1, d2), FFMIN(d3, d4));
247 
248  if (mindiff == d4) {
249  return av_clip(c, mi4, ma4);
250  } else if (mindiff == d2) {
251  return av_clip(c, mi2, ma2);
252  } else if (mindiff == d3) {
253  return av_clip(c, mi3, ma3);
254  }
255 
256  return av_clip(c, mi1, ma1);
257 }
258 
259 static int mode10(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
260 {
261  const int d1 = FFABS(c - a1);
262  const int d2 = FFABS(c - a2);
263  const int d3 = FFABS(c - a3);
264  const int d4 = FFABS(c - a4);
265  const int d5 = FFABS(c - a5);
266  const int d6 = FFABS(c - a6);
267  const int d7 = FFABS(c - a7);
268  const int d8 = FFABS(c - a8);
269 
270  const int mindiff = FFMIN(FFMIN(FFMIN(d1, d2), FFMIN(d3, d4)),
271  FFMIN(FFMIN(d5, d6), FFMIN(d7, d8)));
272 
273  if (mindiff == d7) return a7;
274  if (mindiff == d8) return a8;
275  if (mindiff == d6) return a6;
276  if (mindiff == d2) return a2;
277  if (mindiff == d3) return a3;
278  if (mindiff == d1) return a1;
279  if (mindiff == d5) return a5;
280 
281  return a4;
282 }
283 
284 static int mode1112(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
285 {
286  const int sum = 4 * c + 2 * (a2 + a4 + a5 + a7) + a1 + a3 + a6 + a8;
287  const int val = (sum + 8) >> 4;
288 
289  return val;
290 }
291 
292 static int mode1314(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
293 {
294  const int d1 = FFABS(a1 - a8);
295  const int d2 = FFABS(a2 - a7);
296  const int d3 = FFABS(a3 - a6);
297 
298  const int mindiff = FFMIN(FFMIN(d1, d2), d3);
299 
300  if (mindiff == d2) {
301  return (a2 + a7 + 1) >> 1;
302  }
303  if (mindiff == d3) {
304  return (a3 + a6 + 1) >> 1;
305  }
306 
307  return (a1 + a8 + 1) >> 1;
308 }
309 
310 static int mode1516(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
311 {
312  const int d1 = FFABS(a1 - a8);
313  const int d2 = FFABS(a2 - a7);
314  const int d3 = FFABS(a3 - a6);
315 
316  const int mindiff = FFMIN(FFMIN(d1, d2), d3);
317  const int average = (2 * (a2 + a7) + a1 + a3 + a6 + a8 + 4) >> 3;
318 
319  if (mindiff == d2) {
320  return av_clip(average, FFMIN(a2, a7), FFMAX(a2, a7));
321  }
322  if (mindiff == d3) {
323  return av_clip(average, FFMIN(a3, a6), FFMAX(a3, a6));
324  }
325 
326  return av_clip(average, FFMIN(a1, a8), FFMAX(a1, a8));
327 }
328 
329 static int mode17(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
330 {
332 
333  const int l = FFMAX(FFMAX(mi1, mi2), FFMAX(mi3, mi4));
334  const int u = FFMIN(FFMIN(ma1, ma2), FFMIN(ma3, ma4));
335 
336  return av_clip(c, FFMIN(l, u), FFMAX(l, u));
337 }
338 
339 static int mode18(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
340 {
341  const int d1 = FFMAX(FFABS(c - a1), FFABS(c - a8));
342  const int d2 = FFMAX(FFABS(c - a2), FFABS(c - a7));
343  const int d3 = FFMAX(FFABS(c - a3), FFABS(c - a6));
344  const int d4 = FFMAX(FFABS(c - a4), FFABS(c - a5));
345 
346  const int mindiff = FFMIN(FFMIN(d1, d2), FFMIN(d3, d4));
347 
348  if (mindiff == d4) {
349  return av_clip(c, FFMIN(a4, a5), FFMAX(a4, a5));
350  }
351  if (mindiff == d2) {
352  return av_clip(c, FFMIN(a2, a7), FFMAX(a2, a7));
353  }
354  if (mindiff == d3) {
355  return av_clip(c, FFMIN(a3, a6), FFMAX(a3, a6));
356  }
357 
358  return av_clip(c, FFMIN(a1, a8), FFMAX(a1, a8));
359 }
360 
361 static int mode19(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
362 {
363  const int sum = a1 + a2 + a3 + a4 + a5 + a6 + a7 + a8;
364  const int val = (sum + 4) >> 3;
365 
366  return val;
367 }
368 
369 static int mode20(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
370 {
371  const int sum = a1 + a2 + a3 + a4 + c + a5 + a6 + a7 + a8;
372  const int val = (sum + 4) / 9;
373 
374  return val;
375 }
376 
377 static int mode21(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
378 {
379  const int l1l = (a1 + a8) >> 1;
380  const int l2l = (a2 + a7) >> 1;
381  const int l3l = (a3 + a6) >> 1;
382  const int l4l = (a4 + a5) >> 1;
383 
384  const int l1h = (a1 + a8 + 1) >> 1;
385  const int l2h = (a2 + a7 + 1) >> 1;
386  const int l3h = (a3 + a6 + 1) >> 1;
387  const int l4h = (a4 + a5 + 1) >> 1;
388 
389  const int mi = FFMIN(FFMIN(l1l, l2l), FFMIN(l3l, l4l));
390  const int ma = FFMAX(FFMAX(l1h, l2h), FFMAX(l3h, l4h));
391 
392  return av_clip(c, mi, ma);
393 }
394 
395 static int mode22(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
396 {
397  const int l1 = (a1 + a8 + 1) >> 1;
398  const int l2 = (a2 + a7 + 1) >> 1;
399  const int l3 = (a3 + a6 + 1) >> 1;
400  const int l4 = (a4 + a5 + 1) >> 1;
401 
402  const int mi = FFMIN(FFMIN(l1, l2), FFMIN(l3, l4));
403  const int ma = FFMAX(FFMAX(l1, l2), FFMAX(l3, l4));
404 
405  return av_clip(c, mi, ma);
406 }
407 
408 static int mode23(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
409 {
411 
412  const int linediff1 = ma1 - mi1;
413  const int linediff2 = ma2 - mi2;
414  const int linediff3 = ma3 - mi3;
415  const int linediff4 = ma4 - mi4;
416 
417  const int u1 = FFMIN(c - ma1, linediff1);
418  const int u2 = FFMIN(c - ma2, linediff2);
419  const int u3 = FFMIN(c - ma3, linediff3);
420  const int u4 = FFMIN(c - ma4, linediff4);
421  const int u = FFMAX(FFMAX(FFMAX(u1, u2), FFMAX(u3, u4)), 0);
422 
423  const int d1 = FFMIN(mi1 - c, linediff1);
424  const int d2 = FFMIN(mi2 - c, linediff2);
425  const int d3 = FFMIN(mi3 - c, linediff3);
426  const int d4 = FFMIN(mi4 - c, linediff4);
427  const int d = FFMAX(FFMAX(FFMAX(d1, d2), FFMAX(d3, d4)), 0);
428 
429  return c - u + d; // This probably will never overflow.
430 }
431 
432 static int mode24(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
433 {
435 
436  const int linediff1 = ma1 - mi1;
437  const int linediff2 = ma2 - mi2;
438  const int linediff3 = ma3 - mi3;
439  const int linediff4 = ma4 - mi4;
440 
441  const int tu1 = c - ma1;
442  const int tu2 = c - ma2;
443  const int tu3 = c - ma3;
444  const int tu4 = c - ma4;
445 
446  const int u1 = FFMIN(tu1, linediff1 - tu1);
447  const int u2 = FFMIN(tu2, linediff2 - tu2);
448  const int u3 = FFMIN(tu3, linediff3 - tu3);
449  const int u4 = FFMIN(tu4, linediff4 - tu4);
450  const int u = FFMAX(FFMAX(FFMAX(u1, u2), FFMAX(u3, u4)), 0);
451 
452  const int td1 = mi1 - c;
453  const int td2 = mi2 - c;
454  const int td3 = mi3 - c;
455  const int td4 = mi4 - c;
456 
457  const int d1 = FFMIN(td1, linediff1 - td1);
458  const int d2 = FFMIN(td2, linediff2 - td2);
459  const int d3 = FFMIN(td3, linediff3 - td3);
460  const int d4 = FFMIN(td4, linediff4 - td4);
461  const int d = FFMAX(FFMAX(FFMAX(d1, d2), FFMAX(d3, d4)), 0);
462 
463  return c - u + d; // This probably will never overflow.
464 }
465 
466 static int config_input(AVFilterLink *inlink)
467 {
468  RemoveGrainContext *s = inlink->dst->priv;
469  const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
470  int i;
471 
473 
474  s->planeheight[1] = s->planeheight[2] = FF_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h);
475  s->planeheight[0] = s->planeheight[3] = inlink->h;
476  s->planewidth[1] = s->planewidth[2] = FF_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w);
477  s->planewidth[0] = s->planewidth[3] = inlink->w;
478 
479  for (i = 0; i < s->nb_planes; i++) {
480  switch (s->mode[i]) {
481  case 1: s->rg[i] = mode01; break;
482  case 2: s->rg[i] = mode02; break;
483  case 3: s->rg[i] = mode03; break;
484  case 4: s->rg[i] = mode04; break;
485  case 5: s->rg[i] = mode05; break;
486  case 6: s->rg[i] = mode06; break;
487  case 7: s->rg[i] = mode07; break;
488  case 8: s->rg[i] = mode08; break;
489  case 9: s->rg[i] = mode09; break;
490  case 10: s->rg[i] = mode10; break;
491  case 11: s->rg[i] = mode1112; break;
492  case 12: s->rg[i] = mode1112; break;
493  case 13: s->skip_odd = 1;
494  s->rg[i] = mode1314; break;
495  case 14: s->skip_even = 1;
496  s->rg[i] = mode1314; break;
497  case 15: s->skip_odd = 1;
498  s->rg[i] = mode1516; break;
499  case 16: s->skip_even = 1;
500  s->rg[i] = mode1516; break;
501  case 17: s->rg[i] = mode17; break;
502  case 18: s->rg[i] = mode18; break;
503  case 19: s->rg[i] = mode19; break;
504  case 20: s->rg[i] = mode20; break;
505  case 21: s->rg[i] = mode21; break;
506  case 22: s->rg[i] = mode22; break;
507  case 23: s->rg[i] = mode23; break;
508  case 24: s->rg[i] = mode24; break;
509  }
510  }
511 
512  if (ARCH_X86)
514 
515  return 0;
516 }
517 
518 typedef struct ThreadData {
519  AVFrame *in, *out;
520  int plane;
521 } ThreadData;
522 
523 static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
524 {
525  RemoveGrainContext *s = ctx->priv;
526  ThreadData *td = arg;
527  AVFrame *in = td->in;
528  AVFrame *out = td->out;
529  const int i = td->plane;
530  const int height = s->planeheight[i];
531  const int om = in->linesize[i] - 1;
532  const int o0 = in->linesize[i] ;
533  const int op = in->linesize[i] + 1;
534  int start = (height * jobnr ) / nb_jobs;
535  int end = (height * (jobnr+1)) / nb_jobs;
536  int x, y;
537 
538  start = FFMAX(1, start);
539  end = FFMIN(height-1, end);
540  for (y = start; y < end; y++) {
541  uint8_t *dst = out->data[i];
542  uint8_t *src = in->data[i];
543 
544  src = in->data[i] + y * in->linesize[i];
545  dst = out->data[i] + y * out->linesize[i];
546 
547  if (s->skip_even && !(y & 1)) {
548  memcpy(dst, src, s->planewidth[i]);
549  continue;
550  }
551  if (s->skip_odd && y & 1) {
552  memcpy(dst, src, s->planewidth[i]);
553  continue;
554  }
555 
556  *dst++ = *src++;
557 
558  if (s->fl[i]) {
559  int w_asm = (s->planewidth[i] - 2) & ~15;
560 
561  s->fl[i](dst, src, in->linesize[i], w_asm);
562 
563  x = 1 + w_asm;
564  dst += w_asm;
565  src += w_asm;
566  } else
567  x = 1;
568 
569  for (; x < s->planewidth[i] - 1; x++) {
570  const int a1 = src[-op];
571  const int a2 = src[-o0];
572  const int a3 = src[-om];
573  const int a4 = src[-1 ];
574  const int c = src[ 0 ];
575  const int a5 = src[ 1 ];
576  const int a6 = src[ om];
577  const int a7 = src[ o0];
578  const int a8 = src[ op];
579 
580  const int res = s->rg[i](c, a1, a2, a3, a4, a5, a6, a7, a8);
581 
582  *dst = res;
583  dst++, src++;
584  }
585  dst[0] = src[0];
586  }
587 
588  return 0;
589 }
590 
591 static int filter_frame(AVFilterLink *inlink, AVFrame *in)
592 {
593  AVFilterContext *ctx = inlink->dst;
594  AVFilterLink *outlink = ctx->outputs[0];
595  RemoveGrainContext *s = ctx->priv;
596  ThreadData td;
597  AVFrame *out;
598  int i;
599 
600  out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
601  if (!out) {
602  av_frame_free(&in);
603  return AVERROR(ENOMEM);
604  }
605  av_frame_copy_props(out, in);
606 
607  for (i = 0; i < s->nb_planes; i++) {
608  uint8_t *dst = out->data[i];
609  uint8_t *src = in->data[i];
610 
611  if (s->mode[i] == 0) {
612  av_image_copy_plane(dst, out->linesize[i],
613  src, in->linesize[i],
614  s->planewidth[i], s->planeheight[i]);
615  continue;
616  }
617 
618  memcpy(dst, src, s->planewidth[i]);
619 
620  td.in = in; td.out = out; td.plane = i;
621  ctx->internal->execute(ctx, filter_slice, &td, NULL,
622  FFMIN(s->planeheight[i], ctx->graph->nb_threads));
623 
624  src = in->data[i] + (s->planeheight[i] - 1) * in->linesize[i];
625  dst = out->data[i] + (s->planeheight[i] - 1) * out->linesize[i];
626  memcpy(dst, src, s->planewidth[i]);
627  }
628 
629  av_frame_free(&in);
630  return ff_filter_frame(outlink, out);
631 }
632 
633 static const AVFilterPad removegrain_inputs[] = {
634  {
635  .name = "default",
636  .type = AVMEDIA_TYPE_VIDEO,
637  .filter_frame = filter_frame,
638  .config_props = config_input,
639  },
640  { NULL }
641 };
642 
644  {
645  .name = "default",
646  .type = AVMEDIA_TYPE_VIDEO,
647  },
648  { NULL }
649 };
650 
652  .name = "removegrain",
653  .description = NULL_IF_CONFIG_SMALL("Remove grain."),
654  .priv_size = sizeof(RemoveGrainContext),
656  .inputs = removegrain_inputs,
657  .outputs = removegrain_outputs,
658  .priv_class = &removegrain_class,
660 };
#define NULL
Definition: coverity.c:32
const char const char void * val
Definition: avisynth_c.h:634
void ff_removegrain_init_x86(RemoveGrainContext *rg)
const char * s
Definition: avisynth_c.h:631
AVFrame * out
const AVPixFmtDescriptor * av_pix_fmt_desc_get(enum AVPixelFormat pix_fmt)
Definition: pixdesc.c:2129
This structure describes decoded (raw) audio or video data.
Definition: frame.h:171
#define mi
static int mode08(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
AVOption.
Definition: opt.h:255
#define ma
planar YUV 4:4:4, 24bpp, (1 Cr & Cb sample per 1x1 Y samples)
Definition: pixfmt.h:68
misc image utilities
static const AVFilterPad outputs[]
Definition: af_ashowinfo.c:248
int av_pix_fmt_count_planes(enum AVPixelFormat pix_fmt)
Definition: pixdesc.c:2169
Main libavfilter public API header.
planar GBR 4:4:4 24bpp
Definition: pixfmt.h:188
static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
AVFilter ff_vf_removegrain
#define a1
Definition: regdef.h:47
#define FLAGS
int(* rg[4])(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
Definition: removegrain.h:35
AVFrame * ff_get_video_buffer(AVFilterLink *link, int w, int h)
Request a picture buffer with a specific set of permissions.
Definition: video.c:109
uint8_t log2_chroma_w
Amount to shift the luma width right to find the chroma width.
Definition: pixdesc.h:80
static int mode24(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
AVFilterFormats * ff_make_format_list(const int *fmts)
Create a list of supported formats.
Definition: formats.c:283
struct AVFilterGraph * graph
filtergraph this filter belongs to
Definition: avfilter.h:656
#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:451
const char * name
Pad name.
Definition: internal.h:69
int ff_filter_frame(AVFilterLink *link, AVFrame *frame)
Send a frame of data to the next filter.
Definition: avfilter.c:1158
#define a3
Definition: regdef.h:49
AVFrame * in
planar YUV 4:2:0, 20bpp, (1 Cr & Cb sample per 2x2 Y & A samples)
Definition: pixfmt.h:103
static const AVOption removegrain_options[]
uint8_t
mode
Definition: f_perms.c:27
AVOptions.
static int mode19(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
static av_cold int end(AVCodecContext *avctx)
Definition: avrndec.c:90
int plane
Definition: vf_blend.c:87
planar YUV 4:4:0 full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV440P and setting color_range...
Definition: pixfmt.h:102
int nb_threads
Maximum number of threads used by filters in this graph.
Definition: avfilter.h:1207
static const uint64_t c1
Definition: murmur3.c:49
planar YUV 4:2:2, 16bpp, full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV422P and setting col...
Definition: pixfmt.h:76
static int mode23(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
static int mode04(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
static int query_formats(AVFilterContext *ctx)
A filter pad used for either input or output.
Definition: internal.h:63
planar YUV 4:2:2 24bpp, (1 Cr & Cb sample per 2x1 Y & A samples)
Definition: pixfmt.h:281
int ff_set_common_formats(AVFilterContext *ctx, AVFilterFormats *formats)
A helper for query_formats() which sets all links to the same list of formats.
Definition: formats.c:542
#define td
Definition: regdef.h:70
uint8_t log2_chroma_h
Amount to shift the luma height right to find the chroma height.
Definition: pixdesc.h:89
static int mode20(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
#define AVERROR(e)
Definition: error.h:43
void av_frame_free(AVFrame **frame)
Free the frame and any dynamically allocated objects in it, e.g.
Definition: frame.c:148
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
Definition: internal.h:175
void * priv
private data for use by the filter
Definition: avfilter.h:654
#define AVFILTER_FLAG_SLICE_THREADS
The filter supports multithreading by splitting frames into multiple parts and processing them concur...
Definition: avfilter.h:442
const char * arg
Definition: jacosubdec.c:66
static int mode1112(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
#define FFMAX(a, b)
Definition: common.h:79
planar YUV 4:2:2, 16bpp, (1 Cr & Cb sample per 2x1 Y samples)
Definition: pixfmt.h:67
static int mode1314(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
static int cmp_int(const void *p1, const void *p2)
#define FFMIN(a, b)
Definition: common.h:81
float y
planar YUV 4:2:0, 12bpp, full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV420P and setting col...
Definition: pixfmt.h:75
static int mode09(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
#define FF_CEIL_RSHIFT(a, b)
Definition: common.h:57
#define a2
Definition: regdef.h:48
#define REMOVE_GRAIN_SORT_AXIS
#define FFABS(a)
Absolute value, Note, INT_MIN / INT64_MIN result in undefined behavior as they are not representable ...
Definition: common.h:68
static int mode07(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
float u
static const AVFilterPad removegrain_inputs[]
static int mode01(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
#define a5
Definition: regdef.h:51
AVS_Value src
Definition: avisynth_c.h:482
static int config_input(AVFilterLink *inlink)
int linesize[AV_NUM_DATA_POINTERS]
For video, size in bytes of each picture line.
Definition: frame.h:199
planar YUV 4:4:4 32bpp, (1 Cr & Cb sample per 1x1 Y & A samples)
Definition: pixfmt.h:280
Descriptor that unambiguously describes how the bits of a pixel are stored in the up to 4 data planes...
Definition: pixdesc.h:69
static int mode21(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){if(!*ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) returnNULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method!=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt)>2){ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc){av_free(ac);returnNULL;}returnac;}in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar){ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar?ac->channels:1;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_log(ac->avr, AV_LOG_TRACE,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> in
BYTE int const BYTE int int int height
Definition: avisynth_c.h:676
planar YUV 4:1:0, 9bpp, (1 Cr & Cb sample per 4x4 Y samples)
Definition: pixfmt.h:69
Filter definition.
Definition: avfilter.h:470
static const AVFilterPad inputs[]
Definition: af_ashowinfo.c:239
static const AVFilterPad removegrain_outputs[]
const char * name
Filter name.
Definition: avfilter.h:474
AVFilterLink ** outputs
array of pointers to output links
Definition: avfilter.h:648
static enum AVPixelFormat pix_fmts[]
Definition: libkvazaar.c:209
static int mode06(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
AVFilterInternal * internal
An opaque struct for libavfilter internal use.
Definition: avfilter.h:679
static int flags
Definition: cpu.c:47
static int mode18(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
uint8_t * data[AV_NUM_DATA_POINTERS]
pointer to the picture/channel planes.
Definition: frame.h:182
static int mode10(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
static int op(uint8_t **dst, const uint8_t *dst_end, GetByteContext *gb, int pixel, int count, int *x, int width, int linesize)
Perform decode operation.
Definition: anm.c:78
planar YUV 4:2:0, 12bpp, (1 Cr & Cb sample per 2x2 Y samples)
Definition: pixfmt.h:63
static int mode17(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
Y , 8bpp.
Definition: pixfmt.h:71
static int mode03(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
planar GBRA 4:4:4:4 32bpp
Definition: pixfmt.h:299
static double c[64]
planar YUV 4:4:4, 24bpp, full scale (JPEG), deprecated in favor of AV_PIX_FMT_YUV444P and setting col...
Definition: pixfmt.h:77
static const uint64_t c2
Definition: murmur3.c:50
planar YUV 4:1:1, 12bpp, (1 Cr & Cb sample per 4x1 Y samples)
Definition: pixfmt.h:70
avfilter_execute_func * execute
Definition: internal.h:164
#define a4
Definition: regdef.h:50
#define OFFSET(x)
static int mode1516(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
A list of supported formats for one end of a filter link.
Definition: formats.h:64
planar YUV 4:1:1, 12bpp, (1 Cr & Cb sample per 4x1 Y samples) full scale (JPEG), deprecated in favor ...
Definition: pixfmt.h:302
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){if(!*ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) returnNULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method!=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt)>2){ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc){av_free(ac);returnNULL;}returnac;}in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar){ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar?ac->channels:1;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_log(ac->avr, AV_LOG_TRACE,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> out
An instance of a filter.
Definition: avfilter.h:633
static int mode22(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
planar YUV 4:4:0 (1 Cr & Cb sample per 1x2 Y samples)
Definition: pixfmt.h:101
void INT64 start
Definition: avisynth_c.h:553
static int mode02(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
AVFILTER_DEFINE_CLASS(removegrain)
static int filter_frame(AVFilterLink *inlink, AVFrame *in)
void av_image_copy_plane(uint8_t *dst, int dst_linesize, const uint8_t *src, int src_linesize, int bytewidth, int height)
Copy image plane from src to dst.
Definition: imgutils.c:273
internal API functions
static int mode05(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
AVPixelFormat
Pixel format.
Definition: pixfmt.h:61
int av_frame_copy_props(AVFrame *dst, const AVFrame *src)
Copy only "metadata" fields from src to dst.
Definition: frame.c:553
void(* fl[4])(uint8_t *dst, uint8_t *src, ptrdiff_t stride, int pixels)
Definition: removegrain.h:37