44 #define MAX_CHANNELS 63
47 #define PRE_B0 1.53512485958697
48 #define PRE_B1 -2.69169618940638
49 #define PRE_B2 1.19839281085285
50 #define PRE_A1 -1.69065929318241
51 #define PRE_A2 0.73248077421585
57 #define RLB_A1 -1.99004745483398
58 #define RLB_A2 0.99007225036621
61 #define ABS_UP_THRES 10
62 #define HIST_GRAIN 100
63 #define HIST_SIZE ((ABS_UP_THRES - ABS_THRES) * HIST_GRAIN + 1)
117 #define I400_BINS (48000 * 4 / 10)
118 #define I3000_BINS (48000 * 3)
128 #define OFFSET(x) offsetof(EBUR128Context, x)
129 #define A AV_OPT_FLAG_AUDIO_PARAM
130 #define V AV_OPT_FLAG_VIDEO_PARAM
131 #define F AV_OPT_FLAG_FILTERING_PARAM
154 const int below0 = y > ebur128->
y_zero_lu;
155 const int reached = y >= v;
157 const int colorid = 4*line + 2*reached + below0;
163 v += 2 * ebur128->
meter;
193 for (i = 0; buf[i]; i++) {
197 for (char_y = 0; char_y < font_height; char_y++) {
198 for (mask = 0x80;
mask; mask >>= 1) {
199 if (font[buf[i] * font_height + char_y] & mask)
202 memcpy(p,
"\x00\x00\x00", 3);
215 for (i = 0; i <
len; i++) {
216 memcpy(p,
"\x00\xff\x00", 3);
230 if (ebur128->
w < 640 || ebur128->
h < 480) {
232 "minimum size is 640x480\n", ebur128->
w, ebur128->
h);
235 outlink->
w = ebur128->
w;
236 outlink->
h = ebur128->
h;
242 ebur128->
text.
y = 40;
243 ebur128->
text.
w = 3 * 8;
265 outlink->
w, outlink->
h);
276 memset(outpicref->
data[0], 0, ebur128->
h * outpicref->
linesize[0]);
280 for (i = ebur128->
meter; i >= -ebur128->
meter * 2; i--) {
282 x =
PAD + (i < 10 && i > -10) * 8;
286 "%c%d", i < 0 ? '-' : i > 0 ?
'+' :
' ',
FFABS(i));
293 for (y = 0; y < ebur128->
graph.
h; y++) {
296 for (x = 0; x < ebur128->
graph.
w; x++)
297 memcpy(p + x*3, c, 3);
302 #define DRAW_RECT(r) do { \
303 drawline(outpicref, r.x, r.y - 1, r.w, 3); \
304 drawline(outpicref, r.x, r.y + r.h, r.w, 3); \
305 drawline(outpicref, r.x - 1, r.y, r.h, outpicref->linesize[0]); \
306 drawline(outpicref, r.x + r.w, r.y, r.h, outpicref->linesize[0]); \
321 #define BACK_MASK (AV_CH_BACK_LEFT |AV_CH_BACK_CENTER |AV_CH_BACK_RIGHT| \
322 AV_CH_TOP_BACK_LEFT|AV_CH_TOP_BACK_CENTER|AV_CH_TOP_BACK_RIGHT)
349 #define ENERGY(loudness) (pow(10, ((loudness) + 0.691) / 10.))
350 #define LOUDNESS(energy) (-0.691 + 10 * log10(energy))
370 ebur128->
class = &ebur128_class;
412 #define HIST_POS(power) (int)(((power) - ABS_THRES) * HIST_GRAIN)
420 double relative_threshold;
431 if (!relative_threshold)
432 relative_threshold = 1e-12;
436 return gate_hist_pos;
446 const double *
samples = (
double *)insamples->
data[0];
453 #define MOVE_TO_NEXT_CACHED_ENTRY(time) do { \
454 ebur128->i##time.cache_pos++; \
455 if (ebur128->i##time.cache_pos == I##time##_BINS) { \
456 ebur128->i##time.filled = 1; \
457 ebur128->i##time.cache_pos = 0; \
471 #define FILTER(Y, X, name) do { \
472 double *dst = ebur128->Y + ch*3; \
473 double *src = ebur128->X + ch*3; \
476 dst[0] = src[0]*name##_B0 + src[1]*name##_B1 + src[2]*name##_B2 \
477 - dst[1]*name##_A1 - dst[2]*name##_A2; \
480 ebur128->
x[ch * 3] = *samples++;
484 ebur128->
x[ch * 3 + 2] = ebur128->
x[ch * 3 + 1];
485 ebur128->
x[ch * 3 + 1] = ebur128->
x[ch * 3 ];
488 bin = ebur128->
z[ch * 3] * ebur128->
z[ch * 3];
496 ebur128->
i400.
cache [ch][bin_id_400 ] = bin;
504 double loudness_400, loudness_3000;
505 double power_400 = 1e-12, power_3000 = 1e-12;
507 const int64_t pts = insamples->
pts +
513 #define COMPUTE_LOUDNESS(m, time) do { \
514 if (ebur128->i##time.filled) { \
516 for (ch = 0; ch < nb_channels; ch++) \
517 power_##time += ebur128->ch_weighting[ch] * ebur128->i##time.sum[ch]; \
518 power_##time /= I##time##_BINS; \
520 loudness_##time = LOUDNESS(power_##time); \
527 #define I_GATE_THRES -10 // initially defined to -8 LU in the first EBU standard
530 double integrated_sum = 0;
531 int nb_integrated = 0;
537 for (i = gate_hist_pos; i <
HIST_SIZE; i++) {
539 nb_integrated += nb_v;
547 #define LRA_GATE_THRES -20
548 #define LRA_LOWER_PRC 10
549 #define LRA_HIGHER_PRC 95
558 for (i = gate_hist_pos; i <
HIST_SIZE; i++)
566 for (i = gate_hist_pos; i <
HIST_SIZE; i++) {
577 for (i = HIST_SIZE - 1; i >= 0; i--) {
590 #define LOG_FMT "M:%6.1f S:%6.1f I:%6.1f LUFS LRA:%6.1f LU"
597 const int y_loudness_lu_graph =
lu_to_y(ebur128, loudness_3000 + 23);
598 const int y_loudness_lu_gauge =
lu_to_y(ebur128, loudness_400 + 23);
601 p = pic->data[0] + ebur128->
graph.
y*pic->linesize[0] + ebur128->
graph.
x*3;
602 for (y = 0; y < ebur128->
graph.
h; y++) {
605 memmove(p, p + 3, (ebur128->
graph.
w - 1) * 3);
606 memcpy(p + (ebur128->
graph.
w - 1) * 3, c, 3);
607 p += pic->linesize[0];
611 p = pic->data[0] + ebur128->
gauge.
y*pic->linesize[0] + ebur128->
gauge.
x*3;
612 for (y = 0; y < ebur128->
gauge.
h; y++) {
615 for (x = 0; x < ebur128->
gauge.
w; x++)
616 memcpy(p + x*3, c, 3);
617 p += pic->linesize[0];
623 loudness_400, loudness_3000,
635 loudness_400, loudness_3000,
652 static const int input_srate[] = {48000, -1};
705 " Integrated loudness:\n"
707 " Threshold: %5.1f LUFS\n\n"
710 " Threshold: %5.1f LUFS\n"
711 " LRA low: %5.1f LUFS\n"
712 " LRA high: %5.1f LUFS\n",
749 .priv_class = &ebur128_class,