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

Authors: Bo Chen; Qing-Hua Zou; Wen-Sheng Chen; Yan Li
Journal: Advances in Mathematical Physics
 (2013)

---

## Abstract

By summarizing some classical active contour models 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.


---

## Body

## 1. 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 C(p)=[x(p),y(p)], in which 0≤p≤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 C is zero level set function ϕ(x,t):Rn×[0,∞)→R,n={2,3}; for example, C={x∈Rn:ϕ(x,t)=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 g(x)=g(|∇I|2), where ∇I means the gradient of the given image I, which has a property; that, it equals zero on the border, and others equal one. For example, g(x)=e-(1/σe2)|∇Gσ*I(x)|2, where σe is a scale factor, Gσ=(1/2πσ)e-((x-μ)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 ∇g 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] can be represented as
(1)∂ϕ∂t=αkK|∇ϕ|+VN|∇ϕ|+N→·∇ϕ,
where K is the Euclidean curvature and αk, VN, and N→ 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, [u^,v^] is the GVF, and k(x) is a function based on the normal curvature and GVF [9]. At the same time, I(x) in Table 1 is original gray level image, Cin and Cout are average gray values of I(x) inside and outside of curve, and  e1 and e2 are the weighted average gray values of I(x) inside and outside of curve in the Gaussian window.Table 1
The analysis table of energy function model.
Model
α
k
V
N
N
→
Geometric AC [3]
g
(
x
)
r
g
(
x
)
0
→
Geodesic AC [4]
g
(
x
)
r
g
(
x
)
∇
g
(
x
)
GVF Geo. AC [5]
g
(
x
)
r
k
(
x
)
g
(
x
)
g
(
x
)
(
1
-
|
k
(
x
)
|
[
u
^
,
ν
^
]
)
CV AC [6]
η
μ
+
(
I
-
c
i
n
)
2 -(I-cout)2
0
→
LBF AC [7]
η
λ
1
e
1
-
λ
2
e
2
0
→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 Section4. Some experimental results are shown in Section 5. We conclude the paper in Section 6.


