A Fast Region-Based Segmentation Model with Gaussian Kernel of Fractional Order

By summarizing some classical active contourmodels from the view of level set representation, a simple energy function expression with the Gaussian kernel of fractional order is proposed, and then a novel region-based geometric active contour model is established. In this proposed model, the energy function with value of [−1, 1] is built, the local mean and global mean of the inside and outside of the evolution curve are employed, and the segmentation results are obtained by controlling the expansion and contraction of the evolution curve. The model is simple and easy to implement; it can also protect weak edges because of considering more statistical information. Experimental results on synthetic and natural images show that the proposed model is much more effective in dealing with the images with weak or blurred edges, and it takes less time.


Introduction
Image segmentation is a basic and important topic in the fields of image processing.Accurate image segmentation can provide more important information for the follow-up application, such as machine vision and motion tracking.However, segmental results are always affected by low contrast and the problems of intensity inhomogeneity.The main idea of image segmentation is to extract the concerned regions and their contours from the whole image.There have been thousands of image segmentation algorithms proposed in recent decades.Some researchers put forward the edge detection based on the gradient, derivatives, or Canny edge detection, and so on.Edge detection is good for simple image but not suitable for the clutter target boundary extraction.The main reasons are as follows.Firstly, edge extracted for complex image is often not corresponding to the target boundary.Secondly, the extracted edge is discontinuous, but the goal often needs closed boundary to separate the object from the whole image.In addition, edge detection is dependent on the local information near pixel; it has advantages sometimes, but in many cases overall appearance of the target is the key, so the concepts of the image segmentation and edge detection are not one and the same.
Regional growth is a simple technique to provide segmental region; the algorithm begins with some seed points and found pixels near the seed which has similar image characteristics, such as gray scale and color characteristics.This algorithm has been applied to Mumford-Shah function [1].Another region-based method is active contour (AC) model [2].Active contour model is 2D or 3D surface contour description, which involves the contour evolution under an appropriate energy in order to get a satisfactory segmentation result, such as the target boundary with the closed contour.Over the past decade, researchers have proposed many different active contour models, which are mainly divided into two categories, namely, parametric active contour models and geometric active contour models.In parametric active contour models, the parameter equation of the curve is () = [(), ()], in which 0 ≤  ≤ 1.The parametric active contour model essentially depends on the energy function rather than the geometric figures of the contour.Therefore, this model cannot handle topology changes when it detects multiple targets, but geometric active contour model can deal with topological changes, because it uses the structure of level set, in which the curve  is zero level set function (, ) :   × [0, ∞) → ,  = {2, 3}; for example,  = { ∈   : (, ) = 0}.The first type of geometric active contour model is introduced by Caselles et al. [3]; its main idea is to use curvature and normal direction forcing curve movement, so that it stops on the edge with the edge function () = (|∇| 2 ), where ∇ means the gradient of the given image , which has a property; that, it equals zero on the border, and others equal one.For example, () =  −(1/ 2  ) |∇  * ()| 2 , where   is a scale factor,   = (1/ √ 2) −((−) 2 / 2 ) is a Gaussian kernel, in which  denotes the standard deviation of the given image, and  denotes the expectation of the given image.Another type of geometric active contour model is the Geodesic active contour model [4], which can search for the minimum length of the edge weights under the energy function.This model is similar to the former geometric model, but there is a big difference, a vector filed term is employed in Geodesic active contour model to stop the motion curve on the weak edges.Paragios et al. put forward the famous gradient vector flow (GVF) instead of ∇ to increase the range of results, called GVF geometric AC [5].Chan and Vese [6] proposed a new model CV AC, and Li et al. [7] proposed a new model (LBF) which uses energy function to overcome the problem of nonhomogeneity.LBF model can deal with the image of different gray levels by adding the kernel function, and it can employ local gray level information effectively.
Many structures of different level set evolution models have been summarized before.The level set evolution of the above energy functions [8] where  is the Euclidean curvature and   ,   , and ⃗  are three parameters, which decide the speed and direction of the evolution.The term based on curvature vector is used to smooth the curve.The normal direction is used to control shrinkage and expansion of the curve and force the curve to move along the direction vector.Details are shown in Table 1, where , ,  1 , and  2 are constant, [û, V] is the GVF, and () is a function based on the normal curvature and GVF [9].At the same time, () in Table 1 is original gray level image,  in and  out are average gray values of () inside and outside of curve, and  1 and  2 are the weighted average gray values of () inside and outside of curve in the Gaussian window.
The rest of the paper is organized as follows: in the next two sections, we will review classical existing geometric models, Chan-Vese model and LBF model.The new model is introduced in Section 4. Some experimental results are shown in Section 5. We conclude the paper in Section 6.

Chan-Vese Model
In Chan-Vese (CV) model, we considered the simplest type of segmentation, which divided the image into the target and the background, and the distributions of the gray values of target and background are approximately constant values.CV model is based on the evolution of the level set and can deal with curve topology changes better for the curve which is expressed by the level set function.The energy function of CV model is where  0 is the given image,  is the evolution curve, () : [0, 1] →  is a parameter evolution curve, and  0 is the weight coefficient.In the energy function, the first item is the length of the curve evolution and it can regularize the curve.
The last two items are global binary fitting items.The basic idea of the model is to minimize the fitting item  1 + 2 , where The level set function  is defined as The following are Dirac function and Heaviside function, respectively: Because function () cannot directly take the derivative of , we can replace  with   () in the CV model, where Advances in Mathematical Physics 3 The level set function of the CV model is: where From the definition of the level set function , expressions of  1 ,  2 are, respectively, as follows: According to Euler-Lagrange equation, the level set expression of CV model is obtained as

