809 research outputs found
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
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
Exact reconstruction with directional wavelets on the sphere
A new formalism is derived for the analysis and exact reconstruction of
band-limited signals on the sphere with directional wavelets. It represents an
evolution of the wavelet formalism developed by Antoine & Vandergheynst (1999)
and Wiaux et al. (2005). The translations of the wavelets at any point on the
sphere and their proper rotations are still defined through the continuous
three-dimensional rotations. The dilations of the wavelets are directly defined
in harmonic space through a new kernel dilation, which is a modification of an
existing harmonic dilation. A family of factorized steerable functions with
compact harmonic support which are suitable for this kernel dilation is firstly
identified. A scale discretized wavelet formalism is then derived, relying on
this dilation. The discrete nature of the analysis scales allows the exact
reconstruction of band-limited signals. A corresponding exact multi-resolution
algorithm is finally described and an implementation is tested. The formalism
is of interest notably for the denoising or the deconvolution of signals on the
sphere with a sparse expansion in wavelets. In astrophysics, it finds a
particular application for the identification of localized directional features
in the cosmic microwave background (CMB) data, such as the imprint of
topological defects, in particular cosmic strings, and for their reconstruction
after separation from the other signal components.Comment: 22 pages, 2 figures. Version 2 matches version accepted for
publication in MNRAS. Version 3 (identical to version 2) posted for code
release announcement - "Steerable scale discretised wavelets on the sphere" -
S2DW code available for download at
http://www.mrao.cam.ac.uk/~jdm57/software.htm
3D weak lensing with spin wavelets on the ball
We construct the spin flaglet transform, a wavelet transform to analyze spin
signals in three dimensions. Spin flaglets can probe signal content localized
simultaneously in space and frequency and, moreover, are separable so that
their angular and radial properties can be controlled independently. They are
particularly suited to analyzing of cosmological observations such as the weak
gravitational lensing of galaxies. Such observations have a unique 3D
geometrical setting since they are natively made on the sky, have spin angular
symmetries, and are extended in the radial direction by additional distance or
redshift information. Flaglets are constructed in the harmonic space defined by
the Fourier-Laguerre transform, previously defined for scalar functions and
extended here to signals with spin symmetries. Thanks to various sampling
theorems, both the Fourier-Laguerre and flaglet transforms are theoretically
exact when applied to bandlimited signals. In other words, in numerical
computations the only loss of information is due to the finite representation
of floating point numbers. We develop a 3D framework relating the weak lensing
power spectrum to covariances of flaglet coefficients. We suggest that the
resulting novel flaglet weak lensing estimator offers a powerful alternative to
common 2D and 3D approaches to accurately capture cosmological information.
While standard weak lensing analyses focus on either real or harmonic space
representations (i.e., correlation functions or Fourier-Bessel power spectra,
respectively), a wavelet approach inherits the advantages of both techniques,
where both complicated sky coverage and uncertainties associated with the
physical modeling of small scales can be handled effectively. Our codes to
compute the Fourier-Laguerre and flaglet transforms are made publicly
available.Comment: 24 pages, 4 figures, version accepted for publication in PR
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
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
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
Wavelet reconstruction of E and B modes for CMB polarisation and cosmic shear analyses
We present new methods for mapping the curl-free (E-mode) and divergence-free (B-mode) components of spin 2 signals using spin directional wavelets. Our methods are equally applicable to measurements of the polarisation of the cosmic microwave background (CMB) and the shear of galaxy shapes due to weak gravitational lensing. We derive pseudo and pure wavelet estimators, where E-B mixing arising due to incomplete sky coverage is suppressed in wavelet space using scale- and orientation-dependent masking and weighting schemes. In the case of the pure estimator, ambiguous modes (which have vanishing curl and divergence
simultaneously on the incomplete sky) are also cancelled. On simulations, we demonstrate the improvement (i.e., reduction in leakage) provided by our wavelet space estimators over standard harmonic space approaches. Our new methods can be directly interfaced in a coherent
and computationally-efficient manner with component separation or feature extraction techniques that also exploit wavelets
- …