Abstract and Applied Analysis

Volume 2013 (2013), Article ID 131207, 8 pages

http://dx.doi.org/10.1155/2013/131207

Research Article

Construction of Target Controllable Image Segmentation Model Based on Homotopy Perturbation Technology

Shu-Li Mei

College of Information and Electrical Engineering, China Agricultural University, Postbox 53, East Campus, 17 Qinghua Donglu Road, Haidian District, Beijing 100083, China

Received 31 December 2012; Revised 4 January 2013; Accepted 6 January 2013

Academic Editor: Lan Xu

Copyright © 2013 Shu-Li Mei. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Based on the basic idea of the homotopy perturbation method which was proposed by Jihuan He, a target controllable image segmentation model and the corresponding multiscale wavelet numerical method are constructed. Using the novel model, we can get the only right object from the multiobject images, which is helpful to avoid the oversegmentation and insufficient segmentation. The solution of the variational model is the nonlinear PDEs deduced by the variational approach. So, the bottleneck of the variational model on image segmentation is the lower efficiency of the algorithm. Combining the multiscale wavelet interpolation operator and HPM, a semianalytical numerical method can be obtained, which can improve the computational efficiency and accuracy greatly. The numerical results on some images segmentation show that the novel model and the numerical method are effective and practical.

1. Introduction

In general, choosing different parameters in the most common image segmentation methods usually leads to different image segmentation results [1]. In other words, the object segmentation results are uncontrollable by the common methods. To solve the problem, one of the most common strategies is choosing thresholds using prior knowledge or analyzing the distribution of gray values of an image with the gray value histogram. Another method is image enhancement, which can often destroy the contour of the objects.

The variational method on image segmentation is a new image processing technology, which processes lots of better properties in processing medical images such as MRI and CI [2]. In this method, the pictures are taken as continuous energetic fields, and so the corresponding information in digital images such as gradient, divergence, and the curvature of the object contour can be viewed as the differential operators embedded in the variational model on image processing. The traditional complicated image processing such as denoising with texture preserving and exact segmentation can be done by this model. The outstanding work of this field is the energy function for image segmentation proposed by Mumford and Shah, which has been widely used, and its mathematical properties are well analyzed. This is a general approach on image segmentation, where it is assumed that objects can be characterized by smooth surfaces or volumes in three dimensions. In order to solve the Mumford-Shah model with the Euler-Lagrange method, a simplified model was deduced by Chan and Vese, in which the Euclid length was employed instead of the Hausdorff length [3]. So, the simplified model is also called Chan-Vese model. Similar to other image segmentation methods, Chan-Vese model cannot identify the object as well. Multilevel set approach for solving C-V model can segment all the objects in a picture. But it will obviously lead to oversegmentation [4].

In many cases, the purpose of the image segmentation is to get one special single object instead of all the objects in a multiobject image. Therefore, the purpose of this paper is to construct a target controllable image segmentation model based on the basic idea of homotopy perturbation technology (HPM). Using the variational method, the optimal solution of the energy function can be expressed as a nonlinear partial differential equation. So, another task of this study is to construct an effective numerical method on nonlinear PDEs by combining the multiscale wavelet interpolation operator and the homotopy perturbation method. The homotopy perturbation method (HPM) proposed by He [5, 6] is constantly being developed and applied to solve various nonlinear problems by He [7–15] and by others [16–20]. The better improvement is adding an auxiliary parameter into the homotopy equation, which is helpful to eliminate the secular term in the perturbation solution. This can improve the rate of convergence greatly. Unlike analytical perturbation methods, HPM does not depend on small parameter which is difficult to find. The variational iteration method was another simple and effective method for nonlinear equations proposed by He [21–26], which can provide analytical approximations to a rather wide class of nonlinear equations [27–33] without linearization, perturbation, or discretization which can result in massive numerical computation. In order to solve the nonlinear PDEs, it is necessary to introduce the wavelet numerical algorithm [34–37] into HPM.

