1 |
|
|
/* |
2 |
|
|
* Copyright (c) 2015 Paul B Mahol |
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> |
22 |
|
|
#include <math.h> |
23 |
|
|
|
24 |
|
|
#include "libavcodec/avfft.h" |
25 |
|
|
#include "libavutil/audio_fifo.h" |
26 |
|
|
#include "libavutil/avassert.h" |
27 |
|
|
#include "libavutil/avstring.h" |
28 |
|
|
#include "libavutil/channel_layout.h" |
29 |
|
|
#include "libavutil/intreadwrite.h" |
30 |
|
|
#include "libavutil/opt.h" |
31 |
|
|
#include "libavutil/parseutils.h" |
32 |
|
|
#include "audio.h" |
33 |
|
|
#include "filters.h" |
34 |
|
|
#include "video.h" |
35 |
|
|
#include "avfilter.h" |
36 |
|
|
#include "internal.h" |
37 |
|
|
#include "window_func.h" |
38 |
|
|
|
39 |
|
|
enum DisplayMode { LINE, BAR, DOT, NB_MODES }; |
40 |
|
|
enum ChannelMode { COMBINED, SEPARATE, NB_CMODES }; |
41 |
|
|
enum FrequencyScale { FS_LINEAR, FS_LOG, FS_RLOG, NB_FSCALES }; |
42 |
|
|
enum AmplitudeScale { AS_LINEAR, AS_SQRT, AS_CBRT, AS_LOG, NB_ASCALES }; |
43 |
|
|
|
44 |
|
|
typedef struct ShowFreqsContext { |
45 |
|
|
const AVClass *class; |
46 |
|
|
int w, h; |
47 |
|
|
int mode; |
48 |
|
|
int cmode; |
49 |
|
|
int fft_size; |
50 |
|
|
int fft_bits; |
51 |
|
|
int ascale, fscale; |
52 |
|
|
int avg; |
53 |
|
|
int win_func; |
54 |
|
|
FFTContext *fft; |
55 |
|
|
FFTComplex **fft_data; |
56 |
|
|
float **avg_data; |
57 |
|
|
float *window_func_lut; |
58 |
|
|
float overlap; |
59 |
|
|
float minamp; |
60 |
|
|
int hop_size; |
61 |
|
|
int nb_channels; |
62 |
|
|
int nb_freq; |
63 |
|
|
int win_size; |
64 |
|
|
float scale; |
65 |
|
|
char *colors; |
66 |
|
|
AVAudioFifo *fifo; |
67 |
|
|
int64_t pts; |
68 |
|
|
} ShowFreqsContext; |
69 |
|
|
|
70 |
|
|
#define OFFSET(x) offsetof(ShowFreqsContext, x) |
71 |
|
|
#define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM |
72 |
|
|
|
73 |
|
|
static const AVOption showfreqs_options[] = { |
74 |
|
|
{ "size", "set video size", OFFSET(w), AV_OPT_TYPE_IMAGE_SIZE, {.str = "1024x512"}, 0, 0, FLAGS }, |
75 |
|
|
{ "s", "set video size", OFFSET(w), AV_OPT_TYPE_IMAGE_SIZE, {.str = "1024x512"}, 0, 0, FLAGS }, |
76 |
|
|
{ "mode", "set display mode", OFFSET(mode), AV_OPT_TYPE_INT, {.i64=BAR}, 0, NB_MODES-1, FLAGS, "mode" }, |
77 |
|
|
{ "line", "show lines", 0, AV_OPT_TYPE_CONST, {.i64=LINE}, 0, 0, FLAGS, "mode" }, |
78 |
|
|
{ "bar", "show bars", 0, AV_OPT_TYPE_CONST, {.i64=BAR}, 0, 0, FLAGS, "mode" }, |
79 |
|
|
{ "dot", "show dots", 0, AV_OPT_TYPE_CONST, {.i64=DOT}, 0, 0, FLAGS, "mode" }, |
80 |
|
|
{ "ascale", "set amplitude scale", OFFSET(ascale), AV_OPT_TYPE_INT, {.i64=AS_LOG}, 0, NB_ASCALES-1, FLAGS, "ascale" }, |
81 |
|
|
{ "lin", "linear", 0, AV_OPT_TYPE_CONST, {.i64=AS_LINEAR}, 0, 0, FLAGS, "ascale" }, |
82 |
|
|
{ "sqrt", "square root", 0, AV_OPT_TYPE_CONST, {.i64=AS_SQRT}, 0, 0, FLAGS, "ascale" }, |
83 |
|
|
{ "cbrt", "cubic root", 0, AV_OPT_TYPE_CONST, {.i64=AS_CBRT}, 0, 0, FLAGS, "ascale" }, |
84 |
|
|
{ "log", "logarithmic", 0, AV_OPT_TYPE_CONST, {.i64=AS_LOG}, 0, 0, FLAGS, "ascale" }, |
85 |
|
|
{ "fscale", "set frequency scale", OFFSET(fscale), AV_OPT_TYPE_INT, {.i64=FS_LINEAR}, 0, NB_FSCALES-1, FLAGS, "fscale" }, |
86 |
|
|
{ "lin", "linear", 0, AV_OPT_TYPE_CONST, {.i64=FS_LINEAR}, 0, 0, FLAGS, "fscale" }, |
87 |
|
|
{ "log", "logarithmic", 0, AV_OPT_TYPE_CONST, {.i64=FS_LOG}, 0, 0, FLAGS, "fscale" }, |
88 |
|
|
{ "rlog", "reverse logarithmic", 0, AV_OPT_TYPE_CONST, {.i64=FS_RLOG}, 0, 0, FLAGS, "fscale" }, |
89 |
|
|
{ "win_size", "set window size", OFFSET(fft_size), AV_OPT_TYPE_INT, {.i64=2048}, 16, 65536, FLAGS }, |
90 |
|
|
{ "win_func", "set window function", OFFSET(win_func), AV_OPT_TYPE_INT, {.i64=WFUNC_HANNING}, 0, NB_WFUNC-1, FLAGS, "win_func" }, |
91 |
|
|
{ "rect", "Rectangular", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_RECT}, 0, 0, FLAGS, "win_func" }, |
92 |
|
|
{ "bartlett", "Bartlett", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_BARTLETT}, 0, 0, FLAGS, "win_func" }, |
93 |
|
|
{ "hanning", "Hanning", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_HANNING}, 0, 0, FLAGS, "win_func" }, |
94 |
|
|
{ "hamming", "Hamming", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_HAMMING}, 0, 0, FLAGS, "win_func" }, |
95 |
|
|
{ "blackman", "Blackman", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_BLACKMAN}, 0, 0, FLAGS, "win_func" }, |
96 |
|
|
{ "welch", "Welch", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_WELCH}, 0, 0, FLAGS, "win_func" }, |
97 |
|
|
{ "flattop", "Flat-top", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_FLATTOP}, 0, 0, FLAGS, "win_func" }, |
98 |
|
|
{ "bharris", "Blackman-Harris", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_BHARRIS}, 0, 0, FLAGS, "win_func" }, |
99 |
|
|
{ "bnuttall", "Blackman-Nuttall", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_BNUTTALL}, 0, 0, FLAGS, "win_func" }, |
100 |
|
|
{ "bhann", "Bartlett-Hann", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_BHANN}, 0, 0, FLAGS, "win_func" }, |
101 |
|
|
{ "sine", "Sine", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_SINE}, 0, 0, FLAGS, "win_func" }, |
102 |
|
|
{ "nuttall", "Nuttall", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_NUTTALL}, 0, 0, FLAGS, "win_func" }, |
103 |
|
|
{ "lanczos", "Lanczos", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_LANCZOS}, 0, 0, FLAGS, "win_func" }, |
104 |
|
|
{ "gauss", "Gauss", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_GAUSS}, 0, 0, FLAGS, "win_func" }, |
105 |
|
|
{ "tukey", "Tukey", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_TUKEY}, 0, 0, FLAGS, "win_func" }, |
106 |
|
|
{ "dolph", "Dolph-Chebyshev", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_DOLPH}, 0, 0, FLAGS, "win_func" }, |
107 |
|
|
{ "cauchy", "Cauchy", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_CAUCHY}, 0, 0, FLAGS, "win_func" }, |
108 |
|
|
{ "parzen", "Parzen", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_PARZEN}, 0, 0, FLAGS, "win_func" }, |
109 |
|
|
{ "poisson", "Poisson", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_POISSON}, 0, 0, FLAGS, "win_func" }, |
110 |
|
|
{ "bohman", "Bohman", 0, AV_OPT_TYPE_CONST, {.i64=WFUNC_BOHMAN} , 0, 0, FLAGS, "win_func" }, |
111 |
|
|
{ "overlap", "set window overlap", OFFSET(overlap), AV_OPT_TYPE_FLOAT, {.dbl=1.}, 0., 1., FLAGS }, |
112 |
|
|
{ "averaging", "set time averaging", OFFSET(avg), AV_OPT_TYPE_INT, {.i64=1}, 0, INT32_MAX, FLAGS }, |
113 |
|
|
{ "colors", "set channels colors", OFFSET(colors), AV_OPT_TYPE_STRING, {.str = "red|green|blue|yellow|orange|lime|pink|magenta|brown" }, 0, 0, FLAGS }, |
114 |
|
|
{ "cmode", "set channel mode", OFFSET(cmode), AV_OPT_TYPE_INT, {.i64=COMBINED}, 0, NB_CMODES-1, FLAGS, "cmode" }, |
115 |
|
|
{ "combined", "show all channels in same window", 0, AV_OPT_TYPE_CONST, {.i64=COMBINED}, 0, 0, FLAGS, "cmode" }, |
116 |
|
|
{ "separate", "show each channel in own window", 0, AV_OPT_TYPE_CONST, {.i64=SEPARATE}, 0, 0, FLAGS, "cmode" }, |
117 |
|
|
{ "minamp", "set minimum amplitude", OFFSET(minamp), AV_OPT_TYPE_FLOAT, {.dbl=1e-6}, FLT_MIN, 1e-6, FLAGS }, |
118 |
|
|
{ NULL } |
119 |
|
|
}; |
120 |
|
|
|
121 |
|
|
AVFILTER_DEFINE_CLASS(showfreqs); |
122 |
|
|
|
123 |
|
|
static int query_formats(AVFilterContext *ctx) |
124 |
|
|
{ |
125 |
|
|
AVFilterFormats *formats = NULL; |
126 |
|
|
AVFilterChannelLayouts *layouts = NULL; |
127 |
|
|
AVFilterLink *inlink = ctx->inputs[0]; |
128 |
|
|
AVFilterLink *outlink = ctx->outputs[0]; |
129 |
|
|
static const enum AVSampleFormat sample_fmts[] = { AV_SAMPLE_FMT_FLTP, AV_SAMPLE_FMT_NONE }; |
130 |
|
|
static const enum AVPixelFormat pix_fmts[] = { AV_PIX_FMT_RGBA, AV_PIX_FMT_NONE }; |
131 |
|
|
int ret; |
132 |
|
|
|
133 |
|
|
/* set input audio formats */ |
134 |
|
|
formats = ff_make_format_list(sample_fmts); |
135 |
|
|
if ((ret = ff_formats_ref(formats, &inlink->outcfg.formats)) < 0) |
136 |
|
|
return ret; |
137 |
|
|
|
138 |
|
|
layouts = ff_all_channel_layouts(); |
139 |
|
|
if ((ret = ff_channel_layouts_ref(layouts, &inlink->outcfg.channel_layouts)) < 0) |
140 |
|
|
return ret; |
141 |
|
|
|
142 |
|
|
formats = ff_all_samplerates(); |
143 |
|
|
if ((ret = ff_formats_ref(formats, &inlink->outcfg.samplerates)) < 0) |
144 |
|
|
return ret; |
145 |
|
|
|
146 |
|
|
/* set output video format */ |
147 |
|
|
formats = ff_make_format_list(pix_fmts); |
148 |
|
|
if ((ret = ff_formats_ref(formats, &outlink->incfg.formats)) < 0) |
149 |
|
|
return ret; |
150 |
|
|
|
151 |
|
|
return 0; |
152 |
|
|
} |
153 |
|
|
|
154 |
|
|
static av_cold int init(AVFilterContext *ctx) |
155 |
|
|
{ |
156 |
|
|
ShowFreqsContext *s = ctx->priv; |
157 |
|
|
|
158 |
|
|
s->pts = AV_NOPTS_VALUE; |
159 |
|
|
|
160 |
|
|
return 0; |
161 |
|
|
} |
162 |
|
|
|
163 |
|
|
static int config_output(AVFilterLink *outlink) |
164 |
|
|
{ |
165 |
|
|
AVFilterContext *ctx = outlink->src; |
166 |
|
|
AVFilterLink *inlink = ctx->inputs[0]; |
167 |
|
|
ShowFreqsContext *s = ctx->priv; |
168 |
|
|
float overlap; |
169 |
|
|
int i; |
170 |
|
|
|
171 |
|
|
s->fft_bits = av_log2(s->fft_size); |
172 |
|
|
s->nb_freq = 1 << (s->fft_bits - 1); |
173 |
|
|
s->win_size = s->nb_freq << 1; |
174 |
|
|
av_audio_fifo_free(s->fifo); |
175 |
|
|
av_fft_end(s->fft); |
176 |
|
|
s->fft = av_fft_init(s->fft_bits, 0); |
177 |
|
|
if (!s->fft) { |
178 |
|
|
av_log(ctx, AV_LOG_ERROR, "Unable to create FFT context. " |
179 |
|
|
"The window size might be too high.\n"); |
180 |
|
|
return AVERROR(ENOMEM); |
181 |
|
|
} |
182 |
|
|
|
183 |
|
|
/* FFT buffers: x2 for each (display) channel buffer. |
184 |
|
|
* Note: we use free and malloc instead of a realloc-like function to |
185 |
|
|
* make sure the buffer is aligned in memory for the FFT functions. */ |
186 |
|
|
for (i = 0; i < s->nb_channels; i++) { |
187 |
|
|
av_freep(&s->fft_data[i]); |
188 |
|
|
av_freep(&s->avg_data[i]); |
189 |
|
|
} |
190 |
|
|
av_freep(&s->fft_data); |
191 |
|
|
av_freep(&s->avg_data); |
192 |
|
|
s->nb_channels = inlink->channels; |
193 |
|
|
|
194 |
|
|
s->fft_data = av_calloc(s->nb_channels, sizeof(*s->fft_data)); |
195 |
|
|
if (!s->fft_data) |
196 |
|
|
return AVERROR(ENOMEM); |
197 |
|
|
s->avg_data = av_calloc(s->nb_channels, sizeof(*s->avg_data)); |
198 |
|
|
if (!s->avg_data) |
199 |
|
|
return AVERROR(ENOMEM); |
200 |
|
|
for (i = 0; i < s->nb_channels; i++) { |
201 |
|
|
s->fft_data[i] = av_calloc(s->win_size, sizeof(**s->fft_data)); |
202 |
|
|
s->avg_data[i] = av_calloc(s->nb_freq, sizeof(**s->avg_data)); |
203 |
|
|
if (!s->fft_data[i] || !s->avg_data[i]) |
204 |
|
|
return AVERROR(ENOMEM); |
205 |
|
|
} |
206 |
|
|
|
207 |
|
|
/* pre-calc windowing function */ |
208 |
|
|
s->window_func_lut = av_realloc_f(s->window_func_lut, s->win_size, |
209 |
|
|
sizeof(*s->window_func_lut)); |
210 |
|
|
if (!s->window_func_lut) |
211 |
|
|
return AVERROR(ENOMEM); |
212 |
|
|
generate_window_func(s->window_func_lut, s->win_size, s->win_func, &overlap); |
213 |
|
|
if (s->overlap == 1.) |
214 |
|
|
s->overlap = overlap; |
215 |
|
|
s->hop_size = (1. - s->overlap) * s->win_size; |
216 |
|
|
if (s->hop_size < 1) { |
217 |
|
|
av_log(ctx, AV_LOG_ERROR, "overlap %f too big\n", s->overlap); |
218 |
|
|
return AVERROR(EINVAL); |
219 |
|
|
} |
220 |
|
|
|
221 |
|
|
for (s->scale = 0, i = 0; i < s->win_size; i++) { |
222 |
|
|
s->scale += s->window_func_lut[i] * s->window_func_lut[i]; |
223 |
|
|
} |
224 |
|
|
|
225 |
|
|
outlink->frame_rate = av_make_q(inlink->sample_rate, s->win_size * (1.-s->overlap)); |
226 |
|
|
outlink->sample_aspect_ratio = (AVRational){1,1}; |
227 |
|
|
outlink->w = s->w; |
228 |
|
|
outlink->h = s->h; |
229 |
|
|
|
230 |
|
|
s->fifo = av_audio_fifo_alloc(inlink->format, inlink->channels, s->win_size); |
231 |
|
|
if (!s->fifo) |
232 |
|
|
return AVERROR(ENOMEM); |
233 |
|
|
return 0; |
234 |
|
|
} |
235 |
|
|
|
236 |
|
|
static inline void draw_dot(AVFrame *out, int x, int y, uint8_t fg[4]) |
237 |
|
|
{ |
238 |
|
|
|
239 |
|
|
uint32_t color = AV_RL32(out->data[0] + y * out->linesize[0] + x * 4); |
240 |
|
|
|
241 |
|
|
if ((color & 0xffffff) != 0) |
242 |
|
|
AV_WL32(out->data[0] + y * out->linesize[0] + x * 4, AV_RL32(fg) | color); |
243 |
|
|
else |
244 |
|
|
AV_WL32(out->data[0] + y * out->linesize[0] + x * 4, AV_RL32(fg)); |
245 |
|
|
} |
246 |
|
|
|
247 |
|
|
static int get_sx(ShowFreqsContext *s, int f) |
248 |
|
|
{ |
249 |
|
|
switch (s->fscale) { |
250 |
|
|
case FS_LINEAR: |
251 |
|
|
return (s->w/(float)s->nb_freq)*f; |
252 |
|
|
case FS_LOG: |
253 |
|
|
return s->w-pow(s->w, (s->nb_freq-f-1)/(s->nb_freq-1.)); |
254 |
|
|
case FS_RLOG: |
255 |
|
|
return pow(s->w, f/(s->nb_freq-1.)); |
256 |
|
|
} |
257 |
|
|
|
258 |
|
|
return 0; |
259 |
|
|
} |
260 |
|
|
|
261 |
|
|
static float get_bsize(ShowFreqsContext *s, int f) |
262 |
|
|
{ |
263 |
|
|
switch (s->fscale) { |
264 |
|
|
case FS_LINEAR: |
265 |
|
|
return s->w/(float)s->nb_freq; |
266 |
|
|
case FS_LOG: |
267 |
|
|
return pow(s->w, (s->nb_freq-f-1)/(s->nb_freq-1.))- |
268 |
|
|
pow(s->w, (s->nb_freq-f-2)/(s->nb_freq-1.)); |
269 |
|
|
case FS_RLOG: |
270 |
|
|
return pow(s->w, (f+1)/(s->nb_freq-1.))- |
271 |
|
|
pow(s->w, f /(s->nb_freq-1.)); |
272 |
|
|
} |
273 |
|
|
|
274 |
|
|
return 1.; |
275 |
|
|
} |
276 |
|
|
|
277 |
|
|
static inline void plot_freq(ShowFreqsContext *s, int ch, |
278 |
|
|
double a, int f, uint8_t fg[4], int *prev_y, |
279 |
|
|
AVFrame *out, AVFilterLink *outlink) |
280 |
|
|
{ |
281 |
|
|
const int w = s->w; |
282 |
|
|
const float min = s->minamp; |
283 |
|
|
const float avg = s->avg_data[ch][f]; |
284 |
|
|
const float bsize = get_bsize(s, f); |
285 |
|
|
const int sx = get_sx(s, f); |
286 |
|
|
int end = outlink->h; |
287 |
|
|
int x, y, i; |
288 |
|
|
|
289 |
|
|
switch(s->ascale) { |
290 |
|
|
case AS_SQRT: |
291 |
|
|
a = 1.0 - sqrt(a); |
292 |
|
|
break; |
293 |
|
|
case AS_CBRT: |
294 |
|
|
a = 1.0 - cbrt(a); |
295 |
|
|
break; |
296 |
|
|
case AS_LOG: |
297 |
|
|
a = log(av_clipd(a, min, 1)) / log(min); |
298 |
|
|
break; |
299 |
|
|
case AS_LINEAR: |
300 |
|
|
a = 1.0 - a; |
301 |
|
|
break; |
302 |
|
|
} |
303 |
|
|
|
304 |
|
|
switch (s->cmode) { |
305 |
|
|
case COMBINED: |
306 |
|
|
y = a * outlink->h - 1; |
307 |
|
|
break; |
308 |
|
|
case SEPARATE: |
309 |
|
|
end = (outlink->h / s->nb_channels) * (ch + 1); |
310 |
|
|
y = (outlink->h / s->nb_channels) * ch + a * (outlink->h / s->nb_channels) - 1; |
311 |
|
|
break; |
312 |
|
|
default: |
313 |
|
|
av_assert0(0); |
314 |
|
|
} |
315 |
|
|
if (y < 0) |
316 |
|
|
return; |
317 |
|
|
|
318 |
|
|
switch (s->avg) { |
319 |
|
|
case 0: |
320 |
|
|
y = s->avg_data[ch][f] = !outlink->frame_count_in ? y : FFMIN(avg, y); |
321 |
|
|
break; |
322 |
|
|
case 1: |
323 |
|
|
break; |
324 |
|
|
default: |
325 |
|
|
s->avg_data[ch][f] = avg + y * (y - avg) / (FFMIN(outlink->frame_count_in + 1, s->avg) * y); |
326 |
|
|
y = s->avg_data[ch][f]; |
327 |
|
|
break; |
328 |
|
|
} |
329 |
|
|
|
330 |
|
|
switch(s->mode) { |
331 |
|
|
case LINE: |
332 |
|
|
if (*prev_y == -1) { |
333 |
|
|
*prev_y = y; |
334 |
|
|
} |
335 |
|
|
if (y <= *prev_y) { |
336 |
|
|
for (x = sx + 1; x < sx + bsize && x < w; x++) |
337 |
|
|
draw_dot(out, x, y, fg); |
338 |
|
|
for (i = y; i <= *prev_y; i++) |
339 |
|
|
draw_dot(out, sx, i, fg); |
340 |
|
|
} else { |
341 |
|
|
for (i = *prev_y; i <= y; i++) |
342 |
|
|
draw_dot(out, sx, i, fg); |
343 |
|
|
for (x = sx + 1; x < sx + bsize && x < w; x++) |
344 |
|
|
draw_dot(out, x, i - 1, fg); |
345 |
|
|
} |
346 |
|
|
*prev_y = y; |
347 |
|
|
break; |
348 |
|
|
case BAR: |
349 |
|
|
for (x = sx; x < sx + bsize && x < w; x++) |
350 |
|
|
for (i = y; i < end; i++) |
351 |
|
|
draw_dot(out, x, i, fg); |
352 |
|
|
break; |
353 |
|
|
case DOT: |
354 |
|
|
for (x = sx; x < sx + bsize && x < w; x++) |
355 |
|
|
draw_dot(out, x, y, fg); |
356 |
|
|
break; |
357 |
|
|
} |
358 |
|
|
} |
359 |
|
|
|
360 |
|
|
static int plot_freqs(AVFilterLink *inlink, AVFrame *in) |
361 |
|
|
{ |
362 |
|
|
AVFilterContext *ctx = inlink->dst; |
363 |
|
|
AVFilterLink *outlink = ctx->outputs[0]; |
364 |
|
|
ShowFreqsContext *s = ctx->priv; |
365 |
|
|
const int win_size = s->win_size; |
366 |
|
|
char *colors, *color, *saveptr = NULL; |
367 |
|
|
AVFrame *out; |
368 |
|
|
int ch, n; |
369 |
|
|
|
370 |
|
|
out = ff_get_video_buffer(outlink, outlink->w, outlink->h); |
371 |
|
|
if (!out) |
372 |
|
|
return AVERROR(ENOMEM); |
373 |
|
|
|
374 |
|
|
for (n = 0; n < outlink->h; n++) |
375 |
|
|
memset(out->data[0] + out->linesize[0] * n, 0, outlink->w * 4); |
376 |
|
|
|
377 |
|
|
/* fill FFT input with the number of samples available */ |
378 |
|
|
for (ch = 0; ch < s->nb_channels; ch++) { |
379 |
|
|
const float *p = (float *)in->extended_data[ch]; |
380 |
|
|
|
381 |
|
|
for (n = 0; n < in->nb_samples; n++) { |
382 |
|
|
s->fft_data[ch][n].re = p[n] * s->window_func_lut[n]; |
383 |
|
|
s->fft_data[ch][n].im = 0; |
384 |
|
|
} |
385 |
|
|
for (; n < win_size; n++) { |
386 |
|
|
s->fft_data[ch][n].re = 0; |
387 |
|
|
s->fft_data[ch][n].im = 0; |
388 |
|
|
} |
389 |
|
|
} |
390 |
|
|
|
391 |
|
|
/* run FFT on each samples set */ |
392 |
|
|
for (ch = 0; ch < s->nb_channels; ch++) { |
393 |
|
|
av_fft_permute(s->fft, s->fft_data[ch]); |
394 |
|
|
av_fft_calc(s->fft, s->fft_data[ch]); |
395 |
|
|
} |
396 |
|
|
|
397 |
|
|
#define RE(x, ch) s->fft_data[ch][x].re |
398 |
|
|
#define IM(x, ch) s->fft_data[ch][x].im |
399 |
|
|
#define M(a, b) (sqrt((a) * (a) + (b) * (b))) |
400 |
|
|
|
401 |
|
|
colors = av_strdup(s->colors); |
402 |
|
|
if (!colors) { |
403 |
|
|
av_frame_free(&out); |
404 |
|
|
return AVERROR(ENOMEM); |
405 |
|
|
} |
406 |
|
|
|
407 |
|
|
for (ch = 0; ch < s->nb_channels; ch++) { |
408 |
|
|
uint8_t fg[4] = { 0xff, 0xff, 0xff, 0xff }; |
409 |
|
|
int prev_y = -1, f; |
410 |
|
|
double a; |
411 |
|
|
|
412 |
|
|
color = av_strtok(ch == 0 ? colors : NULL, " |", &saveptr); |
413 |
|
|
if (color) |
414 |
|
|
av_parse_color(fg, color, -1, ctx); |
415 |
|
|
|
416 |
|
|
a = av_clipd(M(RE(0, ch), 0) / s->scale, 0, 1); |
417 |
|
|
plot_freq(s, ch, a, 0, fg, &prev_y, out, outlink); |
418 |
|
|
|
419 |
|
|
for (f = 1; f < s->nb_freq; f++) { |
420 |
|
|
a = av_clipd(M(RE(f, ch), IM(f, ch)) / s->scale, 0, 1); |
421 |
|
|
|
422 |
|
|
plot_freq(s, ch, a, f, fg, &prev_y, out, outlink); |
423 |
|
|
} |
424 |
|
|
} |
425 |
|
|
|
426 |
|
|
av_free(colors); |
427 |
|
|
out->pts = in->pts; |
428 |
|
|
out->sample_aspect_ratio = (AVRational){1,1}; |
429 |
|
|
return ff_filter_frame(outlink, out); |
430 |
|
|
} |
431 |
|
|
|
432 |
|
|
static int filter_frame(AVFilterLink *inlink) |
433 |
|
|
{ |
434 |
|
|
AVFilterContext *ctx = inlink->dst; |
435 |
|
|
ShowFreqsContext *s = ctx->priv; |
436 |
|
|
AVFrame *fin = NULL; |
437 |
|
|
int ret = 0; |
438 |
|
|
|
439 |
|
|
fin = ff_get_audio_buffer(inlink, s->win_size); |
440 |
|
|
if (!fin) { |
441 |
|
|
ret = AVERROR(ENOMEM); |
442 |
|
|
goto fail; |
443 |
|
|
} |
444 |
|
|
|
445 |
|
|
fin->pts = s->pts; |
446 |
|
|
s->pts += s->hop_size; |
447 |
|
|
ret = av_audio_fifo_peek(s->fifo, (void **)fin->extended_data, s->win_size); |
448 |
|
|
if (ret < 0) |
449 |
|
|
goto fail; |
450 |
|
|
|
451 |
|
|
ret = plot_freqs(inlink, fin); |
452 |
|
|
av_frame_free(&fin); |
453 |
|
|
av_audio_fifo_drain(s->fifo, s->hop_size); |
454 |
|
|
|
455 |
|
|
fail: |
456 |
|
|
av_frame_free(&fin); |
457 |
|
|
return ret; |
458 |
|
|
} |
459 |
|
|
|
460 |
|
|
static int activate(AVFilterContext *ctx) |
461 |
|
|
{ |
462 |
|
|
AVFilterLink *inlink = ctx->inputs[0]; |
463 |
|
|
AVFilterLink *outlink = ctx->outputs[0]; |
464 |
|
|
ShowFreqsContext *s = ctx->priv; |
465 |
|
|
AVFrame *in = NULL; |
466 |
|
|
int ret = 0; |
467 |
|
|
|
468 |
|
|
FF_FILTER_FORWARD_STATUS_BACK(outlink, inlink); |
469 |
|
|
|
470 |
|
|
if (av_audio_fifo_size(s->fifo) < s->win_size) |
471 |
|
|
ret = ff_inlink_consume_samples(inlink, s->win_size, s->win_size, &in); |
472 |
|
|
if (ret < 0) |
473 |
|
|
return ret; |
474 |
|
|
if (ret > 0) { |
475 |
|
|
av_audio_fifo_write(s->fifo, (void **)in->extended_data, in->nb_samples); |
476 |
|
|
if (s->pts == AV_NOPTS_VALUE) |
477 |
|
|
s->pts = in->pts; |
478 |
|
|
av_frame_free(&in); |
479 |
|
|
} |
480 |
|
|
|
481 |
|
|
if (av_audio_fifo_size(s->fifo) >= s->win_size) { |
482 |
|
|
ret = filter_frame(inlink); |
483 |
|
|
if (ret <= 0) |
484 |
|
|
return ret; |
485 |
|
|
} |
486 |
|
|
|
487 |
|
|
FF_FILTER_FORWARD_STATUS(inlink, outlink); |
488 |
|
|
FF_FILTER_FORWARD_WANTED(outlink, inlink); |
489 |
|
|
|
490 |
|
|
return FFERROR_NOT_READY; |
491 |
|
|
} |
492 |
|
|
|
493 |
|
|
static av_cold void uninit(AVFilterContext *ctx) |
494 |
|
|
{ |
495 |
|
|
ShowFreqsContext *s = ctx->priv; |
496 |
|
|
int i; |
497 |
|
|
|
498 |
|
|
av_fft_end(s->fft); |
499 |
|
|
for (i = 0; i < s->nb_channels; i++) { |
500 |
|
|
if (s->fft_data) |
501 |
|
|
av_freep(&s->fft_data[i]); |
502 |
|
|
if (s->avg_data) |
503 |
|
|
av_freep(&s->avg_data[i]); |
504 |
|
|
} |
505 |
|
|
av_freep(&s->fft_data); |
506 |
|
|
av_freep(&s->avg_data); |
507 |
|
|
av_freep(&s->window_func_lut); |
508 |
|
|
av_audio_fifo_free(s->fifo); |
509 |
|
|
} |
510 |
|
|
|
511 |
|
|
static const AVFilterPad showfreqs_inputs[] = { |
512 |
|
|
{ |
513 |
|
|
.name = "default", |
514 |
|
|
.type = AVMEDIA_TYPE_AUDIO, |
515 |
|
|
}, |
516 |
|
|
{ NULL } |
517 |
|
|
}; |
518 |
|
|
|
519 |
|
|
static const AVFilterPad showfreqs_outputs[] = { |
520 |
|
|
{ |
521 |
|
|
.name = "default", |
522 |
|
|
.type = AVMEDIA_TYPE_VIDEO, |
523 |
|
|
.config_props = config_output, |
524 |
|
|
}, |
525 |
|
|
{ NULL } |
526 |
|
|
}; |
527 |
|
|
|
528 |
|
|
AVFilter ff_avf_showfreqs = { |
529 |
|
|
.name = "showfreqs", |
530 |
|
|
.description = NULL_IF_CONFIG_SMALL("Convert input audio to a frequencies video output."), |
531 |
|
|
.init = init, |
532 |
|
|
.uninit = uninit, |
533 |
|
|
.query_formats = query_formats, |
534 |
|
|
.priv_size = sizeof(ShowFreqsContext), |
535 |
|
|
.activate = activate, |
536 |
|
|
.inputs = showfreqs_inputs, |
537 |
|
|
.outputs = showfreqs_outputs, |
538 |
|
|
.priv_class = &showfreqs_class, |
539 |
|
|
}; |