supdeq_gauss

% SUpDEq - Spatial Upsampling by Directional Equalization

 function [gridData, Npoints, Nmax] = supdeq_gauss(N)

 This function returns a Gaussian sampling grid with azimuth and elevation
 in degree plus the respective sampling weights.

 Output:
 gridData      - Q x 3 matrix where the first column holds the azimuth, the 
                 second the elevation, and the third the sampling weights.
 Npoints       - Total number of sampling points / nodes [2*(N+1)^2]
 Nmax          - Highest stable grid order  

 Input:
 N             - Spatial order of the gauss grid. Twice as many samples 
                 are always distributed along the azimuth compared to 
                 the elevation.

 Dependencies: SOFiA toolbox

 (C) 2018 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 [gridData, Npoints, Nmax] = supdeq_gauss(N)
0004 %
0005 % This function returns a Gaussian sampling grid with azimuth and elevation
0006 % in degree plus the respective sampling weights.
0007 %
0008 % Output:
0009 % gridData      - Q x 3 matrix where the first column holds the azimuth, the
0010 %                 second the elevation, and the third the sampling weights.
0011 % Npoints       - Total number of sampling points / nodes [2*(N+1)^2]
0012 % Nmax          - Highest stable grid order
0013 %
0014 % Input:
0015 % N             - Spatial order of the gauss grid. Twice as many samples
0016 %                 are always distributed along the azimuth compared to
0017 %                 the elevation.
0018 %
0019 % Dependencies: SOFiA toolbox
0020 %
0021 % (C) 2018 by JMA, Johannes M. Arend
0022 %             TH Köln - University of Applied Sciences
0023 %             Institute of Communications Engineering
0024 %             Department of Acoustics and Audio Signal Processing
0025 
0026 function [gridData, Npoints, Nmax] = supdeq_gauss(N)
0027 
0028 if nargin == 0
0029     error('Please specify the desired spatial order N!');
0030 else
0031     AZnodes = 2*(N+1);
0032     ELnodes = AZnodes/2;
0033     [gridData, Npoints, Nmax] = sofia_gauss(AZnodes,ELnodes,0);
0034     %Convert from rad 2 deg
0035     gridData(:,1:2) = gridData(:,1:2)*180/pi;
0036 end
0037 
0038 end
0039