Adaptive Segmentation Model for Images with Intensity Inhomogeneity based on Local Neighborhood Contrast

Segmentation of images with intensity inhomogeneity is a significant task in the field of image processing, especially in medical image processing and analysis. Some local region-based models work well on handling intensity inhomogeneity, but they are always sensitive to contour initialization and high noise. In this paper, we present an adaptive segmentation model for images with intensity inhomogeneity in the form of partial differential equation. Firstly, a global intensity fitting term and a local intensity fitting term are constructed by employing the global and local image information, respectively. Secondly, a tradeoff function is defined to adjust adaptively the weight between two fitting terms, which is based on the neighborhood contrast of image pixel. Finally, a weighted regularization term related to local entropy is used to ensure the smoothness of evolution curve. Meanwhile, a distance regularization term is added for stable level set evolution. Experimental results show that the proposed model without initial contour can segment inhomogeneous images stably and effectively, which thereby avoiding the influence of contour initialization on segmentation results. Besides, the proposed model works better on noise images comparing with two relevant segmentation models.     

On a Variational Model for Selective Image Segmentation of Features with Infinite Perimeter

Variational models provide reliable formulation for segmentation of features and their boundaries in an image, following the seminal work of Mumford-Shah (1989, Commun. Pure Appl. Math.) on dividing a general surface into piecewise smooth sub-surfaces. A central idea of models based on this work is to minimize the length of feature’s boundaries (i.e., H1 Hausdorff measure). However there exist problems with irregular and oscillatory object boundaries, where minimizing such a length is not appropriate, as noted by Barchiesi et al. (2010, SIAM J. Multiscale Model. Simu.) who proposed to miminize L2 Lebesgue measure of the γ-neighborhood of the boundaries. This paper presents a dual level set selective segmentation model based on Barchiesi et al. (2010) to automatically select a local feature instead of all global features. Our model uses two level set functions: a global level set which segments all boundaries, and the local level set which evolves and finds the boundary of the object closest to the geometric constraints. Using real life images with oscillatory boundaries, we show qualitative results demonstrating the effectiveness of the proposed method.     

A variational formulation for physical noised image segmentation

Image segmentation is a hot topic in image science. In this paper we present a new variational segmentation model based on the theory of Mumford-Shah model. The aim of our model is to divide noised image, according to a certain criterion, into homogeneous and smooth regions that should correspond to structural units in the scene or objects of interest. The proposed region-based model uses total variation as a regularization term, and different fidelity term can be used for image segmentation in the cases of physical noise, such as Gaussian, Poisson and multiplicative speckle noise. Our model consists of five weighted terms, two of them are responsible for image denoising based on fidelity term and total variation term, the others assure that the three conditions of adherence to the data, smoothing, and discontinuity detection are met at once. We also develop a primal-dual hybrid gradient algorithm for our model. Numerical results on various synthetic and real images are provided to compare our method with others,these results show that our proposed model and algorithms are effective.     

Generalized row-action methods for tomographic imaging

Row-action methods play an important role in tomographic image reconstruction. Many such methods can be viewed as incremental gradient methods for minimizing a sum of a large number of convex functions, and despite their relatively poor global rate of convergence, these methods often exhibit fast initial convergence which is desirable in applications where a low-accuracy solution is acceptable. In this paper, we propose relaxed variants of a class of incremental proximal gradient methods, and these variants generalize many existing row-action methods for tomographic imaging. Moreover, they allow us to derive new incremental algorithms for tomographic imaging that incorporate different types of prior information via regularization. We demonstrate the efficacy of the approach with some numerical examples.     

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.     

A new family of conjugate gradient methods

