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.

Poissonian Image Restoration Via the -Based Minimization

Discrete Total Variation Model with Gradient Fidelity Term for Image Restoration

Fast roughness minimizing image restoration under mixed Poisson-Gaussian noise

Restoring Poissonian Images By A Combined First-Order And Second-Order Variation Approach

Image restoration based on the minimax-concave and the overlapping group sparsity

Poisson noisy image restoration via overlapping group sparse and nonconvex second-order total variation priors

High-quality image restoration from partial random samples in spatial domain.

A Nonconvex Model with Minimax Concave Penalty for Image Restoration

ℓ_0 Minimization for Wavelet Frame Based Image Restoration

Novel Deconvolution Approach to Image Restoration

3D Poissonian image deblurring via patch-based tensor logarithmic Schatten- p minimization

Image Smoothing Via Gradient Sparsity and Surface Area Minimization

Efficient Reconstruction of Piecewise Constant Images Using Nonsmooth Nonconvex Minimization

Poisson Image Denoising Using Best Linear Prediction: A Post-processing Framework

A Dual Model for Restoring Images Corrupted by Mixture of Additive and Multiplicative Noise

Deep Block Proximal Linearised Minimisation Algorithm for Non-convex Inverse Problems

Blind multi-Poissonian image deconvolution with sparse log-step gradient prior

Nonconvex Regularization Based Sparse Recovery and Demixing with Application to Color Image Inpainting

A HYBRID ALTERNATING MINIMIZATION ALGORITHM FOR STRUCTURED CONVEX OPTIMIZATION PROBLEMS WITH APPLICATION IN POISSONIAN IMAGE PROCESSING

Hybrid Regularization Image Deblurring in the Presence of Impulsive Noise

An Efficient Inexact Gauss–Seidel-Based Algorithm for Image Restoration with Mixed Noise