1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169
//! Head-Related Transfer Function (HRTF) module. Provides all needed types and methods for HRTF rendering.
//!
//! # Overview
//!
//! HRTF stands for [Head-Related Transfer Function](https://en.wikipedia.org/wiki/Head-related_transfer_function)
//! and can work only with spatial sounds. For each of such sound source after it was processed by HRTF you can
//! definitely tell from which locationsound came from. In other words HRTF improves perception of sound to
//! the level of real life.
//!
//! # HRIR Spheres
//!
//! This library uses Head-Related Impulse Response (HRIR) spheres to create HRTF spheres. HRTF sphere is a set of
//! points in 3D space which are connected into a mesh forming triangulated sphere. Each point contains spectrum
//! for left and right ears which will be used to modify samples from each spatial sound source to create binaural
//! sound. HRIR spheres can be found [here](https://github.com/mrDIMAS/hrir_sphere_builder/tree/master/hrtf_base/IRCAM)
//!
//! # Usage
//!
//! To use HRTF you need to change default renderer to HRTF renderer like so:
//!
//! ```no_run
//! use fyrox_sound::context::{self, SoundContext};
//! use fyrox_sound::renderer::hrtf::{HrtfRenderer};
//! use fyrox_sound::renderer::Renderer;
//! use std::path::Path;
//! use hrtf::HrirSphere;
//!
//! fn use_hrtf(context: &mut SoundContext) {
//! // IRC_1002_C.bin is HRIR sphere in binary format, can be any valid HRIR sphere
//! // from base mentioned above.
//! let hrir_sphere = HrirSphere::from_file("examples/data/IRC_1002_C.bin", context::SAMPLE_RATE).unwrap();
//!
//! context.state().set_renderer(Renderer::HrtfRenderer(HrtfRenderer::new(hrir_sphere)));
//! }
//! ```
//!
//! # Performance
//!
//! HRTF is `heavy`. Usually it 4-5 slower than default renderer, this is essential because HRTF requires some heavy
//! math (fast Fourier transform, convolution, etc.). On Ryzen 1700 it takes 400-450 μs (0.4 - 0.45 ms) per source.
//! In most cases this is ok, engine works in separate thread and it has around 100 ms to prepare new portion of
//! samples for output device.
//!
//! # Known problems
//!
//! This renderer still suffers from small audible clicks in very fast moving sounds, clicks sounds more like
//! "buzzing" - it is due the fact that hrtf is different from frame to frame which gives "bumps" in amplitude
//! of signal because of phase shift each impulse response have. This can be fixed by short cross fade between
//! small amount of samples from previous frame with same amount of frames of current as proposed in
//! [here](http://csoundjournal.com/issue9/newHRTFOpcodes.html)
//!
//! Clicks can be reproduced by using clean sine wave of 440 Hz on some source moving around listener.
use crate::{
context::{self, DistanceModel, SoundContext},
listener::Listener,
renderer::render_source_2d_only,
source::SoundSource,
};
use fyrox_core::{
reflect::prelude::*,
visitor::{Visit, VisitResult, Visitor},
};
use hrtf::HrirSphere;
use std::{
fmt::Debug,
path::{Path, PathBuf},
};
/// See module docs.
#[derive(Clone, Debug, Default, Reflect)]
pub struct HrtfRenderer {
hrir_path: PathBuf,
#[reflect(hidden)]
processor: Option<hrtf::HrtfProcessor>,
}
impl Visit for HrtfRenderer {
fn visit(&mut self, name: &str, visitor: &mut Visitor) -> VisitResult {
let mut region = visitor.enter_region(name)?;
self.hrir_path.visit("ResourcePath", &mut region)?;
drop(region);
if visitor.is_reading() {
self.try_reload();
}
Ok(())
}
}
impl HrtfRenderer {
/// Creates new HRTF renderer using specified HRTF sphere. See module docs for more info.
pub fn new(hrir_sphere: hrtf::HrirSphere) -> Self {
Self {
hrir_path: hrir_sphere.source().to_path_buf(),
processor: Some(hrtf::HrtfProcessor::new(
hrir_sphere,
SoundContext::HRTF_INTERPOLATION_STEPS,
SoundContext::HRTF_BLOCK_LEN,
)),
}
}
/// Tries to load a new HRIR sphere from a given path and re-initialize the renderer to use it.
pub fn set_hrir_sphere_from_path<P: AsRef<Path>>(&mut self, path: P) {
self.hrir_path = path.as_ref().to_path_buf();
self.try_reload();
}
/// Returns currently used HRIR path.
pub fn hrir_sphere_path(&self) -> &Path {
&self.hrir_path
}
fn try_reload(&mut self) {
self.processor = HrirSphere::from_file(&self.hrir_path, context::SAMPLE_RATE)
.ok()
.map(|sphere| {
hrtf::HrtfProcessor::new(
sphere,
SoundContext::HRTF_INTERPOLATION_STEPS,
SoundContext::HRTF_BLOCK_LEN,
)
});
}
pub(crate) fn render_source(
&mut self,
source: &mut SoundSource,
listener: &Listener,
distance_model: DistanceModel,
out_buf: &mut [(f32, f32)],
) {
// Render as 2D first with k = (1.0 - spatial_blend).
render_source_2d_only(source, out_buf);
// Then add HRTF part with k = spatial_blend
let new_distance_gain =
source.spatial_blend() * source.calculate_distance_gain(listener, distance_model);
let new_sampling_vector = source.calculate_sampling_vector(listener);
if let Some(processor) = self.processor.as_mut() {
processor.process_samples(hrtf::HrtfContext {
source: &source.frame_samples,
output: out_buf,
new_sample_vector: hrtf::Vec3::new(
new_sampling_vector.x,
new_sampling_vector.y,
new_sampling_vector.z,
),
prev_sample_vector: hrtf::Vec3::new(
source.prev_sampling_vector.x,
source.prev_sampling_vector.y,
source.prev_sampling_vector.z,
),
prev_left_samples: &mut source.prev_left_samples,
prev_right_samples: &mut source.prev_right_samples,
prev_distance_gain: source.prev_distance_gain.unwrap_or(new_distance_gain),
new_distance_gain,
});
}
source.prev_sampling_vector = new_sampling_vector;
source.prev_distance_gain = Some(new_distance_gain);
}
}