Fully bayesian blind source separation of astrophysical images modelled by mixture of Gaussians

In this work, we address the problem of source separation in the presence of prior information. We develop a fully Bayesian source separation technique which assumes a generic model for the sources, namely Gaussian mixtures with a priori unknown number of components and utilise Markov chain Monte Carlo techniques for model parameter estimation. The development of this methodology is motivated by the need to bring an efficient solution to the separation of components in the microwave radiation maps to be obtained by the satellite mission PLANCK which has the objective of uncovering cosmic microwave background radiation. The proposed algorithm successfully incorporates a rich variety of prior information available to us in this problem in contrast to most of the previous work which assume completely blind separation of the sources. We report results on realistic simulations of expected Planck maps and on WMAP 3rd year results. The technique suggested is easily applicable to other source separation applications by modifying some of the priors.