Home > src > main > matlab > s2let_demo2.m

s2let_demo2

PURPOSE ^

s2let_demo2

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 s2let_demo2
 Compute and plot the harmonic tiling and the wavelet kernels.

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % s2let_demo2
0002 % Compute and plot the harmonic tiling and the wavelet kernels.
0003 %
0004 % S2LET package to perform Wavelets on the Sphere.
0005 % Copyright (C) 2012-2015  Boris Leistedt & Jason McEwen
0006 % See LICENSE.txt for license details
0007 
0008 B = 3;
0009 zoomfactor = 1.4;
0010 J_min = 2;
0011 L = 128;
0012 J = s2let_jmax(L, B);
0013 type = 'colour';
0014 lighting = true;
0015 
0016 ns = ceil(sqrt(2+J-J_min+1)) ;
0017 nx = ns - 1 + rem(2+J-J_min + 1, ns) ;
0018 ny = ns;
0019 
0020 [kappa kappa0] = s2let_transform_axisym_tiling(B, L, J_min);
0021 
0022 s2let_plot_axisym_tiling(B, L, J_min);
0023 
0024 figure('Position',[100 100 1100 700])
0025 
0026 h=subplot(nx, ny, 1);
0027 flm = zeros(L^2,1);
0028 for l = 0:L-1
0029     flm(l^2+l+1,1) = kappa0(l+1);
0030 end
0031 f = ssht_inverse(flm, L, 'Reality', true);
0032 ssht_plot_sphere(f, L, 'Type', type, 'Lighting', lighting);
0033 zoom(zoomfactor)
0034 v = caxis;
0035 temp = max(abs(v));
0036 caxis([-temp temp])
0037 %title('k0')
0038 locate = get(h,'title');
0039 pos = get(locate,'position');
0040 pos(1,2) = pos(1,2)+0.7;
0041 pos(1,1) = pos(1,1)-0.7;
0042 set(locate,'pos',pos);
0043 
0044 for j = J_min:5
0045    h=subplot(nx, ny, j-J_min+2);
0046    flm = zeros(L^2,1);
0047     for l = 0:L-1
0048         flm(l^2+l+1,1) = kappa(j+1,l+1);
0049     end
0050    f = ssht_inverse(flm, L, 'Reality', true);
0051    ssht_plot_sphere(f, L, 'Type', type, 'Lighting', lighting);
0052    v = caxis;
0053    temp = max(abs(v));
0054    caxis([-temp temp])
0055    colormap(jet)
0056    %title(['j',int2str(j)-J_min+1])
0057    locate = get(h,'title');
0058    pos = get(locate,'position');
0059    pos(1,2) = pos(1,2)+0.7;
0060    pos(1,1) = pos(1,1)-0.7;
0061    set(locate,'pos',pos);
0062    zoom(zoomfactor)
0063 end
0064 
0065 
0066

Generated on Fri 11-Nov-2016 11:50:36 by m2html © 2005