274 research outputs found
A Posteriori Error Control for the Binary Mumford-Shah Model
The binary Mumford-Shah model is a widespread tool for image segmentation and
can be considered as a basic model in shape optimization with a broad range of
applications in computer vision, ranging from basic segmentation and labeling
to object reconstruction. This paper presents robust a posteriori error
estimates for a natural error quantity, namely the area of the non properly
segmented region. To this end, a suitable strictly convex and non-constrained
relaxation of the originally non-convex functional is investigated and Repin's
functional approach for a posteriori error estimation is used to control the
numerical error for the relaxed problem in the -norm. In combination with
a suitable cut out argument, a fully practical estimate for the area mismatch
is derived. This estimate is incorporated in an adaptive meshing strategy. Two
different adaptive primal-dual finite element schemes, and the most frequently
used finite difference discretization are investigated and compared. Numerical
experiments show qualitative and quantitative properties of the estimates and
demonstrate their usefulness in practical applications.Comment: 18 pages, 7 figures, 1 tabl
A Total Fractional-Order Variation Model for Image Restoration with Non-homogeneous Boundary Conditions and its Numerical Solution
To overcome the weakness of a total variation based model for image
restoration, various high order (typically second order) regularization models
have been proposed and studied recently. In this paper we analyze and test a
fractional-order derivative based total -order variation model, which
can outperform the currently popular high order regularization models. There
exist several previous works using total -order variations for image
restoration; however first no analysis is done yet and second all tested
formulations, differing from each other, utilize the zero Dirichlet boundary
conditions which are not realistic (while non-zero boundary conditions violate
definitions of fractional-order derivatives). This paper first reviews some
results of fractional-order derivatives and then analyzes the theoretical
properties of the proposed total -order variational model rigorously.
It then develops four algorithms for solving the variational problem, one based
on the variational Split-Bregman idea and three based on direct solution of the
discretise-optimization problem. Numerical experiments show that, in terms of
restoration quality and solution efficiency, the proposed model can produce
highly competitive results, for smooth images, to two established high order
models: the mean curvature and the total generalized variation.Comment: 26 page
Active Mean Fields for Probabilistic Image Segmentation: Connections with Chan-Vese and Rudin-Osher-Fatemi Models
Segmentation is a fundamental task for extracting semantically meaningful
regions from an image. The goal of segmentation algorithms is to accurately
assign object labels to each image location. However, image-noise, shortcomings
of algorithms, and image ambiguities cause uncertainty in label assignment.
Estimating the uncertainty in label assignment is important in multiple
application domains, such as segmenting tumors from medical images for
radiation treatment planning. One way to estimate these uncertainties is
through the computation of posteriors of Bayesian models, which is
computationally prohibitive for many practical applications. On the other hand,
most computationally efficient methods fail to estimate label uncertainty. We
therefore propose in this paper the Active Mean Fields (AMF) approach, a
technique based on Bayesian modeling that uses a mean-field approximation to
efficiently compute a segmentation and its corresponding uncertainty. Based on
a variational formulation, the resulting convex model combines any
label-likelihood measure with a prior on the length of the segmentation
boundary. A specific implementation of that model is the Chan-Vese segmentation
model (CV), in which the binary segmentation task is defined by a Gaussian
likelihood and a prior regularizing the length of the segmentation boundary.
Furthermore, the Euler-Lagrange equations derived from the AMF model are
equivalent to those of the popular Rudin-Osher-Fatemi (ROF) model for image
denoising. Solutions to the AMF model can thus be implemented by directly
utilizing highly-efficient ROF solvers on log-likelihood ratio fields. We
qualitatively assess the approach on synthetic data as well as on real natural
and medical images. For a quantitative evaluation, we apply our approach to the
icgbench dataset
Combinatorial Continuous Maximal Flows
Maximum flow (and minimum cut) algorithms have had a strong impact on
computer vision. In particular, graph cuts algorithms provide a mechanism for
the discrete optimization of an energy functional which has been used in a
variety of applications such as image segmentation, stereo, image stitching and
texture synthesis. Algorithms based on the classical formulation of max-flow
defined on a graph are known to exhibit metrication artefacts in the solution.
Therefore, a recent trend has been to instead employ a spatially continuous
maximum flow (or the dual min-cut problem) in these same applications to
produce solutions with no metrication errors. However, known fast continuous
max-flow algorithms have no stopping criteria or have not been proved to
converge. In this work, we revisit the continuous max-flow problem and show
that the analogous discrete formulation is different from the classical
max-flow problem. We then apply an appropriate combinatorial optimization
technique to this combinatorial continuous max-flow CCMF problem to find a
null-divergence solution that exhibits no metrication artefacts and may be
solved exactly by a fast, efficient algorithm with provable convergence.
Finally, by exhibiting the dual problem of our CCMF formulation, we clarify the
fact, already proved by Nozawa in the continuous setting, that the max-flow and
the total variation problems are not always equivalent.Comment: 26 page
Recommended from our members
A Variational Model for Joint Motion Estimation and Image Reconstruction
The aim of this paper is to derive and analyze a variational model for the joint estimation of motion and reconstruction of image sequences, which is based on a time-continuous Eulerian motion model. The model can be set up in terms of the continuity equation or the brightness constancy equation. The analysis in this paper focuses on the latter for robust motion estimation on sequences of twodimensional images. We rigorously prove the existence of a minimizer in a suitable function space setting. Moreover, we discuss the numerical solution of the model based on primal-dual algorithms and investigate several examples. Finally, the benefits of our model compared to existing techniques, such as sequential image reconstruction and motion estimation, are shown.The work of the first author was also supported by the German
Science Foundation DFG via EXC 1003 Cells in Motion Cluster of Excellence, M¨unster, German
Total variation regularization for manifold-valued data
We consider total variation minimization for manifold valued data. We propose
a cyclic proximal point algorithm and a parallel proximal point algorithm to
minimize TV functionals with -type data terms in the manifold case.
These algorithms are based on iterative geodesic averaging which makes them
easily applicable to a large class of data manifolds. As an application, we
consider denoising images which take their values in a manifold. We apply our
algorithms to diffusion tensor images, interferometric SAR images as well as
sphere and cylinder valued images. For the class of Cartan-Hadamard manifolds
(which includes the data space in diffusion tensor imaging) we show the
convergence of the proposed TV minimizing algorithms to a global minimizer
A class of second-order geometric quasilinear hyperbolic PDEs and their application in imaging science
In this paper, we study damped second-order dynamics, which are quasilinear
hyperbolic partial differential equations (PDEs). This is inspired by the
recent development of second-order damping systems for accelerating energy
decay of gradient flows. We concentrate on two equations: one is a damped
second-order total variation flow, which is primarily motivated by the
application of image denoising; the other is a damped second-order mean
curvature flow for level sets of scalar functions, which is related to a
non-convex variational model capable of correcting displacement errors in image
data (e.g. dejittering). For the former equation, we prove the existence and
uniqueness of the solution. For the latter, we draw a connection between the
equation and some second-order geometric PDEs evolving the hypersurfaces which
are described by level sets of scalar functions, and show the existence and
uniqueness of the solution for a regularized version of the equation. The
latter is used in our algorithmic development. A general algorithm for
numerical discretization of the two nonlinear PDEs is proposed and analyzed.
Its efficiency is demonstrated by various numerical examples, where simulations
on the behavior of solutions of the new equations and comparisons with
first-order flows are also documented
- …