3D Multiphase Piecewise Constant Level Set Method Based on Graph Cut Minimization

Segmentation of three-dimensional (3D) complicated structures is of great importance for many real applications. In this work we combine graph cut minimization method with a variant of the level set idea for 3D segmentation based on the Mumford-Shah model. Compared with the traditional approach for solving the Euler-Lagrange equation we do not need to solve any partial differential equations. Instead, the minimum cut on a special designed graph need to be computed. The method is tested on data with complicated structures. It is rather stable with respect to initial value and the algorithm is nearly parameter free. Experiments show that it can solve large problems much faster than traditional approaches.     

基于直方图初始化的有效的多区域Mumford-Shah图像分割方法（英文）

The Mumford-Shah energy functional is a successful image segmentation model. It is a non-convex variational problem and lacks of good initialization techniques so far. In this paper, motivated by the fact that image histogram is a combination of several Gaussian distributions, and their centers can be considered as approximations of cluster centers, we introduce a histogram-based initialization method to compute the cluster centers. With this technique, we then devise an effective multi-region Mumford-Shah image segmentation method, and adopt the recent proximal alternating minimization method to solve the minimization problem. Experiments indicate that our histogram initialization method is more robust than existing methods,and our segmentation method is very effective for both gray and color images.     

A primal-dual active-set algorithm for bilaterally constrained total variation deblurring and piecewise constant Mumford-Shah segmentation problems

In this paper, we propose a fast primal-dual algorithm for solving bilaterally constrained total variation minimization problems which subsume the bilaterally constrained total variation image deblurring model and the two-phase piecewise constant Mumford-Shah image segmentation model. The presence of the bilateral constraints makes the optimality conditions of the primal-dual problem semi-smooth which can be solved by a semi-smooth Newton’s method superlinearly. But the linear system to solve at each iteration is very large and difficult to precondition. Using a primal-dual active-set strategy, we reduce the linear system to a much smaller and better structured one so that it can be solved efficiently by conjugate gradient with an approximate inverse preconditioner. Locally superlinear convergence results are derived for the proposed algorithm. Numerical experiments are also provided for both deblurring and segmentation problems. In particular, for the deblurring problem, we show that the addition of the bilateral constraints to the total variation model improves the quality of the solutions.     

Simultaneous reconstruction and segmentation for tomography data

We present a Mumford-Shah like approach for the inversion of CT and SPECT-data (Single Photon Emission Computerized Tomography). With this approach we aim at the simultaneous reconstruction and segmentation of activity and density distribution from given tomography data. We assume the functions to be piecewise constant with respect to a set of contours. Shape sensitivity analysis is used to find a descent direction for the cost functional which leads to an update formula for the contour in a level set framework. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Retinex modulated piecewise constant variational model for image segmentation and bias correction

In this paper, we propose a novel Retinex induced piecewise constant variational model for simultaneous segmentation of images with intensity inhomogeneity and bias correction. Firstly, we obtain an additive model by decomposing the original image into a smooth bias component and a structure part based on the Retinex theory. Secondly, the structure part can be modeled by the piecewise constant variational model and thus deduced a new data fidelity term. Finally, we formulate a new energy functional by incorporating the data fidelity term into the level set framework and introducing a GL-regularizer to the level set function and a smooth regularizer to model the bias component. Based on the alternating minimization algorithm and the operator splitting method, we present a numerical scheme to solve the minimization problem efficiently. Experimental results on images from diverse modalities demonstrate the competitive performances of the proposed model and algorithm over other representative methods in term of efficiency and robustness.     

求解分段常数图象分割模型的一个快速算法

针对Xue-ChengTai等提出的分段常数图象分割模型,我们提出了一个新的快速求解算法。通过引进一个函数来选择模型中的正则化参数β的值,并判断在迭代过程中何时求解不含惩罚项的泛函F。此函数的引入有效地加速了算法的收敛速度。结合原始-对偶Newton方法来求解总变差最小化问题。数值试验表明新算法具有很快的收敛速度与良好的分割效果,且算法对初始值的要求不高。     

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.     

Two-Stage Image Segmentation Scheme Based on Inexact Alternating Direction Method

Image segmentation is a fundamental problem in both image processing and computer vision with numerous applications. In this paper, we propose a two-stage image segmentation scheme based on inexact alternating direction method. Specifically, we first solve the convex variant of the Mumford-Shah model to get the smooth solution, and the segmentation is then obtained by applying the K-means clustering method to the solution. Some numerical comparisons are arranged to show the effectiveness of our proposed schemes by segmenting many kinds of images such as artificial images, natural images, and brain MRI images.     

Two-parallel machines scheduling with rate-modifying activities to minimize total completion time

This paper considers the two-parallel machines scheduling problem with rate-modifying activities. In this model, each machine has a rate-modifying activity that can change the processing rate of machine under consideration. Hence the actual processing times of jobs vary depending on whether the job is scheduled before or after the rate-modifying activity. We need to make a decision on when to schedule the rate-modifying activities and the sequence of jobs to minimize some objective function. We provide polynomial and pseudo-polynomial time algorithms to solve the total completion time minimization problem and total weighted completion time minimization problem under agreeable ratio condition.     

Shape and topology optimization for elliptic boundary value problems using a piecewise constant level set method

The aim of this paper is to propose a variational piecewise constant level set method for solving elliptic shape and topology optimization problems. The original model is approximated by a two-phase optimal shape design problem by the ersatz material approach. Under the piecewise constant level set framework, we first reformulate the two-phase design problem to be a new constrained optimization problem with respect to the piecewise constant level set function. Then we solve it by the projection Lagrangian method. A gradient-type iterative algorithm is presented. Comparisons between our numerical results and those obtained by level set approaches show the effectiveness, accuracy and efficiency of our algorithm.     

