ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 1 | /* |
| 2 | * Copyright (c) 2014 The WebRTC project authors. All Rights Reserved. |
| 3 | * |
| 4 | * Use of this source code is governed by a BSD-style license |
| 5 | * that can be found in the LICENSE file in the root of the source |
| 6 | * tree. An additional intellectual property rights grant can be found |
| 7 | * in the file PATENTS. All contributing project authors may |
| 8 | * be found in the AUTHORS file in the root of the source tree. |
| 9 | */ |
| 10 | |
| 11 | #ifndef WEBRTC_MODULES_AUDIO_PROCESSING_INTELLIGIBILITY_INTELLIGIBILITY_ENHANCER_H_ |
| 12 | #define WEBRTC_MODULES_AUDIO_PROCESSING_INTELLIGIBILITY_INTELLIGIBILITY_ENHANCER_H_ |
| 13 | |
| 14 | #include <complex> |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 15 | #include <memory> |
ekm | 35b72fb | 2015-07-10 14:11:52 -0700 | [diff] [blame] | 16 | #include <vector> |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 17 | |
| 18 | #include "webrtc/common_audio/lapped_transform.h" |
ekmeyerson | 60d9b33 | 2015-08-14 10:35:55 -0700 | [diff] [blame] | 19 | #include "webrtc/common_audio/channel_buffer.h" |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 20 | #include "webrtc/modules/audio_processing/intelligibility/intelligibility_utils.h" |
Alejandro Luebs | 18fcbcf | 2016-02-22 15:57:38 -0800 | [diff] [blame] | 21 | #include "webrtc/modules/audio_processing/vad/voice_activity_detector.h" |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 22 | |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 23 | namespace webrtc { |
| 24 | |
| 25 | // Speech intelligibility enhancement module. Reads render and capture |
| 26 | // audio streams and modifies the render stream with a set of gains per |
| 27 | // frequency bin to enhance speech against the noise background. |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 28 | // Details of the model and algorithm can be found in the original paper: |
| 29 | // http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=6882788 |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 30 | class IntelligibilityEnhancer { |
| 31 | public: |
Alejandro Luebs | 18fcbcf | 2016-02-22 15:57:38 -0800 | [diff] [blame] | 32 | IntelligibilityEnhancer(int sample_rate_hz, size_t num_render_channels); |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 33 | |
aluebs | c466bad | 2016-02-10 12:03:00 -0800 | [diff] [blame] | 34 | // Sets the capture noise magnitude spectrum estimate. |
| 35 | void SetCaptureNoiseEstimate(std::vector<float> noise); |
ekm | b7553df | 2015-06-16 18:57:32 -0700 | [diff] [blame] | 36 | |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 37 | // Reads chunk of speech in time domain and updates with modified signal. |
ekmeyerson | 60d9b33 | 2015-08-14 10:35:55 -0700 | [diff] [blame] | 38 | void ProcessRenderAudio(float* const* audio, |
| 39 | int sample_rate_hz, |
Peter Kasting | 6955870 | 2016-01-12 16:26:35 -0800 | [diff] [blame] | 40 | size_t num_channels); |
ekmeyerson | 60d9b33 | 2015-08-14 10:35:55 -0700 | [diff] [blame] | 41 | bool active() const; |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 42 | |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 43 | private: |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 44 | // Provides access point to the frequency domain. |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 45 | class TransformCallback : public LappedTransform::Callback { |
| 46 | public: |
aluebs | c466bad | 2016-02-10 12:03:00 -0800 | [diff] [blame] | 47 | TransformCallback(IntelligibilityEnhancer* parent); |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 48 | |
| 49 | // All in frequency domain, receives input |in_block|, applies |
| 50 | // intelligibility enhancement, and writes result to |out_block|. |
pkasting | b297c5a | 2015-07-22 15:17:22 -0700 | [diff] [blame] | 51 | void ProcessAudioBlock(const std::complex<float>* const* in_block, |
Peter Kasting | 6955870 | 2016-01-12 16:26:35 -0800 | [diff] [blame] | 52 | size_t in_channels, |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 53 | size_t frames, |
Peter Kasting | 6955870 | 2016-01-12 16:26:35 -0800 | [diff] [blame] | 54 | size_t out_channels, |
pkasting | b297c5a | 2015-07-22 15:17:22 -0700 | [diff] [blame] | 55 | std::complex<float>* const* out_block) override; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 56 | |
| 57 | private: |
| 58 | IntelligibilityEnhancer* parent_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 59 | }; |
| 60 | friend class TransformCallback; |
ekm | 35b72fb | 2015-07-10 14:11:52 -0700 | [diff] [blame] | 61 | FRIEND_TEST_ALL_PREFIXES(IntelligibilityEnhancerTest, TestErbCreation); |
| 62 | FRIEND_TEST_ALL_PREFIXES(IntelligibilityEnhancerTest, TestSolveForGains); |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 63 | |
Alejandro Luebs | 3234819 | 2016-02-17 20:04:19 -0800 | [diff] [blame] | 64 | // Updates power computation and analysis with |in_block_|, |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 65 | // and writes modified speech to |out_block|. |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 66 | void ProcessClearBlock(const std::complex<float>* in_block, |
| 67 | std::complex<float>* out_block); |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 68 | |
ekm | 35b72fb | 2015-07-10 14:11:52 -0700 | [diff] [blame] | 69 | // Bisection search for optimal |lambda|. |
| 70 | void SolveForLambda(float power_target, float power_bot, float power_top); |
| 71 | |
| 72 | // Transforms freq gains to ERB gains. |
| 73 | void UpdateErbGains(); |
| 74 | |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 75 | // Returns number of ERB filters. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 76 | static size_t GetBankSize(int sample_rate, size_t erb_resolution); |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 77 | |
| 78 | // Initializes ERB filterbank. |
aluebs | c466bad | 2016-02-10 12:03:00 -0800 | [diff] [blame] | 79 | std::vector<std::vector<float>> CreateErbBank(size_t num_freqs); |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 80 | |
| 81 | // Analytically solves quadratic for optimal gains given |lambda|. |
| 82 | // Negative gains are set to 0. Stores the results in |sols|. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 83 | void SolveForGainsGivenLambda(float lambda, size_t start_freq, float* sols); |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 84 | |
Alejandro Luebs | 18fcbcf | 2016-02-22 15:57:38 -0800 | [diff] [blame] | 85 | // Returns true if the audio is speech. |
| 86 | bool IsSpeech(const float* audio); |
| 87 | |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 88 | const size_t freqs_; // Num frequencies in frequency domain. |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 89 | const size_t chunk_length_; // Chunk size in samples. |
| 90 | const size_t bank_size_; // Num ERB filters. |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 91 | const int sample_rate_hz_; |
Peter Kasting | 6955870 | 2016-01-12 16:26:35 -0800 | [diff] [blame] | 92 | const size_t num_render_channels_; |
ekmeyerson | 60d9b33 | 2015-08-14 10:35:55 -0700 | [diff] [blame] | 93 | |
Alejandro Luebs | 18fcbcf | 2016-02-22 15:57:38 -0800 | [diff] [blame] | 94 | intelligibility::PowerEstimator<std::complex<float>> clear_power_estimator_; |
| 95 | std::unique_ptr<intelligibility::PowerEstimator<float>> |
| 96 | noise_power_estimator_; |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 97 | std::unique_ptr<float[]> filtered_clear_pow_; |
| 98 | std::unique_ptr<float[]> filtered_noise_pow_; |
| 99 | std::unique_ptr<float[]> center_freqs_; |
aluebs | c466bad | 2016-02-10 12:03:00 -0800 | [diff] [blame] | 100 | std::vector<std::vector<float>> capture_filter_bank_; |
| 101 | std::vector<std::vector<float>> render_filter_bank_; |
Peter Kasting | dce40cf | 2015-08-24 14:52:23 -0700 | [diff] [blame] | 102 | size_t start_freq_; |
Alejandro Luebs | 18fcbcf | 2016-02-22 15:57:38 -0800 | [diff] [blame] | 103 | |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 104 | std::unique_ptr<float[]> gains_eq_; // Pre-filter modified gains. |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 105 | intelligibility::GainApplier gain_applier_; |
| 106 | |
ekmeyerson | 60d9b33 | 2015-08-14 10:35:55 -0700 | [diff] [blame] | 107 | // Destination buffers used to reassemble blocked chunks before overwriting |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 108 | // the original input array with modifications. |
ekmeyerson | 60d9b33 | 2015-08-14 10:35:55 -0700 | [diff] [blame] | 109 | ChannelBuffer<float> temp_render_out_buffer_; |
ekm | db4fecf | 2015-06-22 17:49:08 -0700 | [diff] [blame] | 110 | |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 111 | TransformCallback render_callback_; |
kwiberg | 88788ad | 2016-02-19 07:04:49 -0800 | [diff] [blame] | 112 | std::unique_ptr<LappedTransform> render_mangler_; |
Alejandro Luebs | 18fcbcf | 2016-02-22 15:57:38 -0800 | [diff] [blame] | 113 | |
| 114 | VoiceActivityDetector vad_; |
| 115 | std::vector<int16_t> audio_s16_; |
| 116 | size_t chunks_since_voice_; |
| 117 | bool is_speech_; |
ekm | 030249d | 2015-06-15 13:02:24 -0700 | [diff] [blame] | 118 | }; |
| 119 | |
| 120 | } // namespace webrtc |
| 121 | |
| 122 | #endif // WEBRTC_MODULES_AUDIO_PROCESSING_INTELLIGIBILITY_INTELLIGIBILITY_ENHANCER_H_ |