25 July 2017 Automated three-dimensional morphology-based clustering of human erythrocytes with regular shapes: stomatocytes, discocytes, and echinocytes
Author Affiliations +
We present unsupervised clustering methods for automatic grouping of human red blood cells (RBCs) extracted from RBC quantitative phase images obtained by digital holographic microscopy into three RBC clusters with regular shapes, including biconcave, stomatocyte, and sphero-echinocyte. We select some good features related to the RBC profile and morphology, such as RBC average thickness, sphericity coefficient, and mean corpuscular volume, and clustering methods, including density-based spatial clustering applications with noise, k-medoids, and k-means, are applied to the set of morphological features. The clustering results of RBCs using a set of three-dimensional features are compared against a set of two-dimensional features. Our experimental results indicate that by utilizing the introduced set of features, two groups of biconcave RBCs and old RBCs (suffering from the sphero-echinocyte process) can be perfectly clustered. In addition, by increasing the number of clusters, the three RBC types can be effectively clustered in an automated unsupervised manner with high accuracy. The performance evaluation of the clustering techniques reveals that they can assist hematologists in further diagnosis.



Blood cells have different functionalities in the human body and tissue. Red blood cells (RBCs), or erythrocytes, are the most abundant among blood cells. An erythrocyte is a discoid cell with a thick rim and thin sunken center. The main functions of erythrocytes are to absorb oxygen from the lungs, release it into tissues during circulation, and transport carbon dioxide from the tissues to the lungs. The biconcave shape of erythrocytes (doughnut-like) allows them to squeeze through capillaries that are smaller than an RBC. However, different blood abnormalities at different stages alter the original bioconcave shape of erythrocytes into different morphologies.12.3.4 During blood storage in blood banks, the shape of erythrocytes changes from biconcave to flat-disk, and then to sphero-echinocyte when the storage time exceeds a few weeks. It has been proven that the transfusion of damaged RBCs can cause severe problems to body tissue and, in some cases, may lead to death.

A typical human erythrocyte RBC has a diameter of 6.4 to 7.76  μm10 and thickness of 2 to 3  μm at its thickest point.11 The surface area of adult mature RBCs is 140  μm2. Under some circumstances, a mature RBC undergoes deformation into different shapes, such as echinocytes, stomatocytes, spherocytes, elliptocytes, acanthocytes, burr cells, and schizocytes among others.12 There are cases in which different types of RBCs may exhibit similar characteristics, such as mean corpuscular volume (MCV), with tiny differences in surface area and, in some cases, constant surface area and similar shapes, which make them difficult to distinguish easily. Therefore, RBC clustering suffers from common characteristics among different kinds of RBCs. Thus, different RBCs may sometimes be categorized in the same group and result in a significant misclassification. In a conventional RBC investigation, a hematologist manually counts and classifies the cells with assistance from a microscope; this is a procedure that is tiresome, time-consuming, and susceptible to error. More specifically, the accuracy of the counting and diagnosis is affected by subjective circumstances, such as experience and fatigue, due to human exhaustion.1314.15.16 Many efforts have been made in RBC studies to detect abnormalities in samples before transfusion to the patient to prevent future disorders caused by the malfunction of RBCs.1718.19

According to the above discussion, we believe the utilization of clustering techniques can provide us with reasonable results. Therefore, because of the irregularly shaped groups of RBC types distributed in feature space, the density-based clustering method can be effective in this case and for overlapped RBC data. Fuzzy clustering methods are efficient unsupervised methods that can be applied in this field for regularly shaped data points.20,21

The idea of data clustering is based on the concept that the human brain processes information as patterns rather than numerical entities. Clustering is a process that groups data observations that are similar to each other, whereas data observations of different clusters are not similar. In supervised classification methods, the data need to be labeled before applying the classification method. As RBC data from the initial phase image are unlabeled, the unsupervised method needs to be used to label and cluster them automatically. In the unsupervised method, the goal is to intrinsically group unlabeled data without predefined data groups.2021.22.23

In conventional two-dimensional (2-D) microscopic imaging techniques, it is difficult to detect the three-dimensional (3-D) shape of erythrocytes; thus, the overall performance is not acceptable. However, digital holographic microscopy (DHM) is capable of imaging semitransparent or transparent biological cells and provides quantitative detailed information about the cell structure and its contents at a single-RBC level. In addition, it is a noninvasive and label-free method. Therefore, samples can remain untouched for further investigations. The quantitative phase image (QPI) obtained by DHM enables us to measure the 3-D properties of RBCs, which include the volume, surface area, projected surface area (PSA), and dry mass of the biological cell.10,11,24,25

