Restoration of interferometric images. III. Efficient Richardson-Lucy methods for LINC-NIRVANA data reduction

In previous papers we proposed methods and software for the restoration of images provided by Fizeau interferometers such as LINC-NIRVANA (LN), the German-Italian beam combiner for the Large Binocular Telescope (LBT). It will provide multiple images of the same target corresponding to different orientations of the baseline. Therefore LN will require routinely the use of multiple-image deconvolution methods in order to produce a unique high-resolution image. As a consequence of the complexity of astronomical images, two kinds of methods will be required: first a quick-look method , namely a method that is computationally efficient, allowing a rapid overview and identification of the object being observed; second an ad hoc method designed for that particular object and as accurate as possible. In this paper we investigate the possibility of using Richardson-Lucy-like (RL-like) methods, namely methods designed for the maximization of the likelihood function in the case of Poisson noise, as possible quick-look methods. To this purpose we propose new techniques for accelerating the Ordered Subsets - Expectation Maximization (OS-EM) method, investigated in our previous papers; moreover, we analyze approaches based on the fusion of the multiple images into a single one, so that one can use single-image deconvolution methods which are presumably more efficient than the multiple-image ones. The results are encouraging and all the methods proposed in this paper have been implemented in our software package AIRY.

[1]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[2]  Eiichi Tanaka A Fast Reconstruction Algorthm for Stationary Positron Emission Tomography Based on a Modified EM Algorithm , 1987, IEEE Transactions on Medical Imaging.

[3]  Robert J. Hanisch,et al.  The restoration of HST images and spectra - II , 2015 .

[4]  T. Holmes,et al.  Acceleration of Maximum-Likelihood Image-Restoration for Fluorescence Microscopy and Other Noncoherent Imagery , 1991, Quantum Limited Imaging and Information Processing.

[5]  Marcel Carbillet,et al.  Restoration of interferometric images - I. The software package AIRY , 2002 .

[6]  L. Lucy An iterative technique for the rectification of observed distributions , 1974 .

[7]  H. Lantéri,et al.  Penalized maximum likelihood image restoration with positivity constraints:multiplicative algorithms , 2002 .

[8]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[9]  H Lantéri,et al.  A general method to devise maximum-likelihood signal restoration multiplicative algorithms with non-negativity constraints , 2001, Signal Process..

[10]  H. Malcolm Hudson,et al.  Accelerated image reconstruction using ordered subsets of projection data , 1994, IEEE Trans. Medical Imaging.

[11]  W. Wamsteker,et al.  A simple but efficient algorithm for multiple-image deblurring , 2004 .

[12]  Marcel Carbillet,et al.  Restoration of interferometric images II. The case-study of the Large Binocular Telescope , 2002 .

[13]  M. Bertero,et al.  Application of the OS-EM method to the restoration of LBT images , 2000 .

[14]  Mario Bertero,et al.  Image restoration methods for the Large Binocular Telescope (LBT) , 2000 .

[15]  I. Csiszár Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems , 1991 .

[16]  Frank Natterer,et al.  Mathematical methods in image reconstruction , 2001, SIAM monographs on mathematical modeling and computation.

[17]  R. White,et al.  Image recovery from data acquired with a charge-coupled-device camera. , 1993, Journal of the Optical Society of America. A, Optics and image science.

[18]  E. R. Podolyak,et al.  High resolution soft X-ray spectrum reconstruction by MWPC attenuation measurements , 1983 .