In this paper we develop a new class of conjugate gradient methods for unconstrained optimization problems. A new nonmonotone line search technique is proposed to guarantee the global convergence of these conjugate gradient methods under some mild conditions. In particular, Polak–Ribiére–Polyak and Liu–Storey conjugate gradient methods are special cases of the new class of conjugate gradient methods. By estimating the local Lipschitz constant of the derivative of objective functions, we can find an adequate step size and substantially decrease the function evaluations at each iteration. Numerical results show that these new conjugate gradient methods are effective in minimizing large-scale non-convex non-quadratic functions.     

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.     

Split Bregman method for minimization of improved active contour model combining local and global information dynamically

This paper presents an improved active contour model by combining the Chan–Vese model, the region-scalable fitting energy model, the globally convex segmentation method and the split Bregman method. A weight function that varies with the location of a given image is used to control the influence of the local and global information dynamically. We first present our model in a 2-phase level set formulation and then extend it to a multi-phase formulation. By taking the local and global information into consideration together, our model can segment more general images, especially images with intensity inhomogeneity. Our model has been applied to synthetic and real images with promising results. Numerical results show the advantages of our model compared with other models. The accuracy and efficiency are demonstrated by the numerical results. Besides, our model is robust in the presence of noise.     

Statistical shape priors for level set segmentation

Starting in the early 1990's level set methods have become a popular mathematical framework for variational image segmentation. In many applications of segmentation, however, cost functionals which merely take into account the intensity information of the input image will not give rise to the desired segmentation results. To cope with missing or misleading image information, researchers have proposed to impose statistical shape priors into the segmentation process. Such shape priors favor the evolving embedding function to remain similar to embedding functions associated with a collection of training shapes. As a consequence, one can obtain shape-consistent segmentation despite large amounts of noise, background clutter and partial occlusion of the object of interest. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A New Implicit Method for Surface Segmentation by Minimal Paths in 3D Images

We introduce a novel implicit approach for single-object segmentation in 3D images. The boundary surface of this object is assumed to contain two known curves (the constraining curves), given by an expert. The aim of our method is to find the wanted surface by exploiting as much as possible the information given in the supplied curves and in the image. As for active surfaces, we use a cost potential that penalizes image regions of low interest (most likely areas of low gradient or too far from the surface to be extracted). In order to avoid local minima, we introduce a new partial differential equation and use its solution for segmentation. We show that the zero level set of this solution contains the constraining curves as well as a set of paths joining them. We present a fast implementation that has been successfully applied to 3D medical and synthetic images.     

An efficient algorithm for large scale global optimization of continuous functions

A fast descent algorithm, resorting to a “stretching” function technique and built on one hybrid method (GRSA) which combines simulated annealing (SA) algorithm and gradient based methods for large scale global optimizations, is proposed. Unlike the previously proposed method in which the original objective functions remain unchanged during the whole course of optimization, the new method firstly constructs an auxiliary function on one local minimizer obtained by gradient based methods and then SA is executed on this constructed auxiliary function instead of on the original objective function in order that we can improve the jumping ability of SA algorithm to escape from the currently discovered local minimum to a better one from which the gradient based methods restart a new local search. The above procedure is repeated until a global minimum is detected. In addition, corresponding to the adopted “stretching” technique, a new next trial point generating scheme is designed. It is verified by simulation especially on large scale problems that the convergence speed is greatly accelerated, which is its main difference from many other reported methods that mostly cope with functions with less than 50 variables and does not apply to large scale optimization problems. Furthermore, the new algorithm functions as a global optimization procedure with a high success probability and high solution precision.     

The Topological Gradient Method: From Optimal Design to Image Processing

The aim of this article is to review and extend the applications of the topological gradient to major image processing problems. We briefly review the topological gradient, and then present its application to the crack localization problem, which can be solved using the Dirichlet to Neumann approach. A very natural application of this technique in image processing is the inpainting problem, which can be solved by identifying the optimal location of the missing edges. Edge detection is of extreme importance, as edges convey essential information in a picture. A second natural application is then the image reconstruction. A class of image reconstruction problems is considered that includes restoration, demosaicing, segmentation and super-resolution. These problems are studied using a unified theoretical framework which is based on the topological gradient method. This tool is able to find the localization and orientation of the edges for blurred, low sampled, partially masked, noisy images. We review existing algorithms and propose new ones. The performance of our approach is compared with conventional image reconstruction processes.     

