Human red blood cell recognition enhancement with three-dimensional morphological features obtained by digital holographic imaging

Abstract. The classification of erythrocytes plays an important role in the field of hematological diagnosis, specifically blood disorders. Since the biconcave shape of red blood cell (RBC) is altered during the different stages of hematological disorders, we believe that the three-dimensional (3-D) morphological features of erythrocyte provide better classification results than conventional two-dimensional (2-D) features. Therefore, we introduce a set of 3-D features related to the morphological and chemical properties of RBC profile and try to evaluate the discrimination power of these features against 2-D features with a neural network classifier. The 3-D features include erythrocyte surface area, volume, average cell thickness, sphericity index, sphericity coefficient and functionality factor, MCH and MCHSD, and two newly introduced features extracted from the ring section of RBC at the single-cell level. In contrast, the 2-D features are RBC projected surface area, perimeter, radius, elongation, and projected surface area to perimeter ratio. All features are obtained from images visualized by off-axis digital holographic microscopy with a numerical reconstruction algorithm, and four categories of biconcave (doughnut shape), flat-disc, stomatocyte, and echinospherocyte RBCs are interested. Our experimental results demonstrate that the 3-D features can be more useful in RBC classification than the 2-D features. Finally, we choose the best feature set of the 2-D and 3-D features by sequential forward feature selection technique, which yields better discrimination results. We believe that the final feature set evaluated with a neural network classification strategy can improve the RBC classification accuracy.


