Mathematical Problems in EngineeringVolume 2012 (2012), Article ID 474938, 10 pagesdoi:10.1155/2012/474938

Research Article

Incorporating a Local Binary Fitting Model into a Maximum Regional Difference Model for Extracting Microscopic Information under Complex Conditions

Chunyan Yao,1,2 Jianwei Zhang,² Min Chen,³ Qiu Guan,³ and Massimo Scalia⁴

¹Key Laboratory of E&M, Zhejiang University of Technology, Ministry of Education of China, 310014 Hangzhou, China²Department of Informatics, University of Hamburg, 22527 Hamburg, Germany³College of Computer Science, Zhejiang University of Technology, Hangzhou 310023, China⁴Department of Mathematics, Sapien za University of Rome, P.le A. Moro 2, 00185 Rome, Italy

Received 30 March 2011; Accepted 15 May 2011

Academic Editor: Shengyong Chen

Copyright © 2012 Chunyan Yao et al. 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.

This paper presents a novel region-based method for extracting useful information from microscopic images under complex conditions. It is especially used for blood cell segmentation and statistical analysis. The active model detects several inner and outer contours of an object from its background. The method incorporates a local binary fitting model into a maximum regional difference model. It utilizes both local and global intensity information as the driving forces of the contour model on the principle of the largest regional difference. The local and global fitting forces ensure that local dissimilarities can be captured and globally different areas can be segmented, respectively. By combining the advantages of local and global information, the motion of the contour is driven by the mixed fitting force, which is composed of the local and global fitting term in the energy function. Experiments are carried out in the laboratory, and results show that the novel model can yield good performances for microscopic image analysis.

1. Introduction

Observing and analyzing microscopic information has become one of the most important ways for exploring the microcosmic world for mankind. One typical issue is cell image analysis, which is at the forefront of clinics and bioengineering. In recent years, biological discovery and its translation into new clinical therapies have advanced rapidly. Pathologists use histopathological images of biopsy samples removed from patients and examine them under a microscope. While examining such images, a pathologist typically assesses the deviations in the cell structures and the change in the proportion between nucleolus and enchylema across the tissue under examination [1]. It is crucial to discover an automatic method to detect the contour of the cell and its nucleolus [2]. Cell detection is critical to disease diagnosis and prevention.

Dealing with cell images, there may exist many complex conditions such as large density of many cell objects, overlapping, and complex nucleus structures. This paper proposes a hybrid contour model for segmentation of the cell nucleolus and membranes. The proposed algorithm can handle cell overlapping and calculate a number of morphological parameters as well. By adopting the affinity propagation algorithm, a classifier can be designed for cell classification. Local and global intensity information is used as the driving forces of the active contour model under the guidance of the maximum regional difference. The local and global fitting forces ensure that local dissimilarities can be captured and global different areas can be segmented, respectively. This hybrid contour can capture the image details in the nucleus area relative to the traditional contour model, which is important to study the morphological nucleus parameters.

A tremendous amount of research work has been conducted in the area of microscopic image processing [3]. There are some challenges mentioned in these studies. The noise pollution is the most frequently occurring problem. It arises from staining the biopsy samples, and uneven distribution of the stain usually causes problems in processing the stained material. Besides noise, a strong shading effect is apparent in the images so that the cell images do not appear smooth. Accurate cell segmentation often precedes other analyses such as that of the morphologic process and behavior. Active contour models have been among the most successful methods facing the challenge mentioned above [4, 5]. Techniques based on active contour models have the potential to produce better estimates of cell morphologies.

The existing active contour models can be categorized into two classes: edge-based models [5–7] and region-based models [8–10]. The edge-based model directly uses intensity gradient information to attract the contour toward the object boundaries. Therefore, this kind of model has a worse performance for weak object boundaries since the cell image is very fuzzy due to low contrast at the location of the cell membrane. The region-based model aims to identify each region of interest by using a certain region descriptor. It guides the motion of the contour and is to some extent less sensitive to the location of initial contours [11, 12]. It is therefore more suitable for cell segmentation [13].

