Commit | Line | Data |
---|---|---|
2ba45a60 DM |
1 | /* |
2 | * Copyright (C) 2012 Michael Niedermayer <michaelni@gmx.at> | |
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 <float.h> /* FLT_MAX */ | |
22 | ||
23 | #include "libavutil/cpu.h" | |
24 | #include "libavutil/common.h" | |
25 | #include "libavutil/opt.h" | |
26 | #include "internal.h" | |
27 | #include "vf_idet.h" | |
28 | ||
29 | #define OFFSET(x) offsetof(IDETContext, x) | |
30 | #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM | |
31 | ||
32 | static const AVOption idet_options[] = { | |
33 | { "intl_thres", "set interlacing threshold", OFFSET(interlace_threshold), AV_OPT_TYPE_FLOAT, {.dbl = 1.04}, -1, FLT_MAX, FLAGS }, | |
34 | { "prog_thres", "set progressive threshold", OFFSET(progressive_threshold), AV_OPT_TYPE_FLOAT, {.dbl = 1.5}, -1, FLT_MAX, FLAGS }, | |
f6fa7814 DM |
35 | { "rep_thres", "set repeat threshold", OFFSET(repeat_threshold), AV_OPT_TYPE_FLOAT, {.dbl = 3.0}, -1, FLT_MAX, FLAGS }, |
36 | { "half_life", "half life of cumulative statistics", OFFSET(half_life), AV_OPT_TYPE_FLOAT, {.dbl = 0.0}, -1, INT_MAX, FLAGS }, | |
2ba45a60 DM |
37 | { NULL } |
38 | }; | |
39 | ||
40 | AVFILTER_DEFINE_CLASS(idet); | |
41 | ||
42 | static const char *type2str(Type type) | |
43 | { | |
44 | switch(type) { | |
f6fa7814 DM |
45 | case TFF : return "tff"; |
46 | case BFF : return "bff"; | |
47 | case PROGRESSIVE : return "progressive"; | |
48 | case UNDETERMINED : return "undetermined"; | |
49 | } | |
50 | return NULL; | |
51 | } | |
52 | ||
53 | #define PRECISION 1048576 | |
54 | ||
55 | static uint64_t uintpow(uint64_t b,unsigned int e) | |
56 | { | |
57 | uint64_t r=1; | |
58 | while(e--) r*=b; | |
59 | return r; | |
60 | } | |
61 | ||
62 | static int av_dict_set_fxp(AVDictionary **pm, const char *key, uint64_t value, unsigned int digits, | |
63 | int flags) | |
64 | { | |
65 | char valuestr[44]; | |
66 | uint64_t print_precision = uintpow(10, digits); | |
67 | ||
68 | value = av_rescale(value, print_precision, PRECISION); | |
69 | ||
70 | snprintf(valuestr, sizeof(valuestr), "%"PRId64".%0*"PRId64, | |
71 | value / print_precision, digits, value % print_precision); | |
72 | ||
73 | return av_dict_set(pm, key, valuestr, flags); | |
74 | } | |
75 | ||
76 | static const char *rep2str(RepeatedField repeated_field) | |
77 | { | |
78 | switch(repeated_field) { | |
79 | case REPEAT_NONE : return "neither"; | |
80 | case REPEAT_TOP : return "top"; | |
81 | case REPEAT_BOTTOM : return "bottom"; | |
2ba45a60 DM |
82 | } |
83 | return NULL; | |
84 | } | |
85 | ||
86 | int ff_idet_filter_line_c(const uint8_t *a, const uint8_t *b, const uint8_t *c, int w) | |
87 | { | |
88 | int x; | |
89 | int ret=0; | |
90 | ||
91 | for(x=0; x<w; x++){ | |
92 | int v = (*a++ + *c++) - 2 * *b++; | |
93 | ret += FFABS(v); | |
94 | } | |
95 | ||
96 | return ret; | |
97 | } | |
98 | ||
99 | int ff_idet_filter_line_c_16bit(const uint16_t *a, const uint16_t *b, const uint16_t *c, int w) | |
100 | { | |
101 | int x; | |
102 | int ret=0; | |
103 | ||
104 | for(x=0; x<w; x++){ | |
105 | int v = (*a++ + *c++) - 2 * *b++; | |
106 | ret += FFABS(v); | |
107 | } | |
108 | ||
109 | return ret; | |
110 | } | |
111 | ||
112 | static void filter(AVFilterContext *ctx) | |
113 | { | |
114 | IDETContext *idet = ctx->priv; | |
115 | int y, i; | |
116 | int64_t alpha[2]={0}; | |
117 | int64_t delta=0; | |
f6fa7814 | 118 | int64_t gamma[2]={0}; |
2ba45a60 | 119 | Type type, best_type; |
f6fa7814 | 120 | RepeatedField repeat; |
2ba45a60 | 121 | int match = 0; |
f6fa7814 | 122 | AVDictionary **metadata = avpriv_frame_get_metadatap(idet->cur); |
2ba45a60 DM |
123 | |
124 | for (i = 0; i < idet->csp->nb_components; i++) { | |
125 | int w = idet->cur->width; | |
126 | int h = idet->cur->height; | |
127 | int refs = idet->cur->linesize[i]; | |
128 | ||
129 | if (i && i<3) { | |
130 | w = FF_CEIL_RSHIFT(w, idet->csp->log2_chroma_w); | |
131 | h = FF_CEIL_RSHIFT(h, idet->csp->log2_chroma_h); | |
132 | } | |
133 | ||
134 | for (y = 2; y < h - 2; y++) { | |
135 | uint8_t *prev = &idet->prev->data[i][y*refs]; | |
136 | uint8_t *cur = &idet->cur ->data[i][y*refs]; | |
137 | uint8_t *next = &idet->next->data[i][y*refs]; | |
138 | alpha[ y &1] += idet->filter_line(cur-refs, prev, cur+refs, w); | |
139 | alpha[(y^1)&1] += idet->filter_line(cur-refs, next, cur+refs, w); | |
140 | delta += idet->filter_line(cur-refs, cur, cur+refs, w); | |
f6fa7814 | 141 | gamma[(y^1)&1] += idet->filter_line(cur , prev, cur , w); |
2ba45a60 DM |
142 | } |
143 | } | |
144 | ||
145 | if (alpha[0] > idet->interlace_threshold * alpha[1]){ | |
146 | type = TFF; | |
147 | }else if(alpha[1] > idet->interlace_threshold * alpha[0]){ | |
148 | type = BFF; | |
149 | }else if(alpha[1] > idet->progressive_threshold * delta){ | |
f6fa7814 | 150 | type = PROGRESSIVE; |
2ba45a60 DM |
151 | }else{ |
152 | type = UNDETERMINED; | |
153 | } | |
154 | ||
f6fa7814 DM |
155 | if ( gamma[0] > idet->repeat_threshold * gamma[1] ){ |
156 | repeat = REPEAT_TOP; | |
157 | } else if ( gamma[1] > idet->repeat_threshold * gamma[0] ){ | |
158 | repeat = REPEAT_BOTTOM; | |
159 | } else { | |
160 | repeat = REPEAT_NONE; | |
161 | } | |
162 | ||
2ba45a60 DM |
163 | memmove(idet->history+1, idet->history, HIST_SIZE-1); |
164 | idet->history[0] = type; | |
165 | best_type = UNDETERMINED; | |
166 | for(i=0; i<HIST_SIZE; i++){ | |
167 | if(idet->history[i] != UNDETERMINED){ | |
168 | if(best_type == UNDETERMINED) | |
169 | best_type = idet->history[i]; | |
170 | ||
171 | if(idet->history[i] == best_type) { | |
172 | match++; | |
173 | }else{ | |
174 | match=0; | |
175 | break; | |
176 | } | |
177 | } | |
178 | } | |
179 | if(idet->last_type == UNDETERMINED){ | |
180 | if(match ) idet->last_type = best_type; | |
181 | }else{ | |
182 | if(match>2) idet->last_type = best_type; | |
183 | } | |
184 | ||
185 | if (idet->last_type == TFF){ | |
186 | idet->cur->top_field_first = 1; | |
187 | idet->cur->interlaced_frame = 1; | |
188 | }else if(idet->last_type == BFF){ | |
189 | idet->cur->top_field_first = 0; | |
190 | idet->cur->interlaced_frame = 1; | |
f6fa7814 | 191 | }else if(idet->last_type == PROGRESSIVE){ |
2ba45a60 DM |
192 | idet->cur->interlaced_frame = 0; |
193 | } | |
194 | ||
f6fa7814 DM |
195 | for(i=0; i<3; i++) |
196 | idet->repeats[i] = av_rescale(idet->repeats [i], idet->decay_coefficient, PRECISION); | |
197 | ||
198 | for(i=0; i<4; i++){ | |
199 | idet->prestat [i] = av_rescale(idet->prestat [i], idet->decay_coefficient, PRECISION); | |
200 | idet->poststat[i] = av_rescale(idet->poststat[i], idet->decay_coefficient, PRECISION); | |
201 | } | |
202 | ||
203 | idet->total_repeats [ repeat] ++; | |
204 | idet->repeats [ repeat] += PRECISION; | |
205 | ||
206 | idet->total_prestat [ type] ++; | |
207 | idet->prestat [ type] += PRECISION; | |
208 | ||
209 | idet->total_poststat[idet->last_type] ++; | |
210 | idet->poststat [idet->last_type] += PRECISION; | |
211 | ||
212 | av_log(ctx, AV_LOG_DEBUG, "Repeated Field:%12s, Single frame:%12s, Multi frame:%12s\n", | |
213 | rep2str(repeat), type2str(type), type2str(idet->last_type)); | |
214 | ||
215 | av_dict_set (metadata, "lavfi.idet.repeated.current_frame", rep2str(repeat), 0); | |
216 | av_dict_set_fxp(metadata, "lavfi.idet.repeated.neither", idet->repeats[REPEAT_NONE], 2, 0); | |
217 | av_dict_set_fxp(metadata, "lavfi.idet.repeated.top", idet->repeats[REPEAT_TOP], 2, 0); | |
218 | av_dict_set_fxp(metadata, "lavfi.idet.repeated.bottom", idet->repeats[REPEAT_BOTTOM], 2, 0); | |
219 | ||
220 | av_dict_set (metadata, "lavfi.idet.single.current_frame", type2str(type), 0); | |
221 | av_dict_set_fxp(metadata, "lavfi.idet.single.tff", idet->prestat[TFF], 2 , 0); | |
222 | av_dict_set_fxp(metadata, "lavfi.idet.single.bff", idet->prestat[BFF], 2, 0); | |
223 | av_dict_set_fxp(metadata, "lavfi.idet.single.progressive", idet->prestat[PROGRESSIVE], 2, 0); | |
224 | av_dict_set_fxp(metadata, "lavfi.idet.single.undetermined", idet->prestat[UNDETERMINED], 2, 0); | |
225 | ||
226 | av_dict_set (metadata, "lavfi.idet.multiple.current_frame", type2str(idet->last_type), 0); | |
227 | av_dict_set_fxp(metadata, "lavfi.idet.multiple.tff", idet->poststat[TFF], 2, 0); | |
228 | av_dict_set_fxp(metadata, "lavfi.idet.multiple.bff", idet->poststat[BFF], 2, 0); | |
229 | av_dict_set_fxp(metadata, "lavfi.idet.multiple.progressive", idet->poststat[PROGRESSIVE], 2, 0); | |
230 | av_dict_set_fxp(metadata, "lavfi.idet.multiple.undetermined", idet->poststat[UNDETERMINED], 2, 0); | |
2ba45a60 DM |
231 | } |
232 | ||
233 | static int filter_frame(AVFilterLink *link, AVFrame *picref) | |
234 | { | |
235 | AVFilterContext *ctx = link->dst; | |
236 | IDETContext *idet = ctx->priv; | |
237 | ||
238 | if (idet->prev) | |
239 | av_frame_free(&idet->prev); | |
240 | idet->prev = idet->cur; | |
241 | idet->cur = idet->next; | |
242 | idet->next = picref; | |
243 | ||
f6fa7814 DM |
244 | if (!idet->cur && |
245 | !(idet->cur = av_frame_clone(idet->next))) | |
246 | return AVERROR(ENOMEM); | |
2ba45a60 DM |
247 | |
248 | if (!idet->prev) | |
f6fa7814 | 249 | return 0; |
2ba45a60 DM |
250 | |
251 | if (!idet->csp) | |
252 | idet->csp = av_pix_fmt_desc_get(link->format); | |
253 | if (idet->csp->comp[0].depth_minus1 / 8 == 1){ | |
254 | idet->filter_line = (ff_idet_filter_func)ff_idet_filter_line_c_16bit; | |
255 | if (ARCH_X86) | |
256 | ff_idet_init_x86(idet, 1); | |
257 | } | |
258 | ||
259 | filter(ctx); | |
260 | ||
261 | return ff_filter_frame(ctx->outputs[0], av_frame_clone(idet->cur)); | |
262 | } | |
263 | ||
f6fa7814 DM |
264 | static int request_frame(AVFilterLink *link) |
265 | { | |
266 | AVFilterContext *ctx = link->src; | |
267 | IDETContext *idet = ctx->priv; | |
268 | ||
269 | do { | |
270 | int ret; | |
271 | ||
272 | if (idet->eof) | |
273 | return AVERROR_EOF; | |
274 | ||
275 | ret = ff_request_frame(link->src->inputs[0]); | |
276 | ||
277 | if (ret == AVERROR_EOF && idet->cur) { | |
278 | AVFrame *next = av_frame_clone(idet->next); | |
279 | ||
280 | if (!next) | |
281 | return AVERROR(ENOMEM); | |
282 | ||
283 | filter_frame(link->src->inputs[0], next); | |
284 | idet->eof = 1; | |
285 | } else if (ret < 0) { | |
286 | return ret; | |
287 | } | |
288 | } while (!idet->prev); | |
289 | ||
290 | return 0; | |
291 | } | |
292 | ||
2ba45a60 DM |
293 | static av_cold void uninit(AVFilterContext *ctx) |
294 | { | |
295 | IDETContext *idet = ctx->priv; | |
296 | ||
f6fa7814 DM |
297 | av_log(ctx, AV_LOG_INFO, "Repeated Fields: Neither:%6"PRId64" Top:%6"PRId64" Bottom:%6"PRId64"\n", |
298 | idet->total_repeats[REPEAT_NONE], | |
299 | idet->total_repeats[REPEAT_TOP], | |
300 | idet->total_repeats[REPEAT_BOTTOM] | |
301 | ); | |
302 | av_log(ctx, AV_LOG_INFO, "Single frame detection: TFF:%6"PRId64" BFF:%6"PRId64" Progressive:%6"PRId64" Undetermined:%6"PRId64"\n", | |
303 | idet->total_prestat[TFF], | |
304 | idet->total_prestat[BFF], | |
305 | idet->total_prestat[PROGRESSIVE], | |
306 | idet->total_prestat[UNDETERMINED] | |
307 | ); | |
308 | av_log(ctx, AV_LOG_INFO, "Multi frame detection: TFF:%6"PRId64" BFF:%6"PRId64" Progressive:%6"PRId64" Undetermined:%6"PRId64"\n", | |
309 | idet->total_poststat[TFF], | |
310 | idet->total_poststat[BFF], | |
311 | idet->total_poststat[PROGRESSIVE], | |
312 | idet->total_poststat[UNDETERMINED] | |
313 | ); | |
2ba45a60 DM |
314 | |
315 | av_frame_free(&idet->prev); | |
316 | av_frame_free(&idet->cur ); | |
317 | av_frame_free(&idet->next); | |
318 | } | |
319 | ||
320 | static int query_formats(AVFilterContext *ctx) | |
321 | { | |
322 | static const enum AVPixelFormat pix_fmts[] = { | |
323 | AV_PIX_FMT_YUV420P, | |
324 | AV_PIX_FMT_YUV422P, | |
325 | AV_PIX_FMT_YUV444P, | |
326 | AV_PIX_FMT_YUV410P, | |
327 | AV_PIX_FMT_YUV411P, | |
328 | AV_PIX_FMT_GRAY8, | |
329 | AV_PIX_FMT_YUVJ420P, | |
330 | AV_PIX_FMT_YUVJ422P, | |
331 | AV_PIX_FMT_YUVJ444P, | |
332 | AV_PIX_FMT_GRAY16, | |
333 | AV_PIX_FMT_YUV440P, | |
334 | AV_PIX_FMT_YUVJ440P, | |
335 | AV_PIX_FMT_YUV420P10, | |
336 | AV_PIX_FMT_YUV422P10, | |
337 | AV_PIX_FMT_YUV444P10, | |
338 | AV_PIX_FMT_YUV420P16, | |
339 | AV_PIX_FMT_YUV422P16, | |
340 | AV_PIX_FMT_YUV444P16, | |
341 | AV_PIX_FMT_YUVA420P, | |
342 | AV_PIX_FMT_NONE | |
343 | }; | |
344 | ||
345 | ff_set_common_formats(ctx, ff_make_format_list(pix_fmts)); | |
346 | ||
347 | return 0; | |
348 | } | |
349 | ||
350 | static int config_output(AVFilterLink *outlink) | |
351 | { | |
352 | outlink->flags |= FF_LINK_FLAG_REQUEST_LOOP; | |
353 | return 0; | |
354 | } | |
355 | ||
356 | static av_cold int init(AVFilterContext *ctx) | |
357 | { | |
358 | IDETContext *idet = ctx->priv; | |
359 | ||
f6fa7814 | 360 | idet->eof = 0; |
2ba45a60 DM |
361 | idet->last_type = UNDETERMINED; |
362 | memset(idet->history, UNDETERMINED, HIST_SIZE); | |
363 | ||
f6fa7814 DM |
364 | if( idet->half_life > 0 ) |
365 | idet->decay_coefficient = (uint64_t) round( PRECISION * exp2(-1.0 / idet->half_life) ); | |
366 | else | |
367 | idet->decay_coefficient = PRECISION; | |
368 | ||
2ba45a60 DM |
369 | idet->filter_line = ff_idet_filter_line_c; |
370 | ||
371 | if (ARCH_X86) | |
372 | ff_idet_init_x86(idet, 0); | |
373 | ||
374 | return 0; | |
375 | } | |
376 | ||
2ba45a60 DM |
377 | static const AVFilterPad idet_inputs[] = { |
378 | { | |
379 | .name = "default", | |
380 | .type = AVMEDIA_TYPE_VIDEO, | |
381 | .filter_frame = filter_frame, | |
382 | }, | |
383 | { NULL } | |
384 | }; | |
385 | ||
386 | static const AVFilterPad idet_outputs[] = { | |
387 | { | |
388 | .name = "default", | |
389 | .type = AVMEDIA_TYPE_VIDEO, | |
390 | .config_props = config_output, | |
f6fa7814 | 391 | .request_frame = request_frame |
2ba45a60 DM |
392 | }, |
393 | { NULL } | |
394 | }; | |
395 | ||
396 | AVFilter ff_vf_idet = { | |
397 | .name = "idet", | |
398 | .description = NULL_IF_CONFIG_SMALL("Interlace detect Filter."), | |
399 | .priv_size = sizeof(IDETContext), | |
400 | .init = init, | |
401 | .uninit = uninit, | |
402 | .query_formats = query_formats, | |
403 | .inputs = idet_inputs, | |
404 | .outputs = idet_outputs, | |
405 | .priv_class = &idet_class, | |
406 | }; |