annotate ffmpeg/libavfilter/vf_idet.c @ 13:844d341cf643 tip

Back up before ISMIR
author Yading Song <yading.song@eecs.qmul.ac.uk>
date Thu, 31 Oct 2013 13:17:06 +0000
parents 6840f77b83aa
children
rev   line source
yading@10 1 /*
yading@10 2 * Copyright (C) 2012 Michael Niedermayer <michaelni@gmx.at>
yading@10 3 *
yading@10 4 * This file is part of FFmpeg.
yading@10 5 *
yading@10 6 * FFmpeg is free software; you can redistribute it and/or
yading@10 7 * modify it under the terms of the GNU Lesser General Public
yading@10 8 * License as published by the Free Software Foundation; either
yading@10 9 * version 2.1 of the License, or (at your option) any later version.
yading@10 10 *
yading@10 11 * FFmpeg is distributed in the hope that it will be useful,
yading@10 12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@10 13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
yading@10 14 * Lesser General Public License for more details.
yading@10 15 *
yading@10 16 * You should have received a copy of the GNU Lesser General Public
yading@10 17 * License along with FFmpeg; if not, write to the Free Software
yading@10 18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@10 19 */
yading@10 20
yading@10 21 #include <float.h> /* FLT_MAX */
yading@10 22
yading@10 23 #include "libavutil/cpu.h"
yading@10 24 #include "libavutil/common.h"
yading@10 25 #include "libavutil/opt.h"
yading@10 26 #include "libavutil/pixdesc.h"
yading@10 27 #include "avfilter.h"
yading@10 28 #include "internal.h"
yading@10 29
yading@10 30 #define HIST_SIZE 4
yading@10 31
yading@10 32 typedef enum {
yading@10 33 TFF,
yading@10 34 BFF,
yading@10 35 PROGRSSIVE,
yading@10 36 UNDETERMINED,
yading@10 37 } Type;
yading@10 38
yading@10 39 typedef struct {
yading@10 40 const AVClass *class;
yading@10 41 float interlace_threshold;
yading@10 42 float progressive_threshold;
yading@10 43
yading@10 44 Type last_type;
yading@10 45 int prestat[4];
yading@10 46 int poststat[4];
yading@10 47
yading@10 48 uint8_t history[HIST_SIZE];
yading@10 49
yading@10 50 AVFrame *cur;
yading@10 51 AVFrame *next;
yading@10 52 AVFrame *prev;
yading@10 53 int (*filter_line)(const uint8_t *prev, const uint8_t *cur, const uint8_t *next, int w);
yading@10 54
yading@10 55 const AVPixFmtDescriptor *csp;
yading@10 56 } IDETContext;
yading@10 57
yading@10 58 #define OFFSET(x) offsetof(IDETContext, x)
yading@10 59 #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM
yading@10 60
yading@10 61 static const AVOption idet_options[] = {
yading@10 62 { "intl_thres", "set interlacing threshold", OFFSET(interlace_threshold), AV_OPT_TYPE_FLOAT, {.dbl = 1.04}, -1, FLT_MAX, FLAGS },
yading@10 63 { "prog_thres", "set progressive threshold", OFFSET(progressive_threshold), AV_OPT_TYPE_FLOAT, {.dbl = 1.5}, -1, FLT_MAX, FLAGS },
yading@10 64 { NULL }
yading@10 65 };
yading@10 66
yading@10 67 AVFILTER_DEFINE_CLASS(idet);
yading@10 68
yading@10 69 static const char *type2str(Type type)
yading@10 70 {
yading@10 71 switch(type) {
yading@10 72 case TFF : return "Top Field First ";
yading@10 73 case BFF : return "Bottom Field First";
yading@10 74 case PROGRSSIVE : return "Progressive ";
yading@10 75 case UNDETERMINED: return "Undetermined ";
yading@10 76 }
yading@10 77 return NULL;
yading@10 78 }
yading@10 79
yading@10 80 static int filter_line_c(const uint8_t *a, const uint8_t *b, const uint8_t *c, int w)
yading@10 81 {
yading@10 82 int x;
yading@10 83 int ret=0;
yading@10 84
yading@10 85 for(x=0; x<w; x++){
yading@10 86 int v = (*a++ + *c++) - 2 * *b++;
yading@10 87 ret += FFABS(v);
yading@10 88 }
yading@10 89
yading@10 90 return ret;
yading@10 91 }
yading@10 92
yading@10 93 static int filter_line_c_16bit(const uint16_t *a, const uint16_t *b, const uint16_t *c, int w)
yading@10 94 {
yading@10 95 int x;
yading@10 96 int ret=0;
yading@10 97
yading@10 98 for(x=0; x<w; x++){
yading@10 99 int v = (*a++ + *c++) - 2 * *b++;
yading@10 100 ret += FFABS(v);
yading@10 101 }
yading@10 102
yading@10 103 return ret;
yading@10 104 }
yading@10 105
yading@10 106 static void filter(AVFilterContext *ctx)
yading@10 107 {
yading@10 108 IDETContext *idet = ctx->priv;
yading@10 109 int y, i;
yading@10 110 int64_t alpha[2]={0};
yading@10 111 int64_t delta=0;
yading@10 112 Type type, best_type;
yading@10 113 int match = 0;
yading@10 114
yading@10 115 for (i = 0; i < idet->csp->nb_components; i++) {
yading@10 116 int w = idet->cur->width;
yading@10 117 int h = idet->cur->height;
yading@10 118 int refs = idet->cur->linesize[i];
yading@10 119
yading@10 120 if (i && i<3) {
yading@10 121 w >>= idet->csp->log2_chroma_w;
yading@10 122 h >>= idet->csp->log2_chroma_h;
yading@10 123 }
yading@10 124
yading@10 125 for (y = 2; y < h - 2; y++) {
yading@10 126 uint8_t *prev = &idet->prev->data[i][y*refs];
yading@10 127 uint8_t *cur = &idet->cur ->data[i][y*refs];
yading@10 128 uint8_t *next = &idet->next->data[i][y*refs];
yading@10 129 alpha[ y &1] += idet->filter_line(cur-refs, prev, cur+refs, w);
yading@10 130 alpha[(y^1)&1] += idet->filter_line(cur-refs, next, cur+refs, w);
yading@10 131 delta += idet->filter_line(cur-refs, cur, cur+refs, w);
yading@10 132 }
yading@10 133 }
yading@10 134
yading@10 135 if (alpha[0] > idet->interlace_threshold * alpha[1]){
yading@10 136 type = TFF;
yading@10 137 }else if(alpha[1] > idet->interlace_threshold * alpha[0]){
yading@10 138 type = BFF;
yading@10 139 }else if(alpha[1] > idet->progressive_threshold * delta){
yading@10 140 type = PROGRSSIVE;
yading@10 141 }else{
yading@10 142 type = UNDETERMINED;
yading@10 143 }
yading@10 144
yading@10 145 memmove(idet->history+1, idet->history, HIST_SIZE-1);
yading@10 146 idet->history[0] = type;
yading@10 147 best_type = UNDETERMINED;
yading@10 148 for(i=0; i<HIST_SIZE; i++){
yading@10 149 if(idet->history[i] != UNDETERMINED){
yading@10 150 if(best_type == UNDETERMINED)
yading@10 151 best_type = idet->history[i];
yading@10 152
yading@10 153 if(idet->history[i] == best_type) {
yading@10 154 match++;
yading@10 155 }else{
yading@10 156 match=0;
yading@10 157 break;
yading@10 158 }
yading@10 159 }
yading@10 160 }
yading@10 161 if(idet->last_type == UNDETERMINED){
yading@10 162 if(match ) idet->last_type = best_type;
yading@10 163 }else{
yading@10 164 if(match>2) idet->last_type = best_type;
yading@10 165 }
yading@10 166
yading@10 167 if (idet->last_type == TFF){
yading@10 168 idet->cur->top_field_first = 1;
yading@10 169 idet->cur->interlaced_frame = 1;
yading@10 170 }else if(idet->last_type == BFF){
yading@10 171 idet->cur->top_field_first = 0;
yading@10 172 idet->cur->interlaced_frame = 1;
yading@10 173 }else if(idet->last_type == PROGRSSIVE){
yading@10 174 idet->cur->interlaced_frame = 0;
yading@10 175 }
yading@10 176
yading@10 177 idet->prestat [ type] ++;
yading@10 178 idet->poststat[idet->last_type] ++;
yading@10 179 av_log(ctx, AV_LOG_DEBUG, "Single frame:%s, Multi frame:%s\n", type2str(type), type2str(idet->last_type));
yading@10 180 }
yading@10 181
yading@10 182 static int filter_frame(AVFilterLink *link, AVFrame *picref)
yading@10 183 {
yading@10 184 AVFilterContext *ctx = link->dst;
yading@10 185 IDETContext *idet = ctx->priv;
yading@10 186
yading@10 187 if (idet->prev)
yading@10 188 av_frame_free(&idet->prev);
yading@10 189 idet->prev = idet->cur;
yading@10 190 idet->cur = idet->next;
yading@10 191 idet->next = picref;
yading@10 192
yading@10 193 if (!idet->cur)
yading@10 194 return 0;
yading@10 195
yading@10 196 if (!idet->prev)
yading@10 197 idet->prev = av_frame_clone(idet->cur);
yading@10 198
yading@10 199 if (!idet->csp)
yading@10 200 idet->csp = av_pix_fmt_desc_get(link->format);
yading@10 201 if (idet->csp->comp[0].depth_minus1 / 8 == 1)
yading@10 202 idet->filter_line = (void*)filter_line_c_16bit;
yading@10 203
yading@10 204 filter(ctx);
yading@10 205
yading@10 206 return ff_filter_frame(ctx->outputs[0], av_frame_clone(idet->cur));
yading@10 207 }
yading@10 208
yading@10 209 static int request_frame(AVFilterLink *link)
yading@10 210 {
yading@10 211 AVFilterContext *ctx = link->src;
yading@10 212 IDETContext *idet = ctx->priv;
yading@10 213
yading@10 214 do {
yading@10 215 int ret;
yading@10 216
yading@10 217 if ((ret = ff_request_frame(link->src->inputs[0])))
yading@10 218 return ret;
yading@10 219 } while (!idet->cur);
yading@10 220
yading@10 221 return 0;
yading@10 222 }
yading@10 223
yading@10 224 static av_cold void uninit(AVFilterContext *ctx)
yading@10 225 {
yading@10 226 IDETContext *idet = ctx->priv;
yading@10 227
yading@10 228 av_log(ctx, AV_LOG_INFO, "Single frame detection: TFF:%d BFF:%d Progressive:%d Undetermined:%d\n",
yading@10 229 idet->prestat[TFF],
yading@10 230 idet->prestat[BFF],
yading@10 231 idet->prestat[PROGRSSIVE],
yading@10 232 idet->prestat[UNDETERMINED]
yading@10 233 );
yading@10 234 av_log(ctx, AV_LOG_INFO, "Multi frame detection: TFF:%d BFF:%d Progressive:%d Undetermined:%d\n",
yading@10 235 idet->poststat[TFF],
yading@10 236 idet->poststat[BFF],
yading@10 237 idet->poststat[PROGRSSIVE],
yading@10 238 idet->poststat[UNDETERMINED]
yading@10 239 );
yading@10 240
yading@10 241 av_frame_free(&idet->prev);
yading@10 242 av_frame_free(&idet->cur );
yading@10 243 av_frame_free(&idet->next);
yading@10 244 }
yading@10 245
yading@10 246 static int query_formats(AVFilterContext *ctx)
yading@10 247 {
yading@10 248 static const enum AVPixelFormat pix_fmts[] = {
yading@10 249 AV_PIX_FMT_YUV420P,
yading@10 250 AV_PIX_FMT_YUV422P,
yading@10 251 AV_PIX_FMT_YUV444P,
yading@10 252 AV_PIX_FMT_YUV410P,
yading@10 253 AV_PIX_FMT_YUV411P,
yading@10 254 AV_PIX_FMT_GRAY8,
yading@10 255 AV_PIX_FMT_YUVJ420P,
yading@10 256 AV_PIX_FMT_YUVJ422P,
yading@10 257 AV_PIX_FMT_YUVJ444P,
yading@10 258 AV_NE( AV_PIX_FMT_GRAY16BE, AV_PIX_FMT_GRAY16LE ),
yading@10 259 AV_PIX_FMT_YUV440P,
yading@10 260 AV_PIX_FMT_YUVJ440P,
yading@10 261 AV_NE( AV_PIX_FMT_YUV420P10BE, AV_PIX_FMT_YUV420P10LE ),
yading@10 262 AV_NE( AV_PIX_FMT_YUV422P10BE, AV_PIX_FMT_YUV422P10LE ),
yading@10 263 AV_NE( AV_PIX_FMT_YUV444P10BE, AV_PIX_FMT_YUV444P10LE ),
yading@10 264 AV_NE( AV_PIX_FMT_YUV420P16BE, AV_PIX_FMT_YUV420P16LE ),
yading@10 265 AV_NE( AV_PIX_FMT_YUV422P16BE, AV_PIX_FMT_YUV422P16LE ),
yading@10 266 AV_NE( AV_PIX_FMT_YUV444P16BE, AV_PIX_FMT_YUV444P16LE ),
yading@10 267 AV_PIX_FMT_YUVA420P,
yading@10 268 AV_PIX_FMT_NONE
yading@10 269 };
yading@10 270
yading@10 271 ff_set_common_formats(ctx, ff_make_format_list(pix_fmts));
yading@10 272
yading@10 273 return 0;
yading@10 274 }
yading@10 275
yading@10 276 static av_cold int init(AVFilterContext *ctx)
yading@10 277 {
yading@10 278 IDETContext *idet = ctx->priv;
yading@10 279
yading@10 280 idet->last_type = UNDETERMINED;
yading@10 281 memset(idet->history, UNDETERMINED, HIST_SIZE);
yading@10 282
yading@10 283 idet->filter_line = filter_line_c;
yading@10 284
yading@10 285 return 0;
yading@10 286 }
yading@10 287
yading@10 288
yading@10 289 static const AVFilterPad idet_inputs[] = {
yading@10 290 {
yading@10 291 .name = "default",
yading@10 292 .type = AVMEDIA_TYPE_VIDEO,
yading@10 293 .filter_frame = filter_frame,
yading@10 294 },
yading@10 295 { NULL }
yading@10 296 };
yading@10 297
yading@10 298 static const AVFilterPad idet_outputs[] = {
yading@10 299 {
yading@10 300 .name = "default",
yading@10 301 .type = AVMEDIA_TYPE_VIDEO,
yading@10 302 .request_frame = request_frame,
yading@10 303 },
yading@10 304 { NULL }
yading@10 305 };
yading@10 306
yading@10 307 AVFilter avfilter_vf_idet = {
yading@10 308 .name = "idet",
yading@10 309 .description = NULL_IF_CONFIG_SMALL("Interlace detect Filter."),
yading@10 310
yading@10 311 .priv_size = sizeof(IDETContext),
yading@10 312 .init = init,
yading@10 313 .uninit = uninit,
yading@10 314 .query_formats = query_formats,
yading@10 315 .inputs = idet_inputs,
yading@10 316 .outputs = idet_outputs,
yading@10 317 .priv_class = &idet_class,
yading@10 318 };