2 * Copyright (c) 2003 Rich Felker
4 * This file is part of FFmpeg.
6 * FFmpeg is free software; you can redistribute it and/or
7 * modify it under the terms of the GNU General Public
8 * License as published by the Free Software Foundation; either
9 * version 2 of the License, or (at your option) any later version.
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
14 * GNU General Public License for more details.
16 * You should have received a copy of the GNU General Public License along
17 * with FFmpeg; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 #include "libavutil/avassert.h"
22 #include "libavutil/imgutils.h"
23 #include "libavutil/opt.h"
24 #include "libavutil/pixdesc.h"
29 #include "vf_pullup.h"
31 #define F_HAVE_BREAKS 1
32 #define F_HAVE_AFFINITY 2
37 #define OFFSET(x) offsetof(PullupContext, x)
38 #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
40 static const AVOption pullup_options[] = {
41 { "jl", "set left junk size", OFFSET(junk_left), AV_OPT_TYPE_INT, {.i64=1}, 0, INT_MAX, FLAGS },
42 { "jr", "set right junk size", OFFSET(junk_right), AV_OPT_TYPE_INT, {.i64=1}, 0, INT_MAX, FLAGS },
43 { "jt", "set top junk size", OFFSET(junk_top), AV_OPT_TYPE_INT, {.i64=4}, 1, INT_MAX, FLAGS },
44 { "jb", "set bottom junk size", OFFSET(junk_bottom), AV_OPT_TYPE_INT, {.i64=4}, 1, INT_MAX, FLAGS },
45 { "sb", "set strict breaks", OFFSET(strict_breaks), AV_OPT_TYPE_BOOL,{.i64=0},-1, 1, FLAGS },
46 { "mp", "set metric plane", OFFSET(metric_plane), AV_OPT_TYPE_INT, {.i64=0}, 0, 2, FLAGS, "mp" },
47 { "y", "luma", 0, AV_OPT_TYPE_CONST, {.i64=0}, 0, 0, FLAGS, "mp" },
48 { "u", "chroma blue", 0, AV_OPT_TYPE_CONST, {.i64=1}, 0, 0, FLAGS, "mp" },
49 { "v", "chroma red", 0, AV_OPT_TYPE_CONST, {.i64=2}, 0, 0, FLAGS, "mp" },
53 AVFILTER_DEFINE_CLASS(pullup);
55 static int query_formats(AVFilterContext *ctx)
57 static const enum AVPixelFormat pix_fmts[] = {
58 AV_PIX_FMT_YUVJ444P, AV_PIX_FMT_YUVJ440P,
59 AV_PIX_FMT_YUVJ422P, AV_PIX_FMT_YUVJ420P,
60 AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUV440P,
61 AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUV420P,
62 AV_PIX_FMT_YUV411P, AV_PIX_FMT_YUV410P,
63 AV_PIX_FMT_YUVJ411P, AV_PIX_FMT_GRAY8,
67 AVFilterFormats *fmts_list = ff_make_format_list(pix_fmts);
69 return AVERROR(ENOMEM);
70 return ff_set_common_formats(ctx, fmts_list);
73 #define ABS(a) (((a) ^ ((a) >> 31)) - ((a) >> 31))
75 static int diff_c(const uint8_t *a, const uint8_t *b, ptrdiff_t s)
79 for (i = 0; i < 4; i++) {
80 for (j = 0; j < 8; j++)
81 diff += ABS(a[j] - b[j]);
89 static int comb_c(const uint8_t *a, const uint8_t *b, ptrdiff_t s)
93 for (i = 0; i < 4; i++) {
94 for (j = 0; j < 8; j++)
95 comb += ABS((a[j] << 1) - b[j - s] - b[j ]) +
96 ABS((b[j] << 1) - a[j ] - a[j + s]);
104 static int var_c(const uint8_t *a, const uint8_t *b, ptrdiff_t s)
108 for (i = 0; i < 3; i++) {
109 for (j = 0; j < 8; j++)
110 var += ABS(a[j] - a[j + s]);
114 return 4 * var; /* match comb scaling */
117 static int alloc_metrics(PullupContext *s, PullupField *f)
119 f->diffs = av_calloc(FFALIGN(s->metric_length, 16), sizeof(*f->diffs));
120 f->combs = av_calloc(FFALIGN(s->metric_length, 16), sizeof(*f->combs));
121 f->vars = av_calloc(FFALIGN(s->metric_length, 16), sizeof(*f->vars));
123 if (!f->diffs || !f->combs || !f->vars) {
127 return AVERROR(ENOMEM);
132 static void free_field_queue(PullupField *head)
134 PullupField *f = head;
143 memset(f, 0, sizeof(*f)); // clear all pointers to avoid stale ones
149 static PullupField *make_field_queue(PullupContext *s, int len)
151 PullupField *head, *f;
153 f = head = av_mallocz(sizeof(*head));
157 if (alloc_metrics(s, f) < 0) {
162 for (; len > 0; len--) {
163 f->next = av_mallocz(sizeof(*f->next));
165 free_field_queue(head);
171 if (alloc_metrics(s, f) < 0) {
172 free_field_queue(head);
183 static int config_input(AVFilterLink *inlink)
185 AVFilterContext *ctx = inlink->dst;
186 PullupContext *s = ctx->priv;
187 const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
188 int mp = s->metric_plane;
190 s->nb_planes = av_pix_fmt_count_planes(inlink->format);
192 if (mp + 1 > s->nb_planes) {
193 av_log(ctx, AV_LOG_ERROR, "input format does not have such plane\n");
194 return AVERROR(EINVAL);
197 s->planeheight[1] = s->planeheight[2] = AV_CEIL_RSHIFT(inlink->h, desc->log2_chroma_h);
198 s->planeheight[0] = s->planeheight[3] = inlink->h;
199 s->planewidth[1] = s->planewidth[2] = AV_CEIL_RSHIFT(inlink->w, desc->log2_chroma_w);
200 s->planewidth[0] = s->planewidth[3] = inlink->w;
202 s->metric_w = (s->planewidth[mp] - ((s->junk_left + s->junk_right) << 3)) >> 3;
203 s->metric_h = (s->planeheight[mp] - ((s->junk_top + s->junk_bottom) << 1)) >> 3;
204 s->metric_offset = (s->junk_left << 3) + (s->junk_top << 1) * s->planewidth[mp];
205 s->metric_length = s->metric_w * s->metric_h;
207 av_log(ctx, AV_LOG_DEBUG, "w: %d h: %d\n", s->metric_w, s->metric_h);
208 av_log(ctx, AV_LOG_DEBUG, "offset: %d length: %d\n", s->metric_offset, s->metric_length);
210 s->head = make_field_queue(s, 8);
212 return AVERROR(ENOMEM);
219 ff_pullup_init_x86(s);
223 static PullupBuffer *pullup_lock_buffer(PullupBuffer *b, int parity)
228 if ((parity + 1) & 1)
230 if ((parity + 1) & 2)
236 static void pullup_release_buffer(PullupBuffer *b, int parity)
241 if ((parity + 1) & 1)
243 if ((parity + 1) & 2)
247 static int alloc_buffer(PullupContext *s, PullupBuffer *b)
253 for (i = 0; i < s->nb_planes; i++) {
254 b->planes[i] = av_malloc(s->planeheight[i] * s->planewidth[i]);
256 if (s->nb_planes == 1)
257 b->planes[1] = av_malloc(4*256);
262 static PullupBuffer *pullup_get_buffer(PullupContext *s, int parity)
266 /* Try first to get the sister buffer for the previous field */
267 if (parity < 2 && s->last && parity != s->last->parity
268 && !s->last->buffer->lock[parity]) {
269 alloc_buffer(s, s->last->buffer);
270 return pullup_lock_buffer(s->last->buffer, parity);
273 /* Prefer a buffer with both fields open */
274 for (i = 0; i < FF_ARRAY_ELEMS(s->buffers); i++) {
275 if (s->buffers[i].lock[0])
277 if (s->buffers[i].lock[1])
279 alloc_buffer(s, &s->buffers[i]);
280 return pullup_lock_buffer(&s->buffers[i], parity);
286 /* Search for any half-free buffer */
287 for (i = 0; i < FF_ARRAY_ELEMS(s->buffers); i++) {
288 if (((parity + 1) & 1) && s->buffers[i].lock[0])
290 if (((parity + 1) & 2) && s->buffers[i].lock[1])
292 alloc_buffer(s, &s->buffers[i]);
293 return pullup_lock_buffer(&s->buffers[i], parity);
299 static int queue_length(PullupField *begin, PullupField *end)
307 for (f = begin; f != end; f = f->next)
313 static int find_first_break(PullupField *f, int max)
317 for (i = 0; i < max; i++) {
318 if (f->breaks & BREAK_RIGHT || f->next->breaks & BREAK_LEFT)
326 static void compute_breaks(PullupContext *s, PullupField *f0)
328 PullupField *f1 = f0->next;
329 PullupField *f2 = f1->next;
330 PullupField *f3 = f2->next;
331 int i, l, max_l = 0, max_r = 0;
333 if (f0->flags & F_HAVE_BREAKS)
336 f0->flags |= F_HAVE_BREAKS;
338 /* Special case when fields are 100% identical */
339 if (f0->buffer == f2->buffer && f1->buffer != f3->buffer) {
340 f2->breaks |= BREAK_RIGHT;
344 if (f0->buffer != f2->buffer && f1->buffer == f3->buffer) {
345 f1->breaks |= BREAK_LEFT;
349 for (i = 0; i < s->metric_length; i++) {
350 l = f2->diffs[i] - f3->diffs[i];
358 /* Don't get tripped up when differences are mostly quant error */
359 if (max_l + max_r < 128)
361 if (max_l > 4 * max_r)
362 f1->breaks |= BREAK_LEFT;
363 if (max_r > 4 * max_l)
364 f2->breaks |= BREAK_RIGHT;
367 static void compute_affinity(PullupContext *s, PullupField *f)
369 int i, max_l = 0, max_r = 0, l;
371 if (f->flags & F_HAVE_AFFINITY)
374 f->flags |= F_HAVE_AFFINITY;
376 if (f->buffer == f->next->next->buffer) {
378 f->next->affinity = 0;
379 f->next->next->affinity = -1;
380 f->next->flags |= F_HAVE_AFFINITY;
381 f->next->next->flags |= F_HAVE_AFFINITY;
385 for (i = 0; i < s->metric_length; i++) {
387 int lv = f->prev->vars[i];
388 int rv = f->next->vars[i];
389 int lc = f-> combs[i] - 2*(v < lv ? v : lv);
390 int rc = f->next->combs[i] - 2*(v < rv ? v : rv);
402 if (max_l + max_r < 64)
405 if (max_r > 6 * max_l)
407 else if (max_l > 6 * max_r)
411 static int decide_frame_length(PullupContext *s)
413 PullupField *f0 = s->first;
414 PullupField *f1 = f0->next;
415 PullupField *f2 = f1->next;
419 if (queue_length(s->first, s->last) < 4)
423 n = queue_length(f, s->last);
424 for (i = 0; i < n - 1; i++) {
426 compute_breaks(s, f);
428 compute_affinity(s, f);
433 if (f0->affinity == -1)
436 l = find_first_break(f0, 3);
438 if (l == 1 && s->strict_breaks < 0)
443 return 1 + (s->strict_breaks < 1 && f0->affinity == 1 && f1->affinity == -1);
445 /* FIXME: strictly speaking, f0->prev is no longer valid... :) */
447 && (f0->prev->breaks & BREAK_RIGHT) && (f2->breaks & BREAK_LEFT)
448 && (f0->affinity != 1 || f1->affinity != -1) )
450 return 1 + (f1->affinity != 1);
452 return 2 + (f2->affinity != 1);
454 /* 9 possibilities covered before switch */
455 if (f1->affinity == 1)
456 return 1; /* covers 6 */
457 else if (f1->affinity == -1)
458 return 2; /* covers 6 */
459 else if (f2->affinity == -1) { /* covers 2 */
460 return (f0->affinity == 1) ? 3 : 1;
462 return 2; /* the remaining 6 */
467 static PullupFrame *pullup_get_frame(PullupContext *s)
469 PullupFrame *fr = &s->frame;
470 int i, n = decide_frame_length(s);
471 int aff = s->first->next->affinity;
473 av_assert1(n < FF_ARRAY_ELEMS(fr->ifields));
479 fr->parity = s->first->parity;
482 for (i = 0; i < n; i++) {
483 /* We cheat and steal the buffer without release+relock */
484 fr->ifields[i] = s->first->buffer;
485 s->first->buffer = 0;
486 s->first = s->first->next;
490 fr->ofields[fr->parity ] = fr->ifields[0];
491 fr->ofields[fr->parity ^ 1] = 0;
493 fr->ofields[fr->parity ] = fr->ifields[0];
494 fr->ofields[fr->parity ^ 1] = fr->ifields[1];
497 aff = (fr->ifields[0] == fr->ifields[1]) ? -1 : 1;
498 fr->ofields[fr->parity ] = fr->ifields[1 + aff];
499 fr->ofields[fr->parity ^ 1] = fr->ifields[1 ];
502 pullup_lock_buffer(fr->ofields[0], 0);
503 pullup_lock_buffer(fr->ofields[1], 1);
505 if (fr->ofields[0] == fr->ofields[1]) {
506 fr->buffer = fr->ofields[0];
507 pullup_lock_buffer(fr->buffer, 2);
514 static void pullup_release_frame(PullupFrame *f)
518 for (i = 0; i < f->length; i++)
519 pullup_release_buffer(f->ifields[i], f->parity ^ (i & 1));
521 pullup_release_buffer(f->ofields[0], 0);
522 pullup_release_buffer(f->ofields[1], 1);
525 pullup_release_buffer(f->buffer, 2);
529 static void compute_metric(PullupContext *s, int *dest,
530 PullupField *fa, int pa, PullupField *fb, int pb,
531 int (*func)(const uint8_t *, const uint8_t *, ptrdiff_t))
533 int mp = s->metric_plane;
535 int ystep = s->planewidth[mp] << 3;
536 int stride = s->planewidth[mp] << 1; /* field stride */
537 int w = s->metric_w * xstep;
541 if (!fa->buffer || !fb->buffer)
544 /* Shortcut for duplicate fields (e.g. from RFF flag) */
545 if (fa->buffer == fb->buffer && pa == pb) {
546 memset(dest, 0, s->metric_length * sizeof(*dest));
550 a = fa->buffer->planes[mp] + pa * s->planewidth[mp] + s->metric_offset;
551 b = fb->buffer->planes[mp] + pb * s->planewidth[mp] + s->metric_offset;
553 for (y = 0; y < s->metric_h; y++) {
554 for (x = 0; x < w; x += xstep)
555 *dest++ = func(a + x, b + x, stride);
556 a += ystep; b += ystep;
560 static int check_field_queue(PullupContext *s)
564 if (s->head->next == s->first) {
565 PullupField *f = av_mallocz(sizeof(*f));
568 return AVERROR(ENOMEM);
570 if ((ret = alloc_metrics(s, f)) < 0) {
584 static void pullup_submit_field(PullupContext *s, PullupBuffer *b, int parity)
588 /* Grow the circular list if needed */
589 if (check_field_queue(s) < 0)
592 /* Cannot have two fields of same parity in a row; drop the new one */
593 if (s->last && s->last->parity == parity)
598 f->buffer = pullup_lock_buffer(b, parity);
603 compute_metric(s, f->diffs, f, parity, f->prev->prev, parity, s->diff);
604 compute_metric(s, f->combs, parity ? f->prev : f, 0, parity ? f : f->prev, 1, s->comb);
605 compute_metric(s, f->vars, f, parity, f, -1, s->var);
608 /* Advance the circular list */
613 s->head = s->head->next;
616 static void copy_field(PullupContext *s,
617 PullupBuffer *dst, PullupBuffer *src, int parity)
622 for (i = 0; i < s->nb_planes; i++) {
623 ss = src->planes[i] + parity * s->planewidth[i];
624 dd = dst->planes[i] + parity * s->planewidth[i];
626 av_image_copy_plane(dd, s->planewidth[i] << 1,
627 ss, s->planewidth[i] << 1,
628 s->planewidth[i], s->planeheight[i] >> 1);
632 static void pullup_pack_frame(PullupContext *s, PullupFrame *fr)
640 return; /* FIXME: deal with this */
642 for (i = 0; i < 2; i++) {
643 if (fr->ofields[i]->lock[i^1])
646 fr->buffer = fr->ofields[i];
647 pullup_lock_buffer(fr->buffer, 2);
648 copy_field(s, fr->buffer, fr->ofields[i^1], i^1);
652 fr->buffer = pullup_get_buffer(s, 2);
654 copy_field(s, fr->buffer, fr->ofields[0], 0);
655 copy_field(s, fr->buffer, fr->ofields[1], 1);
658 static int filter_frame(AVFilterLink *inlink, AVFrame *in)
660 AVFilterContext *ctx = inlink->dst;
661 AVFilterLink *outlink = ctx->outputs[0];
662 PullupContext *s = ctx->priv;
668 b = pullup_get_buffer(s, 2);
670 av_log(ctx, AV_LOG_WARNING, "Could not get buffer!\n");
671 f = pullup_get_frame(s);
672 pullup_release_frame(f);
676 av_image_copy(b->planes, s->planewidth,
677 (const uint8_t**)in->data, in->linesize,
678 inlink->format, inlink->w, inlink->h);
680 p = in->interlaced_frame ? !in->top_field_first : 0;
681 pullup_submit_field(s, b, p );
682 pullup_submit_field(s, b, p^1);
685 pullup_submit_field(s, b, p);
687 pullup_release_buffer(b, 2);
689 f = pullup_get_frame(s);
694 pullup_release_frame(f);
695 f = pullup_get_frame(s);
699 pullup_release_frame(f);
700 if (!in->repeat_pict)
702 f = pullup_get_frame(s);
706 pullup_release_frame(f);
712 /* If the frame isn't already exportable... */
714 pullup_pack_frame(s, f);
716 out = ff_get_video_buffer(outlink, outlink->w, outlink->h);
718 ret = AVERROR(ENOMEM);
721 av_frame_copy_props(out, in);
723 av_image_copy(out->data, out->linesize,
724 (const uint8_t**)f->buffer->planes, s->planewidth,
725 inlink->format, inlink->w, inlink->h);
727 ret = ff_filter_frame(outlink, out);
728 pullup_release_frame(f);
734 static av_cold void uninit(AVFilterContext *ctx)
736 PullupContext *s = ctx->priv;
739 free_field_queue(s->head);
742 for (i = 0; i < FF_ARRAY_ELEMS(s->buffers); i++) {
743 av_freep(&s->buffers[i].planes[0]);
744 av_freep(&s->buffers[i].planes[1]);
745 av_freep(&s->buffers[i].planes[2]);
749 static const AVFilterPad pullup_inputs[] = {
752 .type = AVMEDIA_TYPE_VIDEO,
753 .filter_frame = filter_frame,
754 .config_props = config_input,
759 static const AVFilterPad pullup_outputs[] = {
762 .type = AVMEDIA_TYPE_VIDEO,
767 const AVFilter ff_vf_pullup = {
769 .description = NULL_IF_CONFIG_SMALL("Pullup from field sequence to frames."),
770 .priv_size = sizeof(PullupContext),
771 .priv_class = &pullup_class,
773 .query_formats = query_formats,
774 .inputs = pullup_inputs,
775 .outputs = pullup_outputs,