基于小波与齐次Besov空间的图像分割算法

从具有全局最优解的几何活动轮廓方法出发,分别提出了两种基于齐次Besov窄间与小波变换的图像分割算法,并给出了解的存在性证明.数值求解利用小波软阈值以及分裂Bregman方法,能够有效提高计算效率.由于小波变换具有多分辨特性,对于包含较多细节信息的图像,采用新算法能够得到更好的分割效果.数值实验表明采用新算法能够获得较...     

Preconditioned Iterative Methods for Algebraic Systems from Multiplicative Half-Quadratic Regularization Image Restorations

<正>Image restoration is often solved by minimizing an energy function consisting of a data-fidelity term and a regularization term.A regularized convex term can usually preserve the image edges well in the restored image.In this paper,we consider a class of convex and edge-preserving regularization functions,i.e.,multiplicative half-quadratic regularizations,and we use the Newton method to solve the correspondingly reduced systems of nonlinear equations.At each Newton iterate,the preconditioned conjugate gradient method,incorporated with a constraint preconditioner,is employed to solve the structured Newton equation that has a symmetric positive definite coefficient matrix. The eigenvalue bounds of the preconditioned matrix are deliberately derived,which can be used to estimate the convergence speed of the preconditioned conjugate gradient method.We use experimental results to demonstrate that this new approach is efficient, and the effect of image restoration is reasonably well.     

Noisy image segmentation based on nonlinear diffusion equation model

This paper addresses the segmentation problem in noisy image based on nonlinear diffusion equation model and proposes a new adaptive segmentation model based on gray-level image segmentation model. This model also can be extended to the vector value image segmentation. By virtue of the prior information of regions and boundary of image, a framework is established to construct different segmentation models using different probability density functions. A segmentation model exploiting Gauss probability density function is given in this paper. An efficient and unconditional stable algorithm based on locally one-dimensional (LOD) scheme is developed and it is used to segment the gray image and the vector values image. Comparing with existing classical models, the proposed approach gives the best performance.     

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.     

Existence and regularity of minimizers of a functional for unsupervised multiphase segmentation

We consider a variational model for image segmentation proposed in Sandberg et al. (2010) [12]. In such a model the image domain is partitioned into a finite collection of subsets denoted as phases. The segmentation is unsupervised, i.e., the model finds automatically an optimal number of phases, which are not required to be connected subsets. Unsupervised segmentation is obtained by minimizing a functional of the Mumford–Shah type (Mumford and Shah, 1989 [1]), but modifying the geometric part of the Mumford–Shah energy with the introduction of a suitable scale term. The results of computer experiments discussed in [12] show that the resulting variational model has several properties which are relevant for applications. In this paper we investigate the theoretical properties of the model. We study the existence of minimizers of the corresponding functional, first looking for a weak solution in a class of phases constituted by sets of finite perimeter. Then we find various regularity properties of such minimizers, particularly we study the structure of triple junctions by determining their optimal angles.     

An improved Perry conjugate gradient method with adaptive parameter choice

In this paper, we present a conjugate gradient method for solving unconstrained optimization problems. Motivated by Perry conjugate gradient method and Dai-Liao method, an improved Perry update matrix is proposed to overcome the non-symmetric positive definite property of the Perry matrix. The parameter in the update matrix is determined by minimizing the condition number of the iterative matrix which can ensure the positive definite property. The obtained method can also be considered as a modified form of CG-DESCENT method with an adjusted term. Under some mild conditions, the presented method is global convergent. Numerical experiments under CUTEst environment show that the proposed algorithm is promising.     

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