We propose a Bayesian approach to joint source separation and restoration for astrophysical diffuse sources. We constitute a prior statistical model for the source images by using their gradient maps. We assume a t-distribution for the gradient maps in different directions, because it is able to fit both smooth and sparse data. A Monte Carlo technique, called Langevin sampler, is used to estimate the source images and all the model parameters are estimated by using deterministic techniques.

The detection and flux estimation of point sources in cosmic microwave background (CMB) maps is a very important task in order to clean the maps and also to obtain relevant astrophysical information. In this paper we propose a maximum a posteriori (MAP) approach detection method in a Bayesian scheme which incorporates prior information about the source flux distribution, the locations and the number of sources. We apply this method to CMB simulations with the characteristics of the Planck satellite channels at 30, 44, 70 and 100 GHz. With a similar level of spurious sources, our method yields more complete catalogues than the matched filter with a 5 sigma threshold. Besides, the new technique allows us to fix the number of detected sources in a non-arbitrary way.

We propose to model the image differentials of astrophysical source maps by Student's t-distribution and to use them in the Bayesian source separation method as priors. We introduce an efficient Markov Chain Monte Carlo (MCMC) sampling scheme to unmix the astrophysical sources and describe the derivation details. In this scheme, we use the Langevin stochastic equation for transitions, which enables parallel drawing of random samples from the posterior, and reduces the computation time significantly (by two orders of magnitude). In addition, Student's t-distribution parameters are updated throughout the iterations. The results on astrophysical source separation are assessed with two performance criteria defined in the pixel and the frequency domains.

We present a new approach to the statistical study and modelling of number counts of faint point sources in astronomical images, i.e. counts of sources whose flux falls below the detection limit of a survey. The approach is based on the theory of alpha-stable distributions. We show that the non-Gaussian distribution of the intensity fluctuations produced by a generic point source population -- whose number counts follow a simple power law -- belongs to the alpha-stable family of distributions. Even if source counts do not follow a simple power law, we show that the alpha-stable model is still useful in many astrophysical scenarios. With the alpha-stable model it is possible to totally describe the non-Gaussian distribution with a few parameters which are closely related to the parameters describing the source counts, instead of an infinite number of moments. Using statistical tools available in the signal processing literature, we show how to estimate these parameters in an easy and fast way. We demonstrate that the model proves valid when applied to realistic point source number counts at microwave frequencies. In the case of point extragalactic sources observed at CMB frecuencies, our technique is able to successfully fitting the P(D) distribution of deflections and to precisely determining the main parameters which describe the number counts. In the case of the Planck mission, the relative errors on these parameters are small either at low and at high frequencies. We provide a way to deal with the presence of Gaussian noise in the data using the empirical characteristic function of the P(D). The formalism and methods here presented can be very useful also for experiments in other frequency ranges, e.g. X-ray or radio Astronomy.