Introduction
Human blood contains different type of cells; however, red blood cells (RBC) or erythrocyte are the most abundant cell type.RBC contains hemoglobin, which binds to either oxygen or carbon dioxide.This allows oxygen to be transported to tissues and organs and carbon dioxide to be taken away during microcirculation.The biconcave shape of the erythrocyte is extremely important in the functionality of RBCs.2][3] According to its importance, it has drawn considerable attentions into the pathology research in the clinical relevant blood diseases.Pathological disorders can modify RBCs and lead to significant alterations in its original shape. 45][6][7][8][9] Also, counting cell types in the blood sample is another important task for investigating clinical status, which can be evaluated by wellknown methods such as complete blood count or RBC distribution width, which are part of cytometry field.Because an automated cell counter samples and counts so many cells, the results are reliable in most of the cases. 10However, certain abnormal cells in the blood may not be identified correctly, requiring manual review and identification of any abnormal RBCs the instrument could not categorize.This information can be very helpful regarding identifying the cause of a patient's anemia.Abnormal increase or decrease in RBC counts as revealed in a complete RBC count may indicate that you have an underlying medical condition that calls for further evaluation.
In the case of RBC, biconcaves are a substantial type in a healthy person, but there are other RBCs types with the different percentage varying between healthy and unhealthy persons.It has been shown that the percentages of different types of RBCs will be distinct according to the RBC diseases type. 10ccordingly, it is essential to measure the percentage of each RBC type in a blood sample consists of multiple RBCs for diagnosis and drug testing subjects.Typically, the diagnosing is performed by a human expert, and it shows some drawbacks, such as time-cost consuming and inaccuracy.Generally, experts visualize the sample in the images through a microscope based on their subjective knowledge from the viewpoint of intensity, morphology, texture, and so on based features.Usually, smallscale differences in the features are overlooked by human eyes especially for the border-line diagnostic scenario.
However, the situation has changed completely by the emergence of automatic classification algorithms.These techniques have been applied to problems in biology and have shown promising results for automatic recognition and classification of various micro-organisms. 11,124][15] The fundamental benefit of artificial neural network (ANN) is nothing, but it does not use any mathematical model since ANN learns from data sets and identifies patterns in a sequence of input and output data without any previous assumptions about their type and interrelations.Also, ANN eliminates the drawbacks of the conventional methods by extracting the wanted information using the input data.
7][18][19] These methods generally have good performance but most of them have a significant number of features since they need to discriminate groups by utilizing 2-D features.However, in the case of RBC, which is transparent or semitransparent, we cannot take advantage of conventional intensity-based microscopes.Therefore, we believe that to obtain a satisfied level of accuracy, the classification and recognition should also take into account the threedimensional (3-D) shape of RBCs.6][27] In studying RBC, DHM enables the measurement of 3-D features, such as mean corpuscular volume, surface area, SAV ratio, functionality factor, sphericity index, and sphericity coefficients. 28,29hemical parameters of MCH and MCHSD can also be obtained due to DHM.Accordingly, we believe that any automated RBC classification that can distinguish different RBC types accurately should take into accounts the benefits of DHM imaging technique.
In this study, four main types of RBC shapes, biconcave (doughnut shape), flat discs, stomatocyte, and echinospherocyte are interested RBCs for the quantitative determination of the percentage of RBC types in multiple human RBCs.The reason for differentiating doughnut-shaped and flat-disc RBCs is that it can help in the applications of separating old cells versus young cells.0][31][32][33][34] Transfusion of these old samples might have critical consequences according to the previous studies. 35,36t first, RBCs are visualized by off-axis DHM and the quantitative phase images are reconstructed by the numerical algorithm. 37,38Then, single RBCs are extracted from images with multiple RBCs using the watershed algorithm. 39At the next step, following 2-D features of projected surface area (PSA), perimeter, radius, elongation, and PSA to perimeter rate are extracted.In this paper, we have ignored extracting 2-D features related to the inner section of RBC, unlike the previous method proposed by Refs.25, 26, since flat-disc and echinospherocyte RBCs do not have the inner section.Also, volume, surface area, SAV ratio, average RBC thickness, sphericity index, sphericity coefficient and functionality factors, and MCH and MCH surface density (chemical properties of RBC) are extracted from single RBCs.The latter feature-set is related to the morphological and biochemical properties of RBC 3-D profile.Along with the 3-D features, two new features related to the ring section of RBC are introduced.These features add significant information to the classification model and increase the discrimination power of the classifier.Then, each feature set is fed into PRNN, separately, and the classification results are compared using 10-fold cross validation (CV) technique.Since we are involved in a classification model with nonlinear decision boundary, we have decided to use PRNN strategy.In PRNN, the training algorithm is Bayesian regulation back-propagation, which updates the weights according to Levenberg-Marquardt optimization technique and the activation function for midlevel layers is hyperbolic tangent sigmoid.
Finally, to propose the best feature set, concerning both 2-D and 3-D features, sequential forward feature selection (SFFS) is utilized here.It has been shown that the best performance of a classification model can be achieved by selecting the most informative features and remove noisy ones that are either redundant or irrelevant. 40Indeed, reducing the number of features can shorten training time, reduce the complexity of classifier, and simplify the model for interpretation goals.SFFS technique tries to select a subset of variables that best predict the data by sequentially selecting features until there is no improvement in the prediction.
In this paper, we have extracted 108 biconcave RBCs from a healthy sample stored for 1 day in the blood bank, 106 samples of stomatocyte shape from a sample with predominantly of stomato cells, 38 samples of flat-disc shape, and 71 samples of echinospherocyte shape for training and testing PRNN.Flat-disc and echinospherocyte cells are extracted from RBC samples stored in the blood bank for 40 days and 57 days, respectively.Performance comparison is evaluated by calculating misclassification rate of 10-fold CV technique.It is often claimed that leave-one-out-cross-validation (LOOCV) has higher variance than k-fold CV, and that it is because the training sets in LOOCV have more overlap.This makes the estimates from different folds more dependent than in the k-fold CV and, hence, increases the overall variance. 41ur experimental results demonstrate that the PRNN trained by 3-D features gives a good performance in classifying and counting RBCs in multiple human RBCs in an automated manner in comparisons with the 2-D features.In addition, we introduce the best set of features that combines 2-D and 3-D features to improve the RBCs classification accuracy.We believe that the final feature set evaluated with the presented neural network classification strategy can provide better discrimination results.
This paper is organized as follows.Section 2 explains the general scheme of the off-axis DHM to image RBCs and RBCs preparation.Section 3 explains 2-D and 3-D features extracted in this research and gives a short review about them.In Sec. 4, we will focus on designing the PRNN and counting different RBCs in multiple RBCs, automatically.The experimental results and discussions are provided in Sec. 5. Finally, the conclusion is presented in Sec. 6.
2 Off-Axis Digital Holographic Microscopy and RBC Preparation

