GB2093586A

GB2093586A - Apparatus for classification of cells

Info

Publication number: GB2093586A
Application number: GB8204902A
Authority: GB
Original assignee: Rush Presbyterian St Lukes Medical Center
Current assignee: Rush Presbyterian St Lukes Medical Center
Priority date: 1978-02-03
Filing date: 1979-02-02
Publication date: 1982-09-02
Also published as: CA1135979A; FR2423205B1; GB2093586B; JPS54119294A; JPH0418267B2; DE2903625A1; GB2013878A; DE2903625C2; FR2423205A1; GB2013878B

Abstract

The invention provides an apparatus for automatically analyzing red blood cells in a sample of a patient's blood, comprising; means 14 for examining the red blood cells in a patient's blood sample; means 16 for measuring characteristics of the red blood cells; means 22 for classifying the cells into a plurality of mutually exclusive subpopulations including a normal and abnormal subpopulation; means 28 for determining subpopulations parameters including the dispersion of the distribution of at least one of said subpopulations parameters; and means for reporting a description of the cells based on said dispersion. The apparatus described is substantially similar to that described in specification 2103878. <IMAGE>

Description

1

SPECIFICATION

Apparatus for classification of cells 1 10 GB 2 093 586 A 1 This invention relates to an apparatus for the classification of cells.

The invention relates to an apparatus for automatically analyzing blood, and through the accumulation of measured properties from individual cells, of thereby classifying the blood specimen according to its close resemblance to the normal or to various pathological conditions. More particularly it is concerned with automatically classifying red blood cells and by accumulating measurements relative to each cell generating characteristic values, which by automated means identify the given blood specimen as typical or either normal or of the pathological condition of a specific type of red cell disorder or anemic condition.

According to present medical practice, the diagnosis of a particular one of more than a dozen major types of anemic used three broad categories of information: 1) mean descriptors of cell number, size and hemoglobin content, 2) subjective microscopic visual evaluation of the stained blood cells by a trained hernotologist, and 3) specific biochemical or other tests to pinpoint the precise cause of the anemia.

With conventional equipment, the most common of the first category or red cell sample descriptors are: 1) the red cell count, or the number of red cells per unit volume of blood, 2) the hemoglobin content, or amount of hemoglobin per unit volume of blood, 3) the packed cell volume, or the percentage of blood occupied by red cells, 4) a mean cell size parameter, usually taken as the mean cell volume, which is derived by dividing the packed cell volume by the red cell count, 5) a mean cell hemoglobin parameter, which is derived by 20 dividing the total blood hemoglobin content by the red cell count, and 6) the mean cell hemoglobin concentration, which is derived by dividing the total cell hemoglobin by the packed cell volume.

In the second category, a subjective visual evaluation of the stained blood speciment relates to the tedious and time consuming process by a hernotologist of examining a blood film under the microscope and identifying characteristic abnormal cells, such as large cells, or macrocytes, small cells such as microcytes, 25 target cells, elongated cells such as sickle cells, and giving estimates of size variation e.g. anisocytosis +1, +2, or +3, or a subjective evaluation of population shape changes, such as poikilocytosis +1, +2, or +3.

Today, in addition to the overall mean red cell parameters, and the visual descriptions, certain biochemical and other sophisticated tests are often performed to further clarify the pathology of the anemia. These include: iron kinetics tests, serum iron level tests, hemoglobin electrophoresis, folic acid level tests, vitamin 30 B1 2 tests, and the extraction of a bone marrow sample for evaluation of maturation changes and stainable iron levels. These are very time consuming and expensive, and as to the extraction of a bone marrow sample, very painful.

In accordance with the above, one aspect of the present invention replaces the tedious visual examination with an automatic classification of the normal and abnormal red blood cells into sub-populations and 35 extracts meaningful red blood cells parameters for the separate sub- populations. These and other parameters are used for the automatic classification of the blood specimen with respect to the categories of anemia. One difficulty encountered in separating the normal and abnormal cells into meaningful and widely recognized sub-populations on an automated basis is that of accurately segregating the cells by their morphology and color, particularly where their respective areas or sizes and shapes overlap and their 40 respective principal distinguishing feature is the configuration of their respective central palors (or a lack of central pallor). Central pallor is the thin, disc-shaped central area of red blood cells which may be circular and particularly pronounced for some cells. For instance, target cells and normocytes may have substantially the same size area and shape, but differ in a central pallor configuration. Thus, to distinguish between these cells, and subsequently to distinguish between anemias, the automated analysis should be able to examine 45 and classify cells on the basis of their interior configuration, as well as their exterior configurations.

Other cells, such as spiculed red blood cells, may have the same general size, area and interior configurations of normocytes or the like, but are distinguished principally by their indented spiculed perimeters. Likewise, adding to the difficulty of classifying abnormal cells such as sickle cells from other elongated cells is thatthey may have similar peripheral measurements, sizes, and areas, but differ principally from one another in the presence of pointed projections, or spicules. Still other abnormal cells may be separately categorized from other morphologically similar cells only by their hemoglobin content, measured in terms of color or density. Therefore, it may be desirable to distinguish the hypochromic cells from those that are normochromic.

It has been further discovered that the classification of the red cells into subpopulations is not always discerning enough with regard to automatic anemia classification, but thatfurther statistical descriptions of these subpopulations are often required. For example, different anemia bloods may have the same percentage of round biconcave cells but show considerable variation in the dispersion of the cells with regard to size, hemoglobin content and central pallor. Other anemias may result in blood with the same percentage of elonated or spiculed cells, but differ significantly with regard to measure of mean size, hemoglobin, or with regard to the total population skewness relative to a shape measure.

As will be explained in greater detail hereinafter, the present invention is described in connection with a microscopic slide, digital image and pattern recognition system. However, the invention is not to be construed as limited to such a system, as the feature analysis of the sample to classify the blood and to test the blood for abnormalities maybe performed using other techniques, such as a coherent optical analysis 2 GB 2 093 586 A 2 technique disclosed in U.S. Patent No. 3,947,123; or a liquid flow process technique such as disclosed in U.S. Patent Nos. 3,819,270 and 3, 822,095. To be commercially feasible, the digital image and pattern recognition process for the blood cells should operate on a real time basis and with suff icient speed and accuracy that it will perform as well as the now commercially accepted leukocyte differential counting system, such as for example, the LARC manufactured Corning Glass Works of Corning, New York, and generally diclosed in U.W. Patent No. 3,883,852.

The Coulter Counter, manufactured by Coulter Electronics, Hialeah, Florida, provides results which are helpful in diagnosing anemia in that it provides a red blood cell count and mean red blood cells parameters characterizing the entire population of cells; more specifically, the Wintrobe indices of mean cells volume, mean cell hemoglobin and mean cell hemoglobin concentration along with the number of red cells per cubic millimeter. However, no differentiation between abnormal or normal red blood cells is achieved with the Coulter Counter. Furthermore, the hemoglobin content for individual cells is not determined and measures of dispersion and skewness are not performed. Finally, no automatic categorization or quantitative direct association with known anemias or other pathologies is automatically included as part of the analysis.

Heretofore, some offline experimental work has been performed on image processing of erythrocytes.

One of these works, "Bentley, S.A. and S. M. Lewis, The Use of an Image Analyzing Computerfor the Quantification of Red Cell Morphological Characteristics, Brit. J. Haemat. 29:81,1975", describes an off line analysis of dried and stained red blood cells of a total cell population measuring three red blood cells parameters by an image analysis technique. This analysis is similarto the Coulter Counter analysis in that the parameters measured were from the total population of cells being analyzed, and were analogous to the 20 Wintrobe Indices. The drying of the red cells introduced artifacts, and there was a lack of central pallor, or internal red cell analysis to provide a highly refined classification. Moreover, there was not disclosed the capability for differentiating between and classifying normal red blood cells from abnormal red blood cells.

Likewise, there was no capability disclosed for the automatic classification of the red blood cells with respect to categories of anemias.

Measurements of normal erythrocytes without differentiation of any abnormal erythrocytes by image processing has been disclosed by J. E. Green and reported in a paper entitled "Green, J. E.,'Computer Methods For Erythrocytes Analysis', Proceedings of Symposium of Feature Extraction and Selection and Pattern Recognition, IEE Catalog No. 70C 51 C pp. 100, Argonne, Illinois, 1970". A similar type of paper reporting measurements on the red cells and how to measure their features without any classification thereof was disclosed in a paper entitled "Eden, N.,'Image Processing Techniques in Relation to Studies of Red Cell Shape', edited by M. Bessis, R. Weed and Leblond, Springer- Verlag. New York, pp. 141,1973".

In U.S. Patent No. 3,851,156 Green provides a technique for scene segmentation of stained red and white blood cells through the use of a color algebra technique. In so doing features of perimeter, size and color are generated for red blood cells. These are measured on the total population of cells and no classification into subpopulations is performed. Further a precise central pallor analysis is not considered, and means are not provided to acquire subpopulation statistical measures, such as the bivariate dispersion of the hemoglobin and size, nor is it indicated that they are importantto achieve an anemia categorization, or a profile of similarity measures to prototype anemias.

In short, none of the aforementioned systems has the ability to analyze cells by theirfeatures, particularly 40 the innerfeatures of cell pallor, to quickly classify the blood sample or report it as similar. As a result of recent special analyses, it has been possible to gain a better quantitative understanding of abnormal red blood cell subpopulations for different anemias and the relationship of the abnormal red blood cell subpopulations to normal subpopulations. It has now been found that a red blood cell sample of blood contains significantly more information concerning the type of anemia present than heretofore known or thought.

Thus, with the present invention, it has been possible to quantify and identify blood from patients with anemia by variables or measures of characteristic values of subpopulations, such as size, hemoglobin content, percentages of cells, and other parameters such as measures of dispersion and skewness for certain single and combined parameters over different subpopulations as well as the population of cells as a whole. 50 These measured properties of the subpopulation of cells and the population as a whole provide a robust description of a patient's blood sample.

