A Framework for White Blood Cell Segmentation in Microscopic Blood Images Using Digital Image Processing
© Sadeghian et al. 2009
Received: 14 April 2009
Accepted: 21 May 2009
Published: 11 June 2009
Evaluation of blood smear is a commonly clinical test these days. Most of the time, the hematologists are interested on white blood cells (WBCs) only. Digital image processing techniques can help them in their analysis and diagnosis. For example, disease like acute leukemia is detected based on the amount and condition of the WBC. The main objective of this paper is to segment the WBC to its two dominant elements: nucleus and cytoplasm. The segmentation is conducted using a proposed segmentation framework that consists of an integration of several digital image processing algorithms. Twenty microscopic blood images were tested, and the proposed framework managed to obtain 92% accuracy for nucleus segmentation and 78% for cytoplasm segmentation. The results indicate that the proposed framework is able to extract the nucleus and cytoplasm region in a WBC image sample.
KeywordsWhite blood cell segmentation Active contours Snake algorithm Zack thresholding
White blood cells (WBC) or leukocytes play a significant role in the diagnosis of different diseases, and therefore, extracting information about that is valuable for hematologists. In the past, digital image processing techniques have helped to analyze the cells that lead to more accurate, standard, and remote disease diagnosis systems. However, there are a few complications in extracting the data from WBC due to wide variation of cells in shape, size, edge, and position. Moreover, since illumination is imbalanced, the image contrast between cell boundaries and the background varies depending on the condition during the capturing process.
This study is focusing on WBC segmentation using L2 microscopic images. Our goal is to segment the WBC nucleuses and cytoplasm using a framework that has been developed using digital image processing. The use of image processing techniques have developed rapidly in the last few years, to the point where hematologists can use blood images to automatically process blood slides for the first screening in detecting diseases. These techniques can help to find cell counts in human blood automatically and also can provide information about ratio of nucleus versus cytoplasm to identify and classify different types of WBCs such as neutrophil, basophil, lymphocyte, etc. Therefore, in this paper, we present a proposed framework that consists of several methods that integrates together for nucleus segmentation and cytoplasm extraction.
Many works have been conducted in the area of general segmentation methods. Among the common segmentation methods are edge and border detection, region growing, filtering, mathematical morphology, and watershed clustering. Ritter et al.  presented a fully automatic method for segmentation and border identification of all objects that do not overlap the boundary in an image taken from a peripheral blood smear slide. In their work, pale tips of protuberances are lost. Ongun et al.  did segmentation by morphological preprocessing followed by the snake-balloon algorithm. Jiang et al.  proposed a WBC segmentation scheme on color space images using feature space clustering techniques, scale-space filtering for nucleus extraction, and watershed clustering for cytoplasm extraction. Leyza et al.  used morphological operators and examined the scale-space properties of toggle operator to improve segmentation accuracy. Scotti  presented the automatic morphological method that is based on the morphological analysis of WBCs. Their proposed system extracts the morphological indexes (lymphocytes). Kumar et al.  used teager energy operator for segmentation, nucleus based on the edges, which are detected effectively by teager energy operator but it required at least a weak edge to exist between red blood cell (RBC) and the background. For cytoplasm segmenting, they used a simple morphological method. Cseke introduced multi-step segmentation scheme , which implements the automatic thresholding method proposed by Otsu .
The remainder of this paper is organized as follows. In Section 2, segmentation algorithms and the framework are explained. In Section 3, results obtained by the proposed framework are presented, and finally, the conclusions are drawn in Section 4.
2. Proposed Framework
2.1. Nucleus Segmentation
Nuclei have variable shapes in different kinds of leukocytes. Finding a significant method for shape modeling and segmenting the nucleus has always been a challenge for scientists. Among segmentation methods, active contour models (snakes) have gained a lot of attention recently . Snakes are deformable curves that can move and change their shapes to deform to boundaries of objects in an image. Curves are defined within an image domain and can move under the influence of internal forces within the curve itself and external forces derived from the image data. The internal and external forces are defined in a way that the snake conforms to an object boundary or other desired features within an image [9, 10].
Where α and β are weighting parameters that control the snake's tension and rigidity, respectively. x'(s) and x"(s) denote the first and second derivatives of x(s) with respect to s. The external energy function E ext is derived from the image so that it takes on its smaller values at features of interest, such as boundaries.
Ridge pixels are thresholded using two thresholds T 1 and T 2 with T 1⟨T 2. Ridge pixels with values between T 1 and T 2 are weak edge pixels, and those with values larger than T 2 are strong edge pixels.
Edges segments in T 2 are linked to form continuous edges. To do so, each segment in T 2 is traced to find its end, and its neighbors in T 1 are searched to find any edge segment in T 1, which can bridge the gap until reaching another edge segment in T 2.
By this edge detection method, central connected object boundaries that represent the nucleus are clearly obtained. In next step, GVF of the images were calculated based on Eq. 4 and used as internal and external forces to guide snakes to deform to nucleus boundary edges. Nucleus is the connected boundary in image and has been filled up  by following instruction to have a clear segmented nucleus:
Where R is a reconstruction of f c from f m.
2.2. Cytoplasm Segmentation
By subtracting the segmented nucleus from the original sub-image, we will obtain the cytoplasm, RBC, and the background. Most of the time, RBCs appear in the image border. Looking at the gray level intensities, the cytoplasm and two other components are having almost uniform areas. Therefore, it justifies the need for segmenting these uniform components using thresholding techniques.
3. Results and Discussion
This section is to assess the performance of the proposed WBC segmentation scheme. In our experiment, 20 blood images from acute leukemia cases type L2 were captured using Microscope Olympus BX51. In more details about this digital microscopy acquisition, whole system is soft imaging system with AnalySIS software. Its camera is CC-12, and magnification used ×400 for pictures. Amount of fields per slide acquired is quiet random, which is 2–3 field/slide and total slides are 10.
GVF deformable contour was done with suitable iterations, and the final results are shown in Figure 3d. Snake algorithm finds the connected boundary that is detected in Figure 3c and it selects the nucleus. The result has been shown in Figure 3e. The connected boundaries have been filled up and shown in Figure 3f representing the nucleus of the WBC.
Results present the advantages of our method instead of others. In nucleus segmentation, we used snake algorithm that is not related to size and color of nucleoli because there are various shapes of nucleus in a different kind of white blood cells. So, it gives high accuracy result in segmenting nucleus in any type of WBCs and in any capture illumination that cause different color space in images. Also, in a cytoplasm segmentation method in which we used the thresholding technique, background is completely segmented from other components. And, based on the difference between RBC color in blood image and cytoplasm area, thresholding value is easily segmented in the cytoplasm part. But, we prefer to use a sub-image that contains individual WBC to get a better result. This method is very simple with high speed and trustable accuracy.
By the results we got from this framework, hematologists can decide on the types of WBC and its maturation and also potentially can calculate the amount of cells in specific blood smear and finally in whole body blood. For some of the diseases like leukemia, knowledge about amount of WBCs and also their maturation is very important. In hematology science, information about size and volume of nucleus and cytoplasm is profitable. Our method gives useful information about WBC maturation status by finding the dimension of WBC components, nucleus and cytoplasm.
The framework has been done on sub-images to have easier implementation; this calls the major limitation in our method. In blood image, there are similar color scales in WBCS with some other blood particles that cause a big error in thresholding method for cytoplasm segmentation, so we individuate the WBCs in sub-image to reduce the errors. In the future, we will try to segment sub-images automatically to have a WBC segmentation process that is fully automated. The method has been applied to 20 images. We can calculate the percentage of the accuracy by evaluating WBCs' component segmentation base on the comparison of our own method and manual segmentation. We get an average accuracy of 92% for nucleus segmentation and 70% for cytoplasm segmentation. Since the cytoplasm segmentation process depends on the result of the nucleus segmentation, hence, the 8% (100 - 92%) error yielded by the nucleus segmentation affects directly on the cytoplasm extraction accuracy. Leaving this fact behind, the accuracy of the cytoplasm alone is 78%.
On average, there are about 55 WBCs in a typical blood smear image as explained before for our acute leukemia, type L2, sample blood images. Based on the results (92% accuracy for nucleus and 78% for cytoplasm), after applying the method on a sample image, there is a chance that some parts of the cytoplasm and the nucleus be missed in each WBC. This comprises the 8% and 22% errors and may show its effect on estimating the ratio of nucleus and cytoplasm only. The results show significant accuracy to be used for further analysis of blood images on detection of acute lymphoblastic leukemia or any other diseases related to WBCs.
This paper has demonstrated a proposed framework for segmenting white blood cells using integration of concepts in digital image processing. The proposed scheme has two parts: The nucleus segmentation part is based on morphological analysis, and the cytoplasm segmentation is based on pixel-intensity thresholding. The results show that the proposed method is able to yield 92% accuracy for nucleus segmentation and 78% for cytoplasm segmentation.
- Ritter N, Cooper J: Segmentation and border identification of cells in images of peripheral blood smear slides. Proceedings of the Thirtieth Australasian Conference on Computer Science. 2007, 62: 161-169.Google Scholar
- Ongun G, Halici U, Leblebicioglu K, Atalay V, Beksac M, Beksac S: Feature extraction and classification of blood cells for an automated differential blood count system. Neural Networks. Proceedings. IJCNN '01. International Joint Conference on. 2001, 4: 2461-2466.Google Scholar
- Jiang Kan, Liao Qing-Min, Dai Sheng-Yang: A novel white blood cell segmentation scheme using scale-space filtering and watershed clustering. Machine Learning and Cybernetics, 2003 International Conference on. 2003, 5: 2820-2825.View ArticleGoogle Scholar
- Dorini LB, Minetto R, Leite NJ: White blood cell segmentation using morphological operators and scale-space analysis. SIBGRAPI '07: Proceedings of the XX Brazilian Symposium on Computer Graphics and Image Processing. 2007Google Scholar
- Scotti F: Automatic morphological analysis for acute leukemia identification in peripheral blood microscope images. 2005 IEEE International Conference on Computational Intelligence for Measurement Systems and Applications, 2005. 2005, CIMSA, 96-101.View ArticleGoogle Scholar
- Kumar BR, Joseph DK, Sreenivasc TV: Teager energy based blood cell segmentation. Digital Signal Processing. DSP 2002 14th International Conference on. 2002, 2: 619-622.Google Scholar
- Cseke I: A fast segmentation scheme for white blood cell images. Pattern Recognition. Conference C: Image, Speech and Signal Analysis, Proceedings. 11th IAPR International Conference on. 1992, III: 530-533.Google Scholar
- Otsu N: A threshold selection method from gray-level histograms. Automatica. 1975, 11: 285-296. 10.1016/0005-1098(75)90044-8.View ArticleGoogle Scholar
- Kass M, Witkin A, Terzopoulos D: Snakes: active contour models. Comput Graph Image Process. 1988, 1: 321-331.Google Scholar
- Xu C, Prince JL: Snakes, shapes, and gradient vector flow. IEEE Trans Image Process. 1998, 7: 359-369. 10.1109/83.661186.View ArticlePubMedGoogle Scholar
- Xu C, Prince J: Gradient vector flow: a new external force for snakes. 1997 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 1997. Proceedings. 1997, 66-71.Google Scholar
- Canny J: A computational approach to edge detection. IEEE Trans. Pattern Anal Mach Intell. 1986, 679-698. 10.1109/TPAMI.1986.4767851.Google Scholar
- Gonzalez RC, Woods RE, Eddins SL: Digital image processing using MATLAB. 2003, Prentice-Hall, Upper Saddle RiverGoogle Scholar
- Roberts LG: Machine perception of 3-D solids. Optical and Electro-Optical Information Processing. 1965, 159-197.Google Scholar
- Sobel I: Neighborhood coding of binary images for fast contour following and general array binary processing. Comput Graph Image Process. 1978, 8: 127-135. 10.1016/S0146-664X(78)80020-3.View ArticleGoogle Scholar
- Prewitt JMS: Object enhancement and extraction. Picture Processing and Psychopictorics. 1970, 75-149.Google Scholar
- Zack G, Rogers W, Latt S: Automatic measurement of sister chromatid exchange frequency. J Histochem Cytochem. 1977, 25: 741-753.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.