2. Construction of Target Controllable Image Segmentation Model

In order to solve the Mumford-Shah model with the Euler-Lagrange method, a simplified model was deduced by Chan and Vese, in which the Euclid length was employed instead of the Hausdorff length. This simplified model can also be called the Chan-Vese model, which can be expressed as follows:

where and are positive constants and and are the average gray level values inside () and outside () of the object contour, respectively. denotes the image to process, is the length of the object contour, and is the weight parameter. According to the level set method, the contour curves of the objects should be embedded into the level set function as follows:

Then, the level set-based C-V model can be rewritten as follows:

Using the variational method, the PDEs with respect to the variable can be obtained as follows:

Obviously, is the curvature of the level set function , and is used to constrain the growth of the level set function.

The solution of (4) is the level set function at time . The zero level set is the object contour curve, which can be obtained by solving .

In the following, what we are talking about is how to construct the target controllable image segmentation model based on the basic idea of HPM. It is easy to understand that the function of the curvature in C-V model is just to preserve the smoothness of the object contour. Neglecting the curvature in (4), the simplified model can be obtained as follows:

In solving the C-V model with HPM and iteration method, the average gray level values inside and outside of the contour curves and vary with the evolution of the level set function. This evolution will end up with that the contour curve coincides with the object boundary. Then, and become constants, and the right hand of (5) should equal zero; that is,

In general, and are constant, which are correlated to and obviously.

It is easy to understand that the segmentation results and the values of are in one-to-one correspondence with each other. So, the object segmentation can be controlled by the value of (). Let , , and substitute , into (5), we can obtain

Let

It is obviously that is the necessary condition for the functional extremum problem about C-V model. The solution of the necessary condition is

In the end of the image segmentation processing, the gray level value of the pixel inside the object contour equals ; that is, . Then, the relation between the parameter and the average gray level value of the image can be expressed as

In the end of the image segmentation procedure, the final and should be coincident with average gray level values inside and outside of the segmentation target, respectively. They can be determined in advance by the priori knowledge. But in the beginning of the image segmentation processing, and are the average gray level values inside and outside of the zero level set, respectively. They are determined by the position of the level set function, which is random in most cases. It is easy to understand that there is a continuous map between the two cases, that is, the connection between the two cases can be set up by the HPM. In other words, the parameter in (10) can be taken as the homotopy parameter; then, a linear homotopy function for (4) can be constructed as

where the homotopy parameter , is determined by the final and based on the priori knowledge in advance and is a weight parameter. is an auxiliary parameter, which can be identified by eliminating the secular term in the perturbation analytical solution.

Equation (11) is the target controllable image segmentation model. It should be pointed out that the auxiliary parameter appearing in this model is set to zero, as there is no any secular term in the perturbation analytical solution.

3. HPM on Nonlinear System Based on the Multilevel Wavelet Analysis

3.1. Wavelet Numerical Discretization Schemes on C-V Model

The definition domain of the image is defined as , which should be divided evenly into ( is the level number) subdomains according to the wavelet collocation method. The connection nodes between two adjoining subdomains are the discretization points defined as , where

In addition, denotes the multiscale wavelet function and the corresponding th and th derivatives with respect to and , respectively. The level set function and the corresponding derivative function can be descretized as follows:

where and are constants, which denote the wavelet scale number and the maximum of the scale number, respectively. , , and are the wavelet coefficients at the discretization point .

According to above definitions, the curvature of the level set can be expressed approximately as

Substituting (14) and (13) into (11), we obtain

Obviously, (15) is a nonlinear ordinary differential equation.

3.2. HPM on Discretization Format of C-V Model

There are various ways to construct a homotopy function. For (15), a linear homotopy function can be constructed as

where

where is the weight parameter. According to the perturbation theory, the solution of (16) can be expressed as the power series expansion of

Substituting (18) into (16) and rearranging based on powers of -terms, we have

where and are functions with respect to and , respectively. It is easy to identify the homotopy parameter as

