| 1 | /* |
| 2 | Copyright 2018 Google Inc. All Rights Reserved. |
| 3 | |
| 4 | Licensed under the Apache License, Version 2.0 (the "License"); |
| 5 | you may not use this file except in compliance with the License. |
| 6 | You may obtain a copy of the License at |
| 7 | |
| 8 | http://www.apache.org/licenses/LICENSE-2.0 |
| 9 | |
| 10 | Unless required by applicable law or agreed to in writing, software |
| 11 | distributed under the License is distributed on an "AS-IS" BASIS, |
| 12 | WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. |
| 13 | See the License for the specific language governing permissions and |
| 14 | limitations under the License. |
| 15 | */ |
| 16 | |
| 17 | #include "ambisonics/ambisonic_binaural_decoder.h" |
| 18 | |
| 19 | #include "ambisonics/utils.h" |
| 20 | #include "base/constants_and_types.h" |
| 21 | |
| 22 | |
| 23 | namespace vraudio { |
| 24 | |
| 25 | AmbisonicBinauralDecoder::AmbisonicBinauralDecoder(const AudioBuffer& sh_hrirs, |
| 26 | size_t frames_per_buffer, |
| 27 | FftManager* fft_manager) |
| 28 | : fft_manager_(fft_manager), |
| 29 | freq_input_(kNumMonoChannels, NextPowTwo(input: frames_per_buffer) * 2), |
| 30 | filtered_input_(kNumMonoChannels, frames_per_buffer) { |
| 31 | CHECK(fft_manager_); |
| 32 | CHECK_NE(frames_per_buffer, 0U); |
| 33 | const size_t num_channels = sh_hrirs.num_channels(); |
| 34 | const size_t filter_size = sh_hrirs.num_frames(); |
| 35 | CHECK_NE(num_channels, 0U); |
| 36 | CHECK_NE(filter_size, 0U); |
| 37 | sh_hrir_filters_.reserve(n: num_channels); |
| 38 | for (size_t i = 0; i < num_channels; ++i) { |
| 39 | sh_hrir_filters_.emplace_back( |
| 40 | args: new PartitionedFftFilter(filter_size, frames_per_buffer, fft_manager_)); |
| 41 | sh_hrir_filters_[i]->SetTimeDomainKernel(sh_hrirs[i]); |
| 42 | } |
| 43 | } |
| 44 | |
| 45 | void AmbisonicBinauralDecoder::Process(const AudioBuffer& input, |
| 46 | AudioBuffer* output) { |
| 47 | |
| 48 | DCHECK(output); |
| 49 | DCHECK_EQ(kNumStereoChannels, output->num_channels()); |
| 50 | DCHECK_EQ(input.num_frames(), output->num_frames()); |
| 51 | DCHECK_EQ(input.num_channels(), sh_hrir_filters_.size()); |
| 52 | |
| 53 | output->Clear(); |
| 54 | |
| 55 | AudioBuffer::Channel* freq_input_channel = &freq_input_[0]; |
| 56 | AudioBuffer::Channel* filtered_input_channel = &filtered_input_[0]; |
| 57 | AudioBuffer::Channel* output_channel_0 = &(*output)[0]; |
| 58 | AudioBuffer::Channel* output_channel_1 = &(*output)[1]; |
| 59 | for (size_t channel = 0; channel < input.num_channels(); ++channel) { |
| 60 | const int degree = GetPeriphonicAmbisonicDegreeForChannel(channel); |
| 61 | fft_manager_->FreqFromTimeDomain(time_channel: input[channel], freq_channel: freq_input_channel); |
| 62 | sh_hrir_filters_[channel]->Filter(input: *freq_input_channel); |
| 63 | sh_hrir_filters_[channel]->GetFilteredSignal(output: filtered_input_channel); |
| 64 | if (degree < 0) { |
| 65 | // Degree is negative: spherical harmonic is asymetric. |
| 66 | // So add contributions to the left channel and subtract from the right |
| 67 | // channel. |
| 68 | *output_channel_0 += *filtered_input_channel; |
| 69 | *output_channel_1 -= *filtered_input_channel; |
| 70 | |
| 71 | } else { |
| 72 | // Degree is zero or positive: spherical harmonic is symetric. |
| 73 | // So add contributions to both left and right channels. |
| 74 | *output_channel_0 += *filtered_input_channel; |
| 75 | *output_channel_1 += *filtered_input_channel; |
| 76 | } |
| 77 | } |
| 78 | } |
| 79 | |
| 80 | } // namespace vraudio |
| 81 | |