This page was created by the IDL library routine
mk_html_help. For more information on
this routine, refer to the IDL Online Help Navigator
or type:
? mk_html_help
at the IDL command line prompt.
Last modified: Tue Sep 22 12:42:56 2015.
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_hpx_wav_analysis
PURPOSE:
Compute wavelet transform of a Healpix map
CALLING SEQUENCE:
f_wav = s2let_axisym_hpx_wav_analysis(f, B, L, J_min)
INPUTS:
f - input healpix nside (number of pixels is 12*nside^2)
B - Wavelet parameter
L - Band-limit to be used for the spherical harmonic transforms
J_min - First wavelet scale to be used
wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines)
OUTPUTS:
f_wav - Struc containing the wavelets
EXAMPLE:
see s2let_hpx_demo
read_fits_map, file, f
f_wav = s2let_axisym_hpx_wav_analysis(f, B, L, J_min)
f_rec = s2let_axisym_hpx_wav_synthesis(f_wav)
J_max = s2let_j_max(L, b)
mollview, f_rec, title='Reconstructed map'
mollview, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
mollview, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
Nside of the input map is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_hpx_wav_analysis.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_hpx_wav_synthesis
PURPOSE:
Reconstruct a healpix map from its wavelet transform
CALLING SEQUENCE:
f = s2let_axisym_hpx_wav_synthesis(f_wav)
INPUTS:
f_wav - Struc containing the wavelets maps
OUTPUTS:
f - The reconstructed Healpix map
EXAMPLE:
see s2let_hpx_demo
read_fits_map, file, f
f_wav = s2let_axisym_hpx_wav_analysis(f, B, L, J_min)
f_rec = s2let_axisym_hpx_wav_synthesis(f_wav)
J_max = s2let_j_max(L, b)
mollview, f_rec, title='Reconstructed map'
mollview, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
mollview, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
B, L and J_min of the input wavelets are automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_hpx_wav_synthesis.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_analysis
PURPOSE:
Compute full resolution exact wavelet transform of a complex MW map
CALLING SEQUENCE:
f_wav = s2let_axisym_mw_wav_analysis(f, B, J_min)
INPUTS:
f - input MW map (number of pixels is L*(2*L-1))
B - Wavelet parameter
J_min - First wavelet scale to be used
wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines)
OUTPUTS:
f_wav - Struc containing the wavelets
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The resolution/bandlimit L of the input map is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_analysis.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_analysis_multires
PURPOSE:
Compute multiresolution exact wavelet transform of a complex MW map
CALLING SEQUENCE:
f_wav = s2let_axisym_mw_wav_analysis_multires(f, B, J_min)
INPUTS:
f - input MW map (number of pixels is L*(2*L-1))
B - Wavelet parameter
J_min - First wavelet scale to be used
wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines)
OUTPUTS:
f_wav - Struc containing the wavelets
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The resolution/bandlimit L of the input map is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_analysis_multires.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_analysis_multires_real
PURPOSE:
Compute multiresolution exact wavelet transform of a real MW map
CALLING SEQUENCE:
f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
INPUTS:
f - input MW map (number of pixels is L*(2*L-1))
B - Wavelet parameter
J_min - First wavelet scale to be used
wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines)
OUTPUTS:
f_wav - Struc containing the wavelets
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The resolution/bandlimit L of the input map is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_analysis_multires_real.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_analysis_real
PURPOSE:
Compute full resolution exact wavelet transform of a real MW map
CALLING SEQUENCE:
f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
INPUTS:
f - input MW map (number of pixels is L*(2*L-1))
B - Wavelet parameter
J_min - First wavelet scale to be used
wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines)
OUTPUTS:
f_wav - Struc containing the wavelets
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The resolution/bandlimit L of the input map is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_analysis_real.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_synthesis
PURPOSE:
Reconstruct a complex MW map from its full resolution exact wavelet transform
CALLING SEQUENCE:
f = s2let_axisym_mw_wav_synthesis(f_wav)
INPUTS:
f_wav - input exact wavelet transfrom
OUTPUTS:
f - Reconstructed MW map
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The wavelet parameters B/J_min are automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_synthesis.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_synthesis_multires
PURPOSE:
Reconstruct a complex MW map from its multiresolution exact wavelet transform
CALLING SEQUENCE:
f = s2let_axisym_mw_wav_synthesis_multires(f_wav)
INPUTS:
f_wav - input exact wavelet transfrom
OUTPUTS:
f - Reconstructed MW map
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The wavelet parameters B/J_min are automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_synthesis_multires.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_synthesis_multires_real
PURPOSE:
Reconstruct a complex MW map from its full resolution exact wavelet transform
CALLING SEQUENCE:
f = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
INPUTS:
f_wav - input exact wavelet transfrom
OUTPUTS:
f - Reconstructed MW map
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The wavelet parameters B/J_min are automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_synthesis_multires_real.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_axisym_mw_wav_synthesis_real
PURPOSE:
Reconstruct a real MW map from its full resolution exact wavelet transform
CALLING SEQUENCE:
f = s2let_axisym_mw_wav_synthesis_real(f_wav)
INPUTS:
f_wav - input exact wavelet transfrom
OUTPUTS:
f - Reconstructed MW map
EXAMPLE:
see s2let_demo and s2let_test
f = s2let_read_mw_real_map(file)
if multires eq 0 then f_wav = s2let_axisym_mw_wav_analysis_real(f, B, J_min)
else f_wav = s2let_axisym_mw_wav_analysis_multires_real(f, B, J_min)
if multires eq 0 then f_rec = s2let_axisym_mw_wav_synthesis_real(f_wav)
else f_rec = s2let_axisym_mw_wav_synthesis_multires_real(f_wav)
L = s2let_get_mw_bandlimit(f)
J_max = s2let_j_max(L, B)
s2let_plot_mollweide, f_rec, title='Band-limited map'
s2let_plot_mollweide, f_wav.scal, title='Scaling map'
for j=0, J_max-J_min do begin
s2let_plot_mollweide, f_wav.(j), title='Wavelet map '
+strtrim(j+1,2)+' on '+strtrim(J_max-J_min+1,2)
endfor
COMMENTS:
The wavelet parameters B/J_min are automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_axisym_mw_wav_synthesis_real.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_demo1 PURPOSE: Demo 1 : test the exactness of all MW spherical harmonics and wavelet transforms for a tomography map of the earth and plot the wavelet maps OPTIONAL KEYWORDS: B - The wavelet parameter for the test J_min - The first wavelet scale to be used for the transform multires - multiresolution flag (1 if yes, else 0) wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines) DEFAULT VALUES: B=3, J_min=2, multires=0, wavtype=1
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_demo1.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_demo2 PURPOSE: Demo : read Healpix map, convert it to MW map, then test the exactness of all MW spherical harmonics and wavelet transforms and plot the wavelet maps OPTIONAL KEYWORDS: L - The bandlimit for the spherical harmonic transforms B - The wavelet parameter for the test J_min - The first wavelet scale to be used for the transform multires - multiresolution flag (1 if yes, else 0) wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines) DEFAULT VALUES: L=128, B=3, J_min=2, multires=0, wavtype=1
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_demo2.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_dylib_exists PURPOSE: Check if the s2let dynamic library exists somewhere in the path CALLING SEQUENCE: status = s2let_dylib_exists() OUTPUT 1 if found, 0 otherwise
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_dylib_exists.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_get_dylib PURPOSE: Get the location/name of the s2let dynamic library CALLING SEQUENCE: location = s2let_get_dylib()
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_get_dylib.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_get_mw_bandlimit PURPOSE: Detect the resolution/bandlimit L of an input MW map CALLING SEQUENCE: L = s2let_get_mw_bandlimit(f) INPUT: f - array/map of size L*(2*L-1) OUTPUT L - the resolution/bandlimit
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_get_mw_bandlimit.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_get_wav_bandlimit PURPOSE: Compute the bandlimit of the j-th wavelet constructed with parameter B CALLING SEQUENCE: bl = s2let_get_wav_bandlimit(B, j) OUTPUT The bandlimit of the j-th wavelet constructed with parameter B
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_get_wav_bandlimit.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_healpix2mw
PURPOSE:
Convert a real Healpix map to an MW map
through the spherical harmonic transform
CALLING SEQUENCE:
mwmap = s2let_healpix2mw(hpxmap, lmax=lmax)
INPUT:
hpxmap - The Healpix map (npix = 12*nside*nside)
lmax - The bandlimit of the SHA transform, thus the
resolution of the MW map
OUTPUT
mwmap - The MW map (npix = lmax*(2*lmax-1))
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_healpix2mw.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_hpx_alm2map_real
PURPOSE:
Reconstruct a real Healpix map from its spherical harmonic
transform (nb: not an exact transform due to healpix)
CALLING SEQUENCE:
f = s2let_hpx_alm2map_real(alm, nside)
INPUT:
alm - a complex array containing the spherical harmonic transform
alm_{el em} can be accessed at alm(ind) with ind = el*el + el + em
nside - resolution for the output map
OUTPUT
f - The Healpix map (npix = 12*nside*nside)
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_hpx_alm2map_real.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_hpx_demo PURPOSE: Demo : test all HPX spherical harmonics and wavelet transforms (not exact due to healpix) for a simulated CMB map and plot the wavelet maps OPTIONAL KEYWORDS: L - The bandlimit for the spherical harmonic transforms B - The wavelet parameter for the test J_min - The first wavelet scale to be used for the transform wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines) DEFAULT VALUES: L=192, B=7, J_min=2, wavtype=1
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_hpx_demo.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_j_max PURPOSE: Compute the maximum wavelet to be used, given a bandlimit L and a wavelet parameter B CALLING SEQUENCE: J_max = s2let_j_max(L, B) OUTPUT Compute the maximum wavelet to be used, given a bandlimit L and a wavelet parameter B;
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_j_max.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_mw_alm2map
PURPOSE:
Reconstruct a complex MW map from its exact spherical harmonic transform
CALLING SEQUENCE:
f = s2let_mw_alm2map(alm)
INPUT:
alm - a complex array containing the spherical harmonic transform
alm_{el em} can be accessed at alm(ind) with ind = el*el + el + em
OUTPUT
f - The MW map (npix = L*(2*L-1), L is detected)
COMMENT:
The resolution/bandlimit L is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_alm2map.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_mw_alm2map_real
PURPOSE:
Reconstruct a real MW map from its exact spherical harmonic transform
CALLING SEQUENCE:
f = s2let_mw_alm2map_real(alm)
INPUT:
alm - a complex array containing the spherical harmonic transform
alm_{el em} can be accessed at alm(ind) with ind = el*el + el + em
OUTPUT
f - The MW map (npix = L*(2*L-1), L is detected)
COMMENT:
The resolution/bandlimit L is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_alm2map_real.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_mw_map2alm_real
PURPOSE:
Compute the spherical harmonic transform of a real Healpix map
(nb: not an exact transform due to healpix)
CALLING SEQUENCE:
alm = s2let_hpx_map2alm_real(f, lmax)
INPUTS
f - The input healpix map (npix = 12*nside*nside, nside is detected)
OUTPUT:
alm - a complex array containing the spherical harmonic transform
alm_{el em} can be accessed at alm(ind) with ind = el*el + el + em
COMMENT:
The resolution nside is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_hpx_map2alm_real.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_mw_map2alm_real
PURPOSE:
Compute exact spherical harmonic transform of a complex MW map
CALLING SEQUENCE:
alm = s2let_mw_map2alm(f)
INPUTS
f - The input MW map (npix = L*(2*L-1), L is detected)
OUTPUT:
alm - a complex array containing the spherical harmonic transform
alm_{el em} can be accessed at alm(ind) with ind = el*el + el + em
COMMENT:
The resolution/bandlimit L is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_map2alm.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_mw_map2alm_real
PURPOSE:
Compute exact spherical harmonic transform of a real MW map
CALLING SEQUENCE:
alm = s2let_mw_map2alm_real(f)
INPUTS
f - The input MW map (npix = L*(2*L-1), L is detected)
OUTPUT:
alm - a complex array containing the spherical harmonic transform
alm_{el em} can be accessed at alm(ind) with ind = el*el + el + em
COMMENT:
The resolution/bandlimit L is automatically detected
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_map2alm_real.pro)
S2LET package - Copyright (C) 2012
Boris Leistedt & Jason McEwen
NAME:
s2let_mw_pixel_edges
PURPOSE:
Get (theta, phi) coordinates for the corners of the i-th pixel in
the MW sampling
CALLING SEQUENCE:
arr = s2let_mw_pixel_edges(L, i)
INPUTS
L - resolution/bandlimit of the MW map
i - pixel index
OUTPUT:
arr - an array containing four numbers: theta1, theta2, phi1 and phi2
which are the locations of the corners of the pixel
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_pixel_edges.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_mw_plot_mollweide PURPOSE: Plot a real MW map (from data of FITS file) under Mollweide projection CALLING SEQUENCE: s2let_mw_plot_mollweide, mapfile, title='My map' INPUTS maporfile - filename for the FITS or map read with s2let_read_mw_real_map
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_plot_mollweide.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_mw_sampling PURPOSE: Compute the coordinates of the nodes of the MW sampling for a given bandlimit/resolution CALLING SEQUENCE: s2let_mw_sampling, L, thetas, phis INPUT: L - resolution/bandlimit OUTPUT thetas - The theta coordinates of the nodes phis - The phi coordinates of the nodes
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_sampling.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_mw_write_real_map PURPOSE: Write a real MW map to a FITS file CALLING SEQUENCE: s2let_mw_write_real_map, map, file INPUTS: map - input MW map (npix=L*(2*L-1), L is detected) file - filename for the FITS
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_write_real_map.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_read_mw_real_map PURPOSE: Read a real MW map from a FITS file CALLING SEQUENCE: f = s2let_mw_read_real_map(file) INPUTS file - filename for the FITS OUTPUT: map - MW map (npix=L*(2*L-1), L is detected)
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_mw_read_real_map.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_test PURPOSE: Test the exactness of all MW spherical harmonics and wavelet transform for real and complex signals on the sphere OPTIONAL KEYWORDS: L - The bandlimit/resolution for the test B - The wavelet parameter for the test J_min - The first wavelet scale to be used for the transform multires - multiresolution flag (1 if yes, else 0) wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines) DEFAULT VALUES: L=64, B=3, J_min=2, multires=0, wavtype=1
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_test.pro)
S2LET package - Copyright (C) 2012 Boris Leistedt & Jason McEwen NAME: s2let_valid_wav_parameters PURPOSE: Check and test a set of wavelet parameters CALLING SEQUENCE: status = s2let_valid_wav_parameters(B, L, J_min) INPUTS: B - Wavelet parameter L - Band-limit to be used for the spherical harmonic transforms J_min - First wavelet scale to be used wavtype - Wavelet type (1: scale-discretised, 2:needlets, 3: cubic splines)
(See /Users/bl/Dropbox/Wavelets/s2let/src/main/idl/s2let_valid_wav_parameters.pro)