Also, with the present invention, it is now possible to compare these red blood cells descriptors relative to characteristic values for the red blood cells in a standard normal blood sample and to those typically found in each of a plurality of recognized kinds of anemia. It is also possible to generate an indice of the relative closeness of the blood sample to one or more standard anemias so that the clinician is given a powerful quantitative relationship to aid in his diagnosis. Additionally, by testing the patient's blood at differenttimes, particularly after successive treatments, one can generate a series or indices of the patient's blood relative to each type of anemia. Thus, it is possible to determine a patient's progress to see if his blood is deteriorating or is progressing towards a more normal blood. With expanded and accepted usage of the present invention, it is thought that some of the othertime consuming, painful and/or expensive tests, above discussed, and now commonly used in the diagnosis of anemia, may be eliminated. Therefore, and as will be explained in greater detail hereinafter, a real time analysis of the red blood cell sample and to determine its makeup and a comparison to standard types of anemia is feasible.

To this end, the present invention uses multiple 10giG systems operating simultaneously and in a A 1 w& 3 GB 2 093 586 A controlled relationship one with another to divide and perform tasks therebetween. Thus, the ability to analyze hundreds of red blood cells and to extract their various features and then to define the parameters for the subpopulations for comparison with anemia standards can be accomplished by the division of the functions and tasks between these simultaneously cooperating logic systems, all as will be explained in greater detail hereinafter.

The invention will now be further illustrated and described with reference to the accompanying drawings, in which:

Figure 1 is a perspective view of an apparatus for practicing the method of blood analysis and embodying novel features of the invention; Figure 2 is a block diagram showing the operation of the apparatus illustrated in Figure 1; Figure 3 is a block diagram of the preferred process for analyzing and classifying blood cells; Figure 4 illustrates a scanning technique for locating cells and determining the boundary points of cells in an image; Figures 5a, 5b, and 5c are flow charts of the preferred classification technique for classifying the blood cells into mutually exclusive subpopulations.

Figure 6 is a diagrammatic view of a model for red blood cell central pallor measurement.

Figure 7 illustrates a chain code description and analysis method for three diagrammatic red blood cell outlines; Figure 8 is a block diagram of the preferred process for determining whether a cehis round; Figures 9a, 9b, and 9c are graphs illustrating thickness/density profile measurements for three different, 20 typically appearing cell types, measured in two orthogonal directions. These profiles are used to measure the cell central pallor features and target cell features, Figure 9a illustrating a "flat" cell having little or no central pallor development; Figures 10a, 10b, and 10c are graphs illustrating the profiles of the cells of Figures 9a, 9b, and 9c with the peaks and valleys of each profile labelled; Figure 11 is a schematic of the preferred process for accumulating red blood cell subpopulations parameters; Figures 12a, 12b, 12c, 12d, and 12e are schematics illustrating the preferred process of computing the subpopulation characteristics from the accumulated values from a plurality of cells; Figures 13a, 13b, 13c and 13dare graphs of bivariate distributions or red blood cells subpopulations exhibiting normal and anemic characteristics; Figure 14 is a schematic block diagram of the preferred process for generating a similarity, or n-space distance, measure between a given set of measurement values describing a blood and a stored set of characteristic values, prototypic of various anemic conditions, or the normal blood.

Figure 15a is a graph of population distributions of the cell circularity shape measure, illustrating the 35 differences in skewness of these distributions, over all subpopulations, for normal blood compared to that of sickle cell anemia; and Figure 16 is a graph of population distributions of individual cell central pallor measurements, forthe two subpopulations, biconcave and spherocyte, illustrating the differences in mean values and dispersion for bloods of spherocytosis, normal and iron deficiency anemia.

Forthe purposes of illustration of the invention, the drawings illustrate a method and apparatus for automatically classifying red blood cells and the analyzing the relationship of the patient's blood sample to at least one recognized category of anemia orto a normal blood, or to a red blood cell disorder other than an anemia. More specifically, individual red blood cells are automatically examined and classified into different 46 subpopulations such asl for example, a spherocytic cell subpopulation, an elongated cell subpopulation, an 45 irregular shape cell subpopulation, a target cell subpopulation and a generally round and biconcave cell subpopulation, and then a plurality of characteristic values are generated for the patients subpopulations and population of cells as a whole for comparison with reference characteristic values which define a recognized anemia. By way of example, the selected characteristic values, which identify a given anemia, have been found and are given hereinafter for the following anemias: iron deficiency, chronic disease, B-thalassemia, megaloblastic, hemoglobin SS, hemoglobin SC and spherocytic; and likewise, reference characteristic values defining a normal blood, i.e., substantially all normocytic cells orthe like, have been developed and will be given hereinafter.

Also, as will be explained in greater detail hereinafter, a population or subpopulation dispersion measure of the red blood cells in a patient's blood relative to characteristics such as, for example, hemoglobin, mean 55 cell size (or area) shape, and central pallor, may be reported to the clinician. For example, broadly speaking, the bivariate red cell distribution of size and hemoglobin content for each cell is generally in the form of an elliptically shaped profile, as best seen in Figures 13a-1 3d and having axes at 45' and 1350. The length and width as measures of the bivariate dispersion, and the location of the profile by measures of the mean values, may be reported to provide the clinician with an impression of the patient's total cell makeup. 60 Similarly, measures of central tendency dispersion and skewness of pallor and shape are provided to further quantify the total cell makeup on the total cell population, or on subpopulations, as illustrated in Figures 15 and 16.

The patientwho has anemia generally is experiencing difficulty in either manufacturing new normocytic red blood cells or his existing red blood cells are being destroyed atan abnormal rate or by an abnormal 65 4 GB 2 093 586 A 4 process. Red blood cells typically have a life of about 120 days and their generation, growth, and death is a continuous process. An anemic disorder generally manifests itself in blood cells having unusual sizes or shapes relative to a normal red blood sample, which predominantly contains round normocytic red blood cells, or in blood cells having hemoglobin characteristics differing from the hemoglobin characteristics for normal blood cells. Thus, since the currently existing cells, for example in a normal patient who has just developed an underlying disease process leading to anemia, have a life span of 120 days, new cells with differing characteristics will tend to produce a wider dispersion of population measurements, as in Figure 13c compared to Figure 13a, when sampled and examined by the precise measurement techniques described herein. Such information provides the clinician with a knowledge of the presence of any previously visually estimated anisocytosis, i.e., a large measure of cell variations in size, and also as well the 10 variations in hemoglobin content. The mean cell hemoglobin and mean-cell size information locate the central tendency of the cell distributions in Figures 13a-13d.

In the embodiment of the invention described herein, each of the specified anemias is identified by 16 stored parameters or properties. The patient's blood is analyzed on an individual cell basis with each cell being classified into a subpopulation and then parameters such as mean cell size, mean cell hemoglobin, 15 and the percentage of cells in the subpopulation of the total cell population are generated to give subpopulation results. Also, a plurality of other measured properties or parameters of the patient's blood are generated from the subpopulation parameters, to total 16 parameters to define a set of reference characteristic values, i.e., an n-space location, and a calculation is made of the closeness of the patient's blood location relative to the eight reference character values or n- space locations for the seven anemias and 20 the normal blood. A report of the closeness of the patient's blood sample relative to these standard anemias provides the clinician with a statement as to what type of anemia, if any, the patient has, or how similar it is to a known type. Then, after the patient's treatment, the clinician is able to make later analyses and achieve new quantified results showing whether the patient is progressing towards a more normal blood or is deteriorating.

To achieve the analysis of the individual cells and the classification of same into subpopulations and the comparison of the blood subpopulations variables to those defining a specific anemia on a real time basis, the preferred equipment employes first and second logic systems which operate simultaneously and in a controlled manner so as to proportion the work and efforts therebetween. Also, as will be explaind in greater detail, the present described apparatus and method include a number of powerful and novel techniques and 30 means of and for cell classifying and analyzing which result in an efficient and less expensive method and apparatus for doing the red blood cell analysis. For instance, it is recognized that a normal blood sample generally will have a very high percentage of round cells with identifiable central pallor which can be grouped into a common subpopulation called a "bioconcave" cell subpopulation and that seven different anemias can be identified when using only four other subpopulations with this biconcave subpopulation. It 35 is to be understood, however, that the application of the invention in this context is not limited to any subpopulations described or defined herein, as the particular names and makeup of supbopulations may be varied and still fall within the purview of such application of the invention.

