Polarized wavelets and curvelets on the sphere

The statistics of the temperature anisotropies in the primordial cosmic microwave background radiation field provide a wealth of information for cosmology and for estimating cosmological parameters. An even more acute inference should stem from the study of maps of the polarization state of the CMB radiation. Measuring the extremely weak CMB polarization signal requires very sensitive instruments. The full-sky maps of both temperature and polarization anisotropies of the CMB to be delivered by the upcoming Planck Surveyor satellite experiment are hence being awaited with excitement. Multiscale methods, such as isotropic wavelets, steerable wavelets, or curvelets, have been proposed in the past to analyze the CMB temperature map. In this paper, we contribute to enlarging the set of available transforms for polarized data on the sphere. We describe a set of new multiscale decompositions for polarized data on the sphere, including decimated and undecimated Q-U or E-B wavelet transforms and Q-U or E-B curvelets. The proposed transforms are invertible and so allow for applications in data restoration and denoising.Comment: Accepted. Full paper will figures available at http://jstarck.free.fr/aa08_pola.pd

Sparsity and morphological diversity for hyperspectral data analysis

Recently morphological diversity and sparsity have emerged as new and effective sources of diversity for Blind Source Separation. Based on these new concepts, novelmethods such as Generalized Morphological Component Analysis have been put forward. The latter takes advantage of the very sparse representation of structured data in large overcomplete dictionaries, to separate sources based on their morphology. Building on GMCA, the purpose of this contribution is to describe a new algorithm for hyperspectral data processing. Large-scale hyperspectral data refers to collected data that exhibit sparse spectral signatures in addition to sparse spatial morphologies, in specified dictionaries of spectral and spatial waveforms. Numerical experiments are reported which demonstrate the validity of the proposed extension for solving source separation problems involving hyperspectral data

Sunyaev-Zel'dovich clusters reconstruction in multiband bolometer camera surveys

We present a new method for the reconstruction of Sunyaev-Zel'dovich (SZ) galaxy clusters in future SZ-survey experiments using multiband bolometer cameras such as Olimpo, APEX, or Planck. Our goal is to optimise SZ-Cluster extraction from our observed noisy maps. We wish to emphasize that none of the algorithms used in the detection chain is tuned on prior knowledge on the SZ -Cluster signal, or other astrophysical sources (Optical Spectrum, Noise Covariance Matrix, or covariance of SZ Cluster wavelet coefficients). First, a blind separation of the different astrophysical components which contribute to the observations is conducted using an Independent Component Analysis (ICA) method. Then, a recent non linear filtering technique in the wavelet domain, based on multiscale entropy and the False Discovery Rate (FDR) method, is used to detect and reconstruct the galaxy clusters. Finally, we use the Source Extractor software to identify the detected clusters. The proposed method was applied on realistic simulations of observations. As for global detection efficiency, this new method is impressive as it provides comparable results to Pierpaoli et al. method being however a blind algorithm. Preprint with full resolution figures is available at the URL: w10-dapnia.saclay.cea.fr/Phocea/Vie_des_labos/Ast/ast_visu.php?id_ast=728Comment: Submitted to A&A. 32 Pages, text onl

SZ and CMB reconstruction using Generalized Morphological Component Analysis

