3,049 research outputs found
Fast directional continuous spherical wavelet transform algorithms
We describe the construction of a spherical wavelet analysis through the
inverse stereographic projection of the Euclidean planar wavelet framework,
introduced originally by Antoine and Vandergheynst and developed further by
Wiaux et al. Fast algorithms for performing the directional continuous wavelet
analysis on the unit sphere are presented. The fast directional algorithm,
based on the fast spherical convolution algorithm developed by Wandelt and
Gorski, provides a saving of O(sqrt(Npix)) over a direct quadrature
implementation for Npix pixels on the sphere, and allows one to perform a
directional spherical wavelet analysis of a 10^6 pixel map on a personal
computer.Comment: 10 pages, 3 figures, replaced to match version accepted by IEEE
Trans. Sig. Pro
On the computation of directional scale-discretized wavelet transforms on the sphere
We review scale-discretized wavelets on the sphere, which are directional and
allow one to probe oriented structure in data defined on the sphere.
Furthermore, scale-discretized wavelets allow in practice the exact synthesis
of a signal from its wavelet coefficients. We present exact and efficient
algorithms to compute the scale-discretized wavelet transform of band-limited
signals on the sphere. These algorithms are implemented in the publicly
available S2DW code. We release a new version of S2DW that is parallelized and
contains additional code optimizations. Note that scale-discretized wavelets
can be viewed as a directional generalization of needlets. Finally, we outline
future improvements to the algorithms presented, which can be achieved by
exploiting a new sampling theorem on the sphere developed recently by some of
the authors.Comment: 13 pages, 3 figures, Proceedings of Wavelets and Sparsity XV, SPIE
Optics and Photonics 2013, Code is publicly available at http://www.s2dw.org
Sparse image reconstruction on the sphere: analysis and synthesis
We develop techniques to solve ill-posed inverse problems on the sphere by
sparse regularisation, exploiting sparsity in both axisymmetric and directional
scale-discretised wavelet space. Denoising, inpainting, and deconvolution
problems, and combinations thereof, are considered as examples. Inverse
problems are solved in both the analysis and synthesis settings, with a number
of different sampling schemes. The most effective approach is that with the
most restricted solution-space, which depends on the interplay between the
adopted sampling scheme, the selection of the analysis/synthesis problem, and
any weighting of the l1 norm appearing in the regularisation problem. More
efficient sampling schemes on the sphere improve reconstruction fidelity by
restricting the solution-space and also by improving sparsity in wavelet space.
We apply the technique to denoise Planck 353 GHz observations, improving the
ability to extract the structure of Galactic dust emission, which is important
for studying Galactic magnetism.Comment: 11 pages, 6 Figure
Complex data processing: fast wavelet analysis on the sphere
In the general context of complex data processing, this paper reviews a
recent practical approach to the continuous wavelet formalism on the sphere.
This formalism notably yields a correspondence principle which relates wavelets
on the plane and on the sphere. Two fast algorithms are also presented for the
analysis of signals on the sphere with steerable wavelets.Comment: 20 pages, 5 figures, JFAA style, paper invited to J. Fourier Anal.
and Appli
A novel sampling theorem on the rotation group
We develop a novel sampling theorem for functions defined on the
three-dimensional rotation group SO(3) by connecting the rotation group to the
three-torus through a periodic extension. Our sampling theorem requires
samples to capture all of the information content of a signal band-limited at
, reducing the number of required samples by a factor of two compared to
other equiangular sampling theorems. We present fast algorithms to compute the
associated Fourier transform on the rotation group, the so-called Wigner
transform, which scale as , compared to the naive scaling of .
For the common case of a low directional band-limit , complexity is reduced
to . Our fast algorithms will be of direct use in speeding up the
computation of directional wavelet transforms on the sphere. We make our SO3
code implementing these algorithms publicly available.Comment: 5 pages, 2 figures, minor changes to match version accepted for
publication. Code available at http://www.sothree.or
Scale-discretised ridgelet transform on the sphere
We revisit the spherical Radon transform, also called the Funk-Radon
transform, viewing it as an axisymmetric convolution on the sphere. Viewing the
spherical Radon transform in this manner leads to a straightforward derivation
of its spherical harmonic representation, from which we show the spherical
Radon transform can be inverted exactly for signals exhibiting antipodal
symmetry. We then construct a spherical ridgelet transform by composing the
spherical Radon and scale-discretised wavelet transforms on the sphere. The
resulting spherical ridgelet transform also admits exact inversion for
antipodal signals. The restriction to antipodal signals is expected since the
spherical Radon and ridgelet transforms themselves result in signals that
exhibit antipodal symmetry. Our ridgelet transform is defined natively on the
sphere, probes signal content globally along great circles, does not exhibit
blocking artefacts, supports spin signals and exhibits an exact and explicit
inverse transform. No alternative ridgelet construction on the sphere satisfies
all of these properties. Our implementation of the spherical Radon and ridgelet
transforms is made publicly available. Finally, we illustrate the effectiveness
of spherical ridgelets for diffusion magnetic resonance imaging of white matter
fibers in the brain.Comment: 5 pages, 4 figures, matches version accepted by EUSIPCO, code
available at http://www.s2let.or
S2LET: A code to perform fast wavelet analysis on the sphere
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-depended
features of signals on the sphere. The scale-discretised wavelet transform was
developed previously and reduces to the needlet transform in the axisymmetric
case. 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.Comment: 8 pages, 6 figures, version accepted for publication in A&A. Code is
publicly available from http://www.s2let.or
Localisation of directional scale-discretised wavelets on the sphere
Scale-discretised wavelets yield a directional wavelet framework on the
sphere where a signal can be probed not only in scale and position but also in
orientation. Furthermore, a signal can be synthesised from its wavelet
coefficients exactly, in theory and practice (to machine precision).
Scale-discretised wavelets are closely related to spherical needlets (both were
developed independently at about the same time) but relax the axisymmetric
property of needlets so that directional signal content can be probed. Needlets
have been shown to satisfy important quasi-exponential localisation and
asymptotic uncorrelation properties. We show that these properties also hold
for directional scale-discretised wavelets on the sphere and derive similar
localisation and uncorrelation bounds in both the scalar and spin settings.
Scale-discretised wavelets can thus be considered as directional needlets.Comment: 28 pages, 8 figures, minor changes to match version accepted for
publication by ACH
- …