| /* |
| * Copyright (c) 2019 Vladimir Panteleev |
| * |
| * This file is part of FFmpeg. |
| * |
| * FFmpeg is free software; you can redistribute it and/or |
| * modify it under the terms of the GNU Lesser General Public |
| * License as published by the Free Software Foundation; either |
| * version 2.1 of the License, or (at your option) any later version. |
| * |
| * FFmpeg is distributed in the hope that it will be useful, |
| * but WITHOUT ANY WARRANTY; without even the implied warranty of |
| * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| * Lesser General Public License for more details. |
| * |
| * You should have received a copy of the GNU Lesser General Public |
| * License along with FFmpeg; if not, write to the Free Software |
| * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
| */ |
| |
| #include <float.h> |
| |
| #include "libavutil/imgutils.h" |
| #include "libavutil/opt.h" |
| #include "libavutil/pixdesc.h" |
| #include "avfilter.h" |
| |
| #include "formats.h" |
| #include "internal.h" |
| #include "video.h" |
| |
| #define MAX_FRAMES 240 |
| #define GRID_SIZE 8 |
| #define NUM_CHANNELS 3 |
| |
| typedef struct PhotosensitivityFrame { |
| uint8_t grid[GRID_SIZE][GRID_SIZE][4]; |
| } PhotosensitivityFrame; |
| |
| typedef struct PhotosensitivityContext { |
| const AVClass *class; |
| |
| int nb_frames; |
| int skip; |
| float threshold_multiplier; |
| int bypass; |
| |
| int badness_threshold; |
| |
| /* Circular buffer */ |
| int history[MAX_FRAMES]; |
| int history_pos; |
| |
| PhotosensitivityFrame last_frame_e; |
| AVFrame *last_frame_av; |
| } PhotosensitivityContext; |
| |
| #define OFFSET(x) offsetof(PhotosensitivityContext, x) |
| #define FLAGS AV_OPT_FLAG_VIDEO_PARAM|AV_OPT_FLAG_FILTERING_PARAM |
| |
| static const AVOption photosensitivity_options[] = { |
| { "frames", "set how many frames to use", OFFSET(nb_frames), AV_OPT_TYPE_INT, {.i64=30}, 2, MAX_FRAMES, FLAGS }, |
| { "f", "set how many frames to use", OFFSET(nb_frames), AV_OPT_TYPE_INT, {.i64=30}, 2, MAX_FRAMES, FLAGS }, |
| { "threshold", "set detection threshold factor (lower is stricter)", OFFSET(threshold_multiplier), AV_OPT_TYPE_FLOAT, {.dbl=1}, 0.1, FLT_MAX, FLAGS }, |
| { "t", "set detection threshold factor (lower is stricter)", OFFSET(threshold_multiplier), AV_OPT_TYPE_FLOAT, {.dbl=1}, 0.1, FLT_MAX, FLAGS }, |
| { "skip", "set pixels to skip when sampling frames", OFFSET(skip), AV_OPT_TYPE_INT, {.i64=1}, 1, 1024, FLAGS }, |
| { "bypass", "leave frames unchanged", OFFSET(bypass), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1, FLAGS }, |
| { NULL } |
| }; |
| |
| AVFILTER_DEFINE_CLASS(photosensitivity); |
| |
| static int query_formats(AVFilterContext *ctx) |
| { |
| static const enum AVPixelFormat pixel_fmts[] = { |
| AV_PIX_FMT_RGB24, |
| AV_PIX_FMT_BGR24, |
| AV_PIX_FMT_NONE |
| }; |
| AVFilterFormats *formats = ff_make_format_list(pixel_fmts); |
| if (!formats) |
| return AVERROR(ENOMEM); |
| return ff_set_common_formats(ctx, formats); |
| } |
| |
| typedef struct ThreadData_convert_frame |
| { |
| AVFrame *in; |
| PhotosensitivityFrame *out; |
| int skip; |
| } ThreadData_convert_frame; |
| |
| #define NUM_CELLS (GRID_SIZE * GRID_SIZE) |
| |
| static int convert_frame_partial(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs) |
| { |
| int cell, gx, gy, x0, x1, y0, y1, x, y, c, area; |
| int sum[NUM_CHANNELS]; |
| const uint8_t *p; |
| |
| ThreadData_convert_frame *td = arg; |
| |
| const int slice_start = (NUM_CELLS * jobnr) / nb_jobs; |
| const int slice_end = (NUM_CELLS * (jobnr+1)) / nb_jobs; |
| |
| int width = td->in->width, height = td->in->height, linesize = td->in->linesize[0], skip = td->skip; |
| const uint8_t *data = td->in->data[0]; |
| |
| for (cell = slice_start; cell < slice_end; cell++) { |
| gx = cell % GRID_SIZE; |
| gy = cell / GRID_SIZE; |
| |
| x0 = width * gx / GRID_SIZE; |
| x1 = width * (gx+1) / GRID_SIZE; |
| y0 = height * gy / GRID_SIZE; |
| y1 = height * (gy+1) / GRID_SIZE; |
| |
| for (c = 0; c < NUM_CHANNELS; c++) { |
| sum[c] = 0; |
| } |
| for (y = y0; y < y1; y += skip) { |
| p = data + y * linesize + x0 * NUM_CHANNELS; |
| for (x = x0; x < x1; x += skip) { |
| //av_log(NULL, AV_LOG_VERBOSE, "%d %d %d : (%d,%d) (%d,%d) -> %d,%d | *%d\n", c, gx, gy, x0, y0, x1, y1, x, y, (int)row); |
| sum[0] += p[0]; |
| sum[1] += p[1]; |
| sum[2] += p[2]; |
| p += NUM_CHANNELS * skip; |
| // TODO: variable size |
| } |
| } |
| |
| area = ((x1 - x0 + skip - 1) / skip) * ((y1 - y0 + skip - 1) / skip); |
| for (c = 0; c < NUM_CHANNELS; c++) { |
| if (area) |
| sum[c] /= area; |
| td->out->grid[gy][gx][c] = sum[c]; |
| } |
| } |
| return 0; |
| } |
| |
| static void convert_frame(AVFilterContext *ctx, AVFrame *in, PhotosensitivityFrame *out, int skip) |
| { |
| ThreadData_convert_frame td; |
| td.in = in; |
| td.out = out; |
| td.skip = skip; |
| ctx->internal->execute(ctx, convert_frame_partial, &td, NULL, FFMIN(NUM_CELLS, ff_filter_get_nb_threads(ctx))); |
| } |
| |
| typedef struct ThreadData_blend_frame |
| { |
| AVFrame *target; |
| AVFrame *source; |
| uint16_t s_mul; |
| } ThreadData_blend_frame; |
| |
| static int blend_frame_partial(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs) |
| { |
| int x, y; |
| uint8_t *t, *s; |
| |
| ThreadData_blend_frame *td = arg; |
| const uint16_t s_mul = td->s_mul; |
| const uint16_t t_mul = 0x100 - s_mul; |
| const int slice_start = (td->target->height * jobnr) / nb_jobs; |
| const int slice_end = (td->target->height * (jobnr+1)) / nb_jobs; |
| const int linesize = td->target->linesize[0]; |
| |
| for (y = slice_start; y < slice_end; y++) { |
| t = td->target->data[0] + y * td->target->linesize[0]; |
| s = td->source->data[0] + y * td->source->linesize[0]; |
| for (x = 0; x < linesize; x++) { |
| *t = (*t * t_mul + *s * s_mul) >> 8; |
| t++; s++; |
| } |
| } |
| return 0; |
| } |
| |
| static void blend_frame(AVFilterContext *ctx, AVFrame *target, AVFrame *source, float factor) |
| { |
| ThreadData_blend_frame td; |
| td.target = target; |
| td.source = source; |
| td.s_mul = (uint16_t)(factor * 0x100); |
| ctx->internal->execute(ctx, blend_frame_partial, &td, NULL, FFMIN(ctx->outputs[0]->h, ff_filter_get_nb_threads(ctx))); |
| } |
| |
| static int get_badness(PhotosensitivityFrame *a, PhotosensitivityFrame *b) |
| { |
| int badness, x, y, c; |
| badness = 0; |
| for (c = 0; c < NUM_CHANNELS; c++) { |
| for (y = 0; y < GRID_SIZE; y++) { |
| for (x = 0; x < GRID_SIZE; x++) { |
| badness += abs((int)a->grid[y][x][c] - (int)b->grid[y][x][c]); |
| //av_log(NULL, AV_LOG_VERBOSE, "%d - %d -> %d \n", a->grid[y][x], b->grid[y][x], badness); |
| //av_log(NULL, AV_LOG_VERBOSE, "%d -> %d \n", abs((int)a->grid[y][x] - (int)b->grid[y][x]), badness); |
| } |
| } |
| } |
| return badness; |
| } |
| |
| static int config_input(AVFilterLink *inlink) |
| { |
| /* const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format); */ |
| AVFilterContext *ctx = inlink->dst; |
| PhotosensitivityContext *s = ctx->priv; |
| |
| s->badness_threshold = (int)(GRID_SIZE * GRID_SIZE * 4 * 256 * s->nb_frames * s->threshold_multiplier / 128); |
| |
| return 0; |
| } |
| |
| static int filter_frame(AVFilterLink *inlink, AVFrame *in) |
| { |
| int this_badness, current_badness, fixed_badness, new_badness, i, res; |
| PhotosensitivityFrame ef; |
| AVFrame *src, *out; |
| int free_in = 0; |
| float factor; |
| AVDictionary **metadata; |
| |
| AVFilterContext *ctx = inlink->dst; |
| AVFilterLink *outlink = ctx->outputs[0]; |
| PhotosensitivityContext *s = ctx->priv; |
| |
| /* weighted moving average */ |
| current_badness = 0; |
| for (i = 1; i < s->nb_frames; i++) |
| current_badness += i * s->history[(s->history_pos + i) % s->nb_frames]; |
| current_badness /= s->nb_frames; |
| |
| convert_frame(ctx, in, &ef, s->skip); |
| this_badness = get_badness(&ef, &s->last_frame_e); |
| new_badness = current_badness + this_badness; |
| av_log(s, AV_LOG_VERBOSE, "badness: %6d -> %6d / %6d (%3d%% - %s)\n", |
| current_badness, new_badness, s->badness_threshold, |
| 100 * new_badness / s->badness_threshold, new_badness < s->badness_threshold ? "OK" : "EXCEEDED"); |
| |
| fixed_badness = new_badness; |
| if (new_badness < s->badness_threshold || !s->last_frame_av || s->bypass) { |
| factor = 1; /* for metadata */ |
| av_frame_free(&s->last_frame_av); |
| s->last_frame_av = src = in; |
| s->last_frame_e = ef; |
| s->history[s->history_pos] = this_badness; |
| } else { |
| factor = (float)(s->badness_threshold - current_badness) / (new_badness - current_badness); |
| if (factor <= 0) { |
| /* just duplicate the frame */ |
| s->history[s->history_pos] = 0; /* frame was duplicated, thus, delta is zero */ |
| } else { |
| res = av_frame_make_writable(s->last_frame_av); |
| if (res) { |
| av_frame_free(&in); |
| return res; |
| } |
| blend_frame(ctx, s->last_frame_av, in, factor); |
| |
| convert_frame(ctx, s->last_frame_av, &ef, s->skip); |
| this_badness = get_badness(&ef, &s->last_frame_e); |
| fixed_badness = current_badness + this_badness; |
| av_log(s, AV_LOG_VERBOSE, " fixed: %6d -> %6d / %6d (%3d%%) factor=%5.3f\n", |
| current_badness, fixed_badness, s->badness_threshold, |
| 100 * new_badness / s->badness_threshold, factor); |
| s->last_frame_e = ef; |
| s->history[s->history_pos] = this_badness; |
| } |
| src = s->last_frame_av; |
| free_in = 1; |
| } |
| s->history_pos = (s->history_pos + 1) % s->nb_frames; |
| |
| out = ff_get_video_buffer(outlink, in->width, in->height); |
| if (!out) { |
| if (free_in == 1) |
| av_frame_free(&in); |
| return AVERROR(ENOMEM); |
| } |
| av_frame_copy_props(out, in); |
| metadata = &out->metadata; |
| if (metadata) { |
| char value[128]; |
| |
| snprintf(value, sizeof(value), "%f", (float)new_badness / s->badness_threshold); |
| av_dict_set(metadata, "lavfi.photosensitivity.badness", value, 0); |
| |
| snprintf(value, sizeof(value), "%f", (float)fixed_badness / s->badness_threshold); |
| av_dict_set(metadata, "lavfi.photosensitivity.fixed-badness", value, 0); |
| |
| snprintf(value, sizeof(value), "%f", (float)this_badness / s->badness_threshold); |
| av_dict_set(metadata, "lavfi.photosensitivity.frame-badness", value, 0); |
| |
| snprintf(value, sizeof(value), "%f", factor); |
| av_dict_set(metadata, "lavfi.photosensitivity.factor", value, 0); |
| } |
| av_frame_copy(out, src); |
| if (free_in == 1) |
| av_frame_free(&in); |
| return ff_filter_frame(outlink, out); |
| } |
| |
| static av_cold void uninit(AVFilterContext *ctx) |
| { |
| PhotosensitivityContext *s = ctx->priv; |
| |
| av_frame_free(&s->last_frame_av); |
| } |
| |
| static const AVFilterPad inputs[] = { |
| { |
| .name = "default", |
| .type = AVMEDIA_TYPE_VIDEO, |
| .filter_frame = filter_frame, |
| .config_props = config_input, |
| }, |
| { NULL } |
| }; |
| |
| static const AVFilterPad outputs[] = { |
| { |
| .name = "default", |
| .type = AVMEDIA_TYPE_VIDEO, |
| }, |
| { NULL } |
| }; |
| |
| AVFilter ff_vf_photosensitivity = { |
| .name = "photosensitivity", |
| .description = NULL_IF_CONFIG_SMALL("Filter out photosensitive epilepsy seizure-inducing flashes."), |
| .priv_size = sizeof(PhotosensitivityContext), |
| .priv_class = &photosensitivity_class, |
| .uninit = uninit, |
| .query_formats = query_formats, |
| .inputs = inputs, |
| .outputs = outputs, |
| }; |