3D deconvolution of hyper-spectral astronomical data

• Published in: Monthly notices of the Royal Astronomical Society Vol. 418; no. 1; pp. 258 - 270
• Main Authors: Bongard, S., Soulez, F., Thiébaut, É., Pecontal, É.
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Publishing Ltd 01.11.2011; Wiley-Blackwell; Oxford University Press
• Subjects: 3-D graphics; Astronomical research; Astronomy; cosmology: observations; Data analysis; Earth, ocean, space; Exact sciences and technology; instrumentation: spectrographs; methods: data analysis; Spectrum analysis; techniques: image processing; techniques: imaging spectroscopy; cosmology: observations; methods: data analysis; instrumentation: spectrographs; techniques: image processing; techniques: imaging spectroscopy; Spectroscopy; Data analysis; Galaxies; Supernovae; Image processing; Image restoration; Deconvolution; Cosmology; Dynamic model; Regularization
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1111/j.1365-2966.2011.19480.x

In this paper we present a general method for multichannel image restoration based on regularized χ2. We introduce separable regularizations that account for the dynamics of the model and take advantage of the continuities present in the data, leaving only two hyper-parameters to tune. We illustrate a practical implementation of this method in the context of host galaxy subtraction for the Nearby SuperNova Factory (SNfactory). We show that the image restoration obtained fulfils the stringent requirements on bias and photometricity needed by this programme. The reconstruction yields sub-per cent integrated residuals in all the synthetic filters considered both on real and simulated data. Even though our implementation is tied to the SNfactory data, the method translates to any hyper-spectral data. As such, it is of direct relevance to several new generation instruments like MUSE. Also, this technique could be applied to multiband astronomical imaging for which image reconstruction is important, for example, to increase image resolution for weak-lensing surveys.