Images from adaptive optics systems are generally affected by significant distortions of the point spread function (PSF) across the field of view, depending on the position of natural and artificial guide stars. Image reduction techniques circumventing or mitigating these effects are important tools to take full advantage of the scientific information encoded in AO images. The aim of this paper is to propose a method for the deblurring of the astronomical image, given a set of samples of the space-variant PSF. The method is based on a partitioning of the image domain into regions of isoplanatism and on applying suitable deconvolution methods with boundary effects correction to each region. The effectiveness of the boundary effects correction is proved. Moreover, the criterion for extending the disjoint sections to partially overlapping sections is validated. The method is applied to simulated images of a stellar system characterized by a spatially variable PSF. We obtain good photometric quality, and therefore good science quality, by performing aperture photometry on the deblurred images. The proposed method is implemented in IDL in the Software Package "Patch", which is available on http://www.airyproject.eu.

In the case of ground-based telescopes equipped with adaptive optics systems, the point spread function (PSF) is only poorly known or completely unknown. Moreover, an accurate modeling of the PSF is in general not available. Therefore in several imaging situations the so-called blind deconvolution methods, aiming at estimating both the scientific target and the PSF from the detected image, can be useful. A blind deconvolution problem is severely ill-posed and, in order to reduce the extremely large number of possible solutions, it is necessary to introduce sensible constraints on both the scientific target and the PSF. In a previous paper we proposed a sound mathematical approach based on a suitable inexact alternating minimization strategy for minimizing the generalized Kullback-Leibler divergence, assuring global convergence. In the framework of this method we showed that an important constraint on the PSF is the upper bound which can be derived from the knowledge of its Strehl ratio. The efficacy of the approach was demonstrated by means of numerical simulations. In this paper, besides improving the previous approach by the use of a further constraint on the unknown scientific target, we extend it to the case of multiple images of the same target obtained with different PSFs. The main application we have in mind is to Fizeau interferometry. As it is known this is a special feature of the Large Binocular Telescope (LBT). The method is applied to realistic simulations of imaging both by single mirrors and Fizeau interferometers. Successes and failures of the method in the imaging of stellar fields are demonstrated in simple cases. These preliminary results look promising at least in specific situations. The IDL code of the proposed method is available on request and will be included in the forthcoming version of the Software Package AIRY (v.6.1).

In this paper we propose a blind deconvolution method which applies to data perturbed by Poisson noise. The objective function is a generalized Kullback-Leibler divergence, depending on both the unknown object and unknown point spread function (PSF), without the addition of regularization terms; constrained minimization, with suitable convex constraints on both unknowns, is considered. The problem is nonconvex and we propose to solve it by means of an inexact alternating minimization method, whose global convergence to stationary points of the objective function has been recently proved in a general setting. The method is iterative and each iteration, also called outer iteration, consists of alternating an update of the object and the PSF by means of fixed numbers of iterations, also called inner iterations, of the scaled gradient projection (SGP) method. The use of SGP has two advantages: first, it allows to prove global convergence of the blind method; secondly, it allows the introduction of different constraints on the object and the PSF. The specific constraint on the PSF, besides non-negativity and normalization, is an upper bound derived from the so-called Strehl ratio, which is the ratio between the peak value of an aberrated versus a perfect wavefront. Therefore a typical application is the imaging of modern telescopes equipped with adaptive optics systems for partial correction of the aberrations due to atmospheric turbulence. In the paper we describe the algorithm and we recall the results leading to its convergence. Moreover we illustrate its effectiveness by means of numerical experiments whose results indicate that the method, pushed to convergence, is very promising in the reconstruction of non-dense stellar clusters. The case of more complex astronomical targets is also considered, but in this case regularization by early stopping of the outer iterations is required.