Off-Axis Digital Holographic Microscopy
Figure 1 shows the off-axis DHM based on the Mach-Zehnder interferometer geometry. 42In this optical setup, light from a coherent source (HeNe laser diode source λ ¼ 682 nm) is split into reference and object beams by a beam splitter.The object beam is transmitted through the RBC sample and magnified by a 40 × ∕0.75 NA microscope objective and interferes with the reference beam.The interference patterns between the diffracted object beam and the reference beam are recorded onto a CCD camera.The reconstruction of the RBC wavefront is obtained from the recorded hologram by using the numerical algorithm described in Refs.37 and 38.
After reconstructing RBC phase images, single RBCs are extracted by the maker-controlled watershed segmentation algorithm. 39Figure 2 shows a reconstructed phase image of a sample consisting of biconcave, flat-disc, and stomatocyte shape RBCs.This image shows that in a single sample it is possible to see RBC with different morphologies.

RBC Preparation
The RBCs of healthy laboratory personnel were obtained through the Laboratoire Suisse d' Analyse Du Dopage, CHUV and stored at 4°C during the storage period.The DHM measurements were performed several days after the blood was collected from the laboratory personnel.A total of 100 to 150 μl of RBC stock solution was suspended in a high-efficiency particulate air (HEPA) buffer at 0.2% hematocrit for predominantly stomatocyte and discocyte-shaped RBCs while at a concentration of ∼0.15% for predominantly echinocyte-shaped RBCs.A total of 4 μl of the erythrocyte suspension was diluted to 150 μl of the HEPA buffer and introduced into the experimental chamber, including two cover slips separated by spacers 1.2 mm thick.The cells were incubated for 30 min at a temperature of 37°C before mounting on the chamber on the DHM stage.All experiments were performed at room temperature (22°C).
3 Feature Extraction

Two-Dimensional Features
After segmentation step and extracting many single RBCs, features can be extracted.We first start with 2-D features.Following features are extracted in 2-D case.
Elongation of the RBC is a measure of the ratio of width to length for oblong RBCs.It can be computed from the chain code by summing the number of each type of elements 0 to 7 and combining 0 and 4, 1 and 5, 2 and 6, and 3 and 7. 16 Average and STD values of the above features are in the agreement with previously reported values (data not shown here). 25,26

Three-Dimensional Features
Since some 3-D features require thickness of the RBC, we first need to convert phase image into thickness image.Accordingly, the thickness value hði; jÞ for each pixel of ði; jÞ with phase value φði; jÞ in a phase image can be expressed as 28,29 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 3 2 6 ; 4 6 8 hði; jÞ where φði; jÞ is the phase value in radians, and the refractive index of RBCs, n rbc , has been measured with a dual-wavelength DHM.Here, n rbc is 1.396 with no significant difference among groups of RBCs.The index of refraction of the HEPA medium, n m , is 1.3334.Surface area is another important property and has the main contribution in different 3-D features.Generally speaking, the surface area of the RBC is the surface area of the membrane mesh plus PSA.The method in this paper splits and divides RBC surfaces into smaller regular areas (triangles here) and adds these smaller areas to give the entire surface area.Understandably, the accuracy of such a calculation is dependent on the smaller area chosen. 29g. 1 Schematic representation of an off-axis DHM setup.3-D features listed in Table 1 are extracted in this paper.We only give a short description of the eight features related to the morphological properties of RBC, but interested readers can refer to Refs.28 and 29 for the detail of each feature.Regarding the calculation of three features (F9-F11), we obtained many points over the ring section of RBC by applying two techniques.First, we estimate the ring section (blue triangles in Fig. 3) by calculating the radius of a circle having the area of the projection of RBC on X-Y plane (the ring is around three-fourth of the RBC radius).Then, we update the position of each point on estimated ring (blue triangles) by finding the thickest point in a 3 × 3 neighbor (red stars in Fig. 3).From now on we name the red stars as RPs.The single green point shows the center of RBC regarding the calculation of sphericity coefficient.In this study, we have decided to obtain 30 RP points.The reason is that the total number of points on the ring section for the stomatocyte RBCs is between 90 and 120 points according to our calculation.Also, since we are looking at a 3 × 3 neighborhood around the blue triangles to find the point on the ring (red stars), there is high possibility of overlapping between red stars.For example, Fig. 4 shows the position of red stars and blue triangles when we increase the number of points.We can see that by increasing the number of red stars many of them overlap.
We believe that earlier 3-D features can distinguish among different RBCs since they are related to the 3-D profile of RBC.The statistic t-test using two-sample Kolmogorov-Smirnov test revealed that some of these features are independent (data not shown).Table 2 shows average and STD values of each feature for each RBC type.

