THE AUDITORY MODELING TOOLBOX

Applies to version: 0.9.5

View the code

Go to function

ZIEGELWANGER2013 - Time of arrival estimates

Usage

[Obj,results]=ziegelwanger2013(Obj,method,model,p0_onaxis)

Input parameters

Obj SOFA object
method (Optional) Select one of the estimation methods
model (Optional) Correct estimated toa, using geometrical TOA-Model. If \(model=0\) use TOA estimated, if \(model=1\) use TOA modeled (default).
p0_onaxis
(optional) Starting values for lsqcurvefit

System Message: WARNING/2 (<string>, line 27)

Block quote ends without a blank line; unexpected unindent.
dim 1: [sphere-radius in m,
azimut of ear in radiants, elevation of ear in radiants, direction-independent delay in seconds]

System Message: WARNING/2 (<string>, line 31)

Definition list ends without a blank line; unexpected unindent.

dim 2: each record channel

Output parameters

Obj SOFA Object
results.toa
data matrix with time of arrival (TOA) for each impulse response (IR):

System Message: WARNING/2 (<string>, line 42)

Block quote ends without a blank line; unexpected unindent.

dim 1: each toa in samples dim 2: each record channel

results.p_onaxis
estimated on-axis model-parameters

System Message: WARNING/2 (<string>, line 46)

Block quote ends without a blank line; unexpected unindent.
dim 1: [sphere-radius in m,
azimut of ear in radiants, elevation of ear in radiants, direction-independent delay in seconds]

System Message: WARNING/2 (<string>, line 50)

Definition list ends without a blank line; unexpected unindent.

dim 2: each record channel

results.p_offaxis
estimated off-axis model-parameters

System Message: WARNING/2 (<string>, line 53)

Block quote ends without a blank line; unexpected unindent.
dim 1: [sphere-radius in m,
xM in m, yM in m, zM in m, direction-independent delay in seconds, channel (starting at 1), azimut of ear in radiants, elevation of ear in radiants]

System Message: WARNING/2 (<string>, line 61)

Definition list ends without a blank line; unexpected unindent.

dim 2: each record channel

Description

Estimates the Time-of-Arrival for each measurement in Obj (SOFA) and corrects the results with a geometrical model of the head.

The value of method is an integer choosing one of the following methods. XXX Explain for each method a little about how they work:

  1. Threshold-Detection
  2. Centroid of squared IR
  3. Mean Groupdelay
  4. Minimal-Phase Cross-Correlation (Max) (default)

Requirements:

  1. SOFA API from http://sourceforge.net/projects/sofacoustics for Matlab (in e.g. thirdparty/SOFA)
  2. Optimization Toolbox for Matlab
  3. Data in hrtf/ziegelwanger2013

Examples:

To calculate the model parameters for the on-axis time-of-arrival model (p_onaxis) and for the off-axis time-of-arrival model (p_offaxis) for a given HRTF set (SOFA object, 'Obj') with the minimum-phase cross-correlation method, use:

[Obj,results]=ziegelwanger2013(Obj,4,1);

References:

P. Majdak and H. Ziegelwanger. Continuous-direction model of the broadband time-of-arrival in the head-related transfer functions. In ICA 2013 Montreal, volume 19, page 050016, Montreal, Canada, 2013. ASA. [ DOI ]

H. Ziegelwanger and P. Majdak. Modeling the broadband time-of-arrival of the head-related transfer functions for binaural audio. In Proceedings of the 134th Convention of the Audio Engineering Society, page 7, Rome, 2013.