LCOV - code coverage report
Current view: top level - libavfilter - vf_removegrain.c (source / functions) Hit Total Coverage
Test: coverage.info Lines: 355 363 97.8 %
Date: 2017-12-10 21:22:29 Functions: 26 26 100.0 %

          Line data    Source code
       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 "libavutil/qsort.h"
      28             : #include "avfilter.h"
      29             : #include "formats.h"
      30             : #include "internal.h"
      31             : #include "removegrain.h"
      32             : #include "video.h"
      33             : 
      34             : #define OFFSET(x) offsetof(RemoveGrainContext, x)
      35             : #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM
      36             : 
      37             : static const AVOption removegrain_options[] = {
      38             :     { "m0", "set mode for 1st plane", OFFSET(mode[0]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
      39             :     { "m1", "set mode for 2nd plane", OFFSET(mode[1]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
      40             :     { "m2", "set mode for 3rd plane", OFFSET(mode[2]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
      41             :     { "m3", "set mode for 4th plane", OFFSET(mode[3]), AV_OPT_TYPE_INT, {.i64=0}, 0, 24, FLAGS },
      42             :     {NULL}
      43             : };
      44             : 
      45             : AVFILTER_DEFINE_CLASS(removegrain);
      46             : 
      47          25 : static int query_formats(AVFilterContext *ctx)
      48             : {
      49             :     static const enum AVPixelFormat pix_fmts[] = {
      50             :         AV_PIX_FMT_GRAY8,
      51             :         AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUV444P,
      52             :         AV_PIX_FMT_YUVA420P, AV_PIX_FMT_YUVA422P, AV_PIX_FMT_YUVA444P,
      53             :         AV_PIX_FMT_YUV440P, AV_PIX_FMT_YUV411P, AV_PIX_FMT_YUV410P,
      54             :         AV_PIX_FMT_YUVJ411P, AV_PIX_FMT_YUVJ420P, AV_PIX_FMT_YUVJ422P,
      55             :         AV_PIX_FMT_YUVJ440P, AV_PIX_FMT_YUVJ444P,
      56             :         AV_PIX_FMT_GBRP, AV_PIX_FMT_GBRAP,
      57             :         AV_PIX_FMT_NONE
      58             :     };
      59             : 
      60          25 :     AVFilterFormats *fmts_list = ff_make_format_list(pix_fmts);
      61          25 :     if (!fmts_list)
      62           0 :         return AVERROR(ENOMEM);
      63          25 :     return ff_set_common_formats(ctx, fmts_list);
      64             : }
      65             : 
      66             : #define REMOVE_GRAIN_SORT_AXIS       \
      67             :     const int ma1 = FFMAX(a1, a8);   \
      68             :     const int mi1 = FFMIN(a1, a8);   \
      69             :     const int ma2 = FFMAX(a2, a7);   \
      70             :     const int mi2 = FFMIN(a2, a7);   \
      71             :     const int ma3 = FFMAX(a3, a6);   \
      72             :     const int mi3 = FFMIN(a3, a6);   \
      73             :     const int ma4 = FFMAX(a4, a5);   \
      74             :     const int mi4 = FFMIN(a4, a5);
      75             : 
      76      149516 : static int mode01(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
      77             : {
      78      149516 :     const int mi = FFMIN(FFMIN(FFMIN(a1, a2), FFMIN(a3, a4)), FFMIN(FFMIN(a5, a6), FFMIN(a7, a8)));
      79      149516 :     const int ma = FFMAX(FFMAX(FFMAX(a1, a2), FFMAX(a3, a4)), FFMAX(FFMAX(a5, a6), FFMAX(a7, a8)));
      80             : 
      81      149516 :     return av_clip(c, mi, ma);
      82             : }
      83             : 
      84     8093199 : static int cmp_int(const void *p1, const void *p2)
      85             : {
      86     8093199 :     int left  = *(const int *)p1;
      87     8093199 :     int right = *(const int *)p2;
      88     8093199 :     return FFDIFFSIGN(left, right);
      89             : }
      90             : 
      91      149516 : static int mode02(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
      92             : {
      93      149516 :     int a[8] = { a1, a2, a3, a4, a5, a6, a7, a8 };
      94             : 
      95      149516 :     AV_QSORT(a, 8, int, cmp_int);
      96             : 
      97      149516 :     return av_clip(c, a[2 - 1 ], a[7 - 1]);
      98             : }
      99             : 
     100      149516 : static int mode03(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     101             : {
     102      149516 :     int a[8] = { a1, a2, a3, a4, a5, a6, a7, a8 };
     103             : 
     104      149516 :     AV_QSORT(a, 8, int, cmp_int);
     105             : 
     106      149516 :     return av_clip(c, a[3 - 1 ], a[6 - 1]);
     107             : }
     108             : 
     109      149516 : static int mode04(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     110             : {
     111      149516 :     int a[8] = { a1, a2, a3, a4, a5, a6, a7, a8 };
     112             : 
     113      149516 :     AV_QSORT(a, 8, int, cmp_int);
     114             : 
     115      149516 :     return av_clip(c, a[4 - 1 ], a[5 - 1]);
     116             : }
     117             : 
     118      149516 : static int mode05(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     119             : {
     120      149516 :     REMOVE_GRAIN_SORT_AXIS
     121             : 
     122      149516 :     const int c1 = FFABS(c - av_clip(c, mi1, ma1));
     123      149516 :     const int c2 = FFABS(c - av_clip(c, mi2, ma2));
     124      149516 :     const int c3 = FFABS(c - av_clip(c, mi3, ma3));
     125      149516 :     const int c4 = FFABS(c - av_clip(c, mi4, ma4));
     126             : 
     127      149516 :     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      149516 :     if (mindiff == c4) {
     131      116368 :         return av_clip(c, mi4, ma4);
     132       33148 :     } else if (mindiff == c2) {
     133       14995 :         return av_clip(c, mi2, ma2);
     134       18153 :     } else if (mindiff == c3) {
     135       14086 :         return av_clip(c, mi3, ma3);
     136             :     }
     137             : 
     138        4067 :     return av_clip(c, mi1, ma1);
     139             : }
     140             : 
     141      149516 : static int mode06(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     142             : {
     143      149516 :     REMOVE_GRAIN_SORT_AXIS
     144             : 
     145      149516 :     const int d1 = ma1 - mi1;
     146      149516 :     const int d2 = ma2 - mi2;
     147      149516 :     const int d3 = ma3 - mi3;
     148      149516 :     const int d4 = ma4 - mi4;
     149             : 
     150      149516 :     const int cli1 = av_clip(c, mi1, ma1);
     151      149516 :     const int cli2 = av_clip(c, mi2, ma2);
     152      149516 :     const int cli3 = av_clip(c, mi3, ma3);
     153      149516 :     const int cli4 = av_clip(c, mi4, ma4);
     154             : 
     155      149516 :     const int c1 = av_clip_uint16((FFABS(c - cli1) << 1) + d1);
     156      149516 :     const int c2 = av_clip_uint16((FFABS(c - cli2) << 1) + d2);
     157      149516 :     const int c3 = av_clip_uint16((FFABS(c - cli3) << 1) + d3);
     158      149516 :     const int c4 = av_clip_uint16((FFABS(c - cli4) << 1) + d4);
     159             : 
     160      149516 :     const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
     161             : 
     162      149516 :     if (mindiff == c4) {
     163       30362 :         return cli4;
     164      119154 :     } else if (mindiff == c2) {
     165       26836 :         return cli2;
     166       92318 :     } else if (mindiff == c3) {
     167       85122 :         return cli3;
     168             :     }
     169             : 
     170        7196 :     return cli1;
     171             : }
     172             : 
     173      149516 : static int mode07(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     174             : {
     175      149516 :     REMOVE_GRAIN_SORT_AXIS
     176             : 
     177      149516 :     const int d1 = ma1 - mi1;
     178      149516 :     const int d2 = ma2 - mi2;
     179      149516 :     const int d3 = ma3 - mi3;
     180      149516 :     const int d4 = ma4 - mi4;
     181             : 
     182      149516 :     const int cli1 = av_clip(c, mi1, ma1);
     183      149516 :     const int cli2 = av_clip(c, mi2, ma2);
     184      149516 :     const int cli3 = av_clip(c, mi3, ma3);
     185      149516 :     const int cli4 = av_clip(c, mi4, ma4);
     186             : 
     187      149516 :     const int c1 = FFABS(c - cli1) + d1;
     188      149516 :     const int c2 = FFABS(c - cli2) + d2;
     189      149516 :     const int c3 = FFABS(c - cli3) + d3;
     190      149516 :     const int c4 = FFABS(c - cli4) + d4;
     191             : 
     192      149516 :     const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
     193             : 
     194      149516 :     if (mindiff == c4) {
     195       30479 :         return cli4;
     196      119037 :     } else if (mindiff == c2) {
     197       27791 :         return cli2;
     198       91246 :     } else if (mindiff == c3) {
     199       84095 :         return cli3;
     200             :     }
     201             : 
     202        7151 :     return cli1;
     203             : }
     204             : 
     205      149516 : static int mode08(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     206             : {
     207      149516 :     REMOVE_GRAIN_SORT_AXIS
     208             : 
     209      149516 :     const int d1 = ma1 - mi1;
     210      149516 :     const int d2 = ma2 - mi2;
     211      149516 :     const int d3 = ma3 - mi3;
     212      149516 :     const int d4 = ma4 - mi4;
     213             : 
     214      149516 :     const int cli1 = av_clip(c, mi1, ma1);
     215      149516 :     const int cli2 = av_clip(c, mi2, ma2);
     216      149516 :     const int cli3 = av_clip(c, mi3, ma3);
     217      149516 :     const int cli4 = av_clip(c, mi4, ma4);
     218             : 
     219      149516 :     const int c1 = av_clip_uint16(FFABS(c - cli1) + (d1 << 1));
     220      149516 :     const int c2 = av_clip_uint16(FFABS(c - cli2) + (d2 << 1));
     221      149516 :     const int c3 = av_clip_uint16(FFABS(c - cli3) + (d3 << 1));
     222      149516 :     const int c4 = av_clip_uint16(FFABS(c - cli4) + (d4 << 1));
     223             : 
     224      149516 :     const int mindiff = FFMIN(FFMIN(c1, c2), FFMIN(c3, c4));
     225             : 
     226      149516 :     if (mindiff == c4) {
     227       31274 :         return cli4;
     228      118242 :     } else if (mindiff == c2) {
     229       25063 :         return cli2;
     230       93179 :     } else if (mindiff == c3) {
     231       83083 :         return cli3;
     232             :     }
     233             : 
     234       10096 :     return cli1;
     235             : }
     236             : 
     237      149516 : static int mode09(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     238             : {
     239      149516 :     REMOVE_GRAIN_SORT_AXIS
     240             : 
     241      149516 :     const int d1 = ma1 - mi1;
     242      149516 :     const int d2 = ma2 - mi2;
     243      149516 :     const int d3 = ma3 - mi3;
     244      149516 :     const int d4 = ma4 - mi4;
     245             : 
     246      149516 :     const int mindiff = FFMIN(FFMIN(d1, d2), FFMIN(d3, d4));
     247             : 
     248      149516 :     if (mindiff == d4) {
     249       31449 :         return av_clip(c, mi4, ma4);
     250      118067 :     } else if (mindiff == d2) {
     251       25061 :         return av_clip(c, mi2, ma2);
     252       93006 :     } else if (mindiff == d3) {
     253       82804 :         return av_clip(c, mi3, ma3);
     254             :     }
     255             : 
     256       10202 :     return av_clip(c, mi1, ma1);
     257             : }
     258             : 
     259      149516 : static int mode10(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     260             : {
     261      149516 :     const int d1 = FFABS(c - a1);
     262      149516 :     const int d2 = FFABS(c - a2);
     263      149516 :     const int d3 = FFABS(c - a3);
     264      149516 :     const int d4 = FFABS(c - a4);
     265      149516 :     const int d5 = FFABS(c - a5);
     266      149516 :     const int d6 = FFABS(c - a6);
     267      149516 :     const int d7 = FFABS(c - a7);
     268      149516 :     const int d8 = FFABS(c - a8);
     269             : 
     270      149516 :     const int mindiff = FFMIN(FFMIN(FFMIN(d1, d2), FFMIN(d3, d4)),
     271             :                               FFMIN(FFMIN(d5, d6), FFMIN(d7, d8)));
     272             : 
     273      149516 :     if (mindiff == d7) return a7;
     274      135259 :     if (mindiff == d8) return a8;
     275      130631 :     if (mindiff == d6) return a6;
     276       82359 :     if (mindiff == d2) return a2;
     277       69646 :     if (mindiff == d3) return a3;
     278       28337 :     if (mindiff == d1) return a1;
     279       24786 :     if (mindiff == d5) return a5;
     280             : 
     281        5870 :     return a4;
     282             : }
     283             : 
     284      299032 : static int mode1112(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     285             : {
     286      299032 :     const int sum = 4 * c + 2 * (a2 + a4 + a5 + a7) + a1 + a3 + a6 + a8;
     287      299032 :     const int val = (sum + 8) >> 4;
     288             : 
     289      299032 :     return val;
     290             : }
     291             : 
     292      149516 : static int mode1314(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     293             : {
     294      149516 :     const int d1 = FFABS(a1 - a8);
     295      149516 :     const int d2 = FFABS(a2 - a7);
     296      149516 :     const int d3 = FFABS(a3 - a6);
     297             : 
     298      149516 :     const int mindiff = FFMIN(FFMIN(d1, d2), d3);
     299             : 
     300      149516 :     if (mindiff == d2) {
     301       29535 :         return (a2 + a7 + 1) >> 1;
     302             :     }
     303      119981 :     if (mindiff == d3) {
     304       99814 :         return (a3 + a6 + 1) >> 1;
     305             :     }
     306             : 
     307       20167 :     return (a1 + a8 + 1) >> 1;
     308             : }
     309             : 
     310      149516 : static int mode1516(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     311             : {
     312      149516 :     const int d1 = FFABS(a1 - a8);
     313      149516 :     const int d2 = FFABS(a2 - a7);
     314      149516 :     const int d3 = FFABS(a3 - a6);
     315             : 
     316      149516 :     const int mindiff = FFMIN(FFMIN(d1, d2), d3);
     317      149516 :     const int average = (2 * (a2 + a7) + a1 + a3 + a6 + a8 + 4) >> 3;
     318             : 
     319      149516 :     if (mindiff == d2) {
     320       29535 :         return av_clip(average, FFMIN(a2, a7), FFMAX(a2, a7));
     321             :     }
     322      119981 :     if (mindiff == d3) {
     323       99814 :         return av_clip(average, FFMIN(a3, a6), FFMAX(a3, a6));
     324             :     }
     325             : 
     326       20167 :     return av_clip(average, FFMIN(a1, a8), FFMAX(a1, a8));
     327             : }
     328             : 
     329      149516 : static int mode17(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     330             : {
     331      149516 :     REMOVE_GRAIN_SORT_AXIS
     332             : 
     333      149516 :     const int l = FFMAX(FFMAX(mi1, mi2), FFMAX(mi3, mi4));
     334      149516 :     const int u = FFMIN(FFMIN(ma1, ma2), FFMIN(ma3, ma4));
     335             : 
     336      149516 :     return av_clip(c, FFMIN(l, u), FFMAX(l, u));
     337             : }
     338             : 
     339      149516 : static int mode18(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     340             : {
     341      149516 :     const int d1 = FFMAX(FFABS(c - a1), FFABS(c - a8));
     342      149516 :     const int d2 = FFMAX(FFABS(c - a2), FFABS(c - a7));
     343      149516 :     const int d3 = FFMAX(FFABS(c - a3), FFABS(c - a6));
     344      149516 :     const int d4 = FFMAX(FFABS(c - a4), FFABS(c - a5));
     345             : 
     346      149516 :     const int mindiff = FFMIN(FFMIN(d1, d2), FFMIN(d3, d4));
     347             : 
     348      149516 :     if (mindiff == d4) {
     349       30555 :         return av_clip(c, FFMIN(a4, a5), FFMAX(a4, a5));
     350             :     }
     351      118961 :     if (mindiff == d2) {
     352       26866 :         return av_clip(c, FFMIN(a2, a7), FFMAX(a2, a7));
     353             :     }
     354       92095 :     if (mindiff == d3) {
     355       85057 :         return av_clip(c, FFMIN(a3, a6), FFMAX(a3, a6));
     356             :     }
     357             : 
     358        7038 :     return av_clip(c, FFMIN(a1, a8), FFMAX(a1, a8));
     359             : }
     360             : 
     361      149516 : static int mode19(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     362             : {
     363      149516 :     const int sum = a1 + a2 + a3 + a4 + a5 + a6 + a7 + a8;
     364      149516 :     const int val = (sum + 4) >> 3;
     365             : 
     366      149516 :     return val;
     367             : }
     368             : 
     369      149516 : static int mode20(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     370             : {
     371      149516 :     const int sum = a1 + a2 + a3 + a4 + c + a5 + a6 + a7 + a8;
     372      149516 :     const int val = (sum + 4) / 9;
     373             : 
     374      149516 :     return val;
     375             : }
     376             : 
     377      149516 : static int mode21(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     378             : {
     379      149516 :     const int l1l = (a1 + a8) >> 1;
     380      149516 :     const int l2l = (a2 + a7) >> 1;
     381      149516 :     const int l3l = (a3 + a6) >> 1;
     382      149516 :     const int l4l = (a4 + a5) >> 1;
     383             : 
     384      149516 :     const int l1h = (a1 + a8 + 1) >> 1;
     385      149516 :     const int l2h = (a2 + a7 + 1) >> 1;
     386      149516 :     const int l3h = (a3 + a6 + 1) >> 1;
     387      149516 :     const int l4h = (a4 + a5 + 1) >> 1;
     388             : 
     389      149516 :     const int mi = FFMIN(FFMIN(l1l, l2l), FFMIN(l3l, l4l));
     390      149516 :     const int ma = FFMAX(FFMAX(l1h, l2h), FFMAX(l3h, l4h));
     391             : 
     392      149516 :     return av_clip(c, mi, ma);
     393             : }
     394             : 
     395      149516 : static int mode22(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     396             : {
     397      149516 :     const int l1 = (a1 + a8 + 1) >> 1;
     398      149516 :     const int l2 = (a2 + a7 + 1) >> 1;
     399      149516 :     const int l3 = (a3 + a6 + 1) >> 1;
     400      149516 :     const int l4 = (a4 + a5 + 1) >> 1;
     401             : 
     402      149516 :     const int mi = FFMIN(FFMIN(l1, l2), FFMIN(l3, l4));
     403      149516 :     const int ma = FFMAX(FFMAX(l1, l2), FFMAX(l3, l4));
     404             : 
     405      149516 :     return av_clip(c, mi, ma);
     406             : }
     407             : 
     408      149516 : static int mode23(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     409             : {
     410      149516 :     REMOVE_GRAIN_SORT_AXIS
     411             : 
     412      149516 :     const int linediff1 = ma1 - mi1;
     413      149516 :     const int linediff2 = ma2 - mi2;
     414      149516 :     const int linediff3 = ma3 - mi3;
     415      149516 :     const int linediff4 = ma4 - mi4;
     416             : 
     417      149516 :     const int u1 = FFMIN(c - ma1, linediff1);
     418      149516 :     const int u2 = FFMIN(c - ma2, linediff2);
     419      149516 :     const int u3 = FFMIN(c - ma3, linediff3);
     420      149516 :     const int u4 = FFMIN(c - ma4, linediff4);
     421      149516 :     const int u = FFMAX(FFMAX(FFMAX(u1, u2), FFMAX(u3, u4)), 0);
     422             : 
     423      149516 :     const int d1 = FFMIN(mi1 - c, linediff1);
     424      149516 :     const int d2 = FFMIN(mi2 - c, linediff2);
     425      149516 :     const int d3 = FFMIN(mi3 - c, linediff3);
     426      149516 :     const int d4 = FFMIN(mi4 - c, linediff4);
     427      149516 :     const int d = FFMAX(FFMAX(FFMAX(d1, d2), FFMAX(d3, d4)), 0);
     428             : 
     429      149516 :     return c - u + d;  // This probably will never overflow.
     430             : }
     431             : 
     432      149516 : static int mode24(int c, int a1, int a2, int a3, int a4, int a5, int a6, int a7, int a8)
     433             : {
     434      149516 :     REMOVE_GRAIN_SORT_AXIS
     435             : 
     436      149516 :     const int linediff1 = ma1 - mi1;
     437      149516 :     const int linediff2 = ma2 - mi2;
     438      149516 :     const int linediff3 = ma3 - mi3;
     439      149516 :     const int linediff4 = ma4 - mi4;
     440             : 
     441      149516 :     const int tu1 = c - ma1;
     442      149516 :     const int tu2 = c - ma2;
     443      149516 :     const int tu3 = c - ma3;
     444      149516 :     const int tu4 = c - ma4;
     445             : 
     446      149516 :     const int u1 = FFMIN(tu1, linediff1 - tu1);
     447      149516 :     const int u2 = FFMIN(tu2, linediff2 - tu2);
     448      149516 :     const int u3 = FFMIN(tu3, linediff3 - tu3);
     449      149516 :     const int u4 = FFMIN(tu4, linediff4 - tu4);
     450      149516 :     const int u = FFMAX(FFMAX(FFMAX(u1, u2), FFMAX(u3, u4)), 0);
     451             : 
     452      149516 :     const int td1 = mi1 - c;
     453      149516 :     const int td2 = mi2 - c;
     454      149516 :     const int td3 = mi3 - c;
     455      149516 :     const int td4 = mi4 - c;
     456             : 
     457      149516 :     const int d1 = FFMIN(td1, linediff1 - td1);
     458      149516 :     const int d2 = FFMIN(td2, linediff2 - td2);
     459      149516 :     const int d3 = FFMIN(td3, linediff3 - td3);
     460      149516 :     const int d4 = FFMIN(td4, linediff4 - td4);
     461      149516 :     const int d = FFMAX(FFMAX(FFMAX(d1, d2), FFMAX(d3, d4)), 0);
     462             : 
     463      149516 :     return c - u + d;  // This probably will never overflow.
     464             : }
     465             : 
     466          25 : static int config_input(AVFilterLink *inlink)
     467             : {
     468          25 :     RemoveGrainContext *s = inlink->dst->priv;
     469          25 :     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
     470             :     int i;
     471             : 
     472          25 :     s->nb_planes = av_pix_fmt_count_planes(inlink->format);
     473             : 
     474          25 :     s->planeheight[1] = s->planeheight[2] = AV_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h);
     475          25 :     s->planeheight[0] = s->planeheight[3] = inlink->h;
     476          25 :     s->planewidth[1]  = s->planewidth[2]  = AV_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w);
     477          25 :     s->planewidth[0]  = s->planewidth[3]  = inlink->w;
     478             : 
     479         100 :     for (i = 0; i < s->nb_planes; i++) {
     480          75 :         switch (s->mode[i]) {
     481           3 :         case 1:  s->rg[i] = mode01;   break;
     482           3 :         case 2:  s->rg[i] = mode02;   break;
     483           3 :         case 3:  s->rg[i] = mode03;   break;
     484           3 :         case 4:  s->rg[i] = mode04;   break;
     485           3 :         case 5:  s->rg[i] = mode05;   break;
     486           3 :         case 6:  s->rg[i] = mode06;   break;
     487           3 :         case 7:  s->rg[i] = mode07;   break;
     488           3 :         case 8:  s->rg[i] = mode08;   break;
     489           3 :         case 9:  s->rg[i] = mode09;   break;
     490           3 :         case 10: s->rg[i] = mode10;   break;
     491           3 :         case 11: s->rg[i] = mode1112; break;
     492           3 :         case 12: s->rg[i] = mode1112; break;
     493           3 :         case 13: s->skip_odd = 1;
     494           3 :                  s->rg[i] = mode1314; break;
     495           3 :         case 14: s->skip_even = 1;
     496           3 :                  s->rg[i] = mode1314; break;
     497           3 :         case 15: s->skip_odd = 1;
     498           3 :                  s->rg[i] = mode1516; break;
     499           3 :         case 16: s->skip_even = 1;
     500           3 :                  s->rg[i] = mode1516; break;
     501           3 :         case 17: s->rg[i] = mode17;   break;
     502           3 :         case 18: s->rg[i] = mode18;   break;
     503           3 :         case 19: s->rg[i] = mode19;   break;
     504           3 :         case 20: s->rg[i] = mode20;   break;
     505           3 :         case 21: s->rg[i] = mode21;   break;
     506           3 :         case 22: s->rg[i] = mode22;   break;
     507           3 :         case 23: s->rg[i] = mode23;   break;
     508           3 :         case 24: s->rg[i] = mode24;   break;
     509             :         }
     510             :     }
     511             : 
     512             :     if (ARCH_X86)
     513          25 :         ff_removegrain_init_x86(s);
     514             : 
     515          25 :     return 0;
     516             : }
     517             : 
     518             : typedef struct ThreadData {
     519             :     AVFrame *in, *out;
     520             :     int plane;
     521             : } ThreadData;
     522             : 
     523         648 : static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
     524             : {
     525         648 :     RemoveGrainContext *s = ctx->priv;
     526         648 :     ThreadData *td = arg;
     527         648 :     AVFrame *in = td->in;
     528         648 :     AVFrame *out = td->out;
     529         648 :     const int i = td->plane;
     530         648 :     const int height = s->planeheight[i];
     531         648 :     const int om = in->linesize[i] - 1;
     532         648 :     const int o0 = in->linesize[i]    ;
     533         648 :     const int op = in->linesize[i] + 1;
     534         648 :     int start = (height *  jobnr   ) / nb_jobs;
     535         648 :     int end   = (height * (jobnr+1)) / nb_jobs;
     536             :     int x, y;
     537             : 
     538         648 :     start = FFMAX(1, start);
     539         648 :     end   = FFMIN(height-1, end);
     540       14328 :     for (y = start; y < end; y++) {
     541       13680 :         uint8_t *dst = out->data[i];
     542       13680 :         uint8_t *src = in->data[i];
     543             : 
     544       13680 :         src = in->data[i] + y * in->linesize[i];
     545       13680 :         dst = out->data[i] + y * out->linesize[i];
     546             : 
     547       13680 :         if (s->skip_even && !(y & 1)) {
     548         570 :             memcpy(dst, src, s->planewidth[i]);
     549         570 :             continue;
     550             :         }
     551       13110 :         if (s->skip_odd && y & 1) {
     552         570 :             memcpy(dst, src, s->planewidth[i]);
     553         570 :             continue;
     554             :         }
     555             : 
     556       12540 :         *dst++ = *src++;
     557             : 
     558       12540 :         if (s->fl[i]) {
     559           0 :             int w_asm = (s->planewidth[i] - 2) & ~15;
     560             : 
     561           0 :             s->fl[i](dst, src, in->linesize[i], w_asm);
     562             : 
     563           0 :             x = 1 + w_asm;
     564           0 :             dst += w_asm;
     565           0 :             src += w_asm;
     566             :         } else
     567       12540 :             x = 1;
     568             : 
     569     3301892 :         for (; x < s->planewidth[i] - 1; x++) {
     570     3289352 :             const int a1 = src[-op];
     571     3289352 :             const int a2 = src[-o0];
     572     3289352 :             const int a3 = src[-om];
     573     3289352 :             const int a4 = src[-1 ];
     574     3289352 :             const int c  = src[ 0 ];
     575     3289352 :             const int a5 = src[ 1 ];
     576     3289352 :             const int a6 = src[ om];
     577     3289352 :             const int a7 = src[ o0];
     578     3289352 :             const int a8 = src[ op];
     579             : 
     580     3289352 :             const int res = s->rg[i](c, a1, a2, a3, a4, a5, a6, a7, a8);
     581             : 
     582     3289352 :             *dst = res;
     583     3289352 :             dst++, src++;
     584             :         }
     585       12540 :         dst[0] = src[0];
     586             :     }
     587             : 
     588         648 :     return 0;
     589             : }
     590             : 
     591          25 : static int filter_frame(AVFilterLink *inlink, AVFrame *in)
     592             : {
     593          25 :     AVFilterContext *ctx = inlink->dst;
     594          25 :     AVFilterLink *outlink = ctx->outputs[0];
     595          25 :     RemoveGrainContext *s = ctx->priv;
     596             :     ThreadData td;
     597             :     AVFrame *out;
     598             :     int i;
     599             : 
     600          25 :     out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
     601          25 :     if (!out) {
     602           0 :         av_frame_free(&in);
     603           0 :         return AVERROR(ENOMEM);
     604             :     }
     605          25 :     av_frame_copy_props(out, in);
     606             : 
     607         100 :     for (i = 0; i < s->nb_planes; i++) {
     608          75 :         uint8_t *dst = out->data[i];
     609          75 :         uint8_t *src = in->data[i];
     610             : 
     611          75 :         if (s->mode[i] == 0) {
     612           6 :             av_image_copy_plane(dst, out->linesize[i],
     613           3 :                                 src, in->linesize[i],
     614             :                                 s->planewidth[i], s->planeheight[i]);
     615           3 :             continue;
     616             :         }
     617             : 
     618          72 :         memcpy(dst, src, s->planewidth[i]);
     619             : 
     620          72 :         td.in = in; td.out = out; td.plane = i;
     621         144 :         ctx->internal->execute(ctx, filter_slice, &td, NULL,
     622          72 :                                FFMIN(s->planeheight[i], ff_filter_get_nb_threads(ctx)));
     623             : 
     624          72 :         src = in->data[i] + (s->planeheight[i] - 1) * in->linesize[i];
     625          72 :         dst = out->data[i] + (s->planeheight[i] - 1) * out->linesize[i];
     626          72 :         memcpy(dst, src, s->planewidth[i]);
     627             :     }
     628             : 
     629          25 :     av_frame_free(&in);
     630          25 :     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             : 
     643             : static const AVFilterPad removegrain_outputs[] = {
     644             :     {
     645             :         .name = "default",
     646             :         .type = AVMEDIA_TYPE_VIDEO,
     647             :     },
     648             :     { NULL }
     649             : };
     650             : 
     651             : AVFilter ff_vf_removegrain = {
     652             :     .name          = "removegrain",
     653             :     .description   = NULL_IF_CONFIG_SMALL("Remove grain."),
     654             :     .priv_size     = sizeof(RemoveGrainContext),
     655             :     .query_formats = query_formats,
     656             :     .inputs        = removegrain_inputs,
     657             :     .outputs       = removegrain_outputs,
     658             :     .priv_class    = &removegrain_class,
     659             :     .flags         = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | AVFILTER_FLAG_SLICE_THREADS,
     660             : };

Generated by: LCOV version 1.13