Home > src > main > matlab > flaglet_demo4.m

flaglet_demo4

PURPOSE ^

flaglet_demo4

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 flaglet_demo4

 Plot wavelet decomposition of geophysics data

 B3LET 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_demo4
0002 %
0003 % Plot wavelet decomposition of geophysics data
0004 %
0005 % B3LET 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 = 50 ;
0010 P = 50 ;
0011 R = 3000;
0012 
0013 B_l = 4;
0014 B_p = 4;
0015 J_min_l = 2;
0016 J_min_p = 2;
0017 downsample = true;
0018 
0019 J_l = ceil(log(L)/log(B_l))
0020 J_p = ceil(log(P)/log(B_p))
0021 
0022 f = gen_geomodel(L, P);
0023 
0024 disp('Flag_analysis...')
0025 f_lmp = flag_analysis(f, 'R', R, 'reality', true);
0026 
0027 % subplot(1,2,1)
0028 % flaglet_plot_f( f_ini )
0029 % caxis([minbound maxbound])
0030 % subplot(1,2,2)
0031 % flaglet_plot_f( f )
0032 % caxis([minbound maxbound])
0033  
0034 disp('Wavelet analysis...')
0035 [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);
0036 
0037 
0038 
0039 disp('Plot...')
0040 figure('Position',[100 100 600 900]) 
0041 
0042 % Layout parameters
0043 minbound = -1.0; 
0044 maxbound = 1.0;
0045 Lplot = 192;
0046 Pplot = 192;
0047 layer = 173;
0048 zoomfactor = 2.2;
0049 
0050 J_l = [ 1 2 ];
0051 J_n = [ 1 2 ];
0052 nx = 2;
0053 ny = 3;
0054 
0055 
0056 ind = 1;
0057 subplot(ny, nx, ind);
0058 flaglet_plot_f( f, 'L', Lplot, 'P', Pplot, 'layer', layer )
0059 zoom(zoomfactor)
0060 colormap(flipud(gray))
0061 caxis([-1.4 1.4])
0062 title('Data')
0063 
0064 ind = 2;
0065 subplot(ny, nx, ind);
0066 flaglet_plot_f( f_scal, 'L', Lplot, 'P', Pplot, 'layer', layer )
0067 zoom(zoomfactor)
0068 colormap(flipud(gray))
0069 caxis([-1.2 1.2])
0070 title('Scaling part')
0071 
0072 for jl = J_l
0073     for jn = J_n
0074         ind = ind + 1; 
0075         subplot(ny, nx, ind);
0076         flaglet_plot_f( f_wav{jl, jn}, 'L', Lplot, 'P', Pplot, 'layer', layer )
0077         zoom(zoomfactor)
0078         %colormap(flipud(gray))
0079         title(['Wavelet : jl=',int2str(jl),' jn=',int2str(jn)])
0080         caxis([minbound maxbound])
0081     end
0082 end

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