85 research outputs found
Geometrical-based algorithm for variational segmentation and smoothing of vector-valued images
An optimisation method based on a nonlinear functional is considered for segmentation and smoothing of vector-valued images. An edge-based approach is proposed to initially segment the image using geometrical properties such as metric tensor of the linearly smoothed image. The nonlinear functional is then minimised for each segmented region to yield the smoothed image. The functional is characterised with a unique solution in contrast with the MumfordâShah functional for vector-valued images. An operator for edge detection is introduced as a result of this unique solution. This operator is analytically calculated and its detection performance and localisation are then compared with those of the DroGoperator. The implementations are applied on colour images as examples of vector-valued images, and the results demonstrate robust performance in noisy environments
Piecewise smooth reconstruction of normal vector field on digital data
International audienceWe propose a novel method to regularize a normal vector field defined on a digital surface (boundary of a set of voxels). When the digital surface is a digitization of a piecewise smooth manifold, our method localizes sharp features (edges) while regularizing the input normal vector field at the same time. It relies on the optimisation of a variant of the Ambrosio-Tortorelli functional, originally defined for denoising and contour extraction in image processing [AT90]. We reformulate this functional to digital surface processing thanks to discrete calculus operators. Experiments show that the output normal field is very robust to digitization artifacts or noise, and also fairly independent of the sampling resolution. The method allows the user to choose independently the amount of smoothing and the length of the set of discontinuities. Sharp and vanishing features are correctly delineated even on extremely damaged data. Finally, our method can be used to enhance considerably the output of state-of- the-art normal field estimators like Voronoi Covariance Measure [MOG11] or Randomized Hough Transform [BM12]
A new proximal method for joint image restoration and edge detection with the Mumford-Shah model
International audienceIn this paper, we propose an adaptation of the PAM algorithm to the minimization of a nonconvex functional designed for joint image denoising and contour detection. This new functional is based on the AmbrosioâTortorelli approximation of the well-known MumfordâShah functional. We motivate the proposed approximation, offering flexibility in the choice of the possibly non-smooth penalization, and we derive closed form expression for the proximal steps involved in the algorithm. We focus our attention on two types of penalization: 1-norm and a proposed quadratic-1 function. Numerical experiments show that the proposed method is able to detect sharp contours and to reconstruct piecewise smooth approximations with low computational cost and convergence guarantees. We also compare the results with state-of-the-art re-laxations of the MumfordâShah functional and a recent discrete formulation of the AmbrosioâTortorelli functional
A discrete approximation of Blake & Zisserman energy in image denoising and optimal choice of regularization parameters
We consider a multi-scale approach for the discrete approximation of a functional proposed by Bake and Zisserman (BZ) for solving image denoising and segmentation problems. The proposed method is based on simple and effective higher order varia-tional model. It consists of building linear discrete energies family which Î-converges to the non-linear BZ functional. The key point of the approach is the construction of the diffusion operators in the discrete energies within a finite element adaptive procedure which approximate in the Î-convergence sense the initial energy including the singular parts. The resulting model preserves the singularities of the image and of its gradient while keeping a simple structure of the underlying PDEs, hence efficient numerical method for solving the problem under consideration. A new point to make this approach work is to deal with constrained optimization problems that we circumvent through a Lagrangian formulation. We present some numerical experiments to show that the proposed approach allows us to detect first and second-order singularities. We also consider and implement to enhance the algorithms and convergence properties, an augmented Lagrangian method using the alternating direction method of Multipliers (ADMM)
Uncertainty quantification in image segmentation using the Ambrosio--Tortorelli approximation of the Mumford--Shah energy
The quantification of uncertainties in image segmentation based on the Mumford-Shah model is studied. The aim is to address the error propagation of noise and other error types in the original image to the restoration result and especially the reconstructed edges (sharp image contrasts). Analytically, we rely on the Ambrosio-Tortorelli approximation and discuss the existence of measurable selections of its solutions as well as sampling-based methods and the limitations of other popular methods. Numerical examples illustrate the theoretical findings
Variational methods and its applications to computer vision
Many computer vision applications such as image segmentation can be formulated in a ''variational'' way as energy minimization problems. Unfortunately, the computational task of minimizing these energies is usually difficult as it generally involves non convex functions in a space with thousands of dimensions and often the associated combinatorial problems are NP-hard to solve. Furthermore, they are ill-posed inverse problems and therefore are extremely sensitive to perturbations (e.g. noise). For this reason in order to compute a physically reliable approximation from given noisy data, it is necessary to incorporate into the mathematical model appropriate regularizations that require complex computations.
The main aim of this work is to describe variational segmentation methods that are particularly effective for curvilinear structures. Due to their complex geometry, classical regularization techniques cannot be adopted because they lead to the loss of most of low contrasted details. In contrast, the proposed method not only better preserves curvilinear structures, but also reconnects some parts that may have been disconnected by noise. Moreover, it can be easily extensible to graphs and successfully applied to different types of data such as medical imagery (i.e. vessels, hearth coronaries etc), material samples (i.e. concrete) and satellite signals (i.e. streets, rivers etc.). In particular, we will show results and performances about an implementation targeting new generation of High Performance Computing (HPC) architectures where different types of coprocessors cooperate. The involved dataset consists of approximately 200 images of cracks, captured in three different tunnels by a robotic machine designed for the European ROBO-SPECT project.Open Acces
Random finite-difference discretizations of the Ambrosio-Tortorelli functional with optimal mesh size
We propose and analyze a finite-difference discretization of the
Ambrosio-Tortorelli functional. It is known that if the discretization is made
with respect to an underlying periodic lattice of spacing , the
discretized functionals -converge to the Mumford-Shah functional only
if , being the elliptic approximation
parameter of the Ambrosio-Tortorelli functional. Discretizing with respect to
stationary, ergodic and isotropic random lattices we prove this
-convergence result also for , a regime at which
the discretization with respect to a periodic lattice converges instead to an
anisotropic version of the Mumford-Shah functional.Comment: 36 pages, 6 figures. Added some numerical example
Joint Image Reconstruction and Segmentation Using the Potts Model
We propose a new algorithmic approach to the non-smooth and non-convex Potts
problem (also called piecewise-constant Mumford-Shah problem) for inverse
imaging problems. We derive a suitable splitting into specific subproblems that
can all be solved efficiently. Our method does not require a priori knowledge
on the gray levels nor on the number of segments of the reconstruction.
Further, it avoids anisotropic artifacts such as geometric staircasing. We
demonstrate the suitability of our method for joint image reconstruction and
segmentation. We focus on Radon data, where we in particular consider limited
data situations. For instance, our method is able to recover all segments of
the Shepp-Logan phantom from angular views only. We illustrate the
practical applicability on a real PET dataset. As further applications, we
consider spherical Radon data as well as blurred data
- âŠ