We discuss a new scale-discretised directional wavelet transform to analyse spin signals defined on the sphere, in particular the polarisation of the cosmic microwave background (CMB).

A new spin wavelet transform on the sphere is proposed to analyse the polarisation of the cosmic microwave background (CMB), a spin $\pm 2$ signal observed on the celestial sphere. The scalar directional scale-discretised wavelet transform on the sphere is extended to analyse signals of arbitrary spin. The resulting spin scale-discretised wavelet transform probes the directional intensity of spin signals. A procedure is presented using this new spin wavelet transform to recover E- and B-mode signals from partial-sky observations of CMB polarisation.

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.

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.

We study the impact of sampling theorems on the fidelity of sparse image reconstruction on the sphere. We discuss how a reduction in the number of samples required to represent all information content of a band-limited signal acts to improve the fidelity of sparse image reconstruction, through both the dimensionality and sparsity of signals. To demonstrate this result we consider a simple inpainting problem on the sphere and consider images sparse in the magnitude of their gradient. We develop a framework for total variation (TV) inpainting on the sphere, including fast methods to render the inpainting problem computationally feasible at high-resolution. Recently a new sampling theorem on the sphere was developed, reducing the required number of samples by a factor of two for equiangular sampling schemes. Through numerical simulations we verify the enhanced fidelity of sparse image reconstruction due to the more efficient sampling of the sphere provided by the new sampling theorem.

A sampling theorem on the sphere has been developed recently, requiring half as many samples as alternative equiangular sampling theorems on the sphere. A reduction by a factor of two in the number of samples required to represent a band-limited signal on the sphere exactly has important implications for compressed sensing, both in terms of the dimensionality and sparsity of signals. We illustrate the impact of this property with an inpainting problem on the sphere, where we show the superior reconstruction performance when adopting the new sampling theorem compared to the alternative.

We discuss a novel sampling theorem on the sphere developed by McEwen & Wiaux recently through an association between the sphere and the torus. To represent a band-limited signal exactly, this new sampling theorem requires less than half the number of samples of other equiangular sampling theorems on the sphere, such as the canonical Driscoll & Healy sampling theorem. A reduction in the number of samples required to represent a band-limited signal on the sphere has important implications for compressive sensing, both in terms of the dimensionality and sparsity of signals. We illustrate the impact of this property with an inpainting problem on the sphere, where we show superior reconstruction performance when adopting the new sampling theorem.

We propose an algorithm for the reconstruction of the signal induced by cosmic strings in the cosmic microwave background (CMB), from radio-interferometric data at arcminute resolution. Radio interferometry provides incomplete and noisy Fourier measurements of the string signal, which exhibits sparse or compressible magnitude of the gradient due to the Kaiser-Stebbins (KS) effect. In this context the versatile framework of compressed sensing naturally applies for solving the corresponding inverse problem. Our algorithm notably takes advantage of a model of the prior statistical distribution of the signal fitted on the basis of realistic simulations. Enhanced performance relative to the standard CLEAN algorithm is demonstrated by simulated observations under noise conditions including primary and secondary CMB anisotropies.

We consider the probe of astrophysical signals through radio interferometers with small field of view and baselines with non-negligible and constant component in the pointing direction. In this context, the visibilities measured essentially identify with a noisy and incomplete Fourier coverage of the product of the planar signals with a linear chirp modulation. In light of the recent theory of compressed sensing and in the perspective of defining the best possible imaging techniques for sparse signals, we analyze the related spread spectrum phenomenon and suggest its universality relative to the sparsity dictionary. Our results rely both on theoretical considerations related to the mutual coherence between the sparsity and sensing dictionaries, as well as on numerical simulations.

An algorithm is proposed for denoising the signal induced by cosmic strings in the cosmic microwave background (CMB). A Bayesian approach is taken, based on modeling the string signal in the wavelet domain with generalized Gaussian distributions. Good performance of the algorithm is demonstrated by simulated experiments at arcminute resolution under noise conditions including primary and secondary CMB anisotropies, as well as instrumental noise.

Radio interferometry probes astrophysical signals through incomplete and noisy Fourier measurements. The theory of compressed sensing demonstrates that such measurements may actually suffice for accurate reconstruction of sparse or compressible signals. We propose new generic imaging techniques based on convex optimization for global minimization problems defined in this context. The versatility of the framework notably allows introduction of specific prior information on the signals, which offers the possibility of significant improvements of reconstruction relative to the standard local matching pursuit algorithm CLEAN used in radio astronomy. We illustrate the potential of the approach by studying reconstruction performances on simulations of two different kinds of signals observed with very generic interferometric configurations. The first kind is an intensity field of compact astrophysical objects. The second kind is the imprint of cosmic strings in the temperature field of the cosmic microwave background radiation, of particular interest for cosmology.

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.