In the last decade, the study of cosmic microwave background (CMB) data has become one of the most powerful tools to study and understand the Universe. More precisely, measuring the CMB power spectrum leads to the estimation of most cosmological parameters. Nevertheless, accessing such precious physical information requires extracting several different astrophysical components from the data. Recovering those astrophysical sources (CMB, Sunyaev-Zel'dovich clusters, galactic dust) thus amounts to a component separation problem which has already led to an intense activity in the field of CMB studies. In this paper, we introduce a new sparsity-based component separation method coined Generalized Morphological Component Analysis (GMCA). The GMCA approach is formulated in a Bayesian maximum a posteriori (MAP) framework. Numerical results show that this new source recovery technique performs well compared to state-of-the-art component separation methods already applied to CMB data.Comment: 11 pages - Statistical Methodology - Special Issue on Astrostatistics - in pres

Sparsity constraints for hyperspectral data analysis: linear mixture model and beyond

The recent development of multi-channel sensors has motivated interest in devising new methods for the coherent processing of multivariate data. An extensive work has already been dedicated to multivariate data processing ranging from blind source separation (BSS) to multi/hyper-spectral data restoration. Previous work has emphasized on the fundamental role played by sparsity and morphological diversity to enhance multichannel signal processing. GMCA is a recent algorithm for multichannel data analysis which was used successfully in a variety of applications including multichannel sparse decomposition, blind source separation (BSS), color image restoration and inpainting. Inspired by GMCA, a recently introduced algorithm coined HypGMCA is described for BSS applications in hyperspectral data processing. It assumes the collected data is a linear instantaneous mixture of components exhibiting sparse spectral signatures as well as sparse spatial morphologies, each in specified dictionaries of spectral and spatial waveforms. We report on numerical experiments with synthetic data and application to real observations which demonstrate the validity of the proposed method

Blind component separation in wavelet space. Application to CMB analysis

It is a recurrent issue in astronomical data analysis that observations are unevenly sampled or incomplete maps with missing patches or intentionaly masked parts. In addition, many astrophysical emissions are non stationary processes over the sky. Hence spectral estimation using standard Fourier transforms is no longer reliable. Spectral matching ICA (SMICA) is a source separation method based on covariance matching in Fourier space which is successfully used for the separation of diffuse astrophysical emissions in Cosmic Microwave Background observations. We show here that wavelets, which are standard tools in processing non stationary data, can profitably be used to extend SMICA. Among possible applications, it is shown that gaps in data are dealt with more conveniently and with better results using this extension, wSMICA, in place of the original SMICA. The performances of these two methods are compared on simulated CMB data sets, demonstrating the advantageous use of wavelets

Blind Source Separation: the Sparsity Revolution

International audienceOver the last few years, the development of multi-channel sensors motivated interest in methods for the coherent processing of multivariate data. Some specific issues have already been addressed as testified by the wide literature on the so-called blind source separation (BSS) problem. In this context, as clearly emphasized by previous work, it is fundamental that the sources to be retrieved present some quantitatively measurable diversity. Recently, sparsity and morphological diversity have emerged as a novel and effective source of diversity for BSS. We give here some essential insights into the use of sparsity in source separation and we outline the essential role of morphological diversity as being a source of diversity or contrast between the sources. This paper overviews a sparsity-based BSS method coined Generalized Morphological Component Analysis (GMCA) that takes advantages of both morphological diversity and sparsity, using recent sparse overcomplete or redundant signal representations. GMCA is a fast and efficient blind source separation method. In remote sensing applications, the specificity of hyperspectral data should be accounted for. We extend the proposed GMCA framework to deal with hyperspectral data. In a general framework, GMCA provides a basis for multivariate data analysis in the scope of a wide range of classical multivariate data restorate. Numerical results are given in color image denoising and inpainting. Finally, GMCA is applied to the simulated ESA/Planck data. It is shown to give effective astrophysical component separation

Wavelets, ridgelets and curvelets on the sphere

We present in this paper new multiscale transforms on the sphere, namely the isotropic undecimated wavelet transform, the pyramidal wavelet transform, the ridgelet transform and the curvelet transform. All of these transforms can be inverted i.e. we can exactly reconstruct the original data from its coefficients in either representation. Several applications are described. We show how these transforms can be used in denoising and especially in a Combined Filtering Method, which uses both the wavelet and the curvelet transforms, thus benefiting from the advantages of both transforms. An application to component separation from multichannel data mapped to the sphere is also described in which we take advantage of moving to a wavelet representation.Comment: Accepted for publication in A&A. Manuscript with all figures can be downloaded at http://jstarck.free.fr/aa_sphere05.pd

First evidence for a gravitational lensing-induced echo in gamma rays with Fermi LAT

Aims. This article shows the first evidence for gravitational lensing phenomena in high energy gamma-rays. This evidence comes from the observation of a gravitational lens induced echo in the light curve of the distant blazar PKS 1830-211. Methods. Traditional methods for the estimation of time delays in gravitational lensing systems rely on the cross-correlation of the light curves of the individual images. In this paper, we use 300 MeV-30 GeV photons detected by the Fermi-LAT instrument. The Fermi-LAT instrument cannot separate the images of known lenses. The observed light curve is thus the superposition of individual image light curves. The Fermi-LAT instrument has the advantage of providing long, evenly spaced, time series. In addition, the photon noise level is very low. This allows to use directly Fourier transform methods. Results. A time delay between the two compact images of PKS 1830-211 has been searched for both by the autocorrelation method and the "double power spectrum" method. The double power spectrum shows a 3 {\sigma} evidence for a time delay of 27.5 $\pm$ 1.3 days, consistent with the result from Lovell et al. (1998). The relative uncertainty on the time delay estimation is reduced from 20% to 5%.Comment: submitted to A&

Structure and Spin Dynamics of La0.85_{0.85}Sr0.15_{0.15}MnO3_3

Neutron scattering has been used to study the structure and spin dynamics of La$_{0.85}$Sr$_{0.15}$MnO$_3$. The magnetic structure of this system is ferromagnetic below T_C = 235 K. We see anomalies in the Bragg peak intensities and new superlattice peaks consistent with the onset of a spin-canted phase below T_{CA} = 205 K, which appears to be associated with a gap at q = (0, 0, 0.5) in the spin-wave spectrum. Anomalies in the lattice parameters indicate a concomitant lattice distortion. The long-wavelength magnetic excitations are found to be conventional spin waves, with a gapless (< 0.02 meV) isotropic dispersion relation $E = Dq^2$. The spin stiffness constant D has a $T^{5/2}$ dependence at low T, and the damping at small q follows $q^4T^{2}$. An anomalously strong quasielastic component, however, develops at small wave vector above 200 K and dominates the fluctuation spectrum as T -> T_C. At larger q, on the other hand, the magnetic excitations become heavily damped at low temperatures, indicating that spin waves in this regime are not eigenstates of the system, while raising the temperature dramatically increases the damping. The strength of the spin-wave damping also depends strongly on the symmetry direction in the crystal. These anomalous damping effects are likely due to the itinerant character of the $e_g$ electrons.Comment: 8 pages (RevTex), 9 figures (encapsulated postscript