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 |
|
|
#if FF_API_INTERLACED_FRAME |
444 |
|
|
FF_DISABLE_DEPRECATION_WARNINGS |
445 |
|
✗ |
out->interlaced_frame = 0; |
446 |
|
|
FF_ENABLE_DEPRECATION_WARNINGS |
447 |
|
|
#endif |
448 |
|
✗ |
out->flags &= ~AV_FRAME_FLAG_INTERLACED; |
449 |
|
✗ |
out->pts = pts; |
450 |
|
✗ |
out->duration = duration; |
451 |
|
|
|
452 |
|
✗ |
td.out = out; td.in = in; |
453 |
|
✗ |
ff_filter_execute(ctx, deinterlace_slice, &td, NULL, |
454 |
|
✗ |
FFMIN(s->planeheight[1] / 2, s->nb_threads)); |
455 |
|
|
|
456 |
|
✗ |
if (s->mode) |
457 |
|
✗ |
s->field = !s->field; |
458 |
|
|
|
459 |
|
✗ |
return ff_filter_frame(outlink, out); |
460 |
|
|
} |
461 |
|
|
|
462 |
|
✗ |
static int config_input(AVFilterLink *inlink) |
463 |
|
|
{ |
464 |
|
✗ |
AVFilterContext *ctx = inlink->dst; |
465 |
|
✗ |
ESTDIFContext *s = ctx->priv; |
466 |
|
✗ |
const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format); |
467 |
|
|
int ret; |
468 |
|
|
|
469 |
|
✗ |
if ((ret = av_image_fill_linesizes(s->linesize, inlink->format, inlink->w)) < 0) |
470 |
|
✗ |
return ret; |
471 |
|
|
|
472 |
|
✗ |
s->planeheight[1] = s->planeheight[2] = AV_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h); |
473 |
|
✗ |
s->planeheight[0] = s->planeheight[3] = inlink->h; |
474 |
|
✗ |
s->planewidth[1] = s->planewidth[2] = AV_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w); |
475 |
|
✗ |
s->planewidth[0] = s->planewidth[3] = inlink->w; |
476 |
|
|
|
477 |
|
✗ |
if (inlink->h < 3) { |
478 |
|
✗ |
av_log(ctx, AV_LOG_ERROR, "Video of less than 3 lines is not supported\n"); |
479 |
|
✗ |
return AVERROR(EINVAL); |
480 |
|
|
} |
481 |
|
|
|
482 |
|
✗ |
s->nb_planes = av_pix_fmt_count_planes(inlink->format); |
483 |
|
✗ |
s->nb_threads = ff_filter_get_nb_threads(ctx); |
484 |
|
✗ |
s->depth = desc->comp[0].depth; |
485 |
|
✗ |
s->interpolate = s->depth <= 8 ? interpolate_8 : interpolate_16; |
486 |
|
✗ |
s->mid_8[0] = mid2_8; |
487 |
|
✗ |
s->mid_8[1] = mid4_8; |
488 |
|
✗ |
s->mid_8[2] = mid6_8; |
489 |
|
✗ |
s->mid_16[0] = mid2_16; |
490 |
|
✗ |
s->mid_16[1] = mid4_16; |
491 |
|
✗ |
s->mid_16[2] = mid6_16; |
492 |
|
✗ |
s->max = (1 << (s->depth)) - 1; |
493 |
|
|
|
494 |
|
✗ |
return 0; |
495 |
|
|
} |
496 |
|
✗ |
static int filter_frame(AVFilterLink *inlink, AVFrame *in) |
497 |
|
|
{ |
498 |
|
✗ |
AVFilterContext *ctx = inlink->dst; |
499 |
|
✗ |
ESTDIFContext *s = ctx->priv; |
500 |
|
|
int ret; |
501 |
|
|
|
502 |
|
✗ |
if (!s->prev) { |
503 |
|
✗ |
s->prev = in; |
504 |
|
✗ |
return 0; |
505 |
|
|
} |
506 |
|
|
|
507 |
|
✗ |
if ((s->deint && !(s->prev->flags & AV_FRAME_FLAG_INTERLACED)) || ctx->is_disabled) { |
508 |
|
✗ |
s->prev->pts *= 2; |
509 |
|
✗ |
s->prev->duration *= 2; |
510 |
|
✗ |
ret = ff_filter_frame(ctx->outputs[0], s->prev); |
511 |
|
✗ |
s->prev = in; |
512 |
|
✗ |
return ret; |
513 |
|
|
} |
514 |
|
|
|
515 |
|
✗ |
ret = filter(ctx, s->prev, s->prev->pts * 2, |
516 |
|
✗ |
s->prev->duration * (s->mode ? 1 : 2)); |
517 |
|
✗ |
if (ret < 0 || s->mode == 0) { |
518 |
|
✗ |
av_frame_free(&s->prev); |
519 |
|
✗ |
s->prev = in; |
520 |
|
✗ |
return ret; |
521 |
|
|
} |
522 |
|
|
|
523 |
|
✗ |
ret = filter(ctx, s->prev, s->prev->pts + in->pts, in->duration); |
524 |
|
✗ |
av_frame_free(&s->prev); |
525 |
|
✗ |
s->prev = in; |
526 |
|
✗ |
return ret; |
527 |
|
|
} |
528 |
|
|
|
529 |
|
✗ |
static int request_frame(AVFilterLink *link) |
530 |
|
|
{ |
531 |
|
✗ |
FilterLink *l = ff_filter_link(link); |
532 |
|
✗ |
AVFilterContext *ctx = link->src; |
533 |
|
✗ |
ESTDIFContext *s = ctx->priv; |
534 |
|
|
int ret; |
535 |
|
|
|
536 |
|
✗ |
if (s->eof) |
537 |
|
✗ |
return AVERROR_EOF; |
538 |
|
|
|
539 |
|
✗ |
ret = ff_request_frame(ctx->inputs[0]); |
540 |
|
|
|
541 |
|
✗ |
if (ret == AVERROR_EOF && s->prev) { |
542 |
|
✗ |
AVFrame *next = av_frame_clone(s->prev); |
543 |
|
|
|
544 |
|
✗ |
if (!next) |
545 |
|
✗ |
return AVERROR(ENOMEM); |
546 |
|
|
|
547 |
|
✗ |
next->pts = s->prev->pts + av_rescale_q(1, av_inv_q(l->frame_rate), |
548 |
|
✗ |
ctx->outputs[0]->time_base); |
549 |
|
✗ |
s->eof = 1; |
550 |
|
✗ |
ret = filter_frame(ctx->inputs[0], next); |
551 |
|
✗ |
} else if (ret < 0) { |
552 |
|
✗ |
return ret; |
553 |
|
|
} |
554 |
|
|
|
555 |
|
✗ |
return ret; |
556 |
|
|
} |
557 |
|
|
|
558 |
|
✗ |
static av_cold void uninit(AVFilterContext *ctx) |
559 |
|
|
{ |
560 |
|
✗ |
ESTDIFContext *s = ctx->priv; |
561 |
|
|
|
562 |
|
✗ |
av_frame_free(&s->prev); |
563 |
|
✗ |
} |
564 |
|
|
|
565 |
|
|
static const AVFilterPad estdif_inputs[] = { |
566 |
|
|
{ |
567 |
|
|
.name = "default", |
568 |
|
|
.type = AVMEDIA_TYPE_VIDEO, |
569 |
|
|
.filter_frame = filter_frame, |
570 |
|
|
.config_props = config_input, |
571 |
|
|
}, |
572 |
|
|
}; |
573 |
|
|
|
574 |
|
|
static const AVFilterPad estdif_outputs[] = { |
575 |
|
|
{ |
576 |
|
|
.name = "default", |
577 |
|
|
.type = AVMEDIA_TYPE_VIDEO, |
578 |
|
|
.config_props = config_output, |
579 |
|
|
.request_frame = request_frame, |
580 |
|
|
}, |
581 |
|
|
}; |
582 |
|
|
|
583 |
|
|
const AVFilter ff_vf_estdif = { |
584 |
|
|
.name = "estdif", |
585 |
|
|
.description = NULL_IF_CONFIG_SMALL("Apply Edge Slope Tracing deinterlace."), |
586 |
|
|
.priv_size = sizeof(ESTDIFContext), |
587 |
|
|
.priv_class = &estdif_class, |
588 |
|
|
.uninit = uninit, |
589 |
|
|
FILTER_INPUTS(estdif_inputs), |
590 |
|
|
FILTER_OUTPUTS(estdif_outputs), |
591 |
|
|
FILTER_PIXFMTS_ARRAY(pix_fmts), |
592 |
|
|
.flags = AVFILTER_FLAG_SUPPORT_TIMELINE_INTERNAL | AVFILTER_FLAG_SLICE_THREADS, |
593 |
|
|
.process_command = ff_filter_process_command, |
594 |
|
|
}; |
595 |
|
|
|