THE AUDITORY MODELING TOOLBOX

Applies to version: 0.9.5

View the code

Go to function

ZIEGELWANGER2014 - Time of arrival estimates

Usage

[Obj,results]=ziegelwanger2014(data,estimation,outlierDetection,model,p0_onaxis)

Input parameters

Obj SOFA object
estimation (optional)
select one of the estimation methods

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

Block quote ends without a blank line; unexpected unindent.

1: Threshold-Detection 2: Centroid of squared IR 3: Mean Groupdelay 4: Minimal-Phase Cross-Correlation (Max) (default) [TOAest]: pre-estimated TOAs

outlierDetection (optional)
detect outliers in estimated TOAs

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

Block quote ends without a blank line; unexpected unindent.

0: off 1: on (default values: [0.05;0.01]) [alpha r]: rejects outliers using the extreme Studentized

System Message: ERROR/3 (<string>, line 33)

Unexpected indentation.
deviance test with the significance level of ALPHA and upper bound of outlier rate R.
model (optional)
correct estimated toa, using geometrical TOA-Model

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

Block quote ends without a blank line; unexpected unindent.

0: TOA estimated 1: off-axis TOA modeled (default) 2: on-axis TOA modeled

p0_onaxis (optional)
startvalues for lsqcurvefit

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

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 48)

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 59)

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 63)

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 67)

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 70)

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, azimut of ear in radiants, elevation of ear in radiants]

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

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.

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/ziegelwanger2014

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 estimation, use:

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

References:

H. Ziegelwanger and P. Majdak. Modeling the direction-continuous time-of-arrival in head-related transfer functions. J. Acoust. Soc. Am., 135:1278-1293, 2014.