Parallel machines scheduling with a deteriorating maintenance activity

In this paper, we consider parallel identical machines scheduling problems with a deteriorating maintenance activity. In this model, each machine has a deteriorating maintenance activity, that is, delaying the maintenance increases the time required to perform it. We need to make a decision on when to schedule the deteriorating maintenance activities and the sequence of jobs to minimize total completion time. We provide a polynomial time algorithm to solve the total completion time minimization problem.     

A Retinex-based total variation approach for image segmentation and bias correction

Image segmentation methods usually suffer from intensity inhomogeneity problem caused by many factors such as spatial variations in illumination (or bias fields of imaging devices). In order to address this problem, this paper proposes a Retinex-based variational model for image segmentation and bias correction. According to Retinex theory, the input inhomogeneous image can be decoupled into illumination bias and reflectance parts. The main contribution of this paper is to consider piecewise constant of the reflectance, and thereby introduce the total variation term in the proposed model for correcting and segmenting the input image. This is different from the existing model in which the spatial smoothness of the illumination bias is employed only. The existence of the minimizers to the variational model is established. Furthermore, we develop an efficient algorithm to solve the model numerically by using the alternating minimization method. Our experimental results are reported to demonstrate the effectiveness of the proposed method, and its performance is competitive with that of the other testing methods.     

Adaptive level set evolution starting with a constant function

In this paper, we propose a novel level set evolution model in a partial differential equation (PDE) formulation. According to the governing PDE, the evolution of level set function is controlled by two forces, an adaptive driving force and a total variation (TV)-based regularizing force that smoothes the level set function. Due to the adaptive driving force, the evolving level set function can adaptively move up or down in accordance with image information as the evolution proceeds forward in time. As a result, the level set function can be simply initialized to a constant function rather than the widely-used signed distance function or piecewise constant function in existing level set evolution models. Our model completely eliminates the needs of initial contours as well as re-initialization, and so avoids the problems resulted from contours initialization and re-initialization. In addition, the evolution PDE can be solved numerically via a simple explicit finite difference scheme with a significantly larger time step. The proposed model is fast enough for near real-time segmentation applications while still retaining enough accuracy; in general, only a few iterations are needed to obtain segmentation results accurately.     

Weak solutions of the forward problem in EEG for different conductivity values

The potential distribution on the scalp produced by current sources in the brain can be measured by an EEG recorder. The relationship between these sources and the scalp potential distribution may be described by a well-known mathematical model where some simplifications are usually introduced. The head is modeled as a multicompartment nested set and the conductivity of the different tissues is approximated by a positive piecewise constant function. This simplified model is used to solve the forward problem (FP), i.e., to calculate the scalp potential for a current source configuration. In this work, we prove that the weak solutions of the FP are continuous with respect to the conductivity values, that is, the difference between the scalp potentials is small if the conductivity values are closed enough. We present numerical examples that illustrates this property.     

Sugar cane transportation in Cuba,a case study

In this study the authors present a mixed integer linear programming model to solve the problem of cost minimization of sugar cane removal and its transport from the fields to the sugar mill at operational level. The complexity of the problem is basically determined by the system approach which results in the generation of a great number of variables and constraints that refer to the following operational dimensions: (a) need for continuous supply to the sugar mill; (b) cutting means used in cane harvesting; (c) transportation vehicles and (d) providing routes, which are characterized by the existence of storage facilities at the beginning of the railroads. The results demonstrate the model is not only useful to minimize transportation cost, but also for scheduling daily cane road transport and harvesting quotas of cutting means.     

基于改进C-V水平集模型的SAR图像分割

本文研究了SAR图像分割的问题.利用一种加入图像边缘信息且无需重新初始化的改进水平集方法,获得了比传统C-V模型分割速度更快、准确度更高的分割结果.推广了C-V水平集模型分割灰度不均匀的SAR图像以及零水平集曲线的初始化等结果.     

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     

Bogolyubov-type theorem with constraints induced by a control system with hysteresis effect

We consider the problem of minimization of an integral functional with nonconvex with respect to the control integrand. We minimize our functional over the solution set of a control system described by two ordinary differential equations subject to a control constraint given by a multivalued mapping with closed nonconvex values. The coefficients of the equations and the constraint depend on the phase variables. One of the equations contains the subdifferential of the indicator function of a closed convex set depending on the unknown phase variable. The equation containing the subdifferential describes an input–output relation of hysteresis type.     

基于小波多分辨率分析和改进窄带法的C-V水平集图像分割模型

本文研究了基于小波分析改进的C-V模型图像分割问题.利用小波多分辨率分析和改进的窄带水平集方法,获得了比传统C-V模型分割速度更快、准确度更高、算法复杂度更低的分割结果.推广了C-V水平集模型如何快速准确地分割灰度不均匀的图像和窄带水平集法等结果.     

求解一类逆鲁棒优化问题的牛顿型扰动方法

逆优化问题是指通过调整目标函数和约束中的某些参数使得已知的一个解成为参数调整后的优化问题的最优解.本文考虑求解一类逆鲁棒优化问题.首先,我们将该问题转化为带有一个线性等式约束,一个二阶锥互补约束和一个线性互补约束的极小化问题;其次,通过一类扰动方法来对转化后的极小化问题进行求解,然后利用带Armijo线搜索的非精确牛顿法求解每一个扰动问题.最后,通过数值实验验证该方法行之有效.