A fractional order derivative based active contour model for inhomogeneous image segmentation

Segmenting intensity inhomogeneous images is a challenging task for both local and global methods. Some hybrid methods have great advantages over the traditional methods in inhomogeneous image segmentation. In this paper, a new hybrid method is presented, which incorporates image gradient, local environment and global information into a framework, called adaptive-weighting active contour model. The energy or level set functions in the framework mainly include two parts: a global term and local term. The global term aims to enhance the image contrast, and it can also accelerate the convergence rate when minimizing the energy function. The local term integrates fractional order differentiation, fractional order gradient magnitude, and difference image information into the well-known local Chan–Vese model, which has been shown to be effective and efficient in modeling the local information. The local term can also enhance low frequency information and improve the inhomogeneous image segmentation. An adaptive weighting strategy is proposed to balance the actions of the global and local terms automatically. When minimizing the level set functions, regularization can be imposed by applying Gaussian filtering to ensure smoothness in the evolution process. In addition, a corresponding stopping criterion is proposed to ensure the evolving curve automatically stops on true boundaries of objects. Dice similarity coefficient is employed as the comparative quantitative measures for the segmented results. Experiments on synthetic images as well as real images are performed to demonstrate the segmentation accuracy and computational efficiency of the presented hybrid method.     

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.     

A Locally-Biased form of the DIRECT Algorithm

In this paper we propose a form of the DIRECT algorithm that is strongly biased toward local search. This form should do well for small problems with a single global minimizer and only a few local minimizers. We motivate our formulation with some results on how the original formulation of the DIRECT algorithm clusters its search near a global minimizer. We report on the performance of our algorithm on a suite of test problems and observe that the algorithm performs particularly well when termination is based on a budget of function evaluations.     

A kriging based method for the solution of mixed-integer nonlinear programs containing black-box functions

In this paper a new methodology is developed for the solution of mixed-integer nonlinear programs under uncertainty whose problem formulation is complicated by both noisy variables and black-box functions representing a lack of model equations. A branch-and-bound framework is employed to handle the integer complexity whereby the solution to the relaxed nonlinear program subproblem at each node is obtained using both global and local information. Global information is obtained using kriging models which are used to identify promising neighborhoods for local search. Response surface methodology (RSM) is then employed whereby local models are sequentially optimized to refine the problem’s lower and upper bounds. This work extends the capabilities of a previously developed kriging-response surface method enabling a wider class of problems to be addressed containing integer decisions and black box models. The proposed algorithm is applied to several small process synthesis examples and its effectiveness is evaluated in terms of the number of function calls required, number of times the global optimum is attained, and computational time.     

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.     

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.     

Square-root dynamics of a giving up smoking model

The paper presents a new giving up smoking model for which interaction term is square root of potential and occasional smokers of model presented in Zaman (2011) [15]. First, we will show formulation of the model. Then we will discuss local and global stability of the model and its general solutions. The non-standard finite difference method (NSFD) is used to solve the new giving up smoking model. Both non-negativity and conservative law for differential equations system are discussed. Numerical results are presented graphically and compared well with those obtained by Runge–Kutta fourth-order method (RK4) and ODE45.     

A computational methodology to study coupled physical processes over partitioned domains

In this paper, we describe a computational methodology to couple physical processes defined over independent subdomains, that are partitions of a global domain in three-dimensions. The methodology presented helps to compute the numerical solution on the global domain by appropriately piecing the local solutions from each subdomain. We discuss the mixed method formulation for the technique applied to a model problem and derive an error estimate for the finite element solution. We demonstrate through numerical experiments that the method is robust and reliable in higher dimensions.     

Nonlinear multigrid method for solving the anisotropic image denoising models

In this paper, we study a nonlinear multigrid method for solving a general image denoising model with two L ¹-regularization terms. Different from the previous studies, we give a simpler derivation of the dual formulation of the general model by augmented Lagrangian method. In order to improve the convergence rate of the proposed multigrid method, an improved dual iteration is proposed as its smoother. Furthermore, we apply the proposed method to the anisotropic ROF model and the anisotropic LLT model. We also give the local Fourier analysis (LFAs) of the Chambolle’s dual iterations and a modified smoother for solving these two models, respectively. Numerical results illustrate the efficiency of the proposed method and indicate that such a multigrid method is more suitable to deal with large-sized images.     