## 2. 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(2)E(C,c1,c2)=μ0Length(C)+λ1∫inside(C)|u0(x,y)-c1|2dxdy+λ2∫outside(C)|u0(x,y)-c2|2dxdy,
where u0 is the given image, C is the evolution curve, C(s):[0,1]→R 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 F1+F2, where
(3)F1=∫inside(C)|u0-c1|2dxdy,F2=∫outside(C)|u0-c2|2dxdy.The level set functionϕ is defined as
(4)ϕ(x,y)>0,(x,y)∈in(C),ϕ(x,y)<0,(x,y)∈out(C),ϕ(x,y)=0,(x,y)∈(C).The following are Dirac function and Heaviside function, respectively:(5)δ(ϕ)=dH(ϕ)dϕ,H(ϕ)={1,ϕ≥0,0,ϕ<0.Because functionH(ϕ) cannot directly take the derivative of ϕ, we can replace H with Hε(ϕ) in the CV model, where
(6)Hε(ϕ)=12(1+2πarctan(ϕε)).The level set function of the CV model is:(7)E(ϕ,c1,c2)=μ∫Ωδε(ϕ)|∇ϕ|dxdy+λ1∫Ω|u0(x,y)-c1|2Hε(ϕ)dxdy+λ2∫Ω|u0(x,y)-c2|2(1-Hε(ϕ))dxdy,
where
(8)δε(ϕ)=Hε′(ϕ)=1πεε2+π2.
From the definition of the level set function ϕ, expressions of c1, c2 are, respectively, as follows:
(9)c1=∫Ωu0(x,y)H(ϕ(x,y))dxdy∫ΩH(ϕ(x,y))dxdy,c2=∫Ωu0(x,y)(1-H(ϕ(x,y)))dxdy∫Ω(1-H(ϕ(x,y)))dxdy.
According to Euler-Lagrange equation, the level set expression of CV model is obtained as
(10)∂ϕ∂t=δ(ϕ)[μdiv(∇ϕ|∇ϕ|)-λ1(u0-c1)2+λ2(u0-c2)2].


## 3. LBF Model

LBF model defines a local binary fitting energy item, which is actually a kernel function; the model is as follows:(11)E(C,f1,f2)=λ1∫Ω∫inside(C)Kσ(x-y)|I(y)-f1(x)|2dydx+λ2∫Ω∫outside(C)Kσ(x-y)|I(y)-f2(x)|2dydx,
where I:Ω→R is the original image, (x,y)∈Ω, Kσ=(1/2πσ2)e-((x-μx)2+(y-μy)2/σ2) and is a Gaussian kernel function, μx and μy are its expectancies, and σ is its standard deviation. f1 and f2 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:
(12)∂ϕ∂t=-δ(ϕ)(λ1e1-λ2e2).By introducing a sign distance constraint and length constraint item, the level set evolution equation is(13)∂ϕ∂t=μ(Δϕ-div(∇ϕ|∇ϕ|))+vδ(ϕ)div(∇ϕ|∇ϕ|)-δ(ϕ)(λ1e1-λ2e2),
where functions e1, e2 and f1, f2 are as follows, respectively:
(14)e1(x)=∫ΩKσ(y-x)|I(x)-f1(y)|2dy,e2(x)=∫ΩKσ(y-x)|I(x)-f2(y)|2dy,(15)f1=Kσ*[H(ϕ)I(x)]Kσ*H(ϕ),f2=Kσ*[(1-H(ϕ))I(x)]Kσ*((1-H(ϕ))).Equation (15) shows that f1 and f2 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.


## 4. 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 Table1, we know that the term VN in energy function is very important, and many models made a breakthrough on it. At the same time, the item N→ of the models is often set as 0→, and item αk is only ordinary parameters. Therefore, only keeping item VN 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 αk in Table 1) can be removed, so the key of whole model is the choice of item VN.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:
(16)f(I(x))=(c1+c2+f1+f2)/4-I(x)max(|(c1+c2+f1+f2)/4-I(x)|),x∈Ω,
where
(17)Min(I(x))≤c1,c2,f1,f2≤Max(I(x)),Min(I(x))≤(c1+c2+f1+f2)4≤Max(I(x)),(18)f1=Kσ*[H(ϕ)I(x)]Kσ*H(ϕ),f2=Kσ*[(1-H(ϕ))I(x)]Kσ*(1-H(ϕ)),(19)Kσ(α)=12πσ2e-(∥x-μ∥α/σ2),
where Kσ(α) 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 VN, so we obtain the corresponding variational level set formulation as follows:
(20)∂ϕ∂t=f(I(x))|∇ϕ|,x∈Ω.
By adding a parameter, the final level set equation of the new model is
(21)∂ϕ∂t=ηf(I(x))|∇ϕ|,x∈Ω.
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 curvec1, c2, f1, and f2.Step 3:
compute the evolution of level set function by (21).Step 4:
use Gaussian filter to regularize level set function,ϕ=ϕ*Gσ.Step 5:
repeat Steps 2 to 4 until convergence.


## 5. 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. 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.Figures1(d)–1(f) got with 2000, 1500, 40 iterations separately and consume time t=217.6 s, 232.6 s, 13.8 s in turn. We can see that Chan-Vese model and LBF model cannot get satisfactory result. While here the proposed model gets satisfactory result, meanwhile the proposed model costs less time.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.
(a)
(b)
(c)
(d)
(e)
(f)Figure2 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.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.
(a)
(b)
(c)
(d)
(e)
(f)Figure3 shows the experimental results of natural figure with CV model, LBF model, and the proposed model; gray level of the image is extremely inhomogeneous. The three models share the same initial value. Figures 3(a) and 3(d) show the segmentation result of CV model; Figures 3(b) and 3(e) show the segmentation result of LBF model; Figures 3(c) and 3(f) show the segmentation result of the new model. The segmentation results show that the proposed model gets the satisfactory complete result, while the CV model and LBF model have some redundant and inaccurate segmentation.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.
(a)
(b)
(c)
(d)
(e)
(f)In Figure4, the experiment shows the segmentation results of a gray matter blood vessels image. The first column shows the segmentation results of Chan-Vese model. The second column shows the results of LBF model. The third column shows the result of proposed model. Figures 4(d)–4(f) got with 1000, 5000, 40 iterations separately and consume time t=79.2 s, 174.2 s, 6.1 s in turn. Because the mean information of Chan-Vese model is very sensitive to inhomogeneity image, it fails to extract the accurate contour. The LBF model is better than Chan-Vese model, but it gets  many iterations and it is very sensitive to initial value. As shown in the low left of the segmental contour, the proposed model got a better segmentation result.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.
(a)
(b)
(c)
(d)
(e)
(f)In Figure5, experimental results show the segmentation results of CT bone images. The initialization is a single circle. The first column shows the segmentation result of Chan-Vese model. The second column shows the result of the LBF model. The third row shows the result of proposed model. Figures 5(d)–5(f) are the corresponding contours to Figures 5(a)–5(c). Figures 5(d)–5(f) got with 50, 1000, 80 iterations separately and consume time t=4.0 s, 34.6 s, 12.5 s in turn. As seen from Figure 5, the new method can get more complete contour than Chan-Vese model and the LBF model and take less time.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.
(a)
(b)
(c)
(d)
(e)
(f)


## 6. 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.



---

*Source: 501628-2013-10-20.xml*