| 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 | #if defined(_WIN32) |
| 18 | #define _SCL_SECURE_NO_WARNINGS |
| 19 | #endif |
| 20 | |
| 21 | #include "ambisonics/foa_rotator.h" |
| 22 | |
| 23 | #include <algorithm> |
| 24 | |
| 25 | #include "base/constants_and_types.h" |
| 26 | #include "base/misc_math.h" |
| 27 | |
| 28 | namespace vraudio { |
| 29 | |
| 30 | bool FoaRotator::Process(const WorldRotation& target_rotation, |
| 31 | const AudioBuffer& input, AudioBuffer* output) { |
| 32 | |
| 33 | DCHECK(output); |
| 34 | DCHECK_EQ(input.num_channels(), kNumFirstOrderAmbisonicChannels); |
| 35 | DCHECK_EQ(input.num_channels(), output->num_channels()); |
| 36 | DCHECK_EQ(input.num_frames(), output->num_frames()); |
| 37 | |
| 38 | static const WorldRotation kIdentityRotation; |
| 39 | |
| 40 | if (current_rotation_.AngularDifferenceRad(other: kIdentityRotation) < |
| 41 | kRotationQuantizationRad && |
| 42 | target_rotation.AngularDifferenceRad(other: kIdentityRotation) < |
| 43 | kRotationQuantizationRad) { |
| 44 | return false; |
| 45 | } |
| 46 | |
| 47 | if (current_rotation_.AngularDifferenceRad(other: target_rotation) < |
| 48 | kRotationQuantizationRad) { |
| 49 | // Rotate the whole input buffer frame by frame. |
| 50 | Rotate(target_rotation: current_rotation_, start_location: 0, duration: input.num_frames(), input, output); |
| 51 | return true; |
| 52 | } |
| 53 | |
| 54 | // In order to perform a smooth rotation, we divide the buffer into |
| 55 | // chunks of size |kSlerpFrameInterval|. |
| 56 | // |
| 57 | |
| 58 | const size_t kSlerpFrameInterval = 32; |
| 59 | |
| 60 | WorldRotation slerped_rotation; |
| 61 | // Rotate the input buffer at every slerp update interval. Truncate the |
| 62 | // final chunk if the input buffer is not an integer multiple of the |
| 63 | // chunk size. |
| 64 | for (size_t i = 0; i < input.num_frames(); i += kSlerpFrameInterval) { |
| 65 | const size_t duration = |
| 66 | std::min(a: input.num_frames() - i, b: kSlerpFrameInterval); |
| 67 | const float interpolation_factor = static_cast<float>(i + duration) / |
| 68 | static_cast<float>(input.num_frames()); |
| 69 | slerped_rotation = |
| 70 | current_rotation_.slerp(t: interpolation_factor, other: target_rotation); |
| 71 | // Rotate the input buffer frame by frame within the current chunk. |
| 72 | Rotate(target_rotation: slerped_rotation, start_location: i, duration, input, output); |
| 73 | } |
| 74 | |
| 75 | current_rotation_ = target_rotation; |
| 76 | |
| 77 | return true; |
| 78 | } |
| 79 | |
| 80 | void FoaRotator::Rotate(const WorldRotation& target_rotation, |
| 81 | size_t start_location, size_t duration, |
| 82 | const AudioBuffer& input, AudioBuffer* output) { |
| 83 | |
| 84 | const AudioBuffer::Channel& input_channel_audio_space_w = input[0]; |
| 85 | const AudioBuffer::Channel& input_channel_audio_space_y = input[1]; |
| 86 | const AudioBuffer::Channel& input_channel_audio_space_z = input[2]; |
| 87 | const AudioBuffer::Channel& input_channel_audio_space_x = input[3]; |
| 88 | AudioBuffer::Channel* output_channel_audio_space_w = &(*output)[0]; |
| 89 | AudioBuffer::Channel* output_channel_audio_space_y = &(*output)[1]; |
| 90 | AudioBuffer::Channel* output_channel_audio_space_z = &(*output)[2]; |
| 91 | AudioBuffer::Channel* output_channel_audio_space_x = &(*output)[3]; |
| 92 | |
| 93 | for (size_t frame = start_location; frame < start_location + duration; |
| 94 | ++frame) { |
| 95 | // Convert the current audio frame into world space position. |
| 96 | temp_audio_position_(0) = input_channel_audio_space_x[frame]; |
| 97 | temp_audio_position_(1) = input_channel_audio_space_y[frame]; |
| 98 | temp_audio_position_(2) = input_channel_audio_space_z[frame]; |
| 99 | ConvertWorldFromAudioPosition(audio_position: temp_audio_position_, world_position: &temp_world_position_); |
| 100 | // Apply rotation to |world_position| and return to audio space. |
| 101 | temp_rotated_world_position_ = target_rotation * temp_world_position_; |
| 102 | |
| 103 | ConvertAudioFromWorldPosition(world_position: temp_rotated_world_position_, |
| 104 | audio_position: &temp_rotated_audio_position_); |
| 105 | (*output_channel_audio_space_x)[frame] = |
| 106 | temp_rotated_audio_position_(0); // X |
| 107 | (*output_channel_audio_space_y)[frame] = |
| 108 | temp_rotated_audio_position_(1); // Y |
| 109 | (*output_channel_audio_space_z)[frame] = |
| 110 | temp_rotated_audio_position_(2); // Z |
| 111 | } |
| 112 | // Copy W channel. |
| 113 | std::copy_n(first: &input_channel_audio_space_w[start_location], n: duration, |
| 114 | result: &(*output_channel_audio_space_w)[start_location]); |
| 115 | } |
| 116 | |
| 117 | } // namespace vraudio |
| 118 | |