WO1981003080A1 - Method and apparatus for measuring mean cell volume of red blood cells - Google Patents

Method and apparatus for measuring mean cell volume of red blood cells Download PDF

Info

Publication number
WO1981003080A1
WO1981003080A1 PCT/US1981/000510 US8100510W WO8103080A1 WO 1981003080 A1 WO1981003080 A1 WO 1981003080A1 US 8100510 W US8100510 W US 8100510W WO 8103080 A1 WO8103080 A1 WO 8103080A1
Authority
WO
WIPO (PCT)
Prior art keywords
cells
cell
central
size
representative
Prior art date
Application number
PCT/US1981/000510
Other languages
French (fr)
Inventor
J Bacus
Original Assignee
Rush Presbyterian St Luke
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Rush Presbyterian St Luke filed Critical Rush Presbyterian St Luke
Publication of WO1981003080A1 publication Critical patent/WO1981003080A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/64Analysis of geometric attributes of convexity or concavity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/42Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/698Matching; Classification
    • G01N2015/012
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30024Cell structures in vitro; Tissue sections in vitro

Definitions

  • This invention relates to a method of and an apparatus for measuring and reporting size information about red blood cells and particularly the mean cell volume of red blood cells in a blood specimen.
  • the Coulter blood cell ⁇ analyzing equipment is calibrated using spheres of a known size.
  • the blood cells are flattened and thin and many red blood cells contain thin central areas called central pallors which reduce substantially the volume of the cell from the volume the cell would have if it did not have a central pallor.
  • the change in size of central pallors of red blood cells appears to be a good indicator of changes in the blood disorders, as discussed in the first mentioned copending patent application.
  • Hydrodynamic forces used in focusing the blood cells and passing the cells through the aperture of the Coulter cell sizing apparatus change the shape of the thin flexible cells from the shape the cells have in their natural relaxed state.
  • the Coulter counter is provided with a shape factor, so that the actual measured signal from the measuring aperture is multiplied by this factor to obtain the final value of mean cell volume. It is thought that this shape factor is about 1.4 for today's conventional Coulter equipment. Another reason for preferring to generate mean cell volume data for red blood cells rather than mean cell area, is a better segregation of the data into more meaningful and more discrete patterns for blood order diagnosis. More specifically, cells may be substantially similar in area and differ markedly in volume because of central pallor size or a lack of central pallor.
  • normal blood cells and iron deficiency cells each typically have sizes clustered in about the 40 to 60 square micron range and it is not possible to distinguish such cells from one another on the basis of area.
  • these same normal cells should have volumes clustered in range of about 75 to 100 cubic microns which is substantially different than the range of about 50 to 75 cubic microns for the iron deficiency cells.
  • FIGURE 1 A graphic illustration of the plot of mean cell area by an image analysis technique versus mean cell volume from a Coulter counter type of sizing apparatus is set forth in FIGURE 1 of an article of "Bentley, S.A. and S.M. Lewis, 'The Use of an Image Analyzing Computer for the Quantification of Red Blood Morphological Characteristics', Brit. J. Haemat.
  • the cells used in this study were selected manually from each patient specimen and were processed with off-line general purpose computer equipment. Unlike the apparatus used in that work, a commercially practical image analysis system must be automatic, and competitive in speed and cost with the Coulter system in obtaining information on cells.
  • the Bentley and Lewis technique, the introbe Indices technique, and the Coulter counter technique all provide size information for the total blood cell count and cannot make and correlate size data for a given kind of cell or for a given subpopulation of abnormal cells.
  • the area differentiations are not as separated as volume size characteristics.
  • _ ery specific data as to volume differences for microcytic cells versus normal cells because of the difference in the thicknesses (the microcytic cells being thinner) as well as the differences in areas at various times from a patient undergoing treatment may provide an insight as to the effectiveness of the treatment at a very early stage.
  • the present invention may also be used to correlate the sizes of cells measured into different classes or categories such as microcytic, normocytic, or macrocytic with the cell size measurements obtained with conventional techniques.
  • cell sizes have been measured and classified into these classes by image analysis techniques but the results have been poor and did not correlate sufficiently with the results from conventional equipment.
  • a general object of the invention -i-s to provide a new and improved method and apparatus for measuring the size of red blood cells.
  • Another object of the invention is to provide a method of and an apparatus in using a central pallor analysis on the red blood cells having central pallor as part of the determination of the mean cell volume for a specimen of cells.
  • a further object of the invention is to determine the mean cell volume of a total population or of a particular subpopulation of red blood cells.
  • a further object of the invention is to provide an image analysis of red blood cells and to derive therefrom a size classification or size data which is correlated with data derived with conventional equipment.
  • FIGURE 1 is a block diagram of the preferred embodiment of the invention.
  • FIGURE 2 is a graph illustrating the mean cell- volume data as generated by the present invention and by conventional equipment
  • FIGURE 2A is a graph illustrating a comparison - of mean cell area and mean cell volume for the same blood specimens indicated in FIGURE 2A.
  • FIGURE 3 is a perspective view of an apparatus for practicing the method and embodying novel features of the invention.
  • FIGURE 4 is a block diagram showing the operation of the apparatus illustrated in FIGURE 3;
  • FIGURE 5 is a block diagram of the preferred process for analyzing and classifying blood cells and for determining mean cell of volumes;
  • FIGURE 6 illustrates a scanning technique for locating cells and determining the boundary points of cells in an image
  • FIGURES 7a, 7b, and 7c are flow charts of the preferred techniques for classifying the blood cells into mutually exclusive subpopulations
  • FIGURE 8 is a diagrammatic view of a model for red blood cell central pallor measurement.
  • FIGURE 10 is a block diagram of the preferred process for determining whether a cell is round
  • FIGURES 11a, lib, and lie are graphs illustrating thickness/density profile measurements for three different, 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 lie illustrating a "flat" cell having little or no central pallor development;
  • FIGURES 12a, 12b and 12c are graphs illustrating the profiles of the cells of Figures 11a, lib, and lie with the peaks and valleys of each profile labelled;
  • FIGURE 13 is a schematic of the preferred process for accumulating red blood cell subpopulations parameters
  • FIGURES 14a, 14b, 14c, 14d and 14e are schematics illustrating the preferred process of computing the subpopulation characteristics from the accumulated values from a plurality of cells.
  • FIGURE 15 is a logic section for generating a mean cell volume report.
  • the apparatus 10 comprises a microscopic digital image processing and pattern recognition system which analyzes a monolayer of red blood cells on a microscope 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 digitized 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 analog to digital-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 electronic 26 which are in 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.
  • 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 parameter 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 3 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 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 ca era electronics are housed in a section 49 below the monitor.
  • the next lower section 51 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 examined 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 classified into mutually exclusive subpopulations and reported out so that the presence of a minor number of abnormal cells is not overlooked or forgotten and so that accurate parameters about a given subpopulation may also be provided. 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.
  • the individual red blood cells may be analyzed and classified with less subjectivity into a large number of mutually exclusive subpopulations such as biconcave (round cells with central pallor), elongated cells, targets, and irregular cells (cells not fitting into any of the above classifications) .
  • mutually exclusive subpopulations 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 I of the aforesaid patent application.
  • MCH mean cell hemoglobin
  • the individual cells are counted for each subpopulation to provide their respective G percentages of the total population; and likew ⁇ se mean cell volume (MCV) for each subpopulation may also be reported out in a format such as shown in Table I in the copending application. It has been found to be helpful,. in detecting abnormalities 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.
  • red cell size was measured as the projected area of the red cell in square microns. This is a two-dimensional description of size and does not contain any volumetric information regarding size, such as the thickness of the cell, or a decrease in the volume of the cell due to increased central pallor. It was not evident prior to the invention that a projected area measurement of "size" in square microns, was not equivalent in a diagnostic sense to a volume measurement of "size” in cubic microns, such as that which could be obtained with a Coulter counter or the like, or by determining the hematocrit value of the blood sample and then dividing by the red cell count.
  • FIGURE 2A is compared cell size measurements from the blood of persons with iron deficiency anemia, megaloblastic anemia and normal blood. Size measured by an electrical impedance apparatus, in this case the Coulter Model S (MCV) is compared to the area analysis (MCA) from the image analysis equipment disclosed in the aforesaid co-pending application. Notice that if the results are projected on the MCV axis there are three distinct clusters, i.e. a separation of data; whereas, if the results are projected on the MCA axis these distinctions are not as apparent. This indicates that an inclusion of volumetric information is desirable when reporting a measure of red cell size.
  • MCV Coulter Model S
  • MCA area analysis
  • mean cell volumes are generated which take into account the central pallors of the red blood cells to provide data which may be directly correlated with MCV data generated in the past or presently being generated by conventional equipment using conventional Coulter counter equipment.
  • This is achieved by using central pallor data or central pallor signals in combination with the area and hemoglobin characteristic data or signals to generate an output representative of the mean cell volume which takes into consideration the actual volumes of the individual central pallors (if any) of the cells being measured.
  • representative signals of mean cell area (MCA) , mean cell hemoglobin (MCH) , and mean cell pallor (PAL) are generated with the above described apparatus and are sent to a means which generates an output representative of the mean cell volumes for the blood cells.
  • the accuracy of the present invention in providing mean cell volume as related to similar measurements from a Coulter counter instrument is readily apparent from a consideration of FIGURE 2.
  • the data from the image analysis is substantially similar on the ordinate for the volume on the abscissa for the Coulter counter measurements of the same blood samples.
  • Kl, K2, K and K4 are used in connection with the measured values of MCA, MCH and PAL, with the parameters Kl, K2 and K3 each being a multiplier for these measured values as indicated in FIGURE 15.
  • the values for Kl, K2 and K3 have been determined experimentally as will be explained hereinafter.
  • a fourth factor K4 is added to the sum of MCA(Kl) + MCH(K2) + PAL(-K3) and is an offset factor indicating the amount of offset from the juncture of the- abscissa and ordinate of a plot of the MCV's. This offset is thought to represent a factor due to drying of the blood cells prior to image analysis with the apparatus 10 disclosed herein.
  • the preferred means for determining mean cell volume comprises either a digital logic system of electrical devices or a programmed microprocessor which uses Boolean logic.
  • the cells are classified into subpopulations related to a specific anemia such as set forth in Table I therein.
  • the subpopulations given in Table I were into he atologically recognized subpopulations such as normocytes with central pallor, normocytes without central pallor, spherocytes, etc.
  • the size of the cells being listed as mean cell area in square microns.
  • the mean cell size may be printed out for the entire population as well as for a given subpopulation.
  • the mean cell volumes given hereinafter may also be given for subpopulations as well as for the total population.
  • the present invention is capable of reporting the total population and the average mean cell hemoglobin as well as the average mean cell ⁇ volume for the entire population, the average mean cell hemoglobin may be reported out in the line with average parameters of Table I of the aforesaid copending application.
  • the invention will be described as having the ability to classify red blood cells into the several mutually exclusive subpopulations set forth in Table I of the aforesaid application.
  • the subpopulations listed therein are preferred subpopulations for classifying blood with respect to recognized categories of anemias but there may be other subpopulations defined. It is also possible to provide a mean cell hemoglobin for a subpopulation of cells, such as biconcave cells, as will be explained hereinafter. Additionally, it is possible with modifications of the analysis logic to determine the cell volume for each cell and then to subsequently determine the mean cell volume for the total population or for any given subpopulation of course, the size classification of the cells may be reported out in other manners such as microcytic, normoc tic, or macrocytic.
  • the present invention is also of utility in correlating other and existing equipment which has not performed adequately in reporting out cell size classifications as microcytic, normocytic, or macrocytic for the reason that the reported classification did not match the results obtained, with conventional equipment.
  • a multiple parallel logic architecture has been found to provide the rapid processing necessary for efficient analyzing of cells on a slide.
  • a first processing means the master control logic 28 (FIGURE 4)
  • a second processing means the image processing logic 22 as shown in FIGURE 4.
  • the analysis of the cells on a slide requires a sequence of operations to be 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 5 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.
  • 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 waiting for the start signal (operation 60) from the master control logic.
  • the image processing logic 22 Upon receipt of the start signal, the image processing logic 22 proceeds to operation 62 which includes digitizing the image produced b the vidicon camera 16 (FIGURE 4) .
  • 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 4 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.
  • the image processing logic After transmitting the "digitizing dpne" 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) .
  • 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 secton 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.
  • the master control logic proceeds to operation 72 wherein the master control logic 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.
  • the image processing logic Upon receipt of the acknowledge signal, the image processing logic proceeds to digitize the image of the 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 a means cell volume (MCU) determination (operation 100) and prints the results (operation 102) , as will also be more fully explained below.
  • MCU means cell volume
  • the apparatus may be used to provide an output of an anemia classification as described in the aforementioned copending applications, or the present invention could be made a "stand alone" unit whose only function would be to provide mean cell size (MCV) for a total specimen without having to do any classifying into subpopulations, or anemia classifications.
  • MCV mean cell size
  • the master control logic is free to -16- control the imaging means wherein a new field 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.
  • the image processing logic may simultaneously digitize and analyze a new image provided by the new field which had been brought into view by the master control logic.
  • the optimization of the time of analysis as well as the number of features used in the classification logic is achieved so that the amount of storage and classifying techniques may be
  • classification logic has been invented for a large number of subpopulations such as those shown in Table I in the aforesaid application.
  • the preferred classification features are size, hemoglobin content, spicularity, roundness, elongation, central peak height
  • 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 wave length of about 415 Nanometers. At this optical wave length, the red 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.
  • 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.
  • 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.
  • the preferred classification means comprise either a digital logic system of electrical devices or a programmed microprocessor which uses Boolean logic to classify the red blood cells.
  • the images of the cells are digitized (operaton 62 of FIGURE 5) 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.
  • FIGURE 6 which illustrates in
  • the raster scan of the digitized image is then continued from the bottom pixel 114f to hit the next digitized cell 112 by impacting a pixel 112a which is above the threshold as seen in block 116.
  • the raster scan continues from the bottom pixel 112b, and, as seen in block 118, no more cells are located in the image field.
  • the image processing logic transmits the cell features and subpopulation classifications to the master control logic (operation 80) as shown in FIGURE 5.
  • the initial image processing done by the image processing logic outlined in FIGURE 5 is shown in _ greater detail in FIGURE 7a.
  • the image is scanned to locate a cell (operation 66) and the boundary is traced as explained above.
  • 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 to FIGURE 9, a digital image of cells as* defined by their boundary pixels 120 are illustrated.
  • features F1-F4 are calculated in an operation 124 by the image processing logic as shown in FIGURE 7a.
  • Feature Fl 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 2 perimeter) /area and is of assistance in classifying round and non-round objects. A round object would have a theoretical value of 4 and non-round objects have greater values.
  • 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 is 14.0, and is a better approximation to this reference number as the number of pixels, or size, of an object increases.
  • Feature F3 and F4 relate to the spicularity and elongated 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 9.
  • 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 cell.
  • Feature F6 which is a measure of the pallor 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 cross-sectional orthogonal 3-peaked thickness/density profiles, either having 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 cross-sectional, orthogonal, 2-peaked thickness/density profiles.
  • Feature F9 is a measure of the degree of roundness of the pallor itself, and is also used in distinguishing target cells.
  • FIGURES 7a, 7b, and 7c 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 7a, 7b, and 7c.
  • the logic decision are made using the various features together with threshold values that are identified as Tl through Til.
  • the thresholds Tl-Tll are described in Table V and specific values are also provided.
  • 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 7a, 7b, and 7c.
  • FIGURES 7a, 7b, and 7c 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.
  • an object that is located is examined by logic section 138 to determine if it is 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 microns , 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.
  • the optics and television electronics should be adjusted such that round objects of the following dimensions produce the given number of pixels.
  • the decision logic then operates to determine whether the cell is round or non-round. This is performed by a logic section indicated generally 140.
  • the logic section 140 is shown in FIGURE 10 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 7a) which is the first step in the target cell analysis and central pallor analysis.
  • logic subsection 150 (FIGURE 7a) 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 nonround cell of the type for cells of classes 3 and 4. Thus, it cannot be accurately classified and it is for this reason that the subsection 150 eliminates such cells from further consideration.
  • 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.
  • 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 object is not round. If feature F2 is not greater than threshold T2, it is one indication that the cell 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 7a). .
  • thickness/density profiles are extracted from the cell image. These profiles are illustrated in FIGURES lla-llc and 12a-12c.
  • 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 11a for a biconcave cell determined in two orthogonal or transverse directions, a and b.
  • FIGURES lib Two profiles each are also illustrated in FIGURES lib and lie for a target cell and a spherocyte cell. As seen in FIGURE lib, 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.
  • the profile has two "peaks”, Pi and P2, and one "valley", VI.
  • Pi 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.
  • VI determines the relative minimum point of thickness density.
  • the target cells have three relative maxima, Pi, P2, and P3, with two relative minima, VI and V2, as shown in FIGURE 12b.
  • the spherocyte has one peak. Pi, and no valleys (FIGURE
  • the image processing logic 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 7b.
  • 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, Pi and P3, or "LEVl".
  • a logic subsection 164 determines whether the two valleys, Via and V2a, are less then LEVl 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 in 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.
  • 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 174 wherein the two valleys, V2a and V2b, are compared to LEV2 which is half the average of the two non-center peaks Plb and 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 compared to a threshold T7 in a logic subsection 178, and if larger, the cell is classified as a target cell (C5) .
  • a threshold T7 in a logic subsection 178
  • the cell is classified as a target cell (C5) .
  • 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 case, both profiles are smoothed to two peaks or less in operations 180 and 183.
  • the logic subsection 174 determines whether the valleys of profile b were not less than LEV2
  • 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 F7 as described.
  • a logic subsection 186 determines whether profile a has only one peak. If so, the variables Pla, P2a, and Via 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
  • a feature F8 which is the average value of the two valleys subtracted from the average value of the four peaks of the two profiles of the cell, is determined by subsection 194. Then the cell feature Fl is examined to determine whether the size of the cell is larger than a threshold T8 at a logic subsection 196.
  • a variable LEV3 is set equal to one-half the value of feature F8 (operation 198) .
  • 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 for the central pallor boundary (operation 202) .
  • the pallor circularity feature F9 is then computed in an operation 204. F9.is calculated as the number of pallor boundary pixels squared divided by the area of the central pallor.
  • F9 is then compared to a threshold 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 cell 208.
  • 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 pallor feature 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 8, where T represents the cell height or thickness, and 132 indicates the indented central pallor region.
  • the cell area is known from previous analysis on that cell, i.e., Fl.
  • feature F5 is the sum of the grey levels for... pixels enclosed by the chain code defining the boundary of the cell.
  • the hemoglobin density is related to the thickness of the cell and in this manner the hemoglobin feature F5 defines a volume of the cell.
  • the cylinder height, or thickness (T) is derived by using the average value of the peaks of the two thickness/density profiles of the cell as:
  • the volume of the central pallor may be calculated as: T times the area of the cell (Fl) minus the hemoglobin content.
  • the percentage pallor volume F6 is:
  • 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 (C3) , an irregular cell (C4) , or a target (C5) .
  • the logic subsection 210 compares the percentage pallor volume feature F6 to a threshold value T10 and the pallor depth feature F8 to a threshold Til and if either feature is less than its associated threshold then the cell is deemed a spherocyte cell (C2) , otherwise it is a biconcave cell (Cl) .
  • 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.
  • FIGURE 13 A diagram illustrating the updating operation in greater detail is shown in FIGURE 13.
  • 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 Fl, F2, F4, F5, and F6 for a particular cell are provided as inputs to the accumulators together with the cell classification C ⁇ to which the cell features pertain. After the measurements for the cell 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.
  • the feature F2 (cell circularity feature) is provided at a line 212 to an accumulator 214.
  • the accumulator 214 produces a running total SI, i.e., accumulates the measurement (F2 - 14.1) 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
  • c:.: - 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.
  • accumulated is the elongation feature F 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.
  • (F4) 220 provides a sum or running total (S3) of the elongation feature F4 squared, 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.
  • the cell classification C- is provided to the * inputs of a logic AND gate 224 and an AND gate 226 with subpopulation Cl constant (i.e., a 1) provided to the other input of the AND gate 224 and subpopulation C2 constant (i.e., a 2) provided to the othe input of AND gate 226.
  • the output of these AND gates are provided to on OR gate 228 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.
  • the accumulator 232 accumulates the sum 2 of the feature (F6) but only when enabled.
  • the gates 224, 226, and 228 permit the accumulators 230 and
  • the output of the accumator 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 central pallor volume for the spherocyte and biconcave cells.
  • a logic AND gate 234 enables accumulators 236, 238, and 240 when C- 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 Fl (cell area) which is provided at Sll, which will be used to calculate the mean cell area parameter for the 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 total of all of the cells' areas for the 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 Nl.
  • additional logic gates permit accumulators to discriminate among the class cells.
  • an AND gate 248 enables accumulators 250, 252, and 254 when the features appearing at the lines
  • the accumulator 250 provides the
  • the accumulated S9 and S10 are used to calculate parameters descriptive of the dispersion, or variation of the bivariate distribution which will be further explained hereinafter.
  • the features for each cell examined by the image processing logic provide the inputs to the logic described in FIGURE 13 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 13 may be used as 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 cell features.
  • 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 for the "digitizing done" signal indicating that the image processing logic has completed digitizing the next field.
  • a parameter for the mean central pallor volume (PAL) of the biconcave and spherocyte cells is provided by a logic subsection 264 having inputs Nl (the number of biconcave cells) , N2 (the number of spherocyte cells) , 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 logic subsection 268 which takes the square root of the output provided by the logic subsection 256 to finally produce the parameter PSD in a manner similar to that of the parameter ESD.
  • a logic diagram is shown for the computation of the parameters EVl and EV2.
  • the general formula for computing the variance of a distribution with respect to 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
  • 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 logic sections 270 and 272.
  • a logic section 278 having inputs Nl, S7, S6, and S10 (the summation of the product Fl times F5 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, EVl, 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 (Nl) 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 (Nl) of the biconcave cells.
  • the mean cell area and mean cell hemoglobin parameters are computed for the remaining four classes or subpopulations, i.e., the spherocytes, elongated, irregular, and target cells by either 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 14b.
  • the percentage of biconcave cells (NCI) is provided by logic subsection 300 which divides the number of biconcave cells (Nl) by a total number of cells located by the image processing means (N) and multiplies by 100.
  • a mean cell area parameter is calculated as a weighted average by multiplying the percentage of a subpopulation (i.e., NC1-NC5 being first divided by 100) by the mean cell area for that subpopulation for each subpopulation and adding the products to produce the weighted average.
  • the percentage of biconcave cells is multiplied by the mean -cell area (MCAl) 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.
  • FIGURE 15 the general method of measuring the mean cell size to provide the close correlation to the MCVs achieved by the Coulter counter, is shown in FIGURE 15.
  • the logic section shown in FIGURE 15 thus will receive as inputs the MCA over line 512 from FIGURE 14d, the MCH over line 511 from FIGURE 14e, and the PAL over line 510 from FIGURE 14c. More specifically, the output for the pallor volume PAL may be applied from FIGURE 14c as input to a multiplier logic section 400 (FIGURE 15) which also receives an input factor -K3 to provide a calculation in form of an output on line 401 leading to the accumulator 403. Since the pallor volume is a negative value it will be subtracted in the accumulator 403.
  • the mean cell areas MCA for the entire population from logic section 310 in FIGURE 14d is applied over line 512 as input to a multiplier logic section 405 (FIGURE 15) along with the input factor K2 to provide an output of their product to the accumulator 403.
  • the mean cell hemoglobin (MCH) provided over line 511 from summation logic 318 (FIGURE 14e) serves an input ot a multiplier logic section 411 along with a factor K2 to provide an output to the summation logic 403.
  • the output of the summation logic 403 is applied over a line 413 to an adder 415 to which is also applied an input factor K4 which is the offset factor.
  • the output from the adder logic 415 is the mean cell volume for the total population of cells.
  • the mean cell volume is reported at 430 by printing it out on a form or by displaying it on a cathode ray tube.
  • mean cell volume (MCV) for a subpopulation of cells can be computed for any given subpopulation of cells.
  • the biconcave cells are separately classified in the example given herein and in the aforesaid application. More specifically, the mean cell area (MCA 1) for the biconcave cells is provided by logic section 288 in FIGURE 14a and this may be applied as an input on line 512 of FIGURE 15. Likewise, the mean cell hemoglobin (MCH 1) for the biconcave cells is provided from logic section 286 (FIGURE 14a) and this may be applied as an input on line 511 of FIGURE 15.
  • the mean central pallor volume (PAL) for the biconcave cells is available from logic section 264 in FIGURE 14c and this may be applied as input over line 510 to the multiplier logic section 400 in the logic shown in FIGURE 15.
  • The. respective constants will be applied as Kl, K2 and K3 to the respective multiplier logic sections 405, 411 and 400.
  • the remaining operation of the logic section shown in FIGURE 15 will be as above described above to provide an output which is the mean cell volume (MCV) for the biconcave subpopulation of cells.
  • the logic section shown in FIGURE 14d may also be used with various other subpopulations of cells to provide mean cell volume data for various abnormal cell populations. Because such data has not been heretofore available, there is a whole new data base for diagnosis of cell disorders, blood diseases or of cell morphological changes. For such abnormal cell populations it is preferred to generate for each cell its cell volume, its hemoglobin, and its pallor measurement and then after the cells are classified to accumulate in logic sections these values and then determine the mean cell area (MCA) , mean cell hemoglobin (MCH) , and mean central pallor volume (PAL) . After this, a logic section, such as shown in FIGURE 15 is used to provide a mean cell volume report or output for a given abnormal cell population.
  • MCA mean cell area
  • MCH mean cell hemoglobin
  • PAL mean central pallor volume
  • the constants -Kl, K2, K3 and K4 are exemplary of the constants which may be used and these constants have been derived for red blood cells which have been spread by a spinner and then dried before image analysis, as above described.
  • the factors tested herein automatically take into account the cell distortion due to drying of the cells prior to examination.
  • the constants would be different in order to adjust the mean cell volumes to that achieved by another piece of equipment such as a Coulter Model S counter.
  • the factors Kl, K2, K3, and K4 were obtained by using a standard multiple linear regression technique as described fully in a publication entitled Nunsrical Method For Engineers and Engineers by R.W. Hamming, published by McGraw Hill in 1962, and by using the mean cell volume data for the same blood specimens as had been previously measured with a Coulter Model S counter.
  • Other techniques such as comparison of data obtained emperically, may be used to develop a correlation between the image analysis mean cell volumes and the mean cell volumes obtained with a conventional mean cell volume measuring technique such as the Wintrobe technique or the electrical impedance technique of Coulter.
  • the output from logic section 415 may be used to classify cells as microcytic, normocytic or macrocytic by having the output of logic section 415 applied to a logic section (not shown) having three levels with MCV below a given level being classified as microcytic, with the MCV's in a central range as normocytic, and with the* MCV's above an upperlevel of the normocytic range being classified as macrocytic.
  • the size or volume information already being generated by such equipment may be correlated to that of conventional equipment producing MCV data by using the techniques herein disclosed.
  • the present invention is not to be construed as being limited to the equipment herein described or to equipment that provides an output only in the terms of mean cell volume (MCV) .
  • the present invention provides a new and improved method and apparatus for generating cell size information correlated with cell size information such as mean cell volumes generated with conventional equipment.
  • cell size information may be reported out as placing the cells in a given size category such as microcytic, normocytic, or macrocytic, it preferably provides a mean cell volume output. Additionally, the present invention may be used to provide such size information for a subpopulation of red blood cells.

Abstract

Apparatus and method for producing signals representative of the mean cell volume of red blood cells in a blood specimen. The apparatus includes means (310) for generating signals representative of the area of the cells and means (318) for measuring the optical density of the individual cells and for generating signals representative of the hemoglobin content or mass of the cells. The central pallor is measured for cells having central pallors and means (264) generate a signal representative of the central pallor of these cells. The representative area signals, hemoglobin signals and central pallor signals are sent to a means (403, 415) which generates an output (430) representative of the mean cell volume of the cells.

Description

-1-
METHOD AND APPARATUS FOR MEASURING MEAN" CELL VOLUME OF RED BLOOD CELLS
This invention relates to a method of and an apparatus for measuring and reporting size information about red blood cells and particularly the mean cell volume of red blood cells in a blood specimen.
In the prior applications referred to above, there is described various apparatus for automatically analyzing blood and providing representative output data of the mean cell size, mean cell hemoglobin, and mean cell density. The mean cell size information generated
2 and reported was expressed as area (u ") rather than
3 volume (u ) the latter being the conventional size information generated with conventional measuring techniques such as the Coulter kind of particle counting apparatus. Since physicians are most familiar with mean cell volume data than with mean cell area data, there is a desire to generate a mean cell volume output that can aid physicians and may also aid in automatic diagnosis of anemia or other blood disorders.
Some work has been done with image analysis and pattern recognition techniques to compare mean cell areas with mean cell volumes. Such equipment has not been very accurate in correlating with mean cell volume data generated by conventional Coulter particle sizing equipment used for blood analysis. The fault for this descrepancy may not be entirely with the image analysis equipment. As is known and has been reported in the literature, the Coulter blood counting equipment suffers from several shortcomings such as producing signals which are in error of a true blood cell volume for cells which are (1) tumbling as they pass through the measuring aperture, (2) passing very close to the aperture wall, (3) in reality a pair of closely adjacent cells in the aperture rather than a single cell, (4) being measured while a previously measured cell is closely adjacent the exit end of the aperture and is generating a ghost signal, etc.
The Coulter blood cell~ analyzing equipment is calibrated using spheres of a known size. However, the blood cells are flattened and thin and many red blood cells contain thin central areas called central pallors which reduce substantially the volume of the cell from the volume the cell would have if it did not have a central pallor. The change in size of central pallors of red blood cells appears to be a good indicator of changes in the blood disorders, as discussed in the first mentioned copending patent application. Hydrodynamic forces used in focusing the blood cells and passing the cells through the aperture of the Coulter cell sizing apparatus change the shape of the thin flexible cells from the shape the cells have in their natural relaxed state. To compensate for various ones of these factors, the Coulter counter is provided with a shape factor, so that the actual measured signal from the measuring aperture is multiplied by this factor to obtain the final value of mean cell volume. It is thought that this shape factor is about 1.4 for today's conventional Coulter equipment. Another reason for preferring to generate mean cell volume data for red blood cells rather than mean cell area, is a better segregation of the data into more meaningful and more discrete patterns for blood order diagnosis. More specifically, cells may be substantially similar in area and differ markedly in volume because of central pallor size or a lack of central pallor. For example, normal blood cells and iron deficiency cells each typically have sizes clustered in about the 40 to 60 square micron range and it is not possible to distinguish such cells from one another on the basis of area. However, these same normal cells should have volumes clustered in range of about 75 to 100 cubic microns which is substantially different than the range of about 50 to 75 cubic microns for the iron deficiency cells. A graphic illustration of the plot of mean cell area by an image analysis technique versus mean cell volume from a Coulter counter type of sizing apparatus is set forth in FIGURE 1 of an article of "Bentley, S.A. and S.M. Lewis, 'The Use of an Image Analyzing Computer for the Quantification of Red Blood Morphological Characteristics', Brit. J. Haemat. 29:81, 1975". The cells used in this study were selected manually from each patient specimen and were processed with off-line general purpose computer equipment. Unlike the apparatus used in that work, a commercially practical image analysis system must be automatic, and competitive in speed and cost with the Coulter system in obtaining information on cells.
The Bentley and Lewis technique, the introbe Indices technique, and the Coulter counter technique all provide size information for the total blood cell count and cannot make and correlate size data for a given kind of cell or for a given subpopulation of abnormal cells. With the equipment described in the copending applications, it is possible to segregate and to measure the size of individual cells as well as the mean size for various abnormal cell subpopulations such as spherocytes, macrocytes, microcytes, etc. However, as above explained, the area differentiations are not as separated as volume size characteristics. By providing mean cell volume for subpopulations of abnormal cells on a large scale basis, new insights should be gained into the cell volume characteristics of a given blood disorder and should lead to more substantial and more accurate information on the volume relationship of abnormal cells to normal cells. For example,_ ery specific data as to volume differences for microcytic cells versus normal cells because of the difference in the thicknesses (the microcytic cells being thinner) as well as the differences in areas at various times from a patient undergoing treatment may provide an insight as to the effectiveness of the treatment at a very early stage.
The present invention may also be used to correlate the sizes of cells measured into different classes or categories such as microcytic, normocytic, or macrocytic with the cell size measurements obtained with conventional techniques. Heretofore, cell sizes have been measured and classified into these classes by image analysis techniques but the results have been poor and did not correlate sufficiently with the results from conventional equipment.
Accordingly, a general object of the invention -i-s to provide a new and improved method and apparatus for measuring the size of red blood cells.
Another object of the invention is to provide a method of and an apparatus in using a central pallor analysis on the red blood cells having central pallor as part of the determination of the mean cell volume for a specimen of cells.
A further object of the invention is to determine the mean cell volume of a total population or of a particular subpopulation of red blood cells.
A further object of the invention is to provide an image analysis of red blood cells and to derive therefrom a size classification or size data which is correlated with data derived with conventional equipment.
These and other objects and advantages of the present invention will become apparent from the following detailed description taken in connection with
c. the accompanied drawings in which:
FIGURE 1 is a block diagram of the preferred embodiment of the invention;
FIGURE 2 is a graph illustrating the mean cell- volume data as generated by the present invention and by conventional equipment;
FIGURE 2A is a graph illustrating a comparison - of mean cell area and mean cell volume for the same blood specimens indicated in FIGURE 2A. FIGURE 3 is a perspective view of an apparatus for practicing the method and embodying novel features of the invention;
FIGURE 4 is a block diagram showing the operation of the apparatus illustrated in FIGURE 3; FIGURE 5 is a block diagram of the preferred process for analyzing and classifying blood cells and for determining mean cell of volumes;
FIGURE 6 illustrates a scanning technique for locating cells and determining the boundary points of cells in an image;
FIGURES 7a, 7b, and 7c are flow charts of the preferred techniques for classifying the blood cells into mutually exclusive subpopulations;
FIGURE 8 is a diagrammatic view of a model for red blood cell central pallor measurement.
FIGURE 9 illustrates a chain code description and analysis method for three diagrammatic red blood cell outlines;
FIGURE 10 is a block diagram of the preferred process for determining whether a cell is round;
FIGURES 11a, lib, and lie are graphs illustrating thickness/density profile measurements for three different, 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 lie illustrating a "flat" cell having little or no central pallor development;
FIGURES 12a, 12b and 12c are graphs illustrating the profiles of the cells of Figures 11a, lib, and lie with the peaks and valleys of each profile labelled;
FIGURE 13 is a schematic of the preferred process for accumulating red blood cell subpopulations parameters;
FIGURES 14a, 14b, 14c, 14d and 14e are schematics illustrating the preferred process of computing the subpopulation characteristics from the accumulated values from a plurality of cells; and
FIGURE 15 is a logic section for generating a mean cell volume report.
As shown in the drawings for purposes of illustration, the invention is embodied in an apparatus 10, such as disclosed in U.S. Patent No. 4,199,748 which is incorporated by reference as if fully reproduced herein. In this equipment, as described fully in U.S.
Patent No. 4,199,748, and as shown in FIGURES 3-6 herein the apparatus 10 comprises a microscopic digital image processing and pattern recognition system which analyzes a monolayer of red blood cells on a microscope 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 digitized 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 analog to digital-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 electronic 26 which are in 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 parameter 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 3 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 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 ca era electronics are housed in a section 49 below the monitor. The next lower section 51 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 examined 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 classified into mutually exclusive subpopulations and reported out so that the presence of a minor number of abnormal cells is not overlooked or forgotten and so that accurate parameters about a given subpopulation may also be provided. 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 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 I of the aforesaid patent application. In addition to the hemoglobin parameters, the individual cells are counted for each subpopulation to provide their respective G percentages of the total population; and likewάse mean cell volume (MCV) for each subpopulation may also be reported out in a format such as shown in Table I in the copending application. It has been found to be helpful,. in detecting abnormalities 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.
According to the method described in U.S. Patent No. 4,199,748, red cell size was measured as the projected area of the red cell in square microns. This is a two-dimensional description of size and does not contain any volumetric information regarding size, such as the thickness of the cell, or a decrease in the volume of the cell due to increased central pallor. It was not evident prior to the invention that a projected area measurement of "size" in square microns, was not equivalent in a diagnostic sense to a volume measurement of "size" in cubic microns, such as that which could be obtained with a Coulter counter or the like, or by determining the hematocrit value of the blood sample and then dividing by the red cell count.
Experimentation with normal and several types of anemic blood specimens have recently indicated that the diagnostic information relative to size is better preserved as a volume measurement rather than an area measurement. This can be better understood from FIGURE 2A, in which are compared cell size measurements from the blood of persons with iron deficiency anemia, megaloblastic anemia and normal blood. Size measured by an electrical impedance apparatus, in this case the Coulter Model S (MCV) is compared to the area analysis (MCA) from the image analysis equipment disclosed in the aforesaid co-pending application. Notice that if the results are projected on the MCV axis there are three distinct clusters, i.e. a separation of data; whereas, if the results are projected on the MCA axis these distinctions are not as apparent. This indicates that an inclusion of volumetric information is desirable when reporting a measure of red cell size.
In accordance with the present invention, mean cell volumes are generated which take into account the central pallors of the red blood cells to provide data which may be directly correlated with MCV data generated in the past or presently being generated by conventional equipment using conventional Coulter counter equipment. This is achieved by using central pallor data or central pallor signals in combination with the area and hemoglobin characteristic data or signals to generate an output representative of the mean cell volume which takes into consideration the actual volumes of the individual central pallors (if any) of the cells being measured. As diagramatically illustrated in FIGURE 15, representative signals of mean cell area (MCA) , mean cell hemoglobin (MCH) , and mean cell pallor (PAL) are generated with the above described apparatus and are sent to a means which generates an output representative of the mean cell volumes for the blood cells. The accuracy of the present invention in providing mean cell volume as related to similar measurements from a Coulter counter instrument is readily apparent from a consideration of FIGURE 2. The data from the image analysis is substantially similar on the ordinate for the volume on the abscissa for the Coulter counter measurements of the same blood samples.
To achieve this MCV measurement four parameters Kl, K2, K and K4 are used in connection with the measured values of MCA, MCH and PAL, with the parameters Kl, K2 and K3 each being a multiplier for these measured values as indicated in FIGURE 15. The values for Kl, K2 and K3 have been determined experimentally as will be explained hereinafter. A fourth factor K4 is added to the sum of MCA(Kl) + MCH(K2) + PAL(-K3) and is an offset factor indicating the amount of offset from the juncture of the- abscissa and ordinate of a plot of the MCV's. This offset is thought to represent a factor due to drying of the blood cells prior to image analysis with the apparatus 10 disclosed herein. As will be explained, in the described embodiment of the invention, these values are Kl=.43; K2=1.94; K3=-.84 and K4=27. The preferred means for determining mean cell volume comprises either a digital logic system of electrical devices or a programmed microprocessor which uses Boolean logic. In the analysis given in aforesaid copending application,, the cells are classified into subpopulations related to a specific anemia such as set forth in Table I therein. In U.S. Patent 4,097,845, the subpopulations given in Table I were into he atologically recognized subpopulations such as normocytes with central pallor, normocytes without central pallor, spherocytes, etc. with the size of the cells being listed as mean cell area in square microns. In this same patent several examples were printed out. Manifestly, the mean cell size may be printed out for the entire population as well as for a given subpopulation. Likewise, with, only slight modification of the analysis logic described herein the mean cell volumes given hereinafter may also be given for subpopulations as well as for the total population.
To aid in understanding classification of the cells as well as the measurements used to classify the cells prior to determining the mean cell volumes therefor, some of the description given in the aforesaid U.S. Patent No. 4,199,748 will be repeated. As will be explained, the present invention is capable of reporting the total population and the average mean cell hemoglobin as well as the average mean cell^volume for the entire population, the average mean cell hemoglobin may be reported out in the line with average parameters of Table I of the aforesaid copending application. Thus, as indicated above, herein the invention will be described as having the ability to classify red blood cells into the several mutually exclusive subpopulations set forth in Table I of the aforesaid application. The subpopulations listed therein are preferred subpopulations for classifying blood with respect to recognized categories of anemias but there may be other subpopulations defined. It is also possible to provide a mean cell hemoglobin for a subpopulation of cells, such as biconcave cells, as will be explained hereinafter. Additionally, it is possible with modifications of the analysis logic to determine the cell volume for each cell and then to subsequently determine the mean cell volume for the total population or for any given subpopulation of course, the size classification of the cells may be reported out in other manners such as microcytic, normoc tic, or macrocytic. The present invention is also of utility in correlating other and existing equipment which has not performed adequately in reporting out cell size classifications as microcytic, normocytic, or macrocytic for the reason that the reported classification did not match the results obtained, with conventional equipment. AS disclosed in the aforesaid U.S. Patent No.
4,199,748, a multiple parallel logic architecture has been found to provide the rapid processing necessary for efficient analyzing of cells on a slide. Herein, there is provided a first processing means, the master control logic 28 (FIGURE 4) , and a second processing means, the image processing logic 22 as shown in FIGURE 4. The analysis of the cells on a slide requires a sequence of operations to be 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 5 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 5, 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 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 b the vidicon camera 16 (FIGURE 4) . 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 4, 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 dpne" 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) .
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 secton 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 thes_e___ operations have already been completed then the master control logic proceeds to operation 72 wherein the master control logic 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 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 a means cell volume (MCU) determination (operation 100) and prints the results (operation 102) , as will also be more fully explained below. The apparatus may be used to provide an output of an anemia classification as described in the aforementioned copending applications, or the present invention could be made a "stand alone" unit whose only function would be to provide mean cell size (MCV) for a total specimen without having to do any classifying into subpopulations, or anemia classifications.
Thus, because of the dual processor architecture, the master control logic is free to -16- control the imaging means wherein a new field 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 accumlating 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 had been brought into view by the master control logic. It
10 should be notd 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 different images.
With the present invention, the optimization of the time of analysis as well as the number of features used in the classification logic is achieved so that the amount of storage and classifying techniques may be
20 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 calssification are comprised. Despite this, a relatively foolproof feature set and
25 classification logic has been invented for a large number of subpopulations such as those shown in Table I in the aforesaid application. The preferred classification features are size, hemoglobin content, spicularity, roundness, elongation, central peak height
30 (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
35 elongated cells.
c?: 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 wave length of about 415 Nanometers. At this optical wave length, the red 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. Furthermore, 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. A particular manner of vapor fixing of cells before they dry without staining thereof to prevent the formation of artifacts by distortion of the central pallor is disclosed in U.S. Patent No. 4,209,548 entitled "Method for the Preparation of Blood Samples for Automated Analysis", issued June 24, 1980, and hereby incorporated by reference as fully reproduced herein. Thus, by rapidly preparing the specimens to a monolayer and fixing with a formaldehyde vapor prior to the drying of red blood cells, as disclosed in the aforementioned co-pending patent application, and by not employing a time consuming staining to contrast enhance the cells, as in white blood cell analysis, these specimens may be quickly prepared and analyzed accurately.
The location of the cell image and the identification and feature extraction has been greatly simplified as described below to 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 measurement. Finally, inner central pallor boundaries are determined and features analyzed for more precise target cell identification.
-After having extracted these identifying features, the cells are then categorized by a classification means. The preferred classification means (FIGURES 7a, 7b, and 7c) 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 to the specific features of the illustrated embodiment of the invention, the images of the cells are digitized (operaton 62 of FIGURE 5) 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 6 which illustrates in
- mTVmΣ greater detail the operation 66 (FIGURE 5) by the image processing logic, an original microscopic image which had been digitized is stored as represented by the 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 5) . In this preferred embodiment of the invention, individual cells 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 114a, which is usually the pixel located first, and the one at the "bottom" of the cell, here pixel 114f, 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 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 114f to hit the next digitized cell 112 by impacting a pixel 112a 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 112b, and, as seen in block 118, no more cells are located in the image field. At this time, the image processing logic transmits the cell features and subpopulation classifications to the master control logic (operation 80) as shown in FIGURE 5.
The initial image processing done by the image processing logic outlined in FIGURE 5 is shown in _ greater detail in FIGURE 7a. 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 to FIGURE 9, 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. and J. H. WEENS, 'An Automated Method of Diffential Red Blood Cell Classification with
Application to the Diagnosis of Anemia', Journal of Histochemistry 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 other features for use in the classification of the cells. In this regard, the following features are used herein:
O:: Table II
Feature Description How Determined
Fl Area size Number of pixels enclosed by cell boundary
F2 Shape (circularity) .(Number of perimeter
2 pixels) /area
F3 Shape (spicularity) ' Number of "spicules"
.on boundary
F4 Shape (elongation) Comparison of orthogonal boundary chain code orientations
F5 Grey levels Sum of grey levels as a measure of Cell Hemoglobin
F6 Pallor (volume) The percentage volume of the central pallor
F7 Central peak The height of the central peak of a 3-peaked profile of a cell
F8 Pallor (depth) For a 2-peaked profile, the difference of the valley from the peak heights F9 Pallor (circularity) (Number of pallor
2 boundary pixels) / area of pallor
As indicated above, features F1-F4 are calculated in an operation 124 by the image processing logic as shown in FIGURE 7a. Feature Fl 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 2 perimeter) /area and is of assistance in classifying round and non-round objects. A round object would have a theoretical value of 4 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 is 14.0, and is a better approximation to this reference number as the number of pixels, or size, of an object increases.
Features F3 and F4 relate to the spicularity and elongated 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 9. 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 cell. Feature F6, which is a measure of the pallor 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 cross-sectional orthogonal 3-peaked thickness/density profiles, either having 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 cross-sectional, orthogonal, 2-peaked thickness/density profiles. Feature F9 is a measure of the degree of roundness of the pallor itself, and is also 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 7a, 7b, and 7c. The logic decision are made using the various features together with threshold values that are identified as Tl through Til. The thresholds Tl-Tll 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 7a, 7b, and 7c. In this regard, FIGURES 7a, 7b, and 7c 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 7a, an object that is located is examined by logic section 138 to determine if it is 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 microns , 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 _gf the chain coded cell and then dividing by a conversion factor 1290 or thereabout to 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 adjusted to produce grey levels corresponding to the following optical.density at 418 nanometers:
Table III
Optical Density Grey Level
.134 17
.294 35
.403 52
.505 43
.605 57
'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 2 Pixels
111 1850
93 1550
77 1283 58 967
34 567
23 383
17 283
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 140. the logic section 140 is shown in FIGURE 10 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 7a) which is the first step in the target cell analysis and central pallor analysis. Similarly, if it is determined that the cell is not round, then logic subsection 150 (FIGURE 7a) 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 nonround cell of the type for cells of classes 3 and 4. Thus, it cannot be accurately classified and it is for this 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 object is not round. If feature F2 is not greater than threshold T2, it is one indication that the cell 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 7a). .
In operation 148 thickness/density profiles are extracted from the cell image. These profiles are illustrated in FIGURES lla-llc and 12a-12c. 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 11a for a biconcave cell determined in two orthogonal or transverse directions, a and b. Two profiles each are also illustrated in FIGURES lib and lie for a target cell and a spherocyte cell. As seen in FIGURE lib, 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 11 is discussed more fully in connection with FIGURES
12a-12c. As seen in FIGURE 12a, the profile has two "peaks", Pi and P2, and one "valley", VI. Pi 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. VI determines the relative minimum point of thickness density. Similarly, the target cells have three relative maxima, Pi, P2, and P3, with two relative minima, VI and V2, as shown in FIGURE 12b. The spherocyte has one peak. Pi, and no valleys (FIGURE
C 12c). These profiles are utilized 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 7b. 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, Pi and P3, or "LEVl". A logic subsection 164 determines whether the two valleys, Via and V2a, are less then LEVl 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 in 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 174 wherein the two valleys, V2a and V2b, are compared to LEV2 which is half the average of the two non-center peaks Plb and 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 compared 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 of the 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 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 F7 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 respectively. 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 above, a logic subsection 186 determines whether profile a has only one peak. If so, the variables Pla, P2a, and Via 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 Plb, P2b, and Vlb are set equal to each other in an operation 192.
Continuing with FIGURE 7c therein, a feature F8, which is the average value of the two valleys subtracted from the average value of the four peaks of the two profiles of the cell, is determined by subsection 194. Then the cell feature Fl 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., Fl is greater than T8, it is possible that the cell is a target cell despite the previous target cell analysis and therefore another target cell analysis will be performed beginning in operation 198. 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 for the central pallor boundary (operation 202) . The pallor circularity feature F9 is then computed in an operation 204. F9.is calculated as the number of pallor boundary pixels squared divided by the area of the central pallor. F9 is then compared to a threshold 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 cell 208. 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 pallor feature 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 8, where T represents the cell height or thickness, 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 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 of the cell. The cylinder height, or thickness (T) , is derived by using the average value of the peaks of the two thickness/density profiles of the cell as:
T = Pla = P2a + Plb + P2b *
Thus, the volume of the central pallor may be calculated as: T times the area of the cell (Fl) minus the hemoglobin content. Finally, the percentage pallor volume F6 is:
F6 = (T x Fl - F5) x 100% T x Fl 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 (C3) , an irregular cell (C4) , or a target (C5) . The logic subsection 210 compares the percentage pallor volume feature F6 to a threshold value T10 and the pallor depth feature F8 to a threshold Til and if either feature is less than its associated threshold then the cell is deemed a spherocyte cell (C2) , otherwise it is a biconcave cell (Cl) .
Referring back to FIGURE 5, the feature extraction operation 76 and the cell subpopulation classification operation 78 have been completed for the
c 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 7a, 7b, and 7c 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.
Table V
Threshold Value Description
Tl Size threshold for artifact
T2 54 Size threshold for double cells
T3 25 Elongated threshold
T4 16 Cell cirularity threshold
T5 7 Spiculed threshold
T6 25 Elongation threshold
T7 5 grey levels Target center peak height threshold
T8 47 2 Size threshold for target cells
T9 20 Pallor circularity threshold T10 11% Pallor volume threshold Til 8 grey levels Depth of pallor threshold
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 5. 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 analysed (operation 84) . A diagram illustrating the updating operation in greater detail is shown in FIGURE 13. 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 Fl, F2, F4, F5, and F6 for a particular cell are provided as inputs to the accumulators together with the cell classification C^ to which the cell features pertain. After the measurements for the cell 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 SI, i.e., accumulates the measurement (F2 - 14.1) 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
c:.: - 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 F 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 feature F4 squared, 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.
In the aforesaid U.S. Patent No. 4,199,748, not all feature measurements were accumulated for each subpopulation. For example, in that patent the feature F6 (pallor volume) was only accumulated for the biconcave cells (subpopulation Cl) and the spherocyte cells (subpopulation C2) . Therefore, in addition to the features for a particular cell, the subpopulation classification for the particular cell to which the features pertain was provided as shown as C- at line 222. A plurality of logic utilize the input C- to discriminate among the cell subpopulations. Thus, the cell classification C- is provided to the* inputs of a logic AND gate 224 and an AND gate 226 with subpopulation Cl constant (i.e., a 1) provided to the other input of the AND gate 224 and subpopulation C2 constant (i.e., a 2) provided to the othe input of AND gate 226. The output of these AND gates are provided to on OR gate 228 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 2 of the feature (F6) but only when enabled. Thus, the gates 224, 226, and 228 permit the accumulators 230 and
232 to accumate the measurements derived from the feature F6 only when the feature had been extracted from a Cl or C2 biconcave or spherocyte class cell. The output of the accumator 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 central pallor volume for the spherocyte and biconcave cells. Similarly, a logic AND gate 234 enables accumulators 236, 238, and 240 when C- 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 Fl (cell area) which is provided at Sll, which will be used to calculate the mean cell area parameter for the 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. In a similar manner the total cell area for the elongated (C3) , the irregular (C4) , and target (C5) cells are provided at S13, S15, and S17, respectively. The total of all cells' hemoglobin content for the elongated, irregular, and target cells is provided at S14, S16, and S18, respectively. The total number of cells in each of the above subpopulations is provided at
c:-: N3, N4, and N5.
Likewise, the total of all of the cells' areas for the 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 Nl. For additional accumulated measurements on the biconcave subpopulation, additional logic gates permit accumulators to discriminate among the class cells. Thus, an AND gate 248 enables accumulators 250, 252, and 254 when the features appearing at the lines
244 and 246 have been extracted from a Cl, i.e., a biconcave cell. The accumulator 250 provides the
2 accumulated sum of the measurement (Fl) at S8. The accumulator 252 similarly provides the accumulated total of the measurement <F5) at S9. Finally, the accumulator 254 provides the accumulated sum of the product of the feature Fl times the feature F5 (Fl x F5) . The accumulated S9 and S10 are used to calculate parameters descriptive of the dispersion, or variation 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 13 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 13 may be used as 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 cell features. Referring back to FIGURE 5, 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 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 13 for those N cells are used to calculate the parameters descriptive of the subpopulations (operation 90) which is illustrated in greater detail in FIGURES 14a through 14e.
A parameter for the mean central pallor volume (PAL) of the biconcave and spherocyte cells is provided by a logic subsection 264 having inputs Nl (the number of biconcave cells) , N2 (the number of spherocyte cells) , 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 logic subsection 268 which takes the square root of the output provided by the logic subsection 256 to finally produce the parameter PSD in a manner similar to that of the parameter ESD. Referring to FIGURE 14a, a logic diagram is shown for the computation of the parameters EVl and EV2. The general formula for computing the variance of a distribution with respect to 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
2 biconcave cells) , S8 (the summation of (Fl) for each biconcave cell) , and S6 (the summation of Fl for each biconcave cell) . The variance of the distribution with respect to hemoglobin content is provided by a logic section 272 which has inputs Nl, S9 (the summation of
2 (F5) ), and S7 (the summation of (F5) ) . 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 logic 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 Nl, S7, S6, and S10 (the summation of the product Fl times F5 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, EVl, 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 (Nl) 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 (Nl) of the biconcave cells.
In a similar manner, as shown in FIGURE 14b, the mean cell area and mean cell hemoglobin parameters are computed for the remaining four classes or subpopulations, i.e., the spherocytes, elongated, irregular, and target cells by either 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 14b. For example, the percentage of biconcave cells (NCI) is provided by logic subsection 300 which divides the number of biconcave cells (Nl) 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 calculated which describe the entire population of cells analyzed as illustrated in FIGURES 14d and 14e. First, a mean cell area parameter (MCA) is calculated as a weighted average by multiplying the percentage of a subpopulation (i.e., NC1-NC5 being first divided by 100) by the mean cell area for that subpopulation for each subpopulation and adding the products to produce the weighted average. For example, the percentage of biconcave cells (NCI) is multiplied by the mean -cell area (MCAl) 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.
As explained hereinbefore, the general method of measuring the mean cell size to provide the close correlation to the MCVs achieved by the Coulter counter, is shown in FIGURE 15. The logic section shown in FIGURE 15 thus will receive as inputs the MCA over line 512 from FIGURE 14d, the MCH over line 511 from FIGURE 14e, and the PAL over line 510 from FIGURE 14c. More specifically, the output for the pallor volume PAL may be applied from FIGURE 14c as input to a multiplier logic section 400 (FIGURE 15) which also receives an input factor -K3 to provide a calculation in form of an output on line 401 leading to the accumulator 403. Since the pallor volume is a negative value it will be subtracted in the accumulator 403. Likewise, the mean cell areas MCA for the entire population from logic section 310 in FIGURE 14d is applied over line 512 as input to a multiplier logic section 405 (FIGURE 15) along with the input factor K2 to provide an output of their product to the accumulator 403. The mean cell hemoglobin (MCH) provided over line 511 from summation logic 318 (FIGURE 14e) serves an input ot a multiplier logic section 411 along with a factor K2 to provide an output to the summation logic 403. The output of the summation logic 403 is applied over a line 413 to an adder 415 to which is also applied an input factor K4 which is the offset factor. The output from the adder logic 415 is the mean cell volume for the total population of cells. Typically, the mean cell volume is reported at 430 by printing it out on a form or by displaying it on a cathode ray tube.
In a like manner the computation of mean cell volume (MCV) for a subpopulation of cells can be computed for any given subpopulation of cells. For instance, the biconcave cells are separately classified in the example given herein and in the aforesaid application. More specifically, the mean cell area (MCA 1) for the biconcave cells is provided by logic section 288 in FIGURE 14a and this may be applied as an input on line 512 of FIGURE 15. Likewise, the mean cell hemoglobin (MCH 1) for the biconcave cells is provided from logic section 286 (FIGURE 14a) and this may be applied as an input on line 511 of FIGURE 15. The mean central pallor volume (PAL) for the biconcave cells is available from logic section 264 in FIGURE 14c and this may be applied as input over line 510 to the multiplier logic section 400 in the logic shown in FIGURE 15. The. respective constants will be applied as Kl, K2 and K3 to the respective multiplier logic sections 405, 411 and 400. The remaining operation of the logic section shown in FIGURE 15 will be as above described above to provide an output which is the mean cell volume (MCV) for the biconcave subpopulation of cells.
The logic section shown in FIGURE 14d may also be used with various other subpopulations of cells to provide mean cell volume data for various abnormal cell populations. Because such data has not been heretofore available, there is a whole new data base for diagnosis of cell disorders, blood diseases or of cell morphological changes. For such abnormal cell populations it is preferred to generate for each cell its cell volume, its hemoglobin, and its pallor measurement and then after the cells are classified to accumulate in logic sections these values and then determine the mean cell area (MCA) , mean cell hemoglobin (MCH) , and mean central pallor volume (PAL) . After this, a logic section, such as shown in FIGURE 15 is used to provide a mean cell volume report or output for a given abnormal cell population.
The constants -Kl, K2, K3 and K4 are exemplary of the constants which may be used and these constants have been derived for red blood cells which have been spread by a spinner and then dried before image analysis, as above described. The factors tested herein automatically take into account the cell distortion due to drying of the cells prior to examination. On the other hand, if the red blood ells where kept wet and analyzed, e.g. while in a liquid stream, the constants would be different in order to adjust the mean cell volumes to that achieved by another piece of equipment such as a Coulter Model S counter. Of course, for the MCV of the present invention to accurately correspond to the MCV obtained by a Wintrobe process, the constants would be different as the Coulter counter MCV and the Wintrobe MCV for a given blood specimen will vary. Whether the Coulter counter or the Wintrobe process gives a more true and more accurate depictation of the MCV is not known. It is clear however, that neither the Coulter nor the Wintrobe processes analyze the mean cell pallor of individual cells, as can be done with the present invention, and use such central pallor data in calculating the MCV for a given red blood cell specimen. The factors Kl, K2, K3, and K4 were obtained by using a standard multiple linear regression technique as described fully in a publication entitled Nunsrical Method For Scientists and Engineers by R.W. Hamming, published by McGraw Hill in 1962, and by using the mean cell volume data for the same blood specimens as had been previously measured with a Coulter Model S counter. Other techniques such as comparison of data obtained emperically, may be used to develop a correlation between the image analysis mean cell volumes and the mean cell volumes obtained with a conventional mean cell volume measuring technique such as the Wintrobe technique or the electrical impedance technique of Coulter.
It is to be recognized that equipment has been developed other than that described herein, which measures by image analysis and which provides a cell size output in terms such as microcytic, normocytic, or macrocytic rather than the preferred size output report of a mean cell volume measurement. It is understood however that such equipment has not been accurate in that the results obtained were not consistent _with the results obtained by conventional equipment which provides an MCV output. The present invention may easily provide such a classification of cells. For instance, the output from logic section 415 may be used to classify cells as microcytic, normocytic or macrocytic by having the output of logic section 415 applied to a logic section (not shown) having three levels with MCV below a given level being classified as microcytic, with the MCV's in a central range as normocytic, and with the* MCV's above an upperlevel of the normocytic range being classified as macrocytic.
On the other hand, the size or volume information already being generated by such equipment may be correlated to that of conventional equipment producing MCV data by using the techniques herein disclosed. Thus, the present invention is not to be construed as being limited to the equipment herein described or to equipment that provides an output only in the terms of mean cell volume (MCV) .
In both U.S. Patent 4,097,845 and U.S. Patent 4,199,748 it is pointed out that hard wired logic can be used or that a computer could be used and a specific computer was identified and a long computer program was attached as part of the specification. The computer programs already provided in those disclosures and the information in this disclosure will provide a description sufficient to one skilled in the art to enable the making of a program without undue additional work or experimentation. Hence, the inclusion of another program is submitted not to be warranted and to only result in additional and superfluous material. The present invention likewise may be made in hard wired form without such a computer program and hence the need for a computer program is superfluous for that reason also.
From the foregoing, it will be seen that the present invention provides a new and improved method and apparatus for generating cell size information correlated with cell size information such as mean cell volumes generated with conventional equipment. Although the cell size information may be reported out as placing the cells in a given size category such as microcytic, normocytic, or macrocytic, it preferably provides a mean cell volume output. Additionally, the present invention may be used to provide such size information for a subpopulation of red blood cells.

Claims

CLAIMS :
1. Apparatus for producing a signal representing the mean cell volume from a specimen of red blood cells, some of which have central pallors, said apparatus comprising: means for generating signals representative of the area of the cells; means for measuring the optical density of the individual cells and for generating signals representative of the hemoglobin content or mass of said cells; means for generating signals representative of the central pallor of the cells having central pallors; and means for generating an output representative of the mean cell volume of said cells using said representative area signals, hemoglobin signals, and central pallor signals.' 2. A method for producing a signal representing the mean cell volume from a specimen of red blood cells, some of which have central pallors, said method comprising the steps of: generating signals representative of the area of the cells; measuring the optical density of the individual cells and generating signals representative of the hemoglobin content or mass of said cells, generating signals representative of the central pallor of the cells having central pallors; and generating an output representative of the mean cell volume of said cells using said representative area signals, hemoglobin signals, and central pallor signals. 3. A method of determining the mean cell volume of red blood cells of a particular subpopulation from a blood specimen, comprising the steps of: examining a plurality of red blood cells in the blood, in blood specimen; classifying individual red blood cells by multiple respective features thereof into a plurality subpopulations; and determining the mean cell
c: volume for a given subpopulation of red blood cells.
4. An apparatus for determining the mean cell volume of red blood cells of a particular subpopulation from a blood specimen, said apparatus comprising means
' 5 for examining a plurality of red blood cells in blood specimen; means for classifying individual blood cells by their individual shape and pallor features into a plurality of subpopulations; and means for determining the mean cell volume for a given subpopulation of red
10 blood cells.
5. Apparatus for producing a signal representing the mean cell size from a specimen of red blood cells some of which have central pallors, said apparatus comprising: means for examining said cells by
15 image analysis; means for generating' signals representative of the size of the cells; means for modifying the size information signals of said cells to produce size information signals corresponding to size information from conventional equipment; and means for
20 generating an output representative of the size range of said cells.
6. A method of determining a cell size classification of red blood cells from a blood specimen, comprising the steps of: examining a plurality of red
25 blood cell in the blood specimen by image analysis; extracting size features and generating signals representative of the size of said cells; generating signals for correlating the measured size to that achieved by conventional apparatus; determining a cell
30 size classification for said cells and correlated with the cell size from conventional equipment; and reporting out a cell size classification for said cells.
7. An apparatus in accordance with Claim 6 in which the step of determining a cell size classification
35 comprises the further steps of classifying the cell size information as being representative of microcytic, normocytic or macrocytic cells.
8. A method of determining the mean cell volume of red blood cells from a blood specimen _ comprising the steps of: examining a plurality of red blood cells in the blood in blood specimen; examining individual blood cells by respective features thereof including central pallors, determining a central pallor measurement indicating the percent of volume indentation due to the central pallor, generating size information for said cells, generating hemoglobin information for said cells, and determining the mean cell volume for a given subpopulation of cells by calculating the volume using said hemoglobin and said size information and deducting the central pallor measurement due to central pallor indentation.
"9. A method in accordance ;with Claim 8, including determining the mean cell area for individual cells and determining the mean cell hemoglobin for individual cells.
10. Apparatus for producing a signal representing the mean cell size from a specimen of red blood cells some of which have central pallors, said apparatus comprising: means for generating signals representative of the area of the cells; means for measuring the optical density of the individual cells and for generating signals representative of the hemoglobin content or mass of said cells; means for generating signals representative of the central pallor of the cells having central pallors; and means for generating an output size classification representative of the mean cell size of said cells using said representative are signals, hemoglobin signals, and central pallor signals.
PCT/US1981/000510 1980-04-21 1981-04-20 Method and apparatus for measuring mean cell volume of red blood cells WO1981003080A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US06/142,124 US4453266A (en) 1980-04-21 1980-04-21 Method and apparatus for measuring mean cell volume of red blood cells
US142124 1980-04-21

Publications (1)

Publication Number Publication Date
WO1981003080A1 true WO1981003080A1 (en) 1981-10-29

Family

ID=22498630

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1981/000510 WO1981003080A1 (en) 1980-04-21 1981-04-20 Method and apparatus for measuring mean cell volume of red blood cells

Country Status (5)

Country Link
US (1) US4453266A (en)
EP (1) EP0051633A1 (en)
JP (1) JPH0475463B2 (en)
CA (1) CA1161271A (en)
WO (1) WO1981003080A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2254414A (en) * 1991-03-21 1992-10-07 Univ London Volume measurement of microbial organisms.
CN103823051A (en) * 2008-03-21 2014-05-28 艾博特健康公司 Method and apparatus for determining red blood cell indices of blood sample by utilizing the intrinsic pigmentation of hemoglobin contained within the red blood cells
US9588033B2 (en) 2011-04-15 2017-03-07 Roche Diagnostics Hematology, Inc. Measuring volume and constituents of cells

Families Citing this family (56)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4596464A (en) * 1983-10-14 1986-06-24 Ortho Diagnostic Systems, Inc. Screening method for red cell abnormality
US4593406A (en) * 1984-01-16 1986-06-03 The United States Of America As Represented By The United States Department Of Energy Automatic image acquisition processor and method
US4998284A (en) * 1987-11-17 1991-03-05 Cell Analysis Systems, Inc. Dual color camera microscope and methodology for cell staining and analysis
US4741043B1 (en) * 1985-11-04 1994-08-09 Cell Analysis Systems Inc Method of and apparatus for image analyses of biological specimens
US5109429A (en) * 1985-11-04 1992-04-28 Cell Analysis Systems,Inc. Apparatus and method for analyses of biological specimens
US5016283A (en) * 1985-11-04 1991-05-14 Cell Analysis Systems, Inc. Methods and apparatus for immunoploidy analysis
US5086476A (en) * 1985-11-04 1992-02-04 Cell Analysis Systems, Inc. Method and apparatus for determining a proliferation index of a cell sample
US5134662A (en) * 1985-11-04 1992-07-28 Cell Analysis Systems, Inc. Dual color camera microscope and methodology for cell staining and analysis
US4794450A (en) * 1986-08-29 1988-12-27 Fujirebio Kabushiki Kaisha Method and apparatus for judging agglutination pattern
DE3718066A1 (en) * 1987-05-29 1988-12-08 Zeiss Carl Fa METHOD FOR MICROINJECTION IN CELLS OR. FOR SUCTION FROM SINGLE CELLS OR WHOLE CELLS FROM CELL CULTURES
US4907156A (en) * 1987-06-30 1990-03-06 University Of Chicago Method and system for enhancement and detection of abnormal anatomic regions in a digital image
US5121436A (en) * 1987-08-14 1992-06-09 International Remote Imaging Systems, Inc. Method and apparatus for generating a plurality of parameters of an object in a field of view
JP2686274B2 (en) * 1988-03-24 1997-12-08 東亜医用電子株式会社 Cell image processing method and apparatus
US4918739A (en) * 1988-08-12 1990-04-17 Maraven, S.A. Process and system for digital analysis of images applied to stratigraphic data
US5061632A (en) * 1989-01-31 1991-10-29 Board Of Regents, The University Of Texas System Capillary tube hemoglobinometer and oximeter
US6262798B1 (en) 1992-09-29 2001-07-17 Board Of Regents, The University Of Texas System Method and apparatus for direct spectrophotometric measurements in unaltered whole blood
JP2955873B2 (en) * 1989-08-10 1999-10-04 富士写真フイルム株式会社 Image processing device
US5212637A (en) * 1989-11-22 1993-05-18 Stereometrix Corporation Method of investigating mammograms for masses and calcifications, and apparatus for practicing such method
JP2874008B2 (en) * 1990-05-18 1999-03-24 スズキ株式会社 Particle aggregation pattern judgment device
US5218645A (en) * 1991-03-29 1993-06-08 Cell Analysis Systems, Inc. Method and apparatus for separating cell objects for analysis
US6058209A (en) * 1991-09-27 2000-05-02 E. I. Du Pont De Nemours And Company Method for resolving redundant identifications of an object
US5375177A (en) * 1991-09-27 1994-12-20 E. I. Du Pont De Nemours And Company Method of identifying and characterizing a valid object by color
US5371810A (en) * 1991-09-27 1994-12-06 E. I. Du Pont De Nemours And Company Method of determining the interior points of an object in a background
US5481620A (en) * 1991-09-27 1996-01-02 E. I. Du Pont De Nemours And Company Adaptive vision system
US5579409A (en) * 1991-09-27 1996-11-26 E. I. Du Pont De Nemours And Company Methods for determining the exterior points of an object in a background
US5982944A (en) * 1991-09-27 1999-11-09 E. I. Du Pont De Nemours And Company Adaptive vision system using dual thresholding
WO1993018182A1 (en) * 1992-03-09 1993-09-16 Difco Laboratories Diagnostic microbiological testing apparatus and method
US5889881A (en) * 1992-10-14 1999-03-30 Oncometrics Imaging Corp. Method and apparatus for automatically detecting malignancy-associated changes
US6026174A (en) * 1992-10-14 2000-02-15 Accumed International, Inc. System and method for automatically detecting malignant cells and cells having malignancy-associated changes
US5754688A (en) * 1993-06-04 1998-05-19 Eli Lilly And Company Method and apparatus for stereologic analysis of two-dimensional images
US5598842A (en) * 1993-09-03 1997-02-04 Toa Medical Electronics Co., Ltd. Non-invasive blood analyzer and method using the same
US5791345A (en) * 1993-09-03 1998-08-11 Toa Medical Electronics Co., Ltd. Non-invasive blood analyzer
JP3364323B2 (en) * 1994-05-17 2003-01-08 謙 石原 Non-invasive blood analyzer
US5848177A (en) * 1994-12-29 1998-12-08 Board Of Trustees Operating Michigan State University Method and system for detection of biological materials using fractal dimensions
CA2235772C (en) * 1995-10-23 2002-12-31 Cytometrics, Inc. Method and apparatus for reflected imaging analysis
US7308126B2 (en) * 1997-08-28 2007-12-11 Icad, Inc. Use of computer-aided detection system outputs in clinical practice
US6970587B1 (en) 1997-08-28 2005-11-29 Icad, Inc. Use of computer-aided detection system outputs in clinical practice
US5999639A (en) * 1997-09-04 1999-12-07 Qualia Computing, Inc. Method and system for automated detection of clustered microcalcifications from digital mammograms
US6137898A (en) * 1997-08-28 2000-10-24 Qualia Computing, Inc. Gabor filtering for improved microcalcification detection in digital mammograms
US6801645B1 (en) 1999-06-23 2004-10-05 Icad, Inc. Computer aided detection of masses and clustered microcalcifications with single and multiple input image context classification strategies
JP2003263466A (en) * 2002-03-08 2003-09-19 Umc Japan Electronic cad device, layout data preparing method therefor, and computer program
US7274809B2 (en) * 2002-08-29 2007-09-25 Perceptronix Medical, Inc. And British Columbia Cancer Agency Computerized methods and systems related to the detection of malignancy-associated changes (MAC) to detect cancer
US6981083B2 (en) * 2002-12-05 2005-12-27 International Business Machines Corporation Processor virtualization mechanism via an enhanced restoration of hard architected states
US20040202357A1 (en) 2003-04-11 2004-10-14 Perz Cynthia B. Silhouette image acquisition
JP4328586B2 (en) * 2003-09-02 2009-09-09 キヤノン株式会社 Radiation image connection processing method, radiation image connection processing apparatus, computer program, and computer-readable recording medium
ES2929039T3 (en) * 2006-04-08 2022-11-24 Hoffmann La Roche Analysis of optical data with the help of histograms
WO2009126800A1 (en) * 2008-04-09 2009-10-15 Abbott Point Of Care, Inc. Method for measuring the area of a sample disposed within an analysis chamber
US9459196B2 (en) 2011-07-22 2016-10-04 Roche Diagnostics Hematology, Inc. Blood analyzer calibration and assessment
CN102819765B (en) * 2012-02-28 2016-01-20 浙江工业大学 A kind of milk somatic cell method of counting based on computer vision
CN104655546B (en) * 2015-03-12 2017-07-07 重庆大学 A kind of milk somatic cell method of counting and system based on mobile device
CN108474934B (en) 2015-09-17 2022-01-18 思迪赛特诊断有限公司 Method and apparatus for detecting entities in a body sample
CA3018536A1 (en) 2016-03-30 2017-10-05 S.D. Sight Diagnostics Ltd Distinguishing between blood sample components
EP3455610B1 (en) 2016-05-11 2023-01-04 S.D. Sight Diagnostics Ltd. Sample carrier for optical measurements
WO2018211418A1 (en) * 2017-05-15 2018-11-22 Sigtuple Technologies Private Limited Method and system for determining area to be scanned in peripheral blood smear for analysis
JP7214729B2 (en) 2017-11-14 2023-01-30 エス.ディー.サイト ダイアグノスティクス リミテッド Sample container for optical measurement
US11187644B2 (en) * 2018-04-09 2021-11-30 Sigtuple Technologies Private Limited Method and system for determining total count of red blood cells in peripheral blood smear

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3851156A (en) * 1972-09-05 1974-11-26 Green James E Analysis method and apparatus utilizing color algebra and image processing techniques
US3982183A (en) * 1975-02-20 1976-09-21 Coulter Electronics, Inc. Particle sizing apparatus
US4063309A (en) * 1976-09-16 1977-12-13 Hycel, Inc. Mean corpuscular volume measuring apparatus and method
US4086631A (en) * 1974-12-23 1978-04-25 Hycel, Inc. Hematocrit measurement in an automatic hematology analyzer
US4097845A (en) * 1976-11-01 1978-06-27 Rush-Presbyterian-St. Luke's Medical Center Method of and an apparatus for automatic classification of red blood cells
US4159895A (en) * 1977-03-30 1979-07-03 Coulter Electronics, Inc. Method and apparatus for the detection of beta thalassaemia minor

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2656508A (en) * 1949-08-27 1953-10-20 Wallace H Coulter Means for counting particles suspended in a fluid
US3259842A (en) * 1959-08-19 1966-07-05 Coulter Electronics Particle analyzing device
US3944797A (en) * 1974-08-06 1976-03-16 Coulter Electronics, Inc. Method and apparatus for determining the correct percentiles of the size distribution of a particulate system
GB1538225A (en) * 1975-02-28 1979-01-10 Nat Res Dev Methods and apparatus for analysing frequency distributions

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3851156A (en) * 1972-09-05 1974-11-26 Green James E Analysis method and apparatus utilizing color algebra and image processing techniques
US4086631A (en) * 1974-12-23 1978-04-25 Hycel, Inc. Hematocrit measurement in an automatic hematology analyzer
US3982183A (en) * 1975-02-20 1976-09-21 Coulter Electronics, Inc. Particle sizing apparatus
US4063309A (en) * 1976-09-16 1977-12-13 Hycel, Inc. Mean corpuscular volume measuring apparatus and method
US4097845A (en) * 1976-11-01 1978-06-27 Rush-Presbyterian-St. Luke's Medical Center Method of and an apparatus for automatic classification of red blood cells
US4159895A (en) * 1977-03-30 1979-07-03 Coulter Electronics, Inc. Method and apparatus for the detection of beta thalassaemia minor

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2254414A (en) * 1991-03-21 1992-10-07 Univ London Volume measurement of microbial organisms.
CN103823051A (en) * 2008-03-21 2014-05-28 艾博特健康公司 Method and apparatus for determining red blood cell indices of blood sample by utilizing the intrinsic pigmentation of hemoglobin contained within the red blood cells
CN103823051B (en) * 2008-03-21 2015-11-18 艾博特健康公司 Utilize the intrinsic pigmentation of the haemoglobin contained in red blood cell to determine the method and apparatus of the red cell index of blood sample
US9588033B2 (en) 2011-04-15 2017-03-07 Roche Diagnostics Hematology, Inc. Measuring volume and constituents of cells
US10281382B2 (en) 2011-04-15 2019-05-07 Roche Diagnostics Hematology, Inc. Measuring volume and constituents of cells

Also Published As

Publication number Publication date
JPS57500625A (en) 1982-04-08
JPH0475463B2 (en) 1992-11-30
US4453266A (en) 1984-06-05
CA1161271A (en) 1984-01-31
EP0051633A1 (en) 1982-05-19

Similar Documents

Publication Publication Date Title
US4453266A (en) Method and apparatus for measuring mean cell volume of red blood cells
US4199748A (en) Automated method and apparatus for classification of cells with application to the diagnosis of anemia
US4097845A (en) Method of and an apparatus for automatic classification of red blood cells
US5123055A (en) Method and an apparatus for differentiating a sample of biological cells
US5235522A (en) Method and apparatus for automated analysis of biological specimens
CA1135979A (en) Automated method and apparatus for classification of cells with application to the diagnosis of anemia
US5848177A (en) Method and system for detection of biological materials using fractal dimensions
US4596464A (en) Screening method for red cell abnormality
US9170256B2 (en) Device and method for erythrocyte morphology analysis
KR100303608B1 (en) Method and device for automatically recognizing blood cell
CN108682008B (en) A kind of Leukocyte Image clarity evaluation method and device
JP3581149B2 (en) Method and apparatus for identifying an object using a regular sequence of boundary pixel parameters
US6184978B1 (en) Method and apparatus for verifying uniform flow of a fluid sample through a flow cell and distribution on a slide
JP2001502414A (en) Method and apparatus for assessing slide and sample preparation quality
US20110057946A1 (en) Image processing method, image processing apparatus, and computer-readable recording medium storing image processing program
US20100279341A1 (en) Methods and system for analyzing cells
Zahniser et al. BioPEPR: a system for the automatic prescreening of cervical smears.
Rolke et al. Size structure analysis of zooplankton samples by means of an automated image analyzing system
CN111274949B (en) Blood disease white blood cell scatter diagram similarity analysis method based on structural analysis
KR101106386B1 (en) System for classifying slides using scatter plot distributions
CA2185511C (en) Cytological specimen analysis system with individualized patient data
CN114152557A (en) Image analysis based blood cell counting method and system
EP0138591A2 (en) Screening method for red blood cell abnormality
CA1180811A (en) Method and apparatus for measuring mean cell volume of red blood cells
US11788950B2 (en) Method and device for analysing cells

Legal Events

Date Code Title Description
AK Designated states

Designated state(s): JP

AL Designated countries for regional patents

Designated state(s): CH DE FR GB NL SE