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:

SOFA API from http://sourceforge.net/projects/sofacoustics for Matlab (in e.g. thirdparty/SOFA)
Optimization Toolbox for Matlab
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.