s2let_demo7 - plotting curvelets 1) Plot the scaling function and curvelet functions in real space (s2let_plot_curvelet_tiling). 2) Plot the curvelets on the Sphere (s2let_plot_curvelet_on_sphere). 3) Plot the curvelets parametrically (s2let_plot_curvelet_parametric). wherein curvelets and the scaling functions are generated via the matlab function "s2let_curvelet_tiling(B, L, J_min, <options>)". ----------------------------------------------------------- S2LET package to perform Wavelet Transform on the Sphere. Copyright (C) 2012-2016 Boris Leistedt, Jennifer Chan & Jason McEwen See LICENSE.txt for license details -----------------------------------------------------------
0001 % s2let_demo7 - plotting curvelets 0002 % 0003 % 1) Plot the scaling function and curvelet functions in real space (s2let_plot_curvelet_tiling). 0004 % 2) Plot the curvelets on the Sphere (s2let_plot_curvelet_on_sphere). 0005 % 3) Plot the curvelets parametrically (s2let_plot_curvelet_parametric). 0006 % wherein curvelets and the scaling functions are generated via 0007 % the matlab function "s2let_curvelet_tiling(B, L, J_min, <options>)". 0008 % 0009 % ----------------------------------------------------------- 0010 % S2LET package to perform Wavelet Transform on the Sphere. 0011 % Copyright (C) 2012-2016 Boris Leistedt, Jennifer Chan & Jason McEwen 0012 % See LICENSE.txt for license details 0013 % ----------------------------------------------------------- 0014 0015 clear all; 0016 0017 % --------------- 0018 % Define curvelet parameters: 0019 % --------------- 0020 Spin = 0; % Spin value of curvelet 0021 B = 2; % B = 2 for dyadic sampling 0022 L = 128; % Angular band-limit 0023 J_min = 2; % Minimum scale probed by curvelets 0024 J =s2let_jmax(L, B); % Maximum scale probed by curvelets =ceil(log L/ log B); 0025 0026 0027 % --------------- 0028 % Plot the tiling of scaling function and curvelets in harmonic space: 0029 % And plot the scaling function and curvelet functions in real space: 0030 % --------------- 0031 disp(' - Plot the tiling of curvelets'); 0032 s2let_plot_curvelet_tiling(B, L, J_min,... 0033 'Spin', Spin); 0034 0035 % Define Euler angles for rotating the curvelets functions 0036 % such that the curvelets centred on the North pole of the sphere: 0037 alpha = 0 ; 0038 beta = acos(-Spin/L); 0039 gamma = 0 ; 0040 0041 % --------------- 0042 % Plot the curvelets on the sphere: 0043 % --------------- 0044 disp(' - Plot the curvelets on the sphere'); 0045 s2let_plot_curvelet_on_sphere(alpha, beta, gamma, B, L, J_min, ... 0046 'Spin', Spin); 0047 0048 %{ 0049 disp(' - Plot the curvelet parametrically'); 0050 s2let_plot_curvelet_parametric(alpha, beta, gamma,... 0051 B, L, J_min, ... 0052 'Spin', Spin); 0053 %}