The Chan-Vese model [14] is one of the most popular region-based active contour models. This model has been successfully used for segmenting images. However, the regions of interest in images are often not statistically smooth and, therefore, the models are not applicable to cell images. Li et al. mentioned the intensity inhomogeneity as the biggest problem in low-signal-noise-ratio images such as biopsy samples and medical images [15].

In this paper, a novel model based on the region-based active contour model is proposed, which is able to detect the local unclear boundaries and segment the main differing regions successfully according to the principle of the maximum interregional difference. An energy function with two main terms, the local fitting term and the global fitting term, is defined in this model. The local fitting term produces a strong force to attract the contours and makes them stop at the boundaries although the object boundary is not clear or even blurry. The global fitting term ensures that the curve extracts the main object regions on the principle of the maximal regional difference. In addition, a strategic weight parameter, which is defined using the image intensity gradient information, is introduced to make the two forces work together as a hybrid model. This parametric model induces a mixed force comprising local and global information as driving force.

2. Maximum Regional Difference Model

Chan and Vese [14] proposed an active contour model that segments an image into two sets of possibly disjoint regions, by minimizing a simplified Mumford-Shah function. Assume that are the image domain and are the given image. Mumford and Shah consider segmentation as a problem of seeking an optimal contour that divides the image domain into two approximated piecewise-constant regions with intensities and . Denote as its boundary. Thus, the global data fitting term in the Chan-Vese model is defined as follows: where and represent the regions outside and inside the contour , respectively, and and are two constants that fit the image intensities outside and inside .

This model considers pixels within the same region as having the most similarity and makes up for the shortcomings of the edge detector. When the contour accurately represents the object boundary, the two fitting terms minimize the fitting energy value. In each segmented area, the clustered pixels’ mean value approximately equals and , respectively. Thus the fitting terms with respect to and are the driving forces that evolve the curve motion on the principle of inner region homogeneity.

Since the regional difference is the guideline in image segmentation, the interregional differences should be considered as the model’s driving force as follows:

This kind of region-based active contour model energy is characterized by the maximum dissimilarity between regions. To minimize the energy in (2.2) is the same as maximizing the difference between different regions. Equation (2.2) formulates the global instructive guidance term. This paper revises the guidance term which is defined as follows: where represents the entire image area. The curve evolution of the contour is defined as where and , denoting the areas outside of the contour and inside of the contour, respectively, and denotes the outward unit normal of . In order to show a stable value of (2.4) rather than an oscillating one when the contour divides the area with a small value, that is, to avoid the denominator in (2.4) from being zero, we approximately evolve the as follows: where represents the ratio between and and represents the ratio between and .

3. Incorporation of the Local Binary Fitting Model

Though the in (2.5) regards the regional difference as the guideline of curve evolution, the region-based model has little ability to deal with the intensity inhomogeneity or shading effect caused in the imaging process. The region-based model considers the global difference as the driving force rather than local detail information, and thus it needs a term of local fitting energy to improve the performance of segmentation in the inhomogeneous area.

Li et al. proposed a local binary fitting (LBF) model, which utilizes the local intensity information [15, 16]. An important parameter called the kernel function is introduced to denote the range of the local area. Additionally, and as the two spatially varying fitting functions are used to approximate the local intensities [17, 18]. In the LBF model, the local fitting term is defined as where is the Heaviside function and is a Gaussian kernel that decreases and approaches zero as increases. The function calculates the fitting degree around point , and can be considered as the center point of the local area. One of the greatest advantages of this LBF model is that it discovers the object boundaries more precisely than the Chan-Vese model.

The proposed method makes use of the advantages of the region-based model, as in (2.5), and the local binary fitting model by taking the local and global intensity terms into account. The entire energy model is defined as where is the strategic weight parameter, is defined in (3.1), and is defined in (3.2). The strategic weight parameter is defined as where represents the gradient information in image I and α is a positive constant.