The Richardson-Lucy method is the most popular deconvolution method in astronomy because it preserves the number of counts and the non-negativity of the original object. Regularization is, in general, obtained by an early stopping of Richardson-Lucy iterations. In the case of point-wise objects such as binaries or open star clusters, iterations can be pushed to convergence. However, it is well-known that Richardson-Lucy is an inefficient method. In most cases, acceptable solutions are obtained at the cost of hundreds or thousands of iterations. A general optimization method, referred to as the scaled gradient projection method, has been proposed for the constrained minimization of continuously differentiable convex functions. It is applicable to the non-negative minimization of the Kullback-Leibler divergence. If the scaling suggested by Richardson-Lucy is used in this method, then it provides a considerable increase in the efficiency of Richardson-Lucy. Therefore the aim of this paper is to apply the scaled gradient projection method to a number of imaging problems in astronomy such as single image deconvolution, multiple image deconvolution, and boundary effect correction. The corresponding algorithms are derived and implemented in interactive data language. To attempt to achieve a further increase in efficiency, we also consider an implementation on graphic processing units. The proposed algorithms are tested on simulated images. The acceleration of scaled gradient projection methods achieved with respect to the corresponding Richardson-Lucy methods strongly depends on both the problem and the specific object to be reconstructed, and in our simulations the improvement achieved ranges from about a factor of 4 to more than 30. Moreover, significant accelerations of up to two orders of magnitude have been observed between the serial and parallel implementations of the algorithms.

We report the results of a simulation and reconstruction of observations of a young stellar object (YSO) jet with the LINC-NIRVANA (LN) interferometric instrument, which will be mounted on the Large Binocular Telescope (LBT). This simulation has been performed in order to investigate the ability of observing the weak diffuse jet line emission against the strong IR stellar continuum through narrow band images in the H and K atmospheric windows. In general, this simulation provides clues on the image quality that could be achieved in observations with a high dynamic range. In these cases, standard deconvolution methods, such as Richardson-Lucy, do not provide satisfactory results: we therefore propose here a new method of image reconstruction. It consists in considering the image to be reconstructed as the sum of two terms: one corresponding to the star (whose position is assumed to be known) and the other to the jet. A regularization term is introduced for this second component and the reconstruction is obtained with an iterative method alternating between the two components. An analysis of the results shows that the image quality obtainable with this method is significantly improved with respect to standard deconvolution methods, reducing the number of artifacts and allowing us to reconstruct the original jet intensity distribution with an error smaller than 10%.

We present two wide-field (5'x5'), diffraction limited, broad-band 10 and 20 micron images of the Orion Nebula, plus six 7-13 micron narrow-band images of BN/KL taken at the 3.8m UKIRT telescope. The wide-field images are mosaics of 35"x35" frames reconstructed using a new restoration method developed for this project. They show the filamentary structure of the dust emission and reveal a new remarkable group of arc-like structures ~1' to the South of the Trapezium. The morphology of the Ney-Allen nebula suggests a complex kinematical status at the center of the Cluster. We find indications that the B0.5V star theta1 Ori-D, is surrounded by a photoevaporated circumstellar disk. Among the four historic Trapezium OB stars, this is the only one without a binary companion, suggesting that stellar multiplicity and the presence of massive circumstellar disks may be mutually exclusive. In the BN/KL complex we detect optically thin silicate emission on top of the deep 10 micron absorption feature. We map the foreground optical depth, color temperature and mid-IR luminosity of the embedded sources. We resolve a conspicuous point source at the location of the IRc2-A knot. The spectral profile of the 10 micron feature is compatible with grain crystallization. We detect several point sources In the OMC-1 South region and discuss their association with the mass loss observed at optical and mm wavelengths. Finally, we list the position and photometry of 177 point sources. With two possible exceptions, all point sources detected at shorter wavelengths display 10 micron emission well above the photospheric level. The model of Robberto et al. (2002, ApJ, 578, 897) provides the simplest explanation for the mid-IR excess.

Ground-based astronomical observations at thermal infrared wavelengths face the problem of extracting the weak astronomical signal from the large and rapidly variable background flux. The observing strategy most commonly used, the so-called ``chopping and nodding'' differential technique, provides reliable representations of the target uniquely in the case of compact sources while extended and complex sources can be easily distorted by their negative counterparts. A restoration method, designed to remove the negative values and to provide reliable representations of extended sources has been proposed. We show that the algorithm successfully removes the distortions due to the negative counterparts and, in addition, provides noise reduction.