In this study, we apply several clustering methods to cluster the different shapes of RBCs. Several RBC samples with three major morphologies, biconcave, stomatocyte, and sphero-echinocyte, are visualized by the DHM technique and are combined together. RBCs are obtained from the reconstructed phase image from the DHM technique using the watershed segmentation algorithm.24,25 After feature extraction, similar to our previous work,11 we select some good features that can efficiently discriminate between the RBC types and evaluate the clustering power of the selected features against 2-D features only. To decrease the dimensions of the features dataset, principal component analysis (PCA) is applied to the extracted features, and only three PCAs are retained. Our experimental results reveal that three PCAs can represent 90% of the entire variance. This can help us reduce the problem, enhance the clustering speed, and make the solution more efficient. In addition, the PCA technique is useful for visualizing better the 2-D and 3-D space. Eventually, several clustering methods, including density-based spatial clustering applications with noise (DBSCAN), k-medoids, and k-means clustering are applied to the dimension-reduced 2-D and 3-D features, and the clustering performance is evaluated against the 2-D features. Our experimental results show that the combination (2-D and 3-D) of features can obtain high-accuracy clustering results against 2-D features in the automated clustering of RBCs with regular shapes.

The rest of this paper is organized as follows. Section 2 is dedicated to explaining the schematic of the off-axis DHM and RBC preparation process. The image segmentation technique used in this experiment for extracting RBC samples from QPIs is briefly discussed in Sec. 3. The feature extraction process in this experiment is described in Sec. 4. DBSCAN is discussed in detail in Sec. 5. The k-medoids and k-means clustering methods are discussed in depth in Secs. 6 and 7, respectively. The experimental results and discussion on the accuracy ratio of the combination of features against the 2-D features are explained in Sec. 8. Finally, we evaluate clustering results using silhouette index (SI) in Sec. 9. We conclude the paper in Sec. 10.


Off-Axis Digital Holographic Microscopy and Red Blood Cell Preparation


Off-Axis Digital Holographic Microscopy

The off-axis DHM system uses a laser diode source of wavelength λ=682  nm. The laser beam is divided into two waves, the object wave and reference wave. The object wave passes through the RBC sample and gets diffracted and magnified by a microscope objective (magnification: 40× and numerical aperture: 0.75); then, it interferes with the reference wave in the off-axis geometry. The interference pattern between the object and reference waves is recorded via a charge-coupled device. The QPI of the RBCs is reconstructed from the recorded interference pattern using a specific numerical algorithm.26 A schematic of the off-axis DHM system is shown in Fig. 1.

Fig. 1

Schematic representation of an off-axis DHM setup.



Red Blood Cell Preparation

RBC samples were collected from laboratory personnel of the Laboratoire Suisse d’Analyse du Dopage, Centre Hospitalier Universitaire Vaudois, and for further investigations on the RBC deformation during storage, they were stored at 4°C for a period of time. The total amount of RBCs in the 100 to 150  μl of mainly stomatocytes and discocytes was contained in a high-efficiency particulate air (HEPA) buffer at 0.2% hematocrit, whereas for the echinocyte morphology, the hematocrit concentration was almost 0.15%. To prepare the erythrocytes to be mounted on the DHM stage, they were diluted; 4  μl of suspended erythrocytes were diluted to 150  μl of the HEPA buffer and then carried to the experiment room where the erythrocytes were covered by two cover slides divided by a 1.2-mm-thick splitter. To conduct the RBC experiments, the temperature of the experiment room was 22°C. Before placing the erythrocytes on the DHM stage, the cells were maintained at a temperature of 37°C for 30 min.


Quantitative Phase Image Segmentation of Red Blood Cells

Once we obtain an RBC QPI by the off-axis DHM system, several image-processing algorithms are applied to extract the RBCs from the reconstructed RBC QPI. The first step is to detect the correct RBC samples from the others to be extracted. Next, we remove the noise and background from the RBC QPI by applying the marker-controlled watershed segmentation algorithm to obtain segmented RBC images.24,25 Each RBC sample is detected and extracted for analysis, and more than 14 different characteristics of every RBC sample are automatically measured.11 We further explain the features we used in this paper.

Figure 2 shows the segmentation results for the automatic extraction of RBC samples from the RBC QPI obtained by the DHM system. As we mentioned, the samples we used in this study were extracted from the three main morphologies. The first sample contains the normal RBCs with a biconcave morphology [Fig. 3(a)]. The second sample contains RBCs that are suffering from the sphero-echinocyte process (other morphologies, such as biconcave, can also be found). They were stored for 57 days and then imaged by DHM. The last sample contains cells that are mostly of the stomatocyte morphology. In total, 275 single RBCs were extracted for the feature extraction section.

Fig. 2

Segmentation step: (a) original QPI obtained by off-axis DHM, (b) corresponding segmented RBC image using the marker-controlled watershed algorithm, and (c) binary-segmented RBC image (white bar indicates 10  μm).


Fig. 3

(a) Biconcave sample, (b) RBCs suffering from sphero-echinocytes, and (c) stomatocyte sample.



Feature Extraction and Selection

In this experiment, we extract over 14 2-D and 3-D features related to the RBC profile. We select the best combination of 2-D and 3-D features that can efficiently distinguish between different RBC types and another six 2-D features to compare the clustering method performance.11 The first selected 3-D feature to be used in this experiment is the average cell thickness (ACT), which can be calculated by the following equation:


where h(i,j) is the thickness at the (i,j)’th pixel. For each pixel of (i,j) in the QPI, h(i,j) is calculated by the following equation:


where φ(i,j) is the phase value in radians and the refractive index of RBCs, nrbc, is calculated with dual-wavelength DHM. The refraction index of the HEPA medium, nm, is 1.3334, and k and l are the width and length of the RBC image, respectively. The PSA is another feature used in this experiment. The PSA can be calculated from the following equation:


where N is the number of pixels of the projected cell and p is the size (p=0.159  μm, here) of each pixel in the image. The volume, or MCV, of the RBC can be calculated by the following equation:



The fourth feature is the sphericity coefficient, which can be obtained using the following equation:


where dc and dr are the thickness values in the center and ring of the RBC, respectively. The last feature is the perimeter of the projected RBC profile on the XY plane. The RBC perimeter is the length of the RBC boundary.13 The descriptions of the five selected features and their 2-D and 3-D types are summarized in Table 1.

Table 1

Descriptions and feature type divisions of selected 2-D and 3-D features (best selected feature set).

Segmented single RBC sample
FeatureDescriptionFeature type
ACTThickness of (i,j) pixel for every phase value3-D
PSANumber of pixels within single RBC × one pixel area2-D
Sphericity coefficientCenter part phase value/maximum phase value3-D
PerimeterLength of boundary part of cell2-D
MCVMean corpuscular volume3-D

In our experiment, we evaluated the clustering power of the 2-D features against the best 2-D and 3-D features. We believe that to obtain the best clustering results, we must select the major 2-D and 3-D features. The selected 2-D RBC features and their descriptions are listed in Table 2.

Table 2

Description and feature type divisions of selected 2-D RBC features.

Segmented single RBC sample
PSANumber of pixels within single RBC × one pixel area
Elongation (Ei)Direction of chain code in cell membrane
Perimeter (Pr)Length of RBC cell boundary
Circularity (Ci)Pr2/PSA
Radius (R)Estimation of radius of circle with the area of R=PSA/π


Density-Based Spatial Clustering Applications with Noise

DBSCAN is a clustering method proposed by Ester et al. in 1996, which can identify clusters in large spatial data using the density of data elements. If we consider a point of data as some point distributed in space, the method groups the data points that are in close proximity to each other.27 For a set of samples distributed in feature space to be clustered, DBSCAN has two main parameters, MinPts (number of p points) and ϵ, where p is the number of data points and ϵ (epsilon) is the maximum radius of neighboring p points. It is a nonparametric approach and considers one point as a core point if at least a minimum number of p points is within a distance ϵ of the core point neighborhood. In every step, the core point is changed, and the number of points reachable from the core point with the distance ϵ is checked again; the core point is changed until all reachable points are met, and the unreachable points are marked as noise points.28

Figure 4 shows a schematic representation of the density-based clustering method. The red points are considered as noise points. This implies that the noise data points are far from the other data points. In the DBSCAN clustering method, we need not define the number of clusters, such as in the traditional method. When a cluster is surrounded by another cluster, DBSCAN can cluster the inner and outer clusters effectively. DBSCAN selects two values as input parameters, which can effectively be adjusted according to the density of data points. This parameter can be set by an expert and also be based on data observations.28

Fig. 4

Schematic representation of density-based clustering method. Epsilon and noise points are marked.



k-Medoids Clustering Method

The k-medoids clustering method is a combination of two main algorithms of k-means and medoid shift. Both k-medoids and k-means group similar data and separate them from the other clusters.29 The k-means clustering method is based on the minimization of the total square error, whereas k-medoids attempts to minimize the sum of differences between different samples in the same cluster with the medoid point, which is the center point of every cluster. k-medoids considers one data point as the center of a cluster. It is more powerful than the k-means algorithm because it is not sensitive to outliers. k-means is easily influenced by high-value data points in the dataset because it considers the mean of all data points, whereas k-medoids considers centered data points, which are more reliable.30 When we have an infinite dataset, the k-medoids clustering method considers a data point that exhibits a minimum average dissimilarity to all the other data points. The steps of the k-medoids clustering algorithm are as follows:

  • 1. choose a random k-point among all data points as the medoid point,

  • 2. assign each n data point to the medoid point that is the close distance to medoid point,

  • 3. for all medoid and data points, compute the cost to swap a medoid point between data points and select a medoid point that has the lowest cost, and

  • 4. repeat steps 2 and 3 until all medoid points are fixed.

