2026 |
M. R. P. Thomas; A. Dhekne Are We There Yet? Indoor Localization and Tracking Online 2026, visited: 02.07.2026. @online{Thomas2026, |
M. R. P. Thomas; N. Tuck 2026. @online{nokey, |
2025 |
S. Koyama; E. de Sena; P. Samarasinghe; M. R. P. Thomas; F Antonacci Past, Present, and Future of Spatial Audio and Room Acoustics Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech, and Signal Processing (ICASSP), 2025. @inproceedings{Koyama2025, |
S. S. Kushwaha; J. Ma; M. R. P. Thomas; Y. Tian; A. Bruni Diff-SAGe: End-to-End Spatial Audio Generation Using Diffusion Models Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech, and Signal Processing (ICASSP), 2025. @inproceedings{Kushwaha2025,Spatial audio is a crucial component in creating immersive experiences. Traditional simulation-based approaches to generate spatial audio rely on expertise, have limited scalability, and assume independence between semantic and spatial information. To address these issues, we explore end-to-end spatial audio generation. We introduce and formulate a new task of generating first-order Ambisonics (FOA) given a sound category and sound source spatial location. We propose Diff-SAGe, an end-to-end, flow-based diffusion-transformer model for this task. Diff-SAGe utilizes a complex spectrogram representation for FOA, preserving the phase information crucial for accurate spatial cues. Additionally, a multi-conditional encoder integrates the input conditions into a unified representation, guiding the generation of FOA waveforms from noise. Through extensive evaluations on two datasets, we demonstrate that our method consistently outperforms traditional simulation-based baselines across both objective and subjective metrics. |
2024 |
A. J Seefeldt; J. B Lando; D. Arteaga; M. R. P. Thomas; G. N Dickins Managing playback of multiple streams of audio over multiple speakers Patent 2024, (US Patent 12,170,875 B2). @patent{Seefeldt2024c,A multi-stream rendering system and method may render and play simultaneously a plurality of audio program streams over a plurality of arbitrarily placed loudspeakers. At least one of the program streams may be a spatial mix. The rendering of said spatial mix may be dynamically modified as a function of the simultaneous rendering of one or more additional program streams. The rendering of one or more additional program streams may be dynamically modified as a function of the simultaneous rendering of the spatial mix. |
M. R. P. Thomas; R. J. Cartwright Acoustic zoning with distributed microphones Patent 2024, (US Patent 12,112,750 B2). @patent{Thomas2024,A method for estimating a user's location in an environment may involve receiving output signals from each microphone of a plurality of microphones in the environment. At least two microphones of the plurality of microphones may be included in separate devices at separate locations in the environment and the output signals may correspond to a current utterance of a user. The method may involve determining multiple current acoustic features from the output signals of each microphone and applying a classifier to the multiple current acoustic features. Applying the classifier may involve applying a model trained on previously-determined acoustic features derived from a plurality of previous utterances made by the user in a plurality of user zones in the environment. The method may involve determining, based at least in part on output from the classifier, an estimate of the user zone in which the user is currently located. |
J. G. Stoddard; D. McGrath; M. R. P. Thomas Robust Bayesian Estimation for Spherical Functions Proceedings Article In: AES Intl. Conf. on Audio for Virtual and Augmented Reality, 2024. @inproceedings{Stoddard2024,Spherical functions are continuous functions defined on the surface of a sphere. They are relevant in modelling a range of spherical phenomena that appear across many disciplines of science and engineering. This paper considers the problem of robust estimation of spherical functions from acoustical datasets that may be noisy and incomplete. Employing a spherical harmonic model representation, and using a Bayesian formulation, prior covariance structures are designed for the spherical harmonic coefficients to reflect the prior knowledge which exists in this estimation setting. Taking an empirical Bayes approach, the prior distributions are optimized based on the available data to maintain a robust and flexible framework, and the final model estimate and confidence bounds can be obtained analytically from a well-defined posterior distribution. Numerical examples are presented for the estimation of head related transfer functions (HRTFs) and a sound source directivity pattern, demonstrating the favourable performance of the proposed method unde |
G. N Dickins; C. G. Hines; D. Gunawan; R. J. Cartwright; A. J. Seefeldt; D. Arteaga; M. R. P. Thomas; J. B Lando Acoustic echo cancellation control for distributed audio devices Patent 2024, (US Patent US 12,003,673 B2). @patent{Dickins2024,An audio processing method may involve receiving output signals from each microphone of a plurality of microphones in an audio environment, the output signals corresponding to a current utterance of a person and determining, based on the output signals, one or more aspects of context information relating to the person, including an estimated current proximity of the person to one or more microphone locations. The method may involve selecting two or more loudspeaker-equipped audio devices based, at least in part, on the one or more aspects of the context information, determining one or more types of audio processing changes to apply to audio data being rendered to loudspeaker feed signals for the audio devices and causing one or more types of audio processing changes to be applied. In some examples, the audio processing changes have the effect of increasing a speech to echo ratio at one or more microphones. |
A. J Seefeldt; J. B Lando; D. Arteaga; G. N Dickins; M. R. P. Thomas Adaptable spatial audio playback Patent 2024, (US Patent 12,003,946 B2). @patent{Seefeldt2024b,A rendering mode may be determined for received audio data, including audio signals and associated spatial data. The audio data may be rendered for reproduction via a set of loudspeakers of an environment according to the rendering mode, to produce rendered audio signals. Rendering the audio data may involve determining relative activation of a set of loudspeakers in an environment. The rendering mode may be variable between a reference spatial mode and one or more distributed spatial modes. The reference spatial mode may have an assumed listening position and orientation. In the distributed spatial mode(s), one or more elements of the audio data may each be rendered in a more spatially distributed manner than in the reference spatial mode and spatial locations of remaining elements of the audio data may be warped such that they span a rendering space of the environment more completely than in the reference spatial mode. |
G. N Dickins; M. R. P. Thomas; A. J. Seefeldt; J. B. Lando; D. Arteaga; C. M. Dyonisio; D. Gunawan; R. J. Cartwright; C. G. Hines Coordination of audio devices Patent 2024, (US Patent US 11,968,268 B2). @patent{Dickins2024a,An audio session management method may involve: determining, by an audio session manager, one or more first media engine capabilities of a first media engine of a first smart audio device, the first media engine being configured for managing one or more audio media streams received by the first smart audio device and for performing first smart audio device signal processing for the one or more audio media streams according to a first media engine sample clock; receiving, by the audio session manager and via a first application communication link, first application control signals from the first application; and controlling the first smart audio device according to the first media engine capabilities, by the audio session manager, via first audio session management control signals transmitted to the first smart audio device via a first smart audio device communication link and without reference to the first media engine sample clock. |
N. R. Tsingos; M. R. P. Thomas; C. Fersch Methods, apparatus and systems for encoding and decoding of directional sound sources Patent 2024. @patent{Tsingos2024,Some disclosed methods involve encoding or decoding directional audio data. Some encoding methods may involve receiving a mono audio signal corresponding to an audio object and a representation of a radiation pattern corresponding to the audio object. The radiation pattern may include sound levels corresponding to plurality of sample times, a plurality of frequency bands and a plurality of directions. The methods may involve encoding the mono audio signal and encoding the source radiation pattern to determine radiation pattern metadata. Encoding the radiation pattern may involve determining a spherical harmonic transform of the representation of the radiation pattern and compressing the spherical harmonic transform to obtain encoded radiation pattern metadata. |
2023 |
J-H. Hanschke; D. Arteaga; G. Cengarle; J. B. Lando; M. R. P. Thomas; A. J. Seefeldt Improved Panning on Non-Equidistant Loudspeakers with Direct Sound Level Compensation Proceedings Article In: Proc. Audio Engineering Soc. Conv., New York, USA, 2023. @inproceedings{Hanschke2023,Loudspeaker rendering techniques that create phantom sound sources often assume an equidistant loudspeaker layout. Typical home setups might not fulfill this condition as loudspeakers deviate from canonical positions, thus requiring a corresponding calibration. The standard approach is to compensate for delays and to match the loudness of each loudspeaker at the listener's location. It was found that a shift of the phantom image occurs when this calibration procedure is applied and one of a pair of loudspeakers is significantly closer to the listener than the other. In this paper, a novel approach to panning on non-equidistant loudspeaker layouts is presented whereby the panning position is governed by the direct sound and the perceived loudness is governed by the full impulse response. Subjective listening tests are presented that validate the approach and quantify the perceived effect of the compensation. In a setup where the standard calibration leads to an average error of 10 degrees, the proposed direct sound compensation largely returns the phantom source to its intended position. |
M. R. P. Thomas; J-H. Hanschke Inverted Cardioid Topology for Multi-Radius Spherical Microphone Arrays Proceedings Article In: Proc. Workshop on App. of Signal Processing to Audio and Acoustics (WASPAA), New Paltz, New York, USA, 2023. @inproceedings{Thomas2023,The design of microphone arrays yielding constant directivity over wide bandwidths produces a tradeoff between narrow element spacing for aliasing performance at high frequencies, and wide spacing for white noise gain performance at low frequencies. In the design of spherical microphone arrays, free-space cardioid microphones allow embedding of spheres at multiple radii, providing constant directivity over wider bandwidths. The physical size of the microphone capsules and cabling limits the practical lower bound on sphere radius and consequentially how high the resulting aliasing frequency can be pushed. This paper proposes an inverted cardioid topology in which the main lobes on the inner sphere capsules face radially inward, instead of outward, allowing denser packing of elements. Practical measurements on a dual-sphere design, targeting third-order soundfield capture, investigate the self-shadowing behavior and show near constant modal beamformer directivity between 400 Hz and 8 kHz. |
2022 |
K. Kjoerling; D. S. McGrath; H. Purnhagen; M. R. P. Thomas Methods and devices for coding soundfield representation signals Patent 2022. @patent{Kjoerling2022,The present document describes a method (40) for encoding a soundfield representation (SR) input signal(101, 301) describing a soundfield at a reference position, wherein the SR input signal (101, 301) comprises a plurality of channels for a plurality of different directivity paterns of the sound field at the reference position. The method (40)comprises extracting (401) one or more audio objects (103, 303) from the SR input signal (101,301). Furthermore, the method (40) comprises determining (402) a residuals signal (102, 302) based on the SR input signal (101 ,301) and based on the one or more audio objects (103, 303). The method(40) also comprises performing joint coding of the one or more audio objects (103, 303) and/or the residual signal (102, 302). In addition, the method (40) comprises generating (403) a bitstream (701) based on data generated in the context of joint coding of the one or more audio objects (103, 303) and/or the residual signal(102,302). |
N. R. Tsingos; M. R. P. Thomas; C. Fersch Methods, apparatus and systems for encoding and decoding of directional sound sources Patent 2022. @patent{Tsingos2022,Some disclosed methods involve encoding or decoding directional audio data. Some encoding methods may involve receiving a mono audio signal corresponding to an audio object and a representation of a radiation pattern corresponding to the audio object. The radiation pattern may include sound levels corresponding to a plurality of sample times, a plurality of frequency bands and a plurality of directions. The methods may involve encoding the mono audio signal and encoding the source radiation pattern to determine radiation pattern metadata. Encoding the radiation pattern may involve determining a spherical harmonic transform of the representation of the radiation pattern and compressing the spherical harmonic transform to obtain encoded radiation pattern metadata. |
2021 |
C. Q. Robinson; M. R. P. Thomas; M. J. Smithers Methods and devices for bass management Patent 17/286,313, 2021. @patent{Robinson2021,Some disclosed methods involve multi-band basa management. Some such examples may involve aplying multiple high-pasand low-pas filter frequencies for the purpose of bass input management. Some disclosed methods treat at least some low-frequency signals asaudio objects that can be panned. Some disclosed methods involve paninng low and high frequencies separately. Folowing high-pass rendering, a power audit may determine a low-frequency deficit factor that is to be reproduced by subwofers or other low-frequency-capable loudspeakers. |
N. Akbar; G. N. Dickins; M. R. P. Thomas; P. Samarasinghe; T. Abhayapala Reducing Modal Error Propagation through Correcting Mismatched Microphone Gains Using Rapid Proceedings Article In: Proc. International Conf. Acoustics, Speech, and Signal Process. (ICASSP) , 2021. @inproceedings{Akbar2021, |
2020 |
N. Akbar, G. N. Dickins, M. R. P. Thomas A Practical Approach for Microphone Array Calibration in Augmented and Virtual Reality Applications Proceedings Article In: Proc. International Conf. on 3D Immersion (IC3D), 2020. @inproceedings{Akbar2020b, |
M. R. P. Thomas; J-H. Hanschke Methods, apparatus and systems for audio sound field capture Patent US 10,721,559, 2020. @patent{Thomas2020, |
N. Akbar; G. Dickins; M. R. P. Thomas; P. Samarasinghe; T. Abhayapala A Novel Method for Obtaining Diffuse Field Measurements for Microphone Calibration Proceedings Article In: Proc. International Conf. Acoustics, Speech, and Signal Process. (ICASSP), Barcelona, Spain, 2020. @inproceedings{Akbar2020,We propose a straightforward and cost-effective method to perform diffuse soundfield measurements for calibrating the magnitude response of a microphone array. Typically, such calibration is performed in a diffuse soundfield created in reverberation chambers, an expensive and time-consuming process. A method is proposed for obtaining diffuse field measurements in untreated environments. First, a closed-form expression for the spatial correlation of a wideband signal in a diffuse field is derived. Next, we describe a practical procedure for obtaining the diffuse field response of a microphone array in the presence of a non-diffuse soundfield by the introduction of random perturbations in the microphone location. Experimental spatial correlation data obtained is compared with the theoretical model, confirming that it is possible to obtain diffuse field measurements in untreated environments with relatively few loudspeakers. A 30 second test signal played from 4-8 loudspeakers is shown to be sufficient in obtaining a diffuse field measurement using the proposed method. An Eigenmike ® is then successfully calibrated at two different geographical locations. |
2019 |
P T Bilinski; J Ahrens; M R P Thomas; I J Tashev; J C Platt; D E Johnston HRTF personalization based on anthropometric features Patent 9,900,722, 2019. @patent{Bilinski2019,The derivation of personalized HRTFs for a human subject based on the anthropometric feature parameters of the human subject involves obtaining multiple anthropometric feature parameters and multiple HRTFs of multiple training subjects. Subsequently, multiple anthropometric feature parameters of a human subject are acquired. A representation of the statistical relationship between the plurality of anthropometric feature parameters of the human subject and a subset of the multiple anthropometric feature parameters belonging to the plurality of training subjects is determined. The representation of the statistical relationship is then applied to the multiple HRTFs of the plurality of training subjects to obtain a set of personalized HRTFs for the human subject. |
P T Bilinski; J Ahrens; M R P Thomas; I J Tashev; J C Platt; D E Johnston HRTF Personalization Based on Anthropometric Features Patent 10,313,818, 2019. @patent{Bilinski2019b,The derivation of personalized HRTFs for a human subject based on the anthropometric feature parameters of the human subject involves obtaining multiple anthropometric feature parameters and multiple HRTFs of multiple training subjects. Subsequently, multiple anthropometric feature parameters of a human subject are acquired. A representation of the statistical relationship between the plurality of anthropometric feature parameters of the human subject and a subset of the multiple anthropometric feature parameters belonging to the plurality of training subjects is determined. The representation of the statistical relationship is then applied to the multiple HRTFs of the plurality of training subjects to obtain a set of personalized HRTFs for the human subject. |
M. R. P. Thomas Practical Concentric Open Sphere Cardioid Microphone Array Design for Higher Order Sound Field Capture Proceedings Article In: Proc. {IEEE} Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brighton, UK, 2019. @inproceedings{Thomas2019,The problem of higher order sound field capture with spherical microphone arrays is considered. While A-format cardioid designs are commonplace for first order capture, interest remains in the increased spatial resolution delivered by higher order arrays. Spherical arrays typically use omnidirectional microphones mounted on a rigid baffle, from which higher order spatial components are estimated by accounting for radial mode strength. This produces a design trade-off between with small arrays for spatial aliasing performance and large arrays for reduced amplification of instrument noise at low frequencies. A practical open sphere design is proposed that contains cardioid microphones mounted at multiple radii to fulfill both criteria. A design example with a two spheres of 16-channel cardioids at 42 mm and 420 mm radius produces white noise gain above unity on third order components down to 200 Hz, a decade lower than a rigid 32-channel 42 mm sphere of omnidirectional microphones. |
P T Bilinski; J Ahrens; M R P Thomas; I J Tashev; J C Platt; D E Johnston HRTF personalization based on anthropometric features Patent 10,284,992, 2019. @patent{Bilinski2019a,The derivation of personalized HRTFs for a human subject based on the anthropometric feature parameters of the human subject involves obtaining multiple anthropometric feature parameters and multiple HRTFs of multiple training subjects. Subsequently, multiple anthropometric feature parameters of a human subject are acquired. A representation of the statistical relationship between the plurality of anthropometric feature parameters of the human subject and a subset of the multiple anthropometric feature parameters belonging to the plurality of training subjects is determined. The representation of the statistical relationship is then applied to the multiple HRTFs of the plurality of training subjects to obtain a set of personalized HRTFs for the human subject. |
H. Gamper; D. Johnston; I. J. Tashev; A. Politis; M. R. P. Thomas Systems and Methods for Non-Parametric Processing of Head Geometry for HRTF Personalization Patent US 10,278,002, 2019. @patent{Gamper2018,Systems and methods for HRTF personalization are provided. More specifically, the systems and methods provide HRTF personalization utilizing non-parametric processing of three-dimensional head scans. Accordingly, the systems and methods for HRTF personalization generate a personalized set of HRTFs for a user without having to extract specific geometric and/or anthropometric features from a three dimensional head scan of a user and/or from the three dimensional head scans of training subjects in a database. |
2018 |
M. R. P. Thomas; C. Q. Robinson Amplitude Panning and the Interior Pan Proceedings Article In: Proc. Tonmeistertagung, Cologne, Germany, 2018. @inproceedings{Thomas2018,The perception of source location using multi-loudspeaker amplitude panning is considered. While there have been many studies on the localization of pairwise panned sound sources, relatively few studies investigate the multi-loudspeaker case. This paper evaluates panning scenarios in which a source is panned on the boundary or within the area bounded by discrete loudspeakers, referred to as boundary and interior pans respectively. Subjective testing of a variety of pan locations reveal the following: (1) as expected, pans to a single loudspeaker yield lowest localization error, (2) pairwise pans tend to be consistently localized closer to the listener than single loudspeaker pans, (3) largest errors occur when the virtual source is panned close to the listener, and (4) interior pans are accurately perceived and, surprisingly, in some cases more accurately than pairwise boundary pans. |
M. R. P. Thomas Practical Open-Concentric Spherical Microphone Array Design Proceedings Article In: J. Acoust. Soc. Am., 2018. @inproceedings{Thomas2018a, |
S Kordon; J-M Batke; A Kreuger; M R P Thomas 10,021,508, 2018. @patent{Kordon2018, |
2017 |
M R P Thomas; C Q Robinson Amplitude Panning and the Interior Pan Proceedings Article In: Proc. Audio Eng. Soc. Convention, New York, USA, 2017. @inproceedings{Thomas2017a, |
M R P Thomas Fast Computation of Cubature Formulae for the Sphere Proceedings Article In: Proc. Joint Workshop on Hands-Free Speech Communication and Microphone Arrays (HSCMA), San Francisco, CA, 2017. @inproceedings{Thomas2017,The near-uniform distribution of nodes on the surface of a sphere has found many uses in numerical integration, physics, chemistry, crystallography, and more recently in the capture, representation and reproduction of spatial audio. A popular solution posed by Fliege and Meyer treats nodes as charged particles that are constrained to lie on the surface of a sphere that undergo mutual repulsion. The potential energy contained within the system forms a cost function that is minimized by numerical optimization using simulated annealing followed by Limited memory-BFGS (L-BFGS). In this work the cost function’s numerical gradient is replaced with an analytical gradient and it is shown that this is sufficient for a L-BFGS alone to achieve similar results. It is also proposed to limit iterations once the solution has converged beyond a threshold for use cases in mechanical designs with much higher error tolerances. A single-core C implementation on a modern machine shows up to 146 times faster convergence times compared with the Fliege implementation. In the limited precision case, all numbers of nodes under 1600 are calculated in under one minute and under 380 in under a second. |
H Gamper; M R P Thomas; I J Tashev; D E Johnston Systems and Methods for Audio Creation and Delivery Patent US 9,609,436, 2017. @patent{Gamper2017, |
2016 |
H Gamper; M R P Thomas; L Corbin; I J Tashev Synthess of Device-Independent Noise Corpora for Realistic ASR Evaluation Proceedings Article In: Proc. Interspeech Conf., San Francisco, CA, 2016. @inproceedings{Gamper2016,In order to effectively evaluate the accuracy of automatic speech recognition (ASR) with a novel capture device, it is important to create a realistic test data corpus that is representative of real-world noise conditions. Typically, this involves either recording the output of a device under test (DUT) in a noisy environment, or synthesizing an environment over loudspeakers in a way that simulates realistic signal-to-noise ratios (SNRs), reverberation times, and spatial noise distributions. Here we propose a method that aims at combining the realism of in-situ recordings with the convenience and repeatability of synthetic corpora. A device-independent spatial recording containing noise and speech is combined with the measured directivity pattern of a DUT to generate a synthetic test corpus for evaluating the performance of an ASR system. This is achieved by a spherical harmonic decomposition of both the sound field and the DUT’s directivity patterns. Experimental results suggest that the proposed method can be a viable alternative to costly and cumbersome device-dependent measurements. The proposed simulation method predicted the SNR of the DUT response to within about 3 dB and the word error rate (WER) to within about 20%, across a range of test SNRs, target source directions, and noise types. |
A Politis; M R P Thomas; H Gamper; I J Tashev Application of 3D Spherical Transforms To Personalization Of Head-Related Transfer Functions Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 2016. @inproceedings{Politis2016,Head-related transfer functions (HRTFs) depend on the shape of the human head and ears, motivating HRTF personalization methods that detect and exploit morphological similarities between subjects in an HRTF database and a new user. Prior work determined similarity from sets of morphological parameters. Here we propose a non-parametric morphological similarity based on a harmonic expansion of head scans. Two 3D spherical transforms are explored for this task, and an appropriate shape similarity metric is defined. A case study focusing on personalisation of interaural time differences (ITDs) is conducted by applying this similarity metric on a database of 3D head scans. |
M R P Thomas; H Gamper; I J Tashev BFGUI: An Interactive Tool for the Synthesis and Analysis of Microphone Array Beamformers Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 2016. @inproceedings{Thomas2016,Microphone arrays are beneficial for distant speech capture because the signals they capture can be exploited with beamforming to suppress noise and reverberation. The theory for the design and analysis of microphone arrays is well established, however the performance of a microphone array beamformer is often subject to conflicting criteria that need to be assessed manually. This paper describes BFGUI, a interactive graphical tool for MATLAB, for simulating microphone arrays and synthesizing beamformers, and whose parameters can be modified and performance metrics monitored in real-time. Primarily aimed at teaching and research, this tool provides the user with an intuitive insight into the effects of microphone types, number and geometry, and the influence of design constraints such as regularization and white noise gain on derived metrics. The resulting directivity pattern, directivity index and front-back ratio are examples of such metrics. Multiple analytic microphone models are supported and external measured microphone directivity patterns can also be loaded. The designs can be then exported in a variety of formats for processing of real-world data. |
2015 |
F Lim; M R P Thomas; P A Naylor; I J Tashev Acoustic Blur Kernel with Sliding Window for Blind Estimation of Reverberation Time Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2015. @inproceedings{Lim2015a,Reverberation time, or T60, is a key parameter used for characterizing acoustic spaces. Blind T60 estimation is useful for many applications including speech intelligibility estimation, acoustic scene analysis and dereverberation. In our previous work, a single-channel blind T60 estimator was proposed employing spectral analysis in the modulation frequency domain. It was shown that the estimation accuracy is crucially affected by the window lengths used for transformation to the modulation domain. In this work, we propose the use of a sliding window length that is dynamically updated depending on the length of the detected decay region. Experimental results demonstrated that in the presence of noise, estimation accuracy was improved over our previous work for T60 up to 700 ms. When compared against two alternative algorithms from the literature, the proposed approach demonstrated higher accuracy for T60 between 500 ms and 1 s. Finally, the proposed approach was shown to be more computationally efficient compared to two of the three alternative algorithms. |
H Gamper; M R P Thomas; I J Tashev Anthropometric parameterisation of a spherical scatterer ITD model with arbitrary ear angles Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, NY, 2015. @inproceedings{Gamper2015a,Accurate modelling of the interaural time difference (ITD) is crucial for rendering localised sound. Parametric models allow personalising ITDs using anthropometrics. However, the mapping between anthropometric features and model parameters is not straightforward. Here, we propose deriving personalised ITD model parameters from a sphere fitted to a 3-D head scan. The proposed ITD personalisation is evaluated on an HRTF database containing 181 subjects, for a simple spherical ITD model as well as a frequency and elevation-dependent spherical scatterer model with arbitrary ear angles. |
M R P Thomas; H Gamper; I J Tashev Dereverberation Sweet Spot Dilation with Combined Channel Equalization and Beamforming Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 2015. @inproceedings{Thomas2015,Beamforming and channel equalizers can be formulated as optimal multichannel filter-and-sum operations with different objective criteria. It has been shown in previous studies that the combination of both concepts under a common framework can yield results that combine both the spatial robustness of beamforming and the dereverberation performance of channel equalization. This paper introduces an additional method for leveraging both approaches that exploits channel estimates in a wanted spatial location and derives robustness from knowledge of the array geometry alone. Experiments with an objective assessment of speech quality as a function of source perturbation reveal that the proposed technique can be viewed as a sweet spot dilator when compared with the MINT channel equalizer. |
F Lim; M R P Thomas; I J Tashev Blur Kernel Estimation Approach to Blind Reverberation Time Estimation Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 2015. @inproceedings{Lim2015b,Reverberation time is an important parameter for characterizing acoustic environments. It is useful in many applications including acoustic scene analysis, robust automatic speech recognition and dereverberation. Given knowledge of the acoustic impulse response, reverberation time can be measured using Schroeder's backward integration method. Since it is not always practical to obtain impulse responses, blind estimation algorithms are sometimes desirable. In this work, the reverberation problem is viewed as an image blurring problem. The blur kernel is estimated through spectral analysis in the modulation domain and the T60 is subsequently estimated from the blur kernel's parameters. It is shown through experimental results that the proposed approach is able to improve robustness to higher T60s especially with increasing levels of additive noise up to an signal-to-noise ratio (SNR) of 10 dB. |
H Gamper; M R P Thomas; I J Tashev Estimation of Multipath Propagation Delays and Interaural Time Differences from 3-D Head Scans Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 2015. @inproceedings{Gamper2015,The estimation of acoustic propagation delays from a sound source to a listener's ear entrances is useful for understanding and visualising the wave propagation along the surface of the head, and necessary for individualised spatial sound rendering. The interaural time difference (ITD) is of particular research interest, as it constitutes one of the main localisation cues exploited by the human auditory system. Here, an approach is proposed that employs ray tracing on a 3-D head scan to estimate and visualise the propagation delays and ITDs from a sound source to a subject's ear entrances. Experimental results indicate that the proposed approach is computationally efficient, and performs equally well or better than optimally tuned parametric ITD models, with a mean absolute ITD estimation error of about 14μs. |
2014 |
M R P Thomas; F Lim; I J Tashev; P A Naylor Optimal Beamforming as a Time Domain Equalization Problem with Application to Room Acoustics Proceedings Article In: Proc. Intl. Workshop Acoust. Echo Noise Control (IWAENC), Juan Les Pins, France, 2014. @inproceedings{Thomas2014b,Signals captured by microphone arrays provide spatial diversity that can be exploited by multichannel processing algorithms to suppress noise and reverberation. Beamforming is a class of approaches that treats the problem with respect to the spatial location of wanted and competing sources, leveraging properties of propagation of waves in free space. A related class of algorithms is channel equalization that exploits knowledge of the acoustic impulse response between a source and microphones with a view to near-perfect dereverberation. Beamforming has been shown to be a very powerful and practical tool in a number of domains, whereas channel equalizers are notoriously sensitive to noise and channel mismatch leading to limited practical applicability. This paper investigates some of the common properties of these algorithms and presents a solution incorporating approaches from both disciplines. |
M R P Thomas; J Ahrens; I J Tashev A Method for Converting Between Cylindrical and Spherical Harmonic Representations of Sound Fields Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 2014. @inproceedings{Thomas2014a,Spherical microphone and circular microphone arrays are useful for sampling sound fields that may be resynthesized with loudspeaker arrays. Spherical microphone arrays are desirable because of their ability to capture three-dimensional sound fields, however it is often more practical to construct loudspeaker arrays in the form of a closed circle located in the horizontal plane. This leads to a spatial undersampling as such a circular sampling can only yield a perfect representation of a height-invariant sound field. This paper investigates the consequences of such spatial undersampling by converting between cylindrical and spherical harmonic decompositions of solutions to the wave equation. We show analytically and via numerical simulations that 1) the result of the spatial undersampling is a purely horizontally propagating sound field, and 2) the ratio of travelling and standing components in the undersampled sound field varies depending on the incidence colatitude. The conversion is also used in a beamforming scenario and shows that the beamformer response becomes increasingly omnidirectional as the source moves away from the horizontal plane. |
P Bilinski; J Ahrens; M R P Thomas; I J Tashev; J C Platt HRTF Magnitude Synthesis via Sparse Representation of Anthropometric Features Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 2014. @inproceedings{Bilinski2014,We propose a method for the synthesis of the magnitudes of Head-related Transfer Functions (HRTFs) using a sparse representation of anthropometric features. Our approach treats the HRTF synthesis problem as finding a sparse representation of the subject's anthropometric features w.r.t. the anthropometric features in the training set. The fundamental assumption is that the magnitudes of a given HRTF set can be described by the same sparse combination as the anthropometric data. Thus, we learn a sparse vector that represents the subject's anthropometric features as a linear superposition of the anthropometric features of a small subset of subjects from the training data. Then, we apply the same sparse vector directly on the HRTF tensor data. For evaluation purpose we use a new dataset, containing both anthropometric features and HRTFs. We compare the proposed sparse representation based approach with ridge regression and with the data of a manikin (which was designed based on average anthropometric data), and we simulate the best and the worst possible classifiers to select one of the HRTFs from the dataset. For instrumental evaluation we use log-spectral distortion. Experiments show that our sparse representation outperforms all other evaluated techniques, and that the synthesized HRTFs are almost as good as the best possible HRTF classifier. |
2013 |
F Lim; M R P Thomas; P A Naylor MINTFormer: A Spatially Aware Channel Equalizer Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2013. @inproceedings{Lim2013,Reverberation is a process that distorts a wanted signal and impairs perceived speech quality. In the context of multichannel dereverberation, channel-based methods and beamforming are two common approaches. Channel-based methods such as the multiple input/output inverse theorem (MINT) can provide perfect dereverberation provided the exact acoustic impulse responses (AIRs) are known. However, they have been shown to be very sensitive to AIR estimation errors for which several modifications have consequently been proposed. Conversely, beamformers are significantly more robust but provide comparatively modest dereverberation. While the two approaches are conventionally considered independent, both can be formulated as a filter-and-sum operation with differing filter design criteria. We propose a unified framework, termed MINT-Forming, that exploits this similarity and introduces a mixing parameter to control the tradeoff between the potential performance of MINT and the robustness of beamforming. Empirical results show that the mixing parameter is a monotonic function of channel estimation error, whereby a MINT solution is preferred when channel estimation error is low. |
J Ahrens; M R P Thomas; I J Tashev Efficient Implementation of the Spectral Division Method for Arbitrary Virtual Sound Fields Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2013. @inproceedings{Ahrens2013b,The Spectral Division Method is an analytic approach for sound field synthesis that determines the loudspeaker driving function in the wavenumber domain. Compact expressions for the driving function in time-frequency domain or in time domain can only be determined for a low number of special cases. Generally, the involved spatial Fourier transforms have to be evaluated numerically. We present a detailed description of the computational procedure and minimize the number of required computations by exploiting the following two aspects: 1) The interval for the spatial sampling of the virtual sound field can be selected for each time-frequency bin, whereby low time-frequency bins can be sampled more coarsely, and 2) the driving function only needs to be evaluated at the locations of the loudspeakers of a given array. The inverse spatial Fourier transform is therefore not required to be evaluated at all initial spatial sampling points but only at those locations that coincide with loudspeakers. |
J Ahrens; M R P Thomas; I J Tashev Gentle Acoustic Crosstalk Cancelation Using the Spectral Division Method and Ambiophonics Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2013. @inproceedings{Ahrens2013a,We propose the concept of gentle acoustic crosstalk cancellation, which aims at reducing the crosstalk between a loudspeaker and the listener’s contralateral ear instead of eliminating it completely as aggressive methods intend to do. The expected benefit is higher robustness and a tendency to collapse less unpleasantly. The proposed method employs a linear loudspeaker array and exhibits two stages: 1) Use the Spectral Division Method to illuminate the ipsilateral ear using constructive interference of the loudspeaker signals. This approach provides only little channel separation between the listener’s ears at frequencies below approximately 2000 Hz. 2) There we additionally use destructive interference by Recursive Ambiophonics Crosstalk Elimination (RACE). RACE was chosen because of its tendency to collapse gently. In a sample scenario with realistic parameters, the proposed method achieves around 20 dB of channel separation between 700 Hz and 9000 Hz, which appears to be sufficient to achieve full perceived lateralization when only one ear is illuminated. |
2012 |
D P Jarrett; E A P Habets; M R P Thomas; P A Naylor Rigid sphere room impulse response simulation: algorithm and applications Journal Article In: J. Acoust. Soc. Am., vol. 132, no. 3, pp. 1462–1472, 2012. @article{Jarrett2012,Simulated room impulse responses have been proven to be both useful and indispensable for comprehensive testing of acoustic signal processing algorithms while controlling parameters such as the reverberation time, room dimensions, and source–array distance. In this work, a method is proposed for simulating the room impulse responses between a sound source and the microphones positioned on a spherical array. The method takes into account specular reflections of the source by employing the well-known image method, and scattering from the rigid sphere by employing spherical harmonic decomposition. Pseudocode for the proposed method is provided, taking into account various optimizations to reduce the computational complexity. The magnitude and phase errors that result from the finite order spherical harmonic decomposition are analyzed and general guidelines for the order selection are provided. Three examples are presented: an analysis of a diffuse reverberant sound field, a study of binaural cues in the presence of reverberation, and an illustration of the algorithm’s use as a mouth simulator. |
M R P Thomas; J Ahrens; I Tashev Optimal 3Đ Beamforming Using Measured Microphone Directivity Patterns Proceedings Article In: Proc. Intl. Workshop Acoust. Signal Enhancement, Aachen, Germany, 2012. @inproceedings{Thomas2012b,The design of time-invariant beamformers is often posed as an optimization problem using practical design constraints. In many scenarios it is sufficient to assume that the microphones have an omnidirectional directivity pattern, a flat frequency response in the range of interest, and a 2D environment in which wavefronts propagate as a function of azimuth angle only. In this paper we consider a generalized solution for those cases in which one or more of these assumptions do not hold, yielding a beamformer that is optimized on measured directivity patterns as a function of azimuth, elevation and frequency. A comparative study is made with the 4-element cardioid microphone array employed in Microsoft Kinect for Windows, whose beamformer weights are calculated with directivity patterns using (a) 2D cardioid models, (b) 3D cardioid models and (c) 3D measurements. Results on a recorded noisy speech corpus show similar PESQ and speech recognition accuracy comparing (a) and (b), but a 50% relative improvement in word error rate using measured directivity patterns. |
M R P Thomas; N D Gaubitch; E A P Habets; P A Naylor An Insight into Common Filtering in Noisy SIMO Blind System Identification Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, Kyoto, Japan, 2012. @inproceedings{Thomas2012c,The effect of additive sensor noise on single-input-multiple-output (SIMO) blind system identification (BSI) algorithms based upon cross-relation (CR) error is investigated. Previous studies have shown that additive noise in the observed signal results in systems comprising the true estimated channels convolved with an erroneous `common filter', and additionally that identification and removal of this filter significantly improves estimation error. However, the source of the common filter remained an open question. This paper explains the common filter through a first-order perturbation analysis of the CR matrix, showing that it be estimated from the perturbation and the eigenvectors of the noiseless CR matrix. The analysis given in this paper provides a new insight into the effect of noise on SIMO BSI algorithms and forms the first step towards an overall noise robust solution. |
J Ahrens; M R P Thomas; I J Tashev HRTF Modeling Using a Non-Regularized Least-Squares Fit of Spherical Harmonics Coefficients Proceedings Article In: Proc. Asia-Pacific Signal and Information Process. Assoc. Annu. Summit, Aachen, Germany, 2012. @inproceedings{Ahrens2012,Head-related transfer functions (HRTFs) represent the acoustic transfer function from a sound source at a given location to the ear drums of a human. They are typically measured from discrete source positions at a constant distance. Spherical harmonics decompositions have been shown to provide a flexible representation of HRTFs. Practical constraints often prevent the retrieval of measurement data from certain directions, a circumstance that complicates the decomposition of the measured data into spherical harmonics. A least-squares fit of coefficients is a potential approach to determining the coefficients of incomplete data. However, a straightforward non-regularized fit tends to give unrealistic estimates for the region were no measurement data is available. Recently, a regularized least-squares fit was proposed, which yields well-behaved results for the unknown region at the expense of reducing the accuracy of the data representation in the known region. In this paper, we propose using a lower-order non-regularized least-squares fit to achieve a well-behaved estimation of the unknown data. This data then allows for a high-order non-regularized least-squares fit over the entire sphere. We compare the properties of all three approaches applied to modeling the magnitudes of the HRTFs measured from a manikin. The proposed approach reduces the normalized mean-square error by approximately 7 dB in the known region and 11 dB in the unknown region compared to the regularized fit. |
M R P Thomas; J Ahrens; I Tashev Beamformer Design Using Measured Microphone Directivity Patterns: Robustness to Modelling Error Proceedings Article In: Proc. Asia-Pacific Signal and Information Process. Assoc. Annu. Summit, San Fransisco, USA, 2012. @inproceedings{Thomas2012c,Abstract: The design process for time-invariant acoustic beamformers often assumes that the microphones have an omni-directional directivity pattern, a flat frequency response in the range of interest, and a 2D environment in which wavefronts propagate as a function of azimuth angle only. In this paper we investigate those cases in which one or more of these assumptions do not hold, considering a Minimum Variance Distortionless Response (MVDR)-based solution that is optimized using measured directivity patterns as a function of azimuth, elevation and frequency. Robustness to modelling error is controlled by a regularization parameter that produces a suboptimal but more robust solution. A comparative study is made with the 4-element cardioid microphone array employed in Microsoft Kinect for Windows, whose beamformer weights are calculated with directivity patterns using (a) 2D cardioid models, (b) 3D cardioid models and (c) 3D measurements. Speech recognition and PESQ results are used as evaluation criteria with a noisy speech corpus, revealing empirically optimal regularization parameters for each case and up to a 70% relative improvement in word error rate comparing (a) and (c). |
Publications
2026 |
Are We There Yet? Indoor Localization and Tracking Online 2026, visited: 02.07.2026. |
2026. |
2025 |
Past, Present, and Future of Spatial Audio and Room Acoustics Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech, and Signal Processing (ICASSP), 2025. |
Diff-SAGe: End-to-End Spatial Audio Generation Using Diffusion Models Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech, and Signal Processing (ICASSP), 2025. |
2024 |
Managing playback of multiple streams of audio over multiple speakers Patent 2024, (US Patent 12,170,875 B2). |
Acoustic zoning with distributed microphones Patent 2024, (US Patent 12,112,750 B2). |
Robust Bayesian Estimation for Spherical Functions Proceedings Article In: AES Intl. Conf. on Audio for Virtual and Augmented Reality, 2024. |
Acoustic echo cancellation control for distributed audio devices Patent 2024, (US Patent US 12,003,673 B2). |
Adaptable spatial audio playback Patent 2024, (US Patent 12,003,946 B2). |
Coordination of audio devices Patent 2024, (US Patent US 11,968,268 B2). |
Methods, apparatus and systems for encoding and decoding of directional sound sources Patent 2024. |
2023 |
Improved Panning on Non-Equidistant Loudspeakers with Direct Sound Level Compensation Proceedings Article In: Proc. Audio Engineering Soc. Conv., New York, USA, 2023. |
Inverted Cardioid Topology for Multi-Radius Spherical Microphone Arrays Proceedings Article In: Proc. Workshop on App. of Signal Processing to Audio and Acoustics (WASPAA), New Paltz, New York, USA, 2023. |
2022 |
Methods and devices for coding soundfield representation signals Patent 2022. |
Methods, apparatus and systems for encoding and decoding of directional sound sources Patent 2022. |
2021 |
Methods and devices for bass management Patent 17/286,313, 2021. |
Reducing Modal Error Propagation through Correcting Mismatched Microphone Gains Using Rapid Proceedings Article In: Proc. International Conf. Acoustics, Speech, and Signal Process. (ICASSP) , 2021. |
2020 |
A Practical Approach for Microphone Array Calibration in Augmented and Virtual Reality Applications Proceedings Article In: Proc. International Conf. on 3D Immersion (IC3D), 2020. |
Methods, apparatus and systems for audio sound field capture Patent US 10,721,559, 2020. |
A Novel Method for Obtaining Diffuse Field Measurements for Microphone Calibration Proceedings Article In: Proc. International Conf. Acoustics, Speech, and Signal Process. (ICASSP), Barcelona, Spain, 2020. |
2019 |
HRTF personalization based on anthropometric features Patent 9,900,722, 2019. |
HRTF Personalization Based on Anthropometric Features Patent 10,313,818, 2019. |
Practical Concentric Open Sphere Cardioid Microphone Array Design for Higher Order Sound Field Capture Proceedings Article In: Proc. {IEEE} Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brighton, UK, 2019. |
HRTF personalization based on anthropometric features Patent 10,284,992, 2019. |
Systems and Methods for Non-Parametric Processing of Head Geometry for HRTF Personalization Patent US 10,278,002, 2019. |
2018 |
Amplitude Panning and the Interior Pan Proceedings Article In: Proc. Tonmeistertagung, Cologne, Germany, 2018. |
Practical Open-Concentric Spherical Microphone Array Design Proceedings Article In: J. Acoust. Soc. Am., 2018. |
10,021,508, 2018. |
2017 |
Amplitude Panning and the Interior Pan Proceedings Article In: Proc. Audio Eng. Soc. Convention, New York, USA, 2017. |
Fast Computation of Cubature Formulae for the Sphere Proceedings Article In: Proc. Joint Workshop on Hands-Free Speech Communication and Microphone Arrays (HSCMA), San Francisco, CA, 2017. |
Systems and Methods for Audio Creation and Delivery Patent US 9,609,436, 2017. |
2016 |
Synthess of Device-Independent Noise Corpora for Realistic ASR Evaluation Proceedings Article In: Proc. Interspeech Conf., San Francisco, CA, 2016. |
Application of 3D Spherical Transforms To Personalization Of Head-Related Transfer Functions Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 2016. |
BFGUI: An Interactive Tool for the Synthesis and Analysis of Microphone Array Beamformers Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Shanghai, China, 2016. |
2015 |
Acoustic Blur Kernel with Sliding Window for Blind Estimation of Reverberation Time Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2015. |
Anthropometric parameterisation of a spherical scatterer ITD model with arbitrary ear angles Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, NY, 2015. |
Dereverberation Sweet Spot Dilation with Combined Channel Equalization and Beamforming Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 2015. |
Blur Kernel Estimation Approach to Blind Reverberation Time Estimation Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 2015. |
Estimation of Multipath Propagation Delays and Interaural Time Differences from 3-D Head Scans Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 2015. |
2014 |
Optimal Beamforming as a Time Domain Equalization Problem with Application to Room Acoustics Proceedings Article In: Proc. Intl. Workshop Acoust. Echo Noise Control (IWAENC), Juan Les Pins, France, 2014. |
A Method for Converting Between Cylindrical and Spherical Harmonic Representations of Sound Fields Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 2014. |
HRTF Magnitude Synthesis via Sparse Representation of Anthropometric Features Proceedings Article In: Proc. IEEE Intl. Conf. on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 2014. |
2013 |
MINTFormer: A Spatially Aware Channel Equalizer Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2013. |
Efficient Implementation of the Spectral Division Method for Arbitrary Virtual Sound Fields Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2013. |
Gentle Acoustic Crosstalk Cancelation Using the Spectral Division Method and Ambiophonics Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, New Paltz, New York, 2013. |
2012 |
Rigid sphere room impulse response simulation: algorithm and applications Journal Article In: J. Acoust. Soc. Am., vol. 132, no. 3, pp. 1462–1472, 2012. |
Optimal 3Đ Beamforming Using Measured Microphone Directivity Patterns Proceedings Article In: Proc. Intl. Workshop Acoust. Signal Enhancement, Aachen, Germany, 2012. |
An Insight into Common Filtering in Noisy SIMO Blind System Identification Proceedings Article In: Proc. IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, Kyoto, Japan, 2012. |
HRTF Modeling Using a Non-Regularized Least-Squares Fit of Spherical Harmonics Coefficients Proceedings Article In: Proc. Asia-Pacific Signal and Information Process. Assoc. Annu. Summit, Aachen, Germany, 2012. |
Beamformer Design Using Measured Microphone Directivity Patterns: Robustness to Modelling Error Proceedings Article In: Proc. Asia-Pacific Signal and Information Process. Assoc. Annu. Summit, San Fransisco, USA, 2012. |