New normalized nonlocal hybrid level set method for image segmentation

This article introduces a new normalized nonlocal hybrid level set method for image segmentation. Due to intensity overlapping, blurred edges with complex backgrounds, simple intensity and texture information, such kind of image segmentation is still a challenging task. The proposed method uses both the region and boundary information to achieve accurate segmentation results. The region information can help to identify rough region of interest and prevent the boundary leakage problem. It makes use of normalized nonlocal comparisons between pairs of patches in each region, and a heuristic intensity model is proposed to suppress irrelevant strong edges and constrain the segmentation. The boundary information can help to detect the precise location of the target object, it makes use of the geodesic active contour model to obtain the target boundary. The corresponding variational segmentation problem is implemented by a level set formulation. We use an internal energy term for geometric active contours to penalize the deviation of the level set function from a signed distance function. At last, experimental results on synthetic images and real images are shown in the paper with promising results.     

Using a level set approach for image segmentation under interpolation conditions

In this paper, we propose a new 2D segmentation model including geometric constraints, namely interpolation conditions, to detect objects in a given image. We propose to apply the deformable models to an explicit function using the level set approach (Osher and Sethian [24]); so, we avoid the classical problem of parameterization of both segmentation representation and interpolation conditions. Furthermore, we allow this representation to have topological changes. A problem of energy minimization on a closed subspace of a Hilbert space is defined and introducing Lagrange multipliers enables us to formulate the corresponding variational problem with interpolation conditions. Thus the explicit function evolves, while minimizing the energy and it stops evolving when the desired outlines of the object to detect are reached. The stopping term, as in the classical deformable models, is related to the gradient of the image. Numerical results are given. AMS subject classification 74G65, 46-xx, 92C55     

Image segmentation by level set evolution with region consistency constraint

Image segmentation is a key and fundamental problem in image processing, computer graphics, and computer vision. Level set based method for image segmentation is used widely for its topology flexibility and proper mathematical formulation. However, poor performance of existing level set models on noisy images and weak boundary limit its application in image segmentation. In this paper, we present a region consistency constraint term to measure the regional consistency on both sides of the boundary, this term defines the boundary of the image within a range, and hence increases the stability of the level set model. The term can make existing level set models significantly improve the efficiency of the algorithms on segmenting images with noise and weak boundary. Furthermore, this constraint term can make edge-based level set model overcome the defect of sensitivity to the initial contour. The experimental results show that our algorithm is efficient for image segmentation and outperform the existing state-of-art methods regarding images with noise and weak boundary.     

A variant of the level set method and applications to image segmentation

In this paper we propose a variant of the level set formulation for identifying curves separating regions into different phases. In classical level set approaches, the sign of level set functions are utilized to identify up to phases. The novelty in our approach is to introduce a piecewise constant level set function and use each constant value to represent a unique phase. If phases should be identified, the level set function must approach predetermined constants. We just need one level set function to represent unique phases, and this gains in storage capacity. Further, the reinitializing procedure requested in classical level set methods is superfluous using our approach. The minimization functional for our approach is locally convex and differentiable and thus avoids some of the problems with the nondifferentiability of the Delta and Heaviside functions. Numerical examples are given, and we also compare our method with related approaches.

     

The shape optimization of the arterial graft design by level set methods

The goal of the arterial graft design problem is to find an optimal graft built on an occluded artery, which can be mathematically modeled by a fluid based shape optimization problem. The smoothness of the graft is one of the important aspects in the arterial graft design problem since it affects the flow of the blood significantly. As an attractive design tool for this problem, level set methods are quite efficient for obtaining better shape of the graft. In this paper, a cubic spline level set method and a radial basis function level set method are designed to solve the arterial graft design problem. In both approaches, the shape of the arterial graft is implicitly tracked by the zero-level contour of a level set function and a high level of smoothness of the graft is achieved. Numerical results show the efficiency of the algorithms in the arterial graft design.     

Noninformative priors for the ratio of the shape parameters of two Weibull distributions

The Weibull distribution is one of the most widely used lifetime distributions in reliability engineering. Here, the noninformative priors for the ratio of the shape parameters of two Weibull models are introduced. The first criterion used is the asymptotic matching of the coverage probabilities of Bayesian credible intervals with the corresponding frequentist coverage probabilities. We develop the probability matching priors for the ratio of the shape parameters using the following matching criteria: quantile matching, matching of the distribution function, highest posterior density matching, and matching via inversion of the test statistics. We obtain one particular prior that meets all the matching criteria. Next, we derive the reference priors for different groups of ordering. Our findings show that some of the reference priors satisfy a first-order matching criterion and the one-at-a-time reference prior is a second-order matching prior. Lastly, we perform a simulation study and provide a real-world example.     

