From patchwork Sun Dec 6 22:53:58 2020 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 8bit X-Patchwork-Submitter: Paul B Mahol X-Patchwork-Id: 24372 Return-Path: X-Original-To: patchwork@ffaux-bg.ffmpeg.org Delivered-To: patchwork@ffaux-bg.ffmpeg.org Received: from ffbox0-bg.mplayerhq.hu (ffbox0-bg.ffmpeg.org [79.124.17.100]) by ffaux.localdomain (Postfix) with ESMTP id 81F4E44BC9B for ; Mon, 7 Dec 2020 01:01:48 +0200 (EET) Received: from [127.0.1.1] (localhost [127.0.0.1]) by ffbox0-bg.mplayerhq.hu (Postfix) with ESMTP id 4867568A9C9; Mon, 7 Dec 2020 01:01:48 +0200 (EET) X-Original-To: ffmpeg-devel@ffmpeg.org Delivered-To: ffmpeg-devel@ffmpeg.org Received: from mail-ej1-f46.google.com (mail-ej1-f46.google.com [209.85.218.46]) by ffbox0-bg.mplayerhq.hu (Postfix) with ESMTPS id 65CEF68A6D6 for ; Mon, 7 Dec 2020 01:01:41 +0200 (EET) Received: by mail-ej1-f46.google.com with SMTP id g20so16867278ejb.1 for ; Sun, 06 Dec 2020 15:01:41 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=from:to:subject:date:message-id:mime-version :content-transfer-encoding; bh=pDy+Ah3j+X7C9Dfmoyp+oMPPNyogboBZfwQf1ts39tY=; b=UJhANqzcvPs2MXn4Jz5gf0sF77Av2CGsR2EaVk9FfYGZSepkCiRtixsVZ17PYcdo7F qOfHiDfN87kUqu5GSY1HbwXD/dckDeFOKtMHfRp4W++BW2W5BJmLQmKcLQNfLXu/5q+q jA1D7/yMKfLlo6XW1znxqIAOg03X2bhpXwxsCmf1OkQbXLCXqhP+A5P4eNlF6YPKmil6 SBEyWdpRi4bDRQEOLonDYidy2O8xmi0DOJM4zJH+TYltsUBQfsQVfE8pnJwQjDBA0HqN mvVjb9kaAiqSqGtIaFhsYkwIBnrKu198aqst7XSMXmU6McTNHsSjkn57QwP3vtdQtaHg H9CA== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:from:to:subject:date:message-id:mime-version :content-transfer-encoding; bh=pDy+Ah3j+X7C9Dfmoyp+oMPPNyogboBZfwQf1ts39tY=; b=N4yZq5Hq15dS3nwgDr5QdfMlDBNgLXLiARg30YYZVM77Q3x9JQWCAN9YK+nbW07D/k 1GyZuOJbH+GxqRhlQ/MVRvdLQKyyG7E8xJYDeI3B5z0wzCG1KDrANDVD1qtdrb3k+7Pa oRXwvv0aWXpDYCckv6FzmQyCwyqQCZ8JfVJmohfz13H+hI8qCQoRaOU4E3P8PztNzaZJ fTc84CBLJV1JL7QhVTkYl1kgRTFujA4vyyUkcYNCkRZqjjKKajJL8rUVhofWu2PUL2ty XCYzVXAG85tk8XqIBCfHUbBTG1JDZsBX6nojurjy4GvtXOnlCZFZfKoqS+hR2dxXp0G+ 4lPw== X-Gm-Message-State: AOAM531z/7hdzbEeCrUbRYVjNogSfjiEkPnhuWv4KWK8piIaTpmQGmID VI/jkRRlbOjnfs/JOwtXRbM13F8SfUY= X-Google-Smtp-Source: ABdhPJyPVPnUxHlRON8p1IyHZTgZw1WVtEGftLK+Gyemc67tRbIwcOW8Hc5F55rMmUJz4J8EOOFy4w== X-Received: by 2002:a17:906:3953:: with SMTP id g19mr196956eje.429.1607295254774; Sun, 06 Dec 2020 14:54:14 -0800 (PST) Received: from localhost.localdomain ([77.237.101.115]) by smtp.gmail.com with ESMTPSA id t19sm9590362eje.86.2020.12.06.14.54.13 for (version=TLS1_3 cipher=TLS_AES_256_GCM_SHA384 bits=256/256); Sun, 06 Dec 2020 14:54:14 -0800 (PST) From: Paul B Mahol To: ffmpeg-devel@ffmpeg.org Date: Sun, 6 Dec 2020 23:53:58 +0100 Message-Id: <20201206225358.2129-1-onemda@gmail.com> X-Mailer: git-send-email 2.17.1 MIME-Version: 1.0 Subject: [FFmpeg-devel] [PATCH] avfilter: add asuperpass and asuperstop filter X-BeenThere: ffmpeg-devel@ffmpeg.org X-Mailman-Version: 2.1.20 Precedence: list List-Id: FFmpeg development discussions and patches List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Reply-To: FFmpeg development discussions and patches Errors-To: ffmpeg-devel-bounces@ffmpeg.org Sender: "ffmpeg-devel" Signed-off-by: Paul B Mahol --- doc/filters.texi | 50 ++++++++++++++ libavfilter/Makefile | 2 + libavfilter/af_asupercut.c | 133 +++++++++++++++++++++++++++++++++++-- libavfilter/allfilters.c | 2 + 4 files changed, 183 insertions(+), 4 deletions(-) diff --git a/doc/filters.texi b/doc/filters.texi index 3a31f72d79..5519d90327 100644 --- a/doc/filters.texi +++ b/doc/filters.texi @@ -2733,6 +2733,56 @@ Set input gain level. Allowed range is from 0 to 1. Default value is 1. This filter supports the all above options as @ref{commands}. +@section asuperpass +Apply high order Butterworth band-pass filter. + +The filter accepts the following options: + +@table @option +@item centerf +Set center frequency in Hertz. Allowed range is 2 to 192000. +Default value is 1000. + +@item order +Set filter order. Available values are from 4 to 20. +Default value is 4. + +@item qfactor +Set Q-factor. Allowed range is from 0 to 999. Default value is 1. + +@item level +Set input gain level. Allowed range is from 0 to 1. Default value is 1. +@end table + +@subsection Commands + +This filter supports the all above options as @ref{commands}. + +@section asuperstop +Apply high order Butterworth band-stop filter. + +The filter accepts the following options: + +@table @option +@item centerf +Set center frequency in Hertz. Allowed range is 2 to 192000. +Default value is 1000. + +@item order +Set filter order. Available values are from 4 to 20. +Default value is 4. + +@item qfactor +Set Q-factor. Allowed range is from 0 to 999. Default value is 1. + +@item level +Set input gain level. Allowed range is from 0 to 1. Default value is 1. +@end table + +@subsection Commands + +This filter supports the all above options as @ref{commands}. + @section atempo Adjust audio tempo. diff --git a/libavfilter/Makefile b/libavfilter/Makefile index 77cd55688b..4a96e04999 100644 --- a/libavfilter/Makefile +++ b/libavfilter/Makefile @@ -92,6 +92,8 @@ OBJS-$(CONFIG_ASTREAMSELECT_FILTER) += f_streamselect.o framesync.o OBJS-$(CONFIG_ASUBBOOST_FILTER) += af_asubboost.o OBJS-$(CONFIG_ASUBCUT_FILTER) += af_asupercut.o OBJS-$(CONFIG_ASUPERCUT_FILTER) += af_asupercut.o +OBJS-$(CONFIG_ASUPERPASS_FILTER) += af_asupercut.o +OBJS-$(CONFIG_ASUPERSTOP_FILTER) += af_asupercut.o OBJS-$(CONFIG_ATEMPO_FILTER) += af_atempo.o OBJS-$(CONFIG_ATILT_FILTER) += af_atilt.o OBJS-$(CONFIG_ATRIM_FILTER) += trim.o diff --git a/libavfilter/af_asupercut.c b/libavfilter/af_asupercut.c index 6d0a2b79a9..48738ad083 100644 --- a/libavfilter/af_asupercut.c +++ b/libavfilter/af_asupercut.c @@ -1,4 +1,7 @@ /* + * Copyright (c) 2005 Boðaç Topaktaþ + * Copyright (c) 2020 Paul B Mahol + * * This file is part of FFmpeg. * * FFmpeg is free software; you can redistribute it and/or @@ -33,6 +36,7 @@ typedef struct ASuperCutContext { double cutoff; double level; + double qfactor; int order; int filter_count; @@ -93,10 +97,11 @@ static int get_coeffs(AVFilterContext *ctx) if (s->bypass) return 0; - s->filter_count = s->order / 2 + (s->order & 1); - calc_q_factors(s->order, q); - if (!strcmp(ctx->filter->name, "asubcut")) { + s->filter_count = s->order / 2 + (s->order & 1); + + calc_q_factors(s->order, q); + if (s->order & 1) { BiquadCoeffs *coeffs = &s->coeffs[0]; double omega = 2. * tan(M_PI * w0); @@ -119,7 +124,11 @@ static int get_coeffs(AVFilterContext *ctx) coeffs->a1 = -2.0 * (K * K - 1.0) * norm; coeffs->a2 = -(1.0 - K / q[idx] + K * K) * norm; } - } else { + } else if (!strcmp(ctx->filter->name, "asupercut")) { + s->filter_count = s->order / 2 + (s->order & 1); + + calc_q_factors(s->order, q); + if (s->order & 1) { BiquadCoeffs *coeffs = &s->coeffs[0]; double omega = 2. * tan(M_PI * w0); @@ -142,6 +151,74 @@ static int get_coeffs(AVFilterContext *ctx) coeffs->a1 = -2.0 * (K * K - 1.0) * norm; coeffs->a2 = -(1.0 - K / q[idx] + K * K) * norm; } + } else if (!strcmp(ctx->filter->name, "asuperpass")) { + double alpha, beta, gamma, theta; + double theta_0 = 2. * M_PI * (s->cutoff / inlink->sample_rate); + double d_E; + + s->filter_count = s->order / 2; + d_E = (2. * tan(theta_0 / (2. * s->qfactor))) / sin(theta_0); + + for (int b = 0; b < s->filter_count; b += 2) { + double D = 2. * sin(((b + 1) * M_PI) / (2. * s->filter_count)); + double A = (1. + pow((d_E / 2.), 2)) / (D * d_E / 2.); + double d = sqrt((d_E * D) / (A + sqrt(A * A - 1.))); + double B = D * (d_E / 2.) / d; + double W = B + sqrt(B * B - 1.); + + for (int j = 0; j < 2; j++) { + BiquadCoeffs *coeffs = &s->coeffs[b + j]; + + if (j == 1) + theta = 2. * atan(tan(theta_0 / 2.) / W); + else + theta = 2. * atan(W * tan(theta_0 / 2.)); + + beta = 0.5 * ((1. - (d / 2.) * sin(theta)) / (1. + (d / 2.) * sin(theta))); + gamma = (0.5 + beta) * cos(theta); + alpha = 0.5 * (0.5 - beta) * sqrt(1. + pow((W - (1. / W)) / d, 2.)); + + coeffs->a1 = 2. * gamma; + coeffs->a2 = -2. * beta; + coeffs->b0 = 2. * alpha; + coeffs->b1 = 0.; + coeffs->b2 = -2. * alpha; + } + } + } else if (!strcmp(ctx->filter->name, "asuperstop")) { + double alpha, beta, gamma, theta; + double theta_0 = 2. * M_PI * (s->cutoff / inlink->sample_rate); + double d_E; + + s->filter_count = s->order / 2; + d_E = (2. * tan(theta_0 / (2. * s->qfactor))) / sin(theta_0); + + for (int b = 0; b < s->filter_count; b += 2) { + double D = 2. * sin(((b + 1) * M_PI) / (2. * s->filter_count)); + double A = (1. + pow((d_E / 2.), 2)) / (D * d_E / 2.); + double d = sqrt((d_E * D) / (A + sqrt(A * A - 1.))); + double B = D * (d_E / 2.) / d; + double W = B + sqrt(B * B - 1.); + + for (int j = 0; j < 2; j++) { + BiquadCoeffs *coeffs = &s->coeffs[b + j]; + + if (j == 1) + theta = 2. * atan(tan(theta_0 / 2.) / W); + else + theta = 2. * atan(W * tan(theta_0 / 2.)); + + beta = 0.5 * ((1. - (d / 2.) * sin(theta)) / (1. + (d / 2.) * sin(theta))); + gamma = (0.5 + beta) * cos(theta); + alpha = 0.5 * (0.5 + beta) * ((1. - cos(theta)) / (1. - cos(theta_0))); + + coeffs->a1 = 2. * gamma; + coeffs->a2 = -2. * beta; + coeffs->b0 = 2. * alpha; + coeffs->b1 = -4. * alpha * cos(theta_0); + coeffs->b2 = 2. * alpha; + } + } } return 0; @@ -327,3 +404,51 @@ AVFilter ff_af_asubcut = { .flags = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | AVFILTER_FLAG_SLICE_THREADS, }; + +static const AVOption asuperpass_options[] = { + { "centerf","set center frequency", OFFSET(cutoff), AV_OPT_TYPE_DOUBLE, {.dbl=1000}, 2, 192000, FLAGS }, + { "order", "set filter order", OFFSET(order), AV_OPT_TYPE_INT, {.i64=4}, 4, 20, FLAGS }, + { "qfactor","set Q-factor", OFFSET(qfactor),AV_OPT_TYPE_DOUBLE, {.dbl=1.}, 0., 999., FLAGS }, + { "level", "set input level", OFFSET(level), AV_OPT_TYPE_DOUBLE, {.dbl=1.}, 0., 1., FLAGS }, + { NULL } +}; + +AVFILTER_DEFINE_CLASS(asuperpass); + +AVFilter ff_af_asuperpass = { + .name = "asuperpass", + .description = NULL_IF_CONFIG_SMALL("Apply high order Butterworth band-pass filter."), + .query_formats = query_formats, + .priv_size = sizeof(ASuperCutContext), + .priv_class = &asuperpass_class, + .uninit = uninit, + .inputs = inputs, + .outputs = outputs, + .process_command = process_command, + .flags = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | + AVFILTER_FLAG_SLICE_THREADS, +}; + +static const AVOption asuperstop_options[] = { + { "centerf","set center frequency", OFFSET(cutoff), AV_OPT_TYPE_DOUBLE, {.dbl=1000}, 2, 192000, FLAGS }, + { "order", "set filter order", OFFSET(order), AV_OPT_TYPE_INT, {.i64=4}, 4, 20, FLAGS }, + { "qfactor","set Q-factor", OFFSET(qfactor),AV_OPT_TYPE_DOUBLE, {.dbl=1.}, 0., 999., FLAGS }, + { "level", "set input level", OFFSET(level), AV_OPT_TYPE_DOUBLE, {.dbl=1.}, 0., 1., FLAGS }, + { NULL } +}; + +AVFILTER_DEFINE_CLASS(asuperstop); + +AVFilter ff_af_asuperstop = { + .name = "asuperstop", + .description = NULL_IF_CONFIG_SMALL("Apply high order Butterworth band-stop filter."), + .query_formats = query_formats, + .priv_size = sizeof(ASuperCutContext), + .priv_class = &asuperstop_class, + .uninit = uninit, + .inputs = inputs, + .outputs = outputs, + .process_command = process_command, + .flags = AVFILTER_FLAG_SUPPORT_TIMELINE_GENERIC | + AVFILTER_FLAG_SLICE_THREADS, +}; diff --git a/libavfilter/allfilters.c b/libavfilter/allfilters.c index a7f8d9f864..e90fcfb3bb 100644 --- a/libavfilter/allfilters.c +++ b/libavfilter/allfilters.c @@ -85,6 +85,8 @@ extern AVFilter ff_af_astreamselect; extern AVFilter ff_af_asubboost; extern AVFilter ff_af_asubcut; extern AVFilter ff_af_asupercut; +extern AVFilter ff_af_asuperpass; +extern AVFilter ff_af_asuperstop; extern AVFilter ff_af_atempo; extern AVFilter ff_af_atilt; extern AVFilter ff_af_atrim;