Using local morphological measures on the sphere defined through a steerable wavelet analysis, we examine the three-year WMAP and the NVSS data for correlation induced by the integrated Sachs-Wolfe (ISW) effect. The steerable wavelet constructed from the second derivative of a Gaussian allows one to define three local morphological measures, namely the signed-intensity, orientation and elongation of local features. Detections of correlation between the WMAP and NVSS data are made with each of these morphological measures. The most significant detection is obtained in the correlation of the signed-intensity of local features at a significance of 99.9%. By inspecting signed-intensity sky maps, it is possible for the first time to see the correlation between the WMAP and NVSS data by eye. Foreground contamination and instrumental systematics in the WMAP data are ruled out as the source of all significant detections of correlation. Our results provide new insight on the ISW effect by probing the morphological nature of the correlation induced between the cosmic microwave background and large scale structure of the Universe. Given the current constraints on the flatness of the Universe, our detection of the ISW effect again provides direct and independent evidence for dark energy. Moreover, this new morphological analysis may be used in future to help us to better understand the nature of dark energy.

The decomposition of a signal on the sphere with the steerable wavelet constructed from the second Gaussian derivative gives access to the orientation, signed-intensity, and elongation of the signal's local features. In the present work, the non-Gaussianity of the WMAP temperature data of the cosmic microwave background (CMB) is analyzed in terms of the first four moments of the statistically isotropic random fields associated with these local morphological measures, at wavelet scales corresponding to angular sizes between 27.5 arcminutes and 30 degrees on the celestial sphere. While no detection is made neither in the orientation analysis nor in the elongation analysis, a strong detection is made in the excess kurtosis of the signed-intensity of the WMAP data. The non-Gaussianity is observed with a significance level below 0.5% at a wavelet scale corresponding to an angular size around 10 degrees, and confirmed at neighbour scales. This supports a previous detection of an excess of kurtosis in the wavelet coefficient of the WMAP data with the axisymmetric Mexican hat wavelet (Vielva et al. 2004). Instrumental noise and foreground emissions are not likely to be at the origin of the excess of kurtosis. Large-scale modulations of the CMB related to some unknown systematics are rejected as possible origins of the detection. The observed non-Gaussianity may therefore probably be imputed to the CMB itself, thereby questioning the basic inflationary scenario upon which the present concordance cosmological model relies. Taking the CMB temperature angular power spectrum of the concordance cosmological model at face value, further analysis also suggests that this non-Gaussianity is not confined to the directions on the celestial sphere with an anomalous signed-intensity.

A fast and exact algorithm is developed for the spin +-2 spherical harmonics transforms on equi-angular pixelizations on the sphere. It is based on the Driscoll and Healy fast scalar spherical harmonics transform. The theoretical exactness of the transform relies on a sampling theorem. The associated asymptotic complexity is of order O(L^2 log^2_2(L)), where 2L stands for the square-root of the number of sampling points on the sphere, also setting a band limit L for the spin +-2 functions considered. The algorithm is presented as an alternative to existing fast algorithms with an asymptotic complexity of order O(L^3) on other pixelizations. We also illustrate these generic developments through their application in cosmology, for the analysis of the cosmic microwave background (CMB) polarization data.

Significant alignment and signed-intensity anomalies of local features of the cosmic microwave background (CMB) are detected on the three-year WMAP data, through a decomposition of the signal with steerable wavelets on the sphere. Firstly, an alignment analysis identifies two mean preferred planes in the sky, both with normal axes close to the CMB dipole axis. The first plane is defined by the directions toward which local CMB features are anomalously aligned. A mean preferred axis is also identified in this plane, located very close to the ecliptic poles axis. The second plane is defined by the directions anomalously avoided by local CMB features. This alignment anomaly provides further insight on recent results (Wiaux et al. 2006). Secondly, a signed-intensity analysis identifies three mean preferred directions in the southern galactic hemisphere with anomalously high or low temperature of local CMB features: a cold spot essentially identified with a known cold spot (Vielva et al. 2004), a second cold spot lying very close to the southern end of the CMB dipole axis, and a hot spot lying close to the southern end of the ecliptic poles axis. In both analyses, the anomalies are observed at wavelet scales corresponding to angular sizes around 10 degress on the celestial sphere, with global significance levels around 1%. Further investigation reveals that the alignment and signed-intensity anomalies are only very partially related. Instrumental noise, foreground emissions, as well as some form of other systematics, are strongly rejected as possible origins of the detections. An explanation might still be envisaged in terms of a global violation of the isotropy of the Universe, inducing an intrinsic statistical anisotropy of the CMB.

