Home > src > main > matlab > flaglet_plot_axisym_wavelet_kernel.m

flaglet_plot_axisym_wavelet_kernel

PURPOSE ^

SYNOPSIS ^

function flaglet_plot_axisym_wavelet_kernel(kappa_ln, R, Rplot, P)

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function flaglet_plot_axisym_wavelet_kernel(kappa_ln, R, Rplot, P)
0002 
0003 sz = size(kappa_ln);
0004 N = sz(1);
0005 L = sz(2);
0006 
0007 flmn = zeros(N,L^2);
0008 for l = 0:L-1
0009     for n = 0:N-1
0010         flmn(n + 1, l^2 + l + 1) = kappa_ln(n+1,l+1);
0011     end
0012 end
0013 
0014 h = (R/(P));
0015 rs = 0:h:R ;
0016 [~, thetas, ~] = flag_sampling(L, N, 1.0);
0017 
0018 ind = find(rs <= Rplot);
0019 nind = length(ind);
0020 
0021 thetas = pi/2 - thetas; 
0022 thetas = [ pi/2 thetas -pi-fliplr(thetas(1:L-1)) pi/2 ];
0023 
0024 f = (flag_synthesis(flmn, 'Nodes', rs, 'reality', true));
0025 f(1:3,:,:) = 0.0;
0026 
0027 [rs, thetas] = ndgrid(rs(ind),thetas);
0028 x = rs .* cos(thetas);
0029 y = rs .* sin(thetas);
0030 
0031 v = zeros(nind,2*L+1);
0032 v(:,1) = f(ind,1,1);
0033 v(:,2:(L+1)) = f(ind,:,L);
0034 v(:,(L+2):(2*L)) = fliplr(f(ind,1:L-1,1));
0035 v(:,2*L+1) = f(ind,1,1);
0036 
0037 %size(x)
0038 %size(y)
0039 %size(v)
0040 
0041 ind_neg = find(v < 0);
0042 v(ind_neg) = v(ind_neg) / abs(min(min(min(v))));
0043  
0044 ind_pos = find(v > 0);
0045 v(ind_pos) = v(ind_pos) / max(max(max(v)));
0046 
0047 %[min(min(min(v))), max(max(max(v)))]
0048 %limit = max( [ -min(min(min(v))) max(max(max(v))) ]) ;
0049 %levels = -limit:(2*limit/nlevels):limit ;
0050 %levels = -1:(2/nlevels):1 ;
0051 
0052 p = polar([0 2*pi], [0 Rplot]);
0053 grid off
0054 hold on
0055 surface(x, y, zeros(size(x)), v, 'EdgeColor', 'none');
0056 %colorbar
0057 ph=findall(gca,'type','text'); 
0058 ps=get(ph,'string'); 
0059 %disp([num2cell(1:numel(ps)).',ps]);
0060 ps(1:numel(ps))={''}; 
0061 %ps(16:19)={''};
0062 ps(4:11)={ 
0063           '90'
0064           '-90'
0065           ''
0066           ''
0067           ''
0068           ''
0069           '180'
0070           '0'
0071      }; 
0072 
0073 set(ph,{'string'},ps); 
0074 delete(findall(ancestor(p,'figure'),'HandleVisibility','off','type','line'));
0075 delete(p)
0076 %axis([-Rplot Rplot -Rplot+0.2 Rplot+0.2])
0077 end

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