The cost function is as follows:


where Xi is the data point, Ci is the center point, which is known as the medoid point, and d is the dimension of the data point.


k-Means Clustering Method

The k-means clustering method is a popular method for clustering data observations into k clusters by assigning data observations to each cluster with close distance to the mean value of all similar observations in the same cluster. This process causes the data space to be partitioned according to the mean value of similar data observations.31 If we consider data observations as D={X1,X2,.,Xn} and the number of different clusters as k, the k-means algorithm attempts to find a k centroid point C={C1,..,Cn} to minimize the distance between the data and centroid point.32 The main idea of k-means clustering is to find the best centroid point for each cluster based on the mean value of all data observations in that cluster. During several iterations, k-means attempts to find the best centroid point for each cluster with the closest distance to all observations of k clusters by updating the centroid point. In every iteration, the centroid point will be changed until there are no new points to change. The steps to find the best centroid point in the k-means clustering method are as follows:

  • 1. assign k point to the data space that is going to be clustered,

  • 2. assign each data point of the cluster to the closest centroid point,

  • 3. after assigning all data points to their corresponding centroid points, update the k centroid point, and

  • 4. repeat steps 2 and 3 until there are no points left.

The summation of distances can be calculated by the squared Euclidean distance error function using the following equation:33


where Pi and Qi are the two data points of which their Euclidean distances are going to be calculated.


Experimental Results and Discussion

After feature extraction, since our analysis is in 2-D and 3-D feature space, we applied PCA to reduce the data dimension into 2-D and 3-D feature space and to find a more meaningful basis or coordination of our data instead of original features. Whereas DBSCAN measures the distance of each data sample with neighboring samples using circle radius,28 we used two PCA to be applied to DBSCAN method but for other method, we used three PCA since we analyze results in 3-D feature space as shown in Figs 5 and 10. Figure 5(a) presents the data distribution of the best features according to the first and second PCs. We believe that by varying the number of clusters, we should expect to observe similar morphologies categorized in the same cluster. We evaluated two and three clusters in this study. By utilizing the SI (in k-means clustering), the internal evaluation revealed that increasing the number of clusters decreases the accuracy of the clustering technique (data not shown). Another reason for choosing two and three clusters is that the samples are extracted from three different morphologies. Figure 5(b) shows the original data distribution based on six 2-D features.

Fig. 5

Original data distribution based on the first and second principal components (PCA1 and PCA2) of RBC types of the (a) best features and (b) 2-D features.



Density-Based Spatial Clustering Applications with Noise Clustering Results

As we mentioned before, the density-based clustering approach clusters data based on the density of data observations. Therefore, this method is not based on the shape of data observations, and the number of clusters in a dataset is not predefined. Because of the density of RBC samples distributed in feature space, as shown in Fig. 5, we first applied the DBSCAN clustering algorithm and repeated the clustering experiment while varying the MinPts and epsilon parameters to find the most efficient clustering result. The MinPts parameter denotes the minimum number of data points that can be covered by the radius of a circle, which is known as the epsilon value. Concerning the selection of DBSCAN parameters, there is no special role for MinPts and epsilon parameters. They depend on density and distance of data points around the core point. A low MinPts causes more clusters to build from noise and generates more outliers, and for the epsilon parameter, it is normally considered a number between zero and one on the dataset. Hence, both parameters should not be too small or too large.

DBSCAN considers samples that are far from the other samples as noise or unknown data. This implies that these samples are considered different from the known samples. Therefore, we changed the values of MinPts and epsilon to obtain the best value for these parameters and obtain the best clustering result. In the case of clustering RBCs using DBSCAN based on the 2-D features, according to the data distribution presented in Fig. 5(b), there is no border between different regions of different RBC samples in feature space; therefore, DBSCAN cannot cluster the data points with no border between different regions and considers all data points as one cluster (data not shown). Thus, this proves that the 2-D features are not suitable enough to be clustered by the DBSCAN clustering method. However, for the combination of 2-D and 3-D features, the DBSCAN clustering method can attain a 100% clustering accuracy (see Table 3). The accuracy of the clustering technique can be obtained by comparing an expert’s visual examination (manual clustering) with the automated clustering technique.

Table 3

Clustering performance evaluation results of density-based clustering method using different values for MinPts and epsilon.

EpsilonMinPtsMax iterationNumber of noiseClustering accuracy ratio (%)

Figure 6 shows the graphical representations of the DBSCAN clustering results for a MinPts value of 4 and epsilon values of 0.6 and 0.8. The DBSCAN clustering performance results for MinPts values of 5, 6, and 7, and epsilon values of 0.6 and 0.8 are also graphically represented in Figs. 7Fig. 89, respectively. It is noted that the DBSCAN algorithm that is based on the density of RBC samples in feature space clusters two main RBC types of biconcave and stomatocyte in the same cluster because of their similarities and distinguishes them from the sphero-echinocyte RBC type. In some cases, by changing the epsilon and MinPts value, some of RBCs are considered unknown and marked as noise. It is evident from Fig. 6(b) that we attained a 100% clustering result for the two clusters using a MinPts of 4 and radius of 0.8 and all data points are affected by the radius and MinPts parameters. For the three clusters, according to the density of data points distributed in feature space, the data points of biconcave and sphero-echinocytes are placed in very close proximity to each other. Different radii affect both types and will accordingly be clustered as one cluster. Therefore, we applied the DBSCAN algorithm to cluster RBCs into two clusters.

Fig. 6

Density-based clustering results using MinPts of 4 and epsilon of (a) 0.6 and (b) 0.8.


Fig. 7

Density-based clustering results using MinPts of 5 and epsilon of (a) 0.6 and (b) 0.8.


Fig. 8

Density-based clustering results using MinPts of 6 and epsilon of (a) 0.6 and (b) 0.8.


Fig. 9

Density-based clustering results using MinPts of 7 and epsilon of (a) 0.6 and (b) 0.8.


Different evaluations of the density-based clustering method for different values of epsilon and MinPts, including max iteration, clustering accuracy, and misclassification ratio to cluster all data points in the region of epsilon (radius), as well as number of points far from the other data, which are considered as noise points, are presented in Table 3.


k-Medoids Clustering Results

In this experiment, we applied the k-medoids clustering method. We plotted the original RBC data distribution using three PCAs to analyze the clustering performance in 3-D-space based on the best 2-D and 3-D features (see Fig. 10).

Fig. 10

Original data distribution in 3-D data space based on the three PCAs using (a) best selected feature set and (b) 2-D features.


We applied k-medoids on the 3-D data points based on the three PCAs to cluster RBC data into two and three clusters. For the two clusters, it is observed that the RBC data have been clustered effectively and with high accuracy. As shown in Fig. 11, the k-medoids clustering method found the best central data point in which it can perfectly represent all data points in a specific cluster and marked it as medoid point. The data point that has close distance to the medoid point is allocated to that specific cluster. Because of the similarity between the two main types of biconcave and stomatocyte morphologies, their central point, or so-called medoid point, considers them to be one cluster. Figure 11(a) demonstrates that the sphero-echinocyte RBC samples, which are represented as cluster 1, are fully separated from the biconcave and stomatocyte morphologies using the best 2-D- and 3-D-selected features. Figure 11(b) shows a graphical representation of the k-medoid clustering method using the 2-D features to cluster RBC samples into two clusters. The clustering accuracy for both the best features and 2-D features is presented in Table 4.

Fig. 11

k-medoids clustering results of clustering RBC samples into two clusters using (a) best selected feature set and (b) 2-D features.


Table 4

k-medoids clustering results for two and three clusters on three main types of RBC samples using the best selected features and 2-D features (275 total samples).

Clustering methodNumber of clustersTotal misclassified samplesClustering accuracy rate (%)
k-medoids clustering results using 2-D features22590.9
k-medoids clustering results using best selected features2498.5

Similarly, we also applied the k-medoids clustering method to cluster RBC data into three clusters according to the three main RBC types. A graphical representation of the k-medoids clustering method on the best features and 2-D features is shown in Fig. 12. The medoid points are circled and indicate the center data point for each cluster (see Fig. 12). The experimental results demonstrate that the three main types of RBC samples can effectively be clustered with a high accuracy using the best features. The clustering accuracy of the k-medoids clustering method on 2-D features is presented in Table 4.

Fig. 12

k-medoids clustering results and medoid points for clustering RBC samples into three clusters using (a) best selected feature set and (b) 2-D features.


According to Table 4, the clustering accuracy of the best selected features is much higher than using only the 2-D features. This fact proves that 2-D features cannot suitably discriminate between the different RBC types. As shown in Table 4, since two RBC types of biconcave and stomatocyte have some similarities in shape and features, when we increase the number of clusters, some samples are misclassified and cause a decrease in the clustering accuracy for three clusters.


k-Means Clustering Results

In the first experiment, the k-means clustering method is utilized by defining two numbers of clusters for the RBC data. Figure 13(a) shows the k-means clustering results for clustering RBC types using the best features for the two clusters. As we expected, the clustering results demonstrate that k-means clusters two RBC types of biconcave and stomatocyte into one cluster because of their similarities and separates them from the sphero-echinocyte RBC type. Figure 13(b) shows a graphical representation of the k-means clustering method on 2-D features of the RBC data. The accuracy ratio of both methods is presented in Table 5.

Fig. 13

k-means clustering results on RBC samples for two clusters. Different clusters are represented in different colors using (a) best selected feature set and (b) 2-D features.


Table 5

k-means clustering results on three main types of RBC samples for two and three clusters using the best selected features and 2-D features (275 total samples).

Clustering methodNumber of clustersTotal misclassified samplesClustering accuracy rate (%)
k-means clustering results using 2-D features23786.5
k-means clustering results using best selected features2398.9

Similarly, we also apply the k-means clustering method to the RBC data to cluster data points into three clusters using the best features and 2-D features. Figure 14 indicates that by using the best-selected features, all RBC types of biconcave, discocyte, and stomatocyte can be clustered to a great extent. The clustering accuracy of the 2-D features is significantly lower than that of the best features (see Table 5).

Fig. 14

k-means clustering results on three main RBC samples for three clusters using (a) best selected feature set and (b) 2-D features.


As the k-means clustering method uses the mean of all data points in the same cluster as the centroid point, the experimental results indicate that the k-means method clusters data with a high accuracy level up to 98% for two clusters and 95% for three clusters using the best features. This fact reveals that 3-D features can significantly influence the mean value of data points in the same cluster.

The results obtained by different clustering methods reveal that clustering techniques can be very efficient and accurate if we can choose a good feature set. Specifically, in RBC clustering, the combination of 2-D and 3-D features can significantly increase the accuracy of the clustering results.


Internal Evaluation of Clustering Techniques

In another experiment with internal evaluation of the clustering technique, k-means (almost similar results are obtained for the other clustering techniques) is performed by measuring the SI. SI varies between 1 and +1, and high SI indicates that the input sample is well-matched to its own cluster and poorly matched to neighboring clusters. If most points have a high silhouette value, the plot shows an assessment of how close each sample in one cluster is to samples in the neighboring clusters and using this way, we can measure parameters, such as number of clusters visually.

Silhouette coefficients near +1 mean that the samples are well distinguished from neighboring clusters. Samples that are very close to the neighboring clusters will get zero value, and negative values indicate samples are clustered to the wrong cluster. According to Figs. 15(a) and 15(b), we can see that most of the silhouette values are close to +1. There are a few values below zero that are not well matched to the corresponding cluster.

Fig. 15

SI for (a) two clusters and (b) three clusters.




The quality and functionality of RBCs play major roles in the human health system. Storing blood for long periods can damage the functionality and quality of RBCs. In this study, we applied several unsupervised clustering methods, including DBSCAN clustering, k-medoids, and k-means clustering, for clustering three RBC types of biconcave, stomatocyte, and sphero-echinocyte into two and three clusters. The RBC samples that were visualized by the DHM technique were extracted from the blood sample. DHM provides QPIs of the 3-D profile of RBCs with nanometer accuracy. More than 14 2-D and 3-D features were extracted from every RBC sample. We selected a combinational set of 2-D and 3-D features that can suitably discriminate between the three regular RBC types. The combinational features include the ACT, PSA, sphericity coefficient, perimeter, and MCV. The clustering power of the combinational set of 2-D and 3-D features was compared against a set of six 2-D features, and the clustering results were evaluated for every clustering method. The experimental results and performance of clustering methods indicate that the combinational feature set can yield better RBC clustering. In addition, using the combinational features, we were able to cluster biconcave, stomatocyte, and sphero-echinocyte morphologies to a great extent, which are paramount for RBC abnormality analyses and shape-related diseases.


The authors have no relevant financial interests or conflicts of interest to disclose.


This work was supported by research funds from Chosun University, 2017.


1. M. Bessis, R. Weed and P. Leblond, Red Cell Shape: Physiology, Pathology, Ultrastructure, Springer, New York (2012). Google Scholar

2. P. Canham, “The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell,” J. Theor. Biol. 26, 61–76 (1970).JTBIAP0022-5193 http://dx.doi.org/10.1016/S0022-5193(70)80032-7 Google Scholar

3. C. Uzoigwe, “The human erythrocyte has developed the biconcave disc shape to optimise the flow properties of the blood in the large vessels,” Med. Hypotheses 67, 1159–1163 (2006).MEHYDY0306-9877 http://dx.doi.org/10.1016/j.mehy.2004.11.047 Google Scholar

4. Z. Tu, “Geometry of membranes,” J. Geom. Symmetry Phys. 24, 45–75 (2011).1312-5192 Google Scholar

5. C. Aubron et al., “Age of red blood cells and transfusion in critically ill patients,” Ann Intensive Care 3, 2 (2013). http://dx.doi.org/10.1186/2110-5820-3-2 Google Scholar

6. P. Marik and J. William, “Effect of stored-blood transfusion on oxygen delivery in patients with sepsis,” J. Am. Med. Assoc. 269, 3024–3029 (1993).JAMAAP0098-7484 http://dx.doi.org/10.1001/jama.1993.03500230106037 Google Scholar

7. G. Bosman et al., “Erythrocyte ageing in vivo and in vitro: structural aspects and implications for transfusion,” Transfus. Med. 18, 335–347 (2008).TMEREU0887-7963 http://dx.doi.org/10.1111/tme.2008.18.issue-6 Google Scholar

8. C. Högman and H. Meryman, “Storage parameters affecting red blood cell survival and function after transfusion,” Transfus. Med. Rev. 13, 275–296 (1999).TMEREU0887-7963 http://dx.doi.org/10.1016/S0887-7963(99)80058-3 Google Scholar

9. R. Card, “Red cell membrane changes during storage,” Transfus. Med. Rev. 2, 40–47 (1988).TMEREU0887-7963 http://dx.doi.org/10.1016/S0887-7963(88)70030-9 Google Scholar

10. K. Jaferzadeh and I. Moon, “Quantitative investigation of red blood cell three-dimensional geometric and chemical changes in the storage lesion using digital holographic microscopy,” J. Biomed. Opt. 20 111218 (2015).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.20.11.111218 Google Scholar

11. K. Jaferzadeh and I. Moon, “Human red blood cell recognition enhancement with three-dimensional morphological features obtained by digital holographic imaging,” J. Biomed. Opt. 21, 126015 (2016).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.21.12.126015 Google Scholar

12. L. Theodore, Red Cell Shape, Academic Press, New York (1989). Google Scholar

13. J. Bacus and J. Weens, “An automated method of differential red blood cell classification with application to the diagnosis of anemia,” J. Histochem. Cytochem. 25, 614–632 (1977).JHCYAS0022-1554 http://dx.doi.org/10.1177/25.7.330716 Google Scholar

14. F. Yi, I. Moon and B. Javidi, “Cell morphology-based classification of red blood cells using holographic imaging informatics,” Biomed. Opt. Express 7, 2385–2399 (2016)BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.7.002385 Google Scholar

15. P. Rakshit and K. Bhowmik, “Detection of abnormal findings in human RBC in diagnosing sickle cell anemia using image processing,” Procedia Technol. 10, 28–36 (2013). http://dx.doi.org/10.1016/j.protcy.2013.12.333 Google Scholar

16. M. Buttarello and P. Mario, “Automated blood cell counts,” Am. J. Clin. Pathol. 130, 104–116 (2008).AJCPAI0002-9173 http://dx.doi.org/10.1309/EK3C7CTDKNVPXVTN Google Scholar

17. R. Liu et al., “Recognition and classification of red blood cells using digital holographic microscopy and data clustering with discriminant analysis,” J. Opt. Soc. Am. A 28, 1204–1210 (2011).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.28.001204 Google Scholar

18. J. Dahmen et al., “Automatic classification of red blood cells using Gaussian mixture densities,” Bildverarbeitung für die Medizin 2000, 331–335 (2000). http://dx.doi.org/10.1007/978-3-642-59757-2_62 Google Scholar

19. I. Moon et al., “Automated quantitative analysis of 3D morphology and mean corpuscular hemoglobin in human red blood cells stored in different periods,” Opt. Express 21, 30947–30957 (2013).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.21.030947 Google Scholar

20. A. Vattani, “k-means requires exponentially many iterations even in the plane,” Discrete Comput. Geom. 45, 596–616 (2011). http://dx.doi.org/10.1007/s00454-011-9340-1 Google Scholar

21. M. Ester et al., “A density-based algorithm for discovering clusters in large spatial databases with noise,” Proc. of the Second Int. Conf. on Knowledge Discovery and Data Mining, pp. 226–231 (1996). Google Scholar

22. S. Khanmohammadi, N. Adibeig and S. Shanehbandy, “An improved overlapping k-means clustering method for medical applications,” Expert Syst. Appl. 67, 12–18 (2017).ESAPEH0957-4174 http://dx.doi.org/10.1016/j.eswa.2016.09.025 Google Scholar

