supdeq_optRadius

% SUpDEq - Spatial Upsampling by Directional Equalization

 function r_opt = supdeq_optRadius(headWidth, headHeight, headLength, method)

 This function calculates the optimal radius r (in m) as a linear 
 combination of head width, head height, and head depth according to 
 Algazi et al. [1] or Bahu & Romblom [2]

 Output:
 r             - Optimal radius r in m

 Input:
 headWidth     - Head width in m
 headHeight    - Head height in m
 headLength    - Head length/depth in m
 method        - 'Algazi' or 'Bahu'
                 'Algazi' is default

 Dependencies: -

 References:
 [1] - Algazi VR., Avendano C., Duda RO. - Estimation of a spherical-head 
 model from anthropometry. J. Audio Eng. Soc. 2001; 49(6):472-479.

 [2] H. Bahu and R. David, ?Optimization and prediction of the spherical 
 and ellipsoidal ITD model parameters using offset ears,? 
 in Proceedings of the AES International Conference 
 on Spatial Reproduction - Aesthetics and Science, 2018, pp. 1?11.
   
 (C) 2018-2020 by JMA, Johannes M. Arend
               TH Köln - University of Applied Sciences
               Institute of Communications Engineering
               Department of Acoustics and Audio Signal Processing

0001 %% SUpDEq - Spatial Upsampling by Directional Equalization
0002 %
0003 % function r_opt = supdeq_optRadius(headWidth, headHeight, headLength, method)
0004 %
0005 % This function calculates the optimal radius r (in m) as a linear
0006 % combination of head width, head height, and head depth according to
0007 % Algazi et al. [1] or Bahu & Romblom [2]
0008 %
0009 % Output:
0010 % r             - Optimal radius r in m
0011 %
0012 % Input:
0013 % headWidth     - Head width in m
0014 % headHeight    - Head height in m
0015 % headLength    - Head length/depth in m
0016 % method        - 'Algazi' or 'Bahu'
0017 %                 'Algazi' is default
0018 %
0019 % Dependencies: -
0020 %
0021 % References:
0022 % [1] - Algazi VR., Avendano C., Duda RO. - Estimation of a spherical-head
0023 % model from anthropometry. J. Audio Eng. Soc. 2001; 49(6):472-479.
0024 %
0025 % [2] H. Bahu and R. David, ?Optimization and prediction of the spherical
0026 % and ellipsoidal ITD model parameters using offset ears,?
0027 % in Proceedings of the AES International Conference
0028 % on Spatial Reproduction - Aesthetics and Science, 2018, pp. 1?11.
0029 %
0030 % (C) 2018-2020 by JMA, Johannes M. Arend
0031 %               TH Köln - University of Applied Sciences
0032 %               Institute of Communications Engineering
0033 %               Department of Acoustics and Audio Signal Processing
0034 
0035 function r_opt = supdeq_optRadius(headWidth, headHeight, headLength, method)
0036 
0037 if nargin < 4 || isempty(method)
0038     method = 'Algazi';
0039 end
0040 
0041 
0042 if strcmp(method,'Algazi')
0043     %r_opt = 0.51*X1 + 0.019*X2 + 0.18*X3 + 3.2 cm
0044     %width X1 = head half width, X2 = head half height, X3 = head half length
0045     r_opt = 0.51*(headWidth/2) + 0.019*(headHeight/2) + 0.18*(headLength/2) + 0.032;
0046 end
0047 
0048 if strcmp(method,'Bahu')
0049     %r_opt = 0.44*X1 + 0.23*X3 + 3.2 cm
0050     %width X1 = head half width, X2 = head half height, X3 = head half length
0051     r_opt = 0.44*(headWidth/2) + 0.23*(headLength/2) + 0.032;
0052 end
0053 
0054 end
0055