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).
Go to function
DATA_ZIEGELWANGER2014 - Data from Ziegelwanger and Majdak (2014)
Program code:
function data = data_ziegelwanger2014(varargin)
%DATA_ZIEGELWANGER2014 Data from Ziegelwanger and Majdak (2014)
% Usage: data = data_ziegelwanger2014(flag)
%
% `data_ziegelwanger2014(flag)` returns results for different HRTF
% databases from Ziegelwanger and Majdak (2014).
%
% The flag may be one of:
%
% 'ARI' ARI database. The output has the following
% fields: `data.results` and `data.subjects`.
%
% 'CIPIC' CIPIC database. The output has the following fields:
% `data.results` and `data.subjects`.
%
% 'LISTEN' LISTEN database. The output has the following fields.
% `data.results` and `data.subjects`.
%
% 'SPHERE_ROT' HRTF sets for a rigid sphere placed in the center of
% the measurement setup and varying rotation. The
% output has the following fields: `data.results`,
% `data.subjects`, `data.phi`, `data.theta` and `data.radius`.
%
% 'SPHERE_DIS' HRTF sets for a rigid sphere with various positions in
% the measurement setup. The output has the following fields:
% `data.results`, `data.subjects`, `data.xM`, `data.yM`,
% `data.zM` and `data.radius`.
%
% 'Sphere' HRTF set for a rigid sphere: The
% output has the following fields: `data.hM`,
% `data.meta` and `data.stimPar`.
%
% 'SAT' HRTF set for a rigid sphere combined with a neck and a
% torso: The output has the following fields: `data.hM`,
% `data.meta` and `data.stimPar`.
%
% 'STP' HRTF set for a rigid sphere combined with a neck, a
% torso and a pinna: The output has the following fields:
% `data.hM`,`data.meta` and `data.stimPar`.
%
% 'NH89' HRTF set of listener NH89 of the ARI database: The
% output has the following fields: `data.hM`,
% `data.meta` and `data.stimPar`.
%
% 'reload' Reload previously calculated results
%
% 'recalc' Recalculate the results
%
% The fields are given by:
%
% `data.results` Results for all HRTF sets
%
% `data.subjects` IDs for HRTF sets
%
% `data.phi` Azimuth of ear position
%
% `data.theta` Elevation of ear position
%
% `data.radius` sphere radius
%
% `data.xM` x-coordinate of sphere center
%
% `data.yM` y-coordinate of sphere center
%
% `data.zM` z-coordinate of sphere center
%
% `data` SOFA object
%
% Examples:
% ---------
%
% To get results from the ARI database, use::
%
% data=data_ziegelwanger2014('ARI');
%
% See also: ziegelwanger2014, ziegelwanger2014onaxis,
% ziegelwanger2014offaxis, exp_ziegelwanger2014
%
% References: ziegelwanger2014
% AUTHOR: Harald Ziegelwanger, Acoustics Research Institute, Vienna,
% Austria
%% ------ Check input options --------------------------------------------
% Define input flags
definput.flags.type = {'missingflag','ARI','CIPIC','LISTEN','SPHERE_DIS','SPHERE_ROT','NH89','Sphere','SAT','STP','redo'};
definput.flags.results = {'reload','recalc'};
% Parse input options
[flags,~] = ltfatarghelper({},definput,varargin);
if flags.do_missingflag
flagnames=[sprintf('%s, ',definput.flags.type{2:end-2}),...
sprintf('%s or %s',definput.flags.type{end-1},definput.flags.type{end})];
error('%s: You must specify one of the following flags: %s.',upper(mfilename),flagnames);
else
hpath = which('hrtfinit'); % find local path of hrtf repository
hpath = hpath(1:end-10);
hpath = [hpath 'ziegelwanger2014' filesep];
if ~exist([hpath 'info.mat'],'file')
urlwrite([SOFAdbURL filesep 'ziegelwanger2014/info.mat'],[hpath 'info.mat']);
end
end
%% ARI database
if flags.do_ARI
if ~flags.do_recalc && ~exist([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_ARI.mat'],'file')
urlwrite([SOFAdbURL filesep 'ziegelwanger2014' filesep 'exp_ziegelwanger2014_ARI.mat'],[hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_ARI.mat']);
end
if flags.do_recalc
tmp=load([hpath 'info.mat']);
data=tmp.info.ARI;
for ii=1:length(data.subjects)
disp(['Recalculate data for sujbect ' num2str(ii) filesep num2str(length(data.subjects)) ' of ARI database']);
Obj=SOFAload([hpath 'ARI_' data.subjects{ii} '.sofa']);
if exist([hpath 'ARI_' data.subjects{ii} '.sofa.MCM.mat'],'file')
load([hpath 'ARI_' data.subjects{ii} '.sofa.MCM.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,4,0,0);
toaEst=tmp.toa;
save([hpath 'ARI_' data.subjects{ii} '.sofa.MCM.mat'],'toaEst');
end
[~,results(ii).MCM{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
[~,results(ii).MCM{2}]=ziegelwanger2014(Obj,toaEst,[0.05 0.01],1e-8);
end
save([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_ARI.mat'],'results');
data.results=results;
else
tmp=load([hpath 'info.mat']);
data=tmp.info.ARI;
load([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_ARI.mat']);
data.results=results;
end
end
%% CIPIC database
if flags.do_CIPIC
if ~flags.do_recalc && ~exist([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_CIPIC.mat'],'file')
urlwrite([SOFAdbURL filesep 'ziegelwanger2014' filesep 'exp_ziegelwanger2014_CIPIC.mat'],[hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_CIPIC.mat']);
end
if flags.do_recalc
tmp=load([hpath 'info.mat']);
data=tmp.info.CIPIC;
for ii=1:length(data.subjects)
disp(['Recalculate data for sujbect ' num2str(ii) filesep num2str(length(data.subjects)) ' of CIPIC database']);
Obj=SOFAload([hpath 'CIPIC_' data.subjects{ii} '.sofa']);
if exist([hpath 'CIPIC_' data.subjects{ii} '.sofa.MCM.mat'],'file')
load([hpath 'CIPIC_' data.subjects{ii} '.sofa.MCM.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,4,0,0);
toaEst=tmp.toa;
save([hpath 'CIPIC_' data.subjects{ii} '.sofa.MCM.mat'],'toaEst');
end
[~,results(ii).MCM{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
[~,results(ii).MCM{2}]=ziegelwanger2014(Obj,toaEst,[0.05 0.01],1e-8);
end
save([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_CIPIC.mat'],'results');
data.results=results;
else
tmp=load([hpath 'info.mat']);
data=tmp.info.CIPIC;
load([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_CIPIC.mat']);
data.results=results;
end
end
%% LISTEN database
if flags.do_LISTEN
if ~flags.do_recalc && ~exist([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_LISTEN.mat'],'file')
urlwrite([SOFAdbURL filesep 'ziegelwanger2014' filesep 'exp_ziegelwanger2014_LISTEN.mat'],[hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_LISTEN.mat']);
end
if flags.do_recalc
tmp=load([hpath 'info.mat']);
data=tmp.info.LISTEN;
for ii=1:length(data.subjects)
if ~strcmp(data.subjects{ii},'34')
disp(['Recalculate data for subject ' num2str(ii) filesep num2str(length(data.subjects)) ' of LISTEN database']);
Obj=SOFAload([hpath 'LISTEN_' data.subjects{ii} '.sofa']);
if exist([hpath 'LISTEN_' data.subjects{ii} '.sofa.MCM.mat'],'file')
load([hpath 'LISTEN_' data.subjects{ii} '.sofa.MCM.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,4,0,0);
toaEst=tmp.toa;
save([hpath 'LISTEN_' data.subjects{ii} '.sofa.MCM.mat'],'toaEst');
end
[~,results(ii).MCM{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
[~,results(ii).MCM{2}]=ziegelwanger2014(Obj,toaEst,[0.05 0.01],1e-8);
end
end
save([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_LISTEN.mat'],'results');
data.results=results;
else
tmp=load([hpath 'info.mat']);
data=tmp.info.LISTEN;
load([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_LISTEN.mat']);
data.results=results;
end
end
%% SPHERE (Displacement) database
if flags.do_SPHERE_DIS
if ~flags.do_recalc && ~exist([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_DIS.mat'],'file')
urlwrite([SOFAdbURL filesep 'ziegelwanger2014' filesep 'exp_ziegelwanger2014_SPHERE_DIS.mat'],[hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_DIS.mat']);
end
if flags.do_recalc
tmp=load([hpath 'info.mat']);
data=tmp.info.Sphere_Displacement;
results.p_onaxis=zeros(4,2,length(data.subjects));
results.p_offaxis=zeros(7,2,length(data.subjects));
for ii=1:length(data.subjects)
disp(['Recalculate data for subject ' num2str(ii) filesep num2str(length(data.subjects)) ' of SPHERE_DIS database']);
Obj=SOFAload([hpath 'Sphere_Displacement_' data.subjects{ii} '.sofa']);
if exist([hpath 'Sphere_Displacement_' data.subjects{ii} '.sofa.MCM.mat'],'file')
load([hpath 'Sphere_Displacement_' data.subjects{ii} '.sofa.MCM.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,4,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere_Displacement_' data.subjects{ii} '.sofa.MCM.mat'],'toaEst');
end
[~,results(ii).MCM{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
[~,results(ii).MCM{2}]=ziegelwanger2014(Obj,toaEst,[0.05 0.01],1e-8);
end
save([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_DIS.mat'],'results');
data.results=results;
else
tmp=load([hpath 'info.mat']);
data=tmp.info.Displacement;
load([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_DIS.mat']);
data.results=results;
end
end
%% SPHERE (Rotation) database
if flags.do_SPHERE_ROT
if ~flags.do_recalc && ~exist([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_ROT.mat'],'file')
urlwrite([SOFAdbURL filesep 'ziegelwanger2014' filesep 'exp_ziegelwanger2014_SPHERE_ROT.mat'],[hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_ROT.mat']);
end
if flags.do_recalc
tmp=load([hpath 'info.mat']);
data=tmp.info.Sphere_Rotation;
results.p=zeros(4,2,length(data.phi));
for ii=1:length(data.subjects)
disp(['Recalculate data for subject ' num2str(ii) filesep num2str(length(data.subjects)) ' of SPHERE_ROT database']);
Obj=SOFAload([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa']);
if exist([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.MAX.mat'],'file')
load([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.MAX.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,1,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.MAX.mat'],'toaEst');
end
[~,results(ii).MAX{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
if exist([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.CTD.mat'],'file')
load([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.CTD.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,2,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.CTD.mat'],'toaEst');
end
[~,results(ii).CTD{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
if exist([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.AGD.mat'],'file')
load([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.AGD.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,3,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.AGD.mat'],'toaEst');
end
[~,results(ii).AGD{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
if exist([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.MCM.mat'],'file')
load([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.MCM.mat']);
else
[~,tmp]=ziegelwanger2014(Obj,4,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere_Rotation_' data.subjects{ii} '.sofa.MCM.mat'],'toaEst');
end
[~,results(ii).MCM{1}]=ziegelwanger2014(Obj,toaEst,0,1e-8);
[~,results(ii).MCM{2}]=ziegelwanger2014(Obj,toaEst,[0.05 0.01],1e-8);
end
save([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_ROT.mat'],'results');
data.results=results;
else
tmp=load([hpath 'info.mat']);
data=tmp.info.Rotation;
load([hpath '..' filesep '..' filesep 'experiments' filesep 'exp_ziegelwanger2014_SPHERE_ROT.mat']);
data.results=results;
end
end
%% ARI database (NH89)
if flags.do_NH89
data=SOFAload([hpath 'ARI_NH89.sofa']);
if exist([hpath 'ARI_NH89.sofa.MAX.mat'],'file')
tmp=load([hpath 'ARI_NH89.sofa.MAX.mat']);
data.Data.toaEst{1}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,1,0,0);
toaEst=tmp.toa;
save([hpath 'ARI_NH89.sofa.MAX.mat'],'toaEst');
data.Data.toaEst{1}=toaEst;
end
if exist([hpath 'ARI_NH89.sofa.CTD.mat'],'file')
tmp=load([hpath 'ARI_NH89.sofa.CTD.mat']);
data.Data.toaEst{2}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,2,0,0);
toaEst=tmp.toa;
save([hpath 'ARI_NH89.sofa.CTD.mat'],'toaEst');
data.Data.toaEst{2}=toaEst;
end
if exist([hpath 'ARI_NH89.sofa.AGD.mat'],'file')
tmp=load([hpath 'ARI_NH89.sofa.AGD.mat']);
data.Data.toaEst{3}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,3,0,0);
toaEst=tmp.toa;
save([hpath 'ARI_NH89.sofa.AGD.mat'],'toaEst');
data.Data.toaEst{3}=toaEst;
end
if exist([hpath 'ARI_NH89.sofa.MCM.mat'],'file')
tmp=load([hpath 'ARI_NH89.sofa.MCM.mat']);
data.Data.toaEst{4}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,4,0,0);
toaEst=tmp.toa;
save([hpath 'ARI_NH89.sofa.MCM.mat'],'toaEst');
data.Data.toaEst{4}=toaEst;
end
end
%% Sphere
if flags.do_Sphere
data=SOFAload([hpath 'Sphere.sofa']);
if exist([hpath 'Sphere.sofa.MAX.mat'],'file')
tmp=load([hpath 'Sphere.sofa.MAX.mat']);
data.Data.toaEst{1}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,1,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere.sofa.MAX.mat'],'toaEst');
data.Data.toaEst{1}=toaEst;
end
if exist([hpath 'Sphere.sofa.CTD.mat'],'file')
tmp=load([hpath 'Sphere.sofa.CTD.mat']);
data.Data.toaEst{2}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,2,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere.sofa.CTD.mat'],'toaEst');
data.Data.toaEst{2}=toaEst;
end
if exist([hpath 'Sphere.sofa.AGD.mat'],'file')
tmp=load([hpath 'Sphere.sofa.AGD.mat']);
data.Data.toaEst{3}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,3,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere.sofa.AGD.mat'],'toaEst');
data.Data.toaEst{3}=toaEst;
end
if exist([hpath 'Sphere.sofa.MCM.mat'],'file')
tmp=load([hpath 'Sphere.sofa.MCM.mat']);
data.Data.toaEst{4}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,4,0,0);
toaEst=tmp.toa;
save([hpath 'Sphere.sofa.MCM.mat'],'toaEst');
data.Data.toaEst{4}=toaEst;
end
end
%% Sphere and Torso
if flags.do_SAT
data=SOFAload([hpath 'SAT.sofa']);
if exist([hpath 'SAT.sofa.MAX.mat'],'file')
tmp=load([hpath 'SAT.sofa.MAX.mat']);
data.Data.toaEst{1}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,1,0,0);
toaEst=tmp.toa;
save([hpath 'SAT.sofa.MAX.mat'],'toaEst');
data.Data.toaEst{1}=toaEst;
end
if exist([hpath 'SAT.sofa.CTD.mat'],'file')
tmp=load([hpath 'SAT.sofa.CTD.mat']);
data.Data.toaEst{2}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,2,0,0);
toaEst=tmp.toa;
save([hpath 'SAT.sofa.CTD.mat'],'toaEst');
data.Data.toaEst{2}=toaEst;
end
if exist([hpath 'SAT.sofa.AGD.mat'],'file')
tmp=load([hpath 'SAT.sofa.AGD.mat']);
data.Data.toaEst{3}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,3,0,0);
toaEst=tmp.toa;
save([hpath 'SAT.sofa.AGD.mat'],'toaEst');
data.Data.toaEst{3}=toaEst;
end
if exist([hpath 'SAT.sofa.MCM.mat'],'file')
tmp=load([hpath 'SAT.sofa.MCM.mat']);
data.Data.toaEst{4}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,4,0,0);
toaEst=tmp.toa;
save([hpath 'SAT.sofa.MCM.mat'],'toaEst');
data.Data.toaEst{4}=toaEst;
end
end
%% Sphere, Torso and Pinna
if flags.do_STP
data=SOFAload([hpath 'STP.sofa']);
if exist([hpath 'STP.sofa.MAX.mat'],'file')
tmp=load([hpath 'STP.sofa.MAX.mat']);
data.Data.toaEst{1}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,1,0,0);
toaEst=tmp.toa;
save([hpath 'STP.sofa.MAX.mat'],'toaEst');
data.Data.toaEst{1}=toaEst;
end
if exist([hpath 'STP.sofa.CTD.mat'],'file')
tmp=load([hpath 'STP.sofa.CTD.mat']);
data.Data.toaEst{2}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,2,0,0);
toaEst=tmp.toa;
save([hpath 'STP.sofa.CTD.mat'],'toaEst');
data.Data.toaEst{2}=toaEst;
end
if exist([hpath 'STP.sofa.AGD.mat'],'file')
tmp=load([hpath 'STP.sofa.AGD.mat']);
data.Data.toaEst{3}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,3,0,0);
toaEst=tmp.toa;
save([hpath 'STP.sofa.AGD.mat'],'toaEst');
data.Data.toaEst{3}=toaEst;
end
if exist([hpath 'STP.sofa.MCM.mat'],'file')
tmp=load([hpath 'STP.sofa.MCM.mat']);
data.Data.toaEst{4}=tmp.toaEst;
else
[~,tmp]=ziegelwanger2014(data,4,0,0);
toaEst=tmp.toa;
save([hpath 'STP.sofa.MCM.mat'],'toaEst');
data.Data.toaEst{4}=toaEst;
end
end
Build with Bootstrap