An adaptive level set method for shock-driven fluid-structure interaction

The fluid-structure interaction simulation of shock- and detonation-loaded structures requires numerical methods that can cope with large deformations as well as local topology changes. A robust, level-set-based shock-capturing fluid solver is described that allows coupling to any solid mechanics solver. As computational example, the elastic response of a thin steel panel, modeled with both shell and beam theory, to a shock wave in air is considered. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A fast level set model of multi‐atlas labels fusion for 3D MRI tissues segmentation with application of split Bregman method

In this paper, we propose a new fast level set model of multi‐atlas labels fusion for 3D magnetic resonance imaging (MRI) tissues segmentation. The proposed model is aimed at segmenting regions of interest in MR images, especially the tissues such as the amygdala, the caudate, the hippocampus, the pallidum, the putamen, and the thalamus. We first define a new energy functional by taking full advantage of an image data term, a length term, and a label fusion term. Different from using the region‐scalable fitting image data term and length term directly, we define a new image data term and a new length term, which is also incorporated with an edge detect function. By introducing a spatially weight function into the label fusion term, segmentation sensitivity to warped images can be largely improved. Furthermore, the special structure of the new energy functional ensures the application of the split Bregman method, which is a significant highlight and can improve segmentation efficiency of the proposed model. Because of these promotions, several good characters, such as accuracy, efficiency, and robustness have been exhibited in experimental results. Quantitative and qualitative comparisons with other methods have demonstrated the superior advantages of the proposed model.     

Interval-valued level cut sets of fuzzy set

The connections between Zadeh fuzzy set and three-valued fuzzy set are established in this paper. The concepts of interval-valued level cut sets on Zadeh fuzzy set are presented and new decomposition theorems and representation theorems of Zadeh fuzzy set are established based on new cut sets. Firstly, four interval-valued level cut sets on Zadeh fuzzy set are defined as three-valued fuzzy sets and it is shown that the interval-valued level cut sets of Zadeh fuzzy set are generalizations of normal cut sets on Zadeh fuzzy set, and have the same properties as those of normal cut sets of Zadeh fuzzy set. Secondly, the new decomposition theorems are established based on these new cut sets. It is pointed out that each kind of interval-valued level cut sets corresponds to two decomposition theorems. Thus eight decomposition theorems are obtained. Finally, the definitions of three-valued inverse order nested sets and three-valued order nested sets are presented with eight representation theorems based on new nested sets.     

Asymptotic properties of plug-in level set estimators for right censored data

We assume T1,...,Tn are i.i.d.data sampled from distribution function F with density function f and C1,...,Cn are i.i.d.data sampled from distribution function G.Observed data consists of pairs(Xi,δi),i=1,...,n,where Xi=min{Ti,Ci},δi=I(Ti Ci),I(A)denotes the indicator function of the set A.Based on the right censored data{Xi,δi},i=1,...,n,we consider the problem of estimating the level set{f c}of an unknown one-dimensional density function f and study the asymptotic behavior of the plug-in level set estimators.Under some regularity conditions,we establish the asymptotic normality and the exact convergence rate of theλg-measure of the symmetric difference between the level set{f c}and its plug-in estimator{fn c},where f is the density function of F,and fn is a kernel-type density estimator of f.Simulation studies demonstrate that the proposed method is feasible.Illustration with a real data example is also provided.     

A fast multiphase level set algorithm for solving the Chan-Vese model

In this work, we specifically solve the C-V active contour model by multiphase level set methods. We first develop a fast algorithm based on calculating the variational energy of the C-V model without the length term. We check whether the energy decreases or not when we move a point to another segmented region. Then we draw a connection between this algorithm and the topological derivative, a concept emerged from the shape optimization field. Furthermore, to include the length term of the C-V model, a preprocessing step is taken by using nonlinear diffusion. Numerical experiments have demonstrated the efficiency and the robustness of our algorithm. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A level set approach for computing discontinuous solutions of Hamilton-Jacobi equations

