Multiple image deblurring with high dynamic-range Poisson data

An interesting problem arising in astronomical imaging is the reconstruction of an image with high dynamic range, for example a set of bright point sources superimposed to smooth structures. A few methods have been proposed for dealing with this problem and their performance is not always satisfactory. In this paper we propose a solution based on the representation, already proposed elsewhere, of the image as the sum of a pointwise component and a smooth one, with different regularization for the two components. Our approach is in the framework of Poisson data and to this purpose we need efficient deconvolution methods. Therefore, first we briefly describe the application of the Scaled Gradient Projection (SGP) method to the case of different regularization schemes and subsequently we propose how to apply these methods to the case of multiple image deconvolution of high-dynamic range images, with specific reference to the Fizeau interferometer LBTI of the Large Binocular Telescope (LBT). The efficacy of the proposed methods is illustrated both on simulated images and on real images, observed with LBTI, of the Jovian moon Io. The software is available at http://www.oasis.unimore.it/site/home/software.html