1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
/*
Copyright (C) 2012 Carl Hetherington <cth@carlh.net>
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program 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 General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*/
#include "audio_analysis.h"
#include "analyse_audio_job.h"
#include "compose.hpp"
#include "film.h"
#include "options.h"
#include "decoder_factory.h"
#include "audio_decoder.h"
#include "i18n.h"
using std::string;
using std::max;
using std::min;
using std::cout;
using boost::shared_ptr;
int const AnalyseAudioJob::_num_points = 1024;
AnalyseAudioJob::AnalyseAudioJob (shared_ptr<Film> f)
: Job (f)
, _done (0)
, _samples_per_point (1)
{
}
string
AnalyseAudioJob::name () const
{
return String::compose (_("Analyse audio of %1"), _film->name());
}
void
AnalyseAudioJob::run ()
{
if (!_film->audio_stream () || !_film->length()) {
set_progress (1);
set_state (FINISHED_ERROR);
return;
}
DecodeOptions options;
options.decode_video = false;
Decoders decoders = decoder_factory (_film, options);
assert (decoders.audio);
decoders.audio->set_audio_stream (_film->audio_stream ());
decoders.audio->Audio.connect (bind (&AnalyseAudioJob::audio, this, _1));
int64_t total_audio_frames = video_frames_to_audio_frames (_film->length().get(), _film->audio_stream()->sample_rate(), _film->source_frame_rate());
_samples_per_point = max (int64_t (1), total_audio_frames / _num_points);
_current.resize (_film->audio_stream()->channels ());
_analysis.reset (new AudioAnalysis (_film->audio_stream()->channels()));
while (!decoders.audio->pass()) {
set_progress (float (_done) / total_audio_frames);
}
_analysis->write (_film->audio_analysis_path ());
set_progress (1);
set_state (FINISHED_OK);
}
void
AnalyseAudioJob::audio (shared_ptr<AudioBuffers> b)
{
for (int i = 0; i < b->frames(); ++i) {
for (int j = 0; j < b->channels(); ++j) {
float s = b->data(j)[i];
if (fabsf (s) < 10e-7) {
/* stringstream can't serialise and recover inf or -inf, so prevent such
values by replacing with this (140dB down) */
s = 10e-7;
}
_current[j][AudioPoint::RMS] += pow (s, 2);
_current[j][AudioPoint::PEAK] = max (_current[j][AudioPoint::PEAK], fabsf (s));
if ((_done % _samples_per_point) == 0) {
_current[j][AudioPoint::RMS] = sqrt (_current[j][AudioPoint::RMS] / _samples_per_point);
_analysis->add_point (j, _current[j]);
_current[j] = AudioPoint ();
}
}
++_done;
}
}
|
