A Novel Adaptive Level Set Segmentation Method

The adaptive distance preserving level set (ADPLS) method is fast and not dependent on the initial contour for the segmentation of images with intensity inhomogeneity, but it often leads to segmentation with compromised accuracy. And the local binary fitting model (LBF) method can achieve segmentation with higher accuracy but with low speed and sensitivity to initial contour placements. In this paper, a novel and adaptive fusing level set method has been presented to combine the desirable properties of these two methods, respectively. In the proposed method, the weights of the ADPLS and LBF are automatically adjusted according to the spatial information of the image. Experimental results show that the comprehensive performance indicators, such as accuracy, speed, and stability, can be significantly improved by using this improved method.


Introduction
Since the introduction by Kass et al. [1], active contour models (ACMs) have been widely used in image segmentation [2][3][4]. The existing ACMs based on the level set method initially proposed to handle the topological changes during the curve evolution can be broadly classified as either edge-based models [5][6][7] or region-based models [8][9][10][11][12][13][14][15] according to the type of adopted image features. The basic idea of the level set method is to implicitly embed the moving contour into a higher dimensional level set function and view the contour as its zero level set [16]. ACMs have desirable properties over the conventional image segmentation methods, such as thresholding, edge detection, and region growing. First, ACMs can provide closed and smooth contours in segmentation results, which are necessary for further application such as shape analysis and recognition. Second, ACMs can get the object boundaries with subpixel accuracy.
Edge-based models utilize image gradient to stop evolving contours on the object boundaries. Recently, He et al. proposed adaptive distance preserving level set (ADPLS) [7] evolution for image segmentation, in which the initial curve is no longer required to surround (let in or keep out) the objects to be detected. Moreover, starting with only one initial curve whose position is anywhere in the image, it can automatically detect interior and exterior contours of an object and edges of multiobjects. In addition, a large time step can be used to speed up the curve evolution in the numerical solution of the partial differential equation. However, this method does not contain any local intensity information, which is crucial for segmentation of images with intensity inhomogeneity. As a consequence, the ADPLS method generally fails to segment images with significant intensity inhomogeneity, which is illustrated in the following vessel segmentation in Figure 1. Figure 1(a) is a typical image with intensity inhomogeneity, and Figure 1(b) shows the final segmentation result of the ADPLS method. It is obvious that boundary leakage arises in regions 1 and 3, and less segmentation arises in regions 2, 4, and 5. This example shows the inability of the ADPLS method to segment images with intensity inhomogeneity.
To guide the motion of the active contour, the regionbased models identify each region of interest by using a certain region descriptor such as intensity, color, texture, or motion. Region-based models have better performance than edge-based models in the presence of weak boundaries. Local binary fitting model (LBF) [14] is one of the classical region-based models. It can not only segment the image with weak boundaries well, but also overcome the segmentation error caused by intensity inhomogeneity. The LBF model draws upon spatially varying local region information; thus it is able to deal with intensity inhomogeneity. Figure 1(c) shows the final segmentation result of the LBF model, which shows much better results than that of the ADPLS method. However, compared to the ADPLS method, the LBF model also has its own disadvantages. First, the LBF model is sensitive to the initial contour, and inappropriate initial contours might lead to failure of segmentation. Second, its curve evolution is slow due to some limitation of the time step.
In this paper we proposed a novel adaptive level set method to combine the good properties of both ADPLS and LBF methods. In a simultaneous and automatic way, the proposed method adjusts the proportion of the ADPLS and LBF methods according to spatial image information adaptively. As a consequence, the advantages of the ADPLS method and the LBF method are exploited with minimized disadvantages. The following experiments on both simulated and real images show that the proposed method can achieve segmentation with higher accuracy.