Pattern Recognition Neural Network
ANNs are highly simplified mathematical models of biological neural networks having the ability to learn and provide meaningful solutions to the problems with high-level complexity and nonlinearity.The ANN approach is faster compared to its conventional techniques, robust in noisy environments, and can solve a wide range of problems.4][15] A typical ANN is presented in Fig. 5.An important application of neural networks is pattern recognition that can be implemented by using a feed-forward neural network with a specific training function and specific function in the output layer.During training, the network is trained to associate outputs with input patterns.When the network is used, it identifies the input pattern and tries to output the associated output pattern.
The information processing in the ANN strategy starts from the input layer to the hidden layer (from one hidden layer to another one if there are more than one) and from the last hidden layer to the output layer.A synaptic weight is assigned to each link to represent the relative connection strength of two nodes at both ends in predicting the input-output relationship.y j (the output) of any node j, is given as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 3 2 6 ; 3 3 3 where W i is the connection weight, X i denotes the i'th input of node j, n is the number inputs to node j, and b j is the bias value.Function f or so-called activation function determines the response of a node to the total input signal that is received.Generally, the activation function for hidden layer in PRNN is hyperbolic tangent sigmoid transfer function due to the following advantages: (1) it limits the range of data to values between −1 and 1.This function is almost linear near the mean.It is smooth and has monotonic nonlinearity property at both extremes, (2) it remains differentiable everywhere and the sign of the derivative is unaffected by the normalization.Generally, the differentiable requirement is needed for hidden layers and hyperbolic tangent sigmoid transfer function is often recommended as being more balanced, (3) the 0 for hyperbolic tangent sigmoid transfer function is at the fastest point (highest gradient or gain) and not a trap, (4) the hyperbolic tangent sigmoid transfer function suit the output layer to the competitive outputs of softmax function, (5) since it can be estimated by ð2∕1 þ e −2x Þ − 1 it can be implemented faster in MATLAB ® .For output layer in PRNN, the activation function is softmax transfer function (normalized exponential) that can be interpretable as posterior probabilities for a categorical target variable.It is highly desirable for those outputs to lie between zero and one and to sum to one.The purpose of the softmax activation function is to enforce these constraints on the outputs.Let the input to each output unit be q l , l ¼ 1; : : : ; k, where k is the number of classes.Then the softmax output y l is E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 6 3 ; 1 0 2 y l ¼ e q l P k m¼1 e q m : (3) According to Eq. (3), sum of all "y"s is equal to unity and can be interpreted as the posterior probability for the final decision.
The training algorithm updates the weight and bias values according to Levenberg-Marquardt optimization.It minimizes a combination of squared errors and weights and then determines the correct combination so as to produce a network that generalizes well (Bayesian regularization). 43

Experimental Results and Discussion
In this experiment, 108 RBCs are labeled as biconcaves, 106 RBCs labeled as stomatocytes, 38 RBCs are labeled as

Comparison Between 2-D and 3-D Features
In the case of 2-D features (Table 3), 10-fold CV revealed that the total misclassification rate is significantly high.In detail, misclassification of each group is: flat-disc 64%, stomatocyte 13.4%, biconcave 32.3%, and echinospherocyte are 4.2%.
Only echinospherocyte RBCs can be accurately classified by taking into account 2-D features while other categories have significant error.According to the confusion matrix, PRNN confuses between biconcaves and flat-discs by using 2-D  features (data not shown).In contrast, 3-D features explained in Table 1 generate more accurate and interesting results.
According to the results of 10-fold CV, misclassification rates are 0%, 1.6%, 3.2%, and 0% for flat-disc, stomatocyte, biconcave, and echinospherocyte RBCs, respectively.Table 4 shows   In another experiment, we evaluated the normalized mutual information between each feature in 2-D and 3-D features. 40In the case of 2-D features, it turns out, as shown in Table 5, that there is big mutual information among some features.Mutual information is the amount of information that two features share.If the mutual information between the two features is large (small), the two features are closely (not closely) related.If the mutual information becomes zero, the two features are independent.For example, 2-D-F1 has significant mutual information with features 2-D-F2, 2-D-F4, 2-D-F5, and 2-D-F6 (see first row of Table 5).Therefore, this feature is statistically redundant and cannot add significant information.

Combining 2-D and 3-D Features and Select the Best Feature-Set
We believe that the best classification model should take into account not only 3-D features but also 2-D features.However, it is worth mentioning that not any feature can add significant information to the classification model.Therefore, we believe that the performance of any classification model can be enhanced by taking advantages of feature selection (FS) strategy.In FS, we seek to find the best set of the features that has the strongest ability (or strong) to distinguish each class.Generally speaking, FS preserves the original features intact; features deemed unimportant/irrelevant/redundant are simply eliminated from further consideration while selecting only those features that significantly contribute to the classification problem.Therefore, FS can reduce the number of features (variables) of the classification problem and make the model simpler (or less complex) and shorten training time. 40We also implemented FS in this research by using SFFS technique.Generally speaking, in SFFS features are sequentially added to an empty candidate set until the addition of further features does not decrease the criterion.It has two components of an objective function, called the criterion, and a sequential search algorithm.In former, common criteria are misclassification rate for classification objects (similar in this paper) and mean squared error for regression models.A sequential forward search algorithm adds features from a candidate subset while evaluating the criterion.Since an exhaustive comparison of the criterion value at all 2 n subsets of an n-feature data set is typically infeasible (sometimes feasible but time-consuming), sequential searches grow the candidate set.It turns out that the following features of average RBC thickness (3-D-F1), top view surface (TVS) area volume ratio (3-D-F4), sphericity coefficient (3-D-F9), the upper side of the ring divided by lower side of the ring (3-D-F11), and perimeter (2-D-F2) can better classify multiple RBCs in this research.Our experimental results also show that adding more features does not add significant discrimination ability to the final classification model.As it has been mentioned earlier, SFFS is responsible for adding or removing features from the final feature set.One component of SFFS is the objective function, which here is the misclassification rate.SFFS starts from an empty set and add features one by one to the set and evaluates the misclassification rate.If there is a significant change in the objective function (misclassification rate) then the feature can be added to the final feature set (see Table 6).
According to Table 6 after adding the 7th and 8th features, the misclassification rate never changes.It means that they have no contribution to the final feature set.Adding the 6th feature decreases the misclassification rate marginally but we still did not consider it as part of the feature set since we wanted to keep the final feature set as small as possible (only five features).
Figure 7 shows the data distribution RBCs for each selected feature.We can see that the distribution of average RBC thickness of echinospherocyte RBCs has almost no overlap with other distributions [Fig.7(a)].Table 7 shows misclassification rate of the final feature-set and PRNN approach.
The confusion matrix of the test set shows that PRNN sometimes confuses stomatocytes and echinospherocytes because there are cases in which RBC has a morphology similar to both stomatocyte and echinospherocytes (see Fig. 8).According to Fig. 8, we can see that RBC is similar to both stomatocyte and echinospherocytes morphology and posterior probability for belonging to stomatocyte and echinospherocyte categories are 0.33 and 0.66, respectively.Even for a human examiner, it can be difficult to put it in the correct category.
In another experiment, we tried to count different types of RBCs in five quantitative phase images with multiple RBCs.Images are inspected visually by the human inspector and then   Then, the percentage of the different types of RBCs in the RBC phase images is calculated (see Fig. 9).As expected, the classifier showed that at the first sample stomatocyte RBCs are dominant.In contrast, in the second and fifth figures biconcave and flat-disc RBCs are the major types.Third and fourth figures illustrate that echinospherocytes are dominant RBCs.Although there is a small error in counting nondominant RBCs, the main and important thing is counting and reporting the dominant type for further investigation.
According to the results, the proposed feature set and classifier can automatically count and categorize different types of RBCs in human RBCs by taking advantages of 2-D and 3-D profiles of RBC.The classifier is helpful to assess RBC-related abnormality since the ratio of the different types of RBCs is associated with certain types of diseases.There are some disadvantages in using NN technology regarding the classification such as there is no formula for the number of nodes and hidden layers.Generally, these numbers in the network  are problem dependent and are decided by the trial-and-error method.
We believe that DHM, by providing the quantitative phase images, allows for doing different analysis especially classification goals.For example, someone can easily classify echinoechino-spherocyte from the rest of cells easily by using a conventional binary support vector machine (SVM) classifier and one or two features.Figure 10(a), e.g., shows the density of average cell thickness of echinospherocytes against other RBCs.Any binary classifier can be used to separate these two groups.Figure 10(b) shows the scattering of two groups based on average thickness (AT) value and surface area, and an SVM classifier with its boundary region.

Conclusions
Automatic classification of different types of RBCs in human RBC is a challenging and important problem for pathological diagnosis.In particular, as far as RBCs are concerned, classification based on their 3-D profile is highly relevant.Another issue is that since demands on laboratories are ever increasing  and funding and staffing levels are generally below the desired level, the implementation of a system that shortens staff time, cost-effective, and noninvasive is greatly desired.In this paper, we have presented and assessed the use of PRNN applied to the 2-D and 3-D features of RBCs obtained through DHM in order to categorize and count biconcave, stomatocyte, flat-disc, and echinostomatocyte RBCs in an RBC sample with multiple types.Six 2-D features and 13 3-D features have been extracted and classification results are compared right after.Our experimental results show that the 3-D features have more useful information in RBC classification.In addition, FS shows that average RBC thickness, TVS area volume ratio, sphericity coefficient, the upper side of the ring divided by lower side of the ring, and RBC perimeter can better classify RBCs into the desired categories.The experimental results and the performance imply that the final feature set can help in classifying and counting RBCs, which is substantially important in analyzing RBC abnormality and shape-related diseases.

Other cells Echinospherocyte
Fig. 10 (a) AT value distribution for echinospherocytes and other RBCs and (b) an SVM classifier to separate echinospherocytes and other cells using AT and total surface area.

Fig. 2
Fig. 2 Reconstructed image and three different RBCs inside.A is a flat disc RBC, B is a biconcave RBC, and C is a stomatocyte.White bar is 10 μm.

[(
hði; jÞ is the thickness at ði; jÞ'th pixel] φ is the average phase value, λ is the wavelength of the light source of the setup, and α ¼ 0.00196 dl∕g is a constant known as the specific refraction increment related to the protein concentration) 3-D-F13 MCH surface density (MCHSD) MCHSD ¼ MCH PSA

Fig. 3 3
Fig. 3 3-D representation of four RBC categories and points on the ring section.(a) A typical biconcave sample, (b) a flat disc with center raised, (c) a stomatocyte RBC, and (d) a spherocyte RBC.

Fig. 4 3
Fig. 4 3-D representation of RBC ring section with red stars and blue triangles on it.(a) Number of points is 60 and (b) number of points is 600.

Fig. 5
Fig.5ANN configuration with five nodes, two nodes, and two hidden layers.
the classification error for 2-D and 3-D features.It is noted that the presented neural network classification strategy was used to evaluate the discrimination power of the feature set based on 3-D morphological properties of RBC against 2-D features.The classification results obtained with the neural network demonstrate that the 3-D features can be more effective in RBC classification than the 2-D features.

Fig. 6
Fig. 6 Samples of each RBC group used in this research: (a) four samples of flat disks, (b) four samples of stomatocyte morphology, (c) four samples of biconcave RBC, and (d) four samples of spheroechinocyte RBC.

Fig. 7
Fig. 7 Data distribution of the best feature set.(a) Average thickness value.(b) Sphericity coefficient.(c) Upper side of the ring/lower side of the ring.(d) Top view surface area volume ratio.(e) Perimeter.

Fig. 8
Fig. 8 An RBC sample that confuses neural network, resembles both stomato and spherocyte.(a) 3-D representation and (b) representation on X -Y plane.

Table 4
Misclassification results of 2-D and 3-D features.

Table 5
Normalized mutual information between 2-D features.

Table 6
Misclassification rate after adding each feature to the feature set.
]. Figure 9(e) shows an image with 40 days storage time with many flat-disc RBCs (flat-disc: 16/36, stomatocytes: 7/36, biconcaves: 5/36, and echinospherocytes: 8/36).The numbers in parenthesis show the number of each morphology obtained by the human inspector.At first, each image is segmented into many RBCs since feature extraction should be applied at the single cell level.

Table 7
Misclassification results for the best feature set obtained by sequential features selection technique.