We introduce two types of finite difference methods to compute the L-solution and the proper viscosity solution recently proposed by the second author for semi-discontinuous solutions to a class of Hamilton-Jacobi equations. By regarding the graph of the solution as the zero level curve of a continuous function in one dimension higher, we can treat the corresponding level set equation using the viscosity theory introduced by Crandall and Lions. However, we need to pay special attention both analytically and numerically to prevent the zero level curve from overturning so that it can be interpreted as the graph of a function. We demonstrate our Lax-Friedrichs type numerical methods for computing the L-solution using its original level set formulation. In addition, we couple our numerical methods with a singular diffusive term which is essential to computing solutions to a more general class of HJ equations that includes conservation laws. With this singular viscosity, our numerical methods do not require the divergence structure of equations and do apply to more general equations developing shocks other than conservation laws. These numerical methods are generalized to higher order accuracy using weighted ENO local Lax-Friedrichs methods as developed recently by Jiang and Peng. We verify that our numerical solutions approximate the proper viscosity solutions obtained by the second author in a recent Hokkaido University preprint. Finally, since the solution of scalar conservation law equations can be constructed using existing numerical techniques, we use it to verify that our numerical solution approximates the entropy solution.

     

Tokunaga and Horton self-similarity for level set trees of Markov chains

The Horton and Tokunaga branching laws provide a convenient framework for studying self-similarity in random trees. The Horton self-similarity is a weaker property that addresses the principal branching in a tree; it is a counterpart of the power-law size distribution for elements of a branching system. The stronger Tokunaga self-similarity addresses so-called side branching. The Horton and Tokunaga self-similarity have been empirically established in numerous observed and modeled systems, and proven for two paradigmatic models: the critical Galton–Watson branching process with finite progeny and the finite-tree representation of a regular Brownian excursion. This study establishes the Tokunaga and Horton self-similarity for a tree representation of a finite symmetric homogeneous Markov chain. We also extend the concept of Horton and Tokunaga self-similarity to infinite trees and establish self-similarity for an infinite-tree representation of a regular Brownian motion. We conjecture that fractional Brownian motions are also Tokunaga and Horton self-similar, with self-similarity parameters depending on the Hurst exponent.     

A level set algorithm for a class of reverse convex programs

A new algorithm is presented for minimizing a linear function subject to a set of linear inequalities and one additional reverse convex constraint. The algorithm utilizes a conical partition of the convex polytope in conjuction with its facets in order to remain on the level surface of the reverse convex constraint. The algorithm does not need to solve linear programs on a set of cones which converges to a line segment.     

Robust priors for smoothing and image restoration

The Bayesian method for restoring an image corrupted by added Gaussian noise uses a Gibbs prior for the unknown clean image. The potential of this Gibbs prior penalizes differences between adjacent grey levels. In this paper we discuss the choice of the form and the parameters of the penalizing potential in a particular example used previously by Ogata (1990,Ann. Inst. Statist. Math.,42, 403–433). In this example the clean image is piecewise constant, but the constant patches and the step sizes at edges are small compared with the noise variance. We find that contrary to results reported in Ogata (1990,Ann. Inst. Statist. Math.,42, 403–433) the Bayesian method performs well provided the potential increases more slowly than a quadratic one and the scale parameter of the potential is sufficiently small. Convex potentials with bounded derivatives perform not much worse than bounded potentials, but are computationally much simpler. For bounded potentials we use a variant of simulated annealing. For quadratic potentials data-driven choices of the smoothing parameter are reviewed and compared. For other potentials the smoothing parameter is determined by considering which deviations from a flat image we would like to smooth out and retain respectively.     

Applications of level set methods in computational biophysics

We describe in this paper two applications of Eulerian level set methods to fluid-structure problems arising in biophysics. The first one is concerned with three-dimensional equilibrium shapes of phospholipidic vesicles. This is a complex problem, which can be recast as the minimization of the curvature energy of an immersed elastic membrane, under a constant area constraint. The second deals with isolated cardiomyocyte contraction. This problem corresponds to a generic incompressible fluid-structure coupling between an elastic body and a fluid. By the choice of these two quite different situations, we aim to bring evidence that Eulerian methods provide efficient and flexible computational tools in biophysics applications.     

Estimation of level set trees using adaptive partitions

We present methods for the estimation of level sets, a level set tree, and a volume function of a multivariate density function. The methods are such that the computation is feasible and estimation is statistically efficient in moderate dimensional cases (\(d\approx 8\)) and for moderate sample sizes (\(n\approx \) 50,000). We apply kernel estimation together with an adaptive partition of the sample space. We illustrate how level set trees can be applied in cluster analysis and in flow cytometry.