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 measu
rements 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. Goo
d 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.
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. (20
05). 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.
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 signals local features. In the present work, the non-Gauss
ianity 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.
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 t
he sphere. Two fast algorithms are also presented for the analysis of signals on the sphere with steerable wavelets.
