4 research outputs found
Convergence Rates for Inverse Problems with Impulsive Noise
We study inverse problems F(f) = g with perturbed right hand side g^{obs}
corrupted by so-called impulsive noise, i.e. noise which is concentrated on a
small subset of the domain of definition of g. It is well known that
Tikhonov-type regularization with an L^1 data fidelity term yields
significantly more accurate results than Tikhonov regularization with classical
L^2 data fidelity terms for this type of noise. The purpose of this paper is to
provide a convergence analysis explaining this remarkable difference in
accuracy. Our error estimates significantly improve previous error estimates
for Tikhonov regularization with L^1-fidelity term in the case of impulsive
noise. We present numerical results which are in good agreement with the
predictions of our analysis
Convergence Rates for Exponentially Ill-Posed Inverse Problems with Impulsive Noise
This paper is concerned with exponentially ill-posed operator equations with
additive impulsive noise on the right hand side, i.e. the noise is large on a
small part of the domain and small or zero outside. It is well known that
Tikhonov regularization with an data fidelity term outperforms Tikhonov
regularization with an fidelity term in this case. This effect has
recently been explained and quantified for the case of finitely smoothing
operators. Here we extend this analysis to the case of infinitely smoothing
forward operators under standard Sobolev smoothness assumptions on the
solution, i.e. exponentially ill-posed inverse problems. It turns out that high
order polynomial rates of convergence in the size of the support of large noise
can be achieved rather than the poor logarithmic convergence rates typical for
exponentially ill-posed problems. The main tools of our analysis are Banach
spaces of analytic functions and interpolation-type inequalities for such
spaces. We discuss two examples, the (periodic) backwards heat equation and an
inverse problem in gradiometry.Comment: to appear in SIAM J. Numer. Ana
Efficient Reconstruction of Piecewise Constant Images Using Nonsmooth Nonconvex Minimization
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ε}, 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ε is nonsmooth and is expressed as the sum of an l1 regularization term and
a smooth nonconvex function. Furthermore, the local minimization of each Jε 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.published_or_final_versio