  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).
