tomwalters@284: // Copyright 2010, Thomas Walters
tomwalters@284: //
tomwalters@284: // AIM-C: A C++ implementation of the Auditory Image Model
tomwalters@284: // http://www.acousticscale.org/AIMC
tomwalters@284: //
tomwalters@284: // This program is free software: you can redistribute it and/or modify
tomwalters@284: // it under the terms of the GNU General Public License as published by
tomwalters@284: // the Free Software Foundation, either version 3 of the License, or
tomwalters@284: // (at your option) any later version.
tomwalters@284: //
tomwalters@284: // This program is distributed in the hope that it will be useful,
tomwalters@284: // but WITHOUT ANY WARRANTY; without even the implied warranty of
tomwalters@284: // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
tomwalters@284: // GNU General Public License for more details.
tomwalters@284: //
tomwalters@284: // You should have received a copy of the GNU General Public License
tomwalters@284: // along with this program.  If not, see <http://www.gnu.org/licenses/>.
tomwalters@284: 
tomwalters@284: /*!
tomwalters@284:  * \author Thomas Walters <tom@acousticscale.org>
tomwalters@284:  * \date created 2010/02/19
tomwalters@284:  * \version \$Id$
tomwalters@284:  */
tomwalters@284: 
tomwalters@287: #include <cmath>
tomwalters@287: 
tomwalters@284: #include "Modules/SSI/ModuleSSI.h"
tomwalters@284: 
tomwalters@284: namespace aimc {
tomwalters@284: ModuleSSI::ModuleSSI(Parameters *params) : Module(params) {
tomwalters@284:   module_description_ = "Size-shape image (aka the 'sscAI')";
tomwalters@284:   module_identifier_ = "ssi";
tomwalters@284:   module_type_ = "ssi";
tomwalters@284:   module_version_ = "$Id$";
tomwalters@284: 
tomwalters@305:   // Cut off the SSI at the end of the first cycle
tomwalters@305:   do_pitch_cutoff_ = parameters_->DefaultBool("ssi.pitch_cutoff", false);
tomwalters@305: 
tomwalters@305:   // Weight the values in each channel more strongly if the channel was
tomwalters@305:   // truncated due to the pitch cutoff. This ensures that the same amount of
tomwalters@305:   // energy remains in the SSI spectral profile
tomwalters@305:   weight_by_cutoff_ = parameters_->DefaultBool("ssi.weight_by_cutoff", false);
tomwalters@305: 
tomwalters@305:   // Weight the values in each channel more strongly if the channel was
tomwalters@305:   // scaled such that the end goes off the edge of the computed SSI.
tomwalters@305:   // Again, this ensures that the overall energy of the spectral profile
tomwalters@305:   // remains the same.
tomwalters@305:   weight_by_scaling_ = parameters_->DefaultBool("ssi.weight_by_scaling",
tomwalters@305:                                                 false);
tomwalters@305: 
tomwalters@305:   // Time from the zero-lag line of the SAI from which to start searching
tomwalters@305:   // for a maximum in the input SAI's temporal profile.
tomwalters@305:   pitch_search_start_ms_ = parameters_->DefaultFloat(
tomwalters@305:     "ssi.pitch_search_start_ms", 2.0f);
tomwalters@305: 
tomwalters@305:   // Total width in cycles of the whole SSI
tomwalters@305:   ssi_width_cycles_ = parameters_->DefaultFloat("ssi.width_cycles", 10.0f);
tomwalters@305: 
tomwalters@305:   // Set to true to make the cycles axis logarithmic (ie indexing by gamma
tomwalters@305:   // rather than by cycles) 
tomwalters@305:   log_cycles_axis_ = parameters_->DefaultBool("ssi.log_cycles_axis", true);
tomwalters@305: 
tomwalters@305:   // The centre frequency of the channel which will just fill the complete
tomwalters@305:   // width of the SSI buffer
tomwalters@305:   pivot_cf_ = parameters_->DefaultFloat("ssi.pivot_cf", 1000.0f);
tomwalters@284: }
tomwalters@284: 
tomwalters@284: ModuleSSI::~ModuleSSI() {
tomwalters@284: }
tomwalters@284: 
tomwalters@284: bool ModuleSSI::InitializeInternal(const SignalBank &input) {
tomwalters@284:   // Copy the parameters of the input signal bank into internal variables, so
tomwalters@284:   // that they can be checked later.
tomwalters@284:   sample_rate_ = input.sample_rate();
tomwalters@284:   buffer_length_ = input.buffer_length();
tomwalters@284:   channel_count_ = input.channel_count();
tomwalters@284: 
tomwalters@305:   ssi_width_samples_ = sample_rate_ * ssi_width_cycles_ / pivot_cf_;
tomwalters@287:   if (ssi_width_samples_ > buffer_length_) {
tomwalters@287:     ssi_width_samples_ = buffer_length_;
tomwalters@305:     float cycles = ssi_width_samples_ * pivot_cf_ / sample_rate_;
tomwalters@287:     LOG_INFO(_T("Requested SSI width of %f cycles is too long for the "
tomwalters@287:                 "input buffer length of %d samples. The SSI will be "
tomwalters@287:                 "truncated at %d samples wide. This corresponds to a width "
tomwalters@287:                 "of %f cycles."), ssi_width_cycles_, buffer_length_,
tomwalters@287:                 ssi_width_samples_, cycles);
tomwalters@287:     ssi_width_cycles_ = cycles;
tomwalters@287:   }
tomwalters@287:   output_.Initialize(channel_count_, ssi_width_samples_, sample_rate_);
tomwalters@284:   return true;
tomwalters@284: }
tomwalters@284: 
tomwalters@284: void ModuleSSI::ResetInternal() {
tomwalters@284: }
tomwalters@284: 
tomwalters@305: int ModuleSSI::ExtractPitchIndex(const SignalBank &input) const {
tomwalters@305:   // Generate temporal profile of the SAI
tomwalters@305:   vector<float> sai_temporal_profile(buffer_length_, 0.0f);
tomwalters@305:   for (int i = 0; i < buffer_length_; ++i) {
tomwalters@305:     float val = 0.0f;
tomwalters@305:     for (int ch = 0; ch < channel_count_; ++ch) {
tomwalters@305:       val += input.sample(ch, i);
tomwalters@305:     }
tomwalters@305:     sai_temporal_profile[i] = val;
tomwalters@305:   }
tomwalters@305: 
tomwalters@305:   // Find pitch value
tomwalters@305:   int start_sample = floor(pitch_search_start_ms_ * sample_rate_ / 1000.0f);
tomwalters@305:   int max_idx = 0;
tomwalters@305:   float max_val = 0.0f;
tomwalters@305:   for (int i = start_sample; i < buffer_length_; ++i) {
tomwalters@305:     if (sai_temporal_profile[i] > max_val) {
tomwalters@305:       max_idx = i;
tomwalters@305:       max_val = sai_temporal_profile[i];
tomwalters@305:     }
tomwalters@305:   }
tomwalters@305:   return max_idx;
tomwalters@305: }
tomwalters@305: 
tomwalters@284: void ModuleSSI::Process(const SignalBank &input) {
tomwalters@284:   // Check to see if the module has been initialized. If not, processing
tomwalters@284:   // should not continue.
tomwalters@284:   if (!initialized_) {
tomwalters@285:     LOG_ERROR(_T("Module %s not initialized."), module_identifier_.c_str());
tomwalters@284:     return;
tomwalters@284:   }
tomwalters@284: 
tomwalters@284:   // Check that ths input this time is the same as the input passed to
tomwalters@284:   // Initialize()
tomwalters@284:   if (buffer_length_ != input.buffer_length()
tomwalters@284:       || channel_count_ != input.channel_count()) {
tomwalters@284:     LOG_ERROR(_T("Mismatch between input to Initialize() and input to "
tomwalters@285:                  "Process() in module %s."), module_identifier_.c_str());
tomwalters@284:     return;
tomwalters@284:   }
tomwalters@284: 
tomwalters@287:   output_.set_start_time(input.start_time());
tomwalters@284: 
tomwalters@305:   int pitch_index = buffer_length_ - 1;
tomwalters@305:   if (do_pitch_cutoff_) {
tomwalters@305:     pitch_index = ExtractPitchIndex(input);
tomwalters@305:   }
tomwalters@305: 
tomwalters@287:   for (int ch = 0; ch < channel_count_; ++ch) {
tomwalters@305:     float centre_frequency = input.centre_frequency(ch);
tomwalters@287:     // Copy the buffer from input to output, addressing by h-value
tomwalters@287:     for (int i = 0; i < ssi_width_samples_; ++i) {
tomwalters@305:       float h;
tomwalters@305:       float cycle_samples = sample_rate_ / centre_frequency;
tomwalters@305:       if (log_cycles_axis_) {
tomwalters@305:         float gamma_min = -1.0f;
tomwalters@305:         float gamma_max = log2(ssi_width_cycles_);
tomwalters@305:         float gamma = gamma_min + (gamma_max - gamma_min)
tomwalters@305:                                    * static_cast<float>(i)
tomwalters@305:                                    / static_cast<float>(ssi_width_samples_);
tomwalters@305:         h = pow(2.0f, gamma);
tomwalters@305:       } else {
tomwalters@305:         h = static_cast<float>(i) * ssi_width_cycles_
tomwalters@305:             / static_cast<float>(ssi_width_samples_);
tomwalters@305:       }
tomwalters@284: 
tomwalters@287:       // The index into the input array is a floating-point number, which is
tomwalters@287:       // split into a whole part and a fractional part. The whole part and
tomwalters@287:       // fractional part are found, and are used to linearly interpolate
tomwalters@287:       // between input samples to yield an output sample.
tomwalters@287:       double whole_part;
tomwalters@287:       float frac_part = modf(h * cycle_samples, &whole_part);
tomwalters@305:       int sample = floor(whole_part);
tomwalters@305: 
tomwalters@305:       float weight = 1.0f;
tomwalters@305: 
tomwalters@305:       int cutoff_index = buffer_length_ - 1;
tomwalters@305:       if (do_pitch_cutoff_) {
tomwalters@305:         if (pitch_index < cutoff_index) {
tomwalters@305:           if (weight_by_cutoff_) {
tomwalters@305:             weight *= static_cast<float>(buffer_length_)
tomwalters@305:                       / static_cast<float>(pitch_index);
tomwalters@305:           }
tomwalters@305:           cutoff_index = pitch_index;
tomwalters@305:         }
tomwalters@305:       }
tomwalters@305: 
tomwalters@305:       if (weight_by_scaling_) {
tomwalters@305:         if (centre_frequency > pivot_cf_) {
tomwalters@305:           weight *= (centre_frequency / pivot_cf_);
tomwalters@305:         }
tomwalters@305:       }
tomwalters@287: 
tomwalters@287:       float val;
tomwalters@305:       if (sample < cutoff_index) {
tomwalters@287:         float curr_sample = input.sample(ch, sample);
tomwalters@287:         float next_sample = input.sample(ch, sample + 1);
tomwalters@305:         val = weight * (curr_sample
tomwalters@305:                         + frac_part * (next_sample - curr_sample));
tomwalters@287:       } else {
tomwalters@309:         val = 0.0f;
tomwalters@287:       }
tomwalters@287:       output_.set_sample(ch, i, val);
tomwalters@287:     }
tomwalters@287:   }
tomwalters@284:   PushOutput();
tomwalters@284: }
tomwalters@284: }  // namespace aimc
tomwalters@284: