WO1999040410A9 - System and device for determining particle characteristics - Google Patents
System and device for determining particle characteristicsInfo
- Publication number
- WO1999040410A9 WO1999040410A9 PCT/US1999/002588 US9902588W WO9940410A9 WO 1999040410 A9 WO1999040410 A9 WO 1999040410A9 US 9902588 W US9902588 W US 9902588W WO 9940410 A9 WO9940410 A9 WO 9940410A9
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- particle
- force field
- determining
- cell
- pole piece
- Prior art date
Links
- 239000002245 particle Substances 0.000 title claims abstract description 115
- 230000005291 magnetic effect Effects 0.000 claims abstract description 77
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000000926 separation method Methods 0.000 claims abstract description 10
- 230000001105 regulatory effect Effects 0.000 claims abstract description 3
- 210000004027 cell Anatomy 0.000 claims description 197
- 238000012545 processing Methods 0.000 claims description 27
- 230000004907 flux Effects 0.000 claims description 25
- 238000002372 labelling Methods 0.000 claims description 19
- 238000004891 communication Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 230000005684 electric field Effects 0.000 claims description 4
- 239000007850 fluorescent dye Substances 0.000 claims description 4
- 210000003463 organelle Anatomy 0.000 claims description 4
- 230000001276 controlling effect Effects 0.000 claims description 3
- 230000000242 pagocytic effect Effects 0.000 claims description 2
- 230000005672 electromagnetic field Effects 0.000 claims 2
- 230000006870 function Effects 0.000 description 43
- 230000005298 paramagnetic effect Effects 0.000 description 16
- 238000004458 analytical method Methods 0.000 description 15
- 238000005070 sampling Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 5
- 239000002458 cell surface marker Substances 0.000 description 5
- 230000001413 cellular effect Effects 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000001900 immune effect Effects 0.000 description 4
- 210000004698 lymphocyte Anatomy 0.000 description 4
- 239000003550 marker Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 3
- 230000009087 cell motility Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000011521 glass Substances 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000027455 binding Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000001943 fluorescence-activated cell sorting Methods 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000000037 particle-tracking velocimetry Methods 0.000 description 2
- 229920001296 polysiloxane Polymers 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000012798 spherical particle Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 108010001857 Cell Surface Receptors Proteins 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 206010057249 Phagocytosis Diseases 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- QJVKUMXDEUEQLH-UHFFFAOYSA-N [B].[Fe].[Nd] Chemical compound [B].[Fe].[Nd] QJVKUMXDEUEQLH-UHFFFAOYSA-N 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 239000000427 antigen Substances 0.000 description 1
- 102000036639 antigens Human genes 0.000 description 1
- 108091007433 antigens Proteins 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 239000005388 borosilicate glass Substances 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 150000004696 coordination complex Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- BFMYDTVEBKDAKJ-UHFFFAOYSA-L disodium;(2',7'-dibromo-3',6'-dioxido-3-oxospiro[2-benzofuran-1,9'-xanthene]-4'-yl)mercury;hydrate Chemical compound O.[Na+].[Na+].O1C(=O)C2=CC=CC=C2C21C1=CC(Br)=C([O-])C([Hg])=C1OC1=C2C=C(Br)C([O-])=C1 BFMYDTVEBKDAKJ-UHFFFAOYSA-L 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- -1 for example Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 238000007885 magnetic separation Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 102000006240 membrane receptors Human genes 0.000 description 1
- 239000002923 metal particle Substances 0.000 description 1
- 239000013528 metallic particle Substances 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910001172 neodymium magnet Inorganic materials 0.000 description 1
- 230000009871 nonspecific binding Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008782 phagocytosis Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 230000009870 specific binding Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000000827 velocimetry Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/10—Investigating individual particles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
- G01N27/416—Systems
- G01N27/447—Systems using electrophoresis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/18—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
-
- G01N15/149—
-
- G01N2015/1027—
-
- G01N2015/1028—
Definitions
- the invention relates generally to methods and apparatuses for determining particle characteristics, and more particularly, to methods and apparatuses having particle tracking and image analysis logic and force field devices for determining one or more physical cell characteristics.
- Cell analysis and separation is an increasingly important technique in the diagnosis and treatment of various cancers and diseases.
- primary importance to cell analysis and separation is the ability to identify, or label, cell properties and characteristics of interest.
- the identification, or labeling, of cell properties and characteristics allows them to be used as " handles” which, in turn, can be used to separate " labeled” cells from other cells.
- labels for sorting cells are immunological labels which include, for example, immunofluorescent and immunomagnetic labels.
- Immunofluorescent labels typically include, for example, a fluorescent molecule joined to an antibody.
- Immunomagnetic labels typically include, for example, a paramagnetic compound or molecule joined to either a primary or secondary antibody.
- Cell labeling is performed by attaching the antibody to a marker of interest on the surface of the cell (i.e., cell surface marker) .
- a marker of interest on the surface of the cell
- cell surface marker i.e., cell surface marker
- immunological labels have been developed which allow for the careful labeling of cells, the potential of these labels for cellular analysis and separation has yet to be fully realized. As a result thereof, the cellular properties and characteristics which these labels are capable of identifying have also yet to be fully analyzed. For example, in the case of immunomagnetic labels, the highly accurate quantification of a cell population' s magnetic susceptibility has been impossible to determine. Additionally, in the general case of immunological labels, the cell surface marker and label density has been difficult to accurately determine.
- apparatuses and methods are provided for quantifying at least one of a plurality of particle physical characteristics
- the range of particles includes, for example, cells, cell organelles, platelets, inorganic, organic, biological, and polymeric particles which are optically visible
- the plurality of particle physical characteristics include particle size, shape, magnetic susceptibility, magnetic label density, charge separation, dielectric constant, and derivatives thereof Derivatives include, for example, cell dimensions ⁇ e . g. , diameter, radius, major axis length, minor axis length, etc.
- the method includes generating a region of space having a substantially constant force field, determining the velocity of at least one particle within the region by identifying and locating the particle and ts coordinates in at least two temporally defined digital images, and determining the at least one of a plurality of particle physical characteristics from the determined velocity and a predetermined force field magnitude and direction. These steps include several steps or sub-steps.
- the step of determining the velocity of at least one particle in a force field includes the step or sub-step of processing the at least two temporally defined digital images so that the particle is distinct from the background of each temporally defined digital image
- This step of processing includes one or more of the following steps or sub-steps histogramming, color stretching, filtering, background subtraction, identifying contrast differences.
- the step of determining the velocity of at least one particle in a force field further includes the step or sub-step of tracking the location of the particle through the at least two temporally defined digital images by determining at least one predicted path and the distance the particle has traveled.
- Other steps or sub-steps are more fully set forth in the detailed description.
- a device for determining one or more particle physical characteristics has a force field device for subjecting at least one particle to at least one force field, a substantially transparent flow channel, and a computer system for gathering and analyzing data associated with the at least one particle.
- the force field device has, among other things, a first and second force field producing assembly. Each assembly has a force field producing device and a pole piece for concentrating the force field flux.
- the pole pieces each have a flux concentrating portion having a curved end portion for producing a region of space having a substantially constant force field and wherein the flux concentrating portions are displaced substantially opposite each other with an inter-polar air gap therebetween.
- the substantially transparent flow channel is positioned at least partially within the region of space and provides for the introduction of at least one particle thereinto.
- a flow system having a pump for controlling flow into the channel is also provided.
- a system for determining one or more particle physical characteristics has a force field device for generating at least one force field having a predetermined force field magnitude and direction and for subjecting at least one particle to the at least one force field, a flow system for regulating the introduction of the at least one particle into the force field device, and a computer system for gathering and analyzing data associated with the at least one particle.
- the computer system has, among other things, a digital image system for acquiring at least two temporally defined digital images of the at least one particle, logic for identifying and locating the at least one particle and its coordinates within the at least two temporally defined images , logic for determining the velocity of the at least one particle within the force field, and logic for determining at least one particle physical characteristic from the determined velocity and the predetermined force field magnitude and direction
- a pole piece assembly for producing a region of space having a substantially constant magnetic force field s
- the assembly has, among other things, a first pole piece having a substantially curved flux concentrating portion, a second pole piece having a substantially curved flux concentrating portion, and a device for producing a magnetic flux in flux communication with the first and second pole pieces
- the first and second pole pieces are configured to form an mter-polar air gap
- the substantially curved flux portions of the first and second pole pieces each have a first distal end having a curved portion and a second distal end having a curved portion
- Figure 1 s a high-level functional block diagram of a quantitative analysis device 100 of the present invention
- Figure 2 is a high-level flowchart of the overall Particle Tracking Velocimetry Analysis logic
- Figure 3 is a flow chart illustrating the image processing logic of the present invention.
- Figure 4 is a diagram of the two-dimensional cell tracking function of the present invention
- Figure 5 is a top plan view of one embodiment of a force field device of the present invention
- Figure 6 is a cross-sectional view taken on section line 6-6 of Figure 5 of the force field device of the present invention.
- Figure 7A is a graph showing the magnetic field lines generated by the present invention.
- Figure 7B is a magnified view of a portion of the force field device of Figure 6 ;
- Figure 8 is a graph illustrating the magnetic force and rate of change of magnetic force versus y-axis position of the force field device of Figures 5, 6, and 7 ;
- Figure 9A is a graph showing the distribution of average cell diameters of lymphocyte cells as determined by the present invention and as determined experimentally via a Coulter Multis zer II;
- Figure 9B s a graph showing the distribution of magnetic susceptibility of lymphocyte cells as determined by the present invention.
- the device 100 has a force field device 102, digital or analog camera 104, computer system 106, input device (s) 118, and output device (s) 116.
- the computer system 106 contains a CPU 108, memory 110, Cell Tracking Velocimetry analysis logic 112 (hereinafter CTV logic 112), video board 114, and various other support devices including, for example, floppy, hard, and CD-ROM drives, fax/modem, interface logic, etc. (not shown).
- the CPU 108 is preferably a CPU which s capable of quickly and efficiently processing large amounts of graphical information.
- One such CPU is the Intel MMX ® microprocessor.
- other microprocessors may also be employed such as, for example, the Intel Pentium ® processor, Intel 80486, or other high-performance microprocessors.
- One such computer system having the above-described features, or easily modifiable to have the above-described features is, for example, an IBM APTIVA ® or Gateway 2000 ® personal computer system with an Intel MMX ® microprocessor.
- the video board 114 of the present invention preferably has the ability to capture and process high resolution video and image information from external devices (e.g., digital or analog cameras) or internal devices ( e . g . , built-m video tuner).
- the video board 114 preferably contains a video processor, video memory and various interface logic for interfacing to digital and/or analog cameras, display devices, and the computer system 106.
- One such suitable video board 114 is the P360F Power Grabber manufactured by DPIX Technologies Inc .
- Other high performance video boards having similar characteristics may also be employed.
- the output device (s) 116 mclude(s), for example, display monitors, printers, and external storage devices.
- the input device (s) 118 include (s), for example, a keyboard, mouse, and voice input.
- the digital camera 104 preferably contains the ability to capture high-resolution monochrome and/or color video and still- image information.
- a suitable digital camera is, for example, the CCD 4915 camera, manufactured by Cohu Inc. of San Diego, California, with a 2.5X and a 6.7X magnification insert. Other cameras having similar characteristics may also be employed. Additionally, magnifications other than 2.5X or 6.7X may also be utilized, depending on the size of the cells or particles being analyzed and the required magnification.
- a high-resolution video storage device 120 such as, for example, a sVHS video recorder or optical disk(s) may be employed to store the video and images for subsequent analysis and archiving.
- the light source 124 may be a source coherent light of one or more wavelengths such as, for example laser light, or noncoherent light of one or more wavelengths such as, for example, white light, colored light, ultra-violet light, or infra-red light, or combinations thereof .
- the force field device 102 is responsible for a number of important functions within the present invention. Firstly, the force field device 102 allows the cells or particles to be subjected to and displaced by a force field. Secondly, the force field device 102 allows for the cell or particle displacement to be viewed and/or captured by the digital camera 104. Thirdly, the force field device 102 may, in certain embodiments, provide the force field being applied (e.g., a magnetic field). Other functions will also be apparent from the detailed discussion presented below.
- the CTV logic 112 of the present invention analyzes a plurality of closely-timed cell or particle images to determine, among other things, individual cell or particle locations and velocities through each captured image. All the data collected and determined by the CTV logic, including the captured image information, is preferably stored in a database for subsequent viewing, analysis, or review. As will be described in more detail, once the cell locations and velocities are known, many other characteristics such as, for example, cell acceleration, force, and mass may be determined therefrom. It should also be noted that the CTV logic 112 is not dependent upon the specific type of force field applied to the cells or particles. Therefore, the CTV logic may be employed in systems where the force field device is employing magnetic, electric, or gravitational fields. The present discussion will hereinafter focus on the analysis of cells, however, it should be noted that the discussion is equally applicable to other particles such as, for example, cell organelles, and other metallic/non-metallic particles.
- the CTV logic begins n step 202 where the digital camera 104 (shown m Figure 1) acquires a number of images. Firstly, the digital camera 104 acquires at least one background image . A background image s an image of the observation area of the force field device 102 which typically contains a glass tube or channel through which cells are or will be visible The background image s used m a background subtraction function of the image processing logic See the text associated with Figure 3 for a detailed discussion of the background subtraction function. Secondly, the digital camera 104 acquires a plurality of closely timed digital images of the cells.
- Each digital image is acquired at an image sampling rate generally n the range of 15 to 60 Hz, with a preferable image sampling rate of 30 Hz.
- These digital images are also known as, and commonly referred to as, an image " frame” or " frames.”
- the video board 114 converts each pixel of each image frame into one of a plurality of digital image formats which convey pixel brightness and/or color information.
- the digital formats range from 8 to 24 bits of information per pixel.
- the preferred image format is 8 bits of gray-level information per pixel.
- the 8 bits of gray-level information provide a range of gray-level values from 0 ⁇ i . e.
- Each image frame contains a high-resolution pixel array, preferably m the general range of 600 x 400 to 1280 x 1024 pixels.
- the digital cameral of the illustrated embodiment provides a high- resolution 624 x 450 pixel array
- the digital image sampling rate of 30 Hz may be increased or decreased depending on the required image frame resolution. For example, the sampling rate of 30 Hz produces 30 image frames per second. If a higher image frame resolution is required, the sampling rate may be increased to, for example, 60 Hz (i.e., 60 image frames per second.
- the CTV logic advances to step 204 where the image frames are processed for standardization which facilitates the tasks of locating and identifying individual cells and determining their velocities
- the CTV logic advances to step 204 where the image frames are processed for standardization which facilitates the tasks of locating and identifying individual cells and determining their velocities
- the CTV logic advances to step 204 where the image frames are processed for standardization which facilitates the tasks of locating and identifying individual cells and determining their velocities
- FIG. 3 a flow chart 300 illustrating the image processing logic of the present invention is shown.
- the acquired physical images from the digital camera 104 are not optimal for cell tracking because the gray-level differences between the background image and the cell image (s) are not distinct. Accordingly, some degree of image processing is required.
- the image processing logic executes a plurality of steps including, but not limited to, for example, histogramming, stretching, spatial filtering, background subtraction, pattern filtering, and pixel size matching. Accordingly, the image processing logic starts m step 302 where a histogramming function is executed for each image frame. Histogramming is a statistic measure of the frequency of the gray- level versus the gray-level itself. The range of gray-level is from 0 ( i . e .
- a range of dominant gray-levels may be determined from the frequency of their occurrence.
- the range of dominant gray-levels n turn, indicate whether the cell image s distinct from the background image of the image frame. For example, if the range of dominant gray-levels fall between 175 and 250, it can be said that the cell image is not distinct from the background image. However, if the range of dominant gray-levels comprises two sub-ranges with one sub-range localized near the gray-level of 0 ⁇ i . e.
- the cell image is distinct from the background image of the image frame. It should be noted that whether the cell image is distinct from the background image s a matter of degree. In an ideal image frame, the cell image would appear white m color (i.e., gray-level 255) and the background image would appear black m color (i.e., gray- level 0) . It is a goal of the image processing logic to standardize all image frames as close as possible to this norm.
- step 302 the image processing logic advances to step 304 where a stretching function is executed for each image frame.
- the stretching function is employed when the distribution of the gray-level m the image frame does not cover the full brightness range and produces poor contrast between the cells and the background.
- the stretching function is accomplished by setting a minimum and a maximum gray-level value.
- the minimum gray-level value is set to the low-end and the maximum gray-level value is set to the high-end of the general range of dominant gray-level values as determined by the histogramming function.
- all image frame pixels with a gray-level value less than the minimum gray-level value are set to 0 (i.e., black) .
- All image frame pixels with a gray-level value greater than the maximum gray- level value are set to 255 ( i . e . , white) .
- All image frame pixels with gray-level values in the range between the minimum and maximum are stretched proportionally, to a value between 0 and 255, based on their original gray-level value.
- the result of the stretching function is an image frame that contains distributions of gray- level values which distinctly correspond to the cell images and the background image
- step 306 the image processing logic executes a low-pass spatial filtering function for each image frame.
- the low-pass spatial filtering function removes small details, or noise, from the image frame.
- Low-pass spatial filtering is also known as " neighborhood averaging" and is used to reduce noise by smoothing local gray-levels of the image.
- the image processing logic advances to step 308. Since the CTV logic requires image frames which contain only bright cell images and a dark background image, a background subtraction function is executed for each image m step 308.
- the background subtraction function is particularly useful because most experimental systems include extraneous matter such as , for example, dirt on the surfaces of camera lenses and other observation devices such as, for example, glass tubes or channels.
- each image frame In the background subtraction function, an image frame of the background is subtracted from the image frame containing the cell images and the background image.
- one of the first images acquired n step 202 of Figure 2 is an image of the background of the scene where the cells w ll eventually appear. That is, the background image is an image of the observation area of the force field device which typically contains a glass tube or channel through which cells will be visible. Accordingly, after the background subtraction function n step 308 is performed, each image frame generally contains bright cell images and a dark background image .
- step 310 a cell template function is applied to each cell of each image frame.
- a cell template is a standard complete cell image which has been pre-defined based on an actual physical image of the cell.
- the need to apply the cell template function arises because, after image processing, the cell images may appear as hollow spheres with their interiors having the same gray-level as the background image. Consequently, after application of the cell template function, the cell images no longer appear as hollow spheres but as solid bright spheres.
- application of the cell template function may be optional because after image processing, some cell images may already appear as solid bright spheres. Therefore, in such cases there would be no need to apply the cell template function.
- step 206 two-dimensional cell identification and location s performed for each image frame. This step involves the separation of cell images from the background image and the determination of cell location within the image frame. Cell locations are identified based on the contrast difference between the cell image and the background image
- the cell image is defined as a shape connected set of pixels having an intensity greater than a given threshold gray-level value.
- the threshold gray-level value can fall within a range of possible values such as, for example, 0 ( .e., black) to 79. Therefore, in the given example, all pixels having a gray-level value greater than 79 would have an intensity greater than the given range of threshold gray-level values and would, therefore, be interpreted as belonging to a cell image.
- This result is achieved by scanning the pixels of each image frame and determining whether their intensity is greater than the threshold value .
- the CTV logic determines the location of the center of each cell within each image frame. Th s process is facilitated by using the pixel array of each image frame as a coordinate system. Additionally, since each cell image is a solid sphere, or very nearly thereto, the center of the cell image is easily determined by locating the center of the sphere, which appears as a circle m two-di ensions . These locations are stored m a database for use n the two dimensional cell tracking function of step 208. Additionally, the CTV logic may include a function which allows for the entry of minimum and maximum cell size and aspect ratio data as a criteria for further cell image recognition. In step 208, the CTV logic executes the two-dimensional cell tracking function.
- This function employs a sequence of five image frames to establish, therefrom, the most probable path (as a function of time) taken by each cell.
- the most probable path determination is based on the concept of " path coherence . " Generally, the concept of " path coherence” holds that cell position, velocity, acceleration, change of acceleration, etc . , are internally self-consistent and that these characteristics can be described by mathematical functions that are smooth.
- the two-dimensional cell tracking function of step 208 will now be described with reference to Figure 4. Illustrated in Figure 4 are the positions of a single cell as determined from three consecutive image frames. The position of the cell in the first image frame is indicated by Cl, in the second image frame by C2 , and in the third image frame by C3.
- a search radius " r" is established, thereby establishing search area 402, and will be used for identification of cell images in the second image frame.
- the value of the search radius " r" is dependent on a number of factors including, for example, image sampling frequency and cell displacement between image frames. Therefore, if the cell displacement between image frames is relatively small, a small search radius " r" can be used. For example, " r" can be somewhere in the general range of 1 to 3 times the cell diameter. The ascertained cell diameter is preferably in the range of l ⁇ m or greater. However, if the cell displacement between image frames is relatively large, a larger search radius " r" would be required ⁇ e . g.
- the same search radius " r" is established for every cell of an image frame. Accordingly, once the search radius " r" has been established in the first image frame, all cells in the second image frame within the search radius are tagged and identified. As shown in Figure 4 , cell image C2 is within the search radius and search area 402. Since the coordinates of the center of cell image C2 are known, the distance Dl and vector direction from cell image Cl can be determined. These values are used to determine a candidate path 404 for the cell. Therefore, from the first and second image frame, the path and location of the cell in the third frame is predicted (shown at point " P" on predicted path 404) .
- the path and location of the cell is compared to the actual position of the cell image C3 in the third frame.
- the error between the predicted (i.e., point " P" ) and actual cell image position C3 is determined and added to a penalty function.
- the penalty function is a measure of path coherence. That is, the smaller the penalty function, the greater the path coherence. Therefore, the error between cell images Cl and C2 is zero by default since at least 2 points are required for a path.
- an error may exist, as shown in Figure 4, between the predicted third image frame cell location (i.e., point " P" on predicted path 404) and the actual third image frame cell location C3. Since the coordinates for the predicted cell location " P" and the actual cell location C3 are known, an error value may be determined therefrom.
- This procedure is repeated with the fourth image frame using the data acquired from the first three image frames and, similarly, for the fifth image frame using the data acquired from the first four image frames.
- the CTV logic may have to analyze several possible paths if more than one cell image has been tagged and identified within a search radius. In such situations, the values of the penalty functions for each separate possible path are compared and the path with the smallest penalty function value is chosen as the correct path.
- the bias value is a measure of reliability that a determined cell path is the correct cell path. For example, if drastic differences between the forward analysis penalty function value and the reverse analysis penalty function value are present, a large bias in one of the two directions of analysis exists -- thereby indicating that the presently determined cell path may be less reliable than other potential candidate cell paths . Similarly, a small bias value would indicate that the presently determined cell path is more reliable than other potential candidate cell paths. Accordingly, cell paths having small penalty function values and small, or no, bias values are reliable as determined or most probable cell paths.
- a number of cell characteristics can be determined by the CTV logic 112 such as, for example, the velocity of the cell from image frame to image frame. Furthermore, since the velocity of the cell is known at a plurality of different times, the acceleration of the cell can also be determined. Still further, if the mass of the cell is known, the force acting on the cell can be determined or vice-versa. Therefore, the CTV logic of the present invention can determine a wide range of physical parameters including, for example, all of the above-mentioned parameters and values.
- particle physical characteristics such as, for example, magnetic susceptibility, charge cell separation and size, may be determined from the above physical parameters and values.
- the CTV logic 112 can determine an average value, the standard deviation, and the range for the cell population analyzed through a plurality of known means. These values are the ultimate outputs of the CTV logic 112 and may be output to a display monitor, printer, or storage device.
- the force field device 102 provides a number of important functions within the present invention such as, for example, allowing the cells to be subjected to and displaced by a force field and allowing for cell displacement to be viewed and/or captured by a camera.
- the force field device 102 of the present invention may employ any one of a variety of force fields including, for example, flow, magnetic, electric, and gravitational fields.
- One illustrated embodiment of a force field device 102 employing a magnetic field is shown in Figures 5, 6, 7A, and 7B.
- the present invention can determine, for example, the magnetic susceptibility of cells. Once determined, the magnetic susceptibility can be used to appropriately design magnetic field (s) and the magnetic assemblies of cell sorting devices. This is of particular importance for fractional cell sorting devices which sort cells based on the density of magnetically labeled cell surface markers.
- Force field device 102 for subjecting cells to a magnetic force field is shown.
- Figure 5 is a top plan view and
- Figure 6 is a cross-sectional view taken on section line 6-6 of Figure 5.
- Force field device 102 includes a base 606, two pairs of force field producing magnets 602 and 604, and two pole pieces 502 and 504.
- the magnets 602 and 604 are generally rectangular in shape and preferably in the range of about 2" x 2" x 1" and preferably made from a permanent magnetic material such as, for example, a neodymium-iron-boron alloy.
- the base 606 and the pole pieces 502 and 504 are preferably made from a material that is capable of having a magnetic flux induced therein such as, for example, 1018 low-carbon cold- finished steel.
- the pole pieces 502 and 504 each have a flux concentrating portion 516 and 518, respectively, and top and bottom surfaces 608 and 610, respectively.
- the distance between top and bottom surfaces 608, 610 and 612, 614, is approximately in the range of 12-13 mm, with a preferred distance of approximately 12.5 mm.
- the base 606, magnets 602 and 604, and pole pieces 502 and 504 are assembled as shown in Figures 5 and 6 so as to form an inter-polar air gap 505 and utility spaces 506 and 508.
- utility space 506 may be used for placing the digital camera 104 in close proximity to the inter-polar air gap 505.
- utility space 508 may be used for providing a light source or mirror 510 which directs light through the inter-polar air gap 505.
- a channel 514 is positioned within the inter-polar air gap 505 such that cells therein are subjected to a substantially uniform magnetic field.
- the channel 514 is made from an optically clear inert material such as, for example, borosilicate glass.
- the channel 514 is substantially rectangular in cross-section.
- FIG. 7A a graph showing the magnetic field lines versus x and y position generated by the present invention is illustrated.
- the location of the region of constant magnetic force 702 is shown relative the inter- polar air gap, along with the magnetic field lines B .
- FIG. 7B a magnified view of portion 7 of Figure 6 is shown.
- the flux concentrating portions 516 and 518 have end surfaces 704 and 706.
- End surfaces 704 and 706 each have curved distal ends 712, 714 and 716, 718, respectively.
- Curved distal ends 712 and 716 have a radius in the general range of 3 mm with a preferable radius of 3.18 mm.
- Curved distal ends 714 and 718 are generally hyperbolic and preferably defined by the hyperbolic function:
- distal ends 712 and 714 are configured such that they are continuous with each other.
- Distal ends 716 and 718 are also similarly configured.
- an approximately linear surface may be provided between distal ends 712 and 714 so as to also provide a continuous joinder of the distal ends.
- Distal ends 716 and 718 may also be similarly joined by an approximately linear surface.
- the channel 514 is preferably placed within the inter-polar air gap 505 such that it is very nearly in physical contact with the end surfaces 704 and 706. Alternatively, the channel 514 may be placed in actual physical contact with end surfaces 704 and 706.
- the width of the inter-polar air gap 505 is, at its narrowest, approximately 2 mm So configured, the force field device 104 generates a region 702 of substantially constant magnetic force which is exerted onto the cells present within this region.
- the gradient of the magnetic field energy (V(B 2 /2 ⁇ 0 )) is generally illustrated at 708 and the constant magnetic energy lines (B 2 /2 ⁇ 0 ) are generally illustrated at 710.
- the region 702 of Figure 7B is generally indicated m Figure 8 and falls within the general range of -7 5 to -9.5 mm from the top of the pole pieces 502 and 504. It should be noted that changes in pole piece geometry will affect the y-axis position of the substantially constant magnetic force region 702.
- the force field device 102 of the present invention is n fluid communication with a flow system 122 (shown m Figure 1) which provides for the injection of cells and a carrier medium into the channel 514 and the flow field.
- the flow system preferably contains a disposable injection syringe pump which is m fluid communication with the channel 514 v a silicone tubing.
- the other end of the channel 514 is in flow communication with a waste vessel also via silicone tubing
- the inlet portion of channel 514 may be in fluid communication with an injection device such as, for example, a disposable syringe pump
- the exit portion of the channel 514 may be fluid communication with an aspirating device such as, for example, an aspirating syringe pump
- an injection device such as, for example, a disposable syringe pump
- an aspirating device such as, for example, an aspirating syringe pump
- the flow rate generated by the two pumps being controlled by the computer system 106 or other control logic within the pumps.
- the present invention may be used to determine a plurality of cell characteristics including, for example, cell size and magnetic susceptibility
- cell size including, for example, cell size and magnetic susceptibility
- the present invention can be used to determine data from which equation (1) can be solved for the cell diameter D c .
- the CTV logic analyzes captured images of cells which are subjected to only a gravitational force field (as opposed to an electric or magnetic force field) and determines the cell settling velocities therefrom. Once the cell settling velocities V c are known, along with the other values of Equation (1) , the diameter of the cell D c is determined.
- Illustrated m Figure 9A s a graph showing the distribution of cell diameters of lymphocyte cells as determined by the present invention (i.e., solid line) and as determined experimentally via a commercial cell sizing device ( i . e . , a Coulter Multisizer II)
- the Coulter Multisizer II was used to determine cell diameter D c .
- the primary peak and shoulder generated by each method are in close agreement with each other.
- the large peak around 20 x 10 "6 m corresponds to clumps of cells which can be observed visually.
- the CTV logic of the present invention can also be modified to reject anomalies such as, for example, unreliable or erroneous data Determination of Magnetic Susceptibility based on Cell Velocities.
- a cellular labeling complex generally contains a cell having a surface antigen or marker, e . g. , protem(s) , which serve as a marker to which a magnetic antibody can be attached.
- a surface antigen or marker e . g. , protem(s)
- Other cellular labeling complexes are also possible. For example, one may attach a fluorescent label which includes a primary antibody - fluorescein conjugate to a surface marker and a secondary antibody - magnetic label conjugate to the primary antibody.
- FACS Fluorescence-Activated Cell Sorting
- suitable labeling complexes include, for example, high or low density labels ⁇ e . g.
- the magnetic antibody may be of a plurality of types More particularly, magnetic antibodies can be classified into three broad categories which are based on size: Particulate labels, Colloidal magnetic labels and Molecular magnetic labels Particulate labels are the largest m relative size to the other labels and Molecular magnetic labels are the smallest m terms of relative size.
- cells may be rendered paramagnetic by binding a specific paramagnetic compound to a specific hapten on a cell or the specific or non-specific binding of a paramagnetic metal or metal complex directly to a cell (i.e., metal binding microorganism or by phagocytosis) . Therefore, it should be apparent to those skilled in the art that a cell may be rendered paramagnetic by a number of ways. While under the proper design specifications any of the three types of magnetic labels can be suitable, the present invention preferably employs either colloidal magnetic labels or molecular magnetic labels.
- a homogeneous cell population may be separated from other non- paramagnetic cell populations within a heterogeneous cell population.
- the separation of paramagnetic cells from non- paramagnetic cells is commonly referred to as binary separation.
- magnetic susceptibility the degree to which a paramagnetic label binds to a cell
- the opportunity to fractionally sort the paramagnetically labeled cell population into paramagnetic sub-populations has arisen. Accordingly, knowledge of the average magnetic susceptibility of a homogeneous cell population, along with the standard deviation and range, is required in the design of such fractional cell sorting devices. In particular, it can be shown that the magnetic susceptibility of a homogeneous cell population, along with the standard deviation and range, is required in the design of such fractional cell sorting devices. In particular, it can be shown that the magnetic susceptibility of a homogeneous cell population, along with the standard deviation and range, is required in the design of such fractional cell sorting devices. In particular, it can be shown
- Equation (2) Equation (2)
- TM permeability of free space
- V c velocity of the cell
- D c cell diameter
- a cell surface marker density
- ⁇ number of magnetic labels bound per cell surface receptor site
- B the magnetic flux density
- the CTV logic is capable of
- V determining the velocity of the cell V c y and the cell diameter L> c . Additionally, since the force field device 102 of the present invention provides a substantially uniform magnetic force within a
- Equation (2) is represented by a constant.
- the present invention may alternatively employ cells which have been labeled with differential density labels such as, for example, high density (e . g. , gold) and low density ⁇ e . g. , polymeric particles) labels.
- density labeling of cells modifies the cell motion in response to the various force fields ( e . g. , gravitational field) which can be employed by the present invention.
- electrically-charged labels such as, for example, ions, can also be employed by the present invention.
- electrically-charged labeling of cells also modifies the cell motion in response to the various force fields ⁇ e . g.
- a fluorescent label such as, for example, fluorochrome
- a fluorescent label such as, for example, fluorochrome
- Light sources which are compatible with fluorescent labels include, for example, ultraviolet and visible light sources.
- the fluorescent labels preferably emit or transmit colors m the range of green, yellow, or red. Once labeled, the cells can be classed or sub-classed based on fluorescence.
Abstract
Description
Claims
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU26607/99A AU2660799A (en) | 1998-02-06 | 1999-02-05 | System and device for determining particle characteristics |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/020,327 | 1998-02-06 | ||
US09/020,330 | 1998-02-06 | ||
US09/020,327 US6082205A (en) | 1998-02-06 | 1998-02-06 | System and device for determining particle characteristics |
US09/020,330 US5974901A (en) | 1998-02-06 | 1998-02-06 | Method for determining particle characteristics |
Publications (2)
Publication Number | Publication Date |
---|---|
WO1999040410A1 WO1999040410A1 (en) | 1999-08-12 |
WO1999040410A9 true WO1999040410A9 (en) | 1999-10-28 |
Family
ID=26693317
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US1999/002588 WO1999040410A1 (en) | 1998-02-06 | 1999-02-05 | System and device for determining particle characteristics |
Country Status (2)
Country | Link |
---|---|
AU (1) | AU2660799A (en) |
WO (1) | WO1999040410A1 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020167665A1 (en) * | 2000-05-19 | 2002-11-14 | Yeung Edward S. | High-throughput methods of distinguishing at least one molecule individually in a sample comprising multiple molecules and systems for use therein |
JP2007057360A (en) * | 2005-08-24 | 2007-03-08 | Agilent Technol Inc | Particle detector and particle detecting method used therein |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3434085A (en) * | 1967-05-08 | 1969-03-18 | Varian Associates | Magnets having logarithmic curved pole caps for producing uniform fields above saturation |
US4136950A (en) * | 1976-11-08 | 1979-01-30 | Labrum Engineering, Inc. | Microscope system for observing moving particles |
US5641919A (en) * | 1992-06-12 | 1997-06-24 | Dahneke; Barton E. | Method of characterizing particles by multiple time-of-flight measurements |
-
1999
- 1999-02-05 WO PCT/US1999/002588 patent/WO1999040410A1/en active Application Filing
- 1999-02-05 AU AU26607/99A patent/AU2660799A/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
WO1999040410A1 (en) | 1999-08-12 |
AU2660799A (en) | 1999-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6082205A (en) | System and device for determining particle characteristics | |
US5974901A (en) | Method for determining particle characteristics | |
AU2003219759B2 (en) | Methods and algorithms for cell enumeration in a low-cost cytometer | |
US7282180B2 (en) | Devices and methods to image objects | |
Tibbe et al. | Optical tracking and detection of immunomagnetically selected and aligned cells | |
AU2008249152B2 (en) | Devices and methods to image objects | |
EP2090889B1 (en) | Method and apparatus for imaging target components in a biological sample using permanent magnets | |
US20060094109A1 (en) | Device and method for analytical cell imaging | |
EP0780680A3 (en) | Analyzer for analyzing urine material components | |
Nakamura et al. | Theoretical and experimental analysis of the accuracy and reproducibility of cell tracking velocimetry | |
CA2484730A1 (en) | Device and method for analytical cell imaging | |
EP3913352A1 (en) | Microparticle isolation device, microparticle isolation system, droplet isolation device, droplet control device, and droplet control program | |
AU2003294351B2 (en) | Particle extraction for automatic flow microscope | |
WO1999040410A9 (en) | System and device for determining particle characteristics | |
Tárnok et al. | Immunophenotyping using a laser scanning cytometer | |
AU2002245215A1 (en) | Devices and methods to image objects | |
Varga | Detection of circulating tumor cells in samples from colorectal cancer patients using digital fluorescent microscopy and image cytometry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GE GH GM HR HU ID IL IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG UZ VN YU ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): GH GM KE LS MW SD SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
AK | Designated states |
Kind code of ref document: C2 Designated state(s): AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GE GH GM HR HU ID IL IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG UZ VN YU ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: C2 Designated state(s): GH GM KE LS MW SD SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG |
|
COP | Corrected version of pamphlet |
Free format text: PAGES 1-20, DESCRIPTION, REPLACED BY NEW PAGES 1-23; PAGES 21-29, CLAIMS, REPLACED BY NEW PAGES 24-33; DUE TO LATE TRANSMITTAL BY THE RECEIVING OFFICE |
|
DFPE | Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101) | ||
NENP | Non-entry into the national phase |
Ref country code: KR |
|
REG | Reference to national code |
Ref country code: DE Ref legal event code: 8642 |
|
122 | Ep: pct application non-entry in european phase |