Colour image segmentation by the vector-valued Allen-Cahn phase-field model: a multigrid solution

We propose a new method for the numerical solution of a PDE-driven model for colour image segmentation and give numerical examples of the results. The method combines the vector-valued Allen-Cahn phase field equation with initial data fitting terms. This method is known to be closely related to the Mumford-Shah problem and the level set segmentation by Chan and Vese. Our numerical solution is performed using a multigrid splitting of a finite element space, thereby producing an efficient and robust method for the segmentation of large images.Comment: 17 pages, 9 figure

Functional Liftings of Vectorial Variational Problems with Laplacian Regularization

We propose a functional lifting-based convex relaxation of variational problems with Laplacian-based second-order regularization. The approach rests on ideas from the calibration method as well as from sublabel-accurate continuous multilabeling approaches, and makes these approaches amenable for variational problems with vectorial data and higher-order regularization, as is common in image processing applications. We motivate the approach in the function space setting and prove that, in the special case of absolute Laplacian regularization, it encompasses the discretization-first sublabel-accurate continuous multilabeling approach as a special case. We present a mathematical connection between the lifted and original functional and discuss possible interpretations of minimizers in the lifted function space. Finally, we exemplarily apply the proposed approach to 2D image registration problems.Comment: 12 pages, 3 figures; accepted at the conference "Scale Space and Variational Methods" in Hofgeismar, Germany 201

Discrete stochastic approximations of the Mumford-Shah functional

We propose a $\Gamma$-convergent discrete approximation of the Mumford-Shah functional. The discrete functionals act on functions defined on stationary stochastic lattices and take into account general finite differences through a non-convex potential. In this setting the geometry of the lattice strongly influences the anisotropy of the limit functional. Thus we can use statistically isotropic lattices and stochastic homogenization techniques to approximate the vectorial Mumford-Shah functional in any dimension.Comment: 47 pages, reorganized versio

Sublabel-Accurate Relaxation of Nonconvex Energies

We propose a novel spatially continuous framework for convex relaxations based on functional lifting. Our method can be interpreted as a sublabel-accurate solution to multilabel problems. We show that previously proposed functional lifting methods optimize an energy which is linear between two labels and hence require (often infinitely) many labels for a faithful approximation. In contrast, the proposed formulation is based on a piecewise convex approximation and therefore needs far fewer labels. In comparison to recent MRF-based approaches, our method is formulated in a spatially continuous setting and shows less grid bias. Moreover, in a local sense, our formulation is the tightest possible convex relaxation. It is easy to implement and allows an efficient primal-dual optimization on GPUs. We show the effectiveness of our approach on several computer vision problems

A stochastic-variational model for soft Mumford-Shah segmentation

In contemporary image and vision analysis, stochastic approaches demonstrate great flexibility in representing and modeling complex phenomena, while variational-PDE methods gain enormous computational advantages over Monte-Carlo or other stochastic algorithms. In combination, the two can lead to much more powerful novel models and efficient algorithms. In the current work, we propose a stochastic-variational model for soft (or fuzzy) Mumford-Shah segmentation of mixture image patterns. Unlike the classical hard Mumford-Shah segmentation, the new model allows each pixel to belong to each image pattern with some probability. We show that soft segmentation leads to hard segmentation, and hence is more general. The modeling procedure, mathematical analysis, and computational implementation of the new model are explored in detail, and numerical examples of synthetic and natural images are presented.Comment: 22 page

Phase field approximation of cohesive fracture models

We obtain a cohesive fracture model as a $\Gamma$-limit of scalar damage models in which the elastic coefficient is computed from the damage variable $v$ through a function $f_k$ of the form $f_k(v)=min\{1,\varepsilon_k^{1/2} f(v)\}$, with $f$ diverging for $v$ close to the value describing undamaged material. The resulting fracture energy can be determined by solving a one-dimensional vectorial optimal profile problem. It is linear in the opening $s$ at small values of $s$ and has a finite limit as $s\to\infty$. If the function $f$ is allowed to depend on the index $k$, for specific choices we recover in the limit Dugdale's and Griffith's fracture models, and models with surface energy density having a power-law growth at small openings