Interpolation of room impulse responses in 3d using compressed sensing

In a room, the acoustic transfer between a source and a receiver is described by the so-called “Room Impulse Response”, which depends on both positions and the room characteristics. According to the sampling theorem, directly measuring the full set of acoustic impulse responses within a 3D-space domain would require an unreasonably large number of measurements. Nevertheless, considering that the acoustic wavefield is sparse in some dictionaries, the Compressed Sensing framework allows the recovery of the full wavefield with a reduced set of measurements (microphones), but raises challenging computational and memory issues. In this paper, we exhibit two sparsity assumptions of the wavefield and we derive two practical algorithms for the wavefield estimation. The first one takes advantage of the Modal Theory for the sampling of the Room Impulse Responses in low frequencies (sparsity in frequency), and the second one exploits the Image Source Method for the interpolation of the early reflections (sparsity in time). These two complementary approaches are validated both by numerical and experimental measurements using a 120-microphone 3D array, and results are given as a function of the number of microphones.

[1]  H. Sabine Room Acoustics , 1953, The SAGE Encyclopedia of Human Communication Sciences and Disorders.

[2]  J. Borish Extension of the image model to arbitrary polyhedra , 1984 .

[3]  Nobuhiko Kitawaki,et al.  Common-acoustical-pole and residue model and its application to spatial interpolation and extrapolation of a room transfer function , 1999, IEEE Trans. Speech Audio Process..

[4]  Emmanuel J. Candès,et al.  Decoding by linear programming , 2005, IEEE Transactions on Information Theory.

[5]  Laurent Daudet,et al.  Plate impulse response spatial interpolation with sub-Nyquist sampling , 2011 .

[6]  Manfred R. Schroeder,et al.  Statistical parameters of the frequency response curves of large rooms , 1987 .

[7]  T. Ajdler,et al.  The Plenacoustic Function and Its Sampling , 2006, IEEE Transactions on Signal Processing.

[8]  Michael Elad,et al.  Sparse and Redundant Representations - From Theory to Applications in Signal and Image Processing , 2010 .

[9]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[10]  E.J. Candes Compressive Sampling , 2022 .

[11]  Laurent Daudet,et al.  Compressed sensing for acoustic response reconstruction: Interpolation of the early part , 2011, 2011 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics (WASPAA).

[12]  Laurent Daudet,et al.  Compressively sampling the plenacoustic function , 2011, Optical Engineering + Applications.

[13]  Emmanuel Perrey-Debain,et al.  Plane wave decomposition in the unit disc: convergence estimates and computational aspects , 2006 .

[14]  Laurent Daudet,et al.  Optimal subsampling of multichannel damped sinusoids , 2010, 2010 IEEE Sensor Array and Multichannel Signal Processing Workshop.

[15]  Jean-Marc Jot,et al.  Analysis and synthesis of room reverberation based on a statistical time-frequency model , 1997 .

[16]  Frank Boland,et al.  Acoustic Impulse Response Interpolation for Multichannel Systems Using Dynamic Time Warping , 2009 .

[17]  J. Polack Playing billiards in the concert hall: The mathematical foundations of geometrical room acoustics , 1993 .