106 research outputs found
An algorithm for hybrid regularizers based image restoration with Poisson noise
summary:In this paper, a hybrid regularizers model for Poissonian image restoration is introduced. We study existence and uniqueness of minimizer for this model. To solve the resulting minimization problem, we employ the alternating minimization method with rigorous convergence guarantee. Numerical results demonstrate the efficiency and stability of the proposed method for suppressing Poisson noise
Restoration of Poissonian Images Using Alternating Direction Optimization
Much research has been devoted to the problem of restoring Poissonian images,
namely for medical and astronomical applications. However, the restoration of
these images using state-of-the-art regularizers (such as those based on
multiscale representations or total variation) is still an active research
area, since the associated optimization problems are quite challenging. In this
paper, we propose an approach to deconvolving Poissonian images, which is based
on an alternating direction optimization method. The standard regularization
(or maximum a posteriori) restoration criterion, which combines the Poisson
log-likelihood with a (non-smooth) convex regularizer (log-prior), leads to
hard optimization problems: the log-likelihood is non-quadratic and
non-separable, the regularizer is non-smooth, and there is a non-negativity
constraint. Using standard convex analysis tools, we present sufficient
conditions for existence and uniqueness of solutions of these optimization
problems, for several types of regularizers: total-variation, frame-based
analysis, and frame-based synthesis. We attack these problems with an instance
of the alternating direction method of multipliers (ADMM), which belongs to the
family of augmented Lagrangian algorithms. We study sufficient conditions for
convergence and show that these are satisfied, either under total-variation or
frame-based (analysis and synthesis) regularization. The resulting algorithms
are shown to outperform alternative state-of-the-art methods, both in terms of
speed and restoration accuracy.Comment: 12 pages, 12 figures, 2 tables. Submitted to the IEEE Transactions on
Image Processin
Inexact Bregman iteration with an application to Poisson data reconstruction
This work deals with the solution of image restoration problems by an
iterative regularization method based on the Bregman iteration. Any iteration of this
scheme requires to exactly compute the minimizer of a function. However, in some
image reconstruction applications, it is either impossible or extremely expensive to
obtain exact solutions of these subproblems. In this paper, we propose an inexact
version of the iterative procedure, where the inexactness in the inner subproblem
solution is controlled by a criterion that preserves the convergence of the Bregman
iteration and its features in image restoration problems. In particular, the method
allows to obtain accurate reconstructions also when only an overestimation of the
regularization parameter is known. The introduction of the inexactness in the iterative
scheme allows to address image reconstruction problems from data corrupted by
Poisson noise, exploiting the recent advances about specialized algorithms for the
numerical minimization of the generalized Kullback–Leibler divergence combined with
a regularization term. The results of several numerical experiments enable to evaluat
ACQUIRE: an inexact iteratively reweighted norm approach for TV-based Poisson image restoration
We propose a method, called ACQUIRE, for the solution of constrained
optimization problems modeling the restoration of images corrupted by Poisson
noise. The objective function is the sum of a generalized Kullback-Leibler
divergence term and a TV regularizer, subject to nonnegativity and possibly
other constraints, such as flux conservation. ACQUIRE is a line-search method
that considers a smoothed version of TV, based on a Huber-like function, and
computes the search directions by minimizing quadratic approximations of the
problem, built by exploiting some second-order information. A classical
second-order Taylor approximation is used for the Kullback-Leibler term and an
iteratively reweighted norm approach for the smoothed TV term. We prove that
the sequence generated by the method has a subsequence converging to a
minimizer of the smoothed problem and any limit point is a minimizer.
Furthermore, if the problem is strictly convex, the whole sequence is
convergent. We note that convergence is achieved without requiring the exact
minimization of the quadratic subproblems; low accuracy in this minimization
can be used in practice, as shown by numerical results. Experiments on
reference test problems show that our method is competitive with
well-established methods for TV-based Poisson image restoration, in terms of
both computational efficiency and image quality.Comment: 37 pages, 13 figure
An Augmented Lagrangian Approach to the Constrained Optimization Formulation of Imaging Inverse Problems
We propose a new fast algorithm for solving one of the standard approaches to
ill-posed linear inverse problems (IPLIP), where a (possibly non-smooth)
regularizer is minimized under the constraint that the solution explains the
observations sufficiently well. Although the regularizer and constraint are
usually convex, several particular features of these problems (huge
dimensionality, non-smoothness) preclude the use of off-the-shelf optimization
tools and have stimulated a considerable amount of research. In this paper, we
propose a new efficient algorithm to handle one class of constrained problems
(often known as basis pursuit denoising) tailored to image recovery
applications. The proposed algorithm, which belongs to the family of augmented
Lagrangian methods, can be used to deal with a variety of imaging IPLIP,
including deconvolution and reconstruction from compressive observations (such
as MRI), using either total-variation or wavelet-based (or, more generally,
frame-based) regularization. The proposed algorithm is an instance of the
so-called "alternating direction method of multipliers", for which convergence
sufficient conditions are known; we show that these conditions are satisfied by
the proposed algorithm. Experiments on a set of image restoration and
reconstruction benchmark problems show that the proposed algorithm is a strong
contender for the state-of-the-art.Comment: 13 pages, 8 figure, 8 tables. Submitted to the IEEE Transactions on
Image Processin
Directional TGV-based image restoration under Poisson noise
We are interested in the restoration of noisy and blurry images where the
texture mainly follows a single direction (i.e., directional images). Problems
of this type arise, for example, in microscopy or computed tomography for
carbon or glass fibres. In order to deal with these problems, the Directional
Total Generalized Variation (DTGV) was developed by Kongskov et al. in 2017 and
2019, in the case of impulse and Gaussian noise. In this article we focus on
images corrupted by Poisson noise, extending the DTGV regularization to image
restoration models where the data fitting term is the generalized
Kullback-Leibler divergence. We also propose a technique for the identification
of the main texture direction, which improves upon the techniques used in the
aforementioned work about DTGV. We solve the problem by an ADMM algorithm with
proven convergence and subproblems that can be solved exactly at a low
computational cost. Numerical results on both phantom and real images
demonstrate the effectiveness of our approach.Comment: 20 pages, 1 table, 13 figure
- …