You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

661 lines
20KB

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