| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | /* | ||
| 2 | * Copyright (c) 2021 Paul B Mahol | ||
| 3 | * | ||
| 4 | * This file is part of FFmpeg. | ||
| 5 | * | ||
| 6 | * FFmpeg is free software; you can redistribute it and/or | ||
| 7 | * modify it under the terms of the GNU Lesser General Public | ||
| 8 | * License as published by the Free Software Foundation; either | ||
| 9 | * version 2.1 of the License, or (at your option) any later version. | ||
| 10 | * | ||
| 11 | * FFmpeg is distributed in the hope that it will be useful, | ||
| 12 | * but WITHOUT ANY WARRANTY; without even the implied warranty of | ||
| 13 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | ||
| 14 | * Lesser General Public License for more details. | ||
| 15 | * | ||
| 16 | * You should have received a copy of the GNU Lesser General Public | ||
| 17 | * License along with FFmpeg; if not, write to the Free Software | ||
| 18 | * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA | ||
| 19 | */ | ||
| 20 | |||
| 21 | #include "libavutil/common.h" | ||
| 22 | #include "libavutil/imgutils.h" | ||
| 23 | #include "libavutil/opt.h" | ||
| 24 | #include "libavutil/pixdesc.h" | ||
| 25 | #include "avfilter.h" | ||
| 26 | #include "filters.h" | ||
| 27 | #include "video.h" | ||
| 28 | |||
| 29 | typedef struct ESTDIFContext { | ||
| 30 | const AVClass *class; | ||
| 31 | |||
| 32 | int mode; ///< 0 is frame, 1 is field | ||
| 33 | int parity; ///< frame field parity | ||
| 34 | int deint; ///< which frames to deinterlace | ||
| 35 | int rslope; ///< best edge slope search radius | ||
| 36 | int redge; ///< best edge match search radius | ||
| 37 | int ecost; ///< edge cost for edge matching | ||
| 38 | int mcost; ///< middle cost for edge matching | ||
| 39 | int dcost; ///< distance cost for edge matching | ||
| 40 | int interp; ///< type of interpolation | ||
| 41 | int linesize[4]; ///< bytes of pixel data per line for each plane | ||
| 42 | int planewidth[4]; ///< width of each plane | ||
| 43 | int planeheight[4]; ///< height of each plane | ||
| 44 | int field; ///< which field are we on, 0 or 1 | ||
| 45 | int eof; | ||
| 46 | int depth; | ||
| 47 | int max; | ||
| 48 | int nb_planes; | ||
| 49 | int nb_threads; | ||
| 50 | AVFrame *prev; | ||
| 51 | |||
| 52 | void (*interpolate)(struct ESTDIFContext *s, uint8_t *dst, | ||
| 53 | const uint8_t *prev_line, const uint8_t *next_line, | ||
| 54 | const uint8_t *prev2_line, const uint8_t *next2_line, | ||
| 55 | const uint8_t *prev3_line, const uint8_t *next3_line, | ||
| 56 | int x, int width, int rslope, int redge, | ||
| 57 | int depth, int *K); | ||
| 58 | |||
| 59 | unsigned (*mid_8[3])(const uint8_t *const prev, | ||
| 60 | const uint8_t *const next, | ||
| 61 | const uint8_t *const prev2, | ||
| 62 | const uint8_t *const next2, | ||
| 63 | const uint8_t *const prev3, | ||
| 64 | const uint8_t *const next3, | ||
| 65 | int end, int x, int k, int depth); | ||
| 66 | |||
| 67 | unsigned (*mid_16[3])(const uint16_t *const prev, | ||
| 68 | const uint16_t *const next, | ||
| 69 | const uint16_t *const prev2, | ||
| 70 | const uint16_t *const next2, | ||
| 71 | const uint16_t *const prev3, | ||
| 72 | const uint16_t *const next3, | ||
| 73 | int end, int x, int k, int depth); | ||
| 74 | } ESTDIFContext; | ||
| 75 | |||
| 76 | #define MAX_R 15 | ||
| 77 | #define S (MAX_R * 2 + 1) | ||
| 78 | |||
| 79 | #define OFFSET(x) offsetof(ESTDIFContext, x) | ||
| 80 | #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_RUNTIME_PARAM | ||
| 81 | #define CONST(name, help, val, u) { name, help, 0, AV_OPT_TYPE_CONST, {.i64=val}, 0, 0, FLAGS, .unit = u } | ||
| 82 | |||
| 83 | static const AVOption estdif_options[] = { | ||
| 84 | { "mode", "specify the mode", OFFSET(mode), AV_OPT_TYPE_INT, {.i64=1}, 0, 1, FLAGS, .unit = "mode" }, | ||
| 85 | CONST("frame", "send one frame for each frame", 0, "mode"), | ||
| 86 | CONST("field", "send one frame for each field", 1, "mode"), | ||
| 87 | { "parity", "specify the assumed picture field parity", OFFSET(parity), AV_OPT_TYPE_INT, {.i64=-1}, -1, 1, FLAGS, .unit = "parity" }, | ||
| 88 | CONST("tff", "assume top field first", 0, "parity"), | ||
| 89 | CONST("bff", "assume bottom field first", 1, "parity"), | ||
| 90 | CONST("auto", "auto detect parity", -1, "parity"), | ||
| 91 | { "deint", "specify which frames to deinterlace", OFFSET(deint), AV_OPT_TYPE_INT, {.i64=0}, 0, 1, FLAGS, .unit = "deint" }, | ||
| 92 | CONST("all", "deinterlace all frames", 0, "deint"), | ||
| 93 | CONST("interlaced", "only deinterlace frames marked as interlaced", 1, "deint"), | ||
| 94 | { "rslope", "specify the search radius for edge slope tracing", OFFSET(rslope), AV_OPT_TYPE_INT, {.i64=1}, 1, MAX_R, FLAGS }, | ||
| 95 | { "redge", "specify the search radius for best edge matching", OFFSET(redge), AV_OPT_TYPE_INT, {.i64=2}, 0, MAX_R, FLAGS }, | ||
| 96 | { "ecost", "specify the edge cost for edge matching", OFFSET(ecost), AV_OPT_TYPE_INT, {.i64=2}, 0, 50, FLAGS }, | ||
| 97 | { "mcost", "specify the middle cost for edge matching", OFFSET(mcost), AV_OPT_TYPE_INT, {.i64=1}, 0, 50, FLAGS }, | ||
| 98 | { "dcost", "specify the distance cost for edge matching", OFFSET(dcost), AV_OPT_TYPE_INT, {.i64=1}, 0, 50, FLAGS }, | ||
| 99 | { "interp", "specify the type of interpolation", OFFSET(interp), AV_OPT_TYPE_INT, {.i64=1}, 0, 2, FLAGS, .unit = "interp" }, | ||
| 100 | CONST("2p", "two-point interpolation", 0, "interp"), | ||
| 101 | CONST("4p", "four-point interpolation", 1, "interp"), | ||
| 102 | CONST("6p", "six-point interpolation", 2, "interp"), | ||
| 103 | { NULL } | ||
| 104 | }; | ||
| 105 | |||
| 106 | AVFILTER_DEFINE_CLASS(estdif); | ||
| 107 | |||
| 108 | static const enum AVPixelFormat pix_fmts[] = { | ||
| 109 | AV_PIX_FMT_YUV410P, AV_PIX_FMT_YUV411P, | ||
| 110 | AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUV422P, | ||
| 111 | AV_PIX_FMT_YUV440P, AV_PIX_FMT_YUV444P, | ||
| 112 | AV_PIX_FMT_YUVJ444P, AV_PIX_FMT_YUVJ440P, | ||
| 113 | AV_PIX_FMT_YUVJ422P, AV_PIX_FMT_YUVJ420P, | ||
| 114 | AV_PIX_FMT_YUVJ411P, | ||
| 115 | AV_PIX_FMT_YUVA420P, AV_PIX_FMT_YUVA422P, AV_PIX_FMT_YUVA444P, | ||
| 116 | AV_PIX_FMT_GBRP, AV_PIX_FMT_GBRAP, | ||
| 117 | AV_PIX_FMT_GRAY8, | ||
| 118 | AV_PIX_FMT_GRAY9, AV_PIX_FMT_GRAY10, AV_PIX_FMT_GRAY12, AV_PIX_FMT_GRAY14, AV_PIX_FMT_GRAY16, | ||
| 119 | AV_PIX_FMT_YUV420P9, AV_PIX_FMT_YUV422P9, AV_PIX_FMT_YUV444P9, | ||
| 120 | AV_PIX_FMT_YUV420P10, AV_PIX_FMT_YUV422P10, AV_PIX_FMT_YUV444P10, | ||
| 121 | AV_PIX_FMT_YUV440P10, | ||
| 122 | AV_PIX_FMT_YUV420P12, AV_PIX_FMT_YUV422P12, AV_PIX_FMT_YUV444P12, | ||
| 123 | AV_PIX_FMT_YUV440P12, | ||
| 124 | AV_PIX_FMT_YUV420P14, AV_PIX_FMT_YUV422P14, AV_PIX_FMT_YUV444P14, | ||
| 125 | AV_PIX_FMT_YUV420P16, AV_PIX_FMT_YUV422P16, AV_PIX_FMT_YUV444P16, | ||
| 126 | AV_PIX_FMT_GBRP9, AV_PIX_FMT_GBRP10, AV_PIX_FMT_GBRP12, AV_PIX_FMT_GBRP14, AV_PIX_FMT_GBRP16, | ||
| 127 | AV_PIX_FMT_YUVA444P9, AV_PIX_FMT_YUVA444P10, AV_PIX_FMT_YUVA444P12, AV_PIX_FMT_YUVA444P16, | ||
| 128 | AV_PIX_FMT_YUVA422P9, AV_PIX_FMT_YUVA422P10, AV_PIX_FMT_YUVA422P12, AV_PIX_FMT_YUVA422P16, | ||
| 129 | AV_PIX_FMT_YUVA420P9, AV_PIX_FMT_YUVA420P10, AV_PIX_FMT_YUVA420P16, | ||
| 130 | AV_PIX_FMT_GBRAP10, AV_PIX_FMT_GBRAP12, AV_PIX_FMT_GBRAP16, | ||
| 131 | AV_PIX_FMT_NONE | ||
| 132 | }; | ||
| 133 | |||
| 134 | ✗ | static int config_output(AVFilterLink *outlink) | |
| 135 | { | ||
| 136 | ✗ | FilterLink *outl = ff_filter_link(outlink); | |
| 137 | ✗ | AVFilterContext *ctx = outlink->src; | |
| 138 | ✗ | AVFilterLink *inlink = ctx->inputs[0]; | |
| 139 | ✗ | FilterLink *inl = ff_filter_link(inlink); | |
| 140 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 141 | |||
| 142 | ✗ | outlink->time_base = av_mul_q(inlink->time_base, (AVRational){1, 2}); | |
| 143 | ✗ | if (s->mode) | |
| 144 | ✗ | outl->frame_rate = av_mul_q(inl->frame_rate, (AVRational){2, 1}); | |
| 145 | |||
| 146 | ✗ | return 0; | |
| 147 | } | ||
| 148 | |||
| 149 | typedef struct ThreadData { | ||
| 150 | AVFrame *out, *in; | ||
| 151 | } ThreadData; | ||
| 152 | |||
| 153 | #define MIDL(type, ss) \ | ||
| 154 | static unsigned midl_##ss(const type *const prev, \ | ||
| 155 | const type *const next, \ | ||
| 156 | int end, int x, int k) \ | ||
| 157 | { \ | ||
| 158 | return (prev[av_clip(x + k, 0, end)] + \ | ||
| 159 | next[av_clip(x - k, 0, end)] + 1) >> 1; \ | ||
| 160 | } | ||
| 161 | |||
| 162 | ✗ | MIDL(uint8_t, 8) | |
| 163 | ✗ | MIDL(uint16_t, 16) | |
| 164 | |||
| 165 | #define MID2(type, ss) \ | ||
| 166 | static unsigned mid2_##ss(const type *const prev, \ | ||
| 167 | const type *const next, \ | ||
| 168 | const type *const prev2, \ | ||
| 169 | const type *const next2, \ | ||
| 170 | const type *const prev3, \ | ||
| 171 | const type *const next3, \ | ||
| 172 | int end, int x, int k, int depth) \ | ||
| 173 | { \ | ||
| 174 | return (prev[av_clip(x + k, 0, end)] + \ | ||
| 175 | next[av_clip(x - k, 0, end)] + 1) >> 1; \ | ||
| 176 | } | ||
| 177 | |||
| 178 | ✗ | MID2(uint8_t, 8) | |
| 179 | ✗ | MID2(uint16_t, 16) | |
| 180 | |||
| 181 | #define MID4(type, ss) \ | ||
| 182 | static unsigned mid4_##ss(const type *const prev, \ | ||
| 183 | const type *const next, \ | ||
| 184 | const type *const prev2, \ | ||
| 185 | const type *const next2, \ | ||
| 186 | const type *const prev3, \ | ||
| 187 | const type *const next3, \ | ||
| 188 | int end, int x, int k, int depth) \ | ||
| 189 | { \ | ||
| 190 | return av_clip_uintp2_c(( \ | ||
| 191 | 9 * (prev[av_clip(x + k, 0, end)] + \ | ||
| 192 | next[av_clip(x - k, 0, end)]) - \ | ||
| 193 | 1 * (prev2[av_clip(x + k*3, 0, end)] + \ | ||
| 194 | next2[av_clip(x - k*3, 0, end)]) + 8) >> 4, \ | ||
| 195 | depth); \ | ||
| 196 | } | ||
| 197 | |||
| 198 | ✗ | MID4(uint8_t, 8) | |
| 199 | ✗ | MID4(uint16_t, 16) | |
| 200 | |||
| 201 | #define MID6(type, ss) \ | ||
| 202 | static unsigned mid6_##ss(const type *const prev, \ | ||
| 203 | const type *const next, \ | ||
| 204 | const type *const prev2, \ | ||
| 205 | const type *const next2, \ | ||
| 206 | const type *const prev3, \ | ||
| 207 | const type *const next3, \ | ||
| 208 | int end, int x, int k, int depth) \ | ||
| 209 | { \ | ||
| 210 | return av_clip_uintp2_c(( \ | ||
| 211 | 20 * (prev[av_clip(x + k, 0, end)] + \ | ||
| 212 | next[av_clip(x - k, 0, end)]) - \ | ||
| 213 | 5 * (prev2[av_clip(x + k*3, 0, end)] + \ | ||
| 214 | next2[av_clip(x - k*3, 0, end)]) + \ | ||
| 215 | 1 * (prev3[av_clip(x + k*5, 0, end)] + \ | ||
| 216 | next3[av_clip(x - k*5, 0, end)]) + 16) >> 5, \ | ||
| 217 | depth); \ | ||
| 218 | } | ||
| 219 | |||
| 220 | ✗ | MID6(uint8_t, 8) | |
| 221 | ✗ | MID6(uint16_t, 16) | |
| 222 | |||
| 223 | #define DIFF(type, ss) \ | ||
| 224 | static unsigned diff_##ss(const type *const prev, \ | ||
| 225 | const type *const next, \ | ||
| 226 | int x, int y) \ | ||
| 227 | { \ | ||
| 228 | return FFABS(prev[x] - next[y]); \ | ||
| 229 | } | ||
| 230 | |||
| 231 | ✗ | DIFF(uint8_t, 8) | |
| 232 | ✗ | DIFF(uint16_t, 16) | |
| 233 | |||
| 234 | #define COST(type, ss) \ | ||
| 235 | static unsigned cost_##ss(const type *const prev, \ | ||
| 236 | const type *const next, \ | ||
| 237 | int end, int x, int k) \ | ||
| 238 | { \ | ||
| 239 | const int m = midl_##ss(prev, next, end, x, k); \ | ||
| 240 | const int p = prev[x]; \ | ||
| 241 | const int n = next[x]; \ | ||
| 242 | \ | ||
| 243 | return FFABS(p - m) + FFABS(n - m); \ | ||
| 244 | } | ||
| 245 | |||
| 246 | ✗ | COST(uint8_t, 8) | |
| 247 | ✗ | COST(uint16_t, 16) | |
| 248 | |||
| 249 | #define INTERPOLATE(type, atype, amax, ss) \ | ||
| 250 | static void interpolate_##ss(ESTDIFContext *s, uint8_t *ddst, \ | ||
| 251 | const uint8_t *const pprev_line, \ | ||
| 252 | const uint8_t *const nnext_line, \ | ||
| 253 | const uint8_t *const pprev2_line, \ | ||
| 254 | const uint8_t *const nnext2_line, \ | ||
| 255 | const uint8_t *const pprev3_line, \ | ||
| 256 | const uint8_t *const nnext3_line, \ | ||
| 257 | int x, int width, int rslope, \ | ||
| 258 | int redge, int depth, \ | ||
| 259 | int *K) \ | ||
| 260 | { \ | ||
| 261 | type *dst = (type *)ddst; \ | ||
| 262 | const type *const prev_line = (const type *const)pprev_line; \ | ||
| 263 | const type *const prev2_line = (const type *const)pprev2_line; \ | ||
| 264 | const type *const prev3_line = (const type *const)pprev3_line; \ | ||
| 265 | const type *const next_line = (const type *const)nnext_line; \ | ||
| 266 | const type *const next2_line = (const type *const)nnext2_line; \ | ||
| 267 | const type *const next3_line = (const type *const)nnext3_line; \ | ||
| 268 | const int interp = s->interp; \ | ||
| 269 | const int ecost = s->ecost; \ | ||
| 270 | const int dcost = s->dcost; \ | ||
| 271 | const int mcost = s->mcost; \ | ||
| 272 | atype sd[S], sD[S], di = 0; \ | ||
| 273 | const int end = width - 1; \ | ||
| 274 | atype dmin = amax; \ | ||
| 275 | int id = 0, iD = 0; \ | ||
| 276 | int k = *K; \ | ||
| 277 | \ | ||
| 278 | for (int i = -rslope; i <= rslope && abs(k) > rslope; i++) { \ | ||
| 279 | atype sum = 0; \ | ||
| 280 | \ | ||
| 281 | for (int j = -redge; j <= redge; j++) { \ | ||
| 282 | const int xx = av_clip(x + i + j, 0, end); \ | ||
| 283 | const int yy = av_clip(x - i + j, 0, end); \ | ||
| 284 | sum += diff_##ss(prev_line, next_line, xx, yy); \ | ||
| 285 | sum += diff_##ss(prev2_line, prev_line, xx, yy); \ | ||
| 286 | sum += diff_##ss(next_line, next2_line, xx, yy); \ | ||
| 287 | } \ | ||
| 288 | \ | ||
| 289 | sD[i + rslope] = ecost * sum; \ | ||
| 290 | sD[i + rslope] += mcost * cost_##ss(prev_line, next_line, end, x, i);\ | ||
| 291 | sD[i + rslope] += dcost * abs(i); \ | ||
| 292 | \ | ||
| 293 | if (dmin > sD[i + rslope]) { \ | ||
| 294 | dmin = sD[i + rslope]; \ | ||
| 295 | di = 1; \ | ||
| 296 | iD = i; \ | ||
| 297 | } \ | ||
| 298 | } \ | ||
| 299 | \ | ||
| 300 | for (int i = -rslope; i <= rslope; i++) { \ | ||
| 301 | atype sum = 0; \ | ||
| 302 | \ | ||
| 303 | for (int j = -redge; j <= redge; j++) { \ | ||
| 304 | const int xx = av_clip(x + k + i + j, 0, end); \ | ||
| 305 | const int yy = av_clip(x - k - i + j, 0, end); \ | ||
| 306 | sum += diff_##ss(prev_line, next_line, xx, yy); \ | ||
| 307 | sum += diff_##ss(prev2_line, prev_line, xx, yy); \ | ||
| 308 | sum += diff_##ss(next_line, next2_line, xx, yy); \ | ||
| 309 | } \ | ||
| 310 | \ | ||
| 311 | sd[i + rslope] = ecost * sum; \ | ||
| 312 | sd[i + rslope] += mcost * cost_##ss(prev_line, next_line, end, x, k+i);\ | ||
| 313 | sd[i + rslope] += dcost * abs(k + i); \ | ||
| 314 | \ | ||
| 315 | if (dmin > sd[i + rslope]) { \ | ||
| 316 | dmin = sd[i + rslope]; \ | ||
| 317 | di = 0; \ | ||
| 318 | id = i; \ | ||
| 319 | } \ | ||
| 320 | } \ | ||
| 321 | \ | ||
| 322 | k = di ? iD : k + id; \ | ||
| 323 | \ | ||
| 324 | dst[x] = s->mid_##ss[interp](prev_line, next_line, \ | ||
| 325 | prev2_line, next2_line, \ | ||
| 326 | prev3_line, next3_line, \ | ||
| 327 | end, x, k, depth); \ | ||
| 328 | \ | ||
| 329 | *K = k; \ | ||
| 330 | } | ||
| 331 | |||
| 332 | ✗ | INTERPOLATE(uint8_t, unsigned, UINT_MAX, 8) | |
| 333 | ✗ | INTERPOLATE(uint16_t, uint64_t, UINT64_MAX, 16) | |
| 334 | |||
| 335 | ✗ | static int deinterlace_slice(AVFilterContext *ctx, void *arg, | |
| 336 | int jobnr, int nb_jobs) | ||
| 337 | { | ||
| 338 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 339 | ✗ | ThreadData *td = arg; | |
| 340 | ✗ | AVFrame *out = td->out; | |
| 341 | ✗ | AVFrame *in = td->in; | |
| 342 | ✗ | const int rslope = s->rslope; | |
| 343 | ✗ | const int redge = s->redge; | |
| 344 | ✗ | const int depth = s->depth; | |
| 345 | ✗ | const int interlaced = !!(in->flags & AV_FRAME_FLAG_INTERLACED); | |
| 346 | ✗ | const int tff = (s->field == (s->parity == -1 ? interlaced ? !!(in->flags & AV_FRAME_FLAG_TOP_FIELD_FIRST) : 1 : | |
| 347 | ✗ | s->parity ^ 1)); | |
| 348 | |||
| 349 | ✗ | for (int plane = 0; plane < s->nb_planes; plane++) { | |
| 350 | ✗ | const uint8_t *src_data = in->data[plane]; | |
| 351 | ✗ | uint8_t *dst_data = out->data[plane]; | |
| 352 | ✗ | const int linesize = s->linesize[plane]; | |
| 353 | ✗ | const int width = s->planewidth[plane]; | |
| 354 | ✗ | const int height = s->planeheight[plane]; | |
| 355 | ✗ | const int src_linesize = in->linesize[plane]; | |
| 356 | ✗ | const int dst_linesize = out->linesize[plane]; | |
| 357 | ✗ | const int start = (height * jobnr) / nb_jobs; | |
| 358 | ✗ | const int end = (height * (jobnr+1)) / nb_jobs; | |
| 359 | const uint8_t *prev_line, *prev2_line, *next_line, *next2_line, *in_line; | ||
| 360 | const uint8_t *prev3_line, *next3_line; | ||
| 361 | uint8_t *out_line; | ||
| 362 | int y_out; | ||
| 363 | |||
| 364 | ✗ | y_out = start + (tff ^ (start & 1)); | |
| 365 | |||
| 366 | ✗ | in_line = src_data + (y_out * src_linesize); | |
| 367 | ✗ | out_line = dst_data + (y_out * dst_linesize); | |
| 368 | |||
| 369 | ✗ | while (y_out < end) { | |
| 370 | ✗ | memcpy(out_line, in_line, linesize); | |
| 371 | ✗ | y_out += 2; | |
| 372 | ✗ | in_line += src_linesize * 2; | |
| 373 | ✗ | out_line += dst_linesize * 2; | |
| 374 | } | ||
| 375 | |||
| 376 | ✗ | y_out = start + ((!tff) ^ (start & 1)); | |
| 377 | ✗ | out_line = dst_data + (y_out * dst_linesize); | |
| 378 | |||
| 379 | ✗ | for (int y = y_out; y < end; y += 2) { | |
| 380 | ✗ | int y_prev3_in = y - 5; | |
| 381 | ✗ | int y_next3_in = y + 5; | |
| 382 | ✗ | int y_prev2_in = y - 3; | |
| 383 | ✗ | int y_next2_in = y + 3; | |
| 384 | ✗ | int y_prev_in = y - 1; | |
| 385 | ✗ | int y_next_in = y + 1; | |
| 386 | int k; | ||
| 387 | |||
| 388 | ✗ | while (y_prev3_in < 0) | |
| 389 | ✗ | y_prev3_in += 2; | |
| 390 | |||
| 391 | ✗ | while (y_next3_in >= height) | |
| 392 | ✗ | y_next3_in -= 2; | |
| 393 | |||
| 394 | ✗ | while (y_prev2_in < 0) | |
| 395 | ✗ | y_prev2_in += 2; | |
| 396 | |||
| 397 | ✗ | while (y_next2_in >= height) | |
| 398 | ✗ | y_next2_in -= 2; | |
| 399 | |||
| 400 | ✗ | while (y_prev_in < 0) | |
| 401 | ✗ | y_prev_in += 2; | |
| 402 | |||
| 403 | ✗ | while (y_next_in >= height) | |
| 404 | ✗ | y_next_in -= 2; | |
| 405 | |||
| 406 | ✗ | prev3_line = src_data + (y_prev3_in * src_linesize); | |
| 407 | ✗ | next3_line = src_data + (y_next3_in * src_linesize); | |
| 408 | |||
| 409 | ✗ | prev2_line = src_data + (y_prev2_in * src_linesize); | |
| 410 | ✗ | next2_line = src_data + (y_next2_in * src_linesize); | |
| 411 | |||
| 412 | ✗ | prev_line = src_data + (y_prev_in * src_linesize); | |
| 413 | ✗ | next_line = src_data + (y_next_in * src_linesize); | |
| 414 | |||
| 415 | ✗ | k = 0; | |
| 416 | |||
| 417 | ✗ | for (int x = 0; x < width; x++) { | |
| 418 | ✗ | s->interpolate(s, out_line, | |
| 419 | prev_line, next_line, | ||
| 420 | prev2_line, next2_line, | ||
| 421 | prev3_line, next3_line, | ||
| 422 | x, width, rslope, redge, depth, &k); | ||
| 423 | } | ||
| 424 | |||
| 425 | ✗ | out_line += 2 * dst_linesize; | |
| 426 | } | ||
| 427 | } | ||
| 428 | |||
| 429 | ✗ | return 0; | |
| 430 | } | ||
| 431 | |||
| 432 | ✗ | static int filter(AVFilterContext *ctx, AVFrame *in, int64_t pts, int64_t duration) | |
| 433 | { | ||
| 434 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 435 | ✗ | AVFilterLink *outlink = ctx->outputs[0]; | |
| 436 | AVFrame *out; | ||
| 437 | ThreadData td; | ||
| 438 | |||
| 439 | ✗ | out = ff_get_video_buffer(outlink, outlink->w, outlink->h); | |
| 440 | ✗ | if (!out) | |
| 441 | ✗ | return AVERROR(ENOMEM); | |
| 442 | ✗ | av_frame_copy_props(out, in); | |
| 443 | ✗ | out->flags &= ~AV_FRAME_FLAG_INTERLACED; | |
| 444 | ✗ | out->pts = pts; | |
| 445 | ✗ | out->duration = duration; | |
| 446 | |||
| 447 | ✗ | td.out = out; td.in = in; | |
| 448 | ✗ | ff_filter_execute(ctx, deinterlace_slice, &td, NULL, | |
| 449 | ✗ | FFMIN(s->planeheight[1] / 2, s->nb_threads)); | |
| 450 | |||
| 451 | ✗ | if (s->mode) | |
| 452 | ✗ | s->field = !s->field; | |
| 453 | |||
| 454 | ✗ | return ff_filter_frame(outlink, out); | |
| 455 | } | ||
| 456 | |||
| 457 | ✗ | static int config_input(AVFilterLink *inlink) | |
| 458 | { | ||
| 459 | ✗ | AVFilterContext *ctx = inlink->dst; | |
| 460 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 461 | ✗ | const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format); | |
| 462 | int ret; | ||
| 463 | |||
| 464 | ✗ | if ((ret = av_image_fill_linesizes(s->linesize, inlink->format, inlink->w)) < 0) | |
| 465 | ✗ | return ret; | |
| 466 | |||
| 467 | ✗ | s->planeheight[1] = s->planeheight[2] = AV_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h); | |
| 468 | ✗ | s->planeheight[0] = s->planeheight[3] = inlink->h; | |
| 469 | ✗ | s->planewidth[1] = s->planewidth[2] = AV_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w); | |
| 470 | ✗ | s->planewidth[0] = s->planewidth[3] = inlink->w; | |
| 471 | |||
| 472 | ✗ | if (inlink->h < 3) { | |
| 473 | ✗ | av_log(ctx, AV_LOG_ERROR, "Video of less than 3 lines is not supported\n"); | |
| 474 | ✗ | return AVERROR(EINVAL); | |
| 475 | } | ||
| 476 | |||
| 477 | ✗ | s->nb_planes = av_pix_fmt_count_planes(inlink->format); | |
| 478 | ✗ | s->nb_threads = ff_filter_get_nb_threads(ctx); | |
| 479 | ✗ | s->depth = desc->comp[0].depth; | |
| 480 | ✗ | s->interpolate = s->depth <= 8 ? interpolate_8 : interpolate_16; | |
| 481 | ✗ | s->mid_8[0] = mid2_8; | |
| 482 | ✗ | s->mid_8[1] = mid4_8; | |
| 483 | ✗ | s->mid_8[2] = mid6_8; | |
| 484 | ✗ | s->mid_16[0] = mid2_16; | |
| 485 | ✗ | s->mid_16[1] = mid4_16; | |
| 486 | ✗ | s->mid_16[2] = mid6_16; | |
| 487 | ✗ | s->max = (1 << (s->depth)) - 1; | |
| 488 | |||
| 489 | ✗ | return 0; | |
| 490 | } | ||
| 491 | ✗ | static int filter_frame(AVFilterLink *inlink, AVFrame *in) | |
| 492 | { | ||
| 493 | ✗ | AVFilterContext *ctx = inlink->dst; | |
| 494 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 495 | int ret; | ||
| 496 | |||
| 497 | ✗ | if (!s->prev) { | |
| 498 | ✗ | s->prev = in; | |
| 499 | ✗ | return 0; | |
| 500 | } | ||
| 501 | |||
| 502 | ✗ | if ((s->deint && !(s->prev->flags & AV_FRAME_FLAG_INTERLACED)) || ctx->is_disabled) { | |
| 503 | ✗ | s->prev->pts *= 2; | |
| 504 | ✗ | s->prev->duration *= 2; | |
| 505 | ✗ | ret = ff_filter_frame(ctx->outputs[0], s->prev); | |
| 506 | ✗ | s->prev = in; | |
| 507 | ✗ | return ret; | |
| 508 | } | ||
| 509 | |||
| 510 | ✗ | ret = filter(ctx, s->prev, s->prev->pts * 2, | |
| 511 | ✗ | s->prev->duration * (s->mode ? 1 : 2)); | |
| 512 | ✗ | if (ret < 0 || s->mode == 0) { | |
| 513 | ✗ | av_frame_free(&s->prev); | |
| 514 | ✗ | s->prev = in; | |
| 515 | ✗ | return ret; | |
| 516 | } | ||
| 517 | |||
| 518 | ✗ | ret = filter(ctx, s->prev, s->prev->pts + in->pts, in->duration); | |
| 519 | ✗ | av_frame_free(&s->prev); | |
| 520 | ✗ | s->prev = in; | |
| 521 | ✗ | return ret; | |
| 522 | } | ||
| 523 | |||
| 524 | ✗ | static int request_frame(AVFilterLink *link) | |
| 525 | { | ||
| 526 | ✗ | FilterLink *l = ff_filter_link(link); | |
| 527 | ✗ | AVFilterContext *ctx = link->src; | |
| 528 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 529 | int ret; | ||
| 530 | |||
| 531 | ✗ | if (s->eof) | |
| 532 | ✗ | return AVERROR_EOF; | |
| 533 | |||
| 534 | ✗ | ret = ff_request_frame(ctx->inputs[0]); | |
| 535 | |||
| 536 | ✗ | if (ret == AVERROR_EOF && s->prev) { | |
| 537 | ✗ | AVFrame *next = av_frame_clone(s->prev); | |
| 538 | |||
| 539 | ✗ | if (!next) | |
| 540 | ✗ | return AVERROR(ENOMEM); | |
| 541 | |||
| 542 | ✗ | next->pts = s->prev->pts + av_rescale_q(1, av_inv_q(l->frame_rate), | |
| 543 | ✗ | ctx->outputs[0]->time_base); | |
| 544 | ✗ | s->eof = 1; | |
| 545 | ✗ | ret = filter_frame(ctx->inputs[0], next); | |
| 546 | ✗ | } else if (ret < 0) { | |
| 547 | ✗ | return ret; | |
| 548 | } | ||
| 549 | |||
| 550 | ✗ | return ret; | |
| 551 | } | ||
| 552 | |||
| 553 | ✗ | static av_cold void uninit(AVFilterContext *ctx) | |
| 554 | { | ||
| 555 | ✗ | ESTDIFContext *s = ctx->priv; | |
| 556 | |||
| 557 | ✗ | av_frame_free(&s->prev); | |
| 558 | ✗ | } | |
| 559 | |||
| 560 | static const AVFilterPad estdif_inputs[] = { | ||
| 561 | { | ||
| 562 | .name = "default", | ||
| 563 | .type = AVMEDIA_TYPE_VIDEO, | ||
| 564 | .filter_frame = filter_frame, | ||
| 565 | .config_props = config_input, | ||
| 566 | }, | ||
| 567 | }; | ||
| 568 | |||
| 569 | static const AVFilterPad estdif_outputs[] = { | ||
| 570 | { | ||
| 571 | .name = "default", | ||
| 572 | .type = AVMEDIA_TYPE_VIDEO, | ||
| 573 | .config_props = config_output, | ||
| 574 | .request_frame = request_frame, | ||
| 575 | }, | ||
| 576 | }; | ||
| 577 | |||
| 578 | const FFFilter ff_vf_estdif = { | ||
| 579 | .p.name = "estdif", | ||
| 580 | .p.description = NULL_IF_CONFIG_SMALL("Apply Edge Slope Tracing deinterlace."), | ||
| 581 | .p.priv_class = &estdif_class, | ||
| 582 | .p.flags = AVFILTER_FLAG_SUPPORT_TIMELINE_INTERNAL | AVFILTER_FLAG_SLICE_THREADS, | ||
| 583 | .priv_size = sizeof(ESTDIFContext), | ||
| 584 | .uninit = uninit, | ||
| 585 | FILTER_INPUTS(estdif_inputs), | ||
| 586 | FILTER_OUTPUTS(estdif_outputs), | ||
| 587 | FILTER_PIXFMTS_ARRAY(pix_fmts), | ||
| 588 | .process_command = ff_filter_process_command, | ||
| 589 | }; | ||
| 590 |