plot_axisym_tiling - Plot tiling in harmonic space. -- Axisymmetric wavelets on the sphere. Default usage : plot_axisym_tiling(B, L, J_min) B is the wavelet parameter, L is the angular band-limit, J_min the first wavelet to be used. S2LET package to perform Wavelet transform on the Sphere. Copyright (C) 2012-2015 Boris Leistedt & Jason McEwen See LICENSE.txt for license details
0001 function s2let_plot_axisym_tiling(B, L, J_min) 0002 0003 % plot_axisym_tiling - Plot tiling in harmonic space. 0004 % -- Axisymmetric wavelets on the sphere. 0005 % 0006 % Default usage : 0007 % 0008 % plot_axisym_tiling(B, L, J_min) 0009 % 0010 % B is the wavelet parameter, 0011 % L is the angular band-limit, 0012 % J_min the first wavelet to be used. 0013 % 0014 % S2LET package to perform Wavelet transform on the Sphere. 0015 % Copyright (C) 2012-2015 Boris Leistedt & Jason McEwen 0016 % See LICENSE.txt for license details 0017 0018 [kappa kappa0] = s2let_transform_axisym_tiling(B, L, J_min); 0019 0020 J = s2let_jmax(L, B); 0021 xi = 0:0.01:L-1; 0022 x = 0:L-1; 0023 0024 figure('Position',[100 100 900 450]) 0025 %semilogx(0:L-1, kappa0, 'k', 'LineWidth', 2); 0026 yi = interp1(x,kappa0,xi,'pchip'); 0027 semilogx(xi, yi, 'k', 'LineWidth', 2); 0028 %h = text(2, 1.07, 'k0', 'Color', [0 0 0]); 0029 hold on; 0030 for j = J_min:J 0031 colour = rand(1,3)*0.9; 0032 %plot(0:L-1, kappa(j+1,:), 'LineWidth', 2, 'Color', colour); 0033 yi = interp1(x,kappa(j+1,:),xi,'pchip'); 0034 plot(xi, yi, 'LineWidth', 2, 'Color', colour); 0035 %h = text(B.^j, 1.07, strcat('j',num2str(j+1)), 'Color', colour); 0036 end 0037 %title('Harmonic tiling'); 0038 %xlabel('el'); 0039 axis([0 L -0.05 1.15]); 0040 set(gca,'XTick',2.^[0:(J+2)]); 0041 0042 end