Thus, based on the definition of the Taylor series, can be identified as , and can be identified as

Substituting and into (18) and assuming , the numerical solution of (16) can be obtained subsequently:

Then, the wavelet coefficient can be obtained as follows:

where

is the quasi-Shannon wavelet function; that is,

Substituting the three wavelet coefficients into (16), can be obtained; then we can obtain as follows:

At last, we can obtain the image segmentation result expressed in the level set as follows:

4. Numerical Experiences and Discussion

In this section, we take some multiobject images as examples to illustrate the efficiency of the target controllable image segmentation model compared with the C-V model. The original image showed in Figure 1 consisted of three geometrical solid objects. The color of the background is white, and whole area is 1. The gray level values and the areas of the three objects are showed in Table 1. The image segmentation aims to get the circular and the rectangular objects. In other words, we want to take the black elliptical object as the background.

Table 1: The gray level values and areas of all objects.

Ellipse
Circle
Squareness
Background

Gray level value
0
163
7
255

Area
0.0568
0.0723
0.0733
0.7976

(a) Original image

(b) 10 iteration times

(c) 50 iteration times

(d) 500 iteration times

(e) 750 iteration times

(f) 1000 iteration times

Figure 1: Multiobject image segmentation at (, ).

The segmentation results of C-V model are showed in Figure 1. With the increasing of the iteration times, the rectangular object becomes a part of the background gradually instead of the elliptical object, which does not meet our requirement obviously.

The segmentation results with the target controllable model are showed in Figure 2. The final   , and   . It should be pointed out that the final and can be obtained by priori knowledge in most cases. At the beginning of the segmentation, all the three objects are obtained as the foreground. With the increasing of the iteration times, the black object is gradually pushed into the background and out of the object region.

(a) 10 iteration times

(b) 1000 iteration times

Figure 2: Multitarget image segmentation with target controllable model.

Figure 3 is an enlarged local image of the locust body cavity. The objects have an irregular shape with a slightly serrated border which can introduce over an insufficient segmentation. So it is difficult to segment with other methods. Using the target controllable model, we can get the right object easily. This example shows that the novel model and the corresponding numerical method are practical. Indeed, the novel model has been used to segment the locust’s coelom images (Figure 4).

(a) C-V model

(b) Target controllable model

Figure 3: Special object segmentation.

Figure 4: Locust’s coelom image segmentation results by the target controllable model.

5. Conclusions

C-V model is a kind of the modified Mumford-Shah model which has been widely used in medical images, and its mathematical properties are well analyzed. But the segmentation result is usually uncontrollable. The target controllable image segmentation model proposed in this paper is based on the idea of HPM. The numerical experiences show that the novel model and the corresponding numerical algorithm are effective and practical. It meets the requirement of the medical images segmentation.

Acknowledgment

This work is supported by the National Natural Science Foundation of China under Grant no. 41171337 and National Key Technologies R & D Program of China under Grant no. 2012BAD35B02.

References

O. Wirjadi, “Survey of 3rd image segmentation methods,” Berichte Des Fraunhofer ITWM 123, 2007.

L. A. Vese and T. F. Chan, “A multiphase level set framework for image segmentation using the Mumford and Shah model,” International Journal of Computer Vision, vol. 50, no. 3, pp. 271–293, 2002.

T. F. Chan, S. Esedoglu, and M. Nikolova, “Algorithms for finding global minimizers of image segmentation and denoising models,” UCLA Report, 2004.

T. Goldstein, X. Bresson, and S. Osher, “Geometric applications of the split Bregman method: segmentation and surface reconstruction,” Journal of Scientific Computing, vol. 45, no. 1–3, pp. 272–293, 2010.

J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

J.-H. He, “Asymptotology by homotopy perturbation method,” Applied Mathematics and Computation, vol. 156, no. 3, pp. 591–596, 2004.

J.-H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003.