A Nordtvedt effect at cosmological scales affects the acoustic oscillations imprinted in the cosmic microwave background. The gravitational baryonic mass density of the universe is inferred at the first peak scale from WMAP data. The independent determination of the inertial baryonic mass density through the measurement of the deuterium abundance in the framework of standard big bang nucleosynthesis leads to a new constraint on a possible violation of the strong equivalence principle at the recombination time.

A fast algorithm is developed for the directional correlation of scalar band-limited signals and band-limited steerable filters on the sphere. The asymptotic complexity associated to it through simple quadrature is of order O(L^5), where 2L stands for the square-root of the number of sampling points on the sphere, also setting a band limit L for the signals and filters considered. The filter steerability allows to compute the directional correlation uniquely in terms of direct and inverse scalar spherical harmonics transforms, which drive the overall asymptotic complexity. The separation of variables technique for the scalar spherical harmonics transform produces an O(L^3) algorithm independently of the pixelization. On equi-angular pixelizations, a sampling theorem introduced by Driscoll and Healy implies the exactness of the algorithm. The equi-angular and HEALPix implementations are compared in terms of memory requirements, computation times, and numerical stability. The computation times for the scalar transform, and hence for the directional correlation, of maps of several megapixels on the sphere (L~10^3) are reduced from years to tens of seconds in both implementations on a single standard computer. These generic results for the scale-space signal processing on the sphere are specifically developed in the perspective of the wavelet analysis of the cosmic microwave background (CMB) temperature (T) and polarization (E and B) maps of the WMAP and Planck experiments. As an illustration, we consider the computation of the wavelet coefficients of a simulated temperature map of several megapixels with the second Gaussian derivative wavelet.

This paper reviews the application of a novel methodology for analysing the isotropy of the universe by probing the alignment of local structures in the CMB. The strength of the proposed methodology relies on the steerable wavelet filtering of the CMB signal. One the one hand, the filter steerability renders the computation of the local orientation of the CMB features affordable in terms of computation time. On the other hand, the scale-space nature of the wavelet filtering allows to explore the alignment of the local structures at different scales, probing possible different phenomena. We present the WMAP first-year data analysis recently performed by the same authors (Wiaux et al.), where an extremely significant anisotropy was found. In particular, a preferred plane was detected, having a normal direction with a northern end position close to the northern end of the CMB dipole axis. In addition, a most preferred direction was found in that plane, with a northern end direction very close to the north ecliptic pole. This result synthesised for the first time previously reported anomalies identified in the direction of the dipole and the ecliptic poles axes. In a forthcoming paper (Vielva et al.), we have extended our analysis to the study of individual frequency maps finding first indications for discarding foregrounds as the origin of the anomaly. We have also tested that the preferred orientations are defined by structures homogeneously distributed in the sky, rather than from localised regions. We have also analysed the WMAP 3-year data, finding the same anomaly pattern, although at a slightly lower significance level.

We question the global universe isotropy by probing the alignment of local structures in the cosmic microwave background (CMB) radiation. The original method proposed relies on a steerable wavelet decomposition of the CMB signal on the sphere. The analysis of the first-year WMAP data identifies a mean preferred plane with a normal direction close to the CMB dipole axis, and a mean preferred direction in this plane, very close to the ecliptic poles axis. Previous statistical anisotropy results are thereby synthesized, but further analyses are still required to establish their origin.

Wavelets on the sphere are reintroduced and further developed independently of the original group theoretic formalism, in an equivalent, but more straightforward approach. These developments are motivated by the interest of the scale-space analysis of the cosmic microwave background (CMB) anisotropies on the sky. A new, self-consistent, and practical approach to the wavelet filtering on the sphere is developed. It is also established that the inverse stereographic projection of a wavelet on the plane (i.e. Euclidean wavelet) leads to a wavelet on the sphere (i.e. spherical wavelet). This new correspondence principle simplifies the construction of wavelets on the sphere and allows to transfer onto the sphere properties of wavelets on the plane. In that regard, we define and develop the notions of directionality and steerability of filters on the sphere. In the context of the CMB analysis, these notions are important tools for the identification of local directional features in the wavelet coefficients of the signal, and for their interpretation as possible signatures of non-gaussianity, statistical anisotropy, or foreground emission. But the generic results exposed may find numerous applications beyond cosmology and astrophysics.

We study the effect of a violation of the strong equivalence principle (SEP) on the cosmic microwave background (CMB). Such a violation would modify the weight of baryons in the primordial gravitational potentials and hence their impact in the establishment of the photon-baryon plasma acoustic oscillations before recombination. This cosmological Nordtvedt effect alters the odd peaks height of the CMB temperature anisotropy power spectrum. A gravitational baryonic mass density of the universe may already be inferred at the first peak scale from the analysis of WMAP data. Experimental constraints on a primordial SEP violation are derived from a comparison with the universe's inertial baryonic mass density measured either in a full analysis of the CMB, or in the framework of the standard big bang nucleosynthesis (BBN).