The Auditory Modeling Toolbox

Applies to version: 0.10.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.