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 (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 (0: off, 1: on (default values: [0.05;0.01]), [alpha r]: rejects 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

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

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

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.