Reconstruction of High Dynamic Range Images: Simulations of LBT Observations of a Stellar Jet, a Pathfinder Study for Future AO-Assisted Giant Telescopes

We present simulated Large Binocular Telescope (LBT) infrared narrow-band observations of a star-jet system, in conjunction with improved and optimized deconvolution and image reconstruction algorithms, considering two cases of interest: single-dish direct imaging with an AO-assisted camera and imaging through a Fizeau interferometer that combines the beams of the two mirrors of LBT. We aim at understanding what accuracy can be obtained with the use of present AO-assisted large telescopes (such as LBT) and what improvements an interferometric instrument (such as LINC-NIRVANA) will be able to provide. The proposed deconvolution method is based on the target decomposition as a sum of a point source (the star) and an extended source (the jet). By assuming Poisson noise we add to the negative logarithm of the likelihood a regularization term enforcing smoothness of the jet component. Finally, we use a Richardson-Lucy-like method for the minimization of this function. This approach is an improvement of a method proposed by Lucy in 1994 for accurate photometric restoration of HST images and called two channel photometric restoration. We denote the new method as the multi-component Richardson-Lucy (MC-RL) method. The analysis of the reconstructed objects shows that the MC-RL method applied to the interferometric observations allows us to evaluate the width and the spatial intensity profile of the jet down to 20 mas with an accuracy better than about 20% in the best case of a central star fainter than 10 mag. These limits allow us to obtain a very good reconstruction of the jet acceleration region very close to the exciting source, which would provide fundamental scientific information on the jet collimation degree and eventually on its launching mechanism. As concerns the proposed MC-RL method, it demonstrates a good performance in the reconstruction of images with a very high dynamic range. It can be improved in several directions, by increasing both its efficiency, thanks to recently proposed acceleration techniques, and its accuracy by means of more sophisticated regularization terms. We are also planning to apply the method to simulated observations of upcoming super giant earth-based telescopes.

[1]  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.

[2]  B. Nisini,et al.  1–2.5 $\mu$m spectra of jets from young stars: Strong $\ion{Fe}{ii}$ emission in HH111, HH240-241 and HH120 , 2002 .

[3]  Emiliano Diolaiti,et al.  Analysis of LBT LINC-NIRVANA simulated images of galaxies and young stellar objects , 2008, Astronomical Telescopes + Instrumentation.

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

[5]  Roberto Ragazzoni,et al.  LINC-NIRVANA: a Fizeau beam combiner for the large binocular telescope , 2003, SPIE Astronomical Telescopes + Instrumentation.

[6]  Marcel Carbillet,et al.  Reduction of boundary effects in multiple image deconvolution with an application to LBT LINC-NIRVANA , 2006 .

[7]  Harrison H. Barrett,et al.  Foundations of Image Science , 2003, J. Electronic Imaging.

[8]  M. Bertero,et al.  Image deblurring with Poisson data: from cells to galaxies , 2009 .

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

[10]  Mario Bertero,et al.  Image reconstruction for observations with a high dynamic range: LINC-NIRVANA simulations of a stellar jet , 2012, Other Conferences.

[11]  C. Helstrom,et al.  Compensation for readout noise in CCD images , 1995 .

[12]  C. Vérinaud,et al.  Modelling astronomical adaptive optics – I. The software package caos , 2005 .

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

[14]  J. Giovannelli,et al.  Positive deconvolution for superimposed extended source and point sources , 2005, astro-ph/0507691.

[15]  Mario Bertero,et al.  A simple method for the reduction of boundary effects in the Richardson-Lucy approach to image deconvolution , 2005 .

[16]  Luca Zanni,et al.  A discrepancy principle for Poisson data , 2010 .

[17]  Roberto Ragazzoni,et al.  Layer-oriented simulation tool. , 2004, Applied optics.

[18]  Marcel Carbillet,et al.  AIRY-LN: an ad-hoc numerical tool for deconvolution of images from the LBT instrument LINC-NIRVANA , 2008, Astronomical Telescopes + Instrumentation.

[19]  Jon A. Morse,et al.  Hubble Space Telescope Images of the HH 34 Jet and Bow Shock: Structure and Proper Motions , 2002 .

[20]  Thomas Rodet,et al.  A variational Bayesian approach for unsupervised super-resolution using mixture models of point and smooth sources applied to astrophysical map-making , 2012 .

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

[22]  Johnathan M. Bardsley,et al.  Regularization parameter selection methods for ill-posed Poisson maximum likelihood estimation , 2009 .

[23]  A. Boss,et al.  Protostars and Planets VI , 2000 .

[24]  Valeria Ruggiero,et al.  Inexact Bregman iteration with an application to Poisson data reconstruction , 2013 .

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

[26]  Mario Bertero,et al.  Imaging with LINC-NIRVANA, the Fizeau interferometer of the Large Binocular Telescope: state of the art and open problems , 2011 .

[27]  Marcel Carbillet,et al.  Strehl-constrained iterative blind deconvolution for post-adaptive-optics data , 2009 .

[28]  L. Zanni,et al.  A scaled gradient projection method for constrained image deblurring , 2008 .