s2let_demo_ridgelet_plot

 s2let_demo_ridgelet_plot
 Compute and plot the ridgelet wavelet and scaling functions.

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

0001 
0002 % s2let_demo_ridgelet_plot
0003 % Compute and plot the ridgelet wavelet and scaling functions.
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 % close all;
0010 
0011 L = 256
0012 B = 2;
0013 J_min = 3;
0014 J = s2let_jmax(L, B);
0015 sampling_method = 'MWSS';
0016 % sphere_plot_type = 'parametric';
0017 sphere_plot_type = 'colour';
0018 sphere_plot_lighting = true;
0019 sphere_plot_scale = 1;
0020 save_plots = false;
0021 spin = 2
0022 reality = false
0023 
0024 % plot_component = @(x) (x);
0025 plot_component = @(x) real(x);
0026 % plot_component = @(x) imag(x);
0027 % plot_component = @(x) abs(x);
0028 
0029 [ridgelet_wav, ridgelet_scal] = s2let_ridgelet_compute_wav(L, ...
0030    'B', B, 'J_min', J_min, ...
0031    'Spin', spin, 'Reality', reality, 'Sampling', sampling_method);
0032 
0033 figure
0034 ssht_plot_sphere(plot_component(ridgelet_scal), L, ...
0035    'Lighting', sphere_plot_lighting, ...
0036    'Type', sphere_plot_type, 'ParametricScale', sphere_plot_scale, ...
0037    'Method', sampling_method, ...
0038    'ColourBar', true);
0039 c = caxis;
0040 caxis([-1 1].*max(abs(c)));
0041 view(-37.5,10)
0042 if save_plots, print('-r300', '-dpng', 'ridgelet_scal_Jmin03.png'); end
0043 
0044 for j = J_min:J
0045 
0046    figure
0047    ssht_plot_sphere(plot_component(ridgelet_wav{j-J_min+1}), L, ...
0048       'Lighting', sphere_plot_lighting, ...
0049       'Type', sphere_plot_type, 'ParametricScale', sphere_plot_scale, ...
0050       'ParametricMin', true, ...
0051       'Method', sampling_method, ...
0052       'ColourBar', true);
0053    c = caxis;
0054    caxis([-1 1].*max(abs(c)));
0055    view(-37.5,10)
0056    
0057    plot_filename = sprintf('plots/ridgelet_wav_j%2.2d.png', j);
0058    if save_plots, print('-r300', '-dpng', plot_filename); end
0059    
0060 end
0061