The ADPLS Method.
Let Ω ⊂ 2 be the image domain, and let : Ω → be a given gray level image. The ADPLS method [7] was recently proposed to overcome the disadvantage of the distance preserving level set method proposed by He et al., which requires the initial curve surrounding (let in or keep out) the objects to be detected. It imports a variable weighting coefficient whose sign symbol and size are adjusted by image information. So the zero level set can choose its evolution direction adaptively.
In image segmentation, active contours are dynamic curves that move toward the object boundaries by minimizing a predefined energy functional. Let be the edge indicator function defined by where is the Gaussian kernel with standard deviation and > 0 is a constant.
In the ADPLS method, a variational framework on the level set function is defined as follows: where the first term is the internal energy of that characterizes the deviation of the level set function from a signed distance function, > 0 is the weight of this internal energy term, the second term computes the length of the zero level curve of , is the edge indicator function defined by (1), is the univariate Dirac function, is the Heaviside function, and > 0 is the weight of the second term. The definition of V( ) is given by where > 0 is a constant, sgn( ) is the sign function, and Δ × denotes the image convolved with a Gaussian Computational and Mathematical Methods in Medicine 3 smoothing filter and calculated by Laplace operator. The associated level set evolution equation is given by where Δ ADPLS is the time step; the definition of 1 is given by It is worth mentioning that the variable weighting coefficient V( ) in (4) plays a key role in the evolution of the zero level curves. First, it adaptively guides the zero level curves to the target contour according to the image information. Second, the size of it is adjusted adaptively according to the spatial image information, which can greatly improve the capability of zero level set to detect the edge of multicontour and concavities. Third, by adjusting the size of coefficient , the method can control the zero level curves' ability of capturing the target boundaries. If multicontours exist in the image, we give the coefficient a larger value. On the contrary, if the content of the image is simple, a smaller value is set to coefficient .
Although the ADPLS method has the above advantages, it also has disadvantage as follows. In the segmentation of the images with intensity inhomogeneity, the intensity in one place may vary dramatically from another place, no matter how the edge indicator function ( ) is adjusted; it also has the possibility that the decrease of ( ) is too quickly for one place but too slowly for another one. The situation can be seen from Figure 1(b).

The LBF Model.
Recently Li et al. [14] proposed the LBF model, which can overcome the segmentation error brought by intensity inhomogeneity using two fitting functions 1 ( ) and 2 ( ) which locally approximate the intensities outside and inside the contours. They extracted the object by minimizing the following energy function: where is a Gaussian kernel with standard deviation . The first two terms in (6) are the weighted mean square error of the approximation of the image intensities ( ) outside and inside the contour by the fitting values 1 ( ) and 2 ( ), respectively, with ( − ) as the weight assigned to each intensity ( ) at . The third term | | is the length of the contour, and the last term is the regularization term for level set evolution. 1 > 0, 2 > 0, > 0, and V > 0 are constants. The associated level set evolution equation is given by where Δ LBF is the time step, and the definitions of the edge indicator function 1 , 2 , 1 , and 2 are given by Because of the localization property of the kernel function ( − ), the contribution of the intensity ( ) to the fitting energy LBF decreases to zero as the point goes away from the center point . This localization property plays a key role in segmenting the image with intensity inhomogeneity. And a better result than that of the ADPLS method can be observed from Figure 1(c). However, many local minimums of the energy functional might be simultaneously introduced by this localization property, which means that the LBF model might be sensitive to initial curve, and inappropriate initial curve might lead to the failure of segmentations. Figure 2 shows the failure situations of the segmentations with inappropriate initial contours. Besides, the limitation about time step limits the speed of the LBF model in the evolution of zero level set.

The Proposed Method
The proposed method is a fusion method which combines the advantages of the LBF model and the ADPLS method by taking both local and global intensity information into account. It is built to apply the two methods simultaneously To make the proposed method adjust adaptively to image information, the value of ( ) should consider two aspects. Firstly, the variable weighting coefficient should reflect the degree of intensity homogeneity. Secondly, the proposed method should use ADPLS method mainly to be fast and robust when the zero level set evolves in the intensity homogeneity area and use LBF model mainly to segment image accurately when the zero level set evolves near the objective boundary. According to these points above, the variable weighting coefficient ( ) is designed as follows: where > 0 is a constant and ( ) is the standard deviation of the pixel in a 3 × 3 neighborhood as a matter of experience, which is inversely proportional to the homogeneity of the image in the neighborhood. The more homogeneous the image's intensity level is, the smaller the value of ( ) is. With bigger weight, the ADPLS method will lead the evolution of level set. There are two advantages of the proposed method. First, with the global search ability, the proposed method can avoid being trapped into local minimums. Second, as the ADPLS method, which has a much larger time step, is much faster than the LBF model, the proposed method will be faster than the LBF model. On the other hand, the value of ( ) will be bigger when the pixel is near or on the edge. Accordingly, the weight of the LBF model will increase and be even larger than that of the ADPLS method. So the accuracy of the result will be better than the ADPLS method alone. As discussed above, the proposed method combines the advantages of the ADPLS method and the LBF model. This can be seen from a simple experiment of a simulated image shown in Figure 3. Figures  3(b) and 3(c) show the dominant force of the zero level set evolution on each pixel. As can be seen, the dominant force near the real edge that comes from the LBF model is marked blue (the borders of the image with 3-pixel width are set blue). On the other intensity-homogeneous area, the dominant force that comes from the ADPLS method is marked green.