Inverse Stackelberg Public Goods Game with Multiple Hierarchies Under Global and Local Information Structures

This paper studies the inverse Stackelberg game with multiple hierarchies under global and local information structures, where the players have discrete strategy spaces. For the classic public goods game, we solve the pure-strategy inverse Stackelberg equilibria under three typical hierarchical structures. The results reveal some counterintuitive characteristics within the systems with hierarchies, such as that the cooperation does not increase with the return rate at the equilibria. Furthermore, by defining a local information structure, we give an estimate of the fewest hierarchies required for full cooperation, which can be a constant multiple of the logarithm or square root of the population size or of the population size itself, according to different information structures and return rates. This paper proposes a novel mechanism to play the game and promote cooperation. Both the formulation and analysis method are different from existing works, and the results can find their ample implications in practice, which might help decision making in hierarchical systems.     

A Shared-Constraint Approach to Multi-Leader Multi-Follower Games

Multi-leader multi-follower games are a class of hierarchical games in which a collection of leaders compete in a Nash game constrained by the equilibrium conditions of another Nash game amongst the followers. The resulting equilibrium problem with equilibrium constraints is complicated by nonconvex agent problems and therefore providing tractable conditions for existence of global or even local equilibria has proved challenging. Consequently, much of the extant research on this topic is either model specific or relies on weaker notions of equilibria. We consider a modified formulation in which every leader is cognizant of the equilibrium constraints of all leaders. Equilibria of this modified game contain the equilibria, if any, of the original game. The new formulation has a constraint structure called shared constraints, and our main result shows that if the leader objectives admit a potential function, the global minimizers of the potential function over this shared constraint are equilibria of the modified formulation. We provide another existence result using fixed point theory that does not require potentiality. Additionally, local minima, B-stationary, and strong-stationary points of this minimization problem are shown to be local Nash equilibria, Nash B-stationary, and Nash strong-stationary points of the corresponding multi-leader multi-follower game. We demonstrate the relationship between variational equilibria associated with this modified shared-constraint game and equilibria of the original game from the standpoint of the multiplier sets and show how equilibria of the original formulation may be recovered. We note through several examples that such potential multi-leader multi-follower games capture a breadth of application problems of interest and demonstrate our findings on a multi-leader multi-follower Cournot game.     

Accelerating Branch-and-Bound through a Modeling Language Construct for Relaxation-Specific Constraints

In the tradition of modeling languages for optimization, a single model is passed to a solver for solution. In this paper, we extend BARON’s modeling language in order to facilitate the communication of problem-specific relaxation information from the modeler to the branch-and-bound solver. This effectively results into two models being passed from the modeling language to the solver. Three important application areas are identified and computational experiments are presented. In all cases, nonlinear constraints are provided only to the relaxation constructor in order to strengthen the lower bounding step of the algorithm without complicating the local search process. In the first application area, nonlinear constraints from the reformulation–linearization technique (RLT) are added to strengthen a problem formulation. This approach is illustrated for the pooling problem and computational results show that it results in a scheme that makes global optimization nearly as fast as local optimization for pooling problems from the literature. In the second application area, we communicate with the relaxation constructor the first-order optimality conditions for unconstrained global optimization problems. Computational experiments with polynomial programs demonstrate that this approach leads to a significant reduction of the size of the branch-and-bound search tree. In the third application, problem-specific nonlinear optimality conditions for the satisfiability problem are used to strengthen the lower bounding step and are found to significantly expedite the branch-and-bound algorithm when applied to a nonlinear formulation of this problem.     

On cautious probabilistic inference and default detachment

Conditional probabilities are one promising and widely used approach to model uncertainty in information systems. This paper discusses the DUCK-calculus, which is founded on the cautious approach to uncertain probabilistic inference. Based on a set of sound inference rules, derived probabilistic information is gained by local bounds propagation techniques. Precision being always a central point of criticism to such systems, we demonstrate that DUCK need not necessarily suffer from these problems. We can show that the popular Bayesian networks are subsumed by DUCK, implying that precise probabilities can be deduced by local propagation techniques, even in the multiply connected case. A comparative study with INFERNO and with inference techniques based on global operations-research techniques yields quite favorable results for our approach. Since conditional probabilities are also suited to model nonmonotonic situations by considering different contexts, we investigate the problems of maximal and relevant contexts, needed to draw default conclusions about individuals.     