This strategic weight parameter realizes the trade-off into a hybrid model [19, 20]. Working with the mixed energy restriction terms, the curve’s driving force appears strategically. When the local image area is flat or far away from the object boundaries, the value of is relatively small therefore the right side of (3.3) has a value approximating one. The hybrid model is guided to maximize the interregional difference. would gain a big value when the point is located in the region with intensity inhomogeneity where the object boundaries exist, thus the value of is close to zero. The model possesses a potential force to extract local boundaries. It combines the advantages of the local fitting model and the global regional difference model by the use of intensity gradients as its strategic parameter.

By minimizing the energy function in (3.2), we face solving a partially differential equation. The gradient descent algorithm is chosen to give a numerical solution with respect to the PDE. In addition, we introduce a fitting term mentioned in [5] to regularize the level set function for accurate computation by penalizing the deviation of the level set function from a signed distance function. It is defined as Besides, we introduce a length restriction term that has been used in the Chan-Vese model to smooth the zero level set contour,

Now, the entire model proposed in this paper is defined as The fitting functions and are given by and are constants defined as in [14].

4. Experiments and Results

4.1. Implementation Issues

To minimize the energy function with respect to defined in (3.6), we give the solution using a gradient descent algorithm which is defined as where, and are constants and is the Dirac delta function. and are given by

The iteration scheme is used by discretizing the PDE (4.1). The level set function ultimately sets the image object boundaries at . Furthermore, we set the parameters, , and. equals three in the expression of in the Gaussian kernel equals two.

4.2. Experimental Results

In this section, we use some cell images to demonstrate the mixed driving force in the proposed method for verifying its performance in detecting the contours of cell nucleolus and membrane. Figure 1 shows the process for segmenting neutrophil images. In Figure 1(a), the upper left picture depicts the initial contour and the Figures 1(b), 1(c), and 1(d) show how the curve evolves. It can easily establish the main cell nucleolus boundary and the weak object boundary of the cell membrane. The images resulting after 200 iterations of are shown in Figure 2.

(a)

(b)

(c)

(d)

Figure 1: Neutrophil segmentation and level set function .

Figure 2: Level set function after 200 iterations.

Due to its strong ability to capture local difference, the model can detect the weak boundary of the neutrophil membrane while it accurately distinguishes the nucleolus from the cytoplasm. The global driving force ensures that the model to segment the object has the maximum regional difference and that failure cases are avoided. Moreover, a comparison has also been carried out in our experiments with the Chan-Vese model and the LBF model, using a tumor lymphocyte cell image as test images. The results are shown in Figure 3. Figure 3(a) is the result by our proposed model. Figures 3(b) and 3(c) show the results of the Chan-Vese and LBF model, respectively. The tumor cell in an intensity inhomogeneous condition is detected by the zero level set function. The model successfully extracts the nucleolus, whose size relative to the area of the cell is a significant parameter in biology.

(a)

(b)

(c)

Figure 3: Comparison of different models.

Figure 4 illustrates another example for cell segmentation using the proposed model, where the parameters are, , , , and . Figure 5 shows the corresponding function after 200 iterations.

(a)

(b)

Figure 4: Another example for cell segmentation using the proposed model.

Figure 5: Corresponding function after 200 iterations.

5. Conclusion

The model proposed in this paper yields a relatively more desirable performance than the traditional models, although the results still feature over segmentation in the nucleolus area. However, this over segmentation problem can be overcome by a postmorphological process. For example, we find the locations where function and ( denotes the contour of the object) and substitute them with a binary value of zero and one, respectively. Then, we give a postprocess by setting a threshold value of the area size in the region of the neutrophil nucleolus area to eliminate the oversegmentation.

