s2let_demo3

 s2let_demo3
 Analyse some CMB simulation from a valid HEALPIX map
 Compute the output wavelets as Healpix maps and plot them.

 S2LET package to perform Wavelets on the Sphere.
 Copyright (C) 2012-2015  Boris Leistedt & Jason McEwen
 See LICENSE.txt for license details

0001 % s2let_demo3
0002 % Analyse some CMB simulation from a valid HEALPIX map
0003 % Compute the output wavelets as Healpix maps and plot them.
0004 %
0005 % S2LET package to perform Wavelets on the Sphere.
0006 % Copyright (C) 2012-2015  Boris Leistedt & Jason McEwen
0007 % See LICENSE.txt for license details
0008 
0009 L = 192;
0010 B = 3;
0011 J_min = 2;
0012 nside_recon = 128;
0013 s2let_path = '.'
0014 inputfile = strcat(s2let_path,'/data/somecmbsimu_hpx_128.fits')
0015 
0016 % Read the file
0017 [f_ini, nside] = s2let_hpx_read_real_map(inputfile);
0018 
0019 % Band limit the data
0020 flm = s2let_hpx_map2alm(f_ini, 'L', L);
0021 f = s2let_hpx_alm2map(flm, nside_recon, 'L', L);
0022 
0023 % Perform decomposition
0024 [f_wav, f_scal] = s2let_transform_axisym_analysis_hpx(f,'B',B,'L',L,'J_min',J_min);
0025 
0026 % Plot
0027 J = s2let_jmax(L, B);
0028 zoomfactor = 1.2;
0029 ns = ceil(sqrt(2+J-J_min+1)) ;
0030 ny = ns - 1 ;
0031 nx = ns ;
0032 figure('Position',[100 100 1300 1000])
0033 
0034 subplot(nx, ny, 1);
0035 s2let_hpx_plot_mollweide(f);
0036 campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0037 title('Initial band-limited data')
0038 
0039 subplot(nx, ny, 2);
0040 s2let_hpx_plot_mollweide(f_scal);
0041 campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0042 title('Scaling fct')
0043 
0044 for j = J_min:J
0045    subplot(nx, ny, j-J_min+3);
0046    s2let_hpx_plot_mollweide(f_wav{j-J_min+1});
0047    campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0048    title(['Wavelet scale : ',int2str(j)-J_min+1])
0049 end