X-Git-Url: https://main.carlh.net/gitweb/?p=dcpomatic.git;a=blobdiff_plain;f=src%2Flib%2Faudio_analysis.cc;h=446fcccef1607323d7851dbffe98114cda5fb841;hp=9d708bbfdd07f30e07327781cc328828a4ecce06;hb=8f12e84009d7c2685bb2eeb32665876463d4e6e5;hpb=ab1a666e724911c41dfe08fc96748b38ace839c1 diff --git a/src/lib/audio_analysis.cc b/src/lib/audio_analysis.cc index 9d708bbfd..446fcccef 100644 --- a/src/lib/audio_analysis.cc +++ b/src/lib/audio_analysis.cc @@ -1,95 +1,113 @@ /* - Copyright (C) 2012 Carl Hetherington + Copyright (C) 2012-2018 Carl Hetherington - This program is free software; you can redistribute it and/or modify + This file is part of DCP-o-matic. + + DCP-o-matic 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, + DCP-o-matic 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. + along with DCP-o-matic. If not, see . */ +#include "audio_analysis.h" +#include "cross.h" +#include "util.h" +#include "playlist.h" +#include "audio_content.h" +#include +#include +#include +#include #include #include -#include -#include -#include -#include "audio_analysis.h" +#include +#include +#include using std::ostream; using std::istream; using std::string; -using std::ofstream; -using std::ifstream; using std::vector; using std::cout; using std::max; +using std::pair; +using std::make_pair; using std::list; +using boost::shared_ptr; +using boost::optional; +using boost::dynamic_pointer_cast; +using dcp::raw_convert; +using namespace dcpomatic; -AudioPoint::AudioPoint () -{ - for (int i = 0; i < COUNT; ++i) { - _data[i] = 0; - } -} - -AudioPoint::AudioPoint (istream& s) -{ - for (int i = 0; i < COUNT; ++i) { - s >> _data[i]; - } -} - -void -AudioPoint::write (ostream& s) const -{ - for (int i = 0; i < COUNT; ++i) { - s << _data[i] << "\n"; - } -} - +int const AudioAnalysis::_current_state_version = 3; AudioAnalysis::AudioAnalysis (int channels) { _data.resize (channels); } -AudioAnalysis::AudioAnalysis (string filename) +AudioAnalysis::AudioAnalysis (boost::filesystem::path filename) { - ifstream f (filename.c_str ()); + cxml::Document f ("AudioAnalysis"); + f.read_file (filename); - int channels; - f >> channels; - _data.resize (channels); + if (f.optional_number_child("Version").get_value_or(1) < _current_state_version) { + /* Too old. Throw an exception so that this analysis is re-run. */ + throw OldFormatError ("Audio analysis file is too old"); + } + + BOOST_FOREACH (cxml::NodePtr i, f.node_children ("Channel")) { + vector channel; - for (int i = 0; i < channels; ++i) { - int points; - f >> points; - for (int j = 0; j < points; ++j) { - _data[i].push_back (AudioPoint (f)); + BOOST_FOREACH (cxml::NodePtr j, i->node_children ("Point")) { + channel.push_back (AudioPoint (j)); } + + _data.push_back (channel); + } + + BOOST_FOREACH (cxml::ConstNodePtr i, f.node_children ("SamplePeak")) { + _sample_peak.push_back ( + PeakTime ( + dcp::raw_convert(i->content()), DCPTime(i->number_attribute("Time")) + ) + ); + } + + BOOST_FOREACH (cxml::ConstNodePtr i, f.node_children ("TruePeak")) { + _true_peak.push_back (dcp::raw_convert (i->content ())); } + + _integrated_loudness = f.optional_number_child ("IntegratedLoudness"); + _loudness_range = f.optional_number_child ("LoudnessRange"); + + _analysis_gain = f.optional_number_child ("AnalysisGain"); + _samples_per_point = f.number_child ("SamplesPerPoint"); + _sample_rate = f.number_child ("SampleRate"); + + _leqm = f.optional_number_child("Leqm"); } void AudioAnalysis::add_point (int c, AudioPoint const & p) { - assert (c < channels ()); + DCPOMATIC_ASSERT (c < channels ()); _data[c].push_back (p); } AudioPoint AudioAnalysis::get_point (int c, int p) const { - assert (p < points (c)); + DCPOMATIC_ASSERT (p < points (c)); return _data[c][p]; } @@ -102,24 +120,102 @@ AudioAnalysis::channels () const int AudioAnalysis::points (int c) const { - assert (c < channels ()); + DCPOMATIC_ASSERT (c < channels ()); return _data[c].size (); } void -AudioAnalysis::write (string filename) +AudioAnalysis::write (boost::filesystem::path filename) +{ + shared_ptr doc (new xmlpp::Document); + xmlpp::Element* root = doc->create_root_node ("AudioAnalysis"); + + root->add_child("Version")->add_child_text (raw_convert (_current_state_version)); + + BOOST_FOREACH (vector& i, _data) { + xmlpp::Element* channel = root->add_child ("Channel"); + BOOST_FOREACH (AudioPoint& j, i) { + j.as_xml (channel->add_child ("Point")); + } + } + + for (size_t i = 0; i < _sample_peak.size(); ++i) { + xmlpp::Element* n = root->add_child("SamplePeak"); + n->add_child_text (raw_convert (_sample_peak[i].peak)); + n->set_attribute ("Time", raw_convert (_sample_peak[i].time.get())); + } + + BOOST_FOREACH (float i, _true_peak) { + root->add_child("TruePeak")->add_child_text (raw_convert (i)); + } + + if (_integrated_loudness) { + root->add_child("IntegratedLoudness")->add_child_text (raw_convert (_integrated_loudness.get ())); + } + + if (_loudness_range) { + root->add_child("LoudnessRange")->add_child_text (raw_convert (_loudness_range.get ())); + } + + if (_analysis_gain) { + root->add_child("AnalysisGain")->add_child_text (raw_convert (_analysis_gain.get ())); + } + + root->add_child("SamplesPerPoint")->add_child_text (raw_convert (_samples_per_point)); + root->add_child("SampleRate")->add_child_text (raw_convert (_sample_rate)); + + if (_leqm) { + root->add_child("Leqm")->add_child_text(raw_convert(*_leqm)); + } + + doc->write_to_file_formatted (filename.string ()); +} + +float +AudioAnalysis::gain_correction (shared_ptr playlist) +{ + if (playlist->content().size() == 1 && analysis_gain ()) { + /* In this case we know that the analysis was of a single piece of content and + we know that content's gain when the analysis was run. Hence we can work out + what correction is now needed to make it look `right'. + */ + DCPOMATIC_ASSERT (playlist->content().front()->audio); + return playlist->content().front()->audio->gain() - analysis_gain().get (); + } + + return 0.0f; +} + +/** @return Peak across all channels, and the channel number it is on */ +pair +AudioAnalysis::overall_sample_peak () const +{ + DCPOMATIC_ASSERT (!_sample_peak.empty ()); + + optional pt; + int c = 0; + + for (size_t i = 0; i < _sample_peak.size(); ++i) { + if (!pt || _sample_peak[i].peak > pt->peak) { + pt = _sample_peak[i]; + c = i; + } + } + + return make_pair (pt.get(), c); +} + +optional +AudioAnalysis::overall_true_peak () const { - string tmp = filename + ".tmp"; - - ofstream f (tmp.c_str ()); - f << _data.size() << "\n"; - for (vector >::iterator i = _data.begin(); i != _data.end(); ++i) { - f << i->size () << "\n"; - for (vector::iterator j = i->begin(); j != i->end(); ++j) { - j->write (f); + optional p; + + BOOST_FOREACH (float i, _true_peak) { + if (!p || i > *p) { + p = i; } } - f.close (); - boost::filesystem::rename (tmp, filename); + return p; } +