Abstract:We consider the restoration of piecewise constant images where the number of the regions and their values are not fixed in advance, with a good difference of piecewise constant values between neighboring regions, from noisy data obtained at the output of a linear operator (e.g., a blurring kernel or a Radon transform). Thus we also address the generic problem of unsupervised segmentation in the context of linear inverse problems. The segmentation and the restoration tasks are solved jointly by minimizing an objective function (an energy) composed of a quadratic data-fidelity term and a nonsmooth nonconvex regularization term. The pertinence of such an energy is ensured by the analytical properties of its minimizers. However, its practical interest used to be limited by the difficulty of the computational stage which requires a nonsmooth nonconvex minimization. Indeed, the existing methods are unsatisfactory since they (implicitly or explicitly) involve a smooth approximation of the regularization term and often get stuck in shallow local minima. The goal of this paper is to design a method that efficiently handles the nonsmooth nonconvex minimization. More precisely, we propose a continuation method where one tracks the minimizers along a sequence of approximate nonsmooth energies $\{J_\eps\}$, the first of which being strictly convex and the last one the original energy to minimize. Knowing the importance of the nonsmoothness of the regularization term for the segmentation task, each $J_\eps$ is nonsmooth and is expressed as the sum of an $\ell_1$ regularization term and a smooth nonconvex function. Furthermore, the local minimization of each $J_{\eps}$ is reformulated as the minimization of a smooth function subject to a set of linear constraints. The latter problem is solved by the modified primal-dual interior point method, which guarantees the descent direction at each step. Experimental results are presented and show the effectiveness and the efficiency of the proposed method. Comparison with simulated annealing methods further shows the advantage of our method.

Linearly Constrained Non-Lipschitz Optimization for Image Restoration

A Globally Convergent Algorithm for a Constrained Non-Lipschitz Image Restoration Model

On the Local and Global Minimizers of L(0) Gradient Regularized Model with Box Constraints for Image Restoration

Adaptive Nonconvex Nonsmooth Regularization for Image Restoration Based on Spatial Information

An Iterative Support Shrinking Algorithm for Non-Lipschitz Optimization in Image Restoration

Efficient Convex Optimization for Non-convex Non-smooth Image Restoration

On anisotropic non-lipschitz restoration model: lower bound theory and iterative algorithm

On anisotropic non-Lipschitz restoration model: lower bound theory and convergent algorithm

On a General Smoothly Truncated Regularization for Variational Piecewise Constant Image Restoration: Construction and Convergent Algorithms

Efficient Reconstruction of Piecewise Constant Images Using Nonsmooth Nonconvex Minimization

Single image restoration through l 2 -relaxed truncated l 0 analysis-based sparse optimization in tight frames

On Nonconvex and Nonsmooth Regularization for Contrast-Preserving Variational Signal and Image Restoration

Single Image Restoration Through ℓ2-Relaxed Truncated ℓ0 Analysis-Based Sparse Optimization in Tight Frames

Generalized nonconvex nonsmooth four-directional total variation with overlapping group sparsity for image restoration

A Nonconvex Model with Minimax Concave Penalty for Image Restoration

An Investigation on Semismooth Newton Based Augmented Lagrangian Method for Image Restoration

Tikhonov regularization with conjugate gradient least squares method for large-scale discrete ill-posed problem in image restoration

Image restoration based on Bregman iterative double regularization

Hybrid regularization based restoration algorithm for sparse images

A Novel Variable Exponent Non-Convex TVp,q(x) Model in Image Restoration

Nonlocal Adaptive Biharmonic Regularizer for Image Restoration.