s2let_hpx2mw

 s2let_hpx2mw 
 Converts Healpix map into MW map

 Default usage :

   f_mw = s2let_hpx2mw(f_hpx, <options>)

 f is the input map -- HEALPIX sampling,
 f is the output map -- MW sampling.

 Option :
  'nside'           = { HEALPIX resolution; (default=guessed)}
  'L'               = { Harmonic band-limit; L > 1 (default=2*nside) }

0001 function f_mw = s2let_hpx2mw(f, varargin)
0002 
0003 % s2let_hpx2mw
0004 % Converts Healpix map into MW map
0005 %
0006 % Default usage :
0007 %
0008 %   f_mw = s2let_hpx2mw(f_hpx, <options>)
0009 %
0010 % f is the input map -- HEALPIX sampling,
0011 % f is the output map -- MW sampling.
0012 %
0013 % Option :
0014 %  'nside'           = { HEALPIX resolution; (default=guessed)}
0015 %  'L'               = { Harmonic band-limit; L > 1 (default=2*nside) }
0016 
0017 % S2LET package to perform Wavelets transform on the Sphere.
0018 % Copyright (C) 2012-2015  Boris Leistedt & Jason McEwen
0019 % See LICENSE.txt for license details
0020 
0021 sz = size(f);
0022 nsideguessed = sqrt(max(sz)/12);
0023 Lguessed = 2*nsideguessed;
0024 
0025 p = inputParser;
0026 p.addRequired('f', @isnumeric); 
0027 p.addParamValue('nside', nsideguessed, @isnumeric);   
0028 p.addParamValue('L', Lguessed, @isnumeric);   
0029 p.parse(f, varargin{:});
0030 args = p.Results;
0031 
0032 flm = s2let_hpx_map2alm(f, 'nside', args.nside, 'L', args.L);
0033 %flm = s2let_hpx_map2alm_mex(f, args.nside, args.L);
0034 f_mw = ssht_inverse(flm, args.L, 'Reality', true);
0035 
0036 end