LBF Model
LBF model defines a local binary fitting energy item, which is actually a kernel function; the model is as follows: where and is a Gaussian kernel function,   and   are its expectancies, and  is its standard deviation. 1 and  2 are the image fitting function of the local gray level inside and outside of the contour.The variational level set function of (11) which is got by Euler-Lagrange equation is as follows: By introducing a sign distance constraint and length constraint item, the level set evolution equation is where functions  1 ,  2 and  1 ,  2 are as follows, respectively: Equation (15) shows that  1 and  2 are the weighted average gray values with the Gaussian window inside and outside of contour.Obviously, they share the local characteristics, so that the segmentation of original image by LBF model is more accurate.

Proposed Region-Based Model with Gaussian Kernel of Fractional Order
In order to get better image segmentation results effectively and construct a fast region-based segmentation model, we should keep the energy functional as simple as possible, and energy information must be used effectively.Based on the law of some classical energy function expressions of active contour models summarized in Table 1, we know that the term   in energy function is very important, and many models made a breakthrough on it.At the same time, the item ⃗  of the models is often set as ⃗ 0, and item   is only ordinary parameters.Therefore, only keeping item   in the new model will simplify the expression of energy function.
In order to avoid jumping internally, level set function initialized by symbolic distance function (SDF) in the traditional level set method, but it usually needs to be reinitialized.This will lead to the fact that it is hard to decide when to reinitialize and how to reinitialize, as it is hard to find boundary when the zero level set is away from the inner region.So reinitialization is a very complex operation problem.To solve this problem, we propose a new level set method.At the same time, fractional systems [10,11] gain increasing attention in applied sciences, and functions of fractional order are more flexible, so the new method uses a Gaussian filter with fractional order to regularize binary level set function.The traditional level set method uses curvature item div(∇/|∇|)|∇| to regularize the level set function, and letting |∇| = 1 [12], it can replace the regular items with Laplacian.Based on scale space theory in [13], a function with the Laplacian evolution is equivalent to using a Gaussian filter.Then, we filter the initial conditions of the level set function with Gaussian kernel filters of the level set function, and  controls the regular strength, similar to the item  in Table 1.With Gaussian kernel function, the item div(∇/|∇|)|∇| (similar to the item   in Table 1) can be removed, so the key of whole model is the choice of item   .
From the view of level set function, we need to find a function that can adjust the pressure inside and outside of the interest areas.It drives the curve to contract when the curve is outside the target and expands when the curve is within the target.Based on SPF function with value of [−1, 1] defined in [14], we can construct a function as follow: where where   () is the Gaussian kernel of fractional order .It takes the ordinary Gaussian kernel as its special case when  = 2. Thus, it is more flexible than the ordinary one.We must emphasize that the parameter  may be different from fractal parameters [15,16].In this paper, we call , likely informal, the fractional order of the Gaussian kernel expressed by (19).In the following experiments, we try different values of the fractal order  in the evolution level set function.The value of the function is between −1 and 1.It drives curve to contract externally the target and expands when the curve is within the target.According to the summary of the classic model of general expression (1), we only keep item   , so we obtain the corresponding variational level set formulation as follows: By adding a parameter, the final level set equation of the new model is The main algorithm of the new model is as follows.
Step 1: initialize the level set function .
Step 2: compute the mean and the weighted average values of inside and outside of the curve  1 ,  2 ,  1 , and  2 .
Step 3: compute the evolution of level set function by (21).
Step 4: use Gaussian filter to regularize level set function,  =  *   .

Experimental Results
In this section, we will show some experimental results of the proposed model on synthetic image and nature image; the results also will be compared with those got by the conventional CV model and LBF model.

Conclusion
Inspired by the idea of some classical energy function expressions of active contour model, from the view of level set representation, a novel fast region-based segmentation model with Gaussian kernel of fractional order is proposed.The model is simple and easy to be implementated, and it can protect weak edges because of considering more statistical information.The experimental results on synthetic images and natural images show that the proposed model is superior to the traditional methods.The new model is much more effective in dealing with the images with weak or blurred edges, and it takes less time.

Figure 1 :
Figure 1: The segmentation results of the aircraft composite image: (a) and (d) results using CV model, (b) and (e) results using LBF model, and (c) and (f) results using the new model.

Figure 2 :Figure 3 :
Figure 2: The segmentation results of the natural image: (a) and (d) results using CV model, (b) and (e) results using LBF model, and (c) and (f) results using the new model.

Figure 2
shows the experimental results of natural figure with objects having inhomogeneous background.CV model, LBF model, and the proposed model share the same environment of the initial value.We can find that the gray level of the natural star figure is extremely inhomogeneous easily.The first line shows the segmentation result of CV model.The second line shows the segmentation result of LBF model.The third line shows the segmentation result of new model.The segmentation images reveal that the proposed model gets the most ideal segmentation result.

Figure 3 Figure 4 :
Figure 4: The segmentation results of the blood vessels image: (a) and (d) results using CV model, (b) and (e) results using LBF model, and (c) and (f) results using the new model.

Figure 5 :
Figure 5: The segmentation results of the CT image: (a) and (d) results using CV model, (b) and (e) results using LBF model, and (c) and (f) results using the new model.

Table 1 :
The analysis table of energy function model.
Our algorithm is implemented in Windows 7 Operating System, i3 Dual Core CPU 2.13 GHz and 2 GB RAM.The initial value and parameters such as time step take different values in the specific experiments in this paper.Figures 1(d)-1(f) got with 2000, 1500, 40 iterations separately and consume time  = 217.6 s, 232.6 s, 13.8 s in turn.We can see that Chan-Vese model and LBF model cannot get