Bayesian fusion of hyperspectral astronomical images

André Jalobeanu, Matthieu Petremand, Christophe Collet, Ali Mohammad-Djafari, Jean-François Bercher, Pierre Bessiére
The new integral-field spectrograph MUSE will acquire hyperspectral images of the deep sky, requiring huge amounts of raw data to be processed, posing a challenge to modern algorithms and technologies. In order to achieve the required sensitivity to observe very faint objects, many observations need to be reconstructed and co-added into a single data cube. In this paper, we propose a new fusion method to combine all raw observations while removing most of the instrumental and observational
more » ... observational artifacts such as blur or cosmic rays. Thus, the results can be accurately and consistently analyzed by astronomers. We use a Bayesian framework allowing for optimal data fusion and uncertainty estimation. The knowledge of the instrument allows to write the direct problem (data acquisition on the detector matrix) and then to invert it through Bayesian inference, assuming a smoothness prior for the data cube to be reconstructed. Compared to existing methods, the originality of the new technique is in the propagation of errors throughout the fusion pipeline and the ability to deal with various acquisition parameters for each input image. For this paper, we focus on small-size, simulated astronomical observations with varying parameters to validate the image formation model, the reconstruction algorithm and the predicted uncertainties.
doi:10.1063/1.3573629 fatcat:xsrucjij2nd7dcrkegcgamengi