Issue |
A&A
Volume 558, October 2013
|
|
---|---|---|
Article Number | A128 | |
Number of page(s) | 9 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/201220729 | |
Published online | 18 October 2013 |
S2LET: A code to perform fast wavelet analysis on the sphere
1 Department of Physics and
Astronomy, University College London, London WC1E
6BT, UK
e-mail: boris.leistedt.11@ucl.ac.uk
2 Mullard Space Science Laboratory
(MSSL), University College London, Surrey
RH5 6NT,
UK
e-mail: jason.mcewen@ucl.ac.uk
3 Institute of Electrical Engineering,
Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015
Lausanne,
Switzerland
e-mail: pierre.vandergheynst@epfl.ch
4 Department of Radiology and Medical
Informatics, University of Geneva (UniGE), 1211
Geneva,
Switzerland
5 Department of Medical Radiology,
Lausanne University Hospital (CHUV), 1011
Lausanne,
Switzerland
e-mail: yves.wiaux@epfl.ch
6 Institute of Sensors, Signals
& Systems, Heriot Watt University, Edinburgh
EH14 4AS,
UK
Received:
13
November
2012
Accepted:
27
August
2013
We describe S2LET, a fast and robust implementation of the scale-discretised wavelet transform on the sphere. Wavelets are constructed through a tiling of the harmonic line and can be used to probe spatially localised, scale-dependent features of signals on the sphere. The reconstruction of a signal from its wavelets coefficients is made exact here through the use of a sampling theorem on the sphere. Moreover, a multiresolution algorithm is presented to capture all information of each wavelet scale in the minimal number of samples on the sphere. In addition S2LET supports the HEALPix pixelisation scheme, in which case the transform is not exact but nevertheless achieves good numerical accuracy. The core routines of S2LET are written in C and have interfaces in Matlab, IDL and Java. Real signals can be written to and read from FITS files and plotted as Mollweide projections. The S2LET code is made publicly available, is extensively documented, and ships with several examples in the four languages supported. At present the code is restricted to axisymmetric wavelets but will be extended to directional, steerable wavelets in a future release.
Key words: techniques: high angular resolution / methods: data analysis / techniques: image processing
© ESO, 2013
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.