Home > src > main > matlab > flaglet_demo3.m

flaglet_demo3

PURPOSE ^

flaglet_demo3

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 flaglet_demo3

 Plot wavelet decomposition of LSS dataset

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

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % flaglet_demo3
0002 %
0003 % Plot wavelet decomposition of LSS dataset
0004 %
0005 % flaglet package to perform Wavelets transform on the Solid Sphere.
0006 % Copyright (C) 2012  Boris Leistedt & Jason McEwen
0007 % See LICENSE.txt for license details
0008 
0009 L = 128 ;
0010 P = 128 ;
0011 Lplot = 256;
0012 Pplot = 256;
0013 
0014 B_l = 2;
0015 B_p = 2;
0016 J_min_l = 6;
0017 J_min_p = 6;
0018 downsample = false;
0019 
0020 J_l = ceil(log(L)/log(B_l))
0021 J_p = ceil(log(P)/log(B_p))
0022 
0023 % Load LSS data
0024 load('lss_flmp_312')
0025 R = 420;
0026 
0027 f_lmp = zeros(P, L^2);
0028 f_lmp(:,:) = f_lmp_full(1:P, 1:L^2);
0029 
0030 disp('Flag_synthesis...')
0031 f = flag_synthesis(f_lmp, 'R', R, 'reality', true);
0032 
0033 % subplot(1,2,1)
0034 % flaglet_plot_f( f_ini )
0035 % caxis([minbound maxbound])
0036 % subplot(1,2,2)
0037 % flaglet_plot_f( f )
0038 % caxis([minbound maxbound])
0039  
0040 disp('Wavelet analysis...')
0041 [f_wav, f_scal] = flaglet_axisym_analysis(f, 'B_l', B_l, 'B_p', B_p, 'reality', true, 'downsample', downsample, 'J_min_l', J_min_l, 'J_min_p', J_min_p);
0042 
0043 % nodes = slag_sampling(256, R);
0044 % f_lmp_zeropadded = zeros(P, 256^2);
0045 % f_lmp_zeropadded(1:P, 1:L^2) = f_lmp(:,:);
0046 % disp('flag synthesis')
0047 % f_oversampled = flag_synthesis(f_lmp_zeropadded, 'Nodes', nodes, 'reality', true);
0048 
0049 
0050 disp('Plot...')
0051 figure('Position',[100 100 600 900]) 
0052 
0053 minbound = 10;
0054 maxbound = 100;
0055 layer = 230;
0056 
0057 J_l = [ 1 2 ];
0058 J_n = [ 1 2 ];
0059 nx = 2;
0060 ny = 3;
0061 zoomfactor = 2.1;
0062 
0063 ind = 1;
0064 subplot(ny, nx, ind);
0065 flaglet_plot_f( f, 'L', Lplot, 'P', Pplot, 'layer', layer )
0066 zoom(zoomfactor)
0067 colormap(flipud(gray))
0068 caxis([minbound maxbound])
0069 %title('Data')
0070 
0071 ind = 2;
0072 subplot(ny, nx, ind);
0073 flaglet_plot_f( f_scal, 'L', Lplot, 'P', Pplot, 'layer', layer )
0074 zoom(zoomfactor)
0075 colormap(flipud(gray))
0076 caxis([minbound maxbound])
0077 %title('Scaling part')
0078 
0079 for jl = J_l
0080     for jn = J_n
0081         ind = ind + 1; 
0082         subplot(ny, nx, ind);
0083         flaglet_plot_f( f_wav{jl, jn}, 'L', Lplot, 'P', Pplot, 'layer', layer )
0084         zoom(zoomfactor)
0085         colormap(flipud(gray))
0086         %title(['Wavelet : jl=',int2str(jl),' jn=',int2str(jn)])
0087         caxis([minbound maxbound])
0088     end
0089 end

Generated on Mon 24-Sep-2012 12:26:33 by m2html © 2005