From the experimental results, it follows that the proposed hybrid model can well capture the weak edge boundary by using the local fitting terms, while it observes the rule of regional difference and uses it as guideline to form the fitting terms that drive the contour to accurately capture the boundaries. The experiments demonstrate that it performs well in segmenting the nucleolus and cell membrane.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (51075367, 60870002) and the Science and Technology Department of Zhejiang Province (R1110679, 2010R10006, and 2010C33095).

References

S.-J. Chern and P.-C. Huang, “On the existence of a weak solution of a half-cell model for PEM fuel cells,” Mathematical Problems in Engineering, Article ID 701096, 16 pages, 2010.

S. Y. Chen and Y. F. Li, “Determination of stripe edge blurring for depth sensing,” IEEE Sensors Journal, vol. 11, no. 2, pp. 389–390, 2011.

C. Gunduz, B. Yener, and S. H. Gultekin, “The cell graphs of cancer,” Bioinformatics, vol. 20, no. 1, pp. 145–151, 2004.

C. Li, J. Liu, and M. D. Fox, “Segmentation of external force field for automatic initialization and splitting of snakes,” Pattern Recognition, vol. 38, no. 11, pp. 1947–1960, 2005.

C. Li, C. Xu, C. Gui, and M. D. Fox, “Level set evolution without re-initialization: a new variational formulation,” in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, vol. 1, pp. 430–436, 2005.

A. Vasilevskiy and K. Siddiqi, “Flux maximizing geometric flows,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 24, no. 12, pp. 1565–1578, 2002.

R. Malladi, J. A. Sethian, and B. C. Vemuri, “Shape modeling with front propagation: a level set approach,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 17, no. 2, pp. 158–175, 1995.

S. Y. Chen, Y. F. Li, and J. Zhang, “Realtime structured light vision with the principle of unique color codes,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '07), pp. 429–434, Roma, Italy, April 2007.

A. Tsai, A. Yezzi, and A. S. Willsky, “Curve evolution implementation of the Mumford-Shah functional for image segmentation, denoising, interpolation, and magnification,” IEEE Transactions on Image Processing, vol. 10, no. 8, pp. 1169–1186, 2001.

L. Vese and T. 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.

S. Y. Chen, H. Tong, Z. Wang, S. Liu, M. Li, and B. Zhang, “Improved generalized belief propagation for vision processing,” Mathematical Problems in Engineering, Article ID 416963, 12 pages, 2011.

S. Y. Chen and Y. F. Li, “Vision sensor planning for 3-D model acquisition,” IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 35, no. 5, pp. 894–904, 2005.

M. Chen, S. Chen, and Q. Guan, “Hybrid contour model for segmentation of cell nucleolus and membranes,” in Proceedings of the 2nd International Conference on Biomedical Engineering and Informatics (BMEI '09), Tianjin, China, October 2009.

T. Chan and L. Vese, “Active contours without edges,” IEEE Transactions on Image Processing, vol. 10, no. 2, pp. 266–277, 2001.

C. Li, C. Y. Kao, J. C. Gore, and Z. Ding, “Implicit active contours driven by local binary fitting energy,” in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR '07), June 2007.

C. Li, C.-Y. Kao, J. C. Gore, and Z. Ding, “Minimization of region-scalable fitting energy for image segmentation,” IEEE Transactions on Image Processing, vol. 17, no. 10, pp. 1940–1949, 2008.

S. Y. Chen, “Cardiac deformation mechanics from 4D images,” Electronics Letters, vol. 43, no. 11, pp. 609–611, 2007.

C.-G. Zhu and R.-H. Wang, “Least squares fitting of piecewise algebraic curves,” Mathematical Problems in Engineering, Article ID 78702, 11 pages, 2007.

M. Li, “Fractal time series—a tutorial review,” Mathematical Problems in Engineering, Article ID 157264, 26 pages, 2010.

C. Cattani, “Shannon wavelets for the solution of integrodifferential equations,” Mathematical Problems in Engineering, Article ID 408418, 22 pages, 2010.