A Petrov–Galerkin RBF method for diffusion equation on the unit sphere

This paper concerns a numerical solution for the diffusion equation on the unit sphere. The given method is based on the spherical basis function approximation and the Petrov–Galerkin test discretization. The method is meshless because spherical triangulation is not required neither for approximation nor for numerical integration. This feature is achieved through the spherical basis function approximation and the use of local weak forms instead of a global variational formulation. The local Petrov–Galerkin formulation allows to compute the integrals on small independent spherical caps without any dependence on a connected background mesh. Experimental results show the accuracy and the efficiency of the new method.     

Analysis of finite difference schemes for unsteady Navier-Stokes equations in vorticity formulation

Summary. In this paper, we provide stability and convergence analysis for a class of finite difference schemes for unsteady incompressible Navier-Stokes equations in vorticity-stream function formulation. The no-slip boundary condition for the velocity is converted into local vorticity boundary conditions. Thom's formula, Wilkes' formula, or other local formulas in the earlier literature can be used in the second order method; while high order formulas, such as Briley's formula, can be used in the fourth order compact difference scheme proposed by E and Liu. The stability analysis of these long-stencil formulas cannot be directly derived from straightforward manipulations since more than one interior point is involved in the formula. The main idea of the stability analysis is to control local terms by global quantities via discrete elliptic regularity for stream function. We choose to analyze the second order scheme with Wilkes' formula in detail. In this case, we can avoid the complicated technique necessitated by the Strang-type high order expansions. As a consequence, our analysis results in almost optimal regularity assumption for the exact solution. The above methodology is very general. We also give a detailed analysis for the fourth order scheme using a 1-D Stokes model. Received December 10, 1999 / Revised version received November 5, 2000 / Published online August 17, 2001     

Adaptive Bayesian Nonstationary Modeling for Large Spatial Datasets Using Covariance Approximations

Gaussian process models have been widely used in spatial statistics but face tremendous modeling and computational challenges for very large nonstationary spatial datasets. To address these challenges, we develop a Bayesian modeling approach using a nonstationary covariance function constructed based on adaptively selected partitions. The partitioned nonstationary class allows one to knit together local covariance parameters into a valid global nonstationary covariance for prediction, where the local covariance parameters are allowed to be estimated within each partition to reduce computational cost. To further facilitate the computations in local covariance estimation and global prediction, we use the full-scale covariance approximation (FSA) approach for the Bayesian inference of our model. One of our contributions is to model the partitions stochastically by embedding a modified treed partitioning process into the hierarchical models that leads to automated partitioning and substantial computational benefits. We illustrate the utility of our method with simulation studies and the global Total Ozone Matrix Spectrometer (TOMS) data. Supplementary materials for this article are available online.     

Market Viability and Martingale Measures under Partial Information

We consider a financial market model with a single risky asset whose price process evolves according to a general jump-diffusion with locally bounded coefficients and where market participants have only access to a partial information flow. For any utility function, we prove that the partial information financial market is locally viable, in the sense that the optimal portfolio problem has a solution up to a stopping time, if and only if the (normalised) marginal utility of the terminal wealth generates a partial information equivalent martingale measure (PIEMM). This equivalence result is proved in a constructive way by relying on maximum principles for stochastic control problems under partial information. We then characterize a global notion of market viability in terms of partial information local martingale deflators (PILMDs). We illustrate our results by means of a simple example.     

A nonmonotone filter method for nonlinear optimization

We propose a new nonmonotone filter method to promote global and fast local convergence for sequential quadratic programming algorithms. Our method uses two filters: a standard, global g-filter for global convergence, and a local nonmonotone l-filter that allows us to establish fast local convergence. We show how to switch between the two filters efficiently, and we prove global and superlinear local convergence. A special feature of the proposed method is that it does not require second-order correction steps. We present preliminary numerical results comparing our implementation with a classical filter SQP method.     

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.