As shown in Figures 1 and 2 of the drawings, an apparatus 10 comprises a microscopic digital image processing and pattern recognition system which analyzes a mono layer of red blood cells on a microscopic 40 slide 12 with the cells being spaced from each other to ease the automated classification thereof. Suitable high resolution microscope optics 14 form an optical image for each red blood cell on a vidicon television camera tube or other detector 16 which converts the scanned electronic charged distribution of the optical image point by point into a numerical or digitalized image representing the optical transmission of the points in each image. The output of the vidicon camera is applied to digitizer electronics 20 which includes an 45 anaog to dig ita [-converter which is connected to an image processing logic 22 which controls the digitizer electronics 20 and receives and stores the digitized cell images into a memory store. The image processing logic 22 operates on the digitized cell images in a manner that will be hereinafter described which includes cell feature extraction and cell classification.

A suitable stage motor means 24 is provided and controlled by stage motor electronics 26 which are in 50 turn controlled by a master control logic 28. The stage motor 24 is provided to shift the slide 12 in order to iteratively process different image areas of the blood specimen on the slide. To control the focus of the microscope, a focus control motor means 30 is connected to the microscope and is operated by focus motor electronics 32 which are also controlled by the master control logic 28 by means of the focus prameter electronics 34. Focus control of slides for image analysis is well known in the art, e.g., U.S. Patent No.

3,967,110.

The apparatus 10 shown in Figure 1 includes a housing 38 having a cover 40 enclosing the microscope optics 14 and the television vidicon 16. An upper section 42 of the housing 38 houses the control switches of the apparatus, the next lower section 44 houses the master control logic 28 with the next two lower portions 46 and 47 of the housing containing the memory store for the image processing logic 22 and master control 60 logic 20 and the motor electronics 26 and 32. A terminal 48 is connected to the master control logic 28 and has a keyboard 50 for input of identifying information about the specimen or for other instructions. A monitoring screen 52 provides a visual display of the final report, and preferably a written printout is also made by a printer means 54 to afford a permanent record. A TV monitor 55 provides desired pictorial displays. The TV camera electronics are housed in a section 49 below the monitor. The next lower section 51 65 A l^ GB 2 093 586 A houses the analog to digital converter with the first section 53 housing the image processing logic 22. The results of the red cell analysis may also be transmitted for storage in a medical computer data bank.

Red blood cells may be examiner such that normal cells are distinguished from abnormal cells and classified by the apparatus 10 into subpopulations automatically in a detailed fashion heretofore not possible by a manual/visual examination of cells. Also, each of the red blood cells being examined may be 5 classified into mutually exclusive subpopulations and reported out so that the presence of a minor numbef of abnormal cells is not overlooked or forgotten and so that accurate parameters about a given subpopulation may also be provided. For the first time, the individual red blood cells may be examined individually for the hemoglobin contents. Thus, a report may be made not only of the kind of cells found in the subpopulation but also of their number and their hemoglobin characteristics. Advantageously, the individual red blood cells may be analyzed and classified with less subjectivity into a large number of mutually exclusive subpopulations (Table 1) such as biconcave (round cells with central pallor), elongated cells, targets and irregular cells (cells not fitting into any of the above classifications).

The preferred hemoglobin characteristic gathered from the analysis of the hemoglobin contents of the individual cells within a given subpopulation and reported out is the mean cell hemoglobin (MCH) for a given subpopulation of cells, such as shown in Table 1. In addition, to the hemoglobin parameters, the individual cells are counted for each subpopulation to provide their respective percentages of the total population; and likewise mean cells area (MCA) for each subpopulation may also be reported as shown in Table 1. It has been found to be helpful in detecting subnormalities in blood samples to determine multivariate distributions of the red blood cells in particular subpopulations of a sample with respect to a plurality of quantifiable features. In this regard a bivariate distribution is shown in Figure 13a as a distribution of round biconcave cells, with respect to a preferred quantifiable feature cell area on one axis and the cell hemoglobin content on the other axis.

By means of a measurement and analysis procedure to be described, parameters are reported which describe this distribution with regard to its central disposition, or mean values overthe plurality of variables, 25 and its variability, spread, or dispersion. The mean cell area and mean cell hemoglobin describe the center of the distribution and are reported as shown in Table 1. Two other statistical parameters EV1 and EV2 are reported in Table I and describe the variance of the dispersion of the distribution in the orthogonal directions of its major and minor eliptical spread. EV1 and EV2 stand for eigenvalue 1 and eigenvalue 2, respectively, and describe the dispersion or spread of the distribution. If the points of the distribution are thought of as defining an ellipse, then EV1 and EV2 can be thought of as relating to the length and breadth of the ellipse. Advantages derived from reporting parameters relating to a distribution of

a particular subpopulation will be more fully described hereinafter.

Other parameters reported in Table I include the mean pallor volume (PAL) for the biconcave and spherocyte cells as well as the standard deviation for the distribution of the biconcave and spherocyte cells 35 with respect to central pallor volume. The pallor volume standard deviation (PSD) is a parameter which describes the variance of the distribution of this measure over these subpopulations of cells. Another parameter reported is the skewness (SKW) which measures the skewness of the distribution of all the cells with respect to the quantifiable feature (perimeter of the cell) to area of the cell.

This data has been unavailable prior to this invention from any commercial instrument, or in any other fashion, such as from special research instrumentation. The closest analogous instrument is the Coulter Counter (Coulter Company, Hialeah, Florida) which is unable to classify red blood cells into subpopulations and which reports the mean cell size and mean cell hemoglobin for the entire population of red blood cells.

As seen in Table 1, the described apparatus and method is also capable of reporting the total population, or average mean cell hemoglobin as well as the average mean cell area (which is related to the mean cell size) 45 in addition to the other parameters suggested. In that table these are denoted in the line with AVERAGE parameters.

Thus, as indicated above, the described apparatus and method will be described as having the ability to classify red blood cells into the several mutually exclusive subpopulations set forth in Table 1. The subpopulations listed are the preferred subpopulations for classifying blood with respect to recognized categories of anemias but there may be other subpopulations defined. The mean cell area (MCA) is reported in microns 2 with the mean cell hemoglobin (MCH) reported in picograms (pg).

The several subpopulations described and their associated parameters hereinafter are:

6 GB 2 093 586 A 6 TABLE 1

96.8 BICONCAVE MCA MCH MCA 50 0.5% Spherocytes 47 30 5 MCH 31 0.2% Elongated 5 2 EVII 42 2.3% Irregular 38 23 10 EV2 2 0.2% Targets 57 34 AVERAGE 50 MCA 31 MCH 17 PAL3 PSD 8SKW 0.9 Normal 4.1 Megloblastic 15 4.2 Iron Deficient 6.2 Hemoglobin SS 2.5 Chronic Disease 3.8 fi-Thalassemia 4.8 Hemoglobin SC 4.9 Spherocytic In another technique, samples of blood may be analyzed and thereby classified by "similarity" or "distance" measures being reported to compare said sample to recognized categories of anemia or normal bloods. In the preferred embodiment, 24 parameters are measured for the subpopulations of the sample of 25 blood taken from the patient. Of these, 16 are used for the tested sample of blood and define a point in this 16-space. tonsequently, the typical parameter values for a particular anemia also define a point in the 16 parameter space. The distance, from the point representing the values forthe sample blood taken from the patient, to each of the points representing the typical parameter values for each of the categories of anemia, may be determined. Thus, a physician would be able to determine which of the categories of anemia the sample of blood taken from the patient most closely resembles and could make a diagnosis from that information. Alternatively, simple decision logic could point out the most probable diagnosis. The normalized distance of the parameter values for a sample of blood is shown for a normal category of blood as well as the recognized categories of anemia in Table 1. As seen in Table 1, this particular sample of blood is closest to normal since 0.9 is less than any other distance reported and therefore the blood most closely 35 resembles normal blood.

With reference now to another technique, a multiple parallel logic architecture has been found to provide the rapid processing necessary for eff icient analyzing of cells on a slide. Thus, in the preferred embodiment, there is provided a first processing means, the master control logic 28, and a second processing means, the image processing logic 22 as shown in Figure 3. The analysis of the cells on a slide requires a sequence of 40 operations to performed, and since one operation often requires the results of a previous operation, there are provided synchronizing means for synchronizing the processors so that the results necessary to perform a particular operation are available when that operation is begun.

Figure 3 illustrates the specific interrelationships between the master control logic 28 and the image processing logic 22. Because of this multiple parallel logic or architecture, the master control logic may proceed with one task or operation while the image processing logic is proceeding with another operation.

As seen in Figure 3, the operations carried out by the master control logic 28 are listed in the lefthand column with the operations of the image processing logic 22 in the righthand column. The master control logic, after clearing its associated accumulators, proceeds to operation 56 in which a start signal is sent to the image processing logic and thereafter continues to operation 58. The image processing logic meanwhile is 50 waiting for the start signal (operation 60) from the master control logic. Upon receipt of the start signal, the image processing logic 22 proceeds to operation 62 which includes digitizing the image produced by the vidicon camera 16 (Figure 2). Upon completion of the digitizing, the image processing logic sends a "digitizing done" signal (operation 64) to the master control logic indicating the completion of the digitizing process and proceeds to operation 66. The master control logic operation 58 is currently waiting for the "digitizing done" signal and upon its receipt proceeds to move the stage (operation 60) on which the slide rests so that a new field of cells may be imaged since the previous field has already been digitized by the image processing logic 22. The optics 14, Figure 2, are providing an imaging means of the cells on the slide.

The stage motor drive 24, and the focus motor drive 30, and their associated electronics, are controlled by the master control logic 28. After moving the stage so that a new field may be imaged, the master control logic proceeds to operation 70 wherein the field is focused and then proceeds to operation 72.

After transmitting the "digitizing done" signal, the image processing logic scans the digitized image for a cell boundary point (operation 66). If a cell boundary point is found (operation 74), the image processing logic extracts the cell's boundary and features (operation 76) and classifies the cell as to its proper subpopulation (operation 78).

Z J01 7 GB 2 093 586 A 7 The image processing logic then returns to operation 66 and continues scanning the image for another cell boundary point. The scanning, feature extraction, and cell classification operations will be described in more detail below. If the logic section 74 determines that a new boundary point has not been located, then the image processing logic proceeds to operation 80 wherein the features of each cell located as well as each cell's subpopulation classification is transmitted to the master control logic which will be in the process of executing operations 68, 70, or 72. The transmittal of the information is on an interrupt basis, i.e., should the master control logic be in the process of controlling the imaging means (operations 68 or 70), the master control logic will interrupt these operations and store the information received from the image processing logic before proceeding with moving the stage and focusing the microscope. However, if these operations have already been completed then the master control logic proceeds to operation 72 wherein the master control 10giG waits for the data to be transmitted from the image processing logic. In response to the receipt of the data, the master control logic will transmit an acknowledge signal (operation 82) to the image processing logic and then proceeds to operation 84 wherein the subpopulation data for each subpopulation is updated, as will be more fully explained below.

Upon receipt of the acknowledge signal, the image processing logic proceeds to digitize the image of the 15 new field that has been moved into view by the master control logic. The master control logic, upon completing the update of the subpopulation data, determines at logic section 88 whether N, the total number of cells processed, is equal to 1000. If 1000 cells have not been processed, the master control logic returns to operation 58 and waits for the "digitizing done" signal from the image processing logic, otherwise the master control logic calculates the subpopulation parameters (operation 90), proceeds with an anemia 20 classification (operation 100), and prints the results (operation 102), as will also be more fully explained below.

Thus, because of the dual processor architecture, the master control logic is free to control the imaging means wherein a newfield is brought into view to be imaged while the image processing logic is proceeding with the digitizing and analyzing of the image from the previous field. Similarly, while the master control logic is accumulating the data extracted from the image by the image processing logic, the image processing logic may simultaneously digitize and analyze a new image provided by the new field which has been brought into view by the master control logic. It should be noted that although for purposes of illustration only one image processing logic is described as associated with the master control logic, it is capable of utilizing information from several image processing logics operating in parallel and independently on 30 different images.

The present invention is directed to the optimization of the time of analysis as well as the number of features used in the classification logic so thatthe amount of storage and classifying techniques may be reduced substantially along with equipment requirements therefor. With an optimization of analysis time for classification, there is a danger that the reliability and accuracy of the classification are compromised.

Despite this, a-relatively foolproof feature set and classification logic has been invented for a large number of subpopulations such as those shown in Table 1. The preferred classification features are size, hemoglobin content, spicularity, roundness, elongation, central peak height (if present) from cross-sectional cell scans, and central pallor. By suitable combinations and analyses of such features, it is possible to differentiate from normal blood and to identify biconcave round cells, spherocytes, target cells, irregular-shaped cells, and elongated cells.

In the preferred method and apparatus, the cell classifications are achieved by an image processing and pattern recognition with great accuracy and reliability by rendering white blood cells and other artifacts substantially invisible to the optics 14 by using a light having an optical wavelength of about 415 Nanometers. At this optical wavelength, the red blood cells are relatively contrast enhanced to the utraviolet 45 Vidicon camera without staining, while the white blood cells and other formed elements are substantially invisible. The staining of the red blood cells prior to being analyzed by a microscopic image processing technique has been found to be a time-consuming process, as well as undesirable in that the staining may introduce a number of stained artifacts which detract from the accuracy of the analysis. Futhermore, many of the stains are not stoichiometric in the representation of hemoglobin concentration according to density, thus distorting the quantization of the hemoglobin content of the cell on a per-cell basis. By rapidly preparing the specimens in a monolayer and fixing with a formaldehyde vapor prior to the drying of red blood cells, and by not employing a time consuming staining to contrast enhance the cells, as in white blood cell analysis, specimens may be quickly prepared and analyzed accurately (the formation of artifacts by distortion of the central pallor being avoided).

The location of the cell image and the identification and feature extraction has been greatly simplified as described belowto locate and define the cells by a boundary procedure which defines the cell in the form of an octal chain code. The use of octal chain codes as an image processing technique is described in a paper by H. Freeman, "Computer Processing of Line-Drawing Images", ACM Computing Surveys 6:57,1974. As will be explained in greater detail, the octal chain code allows feature extraction as to: (1) cell size, (2) perimeter length and roundness shape measure (3) irregular shape measure, and (4) elongation shape measure. This is followed by extracting the summed density or hemoglobin feature, and then by extracting cross-sectional scans (thickness/density profiles) for central pallor measurement and target cell measure ment. Finally, inner central pallor boundaries are determined and features analyzed for more precise target cell identification.

8 GB 2 093 586 A 8 After having extracted these identifying features, the cells are then categorized by a classification means. The preferred classification means (Figures 5a, 5b and 5c) comprise either a digital logic system of electrical devices or a programmed microprocessor which uses Boolean logic to classify the red blood cells.

Referring now in greater detail fo the specific features of the illustrated embodiment of the invention, the images of the cells are digitized (operation 62 of Figure 3) in a manner known to the art, e.g., U.S. Patent No. 3,883,852 as a television digitizing system. Magnified blood cell images are obtained by using microscope optics with ultraviolet illumination, arranged to provide a 0. 23 pixel resolution in the image plane. A pixel is a picture element having a specific location in the digitized image stored in the memory analyzer.

Referring now to Figure 4 which illustrates in greater detail the operation 66 (Figure 3) by ihe image processing logic, an original microscopic image which had been digitized is stored as represented by the 10 image 108 for the purpose of further analysis. This analysis is carried out by the image processing logic and is represented by the blocks indicated at 115 which comprise the operations 76 and 78 (Figure 3). In this preferred embodiment of the invention, individual cells and 110 and 112 in a digitized image 108 are located by a technique in which a raster scan is made of the digitized image to locate objects above a critical threshold, such as illustrated for cell 110 in block 113. The boundary of the cell is traced by examining the neighboring pixel elements by a counterclockwise search, by techniques which are well known in the art.

One such technique is disclosed in U.S. Patent No. 3,315,229. During this counterclockwise boundary tracing operation herein, the picture element at the "top" of the cell, pixel 1 14a, which is usually the pixel located first, and the one at the "bottom" of the cell, here pixel 11 4f, are stored for reference in the later analysis. The analysis process then proceeds to extract features and to classify the located cell into one of a plurality of 20 subpopulations, as in block 115, and as described in detail later.

The raster scan of the digitized image is then continued from the bottom pixel 1 14f to hit the next digitized cell 112 by impacting a pixel 11 2a which is above the threshold as seen in block 116. After the boundary is traced and the features for this cell are extracted and the cell is classified, the raster scan continues from the bottom pixel 11 2b, and, as seen in in block 118, no more cells are located in the image field. At this time, the 25 image processing logic transmits the cell features and subpopulation classifications to the master control logic (operation 80) as shown in Figure 4.

The initial image processing done by the image processing logic outlined in Figure 3 is shown in greater detail in Figure 5a. After the image has been digitized (operation 62), the image is scanned to locate a cell (operation 66) and the boundary is traced as explained above.

During this boundary tracing operation, octal chain codes are formed in an operation 119. The outer boundaries, defining a cell, are processed in the following manner. Each pixel element defining the boundary is stored in a list as a series of numbers indicating a line description of the cell. For instance, referring Figure 7, a digital image of cells as defined by their boundary pixels 120 are illustrated.

As is well known in the art, e.g., as described in "Bacus, J. W. andJ. H. Weens,'An Automated Method of 35 Differential Red Blood Cell Classification Application to the Diagnosis of Anemia', Journal of Histochernistry and Cytochemistry, 25:7,1977", a plurality of features F1-F4 can be computed from this chain code. The details of this computation are fully described in the aforementioned publication, which is hereby incorporated by reference as if fully reproduced herein.

The above features are combined with otherfeatures for use in the classification of the cells. In this regard, 40 the following features are used herein:

z 9 TABLE 11

Feature Description

Howdetermined GB 2 093 586 A 9 F1 Area size Number of pixels enclosed 5 by cell boundary F2 Shape (circularity) (Number of perimeter pixels)2/area F3 Shape (spicularity) Number of "spicules" on 10 boundary F4 Shape (elongation) Comparison of orthogonal boundary chain code 15 orientations F5 Grey levels Sum of grey levels as a measure of Cell Hemo globin 20 F6 Pallor (volume) The percentage volume of the central pallor F7 Central peak The height of the central 25 peak of a 3-peaked profile of a cell F8 Pallor (depth) Fora 2-peaked profile, the difference of the 30 valley from the peak heights F9 Pallor (circularity) (Number of pallor boundary pixels)2/area of pallor 35 As indicated above, features F1-F4 are calculated in an operation 124 by the image processing logic as shown in Figure 5a. Feature F1 relates to the area or size of the cell as determined by the number of picture elements or pixels that are enclosed by the cell boundary. Feature F2 is the (boundary perimeter)2 /area and is of assistance in classifying round and non-round objects. A round object would have a theoretical value of 4n 40 and non-round objects have greater values.

In actual practice the value of the perimeter squared divided by the area for round digitized objects varies as a function of the number of pixels, and in addition always involves quantization error, such that in practice for quantized circles the value approximated 14.0, and is a better approximation to this reference number as the number of pixels, or size, of an object increases. For total cell areas above 500 pixels, the quantization 45 error is within .2 units.

Features F3 and F4 relate to the spicularity and elongation shapes, respectively, F3 being a count of the number of spicules in a chain code boundary, and F4 measuring the non- roundness due to elongation of the boundary, as shown in Figure 7. Feature F5 is the integrated optical density of the cell (operation 136). It is the sum of the grey levels within the enclosed boundaries of the Feature F6, which is a measure of the pallor 50 volume, assists in distinguishing cells with large pallors, such as hypochromic cells from normocytes.

Feature F7 is equal to the larger of the two central peaks of two crosssectional orthogonal 3-peaked thickness/density profiles, either having a central peak, and is used to detect target cells. Feature F8 is a measure of the depth of the central pallor, as determined from two crosssectional orthogonal, 2-peaked thickness/density profiles. Feature F9 is a measure of the degree of roundness of the pallor itself, and is also 55 used in distinguishing target cells.

The logic decisions for determining the various features that have been briefly described are carried out by the image processing logic using the logic flow chart shown in Figures 5a, 5b, and 5c. The logic decisions are made using the various features together with threshold values that are identified at T1 through T1 1. The thresholds T1-T1 1 are described in Table V and specific values are also provided. As shown therein, the thresholds are used by the logic with the various features in making logic decisions leading to the classification of the cell of interest in accordance with the flow chart shown in Figures 5a, 5b, and 5c. In this regard, Figures 5a, 5b, and 5c illustrate various decisions that are made on the basis of various features either exceeding or being less than certain threshold values as will be specifically described.

Referring to Figure 5a, an object that is located is examined by logic section 138 to determine if it is 65 GB 2 093 586 A sufficiently large to be a cell, rather than a noise or dirt artifact, and thus is to be further analyzed. If feature Fl, which is the size or area of the object under consideration, is less than the threshold value Tl which may be a value of about 6 micronS2, then the object is not considered by the decision logic and another object will be located for analysis and classification. However, if the area of the cell is greater than the threshold value Tl, feature F5 is computed in operation 136 wherein the hemoglobin content of the cell is determined. This is simply a summing of the grey levels inside the boundary of the chain coded cell and then dividing by a conversion factor 1290 or thereaboutto convert the grey level measurements to picograms of hemoglobin per cell.

For this purpose the electronics generating the television signal and digitizing said signal should be 10 adjusted to produce grey levels corresponding to the following optical density at 418 nanometers:

TABLE Ill

Optical density Greylevel 15 134 17 294 35 20.403 52 20 505 43 605 57 25 Also, for calculation of hemoglobin and the area, the optics and television electronics should be adjusted such that round objects of the following dimensions produce the given number of pixels.

TABLE IV size p? Pixels ill 1850 93 1850 35 77 1283 58 967 40 34 567 23 383 17 283 45 4 67 The decision logic then operates to determine whether the cell is round or non-round. This is performed by a logic section indicated generally at 140. The logic sections 140 is shown in Figure 8 to include logic subsections 142,144, and 146. The subsections 142,144, and 146 are operable to jointly make the roundness determination with the features F2, F3, and F4 being examined with respect to thresholds T4, T5, and T6. If the cell has a small roundness value, a small spiculated value, and a small elongated value, then it is considered to be round and is passed on to the next operation 148 (Figure 5a) which is the first step in the target cell analysis and central pallor analyis. Similarly, if it is determined that the cell is not round, then logic 55 subsection 150 (Figure 5a) operates to determine if the size of the cell exceeds an upper boundary threshold T2, and if it does, the cell is not further analyzed and a new cell will be considered. The effect of the subsection 150 is to eliminate double cells such as that shown in the pictorial representation 152. It should be appreciated from the pictorial representation that such a double cell would not pass the roundness test, but it is also not a non-round cell of the type for cells of classes 3 and 4. Thus, it cannot be accurately classified and 60 it is forthis reason that the subsection 150 eliminates such cells from further consideration.

As previously mentioned, the roundness of the cell is determined by feature F2 which will have a value of 14.0 for a perfect circle and will increase as the shape of the cell departs from circular. Thus, the threshold value T4 is chosen to reflect reasonably good circularity and if the feature F2 exceeds the threshold T4, that is an indication that the shape is not circular, hence the logical flow to subsection 150 indicating that the object65 -W _Z 11 GB 2 093 586 A 11 is not round. If feature F2 is not greater than thresho I d T2, it is one indication that the cel I is round and if the decision from the subsections 144 and 146 also indicate adequate roundness, the logic flow then proceeds to logic subsection 148 (Figure 5a).

In operation 148 thickness/density profiles are extracted from the cell image. These profiles are illustrated in Figures 9a-9c and 10a-10c. A thickness density profile is determined by the grey levels of the pixels along a particular direction across the cell image. As noted earlier, the grey level of a pixel is determined by the hemoglobin density at that point. It has been found that the grey level of the cell at a particular point is related to the hemoglobin density and the cell thickness at that point. Two such thickness/density profiles, profile a and profile b, are shown in Figure 9a for a biconcave cell determined in two orthogonal or transverse directions, a and b. Two profiles each are also illustrated in Figure 9b and 9c for a target cell and a 10 spherocyte cell. As seen in Figure 9b, one direction (direction a) practically missed the center area. Since these profiles are used to distinguish target cells (feature F7), two transverse directions are preferably analyzed. Thus for each cell, two cross-sectional profiles are determined wherein the profile relates to the thickness of the cell along the points of the cross sections.

A profile for each cell of Figure 9 is discussed more fully in connection with Figures 1 Oa-10c. As seen in 15 Figure 1 Oa, the profile has two "peaks", P1 and P2, and one "valley", V1. P1 and P2 are relative maxima of the profile of the cell with respect to the cell thickness and thus determine the two relative maximum thickness density points along the profile. V1 determines the relative minimum point of thickness density.

Similarly, the target cells have three relative maxima, P1, P2, and P3, with two relative minima V1 and V2, as shown in Figure 10b. The spherocyte has one peak, P1, and no valleys (Figure 1 Oc). These profiles are utilized 20 in a target cell analysis and a central pallor analysis as will be more fully explained hereinafter.

After the image processing logic extracts the thickness/density profiles for the cell, it proceeds to the target cell analysis performed by the logic section, referred to generally at 156 of Figure 5b. The first step of the target cell analysis is to smooth the two profiles, profile a and profile b, as shown in operations 156 and 158, which is performed by the image processing logic before proceeding to a logic subsection 160. The logic subsection 160 determines whether a profile has three peaks and if so forwards it to an operation 162 which determines half the average of the two non-center peaks, P1 and P3, or "LEV1 ". A logic subsection 164 determines whether the two valleys, V1 a and V1 a, are less than LEV1 and if so then the cell located might be a target cell and the image processing logic proceeds to examine profile b. If not, then the valleys are not deep enough to profile a to be a target cell, so the center peak, P2a, is set to zero in an operation 166 and profile a is smoothed to two peaks or less in an operation 168.

After profile a is examined, profile b is examined for three peaks in a logic subsection 170. If the logic subsection determines that profile b has three peaks, it is forwarded to an operation 172 and logic subsection 174wherein the two valleys, V2a and V2b, are compared to LEV2 which is half the average of the two non-center peaks P1 band P3b as for profile a. If the two valleys are less than LEV2, then it is forwarded to operation 176 wherein the feature F7 is determined as to which is the larger of the two center peaks, P2a and P2b, of the profiles a and b. Feature F7 is compred to a threshold T7 in a logic subsection 178, and if larger, the cell is classified as a target cell (C5). In other words, if the larger ofthe two center peaks is larger than a certain threshold, then the cell is determined to be a target cell. If not, then the center peaks of the profiles are probably due to "noise" in the image video and digitizing and not due to a center area of a target cell. In that 40 case, both profiles are smoothed to two peaks or less in operations 180 and 183. However, if the logic subsection 174 determined that the valleys of profile b were not less than LEV2, then the profile b is forwarded to a logic subsection 184 which checks whether the center peak of profile a had been set to zero. If not, then profile a may have detected a target cell and thus P2b is set to zero and subsection 176 determines the maximum value for P7 as described.

If the center peak, P2a, had been set to zero, then neither profile has passed the tests at logic subsection 164 and 174 respctively. Thus the cell is probably not a target cell and profile b is also smoothed to two peaks or less at operation 182. However, some target cells might not be detected in this analysis, therefore, other tests are performed on the cell as will be explained later.

After the center peaks of profiles a and b have been examined as explained aboe, a logic subsection 186 determines whether profile a has only one peak. If so, the variables P1 a, P2a, and V1 a are set equal to each other in an operation 188. In either case, the image processing logic then examines profile b to determine whether it has only one peak, at the logic subsection 190. If profile b has only one peak, then the variables P1 b, P2b, and V1 b are set equal to each other in an operation 192.

Continuing with Figure 5c therein, a feature F8, which is the average value of the two valleys subtracted 55 from the average value of the four peaks of the two profiles of the cell, is determined by subsection 194. Then the cell feataure F1 is examined to determine whether the size of the cell is larger than a threshold T8 at a logic subsection 196.

If the cell is large, i.e. F1 is greater than T8, it is possible that the cell is a target cell despite the pevious target cell analysis and therefore another target cell analysis will be performed beginning in operation 198. 60 Therein, a variable LEV3 is set equal to one-half the value of feature F8 (operation 198).

Next, a search for the central pallor of the cell is initiated by searching a direction along the line from the top pixel of the cell through the center of the cell looking for a threshold condition, i.e., hitting a pixel which is below the threshold LEV3, before the center is reached. The chain code is then formed forthe central pallor boundary (operation 202). The pallor circularity feature F9 is then computed in an operation 204. F9 is 65 12 GB 2 093 586 A 12 calculated as the number of pallor boundary pixels squared divided by the area of the central pallor. F9 is then compared to a theshold value T9 at a logic subsection 206 to determine the circularity of the central pallor. This operation is necessary since the two profiles from the previous target cell analysis may have missed the central area as shown for the cel 1208. Thus, if circularity feature F9 is greater than the threshold T9, then the cell is a target cell, otherwise the cell is forwarded to the operation 209 wherein a feature relating to the size of the central pallor of the cell is computed.

The central pallorfeature is defined as the percentage volume of a cylinder, with the height and area of the cell under consideration, not occupied by hemoglobin. This is illustrated in Figure 6, where T represents the cell height orthickness, and 132 indicates the indented central pallor region. The cell area is known from previous analysis on that cell, i.e., Fl. Also, feature F5 is the sum of the grey levels for pixels enclosed by the 10 chain code defining the boundary of the cell. As noted above, the hemoglobin density is related to the thickness of the cell and in this manner the hemoglobin feature F5 defines a volume which is related to the thickness or volume of the cell. The cylinder height, or thickness M, is derived by using the average value of the peaks of the two thickness/density profiles of the cell, as:

T = Pla=P2a+1P1b+M 4 Thus, the volume of the central pallor maybe calculated as: T times the area of the cell (Fl) minus the hemoglobin content. Finally, the percentage pallor volume F6 is:

F6 = (TXF1-F5)x100% Tx171 After this feature has been computed, the image processing logic proceeds to a logic subsection 210 wherein the cell is distinguished between biconcave cells (Cl) and spherocyte cells (C2) as it has already been determined that the cell is not an elongated cell (M), an irregular cell (C4), or a taget cell (C5). The logic subsection 210 compares the percentage pallor volume feature F6 to a threshold value T1 0 and the pallor depth feature F8 to a threshold T1 1 and if either feature is less than its associated threshold then the cell is deemed a spherocyte cell (U), otherwise it is a biconcave cell (Cl).

Referring back to Figure 3. the feature extraction operation 76 and the cell subpopulation classification operation 78 have been completed for the cell that had been located in the image scan. The image processing logic will then continue scanning the image for another cell (operation 66) and if no other cells are found then the features for those cells located as well as the cells'subpopulation classifications will be sent to the master control logic in the operation 80.

While the determination of the various features and decisions contained in the logic diagram of Figures 5a, 5b, and 5c is carried out utilizing the threshold values contained in Table V, it should be understood that the threshold values are based upon empirical and statistical analysis and can be varied somewhat without appreciably affecting the eventual classification of the cells. It should also be appreciated that the threshold values are believed to be optimum values which have been fixed to maximize the accuracy of the classification.

z _f z 13 GB 2 093 586 A 13 TABLE V

Threshold Value Description

T1 6R 2 Size threshold for 5 artifact T2 54R2 Size threshold for double cells 10 T3 25 Elongation threshold T4 16 Cell circularity theshold T5 7 Spiculed threshold 15 T6 25 Elongation threshold T7 5 grey levels Target center peak heightthreshold 20 T8 47 R2 Size threshold for target cells T9 20 Pallor circularity 25 threshold T10 11% Pallor volume threshold T11 8 grey levels Depth of pallor threshold 30 Upon completion of the feature extraction and cell classification analyses for the cells located in the image, these features are transmitted to the master control logic as illustrated in Figure 3. After acknowledging the receipt of the data (operation 82), the master control logic proceeds to update subpopulation measurements for each cell class located in the image just analyzed (operation 84). A diagram illustrating the updating operation in greater detail is shown in Figure 11. A plurality of accumulators are provided to produce a running total of a plurality of measurements for the cell subpopulations or classes. Each accumulation is a function of one or more cell features, such as the cell feature value itself or the value squared, for example.

The cell feature values F1, F2, F4, F5, and F6 for a particular cell are provided as inputs to the accumulators together with the cell classification Ci to which the cell features pertain. After the measurements for the cell 40 have been accumulated, then the other cells in the image are similarly processed to further accumulate the measurements based on all of the cell's features.

Thus, the feature F2 (cell circularity feature) is provided at a line 212 to an accumulator 214. The accumulator 214 produces a running total S1, i.e., accumulates the measurement (F2 - 14.1)3 for all the cells located by the image processing in logic wherein F2 is the cell circularity feature (Table IV). This measurement is used in a later calculation which provides a parameter describing the skewness of the distribution of all the red blood cells located with respect to the circularity feature of the cells.

Also accumulated is the elongation feature F4 which is provided at a line 216 to accumulators 218 and 220.

The accumulator 218 sums the total (S2) of the feature F4 for all the cells which is used to calculate the average elongation for the cells. The accumulator 220 provides a sum or running total (S3) of the elongation 50 featu re F4 sq u a red, i. e., (F4)2, which is used to calculate a parameter describing dispersion, or variation of the distribution of the red blood cells with respect to the mean of the elongation feature F4.

Not all feature mealsurements are accumulated for each subpopulation. For example, the feaure F6 (pallor volume) is only accumulated for the bioconcave cells (subpopulation Cl) and the spherocyte cells (subpopulation C2). Therefore, in addition to the feataures for a particular cell, the subpopulation classification for the particular cell to which the features pertain is provided which is shown as Ci at line 222.

A plurality of logic utilize the input Ci to discriminate among the cell subpopulations. Thus, the cell classification Ci is provided to the inputs of a logic AND gate 224 and an AND gate 226 with subpopulation C1 constant (i.e., a 1) provided to the other input of the AND gate 224 and subpopultion C2 constant (i.e., a 2) provided to the other input of AND gate 226. The output of these AND gates are provided to an OR gate 228 60 which may enable the accumulators 230 and 232. The accumulator 230 provides a summation of the feature F6 (central pallor volume) as indicated by input lines 242, but only when enabled by the logic OR gate 228.

Similarly, the accumulator 232 accumulates the sum of the feature (F6)2 but only when enabled. Thus, the gates 224, 226, and 228 permit the accumulators 230 and 232 to accumulate the measurements derived from the feature F6 only when the feature had been extracted from a C1 or C2 biconcave or spherocyte class cell. 65 14 GB 2 093 586 A 14 The output of the accumulator 232 is provided at S5 which is used to compute the dispersion parameter of the distribution of spherocyte and biconcave cells with respect to the mean volume of the central pallor of the cells. The output of the accumulator 230 is provided at S4 which is also used to calculate the dispersion parameter and also to calculate the mean or average cenral pallor volume for the spherocyte and biconcave 5 cells.

Similarly, a logic AND gate 234 enables accumulators 236,238, and 240 when Ci at line 222 is equal to a 2, i.e., the cell features appearing on the feature lines 244 and 246 were extracted from a class C2 (spherocyte) cell. The accumulator 236 accumulates the feature F2 (cell area) which is provided at S1 1, which will be used to calculate the mean cell area parameter forthe cells in the C2 classification. The accumulator 238 provides at S12 the accumulated total of feature F5 (cell hemoglobin content) which is used to calculate the mean cell hemoglobin content for the class C2. The accumulator 240 provides a total of the number of cells in the C2 class, i.e., N2 equals the number of spherocyte cells located by the image processing logic.

Ina similar manner the total cell area for the elongated (M), the irregular (C4), and target (C5) cells are provided at S1 3, S1 5, and S1 7, respectively. The total of all cells' hemoglobin content for the elongated, irreguir, and target cells is provided at S14, S16, and S18, respectively. The total number of cells in each of 15 the above subpopulations is provided at N3, N4, and N5.

Likewise, the tota of all of the cells'areas forthe biconcave subpopulation is provided at S6, the total of all the cells' hemoglobin contents is provided at S7, and the total number of biconcave cells is provided at N1.

For additional accumulated measurements on the biconcave subpopulation, additional logic gates permit accumulators to dicriminate among the class cells. Thus, an AND gate 248 enables accumulators 250,252, 20 and 254 when the feataures appearing at the lines 244 and 246 have been extracted from a C1, i.e., a biconcave cell. The accumulator 250 provides the accumulated sum of the measurement (F1)2 at S8. The accumulator 252 similarly provides the accumulated total of the measurement (F5)2 at S9. Finally, the accumulator 254 provides the accumulated sum of the product of the feature F1 times the feature F5 (F1 x F5). The accumulated S9 and S10 are used to calculate parameters descriptive of the dispersion, or variation 25 of the bivariate distribution which will be further explained hereinafter.

Thus the features for each cell examined by the image processing logic provide the inputs to the logic described in Figure 11 for updating or accumulating measurements based upon the cell features with the particular measurements updated for each cell depending upon the subpopulation classification to which that particular cell belongs. The measurements updated by the logic of Figure 11 provide an intermediate step for the calculation of parameters which are descriptive of each subpopulation classification as well as parameters which are descriptive of multivariate distributions of cell subpopulations with respect to different cel I features.

Referring back to Figure 3, it is seen that at logic subsection 88 the determination is made whether a preset total of N cells have been processed. If not, the master control logic returns to operation 58 wherein it waits 35 for the "digitizing done" signal indicating that the image processing logic has completed digitizing the next field. If N cells have been processed, e.g., N =one thousand, then the accumulated measurements which had been updated as illustrated in Figure 11 for those N cells are used to calculate the parameters descriptive of the subpopulations (operation 90) which is illustrated in greater detail in Figures 12a through 12e.

The output S1 of the accumulator 214 (Figure 11) is used in the calulation of a dispersion parameter which 40 describes the skewness of a distribution. Herein, a distribution of all the cells with respect to the elongation feature (K. Skewness is calculated a ti -L Thus a logic subsection 256 having inputs S2 and N produces the skewness prameter:

(- 1) J_ SkW % N 3 The calculation of the skewness parameter is quite helpful in describing a population of cells. For example, a distribution of normal cells is shown in Figure 15, generally at 255. The distribution is with respectto feature F2 (circularity). Also shown is a distribution of sickle cell anemia cells, generally at 257. As can be 60 seen there, the distribution of sickle cells is greatly skewed toward the right, indicating a great number of elongated cells. Note, however, that the mode of both distributions is identical. Thus, the skewness parameter is a valuable comparison too[ for indicating anemias.

A logic subsection 258 having inputs S2 (the sum of the elongation measurements forthe cells) and N (the total number of cells) produces the mean cell elongation parameter (ELN).

-i z p A GB 2 093 586 A 15 The general formula for the dispersion in the form of the standard deviation of a distribution with respect to a variable X is given by:

1 Stch bQ4. = -i A logic subsection 260 produces the standard deviation of the elongation distribution of cells with respect 15 to the elongation features. The logic subsection 260 has in input S2 equal to m 1 F- Ftk 20 k%. % 9 _ ( F v.) I,rXIL (Figure 11) and an input S3 equal to and produces the elongation standard deviation (ESD) after the square root of the output has been taken by a logic subsection 262.

A parameter for the mean central pallor volume (PAL) of the biconcave and spherocyte cells is provided by a logic subsection 264 having inputs Ni (the number of biconcave cells), N2 (the number of spherocyte cells), 30 and S4 (the accumulated sum of the volumes of the central pallors of those subclassifications). A parameter of the distribution of the biconcave and spherocyte cells with respect to the central pallor volume, herein, the central pallor volume standard deviation (PSD) is provided by a logic subsection 266 having inputs S4 and S5 and a!ogic subsection 268 which takes the square root of the output provided by the logic subsection 266 to finally produce the parameter PSD in a manner similar to that of the parameter ESD.

A distribution of three different populations of cells, normal, spherocytic, and iron deficient, with respectto the feature F6, the percentage volume of central pallor is shown in Figure 16. It is important to note that the distribution of normal cells at 267 has the same mean value (PAL) as the distribution of iron deficient cells at 269, yet they have a different variation or standard deviation (PSD) in cengtral pallor volume. On the other hand the distribution of normal cells has the same standard deviation as the distribution of spherocyte cells 40 at 271 but a different mean value. Thus both parameters have been found advantageous in the classification of blood with respect to anemias.

Two other parameters, EV1 and EV2, are computed utilizing the accumulated sums S6-S1 0 and N1 and which are descriptive of the amount of dispersions of a bivariate distribution of the biconcave cells. The two variables of the bivariate distribution are the cell size and the cell hemoglobin content. Four such distributions are illustrated in Figures 13a-1 3d wherein the cell area defines the abscissa axis and the cell hemoglobin content defines the ordinate axis. Each "X" represents a biconcave cell with its location within the graph defining the cell's area and hemoglobin content. Thus, as can be seen in the four Figures, the cells are distributed mainly on a 45 line passing through the origin. The mean cell area (MCA) and the mean cell hemoglobin (MCH) define the center of each distribution. The values, EV1 and EV2r define the dispersion or 50 the amount of spread of the distribution in two principal independent axes. In particular, EV1 describes the amount of spread of the cluster or distribution along the direction at essentially 45', or along the line of major dispersion of the ellipse with EV2 describing the dispersion in a direction which is orthogonal or transverse, that is, 90', relative to the dispersion of EV1.

Referring to Figure 12a, a logic diagram is shown for the computation of the parameters EV1 and EV2. The 55 general formula for computing the variance of a distribution with respectto a variable is similar to that given for the standard deviation. The variance of the distribution with respect to cell area is provided by a logic section 270 which has inputs N (the number of biconcave cells), S8 (the summation of (F1)2 foreach biconcave cell), and S6 (the summation of F1 for each biconcave cell). The variance of the distribution with respect to hemoglobin content is provided by a logic section 272 which has inputs N1, S9 (the summation of 60 (F5)2). and S7 (the summation of (M). A logic section 274 provides the sum K of the output of the logic sections 270 and 272 and a logic section 276 provides the product A of the output of the 10giG sections 270 and 272.

The covariance of the distribution with respect to both the cell area and the cell hemoglobin content is provided by a logic section 278 having inputs N1, S7r S6, and S10 (the summation of the product F1 times F5 65 16 GB 2 093 586 A 16 for each biconcave cell). A logic section 280 squares the output of the logic section 278 to produce an output B. A logic section 282 subtracts the output A of the logic section 276 from the output B of the logic section 280 to provide an output D. K and D are coefficients of a quadratic equation wherein a logic section 282 produces the first solution, EV1, to the quadratic equation, and the logic section 284 produces the second solution, EV2, to the equation.

A logic section 286 produces the mean cell hemoglobin parameter for the biconcave cells by dividing the total hemoglobin content S7 for all the biconcave cells by the number (N1) of the biconcave cells. The mean cell area (MCA) of the biconcave cells is produced by a logic section 288 which divides the total cell area (S6) of the biconcave cells by the total number (N1) of the biconcave cells.

In a similar manner, as shown in Figure 12b, the mean cell area and mean cell hemoglobin parameters are 10 computed for the remaining four classes or subpopulations, i.e., the spherocytes, elongated, irregular, and target cells by eight logic sections 290-297. The number of cells in each subpopulation, N1 -N5, are each transformed into a percentage of the total number of cells by five logic subsections 300-304, in Figure 12b. For example, the percentage of biconcave cells (NC1) is provided by logic subsection 300 which divides the number of biconcave cells (N1) by a total number of cells located by the image processing means (N) and multiplies by 100.

Finally, in the preferred embodiment, two other parameters are calulated which describe the entire population of cells analyzed as illustrated in Figures 12d and 12e. First, amean cell area parameter (MCA) is calculated as a weighted average by multiplying the percentage of a subpopulation (i.e., Ni-NC5 being first divided by 100) bythe mean cell area forthat subpopulation for each subpopulation and adding the products 20 to produce the weighted average. For example, the percentage of biconcave cells (NC1) is multiplied by the mean cell area (MCA1) for the biconcave subpopulation by means of a logic section 306 and the percentage of the spherocyte cells (NC2) is multiplied by the mean cell area of the spherocyte cells (MCA2) by means of a logic section 308 and so on for the other subpopulations and adding these five products by means of a summation logic section 310 to produce the mean cell area (MCA) for the entire population. A weighted average of the hemoglobin content for the entire population (MCH) is produced in a similar manner by a plurality of "multiply" logic sections 312-316 and a summation logic section 318.

In the above manner, 24 parameters descriptive of the various subpopulations of red blood cells and the entire population of red cells as a whole may be calculated, 22 of which are listed in Table 1. They are the percentage of the entire population for each subpopulation, the mean cell area (MCA), and the mean cell hemoglobin (MCH) for each subpopulation, the MCA and MCH for the entire population, the mean central pallor volume (PAL)of the distribution of biconcave and spherocyte cells, the standard deviation (PSD) of the central pallor volume distribution, and the skewness (SKW) of the circularity distribution of the entire population. Two parameters, the mean of the elongation distribution (ELN) and the standard deviation of the elongation distribution (ESD), are calulated but in the preferred embodiment are not reported, as in Table 1.

The parameters in Table I show values calculated for a sample of blood taken from a patient. Similarly, a sample of blood may be taken from another that is known to exhibit one of the known categories of anemia such as iron deficiency, for example, and the 16 parameters may be calculated for the known anemic sample.

Subsequently, the parameters calculated from the analysis of the sample of blood taken from the patient may be compared to the parameters of an iron deficient anemic sample to determine if the sample from the 40 patient resembles iron deficient blood. Likewise, parameters may be calculated for a plurality of known anemic samples of blood wherein the parameters of the patient's sample may be compared and in this mannerthe patient's blood may be classified with respect to those recognized categories of anemic blood.

Referring to Table 1, it is seen that the sample, from which the parameters of Table I were calculated, has been compared to eight types of blood which are normal, iron deficient, chronic disease, B-thalassemia, megaloblastic, hemoglobin ss, hemoglobin SC, and spherocytic.

A specified classification technique to produce a similarity measure for the sample blood taken from the patient which is compared to recognized categories of anemia and normal blood is shown in Figure 14.

Sixteen of the 24 parameters can be thought of as defining a sixteen variable space or sixteen-space. Values for the sixteen different prameters would define a vector having sixteen components, one for each parameter. Thus, when a sample of blood taken from a patient is analyzed, the sixteen parameters calculated therefrom would define a vector Y having sixteen components (Yj.... Yj.... Y16)- Similarly, analysis of samples from the eight previously mentioned types of blood, normal and anemic, would define eight vectors, Wj,j to Wj,8. Each component of the vector for a category of anemic or normal blood is determined by obtaining blood with prior knowledge of the anemic condition and measuring mean parameters over a 55 plurality of such bloods forthe sixteen prameters for that particular category. The vector Y representing the parameter values calculated for the sample of blood taken from the patient may be compared to the vectors representing the mean of various categories of anemic and normal blood. The vector that the vector Y most closely resembles, i.e., is the closest to in the sixteen-space, would determine the classification of the patient's blood.

The first step in the anemia classification logic of Figure 14 is to normalize each parameter value to produce the sixteen components of the Y vector. Thus a parameter value X, which represents the mean cell area of biconcave cells in the patient's blood sample is normalized by a logic subsection 320 to produce the first component of the Y vector, Y1. The logic subsection 320 in normalizing the parameter value X, subtracts the mean value a, from X, and divides by the standard deviation b, of the distribution of X with respect to the 65 Q 11 4 17 GB 2 093 586 A 17 biconcave cells. The distribution for each of the sixteen parameters has been determined with the mean ai and the standard deviation bi for each of the sixteen parameters as set forth in Table VI below.

TABLE VI

5 ai bi (mean) (standard deviation) parameter 1 46.677 7.103 MCA, = X, 10 2 26.531 5.775 MCH, = X2 3 54.115 28.972 EV1 = X3 15 4 3.594 2.120 EV2 = X4 20.271 3.821 PAL = X5 6 4.646 1.161 PSD = X6 20 7 76.934 18.331 NC1 = X7 8 4.167 8.193 NC2 = X8 25 9 3.038 5.600 NC3 = Xg 7.459 7.954 NC4 = X10 11 8.402 10.828 NC5 = X11 30 12 46.490 6.609 MCA X12 13 26.490 5.548 MCH = X13 35 14 6.893 3.059 ELN = X14 5.949 2.363 ESID = X15 16 12.218 4.879 SKW = X16 40 There are sixteen logic sections represented by a logic section 322 which normalizes the parameter Xqi to produce one of the sixteen components Yj of the vector Y. The vector Y representing the sixteen parameter values for the sample of blood taken from the patient is compared to the eight vectors Wjj to W,,8 representing the parameter values for each of the eight categories of blood to determine the proper 45 classification for the patient's blood.

Accordingly there are provided eight logic sections represented by a logic section 324 having the sixteen components Y1 -Y16 of the vector Y as inputs. In addition, each of these logic sections has the sixteen components of a vector representing the parameter values of one of the categories of blood. For example, the first logic section 326 compares the parameter values for the patient's blood with the parameter values for normal blood. In this connection, the vector Wjj represents the sixteen components of the vector for normal blood. The two vectors Y and Wj,j are compared by the logic section 326 wherein the standard distance formula is used to calculate the distance between the two vectors to produce a distance D1.

Referring to Table 1, it is seen that the patient's blood parameter vector has a distance of 0.9 to the normal blood parameter vector. In a like manner, the patient's parameter vector is compared to the parameter vectors for the seven categories of anemia, as seen in Table 1. Thepatient's blood parameter vector has a distance of 4.2 to the iron deficient category of anemia, a distance of 2. 5 to the chronic disease category, and so on.

The sixteen components of the parameter vector for each of the eight categories of anemic and normal blood are set forth in Table VII below.

18 GB 2 093 586 A 18 TABLE V11

Thal, Mega.

4 5 6 7 -0.518 1.675 -0.936 1.664 -0.599 0.908 -0.427 1.067 0.731 0.004 -0.144 0.305 0.330 -0.120 -0.303 -0.041 -0.093 -0.338 0.130 -0.509 -0.377 0.782 -0.547 1.806 -0.967 1.688 0.137 -0.281 0.203 -0.392 0.204 -0.451 Referring backto Figure 3, upon completion of the anemia classification (operation 100),themaster 40 control logicproceedsto print the results (operation 102) of the analysis and similarity comparison or classification. One example of a printout by the preferred method and apparatus has already been given as Table 1. The printout in Table 1 indicated that the sample of bloodanalyzed is closest to normal basedonthe features analyzed. Two more examples are given in Tables Vill and IX respectively, with Table Vill indicating hemoglobin SS anemia and Table IX indication P-thalassemia. Two examples of results of the present red blood cell analysis will be listed below in Tables Vill and IX Normal Iron def 2 0.029 -0.594 0.580 -1.037 -0.639 -0.402 -0.835 0.277 -0.009 0.690 -1.317 0.462 1.073 0.356 -0.439 -0.319 -0.483 -0.201 -0.715 0.153 -0.709 -0.370 0.024 -0.652 0.590 -1.039 -0.607 0.082 -0.767 0.080 -0.903 0.171 Chronic 3 -0.198 -0.061 -0.544 -0.237 -0.285 -0.321 0.196 0.143 -0.411 0.154 -0.340 -0.253 -0.072 -0.476 -0.494 -0.360 ss Spherocytic 8 0.588 -0.715 0.304 -0.577 10 0.854 -0.315 0.192 -0.280 0.228 -2.060 0.921 -0.556 -1.065 -0.642 20 -0.162 2.677 -0.336 -0.439 -0.144 -0.033 2.204 -0.687 0.726 -0.679 30 0.272 -0.449 -0.599 -0.750 -0.270 -0.695 -0.013 -0.497 SC 1 2 3 4 6 7 8 9 11 12 13 14 16 i = 0.761 0.442 1.394 1.606 0.125 1.095 -1.028 -0.158 1.297 0.774 0.620 0.556 0.287 1.162 1.610 1.333 Q z Z Z 19 TABLE VIII

GB 2 093 586 A 19 66.51Y. Biconcave MCA MCH MCA 51 2.3% Spherocytes 48 30 5 MCH 29 7.4% Elongated 41 24 EV1 72 14.6% Irregular 42 25 10 EV2 6 9.2% Targets 58 32 Average 49MCA 29MCH18PAL 5PSD 17SKW 5.2 Normal. 4.2 Megaloblastic 15 3.8 Iron Deficient 1.6 Hemoglobin SS 4.0 Disease 3.8 Hemoglobin SC 20 3.9 P-Thalassemia 5.9 Spherocytic TABLE IX

78.2% Biconcave MCA MCH 25 MCA 37 2.4% Spherocytes 31 19 MCH 21 1.0% Elongated 26 15 EVII 36 13.5% Irregular 33 19 EV2 3 4.9% Ta rg ets 39 21 Average 36 MCA 21 MCH17 PAL 5 PSD 12 SKW 4.2 Normal 1.7 Iron Deficient 2.5 Chronic Disease 1.8 P-Thalassemia 6.5 Megaloblastic 5.7 Hemoglobin SS 4.9 Hemoglobin SC 4.1 Spherocytic From the foregoing, it will be seen that a new and improved analysis of red blood cells heretofore not possible is provided, and the discovery that the red blood cell population carries sufficient information to diagnose many anemias without resort to other conventional tests. That is, subtle, slight, and early changes in either cell production or destruction may now be discovered such as the incipiency of an anemia because of the ability to measure accurately the hemoglobin content in individual blood cells, the cell shape, the cell size variations in size of central pallor, as well as the count of individual cell sizes, and an understanding of the total red cell population in blood samples that was heretofore not possible. With the continual process of red blood cell production and destruction over a 120day life span for each cell, the high percentage of old normal red blood would mask the smaller number of new red cells being produced at the incipiency of a particular anemia. For instance, a chronic disease anemia such as caused by an infection, cancer, or tuberculosis, may cause the new red blood cells being produced to be smaller in size with larger central pallors than is normal and with decreased hemoglobin content. Naturally, early detection of a chronic disease kind of anemia would be most helpful in the treatment of this particular anemia.

For the first time, a total spectrum or galaxy of cells (i.e., the red cell population) may be autornatially analyzed on an individual cell basis and over a sufficient number of cells and with sufficient accuracy to detect a dispersion of distribution indicative of a particular kind of anemia. Take, for example, a blood sample from a person suffering from an iron deficiency anemia. Typically, in such an instance, the usual values of dispersion, e.g. EV1, will be increased from a normal blood value of about 30 to values such as 45 or higher, and likewise the usual value of dispersion for EV2 will be increased from about 2 for a normal blood to about 3 or 4 or higher for an anemia such as an iron deficiency anemia. These slightly higher values of dispersion of distribution indicate that the normal cell population is changed because of the addition of GB 2 093 586 A these additional cells having small size, a large central pallor, and low amounts of hemoglobin. Thus, the diagnostician will see thatthe normal closely packed cell population has been expanded by these new cells formed afterthe onset of the anemia. The other values, such as location of the central tendency of the population's dispersion of distribution, which location is defined by the MVA and MCH values, may also have shifted because of the new anemia cells present. It is to be understood that the parameters used herein such as MCA, MCH, EV1, EV2, skewness, etc. , along with the cenral pallor descriptors, have been experimentally found to be most powerful (atthis time) in analysis of the anemias illustrated and described herein. With further investigation, it may be that other parameters and/or measures may be used to describe and define the red cell population. Thus, by way of example, the covariance, or correlation coefficients could be used to describe the red cell population as well as other measures, which are used to define a distribution and dispersion, and which could be used in lieu of the terms used herein in detail to describe the population of red blood cells.

Further, the present description refers principally to the variables of mean cell hemoglobin and mean cell size although other variables such as pallor size and mean cell area have been tried and could'be used. The particular prameters used and the names thereof my also be changed from that described herein.

Of the seven listed anemias, the hardest anemias to distinguish from each other are the iron deficiency, chronic disease, and P-Thalassemia. The accuracy of the diagnosis of these three anemias from one another is thought to be about 80% accurate with the existing equipment described herein. It is thought that the other anemias listed herein can be diagnosed with almost 100% accuracy. Generally speaking, when it has been found that the closest anemia listed was one of the three anemias of iron deficiency, chronic disease, B-Thalassemia and that the closest anemia was not verified, the actual verified anemia will then be the next closest one listed of these three anemias. For this reason and other reasons, it is preferred to quantify the closeness of several anemias so that if the first anemia is not verified then the second closer anemia can be next chosen and examined for verification. 25 It will be recognized by hematologists skilled in the art that the diagnosis of some anemias, such as iron deficiency or chronic disease anemia, are most difficult today with conventional equipment even with all the information of other tests available to the hematologist. The above apparatus and method should provide a very useful too[ for verification of a particular anemia when the other conventional tests have been used and need to be verified. 30 Although the term "anemia" has been used extensively in this description, it should be noted that the term 30 has been used in the general sense, and the present techniques may be used in a detection of other red cell disorders or pathologies such as hereditary elliptocytosis, for example, or others. Although there has been described herein the use,of first and second microcomputers, it is to be understood that only one larger computer could be used, or hard wired logic could be used. On the other hand, more than one additional microcomputer may be added with each simultaneously measuring characteristics of a different red bloodcell and each classifying different red blood cells into subpopulations. Thus it is considered that one or more additional microcomputer may be used than described hereinto expedite the system.

A specific example of the preferred equipment practicing the techniques described is as follows: in the preferred embodiment, the master control logic 28 and the image processing logic 22 which carry out the flow diagram of Figure 3 comprise two microcomputers such as the Digital Equipernt Corporation LSI/1 1 microprocessors.

Claims

1. An apparatus for automatically analyzing red blood cells in a sample of a patient's blood, comprising:

means for examining the red blood cells in a patient's blood sample; means for measuring characteristics of the red blood cells; means for classifying the cells into a plurality of mutually exclusive subpopulations including a normal and abnormal subpopulation; means for determining subpopulations parameters including the dispersion of the distribution of at least one of said subpopulations parameters; and means for reporting a description of the cells based on said dispersion.

2. Apparatus as claimed in claim 1 and substantially as set forth hereinbefore.

3. Apparatus as claimed in claim 1 and substantially as set forth hereinbefore with reference to the 55 accompanying drawings and as illustrated in Figures 1 and 2 of those drawings.

New claims or amendments to claims filed on 5 May 1982 Superseded claim 1 60 Amended claim 1 i Z 21 GB 2 093 586 A 21 CLAIMS 1. An apparatus for automatically analyzing red blood cells in a sample of a patient's blood, comprising: means for examining the red blood cells in a patient's blood sample; means for measuring characteristics having a range of different values for the red blood cells; means for classifying the cells into a plurality of mutually exclusive subpopulations including a normal and abnormal subpopulation; means for determining sub-population parameters with respect to measured characteristics including the dispersion of the distribution of the measured characterictics of at least one of said subpopulations; and means for reporting a description of the cells based on said dispersion.

Printed for Her Majesty's Stationery Office, by Croydon Printing Company Limited, Croydon, Surrey, 1982. Published by The Patent Office, 25 Southampton Buildings, London, WC2A lAY, from which copies may be obtained.

i