s2let_axisym_synthesis Compute axisymmetric wavelet transform, input/outputs as HEALPIX maps. Default usage : f = s2let_transform_axisym_synthesis_hpx(f_wav, f_scal, <options>) f_wav contains the input wavelet contributions -- HEALPIX sampling, f_scal contains the input scaling contributions -- HEALPIX sampling, f is the output field -- HEALPIX sampling, Option : 'nside' = { HEALPIX resolution; (default=guessed)} 'B' = { Dilation factor; B > 1 (default=2) } 'L' = { Harmonic band-limit; L > 1 (default=2*nside) } 'J_min' = { Minimum wavelet scale to consider; 0 <= J_min < log_B(L) (default=0) } S2LET package to perform Wavelets transform on the Sphere. Copyright (C) 2012-2015 Boris Leistedt & Jason McEwen See LICENSE.txt for license details
0001 function f = s2let_transform_axisym_synthesis_hpx(f_wav, f_scal, varargin) 0002 0003 % s2let_axisym_synthesis 0004 % Compute axisymmetric wavelet transform, input/outputs as HEALPIX maps. 0005 % 0006 % Default usage : 0007 % 0008 % f = s2let_transform_axisym_synthesis_hpx(f_wav, f_scal, <options>) 0009 % 0010 % f_wav contains the input wavelet contributions -- HEALPIX sampling, 0011 % f_scal contains the input scaling contributions -- HEALPIX sampling, 0012 % f is the output field -- HEALPIX sampling, 0013 % 0014 % Option : 0015 % 'nside' = { HEALPIX resolution; (default=guessed)} 0016 % 'B' = { Dilation factor; B > 1 (default=2) } 0017 % 'L' = { Harmonic band-limit; L > 1 (default=2*nside) } 0018 % 'J_min' = { Minimum wavelet scale to consider; 0019 % 0 <= J_min < log_B(L) (default=0) } 0020 % 0021 % S2LET package to perform Wavelets transform on the Sphere. 0022 % Copyright (C) 2012-2015 Boris Leistedt & Jason McEwen 0023 % See LICENSE.txt for license details 0024 0025 0026 len = size(f_wav); 0027 temp = f_wav{len}; 0028 sz = size(temp); 0029 nsideguessed = sqrt(max(sz)/12); 0030 Lguessed = 2*nsideguessed; 0031 0032 p = inputParser; 0033 p.addRequired('f_wav'); 0034 p.addRequired('f_scal', @isnumeric); 0035 p.addParamValue('nside', nsideguessed, @isnumeric); 0036 p.addParamValue('B', 2, @isnumeric); 0037 p.addParamValue('L', 2*nsideguessed, @isnumeric); 0038 p.addParamValue('J_min', 0, @isnumeric); 0039 p.parse(f_wav, f_scal, varargin{:}); 0040 args = p.Results; 0041 0042 0043 J = s2let_jmax(args.L, args.B); 0044 f_wav_vec = []; 0045 for j = args.J_min:J 0046 temp = f_wav{j+1-args.J_min}; 0047 f_wav_vec = [f_wav_vec temp]; 0048 end 0049 0050 f = s2let_transform_axisym_synthesis_hpx_mex(f_wav_vec, f_scal, args.nside, args.B, args.L, args.J_min); 0051 0052 end