82 research outputs found
A convergent blind deconvolution method for post-adaptive-optics astronomical imaging
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
Shearlet-based regularized reconstruction in region-of-interest computed tomography
Region of interest (ROI) tomography has gained increasing attention in recent years due to its potential to reducing radiation exposure and shortening the scanning time. However, tomographic reconstruction from ROI-focused illumination involves truncated projection data and typically results in higher numerical instability even when the reconstruction problem has unique solution. To address this problem, both ad hoc analytic formulas and iterative numerical schemes have been proposed in the literature. In this paper, we introduce a novel approach for ROI tomographic reconstruction, formulated as a convex optimization problem with a regularized term based on shearlets. Our numerical implementation consists of an iterative scheme based on the scaled gradient projection method and it is tested in the context of fan-beam CT. Our results show that our approach is essentially insensitive to the location of the ROI and remains very stable also when the ROI size is rather small.Peer reviewe
Preconditioned ADMM with nonlinear operator constraint
We are presenting a modification of the well-known Alternating Direction
Method of Multipliers (ADMM) algorithm with additional preconditioning that
aims at solving convex optimisation problems with nonlinear operator
constraints. Connections to the recently developed Nonlinear Primal-Dual Hybrid
Gradient Method (NL-PDHGM) are presented, and the algorithm is demonstrated to
handle the nonlinear inverse problem of parallel Magnetic Resonance Imaging
(MRI)
- …