tomwalters@305: // Copyright 2010, Thomas Walters
tomwalters@305: //
tomwalters@305: // AIM-C: A C++ implementation of the Auditory Image Model
tomwalters@305: // http://www.acousticscale.org/AIMC
tomwalters@305: //
tomwalters@305: // This program is free software: you can redistribute it and/or modify
tomwalters@305: // it under the terms of the GNU General Public License as published by
tomwalters@305: // the Free Software Foundation, either version 3 of the License, or
tomwalters@305: // (at your option) any later version.
tomwalters@305: //
tomwalters@305: // This program is distributed in the hope that it will be useful,
tomwalters@305: // but WITHOUT ANY WARRANTY; without even the implied warranty of
tomwalters@305: // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
tomwalters@305: // GNU General Public License for more details.
tomwalters@305: //
tomwalters@305: // You should have received a copy of the GNU General Public License
tomwalters@305: // along with this program.  If not, see <http://www.gnu.org/licenses/>.
tomwalters@305: 
tomwalters@305: /*!
tomwalters@305:  * \author Thomas Walters <tom@acousticscale.org>
tomwalters@305:  * \date created 2010/02/24
tomwalters@305:  * \version \$Id$
tomwalters@305:  */
tomwalters@305: 
tomwalters@305: // Use the boost RNGs to generate Gaussian noise
tomwalters@305: #include <boost/random.hpp>
tomwalters@305: #include <math.h>
tomwalters@305: 
tomwalters@305: #include "Modules/SNR/ModuleNoise.h"
tomwalters@305: 
tomwalters@305: namespace aimc {
tomwalters@305: ModuleNoise::ModuleNoise(Parameters *params) :
tomwalters@305:     Module(params),
tomwalters@305:     gaussian_variate_(boost::mt19937(),
tomwalters@305:                       boost::normal_distribution<float>(0.0f, 1.0f)) {
tomwalters@305:   module_description_ = "Adds noise to a signal";
tomwalters@305:   module_identifier_ = "noise";
tomwalters@305:   module_type_ = "snr";
tomwalters@305:   module_version_ = "$Id$";
tomwalters@305: 
tomwalters@305:   // Noise level relative to unit-variance Gaussian noise (ie. 0dB will give a
tomwalters@305:   // noise with an RMS level of 1.0)
tomwalters@305:   float snr_db = parameters_->DefaultFloat("noise.level_db", 0.0f);
tomwalters@305:   multiplier_ = pow(10.0f, snr_db / 20.0f);
tomwalters@305: }
tomwalters@305: 
tomwalters@305: ModuleNoise::~ModuleNoise() {
tomwalters@305: }
tomwalters@305: 
tomwalters@305: bool ModuleNoise::InitializeInternal(const SignalBank &input) {
tomwalters@305:   // Copy the parameters of the input signal bank into internal variables, so
tomwalters@305:   // that they can be checked later.
tomwalters@305:   sample_rate_ = input.sample_rate();
tomwalters@305:   buffer_length_ = input.buffer_length();
tomwalters@305:   channel_count_ = input.channel_count();
tomwalters@305: 
tomwalters@305:   output_.Initialize(input);
tomwalters@305:   return true;
tomwalters@305: }
tomwalters@305: 
tomwalters@305: void ModuleNoise::ResetInternal() {
tomwalters@305: 
tomwalters@305: }
tomwalters@305: 
tomwalters@305: void ModuleNoise::Process(const SignalBank &input) {
tomwalters@305:   // Check to see if the module has been initialized. If not, processing
tomwalters@305:   // should not continue.
tomwalters@305:   if (!initialized_) {
tomwalters@305:     LOG_ERROR(_T("Module %s not initialized."), module_identifier_.c_str());
tomwalters@305:     return;
tomwalters@305:   }
tomwalters@305: 
tomwalters@305:   // Check that ths input this time is the same as the input passed to
tomwalters@305:   // Initialize()
tomwalters@305:   if (buffer_length_ != input.buffer_length()
tomwalters@305:       || channel_count_ != input.channel_count()) {
tomwalters@305:     LOG_ERROR(_T("Mismatch between input to Initialize() and input to "
tomwalters@305:                  "Process() in module %s."), module_identifier_.c_str());
tomwalters@305:     return;
tomwalters@305:   }
tomwalters@305: 
tomwalters@305:   for (int c = 0; c < input.channel_count(); ++c) {
tomwalters@305:     for (int i = 0; i < input.buffer_length(); ++i) {
tomwalters@305:       float s = input[c][i];
tomwalters@305:       s += (multiplier_ * gaussian_variate_());
tomwalters@305:       output_.set_sample(c, i, s);
tomwalters@305:     }
tomwalters@305:   }
tomwalters@305:   PushOutput();
tomwalters@305: }
tomwalters@305: }  // namespace aimc
tomwalters@305: