Novel Sampling Scheme on the Sphere for Head-Related Transfer Function Measurements

This paper presents a novel sampling scheme on the sphere for obtaining head-related transfer function (HRTF) measurements and accurately computing the spherical harmonic transform (SHT). The scheme requires an optimal number of samples, given by the degrees of freedom in the spectral domain, for the accurate representation of the HRTF that is band-limited in the spherical harmonic domain. The proposed scheme allows for the samples to be easily taken over the sphere due to its iso-latitude structure and non-dense sampling near the poles. In addition, the scheme can be used when samples are not taken from the south polar cap region of the sphere as the HRTF measurements are not reliable in south polar cap region due to reflections from the ground. Furthermore, the scheme has a hierarchical structure, which enables the HRTF to be analyzed at different audible frequencies using the same sampling configuration. In comparison to the proposed scheme, none of the other sampling schemes on the sphere simultaneously possess all these properties. We conduct several numerical experiments to determine the accuracy of the SHT associated with the proposed sampling scheme. We show that the SHT attains accuracy on the order of numerical precision (10-14) when samples are taken over the whole sphere, both in the optimal sample placement and hierarchical configurations, and achieves an acceptable level of accuracy (10-5) when samples are not taken over the south polar cap region of the sphere for the band-limits of interest. Simulations are used to show the accurate reconstruction of the HRTF over the whole sphere, including unmeasured locations.

[1]  R. Kennedy,et al.  Hilbert Space Methods in Signal Processing , 2013 .

[2]  Yves Wiaux,et al.  A Novel Sampling Theorem on the Sphere , 2011, IEEE Transactions on Signal Processing.

[3]  Ramani Duraiswami,et al.  Kernel regression for Head-Related Transfer Function interpolation and spectral extrema extraction , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[4]  Russell L. Martin,et al.  Interpolation of Head-Related Transfer Functions , 2007 .

[5]  Ramani Duraiswami,et al.  Sound field decomposition using spherical microphone arrays , 2008, 2008 IEEE International Conference on Acoustics, Speech and Signal Processing.

[6]  Gerald Enzner,et al.  3D reconstruction of HRTF-fields from 1D continuous measurements , 2011, 2011 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics (WASPAA).

[7]  D. Rockmore,et al.  FFTs on the Rotation Group , 2008 .

[8]  V. Rokhlin,et al.  Fast Fourier Transforms for Nonequispaced Data, II , 1995 .

[9]  Sean S. B. Moore,et al.  FFTs for the 2-Sphere-Improvements and Variations , 1996 .

[10]  C. Avendano,et al.  The CIPIC HRTF database , 2001, Proceedings of the 2001 IEEE Workshop on the Applications of Signal Processing to Audio and Acoustics (Cat. No.01TH8575).

[11]  Griffin D. Romigh Individualized Head-Related Transfer Functions: Efficient Modeling and Estimation from Small Sets of Spatial Samples , 2012 .

[12]  Rodney A. Kennedy,et al.  On High-Resolution Head-Related Transfer Function Measurements: An Efficient Sampling Scheme , 2012, IEEE Transactions on Audio, Speech, and Language Processing.

[13]  Anthony I. Tew,et al.  Analyzing head-related transfer function measurements using surface spherical harmonics , 1998 .

[14]  O. V. Verkhodanov,et al.  Gauss - Legendre sky pixelization (GLESP) for CMB maps , 2005 .

[15]  Ramani Duraiswami,et al.  Regularized HRTF fitting using spherical harmonics , 2009, 2009 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics.

[16]  Douglas Brungart,et al.  Spectral HRTF enhancement for improved vertical-polar auditory localization , 2009, 2009 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics.

[17]  K. Gorski,et al.  HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere , 2004, astro-ph/0409513.

[18]  Robert G. Crittenden,et al.  Exactly azimuthal pixelizations of the sky , 1998 .

[19]  W. Shukowsky,et al.  A quadrature formula over the sphere with application to high resolution spherical harmonic analysis , 1986 .

[20]  Bosun Xie,et al.  Recovery of individual head-related transfer functions from a small set of measurements. , 2012, The Journal of the Acoustical Society of America.

[21]  Bill Gardner,et al.  HRTF Measurements of a KEMAR Dummy-Head Microphone , 1994 .

[22]  Wen Zhang,et al.  Binaural sound source localization using the frequency diversity of the head-related transfer function. , 2014, The Journal of the Acoustical Society of America.

[23]  Rodney A. Kennedy,et al.  An Optimal-Dimensionality Sampling Scheme on the Sphere With Fast Spherical Harmonic Transforms , 2014, IEEE Transactions on Signal Processing.

[24]  Akio Ando,et al.  Estimation of individualized head-related transfer function based on principal component analysis , 2009 .

[25]  Ramani Duraiswami,et al.  Flexible and Optimal Design of Spherical Microphone Arrays for Beamforming , 2007, IEEE Transactions on Audio, Speech, and Language Processing.

[26]  Yong Fang,et al.  Interpolation of head-related transfer functions using spherical fourier expansion , 2009 .

[27]  Yôiti Suzuki,et al.  Accuracy of head-related transfer functions synthesized with spherical microphone arrays , 2013 .

[28]  D. Healy,et al.  Computing Fourier Transforms and Convolutions on the 2-Sphere , 1994 .

[29]  Ramani Duraiswami,et al.  Fast head-related transfer function measurement via reciprocity. , 2006, The Journal of the Acoustical Society of America.

[30]  Paweł Strumiłło,et al.  Head related transfer functions measurement and processing for the purpose of creating a spatial sound environment , 2007, Symposium on Photonics Applications in Astronomy, Communications, Industry, and High-Energy Physics Experiments (WILGA).

[31]  R. Duraiswami,et al.  Insights into head-related transfer function: Spatial dimensionality and continuous representation. , 2010, The Journal of the Acoustical Society of America.

[32]  Don Coppersmith,et al.  Matrix multiplication via arithmetic progressions , 1987, STOC.

[33]  Christof Faller,et al.  Sound Field Analysis along a Circle and Its Applications to HRTF Interpolation , 2008 .

[34]  R. Duda,et al.  Range dependence of the response of a spherical head model , 1998 .

[35]  V R Algazi,et al.  Elevation localization and head-related transfer function analysis at low frequencies. , 2001, The Journal of the Acoustical Society of America.

[36]  Bosun Xie,et al.  On the spatial symmetry of head-related transfer functions , 2013 .

[37]  Ramani Duraiswami,et al.  Psychological and Physiological Acoustics Session 1 pPPa : Binaural Hearing and Binaural Techniques I 1 pPPa 3 . Statistical analysis of head related transfer function ( HRTF ) data , 2013 .

[38]  K. Gröchenig A discrete theory of irregular sampling , 1993 .