s2let_demo1

 s2let_demo1
 Analyse Earth tomography data as a real MW map.
 Compute the wavelet maps and plot them.
 Plot 1 : multiresolution wavelet scales
 Plot 2 : full resolution wavelet scales

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

0001 % s2let_demo1
0002 % Analyse Earth tomography data as a real MW map.
0003 % Compute the wavelet maps and plot them.
0004 % Plot 1 : multiresolution wavelet scales
0005 % Plot 2 : full resolution wavelet scales
0006 %
0007 % S2LET package to perform Wavelets on the Sphere.
0008 % Copyright (C) 2012-2015  Boris Leistedt & Jason McEwen
0009 % See LICENSE.txt for license details
0010 
0011 load('EGM2008_Topography_flms_L0128');
0012 f = ssht_inverse(flm, L, 'Reality', true);
0013 
0014 %inputfile = 'data/earth_tomo_mw_128.fits';
0015 %[f, L] = s2let_mw_read_real_map(inputfile);
0016 
0017 B = 3;
0018 J_min = 2;
0019 J = s2let_jmax(L, B);
0020 
0021 zoomfactor = 1.2;
0022 ns = ceil(sqrt(2+J-J_min+1)) ;
0023 ny = ns - 1 + rem(2+J-J_min+1 , ns) ;
0024 nx = ns;
0025 
0026 % Perform decomposition
0027 [f_wav, f_scal] = s2let_transform_axisym_analysis_mw(f, 'B', B, 'J_min', J_min, 'Reality', true, 'Upsample', true);
0028 % FULL RESOLUTION PLOT
0029 figure('Position',[100 100 1300 1000])
0030 subplot(nx, ny, 1);
0031 ssht_plot_mollweide(f, L);
0032 %title('Initial data')
0033 campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0034 v = caxis;
0035 temp = max(abs(v));
0036 caxis([-temp temp])
0037 subplot(nx, ny, 2);
0038 ssht_plot_mollweide(f_scal, L);
0039 campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0040 v = caxis;
0041 temp = max(abs(v));
0042 caxis([-temp temp])
0043 %title('Scaling fct')
0044 for j = J_min:J
0045    subplot(nx, ny, j-J_min+3);
0046    ssht_plot_mollweide(f_wav{j-J_min+1}, L);
0047    campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0048 v = caxis;
0049 temp = max(abs(v));
0050 caxis([-temp temp])
0051    %title(['Wavelet scale : ',int2str(j)-J_min+1])
0052 end
0053 
0054 
0055 % Perform decomposition
0056 [f_wav, f_scal] = s2let_transform_axisym_analysis_mw(f, 'B', B, 'J_min', J_min, 'Reality', true);
0057 
0058 % MULTIRESOLUTION PLOT
0059 figure('Position',[100 100 1300 1000])
0060 subplot(nx, ny, 1);
0061 ssht_plot_mollweide(f, L);
0062 %title('Initial data')
0063 campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0064 v = caxis;
0065 temp = max(abs(v));
0066 caxis([-temp temp])
0067 subplot(nx, ny, 2);
0068 bl = min([ s2let_bandlimit(J_min-1,J_min,B,L) L ]);
0069 ssht_plot_mollweide(f_scal, bl);
0070 campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0071 v = caxis;
0072 temp = max(abs(v));
0073 caxis([-temp temp])
0074 %title('Scaling fct')
0075 for j = J_min:J
0076    subplot(nx, ny, j-J_min+3);
0077    bl =  min([ s2let_bandlimit(j,J_min,B,L) L ]);
0078    ssht_plot_mollweide(f_wav{j-J_min+1}, bl);
0079    campos([0 0 -1]); camup([0 1 0]); zoom(zoomfactor)
0080 v = caxis;
0081 temp = max(abs(v));
0082 caxis([-temp temp])
0083    %title(['Wavelet scale : ',int2str(j)-J_min+1])
0084 end