J. H. He, “Limit cycle and bifurcation of nonlinear problems,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 827–833, 2005.

J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005.

J.-H. He, “Periodic solutions and bifurcations of delay-differential equations,” Physics Letters A, vol. 347, no. 4–6, pp. 228–230, 2005.

J.-H. He, “The homotopy perturbation method nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004.

J. H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 2, pp. 207–208, 2005.

J.-H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012.

J.-H. He, “New interpretation of homotopy perturbation method,” International Journal of Modern Physics B, vol. 20, no. 18, pp. 2561–2568, 2006.

L. Cveticanin, “Homotopy-perturbation method for pure nonlinear differential equation,” Chaos, Solitons and Fractals, vol. 30, no. 5, pp. 1221–1230, 2006.

S. Abbasbandy, “Application of He's homotopy perturbation method for Laplace transform,” Chaos, Solitons and Fractals, vol. 30, no. 5, pp. 1206–1212, 2006.

M. Rafei and D. D. Ganji, “Explicit solutions of Helmholtz equation and fifth-order KdV equation using homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 3, pp. 321–328, 2006.

A. M. Siddiqui, R. Mahmood, and Q. K. Ghori, “Thin film flow of a third grade fluid on a moving belt by he's homotopy perturbation method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 7–14, 2006.

A. M. Siddiqui, M. Ahmed, and Q. K. Ghori, “Couette and poiseuille flows for non-newtonian fluids,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 15–26, 2006.

J.-H. He, “Variational iteration method: a kind of non-linear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.

J.-H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 2-3, pp. 115–123, 2000.

J.-H. He and X.-H. Wu, “Construction of solitary solution and compacton-like solution by variational iteration method,” Chaos, Solitons & Fractals, vol. 29, no. 1, pp. 108–113, 2006.

J.-H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 1-2, pp. 57–68, 1998.

J.-H. He, “Approximate solution of nonlinear differential equations with convolution product nonlinearities,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 1-2, pp. 69–73, 1998.

J.-H. He, “Asymptotic methods: the next frontier towards nonlinear science,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 1907–1908, 2011.

M. A. Abdou and A. A. Soliman, “Variational iteration method for solving Burger's and coupled Burger's equations,” Journal of Computational and Applied Mathematics, vol. 181, no. 2, pp. 245–251, 2005.

A. A. Soliman, “A numerical simulation and explicit solutions of KdV-Burgers' and Lax's seventh-order KdV equations,” Chaos, Solitons and Fractals, vol. 29, no. 2, pp. 294–302, 2006.

E. M. Abulwafa, M. A. Abdou, and A. A. Mahmoud, “The solution of nonlinear coagulation problem with mass loss,” Chaos, Solitons & Fractals, vol. 29, no. 2, pp. 313–330, 2006.

S. Momani and S. Abuasad, “Application of He's variational iteration method to Helmholtz equation,” Chaos, Solitons & Fractals, vol. 27, no. 5, pp. 1119–1123, 2006.

N. Bildik and A. Konuralp, “The use of variational iteration method, differential transform method and adomian decomposition method for solving different types of nonlinear partial differential equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 65–70, 2006.

Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 27–34, 2006.

J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.

Z. Wan-Xie, “On precise integration method,” Journal of Computer and Applied Mathematics, vol. 163, no. 1, pp. 59–78, 2004.

D.-C. Wan and G.-W. Wei, “The study of quasi wavelets based numerical method applied to Burgers' equations,” Applied Mathematics and Mechanics, vol. 21, no. 10, pp. 991–1001, 2000.

G. W. Wei, “Quasi wavelets and quasi interpolating wavelets,” Chemical Physics Letters, vol. 296, no. 3-4, pp. 253–258, 1998.

S.-L. Mei, C. J. Du, and S. W. Zhang, “Asymptotic numerical method for multi-degree-of-freedom nonlinear dynamic systems,” Chaos, Solitons and Fractals, vol. 35, no. 3, pp. 536–542, 2008.