THE AUDITORY MODELING TOOLBOX

This documentation page applies to an outdated major AMT version. We show it for archival purposes only.
Click here for the documentation menu and here to download the latest AMT (1.6.0).

View the code

Go to function

ZIEGELWANGER2014 - Time of arrival estimates

Usage

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

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

Input parameters

'',Obj SOFA object
'',estimation

(optional) TOA estimation method:

[TOAest]: pre-estimated TOAs 1: Built-in threshold detection 2: Built-in centroid of squared IR, 3: Built-in mean group delay, 4: Built-in maximum of the minimum-phase cross-correlation (default) String: An estimator from the 'itdestimator' function ('MaxIACCe' recommended)
'',outlierDetection
(optional) detect outliers in estimated TOAs

0: off 1: on (default values: [0.05;0.01]) [alpha r]: reject outliers using the extreme Studentized deviance test with the significance level of ALPHA and upper bound of outlier rate R

'',model

(optional) correct estimated toa, using geometrical TOA-Model

0: TOA estimated, 1: off-axis TOA modeled (default), 2: on-axis TOA modeled
'',p0_onaxis (optional) startvalues for lsqcurvefit
'',lowpass (optional) bandwidth setting when used with the estimator from 'itdestimator' ('lp' for lowpass, 'bb' for broadband)
'',upper_cutfreq (optional) lowpass cutoff (Hz) when used with the estimator from 'itdestimator'
'',threshlvl (optional) threshold level (dB) when used with the 'Threshold' estimator from 'itdestimator'

Output parameters

'',Obj SOFA Object
'=results.toa'
data matrix with time of arrival (TOA) for each impulse response (IR)
'=results.p_onaxis'
estimated on-axis model-parameters
'=results.p_offaxis'
estimated off-axis model-parameters

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.