Implementation of Level Set.
The implementation steps of the proposed algorithm in this paper can be described as follows.
Step 1. According to the given arbitrarily region in the image domain, initialize the level set function 0 . If the pixel ( , ) belongs to the region , then set 0 ( , ) = −2; otherwise set 0 ( , ) = 2.
Step 2. According to the image intensity and (10), we solve the weighting coefficient ( ) and initialize parameters. Set the variable of iterations iterNum = 0.
Step 3. If the iterNum is less than the maximum number of iterations iterMax, then repeat the following steps.
Step 4. Neumann boundary conditions are applied on the boundary of the image. According to the literature [5,7] and (4) and (7), get Δ ADPLS and Δ LBF and then obtain Δ by solving (9).
Step 5. Update the level set function using the formula = + Δ and set iterNum = iterNum + 1; then return to Step 3.

Experiments
The proposed method has been tested with both synthetic and medical images from different modalities. Unless otherwise specified, we use the following parameters in this paper. The parameters used in LBF model are Δ LBF = 0.1, 2 = 1.0, = 3.0, and the parameter of punishment = 0.002 * 255 * 255; the values of 1 and 2 depend on the actual image to be segmented which is described in the literature [14]. The parameters used in ADPLS method are Δ GIF = 1.0, 1 = 1.5, = 2.0, and the parameter of punishment = 0.2/Δ GIF , = 10; the value of should be small when dealing with simple images and be great with multilayer-complex images, which is referred to in the literature [7]. The constant coefficient in ( ) of (10) is set by experience. It is appropriate to set between 1 and 10 according to a lot of experiments. When the image is intensity inhomogeneous, should be small and be great on the contrary. We do all the experiments with Matlab code run on a Dell Optiplex 210L PC, with Pentium 4 processor, 3.0 GHZ, 1 GB RAM, with Matlab 6.5 on Windows XP. Figure 4 shows the results of the proposed method and the LBF model using the same synthetic image with the same initial contours in the first row and the second row,  respectively. Results in the first row show that the LBF model fails to segment the object correctly after 300 iterations and 26.922 seconds. However, satisfactory segmentations can be obtained by using the proposed method with 10 iterations and 0.672 seconds. In these experiments, we find that, no matter what the initial contours are, the proposed method can lead to better segmentations with less iterations and time costs than the LBF model. Figure 5 shows the results for an X-ray vessel image, which is a typical image with intensity inhomogeneity. In our experiments, the ADPLS method fails to segment the object correctly after 1000 iterations and 24.453 seconds. As can be seen from Figure 5 results (shown in Figure 5(c)). That weak part of the vessel boundaries can be segmented successfully using the proposed method. This demonstrates that, owing to combining LBF model, the proposed method is more accurate than the ADPLS method. Figure 6 shows the segmentation results for a brain MR image using the three methods mentioned above. Column 1 shows the initial contours, and Columns 2, 3, and 4 show the results of the LBF model, the ADPLS method, and the proposed method, respectively; they all have the same initial contours in column 1. Figure 7 gives the enlarged view of segmentation results on local region of row 1 in Figure 6. As to the segmentations in Figures 6 and 7, the stability, iterations, elapsed time, and accuracy are listed in Table 1. Regarding accuracy, the proposed method is close to the LBF model and much more accurate than the ADPLS method. With regard to elapsed time, the proposed method is close to the ADPLS method but much faster than the LBF model. With regard to algorithm stability, the proposed method and the ADPLS method are less sensitive to the initialization than the LBF model. Figure 8 shows the segmentation results for a synthetic image using the three methods mentioned above. Column 1 shows the initial contours, and Columns 2, 3, and 4 show the results of the LBF model, the ADPLS method, and the proposed method, respectively; they all have the same initial contours in Column 1. The extensive experimental results showed the superior performance of the proposed method over the state-of-the-art methods, in terms of both robustness and efficiency.

Conclusions
The fusion method proposed in the paper uses a variable weighting coefficient to combine both ADPLS method and LBF model. The proposed method allows the ADPLS method to be dominant force in intensity-homogeneous area of the image and the LBF model to be dominant force in intensity inhomogeneous area and objective boundaries, so it could take full use of these two algorithms' advantages, respectively, to complement each other. Compared with the LBF model, the fusion method is less susceptible to the initial contours and can always get good segmentation results and can be much faster. On the other hand, compared with the ADPLS, the proposed method can fully avoid boundary leakage and lack of segmentation with high speed when processing images with intensity inhomogeneity. Experiments show that the present method has superior comprehensive performance compared to ADPLS and LBF algorithm, respectively, in segmentation accuracy, speed, and stability.