23. P. Arora and S. Varshne, “Analysis of k-means and k-medoids algorithm for big bata,” Procedia Comput. Sci. 78, 507–512 (2016). http://dx.doi.org/10.1016/j.procs.2016.02.095 Google Scholar

24. F. Yi et al., “Automated segmentation of multiple red blood cells with digital holographic microscopy,” J. Biomed. Opt. 18, 026006 (2013).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.18.2.026006 Google Scholar

25. X. Yang et al., “Automated segmentation and tracking of cells in time-lapse microscopy using watershed and mean shift,” in Proc. of the Int. Symp. on Intelligent Signal Processing Communication Systems, pp. 533–536 (2005). http://dx.doi.org/10.1109/ISPACS.2005.1595464 Google Scholar

26. M. Kim, Digital Holographic Microscopy, Springer, New York (2011). Google Scholar

27. J. Chen and H. He, “A fast density-based data stream clustering algorithm with cluster centers self-determined for mixed data,” Inf. Sci. 345, 271–293 (2016). http://dx.doi.org/10.1016/j.ins.2016.01.071 Google Scholar

28. R. Albarakati, Density Based Data Clustering, California State University, San Bernardino (2015). Google Scholar

29. A. Reynolds, G. Richards and V. Rayward, “The application of k-medoids and pam to the clustering of rules,” in Int. Conf. on Intelligent Data Engineering and Automated Learning, Springer, New York (2004). Google Scholar

30. R. Krishnapuram, A. Joshi and L. Yi, “A fuzzy relative of the k-medoids algorithm with application to web document and snippet clustering,” in IEEE Int. Fuzzy Systems Conf. Proc., Vol. 3 (1999). http://dx.doi.org/10.1109/FUZZY.1999.790086 Google Scholar

31. X. Wu et al., “A hybrid fuzzy K-harmonic means clustering algorithm,” Appl. Math. Modell. 39, 3398–3409 (2015).AMMODL0307-904X http://dx.doi.org/10.1016/j.apm.2014.11.041 Google Scholar

32. M. Capó, A. Pérez and J. Lozanom, “An efficient approximation to the K-means clustering for massive data,” Knowl.-Based Syst. 117, 56–69 (2017). http://dx.doi.org/10.1016/j.knosys.2016.06.031 Google Scholar

33. Z. Huang, “Extensions to the k-means algorithm for clustering large data sets with categorical values,” Data Min. Knowl. Discovery 2, 283–304 (1998). http://dx.doi.org/10.1023/A:1009769707641 Google Scholar


Ezat Ahmadzadeh is a PhD student in the Computer Engineering Department at Chosun University. He received his BS degree in software engineering from Islamic Azad University of Mohabad in 2011 and his MS degree from Islamic Azad University at the Science and Research Branch of Tehran (West Azarbayjan) in 2014. His current research interests include image processing, digital holography, machine vision, and parallel programming.

Keyvan Jaferzadeh is a PhD student in the Computer Engineering Department at Chosun University. He received his BS degree in software engineering and his MS degree in mechatronic engineering in 2006 and 2010, respectively. His current research interests include image processing, digital holography, image compression, and machine vision.

Jieun Lee received her BS degree in computer engineering from Ewha University in 1997, her MS degree from POSTECH in 1999, and her PhD from Seoul National University in 2007. From 1999 to 2002, she worked at LG Electronics Institute of Technology as a research engineer. She is currently an associate professor of the Department of Computer Engineering, Chosun University, Republic of Korea. Her research interests are in geometry processing and computer graphics.

Inkyu Moon received his BS degree in electronics engineering from SungKyunKwan University in Korea in 1996 and his PhD in electrical and computer engineering from the University of Connecticut, USA, in 2007. He joined Chosun University in Korea in 2009 and is currently a professor at the School of Computer Engineering there. His research interests include digital holography, biomedical imaging, and optical information processing. He is a member of IEEE, OSA, and SPIE.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Ezat Ahmadzadeh, Ezat Ahmadzadeh, Keyvan Jaferzadeh, Keyvan Jaferzadeh, Jieun Lee, Jieun Lee, Inkyu Moon, Inkyu Moon, } "Automated three-dimensional morphology-based clustering of human erythrocytes with regular shapes: stomatocytes, discocytes, and echinocytes," Journal of Biomedical Optics 22(7), 076015 (25 July 2017). https://doi.org/10.1117/1.JBO.22.7.076015 . Submission: Received: 29 March 2017; Accepted: 5 July 2017
Received: 29 March 2017; Accepted: 5 July 2017; Published: 25 July 2017

Back to Top