CN1575475A - 用于差分细胞计数的方法以及用于执行该方法的相关装置和软件 - Google Patents

用于差分细胞计数的方法以及用于执行该方法的相关装置和软件 Download PDF

Info

Publication number
CN1575475A
CN1575475A CN02820990.7A CN02820990A CN1575475A CN 1575475 A CN1575475 A CN 1575475A CN 02820990 A CN02820990 A CN 02820990A CN 1575475 A CN1575475 A CN 1575475A
Authority
CN
China
Prior art keywords
survey data
cells
disc
data
points
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
CN02820990.7A
Other languages
English (en)
Inventor
麦克海尔·马特维夫
安德鲁·A·佩尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Burstein Technologies Inc
Original Assignee
Burstein Technologies Inc
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 Burstein Technologies Inc filed Critical Burstein Technologies Inc
Publication of CN1575475A publication Critical patent/CN1575475A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/5005Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells
    • G01N33/5094Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for blood cell populations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429Signal processing
    • G01N15/1433Signal processing using image recognition
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N35/00Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor
    • G01N35/00029Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor provided with flat sample substrates, e.g. slides
    • G01N35/00069Automatic analysis not limited to methods or materials provided for in any single one of groups G01N1/00 - G01N33/00; Handling materials therefor provided with flat sample substrates, e.g. slides whereby the sample substrate is of the bio-disk type, i.e. having the format of an optical disk
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N2015/1486Counting the particles

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Hematology (AREA)
  • Molecular Biology (AREA)
  • Cell Biology (AREA)
  • Urology & Nephrology (AREA)
  • Ecology (AREA)
  • Dispersion Chemistry (AREA)
  • Biotechnology (AREA)
  • Tropical Medicine & Parasitology (AREA)
  • Medicinal Chemistry (AREA)
  • Food Science & Technology (AREA)
  • Signal Processing (AREA)
  • Microbiology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Analysis (AREA)
  • Optical Measuring Cells (AREA)
  • Automatic Analysis And Handling Materials Therefor (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Image Processing (AREA)

Abstract

本发明提供了一种用于成像细胞特别是血细胞的光学方法、系统和软件。在一个实施例中,含有血细胞的实验室样品沉积在光生物盘上,该光生物盘是专门制造的光盘,具有混合室,其含有特效抗原用于锁定血细胞的各种组分。一旦放入光驱动器,光盘便会旋转,且样品和抗原与其它溶液混合。然后将电磁波束导向到光生物盘上,在特殊捕获区与样品相互作用,并通过探测器收集结果波束。然后包含在波束内的信息发送给处理器,由处理器产生数字图象。执行各种图象处理的方法,例如二值化、背景均匀化、标准化和过滤,以强化调查数据中用于精确计数的细胞。设计其它的技术用于修正不规则,例如气泡和模糊细胞。

Description

用于差分细胞计数的方法以及 用于执行该方法的相关装置和软件
相关专利申请参考
本发明要求获得如下专利申请的优先权,它们是美国临时专利申请序列No.60/322,863,其于2001年9月12日提出申请;美国临时专利申请序列No.60/355,300,其于2002年1月31日提出申请;美国临时专利申请序列No.60/353,921,其也是于2002年1月31日提出申请;美国临时专利申请序列No.60/355,644,其于2002年2月5日提出申请;美国临时专利申请序列No.60/355,304,其于2002年2月8日提出申请;美国临时专利申请序列No.60/358,479,其于2002年2月19日提出申请;美国临时专利申请序列No.60/63,949,其于2002年3月12日提出申请;和美国临时专利申请序列No.60/404,921,其于2002年8月21日提出申请。在此全文引用这些专利申请和公开作为参考。
关于已获得版权的材料的声明
本专利文件公开的一部分含有受到版权保护的材料。版权所有者不反对任何人按照其在专利商标局的文档或记录中出现的形式复制本专利文件和专利公开的副本,但保留其它所有的版权。
技术领域
本发明涉及信号探测装置、数据处理方法和相关计算机软件及化验算法(assay algorithms)。本发明更特别地致力于成像生物样品,例如细胞样品,并分析所收集的数据。更特别地,但不仅限于下文根据实践的最佳模式说明的特殊实施例,本发明涉及用于差分细胞计数(differential cell counting)的方法,包括白细胞,以及用于执行该细胞计数的光生物盘的使用。
背景技术
大量的研究和诊断情况都需要从细胞混合物中分离和分析特殊的细胞。这种混合物的来源可以包括血液、脊髓液、骨髓、肿瘤匀浆、淋巴组织和其它含有细胞物质的样品。
完全血细胞计数(complete blood count)(CBC)是检测包括血红蛋白、血细胞比容、平均血细胞血红蛋白、平均血细胞血红蛋白浓度、平均血细胞体积、血小板总数和白细胞总数的集合。最常用的临床检测是总CBC计数,其常规地用于评估健康状况和临床诊断、治疗及跟踪调查(follow-up)。
白细胞(WBC)通过抵抗感染和攻击外来物质(foreign material)保护肌体。差分白细胞计数确定白细胞的数量和人血中每种类型白细胞的百分比。WBC或白细胞计数提供了疾病存在的线索。这些检测包含在通常的健康检测之中,并有助于调查各种疾病,包括感染、过敏症和白血病。在需要额外的白细胞时,骨髓将增加产量。
有五种类型的白细胞,每一种具有不同的功能:嗜中性粒细胞、淋巴细胞、单核细胞、嗜伊红细胞和嗜碱细胞。差分显示出这些细胞是否以正常的分布存在,或者是否有一种细胞类型增加了或减少了。在正常健康人中,典型的WBC总数为每微升(μl)4,000-10,800个细胞。如下因素,例如锻炼、应激(stress)和疾病等能够影响这些数值。这一信息有助于诊断特殊类型的疾病。高WBC可能表明感染、白血病或组织损伤。如果降到每微升低于1,000个细胞,则感染的危险增加。减弱免疫系统的疾病(condition)和药物,例如AIDS或化疗,会导致白细胞减少。从疾病中恢复可以通过白细胞加以监测。总数持续增加或降低到异常的水平表示疾病正在恶化;总数返回到正常表示状况改善。
白细胞差分检测(differential testing)是收集超出可以从白细胞计数自身获得的信息之外的信息不可或缺的。白细胞差分计数用于评估新近的疑似感染或感冒(如果CBC正常)、伴随异常的疑似紊乱、白细胞总数异常、疑似白血病和其它异常,例如嗜酸粒细胞增多、单核细胞增多、嗜碱细胞增多。白细胞或白细胞差分的重复测试可以在出现严重白血球减少症的情况下进行(例如在药物治疗之后)。在治疗期间,例如化疗或放疗,血液计数对于确定治疗是否除了消灭了癌细胞之外也消灭了正常的血细胞是非常重要的。
差分白细胞计数通过计算机细胞计数设备加以确定。机器确定五种主要白细胞类型的总数和百分比。在正常个体中,大多数是嗜中性粒细胞(50-60%),其次是淋巴细胞(20-40%),然后是单核细胞(2-9%),有少量的嗜伊红细胞(1-4%)和嗜碱细胞(0.5-2%)。
在淋巴细胞分类中,有进一步的细胞亚型。例如,淋巴细胞能够广义地分为T细胞(胸腺产生的淋巴细胞)和B细胞(囊等同物淋巴细胞(bursal-equivalent lymphocyte)),它们分别大量响应细胞介导的免疫和体液免疫。尽管形态学特征已经用于划分白细胞的类群,但是已经证明,形态学自身不足以区分淋巴细胞亚型的许多功能。为了区分具有各种功能的淋巴细胞,已经开发出了许多技术,包括玫瑰花结分析(analysis by rosetting)、免疫荧光显微镜、酶组织化学和最近的抗独特细胞表面标记物的单克隆抗体。
嗜中性粒细胞对于抗击感染很重要。当嗜中性粒细胞数目下降到低于每微升1,000个细胞时,该状况称为嗜中性白血球减少症。淋巴瘤治疗能够导致嗜中性白血球减少症。肥胖和吸烟会增加嗜中性粒细胞总数。淋巴细胞分为B淋巴细胞(骨髓中化脓)和T淋巴细胞(胸腺中化脓)。当成年人的淋巴细胞总数下降到低于每微升1,500个细胞,或者在儿童体内低于每微升3,000个细胞时,该状况称为淋巴细胞减少症。淋巴瘤能够导致淋巴细胞减少症。
血小板(凝血细胞)是类似细胞的颗粒,其通过聚集在发生出血的部位而停止出血。然后,它们激活并凝集在一起,以停止出血和促进凝血。如果患者患有骨髓增殖紊乱(myleoproliferative disorder),包括感染、炎症、恶性肿瘤,以及如果切除了脾,则在剧烈激活期间,血小板数量增加。血小板数量过多称为血小板增多症。
标准血液样品内血小板的数量典型地为每微升(μl)133,000-333,000个血小板。血小板数量过高称作血小板增多症(thrombocythemia)。高于正常的血小板总数可能是由于反应性响应(reactive response)或骨髓衰竭(bone marrow failure)。反应性响应典型地由出血、感染、肿瘤形成和骨髓增生紊乱(myleoproliferative disorder)导致。骨髓衰竭通常涉及血细胞损失,称作全血细胞减少症。另一方面,血小板总数的减少是因为免疫血小板减少症。如果血小板总数下降到低于30,000则发生血小板减少症,其会导致异常出血。总数低于5,000则认为有生命危险。
CBC可以通过商业上可获得的手动或电子仪器执行,它们测量血红蛋白水平、血细胞比容、总白细胞和红细胞总数。变化可以包括血小板计数、白细胞差分计数和细胞指数(cell indices)。血液分析仪完全自动化,且细胞计数、体液中的细胞类型和胃排空(gastricaspiration)的结果精确,其中体液包括例如CSF、胸腔液、禁欲液(ascetic fluid)、心包液。
与先前的方法和系统相比,我们开发出了一种用于对成像和分析细胞及其成分的简单、微型、极其灵敏、低廉的系统。该系统使用光学光生物盘、相关探测组件以及信息和信号的处理方法和软件。
发明内容
本发明涉及用于成像和计数实验室样品中细胞物质的方法、装置和软件。本发明的实施例生成样品中细胞的数字图象,并对图象进行计算机分析。本发明成像细胞,特别是血细胞,包括在血液和其它生物液中滋生的寄生虫和病原体。在其它化验中,对珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记(biologicalreporters)进行成像。该系统使用光学光生物盘、相关探测组件,以及信息和信号处理的方法和软件。
本发明还致力于光生物盘、生物驱动和相关方法。本发明或其不同的方面可以容易地在如下的光盘、化验和系统中实现,或者能够容易地适用于或与它们结合使用,这些光盘、化验和系统在如下共同转让的并且同在审理中(commonly assigned and co-pending)的专利申请中公开,这些专利申请包括:美国专利申请序列No.09/378,878,题目为“用于分析从光盘获得的可操作和非可操作数据的方法和装置”(Methods and Apparatus for Analyzing Operational andNon-operational Data Acquired from Optical Discs),其于1999年8月23日提出申请;美国临时专利申请序列No.60/150,288,题目为“用于使用物理同步化标记物获取光盘数据的方法和装置”(Methods andApparatus for Optical Disc Data Acquisition Using PhysicalSynchronization Markers),其于1999年8月23日提出申请;美国专利申请序列No.09/421,870,题目为“具有并发可读分析物材料的可跟踪光盘”(Trackable Optical Discs with Concurrently ReadableAnalyte Material),其于1999年10月26日提出申请;美国专利申请序列No.09/643,106,题目为“用于使用物理同步化标记物获取光盘数据的方法和装置”(Methods and Apparatus for Optical Disc DataAcquisition Using Physical Synchronization Markers),其于2000年8月21日提出申请;美国专利申请序列No.09/999,274,题目为“具有反射层的光生物盘”(Optical Biodiscs with Reflective Layers),其于2001年11月15日提出申请;美国专利申请序列No.09/988,728,题目为“使用光生物盘探测和定量淋巴细胞的方法和装置”(Methods andApparatus for Detecting and Quantifying Lymphocytes with OpticalBiodiscs),其于2001年11月20日提出申请;美国专利申请序列No.09/988,850,题目为“使用光生物盘进行血液分类的方法和装置”(Methods and Apparatus for Blood Typing with Optical Bio-discs),其于2001年11月19日提出申请;美国专利申请序列No.09/989,684,题目为“用于分离粘聚物和分散颗粒的装置和方法”(Apparatus andMethods for Separating Agglutinants and Disperse Particles),其于2001年11月20日提出申请;美国专利申请序列No.09/977,741,题目为“包括光生物盘的双珠化验及其相关方法”(Dual Bead AssaysIncluding Optical Biodiscs and Methods Relating Thereto),其于2001年11月27日提出申请;美国专利申请序列No.09/977,895,题目为“用于分离特殊悬浊液组分的装置和方法”(Apparatus and Methods forSeparating Component of Particulate Suspension),其于2001年11月30日提出申请;美国专利申请序列No.10/005,313,题目为“用于测量分析物的光盘”(Optical Discs for Measuring Analytes),其于2001年12月7日提出申请;美国专利申请序列No.10/006,371,题目为“使用光盘和光盘阅读器探测分析物的方法”(Methods for DetectingAnalytes Using Optical Discs and Optical Disc Readers),其于2001年12月10日提出申请;美国专利申请序列No.10/006,620,题目为“用于探测分析物的多数据层光盘”(Multiple Data Layer Optical Discsfor Detecting Analytes),其于2001年12月10日提出申请;美国专利申请序列No.10/006,619,题目为“用于执行化验的光盘组件”(Optical Disc Assemblies for Performing Assays),其于2001年12月10日提出申请;美国专利申请序列No.10/020,140,题目为“用于基于光盘的实验室的探测系统和包括该系统的改良光生物盘”(Detection System For Disk-Based Laboratory and ImprovedOptical Bio-Disc Including Same),其于2001年12月14日提出申请;美国专利申请序列No.10/035,836,题目为“用于固定DNA捕获探针的表面组件和包括光生物盘的基于珠子的化验及其相关方法”(SurfaceAssembly for Immobilizing DNA Capture Probes and Bead-BasedAssay Including Optical Bio-Discs and Methods Relating Thereto),其于2001年12月21日提出申请;美国专利申请序列No.10/038,297,题目为“包括用于改进特异性的共价键合的双珠化验及相关分析光盘”(Dual Bead Assays Including Covalent Linkages for ImprovedSpecificity and Related Optical Analysis Discs),其于2002年1月4日提出申请;美国专利申请序列No.10/043,688,题目为“光盘分析系统及用于生物学和医学成像的相关方法”(Optical Disc AnalysisSystem Including Related Methods for Biological and MedicalImaging),其于2002年1月10日提出申请;美国临时专利申请序列No.60/348,767,题目为“光盘分析系统及相关信号处理的方法和软件”(Optical Disc Analysis System Including Related Signal ProcessingMethods and Software),其于2002年1月14日提出申请;美国专利申请序列No.10/086,941,题目为“用于将DNA结合在固相上的方法及相关光生物盘和光驱系统”(Methods for DNA Conjugation OntoSolid Phase Including Related Optical Biodiscs and Disc DriveSystems),其于2002年2月26日提出申请;美国专利申请序列No.10/087,549,题目为“用于减少双珠化验中珠子非特异键合的方法及相关光生物盘和光驱系统”(Methods for Decreasing Non-SpecificBinding of Beads in Dual Bead Assays Including Related OpticalBiodiscs and Disc Drive Systems),其于2002年2月28日提出申请;美国专利申请序列No.10/099,256,题目为“使用可断开连接臂和/或结扎以提高特异性和敏感性的双珠化验及相关方法和装置”(Dual BeadAssays Using Cleavable Spacers and/or Ligation to ImproveSpecificity and Sensitivity Including Related Methods andApparatus),其于2002年3月14日提出申请;美国专利申请序列No.10/099,266,题目为“使用限制性酶和其它化学方法减少双珠化验中的非特异键合及用于探测药物靶标的相关光生物盘、方法和系统装置”(Use of Restriction Enzymes and Other Chemical Methods toDecrease Non-Specific Binding in Dual Bead Assays and RelatedBio-Discs,Methods,and System Apparatus for Detecting MedicalTargets),其也于2002年3月14日提出申请;美国专利申请序列No.10/121,281,题目为“多参数化验及其相关分析光盘和方法”(Multi-Parameter Assays Including Analysis Discs and MethodsRelating Thereto),其于2002年4月11日提出申请;美国专利申请序列No.10/150,575,题目为“用于在光生物盘组件内提供分析结果像素的可变取样控制及其相关装置”(Variable Sampling Control forRendering Pixelization of Analysis Results in a Bio-Disc Assembly andApparatus Relating Thereto),其于2002年5月16日提出申请;美国专利申请序列No.10/150,702,题目为“在使用酶反应在光生物盘内产生信号的基因化验中用于固定DNA探针的表面组件及其相关方法”(Surface Assembly for Immobilizing DNA Capture probes in GeneticAssays Using Enzymatic Reactions to Generate Signals in OpticalBio-Discs and Methods Relating Thereto),其于2002年5月17日提出申请;美国专利申请序列No.10/194,418,题目为“用于分析显微结构的光盘系统及相关探测和解码方法”(Optical Disc System andRelated Detecting and Decoding Methods for Analysis of MicroscopicStructures),其于2002年7月12日提出申请;美国专利申请序列No.10/194,396,题目为“用于指导化验的多目的分析光盘及用于使用该光盘的各种指示剂”(Multi-Purpose Optical Analysis Disc forConducting Assays and Various Reporting Agents for UseTherewith),其于2002年7月12日提出申请;美国专利申请序列No.10/199,973,题目为“用于执行体质测量的可传送光盘组件及其相关方法”(Transmissive Optical Disc Assemblies for Performing PhysicalMeasurements and Methods Relating Thereto),其于2002年7月19日提出申请;美国专利申请序列No.10/201,591,题目为“用于执行反应离心的分析光盘及相关驱动组件”(Optical Analysis Disc andRelated Drive Assembly for Performing Interactive Centrifugation),其于2002年7月22日提出申请;美国专利申请序列No.10/205,011,题目为“用于光生物盘键合流体线路的方法和装置”(Methods andApparatus for Bonded Fluidic Circuit for Optical Bio-Disc),其于2002年7月24日提出申请;美国专利申请序列No.10/205,005,题目为“使用光盘驱动的磁珠磁辅助探测”(Magnetic Assisted Detection ofMagnetic Beads Using Optical Disc Drives),其也于2002年7月24日提出申请。本文全文引用所有这些专利申请作为参考。因此,它们提供了背景和相关公开,作为本文所重复内容的支持。
本发明用于执行化验的一个方法是基于位于在光生物盘特殊沟道内的血细胞的光学成像原理。大约数微升的全血注射到光盘上特殊设计的沟道内。使用细胞识别软件分析图象,辨认各种白细胞亚型并生成白细胞差分计数。该方法以使用抗特定细胞的细胞特异抗体捕获特定细胞为基础。在特殊实例中,抗体定向地抗淋巴细胞(CD2、CD19)、单核细胞(CD14)、嗜伊红细胞(CD15)等。这些白细胞亚型特异抗体组装并附着在光生物盘内的固体表面上,该光生物盘包括流动室(flow chamber)。
光生物盘驱动组件用于旋转光盘,阅读和处理任何存储在光盘内的编码信息,并分析光生物盘流动室内的细胞捕获区。光生物盘驱动安装有用于旋转光生物盘的马达,用于控制光盘旋转速度的控制器,用于处理从光盘返回的信号的处理器,和用于分析已处理信号的分析器。旋转速度可变,并可以密切地控制速度、方向和旋转时间。光生物盘也可以用于将信息写入到光生物盘上,其可以发生在利用驱动的阅读波束询问流动室和靶区内的检测材料之前、之中或之后,并用分析器加以分析。光生物盘可以含有用于控制光盘旋转的编码信息,提供特定于(specific to)待执行免疫分类化验类型的处理信息,和用于在与生物驱动关联的监视器上显示任何期望的结果。
差分细胞计数方案(protocol),通常上且特别地,差分白细胞计数方案是为CD、CD-R或DVD格式、这些格式的修改版本及其替代物开发的。驱动的阅读或询问波束探测分析样品中的各种细胞,并产生能够用差分细胞计数器软件加以分析的图象。
显微方法或精致的细胞计数器对于执行这些费事的细胞计数化验是至关重要的。在其它根据本发明进行的化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物质、沉淀(酶反应),或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。
本方法使用光生物盘和相关的光盘组件。通过细胞识别软件程序产生和分析在分析室内游离的或者通过特异抗体方法捕获的各种白细胞亚型的光学图象,该程序通过它们的光散射性鉴别血液或其它体液中的各种细胞基本组分(cellular element)。本发明在分析之前不需要对样品进行任何的处理,例如细胞染色、RBC清除和其它艰苦的规程。这些方法包括在CD型阅读器、DVD型阅读器或其它光盘阅读器中使用顶部探测器(top detector)、底部探测器(bottom detector)、事件计数器(event counter)或细胞计数器(cell counter)进行显微分析或细胞探测。
下述关于聚类指示分析(cluster designation analysis),例如获得CD4/CD8比率代表了能够应用本方法、装置、系统和光盘系统的相关化验的特定组群。
光盘制备:金反射光盘或透射光盘用空气枪(air gun)清洁用以除去任何的灰尘颗粒。使用旋转涂布机(spin coater)用异丙二醇漂洗光盘两次。将2%聚苯乙烯旋转涂布在光盘上从而提供相对厚的涂布(coating throughout)。
化学沉积:一个实施例包括三步培育(incubate)沉积方案:链霉素(streptavidin),培育30分钟;生物素标记的(biotinylated)第一抗体,培育60分钟;和第二捕获抗体,培育30分钟。所有的步骤都适宜在潮湿箱(humidity chamber)内在室温下进行,在沉积之间进行严格的冲洗和干燥步骤。
简而言之,磷酸盐缓冲液中1μl 1mg/ml的链霉素层积(layered)在每个窗口上,并培育30分钟。用蒸馏水冲洗掉多余的链霉素,并干燥光盘。通过合并等体积的生物素标记的lgG(125μg/ml置于PBS中)和乙醛激活的葡聚糖(200μg/ml)制备生物素标记的lgG-葡聚糖复合物。葡聚糖-乙醛生物素标记-lgG复合物层积在每个捕获窗口中的链霉素上,并培育60分钟或者在冰箱中过夜。冲洗掉多余的试剂并旋转干燥光盘。通过在光生物盘狭槽上的指定点上层积捕获抗体产生特殊的条码捕获模式。为了差分计数,抗嗜中性粒细胞(CD128或其它)、抗淋巴细胞(CD2、CD19、CD56及其它)、抗嗜伊红细胞(CD15)、抗单核细胞(CD14)、抗嗜碱细胞(CD63)和抗血小板(CD32和CD151)的抗体层积在每个狭槽的指定点内。下面的表1列举了用于捕获层组件的各种捕获模式的实例。培育30分钟或在冰箱中过夜。组装光盘使用25μm、50μm或100μm(50μm沟道需要的样品体积是25μm室所需要的2倍)、直的、U形的或其它的沟道格式,以及透明(用于顶部探测器)或反射(用于底部探测器)顶盖光盘(cover disc)。
表1:捕获层组件及变型
  窗口   1     2     3     4     5     6
第一层(激活层) 聚苯乙烯 聚苯乙烯 聚苯乙烯 聚苯乙烯 聚苯乙烯 聚苯乙烯
第二层 链霉素 链霉素 链霉素 链霉素 链霉素
第二抗体 B抗鼠lgG+DCHO B抗鼠lgG+DCHO B抗鼠lgG+DCHO B抗鼠lgG+DCHO B抗鼠lgG+DCHO
第一抗体 参考点 淋巴细胞特异抗体 嗜中性粒细胞特异抗体 嗜伊红细胞特异抗体 嗜碱细胞特异抗体 单核细胞特异抗体
光盘泄漏检查:因为典型分析的是血液,一种生物危害性物质,所以这些光盘要进行泄漏检查以保证含有样品的光盘在原位旋转期间,没有室发生泄漏。每个沟道用阻塞剂(blocking agent)StabilGufard填充,并阻塞一小时。光盘以5,000rpm旋转5分钟,并检查泄漏和光盘稳定性。检查泄漏之后,光盘在真空室内放置24小时。抽真空之后,充满磷酸缓冲盐(PBS)缓冲液的,或者选择地,空的室放置在真空袋内,并冷冻保存直至使用。
从全血中分离血沉棕黄层(buffy coat):通过在1,500xg下离心离心管内的静脉血15分钟制备血沉棕黄层,该静脉血含有抗凝剂,如乙二胺四乙酸(EDTA)或酸性柠檬酸葡聚糖(acid citrate dextran)(ACD)。白细胞在血浆与红细胞之间形成层,称为血沉棕黄层。用精确的移液管小心地除去血浆,然后收集血沉棕黄层。一种无离心地从血液中获得血沉棕黄层的可选择方法是允许血液与沉淀作用增强剂,例如纤维蛋白原、葡聚糖、阿拉伯树胶、Ficoll或甲基纤维素,一起沉淀。Boyum试剂(甲基纤维素和甲泛影酸钠)特别适合于获得没有任何红细胞污染的白细胞制品。
光盘上化验——基本技术说明:差分白细胞计数光盘检测的一个优选实施例包括三个单独的部件,(1)含有化学试剂的基础光盘(basedisc),(2)沟道层,和(3)顶盖光盘(cover disc)。
血沉棕黄层或白细胞(7微升置于PBS中)注射到光盘室内,室的入口和出口部分用封闭标签密封,并且光盘在室温下培育15分钟。对于第一方法,光盘上的给定区域(例如1平方毫米的面积)用具有顶部或底部探测器的光驱动器的标准780nm激光加以扫描。根据本发明的细胞识别软件自动地从1平方毫米的捕获图象中给出差分总数,并且外推所获得的值确定每毫升全血的总数。对于第二条形码方法,用标准780nm激光扫描光盘以成像捕获区(淋巴细胞、嗜中性粒细胞、嗜碱细胞、嗜伊红细胞、单核细胞和血小板)。本发明的细胞识别软件除其他事务之外执行如下程序:(a)离心光盘以旋出多余的未键合细胞,(b)成像每个特定细胞捕获区内的限定区域,(c)处理数据,其包括计数每个捕获区中特定的捕获细胞,和(d)获得每毫升全血中不同亚型白细胞的数量。
根据本发明的一个方面,在处理步骤期间,识别软件通读每一个捕获区,并在遇到细胞的时候标记细胞。在其它化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。随后处理来自每个捕获区的数据,软件显示每微升或毫升血液体积中淋巴细胞、嗜中性粒细胞、嗜碱细胞、嗜伊红细胞、单核细胞和血小板区的数量。从将光盘插入光驱动器到获得并显示期望的总数或比率,整个处理花费大约10-15分钟。在本发明这一方面的另一个实施例中,从捕获区读出电响应(electrical response),并存储在光盘上或存储器内,形成可以进行后处理的数据文件以实现在下文中更详细说明的识别目的。
光盘详细说明:下面的分部(subsection)用于概述某些光生物盘的特殊实施例,其可以有利地与本发明一起使用。
(A)跟踪设计(tracking design):在本发明的一个优选实施例中,光盘是涂布有300nm的金的前摆动序列(forward Wobble Set)FDL21:13707或FDL21:1270。在该反射光盘上,通过平版印刷术(lithography)腐蚀出尺寸为2×1mm的椭圆形数据窗口。U-形沟道用于产生高度为25-100微米的室。充满整个室,包括入口和出口部分,需要大约7μl的样品。可以优选地使用4窗口/4沟道格式。然而在透射光盘上不腐蚀数据窗口,并且整个光盘都可以使用。
(B)粘着与键合:使用磨损锁定U-形粘合剂(Fraylock U-shapedadhesive)DBL 201 Rev C 3M94661或直沟道形成室。
(C)顶盖光盘:使用干净光盘(clear disc),其完全反射,具有48个直径为0.040英寸的在半径26mm上等距离放置的样品入口。
数据捕获与处理:用本发明的软件使用特定细胞识别方法以x4的速度和2.67MHz的取样速度扫描和阅读数据光盘。
软件:本发明进一步包括处理方法和相关细胞识别和成像的软件。该软件致力于实施和显示细胞计数和差分细胞计数。在其它化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。本软件可以存储在光盘驱动阅读器中的光生物盘上,或者选择地,只能够通过可靠服务器上的光学阅读器加以访问。该服务器可以在计算网络,例如局域网(LAN)、广域网(WAN),中实现,或者可以在指定的期限和条件下在整个互联网上都能够获得。该分配方法在如下专利申请中公开,即共同转让的美国临时专利申请No.60/246,824,题目为“用于分析生物样品的交互方法和系统和使用专门制备的光生物盘、光盘驱动和网络连接处理相关的医学信息”,其于2000年11月8日提出申请,并对应于美国专利申请序列No.09/986,078。
用于实践本文所公开的不同优选实施例的材料包括前摆动金喷涂光致抗蚀剂光盘(forward wobble gold metalized photo-resist disc),a_透射金喷涂光盘(a_transmissive gold metalized disc),移液管和吸头(tip),旋转涂布器、离心机、回旋转子(swing-out rotor),含有抗凝剂,例如柠檬酸钠或乙二胺四乙酸(EDTA),的VacutainerTMCPT管,潮湿室,绞干机(wringer),粘合剂,顶盖光盘,干净顶盖光盘,磁带或等价物,真空装置,黄吸头和真空室。
在本发明的一个实施例中,含有血细胞的实验室样品沉积在光生物盘上或者沉积在形成于光盘组件内的流体沟道内。在其它化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。光生物盘是专门制造的CD尺寸的光盘,具有流体沟道和/或混合室,其含有特效的用于原位锁定血细胞成分的抗原。因为它们用仔细层积(layered)的金属和物质制成,所以光生物盘具有特殊的光学性能,其允许电磁波束与沉积的测试样品相互作用。一旦光生物盘插入到光驱动器内,驱动便旋转光盘,并且在一些实施例中可以将样品和其它所需溶液混合在一起。电磁波束被导向到驱动内的光生物盘上。在一个称作反射光盘的实施例中,波束从光生物盘的反射表面上反射,且光驱动器内的探测器收集反射的波束。在另一个称作透射光盘的实施例中,部分波束通过光生物盘,并透射到光驱动器内的另一类型探测器上。在上述任一实例中,探测器收集的波束都含有关于光生物盘上实验室样品的信息。然后将该信息发送到模拟-数字处理器,在该处产生表示来自探测器的电信号的数字数据。该数字数据可以实时地加以处理,存储在存储器内或光盘上,然后整体地或部分地作为粗数据加以处理,或者输出成各种格式,包括图象格式。这些格式的任何一种都可以进一步通过其它装置或设备加以处理从而生成预期的结果。该数字数据对于自动的、电的、计算机控制的,和/或机械的计数或分析是有用的。该数字数据也可以用于产生适合于熟练手工计数(experthand counts)、识别或其它手动分析的可视图象。在本发明的另一个实施例中,粗数据、数字数据、可能含有图象的输出数据存储在存档(archive)中。这样,本发明的方法通常可以应用于如本文所使用的“调查数据”(investigational data),其包括但不仅限于,粗探测器输出数据、粗信号数据、数字数据、输出数据、或含有图象或图象数据的输出数据。存档提供一个位置,在该处调查数据能够被分类,并且如果期望的话,与其它辨识信息,诸如例如人口统计学、地理学、医学、历史学或个人数据相结合。随后,能够分析调查数据的组从而进行例如不同人群的健康趋势研究。
本发明的一个实施例致力于计数血细胞的需要。本发明包括处理方法和相关细胞识别和成像的软件。该软件致力于执行细胞计数和显示相应的结果。在本发明的一个实施例中,执行各种图象处理方法,例如二值化(binarization)、背景均匀化、标准化(normalization)和过滤,以强化调查数据(investigational data)中细胞的显现,从而辅助细胞计数的处理。执行其它的技术以修正细胞计数的不规则,例如调查数据中的截留气泡、裂缝和模糊细胞(dim cell)。
本发明的实施例将软件存储在光生物盘上、光盘驱动阅读设备内,或者选择地,仅能通过可靠服务器的光学阅读器加以访问。该服务器可以在计算机网络中实现,例如局域网(LAN)、广域网(WAN)或者整个互联网在特定期限和条件下都可以获得。该分配方法在如下专利申请中公开,即共同转让的美国专利申请No.09/898,078,题目为“用于分析生物样品和处理相关信息的交互系统及其使用”,其于2001年11月7日提出申请,在这里引用作为参考。
更特别地,本发明涉及计数细胞或其它调查特征(investigationalfeature)的方法。该方法包括如下步骤:获得含有细胞的样品的调查数据,在调查数据中选择估计矩形(evaluation rectangle),强化估计矩形内的调查数据,和计数估计矩形中显示的细胞。在本发明的一个特殊实施例中,细胞计数通过识别明亮中心或者选择地黑暗边缘加以执行。
本发明的另一个方面致力于为估计矩形选择自定义尺寸(customsize)的方法。
本发明的再一个方面致力于选择多个估计矩形的方法。
本发明进一步的方面致力于通过如下步骤强化估计矩形内的调查数据,即对调查数据执行背景照明均匀化,对调查数据执行标准化,和过滤调查数据。
本发明附加的方面致力于通过如下步骤对调查数据执行背景照明均匀化,即为相邻矩形(neighborhood rectangle)选择一个尺寸,在调查数据中选取一个点,执行水平扫描以便为定位在将中心定于该点的相邻矩形内的所有相邻点计算第一滑动平均值(sliding average),执行垂直扫描以便为定位在将中心定于该点的相邻矩形内的所有相邻点计算第二滑动平均值,组合第一滑动平均值与第二滑动平均值产生整体平均值,将点的原始值再分配为通过获得整体平均值与原始值之间的差值并将该差值加和到背景值上计算出来的结果值(resultantvalue),和重复执行水平扫描、执行垂直扫描、组合两个平均值和为调查数据内的所有点再分配原始值的步骤。
在本发明的另一个方面中,对调查数据执行标准化的步骤进一步包括如下步骤:计算调查数据中所有点的数值的平均值和标准偏差,使用该平均值和标准偏差使调查数据中所有点的数值标准化,和如果需要的话,截去一些点的数值。
根据本发明的另一个方面,提供了一种过滤调查数据的方法,其包括如下步骤:为相邻矩形选择尺寸,在调查数据中选取一个点,发现所有定位在将中心定于该点的相邻矩形内的足够区别的点,如果足够区别点的数量大于预先确定的标准则再分配点的值,和重复发现所有足够区别的点和为调查数据内的所有点再分配数值的步骤。
根据本发明更进一步的方面,提供了一种处理方法,包括如下步骤:从调查数据中除去不合需要的成分,该步骤包括如下步骤:选择阈值,使用阈值对调查数据执行二值化,对调查数据执行规则化(regularization),提取关联成分(connected components),选择尺寸阈值(size threshold)和除去不能满足尺寸阈值的成分。
本发明的另一个附加方面致力于通过明亮中心计数调查数据内的细胞的方法。该方法包括如下步骤:对调查数据进行卷积,寻找多个局部最大值,从多个局部最大值中除去冗余的局部最大值,断言(declaring)剩余的最大值为细胞中心,和计数细胞中心。
根据本发明另一个方面,提供了另一种通过明亮中心计数调查数据中的细胞的方法。这种可选择的方法包括对所述调查数据进行反转(inversion),进行多次使用移动环(shifted rings)的卷积,总计所述多次卷积的结果,发现局部最大值,断言局部最大值为细胞中心,和计数所述细胞中心。
进一步,本发明的另一个实施例致力于一种方法,其包括如下步骤:从已经通过明亮中心加以计数的调查数据中除去已计数的细胞,通过它们的黑暗边缘计数细胞,和将通过识别明亮中心计数细胞的步骤的总数加和于通过识别黑暗边缘计数细胞的步骤的总数。
本发明的再一个方面包括通过调查数据中的黑暗边缘计数细胞的方法。该方法包括如下步骤:对调查数据进行反转,执行多次使用移动环的卷积,总计该多次卷积的结果,断言最大值是细胞中心,和计数细胞。
在本发明的另一个实施例中,出于细胞计数的目的强化调查数据的方法进一步包括如下步骤:对调查数据进行标准化,对调查数据进行过滤,选择阈值数量,通过用该阈值数量确定调查数据是否与已设定的背景值不同对调查数据进行二值化,对调查数据进行规则化,在调查数据中提取一个像素宽度的边界,填充由一个像素边界限定的区域,和在所填充区域内进行卷积。
本发明的另一个方面涉及获得含细胞样品数字数据的方法。该方法包括如下步骤:(1)在光盘表面提供血液样品,(2)将光盘装载到光学阅读器内,(3)旋转光盘,和(4)将电磁辐射的入射波束导向到光盘上的一个捕获区。该表面上设置有具有一种或多种捕获剂的一个或多个捕获区。该方法接着包括步骤(5)用探测器(detector)探测与光盘在捕获区相互作用之后形成的电磁辐射波束,(6)将探测到的波束转换为模拟输出信号,和(7)将模拟输出信号转换为含有在捕获区捕获到的细胞的数字数据。
根据本发明的另一个方面提供了一种将模拟输出转换为数字数据的方法。该转换方法包括如下步骤:以固定的间隔采样模拟信号的振幅,将采样振幅记录在一维阵列内,利用采样和记录步骤产生多个一维阵列,和组合该多个一维阵列生成含有样品的数字数据的二维阵列。
本发明再一个方面涉及获得含细胞样品的数字数据的方法。该方法包括如下步骤:在光盘表面上提供血液样品(该表面包括含有一种或多种捕获剂的一个或多个捕获区),将光盘装载到光学阅读器内,旋转光盘,将电磁辐射的入射波束导向到光盘的一个捕获区上,用探测器探测与光盘在捕获区相互作用之后形成的电磁辐射波束,和将探测到的波束转换为模拟输出信号。本发明的该特殊实施例包括如下步骤:将模拟输出信号转换成含有在捕获区内捕获的细胞的数字数据。该光盘用反射层构建,从而导向到捕获区的光线被反射到探测器,且探测器是底部探测器。在本发明的另一个方面中,使用顶部探测器或分波探测器(split detector)。
本发明的另一个方面致力于选择估计矩形的方法。该方法包括如下步骤:(1)在调查数据中发现多个窗口中的一个,和(2)在窗口内修剪(cropping)标准尺寸的估计矩形,该步骤包括(a)对调查数据进行压缩,(b)对图象进行阈值估计,(c)对调查数据进行二值化,(d)对调查数据进行规则化,(e)从调查数据中提取关联成分,和(f)从相应于窗口的关联成分中发现成分。
本发明再进一步的方面致力于从调查数据中提取关联成分的方法。该方法包括如下步骤:向调查数据上所有黑色点(black point)的成分分配初始成分数目,设定初始扫描方向,和扫描调查数据重新分配成分数目,从而使关联黑色点的成分数目相同。
且本发明再附加的一个方面致力于在从使用暗点的光盘实施例中获得的调查数据中选择估计矩形的方法。该方法包括如下步骤:在调查数据中发现至少一个暗点,和生成标准尺寸的估计矩形,其中心定位在通过从暗点移动预先确定的距离所发现的点处。
根据本发明的显示方案,提供了一种强化调查数据图象显示的方法。该方法包括如下步骤:对调查数据进行快速傅立叶变换,除去频域(frequency domain)中数据频谱的一些部分,和执行逆变换以恢复调查数据的已修改版本。
同样根据本发明的显示方案,提供了另一种强化调查数据图象显示的方法。该方法包括如下步骤:确定图象是否歪斜(skewed),发现歪斜方向,和修正图象的歪斜。
本发明的另一个方面涉及从存档中获取先前存储的调查数据并对调查数据进行分析的方法。该存档能够根据提供检验样品的患者的特征对存储的调查数据进行分类。在本发明的一个方面中,符合从患者特征中选择出的多个标准的样品被选择进行人口健康趋势研究。
本发明的另一个方面涉及在白细胞计数中计数不同的组分,并显示CD4+细胞和CD8+细胞的数量和CD4与CD8细胞的比率。
附图说明
本发明进一步的目的、方面和方法以及有助于它的附加特点及由其产生的优点,从下文对附图所显示的本发明优选实施例的说明中将显而易见,其中:
图1是根据本发明的光生物盘系统的图示;
图2是结合本发明使用的反射光生物盘的分解透视图;
图3是图2所示光盘的顶视平面图。
图4是图2和3所示出的光盘的透视图,其具有显示光盘不同层的切除部分;
图5是结合本发明使用的透射光生物盘的分解透视图;
图6是图5所示光盘的顶视平面图;
图7是图5和6所示出的透射光盘的透视图,其具有显示光盘不同层,包括图8所示半反射层类型,的切除部分;
图8是描绘图7所示光盘的透视图,其具有示出光盘半反射层功能方面(functional aspect)的切除部分;
图9是显示薄金膜厚度与透射性之间关系的图示;
图10A是示出根据本发明一个实施例的系统操作的透视图和框图;
图10B显示了分波探测器(split detector)和根据本发明一个实施例的透射光生物盘的横截面;
图11是垂直于图2、3和4中示出的反射光生物盘的半径获得的部分剖面图,其显示了在其中形成的流动沟道;
图12是垂直于图5、6和7中示出的透射光生物盘的半径获得的部分剖面图,其显示了在其中形成的流动沟道和单一的顶部探测器;
图13是图2、3和4中显示的反射光生物盘的部分纵向剖面图,其示出了在其中形成的摆动凹槽(wobble groove);
图14是图5、6和7中显示的透射光生物盘的部分纵向剖面图,其示出了在其中形成的摆动凹槽(wobble groove)和顶部探测器;
图15是类似于图11的简图,显示了反射光盘的全部厚度及其初始折射性能;
图16是类似于图12的简图,显示了透射光盘的全部厚度及其初始折射性能;
图17是显示使用本发明的方法处理从光生物盘上收集的信息的流程图;
图18是将取样模拟信号转变为存储成一维阵列的相应数字信号的图示;
图19是光盘的透视图,其具有象征性显示相对于光生物盘轨道定位的被捕获白细胞的放大细节图,该光生物盘在与入射波束相互作用之后产生含有信号的波束;
图20A是根据本发明相对于光生物盘轨道定位的白细胞的图示;
图20B是根据本发明从图20A的白细胞获得的一系列信号迹线(trace);
图21是示出图21A、21B、21C和21D之间关系的图示;
图21A、21B、21C和21D组合在一起是将图20B的信号迹线转变为数字信号的图示,其中数字信号存储成一维阵列并组合成用于数据输出的二维阵列;
图22是描述用于根据本发明的处理方法和计算法则进行数据估计的步骤的流程图;
图23是显示在根据本发明的一个实施例选择估计矩形时所含步骤的流程图;
图24是具有如根据本发明一个特殊实施例的软件所显示的窗口的光生物盘的图示;
图25是示出根据本发明一个实施例在从具有窗口的光生物盘上收集的调查数据中发现窗口时所含步骤的流程图;
图26显示了调查数据阵列的一个示例行,该调查数据阵列出于根据本发明的一个方案发现阈值的目的进行了扫描处理;
图27是显示在根据本发明的另一个方案从调查数据中提取关联成分时所含分步的流程图;
图28描述了根据本发明的一个实施例在软件显示上发现窗口之后修剪估计矩形的结果;
图29显示了无窗口光盘上的示例暗点和含有被捕获细胞的靶区;
图30是类似于图29的简图,其显示了如何根据本发明的一个实施例在无窗口的光盘中利用示例暗点发现细胞;
图31是显示根据本发明的某个方案对调查数据进行背景照明均匀化时所含步骤的流程图;
图32显示了在进行背景照明均匀化之前调查数据的实例,如本发明的软件所显示;
图33显示了在进行背景照明均匀化之后调查数据的实例,如本发明的软件所显示;
图34是示出根据本发明的特殊执行对调查数据进行标准化时所含步骤的流程图;
图35显示了逐步标准化期间如软件所显示的示例调查数据的图示;
图36显示了在标准化之后示例调查数据的图示,如本发明的软件所显示;
图37是显示根据本发明的优选实施例对示例调查数据进行过滤时所含步骤的流程图;
图38显示了在过滤步骤之后示例调查数据的图示,如本发明的软件所显示;
图39是图38所示图示的近视图(close-up view),其具有相应的点值图迹线(point value graph trace);
图40是显示在根据本发明某些方案的特殊实施例从调查数据中除去不需要的成分时所含步骤的流程图;
图41显示了在除去裂缝之前示例调查数据的图示,如本发明的软件所显示;
图42是在除去裂缝之后图41中示例调查数据的图示,如本发明的软件所显示;
图43是显示在根据本发明的明亮中心法标记和计数细胞时所含步骤的流程图;
图44显示了充满由明亮中心法计数的细胞的示例调查数据的图示;
图45显示了图44所示图示一部分的近视图(up-close view)和数值迹线图;
图46A是显示在根据本发明的黑暗边缘法标记和计数细胞时所含步骤的流程图;
图46B是使用移动环的卷积的图示;
图47是示例调查数据的图示,该调查数据中已计数细胞用十字标记,如本发明的软件所显示;
图48提供了显示在使用本发明不同方法所利用的算法提取红细胞时所含步骤的示例流程图;
图49显示了在执行图48所概述的算法之前含有红细胞的调查数据的图示;
图50示出了在执行了图48所概述算法的第一步之后图49的示例调查数据的图示;
图51描述了在应用了图48所概述算法的第二步之后图49的示例调查数据的图示;
图52可视地提供了在执行了图48所概述算法的第三步之后图49的示例调查数据;
图53显示了在执行了图48所概述算法的第四步之后图49的示例调查数据的图示;
图54示出了在应用执行了图48所概述算法的第五步之后图49的示例调查数据的图示;
图55是显示用图48所概述的算法计数的红细胞的近视图;
图56A是在使用本发明的绝对值计数法计数分散细胞(discretecells)之前它们的图示屏幕拍摄(pictorial screen shot);
图56B是显示在利用本发明绝对值计数法的一个实施例计数细胞时所含步骤的流程图;
图57是在执行了根据本发明绝对值计数法的该实施例的标准化和过滤步骤之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图58是在应用了根据本发明绝对值计数法的背景去除和二值化步骤之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图59是在执行了根据本发明绝对值计数法的规则化步骤之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图60是在应用了根据本发明绝对值计数法的所示出实施例的一个象素宽度边界提取(one-pixel wide boundary extraction)步骤之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图61是在执行了根据本发明绝对值计数法的填充成分步骤之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图62是在应用了根据本发明绝对值计数法的该特殊实施例的填充调查数据步骤之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图63是在根据本发明的方法对分散细胞进行计数并用十字标记之后最初在图56A中显示的分散细胞的图示屏幕拍摄;
图64显示了用于计数凝集和分散红细胞的绝对值计数法的结果;
图65是显示在根据本发明的可选择实施例对图象进行快速傅立叶变换时所含步骤的流程图;
图66是在根据本发明进行快速傅立叶变换之前示例调查数据的图示;
图67显示了在执行了快速傅立叶变换之后图66的示例调查数据的图示;
图68示出了在重新对准(realignment)之前调查数据歪斜图象表示(skewed graphical representation)的实例;
图69描述了图68所示图象表示的歪斜方向;
图70显示了在重新对准之后图68的图象表示;
图71A是描述在为气泡痕迹情况修正细胞计数时所含步骤的流程图;
图71B是根据本发明的另一个方案为气泡痕迹情况修正细胞计数的图象表示;
图71C是通过被捕获红细胞靶区的气泡痕迹的实例,如在5X显微镜下所见到的;
图71D是图71中的气泡痕迹和被捕获红细胞周围的放大图,如在40X显微镜下所见到的;和
图72是显示使用本发明的方法分析血液样品的图示流程图。
具体实施方式
本发明致力于成像和计数实验室样品中的细胞物质的方法和装置。这些方法和装置可以用于成像和计数光盘上或光盘内的目标(interest)的任何类型调查特征。本发明的实施例产生样品中调查特征或细胞的调查数据,并对调查数据进行计算机分析。在其它化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。
在随后的描述中,提出了大量的特殊细节以便为本发明的实施例提供更完全的说明。然而本领域的技术人员应当意识到,本发明的实现可以没有这些特殊细节。在其它的例子中,没有详细描述众所周知的特征以便不掩蔽本发明。
这里更详细地讨论使用光盘数据进行白细胞计数的大量实施例。这些实施例并不仅限于成像和计数白细胞,而是可以容易地用于进行任何类型细胞物质的计数。这可以包括,但不仅限于,红细胞、白细胞、珠子和任何其它产生类似的能够被光学阅读器探测的光学信号的物质,包括生物的和非生物的。在其它化验中,目标的调查特征,除了细胞以外,可以替换为珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。下面在讨论细胞计数时进一步详细说明了将本发明用于除白细胞之外的物质所需要的一些修改。
在下面的讨论中,提供了两个主要部分,用于示出关于本发明的数据收集和数据分析方案。第一部分提供了用于从实验室样品中收集调查数据并将该调查数据转换成基于阵列的存储装置、方法和算法的详细说明。第二部分提供了用于分析调查数据的方法和算法的详细说明。这两个部分之后,提供了关于执行白细胞计数化验的方法的部分。
1.数据收集
本发明的实施例包括检索(retrieval)调查数据,该调查数据来自实验室样品中的细胞物质。图1是根据本发明的光生物盘110的透视图。本光生物盘110结合光盘驱动112和显示器114一起显示。测试样品沉积在光生物盘110的指定区域上。一旦光生物盘插入到光盘驱动112中,光盘驱动便使用电磁辐射波束从样品收集信息,该波束通过与测试样品相互作用被调整或调制。在对信息进行分析和处理之后,计算机监视器114显示结果。
能够在本发明中使用的光生物盘110有两个主要的实施例。图2、3和4示出了光生物盘110的反射实施例,而图5、6和7示出了光生物盘110的透射实施例。
A.反射实施例
图2是光生物盘110一个实施例的结构元件的分解透视图。图2是可以在本发明中使用的反射区光生物盘110(下文称为“反射光盘”)的实例。结构元件包括帽部分(cap portion)116、粘合剂或沟道部件118和基片120。帽部分116包括一个或多个入口122和一个或多个出口(vent port)124。帽部分116可以由聚碳酸酯形成,并优选地在其底部上涂布反射表面146(如图4中更详细示出的),如透视图2所示。在优选实施例中,在反射层142的表面上含有触发标记(triggermarking)126,见图4。触发标记126可以包括位于光生物盘所有三个层内的干净窗口(clear window)、不透光区域、或者用信息编码的反射或半反射区域。编码信息用于向处理器166(如图10A所示)发送数据,处理器166顺次与图8和10A所示的询问或入射波束152的操作功能相互作用。图2所示的第二元件是粘合剂或沟道部件118,在其中形成了流体线路128或U形沟道。流体线路128优选地通过冲压或切除膜以除去塑料膜并形成所示的形状而形成。每个流体线路128都包括流动沟道130和返回沟道(return channel)132。图2中示出的一些流体线路128包括混合室134。示出了两种不同类型的混合室134。第一种是对称混合室136,其相对于流动沟道130对称地形成。第二种是偏置混合室138。偏置混合室138形成于图示流动沟道130的一侧。图2中示出的第三元件是基片120,其包括靶或捕获区140。基片120优选地用聚碳酸酯制成,并在其顶部沉积有反射层142,见图4。靶区140通过除去呈图示形状或者选择地呈任何期望形状的反射层142而形成。选择地,靶区140可以通过掩模技术形成,该技术包括在施加反射层142之前掩蔽靶区140。反射层142可以由金属例如铝或金形成。
图3是图2中示出的光生物盘110的顶视平面图,其在帽部分116上具有反射层142,帽部分116被透明地显示以显现位于光盘内的流体线路128、靶或捕获区140和触发标记126。因为每个捕获区都具有一种或多种特异抗原以锁定(lock down)样品中的不同组分(或不同细胞),所以在化验处理之后,室内的每个捕获区都含有一类细胞或细胞组分。锁定或捕获通过使一种或多种抗原具有能够“锁定”在血细胞的特定组分上从而捕获特定细胞的化学结构而实现。分离细胞组分对于对血细胞,例如白细胞,进行差分计数是至关重要的。在其它化验中,除了细胞以外,调查特征能够是珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。靶或捕获区140限定了电磁询问波束与测试样品相互作用的位置。
图4是根据本发明一个实施例的反射区类型光生物盘110的放大透视图。该图包括其各种层的部分,并经过切除从而示出每个层、基片、涂层或膜的部分剖面图。图4显示了涂布有反射层142的基片120。活性层(active layer)144施加在反射层142上面。在优选实施例中,活性层144可以由聚苯乙烯制成。选择地,可以使用聚碳酸酯、金、活化玻璃、经过改性的玻璃、或者经过改性的聚苯乙烯例如聚苯乙烯交联顺丁烯二酸酐(polystyrene-co-maleic anhydride)。此外,能够使用水凝胶。如本特殊实施例所例证的,塑料粘合剂部件118施加在活性层144上。塑料粘合剂部件118的暴露部分示出了被切除或冲压的U形形状,其产生了流体线路128。本光生物盘该反射区实施例中的最终结构层是帽部分116。帽部分116包括在其顶部上的反射表面146。反射表面146可以由金属如铝或金制成。
B.透射实施例
图5是根据本发明的透射类型光生物盘110结构元件的分解透视图。透射型光生物盘110的结构元件类似地包括帽部分116、粘合剂或沟道部件118和基片120层。帽部分116包括一个或多个入口122和一个或多个出口124。帽部分116可以由聚碳酸酯层形成。可选择的触发标记126可以包含在薄半反射层143的表面上,如图7和8所详细示出的。触发标记126可以包括位于光生物盘所有三个层内的干净窗口、不透明区、或者用信息编码的反射或半反射区。编码信息用于向处理器166(图10A所示)发送数据,处理器166然后与图8和10A所示的询问或入射波束152的操作功能相互作用。
图5所示的第二元件是粘合剂或沟道部件118,在其中形成了流体线路128或U形沟道。流体线路128通过冲压或切除膜以除去塑料膜并形成所示的形状而形成。每个流体线路128都包括流动沟道130和返回沟道(return channel)132。图5中示出的一些流体线路128包括混合室134。示出了两种不同类型的混合室134。第一种是对称混合室136,其相对于流动沟道130对称地形成。第二种是偏置混合室138。偏置混合室138形成于图示流动沟道130的一侧。
图5中示出的第三元件是基片120,其包括靶或捕获区140。基片120优选地用聚碳酸酯制成,并在其顶部沉积有反射层143,见图8。与图5中示出的光盘110的基片120连接的半反射层143显著薄于图2、3和4中示出的反射光盘110基片120上的反射层142。更薄的反射层143允许一些询问波束152透射通过透射光盘的结构层,如图11所示。薄半反射层143可以由金属如铝或金形成。
图6是图4中示出的透射型光生物盘的顶视平面图,其具有显示位于光盘内的流体沟道、触发标记126和靶或捕获区140的透明帽部分116。靶或捕获区140是电磁波束与测试样品相互作用的位置。在光盘旋转之后,样品内细胞的特定组分被预先装载在室内的各种捕获剂或抗原捕获在不同捕获区内。
图7是根据本发明透射光盘的光生物盘110的放大透视图。光盘110示出了具有各种层的部分,它们被切除掉以显示每个层、基片、涂层或膜的部分剖面图。图7示出了具有透明帽部分116、位于基片120上的薄半反射层143和触发标记126的透射光盘格式。触发标记126包括布置在帽顶部上的不透明材料。选择地,触发标记126可以由干净的、在光盘的薄反射层143上被腐蚀的非反射窗口,或者任何吸收或不反射来自触发探测器160的信息的标记形成,见图10A。图7还显示,靶区140通过标记呈图示形状,或者选择地呈任何期望形状的指定区域形成。指示靶区140的标记可以在位于基片120上或者位于基片120底部上(在光盘下)的薄半反射层143上制成。选择地,靶区140可以通过掩模技术形成,该技术包括掩蔽除了靶区140之外的整个薄半反射层143。在这一实施例中,靶区140可以通过丝网油墨(silk screening ink)建立到薄半反射层143上产生。活性层144施加在薄半反射层143上。在优选实施例中,活性层144是40-200μm的2%聚苯乙烯厚层。选择地,可以使用聚碳酸酯、金、活化玻璃、经过改性的玻璃、或者经过改性的聚苯乙烯例如聚苯乙烯交联顺丁烯二酸酐(polystyrene-co-maleic anhydride)。此外,能够使用水凝胶。如本实施例所示出的,塑料粘合剂部件118施加在活性层144上。塑料粘合剂部件118的暴露部分示出了产生流体线路128的被切除或冲压掉的U形形状。本光生物盘110的这一透射实施例的最终结构层是干净、非反射帽部分116,其包括入口122和出口124。
C.光盘实施例的光学性能
两个光盘实施例之间的主要差异之一是光盘顶层涂层的厚度。在透射光盘的实例中,在基片层120的顶部上沉积薄的半反射层143。在反射光盘实例中,在其基片层120的顶部上沉积基本上更厚的半反射层。在图8中示出的优选实施例中,透射光盘的薄半反射层143的厚度大约为100-300,不超过400。这是因为金膜层在厚度大于800时完全反射,而在厚度低于大约400时允许光线透过金膜。如下所示,表2提供了金膜相对于膜厚度的反射和透射性能。
表2 Au膜反射率和透射率(绝对值)
厚度(埃) 厚度(nm)   反射率   透射率
    0     0     0.0505     0.9495
    50     5     0.1683     0.7709
    100     10     0.3981     0.5169
    150     15     0.5873     0.3264
    200     20     0.7142     0.2057
    250     25     0.7959     0.1314
    300     30     0.8488     0.0851
    350     35     0.8836     0.0557
    400     40     0.9067     0.0368
    450     45     0.9222     0.0244
    500     50     0.9328     0.0163
    550     55     0.9399     0.0109
    600     60     0.9488     0.0073
    650     65     0.9482     0.0049
    700     70     0.9505     0.0033
    750     75     0.9520     0.0022
    800     80     0.9531     0.0015
光线透过金膜的阈值密度大约是400。除表2之外,图9提供了薄半反射层143根据金的厚度其反射和透射性质成反比例图示。图9中示出的曲线图所使用的反射和透射值是绝对值。如图8所示,更薄的半反射层143允许部分的入射或询问波束152穿透和通过。从而入射或询问波束152能够由如图10A所示的顶部探测器158加以探测。
在反射光生物盘的实例中,返回波束154携带关于生物样品的信息。如上面所讨论的,基本上只有在入射波束处于流动沟道130或靶(或捕获)区140内从而与样品接触时,该关于生物样品的信息才包含在返回波束中。返回波束154也可以携带编码在反射层142内或上或者编码在摆动凹槽170内的信息,如图13和14所示出的。正如对于本领域技术人员所显而易见的,只有在相应的入射波束与反射层142接触时,预先记录的信息才包含在具有靶或捕获区的反射光盘的返回波束154中。当入射波束152位于除去了或者不存在信息携带反射层(information bearing reflective layer)142的区域内时,该信息不包含在返回波束154中。
本发明的方法也可以容易地用于包括等半径沟道(equi-radiochannels)的光生物盘,例如在共同转让的美国临时专利申请序列No.60/353,014中所公开的,其题目为“包括等半径和/或螺旋分析区的光盘及相关的光盘驱动系统和方法”(Optical Discs IncludingEqui-Radial and/or Spiral Analysis Zones and Related Disc DriveSystems and Methods),其于2002年1月29日提出申请,在这里引用作为参考。
D.系统装置
图10A是示出系统装置操作的透视图和框图,该系统装置包括光学组件148,和产生入射或询问波束152、返回波束154和透射波束156的光源150。在反射光生物盘实例中,返回波束154从光生物盘110帽部分116的反射表面146反射。在本光生物盘110的该反射实施例中,返回波束154由底部探测器157加以探测和用于分析信号媒介的存在。在透射光生物盘实施例中,透射波束通过顶部探测器158加以探测,也用于分析信号媒介的存在。在透射实施例中,光电探测器可以用作顶部探测器158。
图10A还显示了硬件触发机构,其包括光盘上的触发标记126和触发探测器160。硬件触发机构在反射光生物盘和透射光生物盘中都使用。触发机构允许处理器166只有在询问波束152处于各自的靶或捕获区140上时收集数据。而且,在透射光生物盘系统中,还使用软件触发器。软件触发器使用底部探测器向处理器166发送信号以便在询问波束152到达各自靶或捕获区140的边缘时收集数据。图10A还示出了驱动马达162和用于控制光生物盘110旋转的控制器164。图10A进一步显示了在可选择实施例中实现的处理器166和分析器168,它们用于处理返回波束154和与透射光生物盘相伴的透射波束156。在透射光生物盘的例子中,透射波束156携带关于生物样品的信息。在该实施例中,光盘上有预先记录的信息。探测器158收集该波束。
在本发明的另一个实施例中,分波顶部探测器用于收集透射波束156。图10B显示了根据本发明一个实施例的分波探测器。探测器170具有两个探测器组件172和174。这两个探测器组件收集被物体186(例如细胞)折射的透射波束156并产生两个信号A和B。物体186能够是调查特征,例如生物细胞如红或白细胞。通过从一个信号中减去另一个信号(也就是A-B或B-A)能够获得差分信号。当探测器组件位于含有散射入射波束152的物体的区域上方时,它们探测信号中的变化。每个探测器察看到与另一个探测器的所见相反的变化。换言之,当光线向一个探测器弯曲时,它看到信号增加,而另一个探测器看到信号减小。由于这个性质,能够通过产生如下信号显著地增加信噪比,该信号是两个探测器的每一个所产生信号的差分。该差分信号具有两个优点。第一,差分信号中消除了系统中的任何等量影响两个探测器的噪音。第二,光盘上能够折射光线而不只是吸收光线的目标物体会在差分信号中引起较大而且容易探测的变化。这有助于分析任务,其需要分离从背景噪音中由目标物体产生的信号。
分波探测器更全面的讨论存在于共同拥有的美国临时专利申请No.60/355,090中,其于2002年2月14日提出申请,题目为“用于生物驱动的分割区域探测器及其相关方法”(Segmented Area Detectorfor BioDrive and Methods Relating Thereto),和相同题目的相关临时专利申请中,其序列编号为No.60/335,123;60/352,649;60/353,739和60/355,090,它们分别于2001年10月10日;2002年1月28日;2002年1月30日和2002年4月7日提出申请,它们全部引用在这里作为参考。能够与本发明联合使用的不同类型探测器更全面的讨论存在于共同拥有的美国专利申请No.10/043,688,其与2002年1月10日提出申请,题目为“光盘分析系统及用于生物和医学成像的相关方法”(optical Disc Analysis System Including Related Methods ForBiological and Medical Imaging),其也在这里引用作为参考。
图11-16显示了反射和透射实施例的剖面图,用于示出光盘的光学性质和探测器如何用于从光盘收集携带信息的波束。
现在更详细地参考图11,其显示了根据本发明的光生物盘110反射光盘实施例的部分剖面图。图11示出了基片120和反射层142。如上所示,反射层142可以由材料例如铝、金或其它合适的反射材料制成。在该实施例中,基片120的顶表面光滑。图11还显示了施加在反射层142上的活性层144。如图11所示,靶区140通过除去反射层142期望位置处的区域或部分而形成,或者选择地,通过在施加反射层142之前掩蔽期望区域而形成。如图11进一步示出的,塑料沟道部件118施加在活性层144上。图11还显示了帽部分116和与之连接的反射表面146。这样,当帽部分116施加在含有期望切除形状的塑料沟道部件118上时,借此形成流动沟道130。如图11显示的箭头所示,入射波束152的路径从光盘110的下方朝基片120初始导向。然后入射波束聚焦在紧邻反射层142的点处。因为这一聚焦发生在靶区140内,在该处除去了反射层142的一部分,所以入射波束继续沿着路径通过活性层144并且进入到流动沟道130内。然后入射波束152继续向上传播通过流动沟道最终入射(fall incident)到反射表面146上。在该点处,入射波束152沿着入射路径返回或反射从而形成返回波束154。
图12是根据本发明的光生物盘110透射实施例的部分剖面图。图12示出了具有干净帽部分116和位于基片120上的半反射层143的透射光盘格式。图12还显示了施加在薄半反射层143上的活性层144。在优选实施例中,透射光盘具有用金属例如铝或金制成的薄半反射层143,其厚度为大约100-300埃,优选地不超过400埃。该薄半反射层143允许部分来自光盘150(见图10A)的入射或询问波束152穿过并向上通过待被顶部探测器158探测的光盘,同时一些光线被沿着与入射波束相同的路径但沿着相反的方向反射回来。在这一布置中,返回或反射波束154从半反射层143反射。这样在该方式中,返回波束154不进入流动沟道130。反射光线或返回波束154可以用于跟踪预先记录信息轨道上的入射波束152,预先记录信息轨道形成于半反射层143内或上,如图13和14更详细描述的。在图12示出的光盘实施例中,可以存在或者可以不存在所限定的靶区140。靶区140可以由在基片120的薄半反射层143上制成的直接标记产生。这些标记可以使用丝网印刷或者任何等效的方法实现。在可选择实施例中,其不使用物理标记限定靶区,流动沟道130有效地用作受限制的(confined)靶区域,在其中对调查特征进行检查。
图13是根据本发明的光生物盘110反射光盘实施例轨道的纵向(taken across)剖面图。该图是纵向沿着光盘的半径和流动沟道获取的。图13包括基片120和反射层142。在该实施例中,基片120包括凹槽170系列。凹槽170呈从光盘的中心附近向外周延伸的螺旋形状。实现凹槽170使得询问波束152可以沿着光盘上的螺旋凹槽寻轨(track)。这种类型的凹槽170称作“摆动凹槽”。具有波浪或波形侧壁的底部分形成凹槽170,同时被提高或提升的部分隔离螺旋中的相邻凹槽170。施加在本实施例中的凹槽170上的反射层142,如图所示,实际上(in nature)是共形的(conformal)。图13还显示了施加在反射层142上的活性层144。如图13所示,靶区140通过除去反射层142期望位置处的区域或部分形成,或者选择地,通过在施加反射层142之前掩蔽期望区域形成。如图13中进一步示出的,塑料粘合剂或沟道部件118施加在活性层144上。图13还显示了帽部分116和与之连接的反射表面146。这样,当帽部分116施加在含有期望切除形状的塑料粘合剂部件118上时,借此形成流动沟道130。
图14是根据本发明的光生物盘110透射光盘实施例,如图12所描述的轨道的纵向(taken across)剖面图。该图是纵向沿着光盘的半径和流动沟道获取的。图14示出基片120和薄半反射层143。该薄半反射层143允许来自光源150的入射或询问波束152穿过和通过待被顶部探测器158探测的光盘,同时一些光线以返回波束154的形式反射回来。薄半反射层143的厚度由光盘阅读器维持其寻轨能力(trackability)所需的最小量反射光线决定。本实施例中的基片120与图13中所讨论的相似,包括凹槽170系列。本实施例中的凹槽170也优选地呈从光盘的中心附近向外周延伸的螺旋形状。实现凹槽170使得询问波束152可以沿着凹槽寻轨。图14还显示了施加在薄半反射层143上的活性层144。如图14中进一步示出的,塑料粘合剂或沟道部件118施加在活性层144上。图14还显示了无反射表面146的帽部分116。这样,当帽部分116施加在含有期望切除形状的塑料粘合剂部件118上时,借此形成流动沟道130并允许部分入射波束152基本上无反射地通过其间。
图15是与图11类似的简图,显示了反射光盘的全部厚度及其初始折射性能。图16是与图12类似的简图,显示了透射光盘的全部厚度及其初始折射性能。凹槽170在图15和16中没有显示,因为这些部分是沿着凹槽170切割的。图15和16显示了狭窄流动沟道130的存在,其在这些实施例中垂直于凹槽170定位。图13、14、15和16显示了各自反射和透射光盘的全部厚度。在这些图中,入射波束152被示出与基片120初始相互作用,基片120具有改变入射波束路径的折射性能,如图所示,从而将波束152聚焦在反射层142上或者薄半反射层143上。
E.模拟-数字处理
无论是从反射光盘的返回波束154获得的还是从透射光盘的透射波束获得的,关于生物测试样品的信息都被导向到处理器166(见图10A)用于信号处理。该处理包括将由底部探测器157(反射光盘)或顶部探测器158(透射光盘)探测的模拟信号转换为离散的数字形式。
图17是与图10B所示装置相关的信息检索处理的概要流程图。在步骤270中,如果实施例是透射光生物盘,那么通过探测器158探测携带关于生物样品信息的透射波束。如果实施例是反射光生物盘,那么在步骤272中通过探测器157探测反射波束154。在任何一个实例中,在步骤274将信息发送到模拟-数字转换。在步骤276中,结果数字数据是阵列。
图18显示了通过处理器166执行的模拟-数字转换。该转换包括以固定的时间间隔212取样模拟信号210,和将信号的相应瞬时模拟振幅214编码成为离散的二进制整数216。取样在某一开始时间218开始,并在某一结束时间220停止。与任何模拟-数字转换处理相关的两个通用值是取样频率和比特深度(bit depth)。取样频率,也称作取样速度,是每单位时间所取的样品数目。较高的取样频率在连续样品之间产生较小的时间间隔212,这导致数字信号222与原始模拟信号210相比具有更高的保真度。比特深度是每个指向(point to)编码模拟信号210的采样振幅214的样品中使用的比特数目。比特深度越大,二进制整数与原始模拟振幅214的近似度越好。在本发明的一个实施例中,取样速度是8MHz,每个样品的比特深度是12比特,其允许的整数样品范围是0-4,095(0至2n-1,其中n是比特深度)。
在其它实施例中,比特深度和取样频率的组合能够被自定义(customized)以适应特殊的精度需要。通过非限制性的实施例,可以期望提高实施例中的取样频率,该实施例包括用于计数通常比细胞小的珠子的方法。在模拟-数字转换期间,沿激光路径的每个连续样品点224都连续存储在光盘上或存储器内成为一维阵列226。每个连续轨道提供一个独立的一维阵列。所有的一维阵列都组合在一起形成二维阵列228(如图21B所示),其类似于普通图象表示。
这里提供数据收集实例以进一步示出在从光生物盘收集数据时涉及的细节。图19显示了本发明的光生物盘110的透视图。图19包括被指定显示相对于光生物盘轨道232定位的被捕获白细胞230的部分的放大详细透视图。如图所示,入射波束152与白细胞230的相互作用产生了含有信号的波束,其或者处于反射光盘返回波束154的形式,或者处于透射光盘透射波束156的形式,这些光盘由底部探测器157或顶部探测器158加以探测。
图20A、20B和图21A-21D示出了细胞如何被捕获进入数字数据。在其它化验中,目标的调查特征不是细胞,而是珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。图20A是相对于光生物盘110的轨道232定位的白细胞230的图示。细胞230定位在类似于图19所示光盘的光盘上。图20B是根据本发明源自于图20A中的白细胞230的信号轨迹系列。图20B描绘了标志有A、B、C和D的相应轨迹。然后模拟信号轨迹(信号)210被导向到处理器166用于转换成相应的数字信号222(如图21A-21D所示)。图20B进一步显示出,对白细胞230的扫描产生了能够被探测和处理的入射波束的扰动(perturbation)231。
图21是示出图21A、21B、21C和21D之间布置关系的图示,这些图组合在一起示出了来自图20B的四个轨道A、B、C和D如何转换为单一的二维数字数据阵列228。
现在特别参考图21A,其显示了来自图20A所示光生物盘轨道A和B的取样模拟信号210。然后处理器166将模拟信号210的相应瞬时模拟振幅214编码成离散的二进制整数216(见图12)。数据点的结果系列是类似于取样模拟信号210的数字信号222。
现在转移到图21B,来自轨道A和B(图21A)的数字信号存储为独立的一维存储阵列226。每个连续轨道都提供相应的一维存储阵列,其在与先前的一维阵列组合时产生数字数据的二维阵列228。然后数字数据存储在存储器中或光盘上作为样品点224(图18)的二维阵列228,该样品点表示样品区域中特殊点处的返回波束154或透射波束156的相对密度。然后二维阵列以原文件(raw file)或图象文件240的形式存储在存储器中或光盘上。然后从存储器242检索存储在文件240中的数据,并用作输入到分析器168(图10A)的数据输出244。
图21C显示了来自图20A所示光生物盘轨道C和D的取样模拟信号210。然后处理器166将模拟信号210的相应瞬时模拟振幅214编码成离散的二进制整数216(图18)。数据点的结果系列是与取样模拟信号210类似的数字信号222。
现在参考图21D,来自轨道C和D(图21C)的数字信号222存储为独立的一维存储阵列226。每个连续轨道提供一个相应的一维阵列,其在与先前的一维阵列组合时产生与图象类似的二维阵列228(图21B)。如上,数字数据而后存储在存储器内或光盘上作为样品点224(图18)的二维阵列228,样片点表示样品区域特殊点处返回波束154或透射波束156(图19)的相对强度。而后二维阵列以原文件、数据文件或图象文件240的形式存储在存储器内或光盘上。然后从存储器240检索存储在文件240内的数据,并用作向分析器168(图10A)的数据输出。
从光生物盘捕获数据并将该数据转换成二维整数阵列的附加方法和算法具有广泛的应用性,并且在如下专利申请中公开,即美国临时专利申请序列No.60/291,233,题目为“用于在光生物盘组件中提供分析结果像素的可变取样控制及相关装置”(Variable SamplingControl for Rendering Pixelation of Analysis Results in OpticalBio-Disc Assembly and Apparatus Relating Thereto),其于2001年5月16日提出申请,并在这里引用作为参考。
本发明的另一个实施例将调查数据存储在存档内。该存档提供能够对调查数据加以分类的空间。随后,能够分析各组数据指导不同人类群体的健康趋势研究。例如,调查数据能够联系患者的信息产生能够通过患者的属性(attributes)加以分类的调查数据目录。诸如年龄、性别、种族和血型的信息例如能够用于分类调查数据。该存档能够利用可搜索的相关数据库。一旦建立这样的存档,就能够对特定分类的调查数据进行分析。例如,人口健康趋势研究可以通过检索调查数据并分析这些调查数据加以执行,调查数据由特定城市的患者捐赠。该优势在于研究的执行不需要人群本身到场。一段时间之后,能够建立历史存档,并能够对一段时期的特殊人口进行研究以及能够分析一段时间的趋势。
II.数据分析
下述部分涉及本发明的数据分析方案。这里特别地结合图22-71D并一般地参考先前的图1-21进行讨论。
A.数据收集和处理
在一个实施例中,对来自光生物盘的调查数据加以分析用于细胞计数,其中调查数据以数字数据的形式存储。在另一个实施例中,使用其它形式的调查数据进行分析。在其它化验中,调查数据能够含有用于计数珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记的信息。本发明的一个实施例将该调查数据实时发送到数据分析器。在另一个实施例中,存储调查数据并在后来加以检索进行分析。在两个实施例中,本发明的计算和处理算法存储在分析器168(图10A)内,并应用于输入调查数据244以产生有用的输出结果262(图22),该结果可以显示在显示器114(图10A)上。作为非限制性实例,下文说明了用于处于数字数据阵列形式的调查数据的分析方法。考虑到本发明,本领域技术人员可以意识到,除了数字数据阵列格式以外,本发明的方法能够用于各种形式的调查数据。
现在继续到图22,其显示的流程图提供了根据本发明的处理方法和计算算法进行数据分析的步骤的一般概览。本处理方法的第一步包括收到输入调查数据244。如上所述,数据分析以二维整数阵列开始,其范围为0-4,095。下一步246是选择用于计数的光盘估计矩形。一旦限定该矩形,目标便成为得出矩形内所含全部白细胞的实际数量。这一区域称作“调查数据区域”。步骤246的实现取决于光盘的构型。两种可能的光盘构型包括具有窗口的光盘和没有窗口的光盘。
作为非限制性的实例,在本发明使用的具有窗口,例如图2和4所示靶或捕获区140,的光盘实施例中,软件识别窗口并切割其一部分用于分析和计数。在一个优选实施例中,如图2所示出的,每个窗口(或捕获区)都具有矩形的形状,该矩形为1×2mm并在其每端(end)上具有半圆形部分。在该实施例中,软件在各自窗口内切割1×2mm的标准尺寸估计矩形。在该实施例的方案中,阅读器可以获取数个连续的样品值,以比较数个不同窗口内的细胞数目。
在本发明使用的无窗口透射光盘的实施例中,如图4所示,步骤246以一种或两种方式执行。标准矩形位置的选择或者通过相对于具有固定坐标的点定位其中心或者通过发现校准点执行,校准点优选地是具有特殊特性的暗染色点。在采用校准点的实例中,具有期望对比度的染料相对于两个细胞群沉积在光盘的特殊位置处。然后光盘阅读器被引导跳跃到其中一个细胞群的中心,而后标准尺寸(1×2mm)估计矩形居中在所选择群上。在其它化验中,细胞可以替换为珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。
除了提供两种类型的光盘之外,步骤246还允许用户进行选择。用户可以通过直接用鼠标选择或者其它方法直接交互地指定期望的用于细胞计数的样品区域形状,例如矩形区域。在软件的这一实施例中,涉及使用鼠标在显示于监视器114(图1)上的调查数据图示的期望部分上点击并拖动一个形状。不管估计区域的选择方法,在下一步骤248对各个矩形区域进行估计用以计数。
图22的第三步是步骤248,其致力于背景照明均匀化。该处理通过一硬件构型修正可能的背景均匀性波动。背景照明均匀化补偿每个样品点的强度水平,从而使整个背景或者非细胞的调查数据的一部分接近一个具有任意背景值V背景(Vbackground)的平面。尽管V背景可以通过许多方式决定,例如采用整个标准矩形样品区的平均值,但在本实施例中V背景被设定为2,000。所选矩形样品区每个点P处的值V用数值(V背景+(P附近区域的V平均值))代替。如果需要,结果V可以缩小以适合数值实际可能的范围,其在本发明优选实施例中为0-4,095。相邻矩形的尺寸选择得比细胞的尺寸充分地大,比标准矩形的尺寸充分地小。在其它化验中,为珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记选择相邻矩形。
图22中流程图的下一步是标准化步骤250。在执行标准化步骤250中,对标准矩形样品区域内的数据进行线性转换,从而使平均值为2,000,标准偏差为600。如果需要,缩小该数值以适应0-4,095的范围。这一步骤250以及背景照明均匀化步骤248,使软件对硬件修改和调谐较不敏感。作为非限制性实例,探测电路例如顶部探测器158(图13)内的信号增益可以改变而不会显著影响结果细胞计数。
如图22所示,接下去执行过滤步骤252。对于标准矩形内的每个点P,计算P附近区域内点的数目,该区域的尺寸小于步骤248所示的尺寸,数值与V背景显著不同。被计算的点的数目应当接近调查数据中细胞的尺寸。在其它化验中,被计算点的数目应当接近于目标物,目标物包括珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。如果所发现区别点的数目足够大,P处的数值保持不变;否则将其指定为V背景。执行该过滤操作是为了除去噪音,并且在最佳情况下,只有细胞保留在调查数据中,而背景均一地等于V背景
可以执行致力于除去坏成分的可选择步骤254,如图22所示。缺陷例如划痕、气泡、灰尘以及其它类似的不规则,可以经过过滤步骤252。这些缺陷可以通过直接或间接地影响调查数据直方图中的整体分布导致细胞计数误差。该步骤利用如下事实,即这些缺陷的尺寸显著大于细胞的尺寸,并通过适当的算法除去它们。在其它化验中,其可以实际地应用到其它调查特征中,例如珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。在确定除去缺陷有多好时要考虑它们的尺寸。在可选择步骤254之后,优选地重复步骤248、250和252。
如图22所示,下一处理步骤是步骤256,其致力于通过明亮中心计数细胞。计数步骤256包括数个分步骤。它们是(1)利用卷积使细胞中心更加可见,(2)标记这些中心,和(3)对细胞进行真实计数。在一些硬件构型中,一些细胞可以显示出没有明亮中心。在这些例子中,只有黑暗边缘可见,从而下述的两个可选择步骤258和260是有用的。
步骤258致力于从图象中除去所发现的细胞。在步骤258中,围绕每个已发现细胞的圆形区域填充数值2,000(背景缺省值)从而使同时具有明亮中心和黑暗边缘的细胞不会被发现两次。步骤260致力于通过识别黑暗边缘计数额外的细胞。步骤258之后对调查数据进行两种转换。目的在于使黑暗边缘更加明显从而能够计数在步骤256中未通过明亮中心方法计数的细胞。在另一个实施例中,通过黑暗边缘计数细胞的方法能够代替通过明亮中心计数细胞的方法。
在计数步骤256之后,或者在选择地采用计数步骤260之后,图2所示的最后一个步骤是结果输出步骤262。在标准矩形中发现的细胞数目显示在图1所示的监视器114上。所识别的每个细胞在所显示的得自光生物盘的调查数据上用十字标记。
下面联系详细概述分步骤的框图给出图22所示每一步骤的详细说明。
步骤1:数据输入
步骤244检索存储在二维阵列内的范围为0-4,095的数据。黑色部分填充数值常数0,而光线探测的面积范围为1-4,095。在结果数据估计中,忽略0。
步骤2a:在具有窗口的光盘中选择估计矩形
图23显示了选择估计矩形的详细处理过程(图22的步骤246)。在步骤300,确定选择的类型。第一选择302包括从光生物盘选择估计矩形,该光生物盘具有物理嵌入的窗口,例如图2和4所示的靶或捕获区140。图24显示了具有窗口326和328的光盘的调查数据实例的图示,其通过根据本发明实施例的软件加以显示。
根据一个实施例的步骤304的详细处理如图25所示。在步骤330,图25,调查数据阵列以如下方式压缩,即只考虑每个第n行和每个第n列。然后在步骤332,以逐行的方式扫描经过压缩的调查数据以确定二值化步骤中使用的阈值。图26提供了用于示出的示例。在示例行342中,每个细胞代表调查数据上的一个点,且每个细胞内的数值表示该点处探测到的光线强度。扫描通过为调查数据阵列的每一行选择所有可能的段(segment)长度L开始。从而,图26为来自调查数据阵列的该示例行显示了所有可能的段长度L。实际阵列将具有许多这样的示例行。L被选择为略小于光盘上窗口的宽度。然后,使用所考虑段内的所有整数值为每个段计算平均值。因为段沿着行“滑动”,所以该处理称作发现“滑动平均值”。一旦发现了所有行段长度L的平均值,便确定了平均值的最小值和最大值。阈值T计算为(最小(平均值)+最大(平均值))/2。该处理减少了搜索窗口以发现覆盖窗口区域的段的任务量。因为窗口区域比非窗口区域亮,所以这些段的平均值高于阈值T,并且在后面的步骤能够更容易地识别。
一个实施例加速了段L平均值的计算,如下。在计算段中所有数值的总和a(n)时,其中n从K+1到K+L(K+1是段的开始点),a(K+L)加到总和a(n)上并从总和a(n)减去a(K)。对K的所有数值重复该步。这使得算法不必在从行中移下一个单元时将段所有数值的总和再重新全加一遍。对每一行重复该发现阈值的整个处理直到以这种方法扫描完所有行为止。
返回图25,在步骤334中执行二值化。在该步中,数值超过阈值T的点显示为黑色,而其余的则显示为白色。所以现在调查数据能够象是黑色和白色的图象一样加以处理。随后的二值化,在步骤336对调查数据进行规则化。规则化包括两个部分:侵蚀(erosion)和扩展(expansion)。侵蚀的执行如下。对于图象P,构建相应的图象P’。如果(1)P中相应点X为白色,或者(2)X的任何相邻点为白色,则P’中的点X’显示为白色。如果每个条件都未满足,则X’显示为黑色。P’是侵蚀的结果图象。以相反的方式执行扩展。对于图象R,构建相应的图象R’。如果(1)R中相应点X为黑色,或者(2)Y的任何相邻点为黑色,则P’中的点Y’显示为黑色。如果每个条件都未满足,则Y’显示为白色。R’是侵蚀的结果图象。多个侵蚀和扩展的组合使得二值化图象更加规则(单一黑色和单一白色消失)。
在执行规则化之后,结果调查数据传递到步骤338用于提取关联成分。在该步骤中,对调查数据加以扫描,从而定义关联成分。对于调查数据中任何给定的黑色点对,如果两个点被它们之间一连串的黑点连接,则该对定义为处于同一成分。这一步骤的主要目的在于将调查数据分解成用白色空间分隔的关联黑色成分的集合。
图27显示了提取关联成分的分步。第一步骤350包括分配初始成分数字。以如下方式扫描调查数据,即给被扫描的第一黑色点分配“0”,下一个分配“1”,如此继续。所有白色点分配“-1”。在步骤352,设定初始扫描方向。有四个方向:
(1)“++”表示从上到下,从左到右,
(2)“+-”表示从上到下,从右到左,
(3)“-+”表示从下到上,从左到右,
(4)“--)表示从下到上,从右到左。
初始扫描方向设定为“++”,意味着从上到下,从左到右。在步骤354,对调查数据的扫描如下。具有如下黑色相邻点P’的每个黑色点P获得P’的分配数字,该黑色相邻点P’的分配数字(在步骤350中分配)小于P的分配数字。例如,如果P的数字为7,其相邻点P’的分配数字为6,则P的新数字为6。
接着改变确定步骤356。算法检查,看是否调查数据点中的任何黑色点都通过步骤354分配了新的数字。如果是,则在步骤358中改变扫描方向。方向的改变遵循如下法则:
(1)如果当前方向为“++”,则新方向为“+-”。
(2)如果当前方向为“+-”,则新方向为“-+”。
(3)如果当前方向为“-+”,则新方向为“--”。
(4)如果当前方向为“--”,则新方向为“++”。
在步骤354中再次以新方向开始扫描。该循环持续到扫描完成而在成分数目中检测不到任何改变时为止。根据成分数目变化改变扫描方向减小了处理所需的扫描经过次数。
步骤356之后,关联成分内的所有点都应当具有相同的数字。在步骤360,对成分中的点重新编号。这一步骤是需要的,因为一些数字可能已经消失了。例如,如果初始调查数据中有20个黑色点,那么每个点都会获得0-19的数字。如果扫描之后,发现有5个成分,则20个初始数字中有5个会消失。在5个数字为例如1,4,9,16和18的成分中将留下黑色的点。重新编号步骤将这5个成分中的点重新编号为0-4。大体上,调查数据中N个成分的编号(enumeration)从0到N-1。这便完成了关联成分提取的处理。
返回图25,在提取关联成分之后(步骤338)进行发现落入窗口内的成分的步骤(步骤340)。选择满足特定逻辑限制的最大的黑色成分。逻辑限制包括能够位于光盘上的窗口的数目以及窗口之间的近似距离。这便完成了发现窗口的处理,其将整个矩形选择处理(图23中)带到下一在窗口内切割标准矩形的步骤(步骤306)。
为了发现标准矩形,其具有(1)最大的总和照明(summaryillumination)和(2)位于相应于窗口的成分的一个点处的中心,包括该成分的区域的扫描如下。首先,算法沿着水平方向扫描并计算滑动平均值,如图26所示。然后算法沿着垂直方向扫描并计算滑动平均值。这包括产生如下段,其垂直地穿过调查数据阵列中的多个行。段的长度近似于窗口的高度。水平滑动平均值与垂直滑动平均值的最大平均值的交点显示为估计矩形的中心。该点具有最大值,因为它是窗口中最亮的点。将该点选择作为估计矩形的中心。一旦定义了该中心,在步骤312中,通过从中心点进行测量产生1×2mm(标准尺寸)的估计矩形。图28显示了在根据本发明实施例的软件显示上发现窗口之后,切割估计矩形的结果。
其它发现中心点的技术包括使用边缘跟踪(edge tracing)发现窗口,发现窗口的中心(使用点值),和通过检查代表调查数据的图象人工选择点。
步骤2b:在无窗口的光盘内选择估计矩形
选择估计矩形的第二个选择是图23中的步骤308,为无窗口的光生物盘选择。一个实施例使用光盘上的暗点定位估计矩形的期望位置。算法开始于步骤310,图23,发现用作用于标定光盘上样品不同区域的标记的暗点。暗点能够以如下的方式发现,该方式与结合图26讨论的在发现窗口中使用的滑动平均值相同。因为现在靶标是非常小的暗点(而不是窗口),所以发现暗点中使用的段要短得多。它们的长度接近于暗点的尺寸。然而,发现滑动平均值的操作原理相同。具有最低滑动平均值的段鉴别为暗点。
图29示出了暗点366的实例。在步骤312,图23,一旦发现了暗点,算法移离暗点以产生估计矩形。切割标准尺寸的估计矩形,其中该矩形具有定位在通过从所发现的暗点移离预先确定距离而发现的点处。图30显示了在从所发现的暗点366移离预先确定距离后切割的估计矩形368的实例。虚椭圆368、370和372鉴别了其它的细胞区域。在本发明的优选实施例中,估计矩形的位置信息能够嵌入在光盘上。代替发现暗点,系统能够从光盘阅读位置信息以定位用于放置估计矩形的区域。
步骤2c:选择估计矩形,包括用户选择的选项
选择估计矩形中的第三个且是最后一个选项涉及用户通过软件用户界面的输入,如图23的步骤316所示。在屏幕上,将根据调查数据生成的光生物盘的图象显示给用户,用户通过在图象上定义(defining)矩形选择估计矩形。在步骤316,确认用户所选的矩形是否大于标准尺寸。如果否,则意味着用户的矩形小于标准尺寸(步骤318,图23)。在这种情况下,在步骤322使用用户所选择的矩形进行计数。如果用户的矩形大于标准尺寸(步骤320,图23),在步骤324中切出标准尺寸估计矩形。
步骤3:背景照明均匀化
图31给出了图22中步骤248的更详细说明。选择出估计矩形之后,在以估计矩形为边界的区域内(称作“调查数据区域”)执行背景照明均匀化。这一步骤的主要目的是消除背景噪音,从而使背景更加均匀。为了实现它,在电执行(electrical implementation)中,背景照明均匀化使用软件算法模拟增益控制(gain control)的效果。
在图31的步骤380,选择相邻矩形(标准矩形内)的标准尺寸。注意相邻矩形不要与估计矩形混淆。相邻矩形设定在单一估计点的周围。其尺寸选择成比细胞足够地大,但受背景照明不均匀的影响足够地小。在其它化验中,能够替代细胞的调查特征是珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。从而,根据所进行化验的类型选择相邻矩形的尺寸。相邻矩形的尺寸为给定点P确定了在P附近有多少点被估计用于背景照明均匀化的处理。
在步骤382和384中分别进行垂直扫描和水平扫描以计算调查数据中每个点的K均值。
K以如下方式导出。首先,执行垂直扫描。在一个实施例中,首先计算调查数据区域中所有点的垂直平均值。点(x,y)的垂直平均值Kvert是(x,y-dy)-(x,y+dy)的范围中所有点的平均值。实际上,所有点列均沿着垂直方向加以扫描。dy项是相邻矩形的半高,该相邻矩形的尺寸在步骤380中确定。这里应用图26所描述的滑动平均值计算技术,除了现在的处理是沿着垂直方向进行之外。
一旦发现所有点的Kvert,便如下地执行水平扫描。对于点(x,y),该点的最终平均值通过获取(x-dx,y)-(x+dx,y)范围内所有点的Kvert给出。实际上,所有行均沿着水平方向加以扫描。dx项是相邻矩形的半宽,该相邻矩形的尺寸在步骤380中确定。总体结果是,对于特定点P,将P的相邻矩形内与P在同一行的所有预先计算的Kvert值加以平均以获得P的最终平均值。预先计算Kvert值减少了计算时间。计算时间不是与调查数据区域的尺寸乘以相邻矩形的尺寸成比例,而是计算时间仅与调查数据区域成比例。
继续参考图31,在步骤386,通过将每个点P的值V再分配给Vbackgroud+(V-Kneighbor),其中Kneighbor是以P的相邻矩形为边界的所有点的平均值。在一个实施例中,背景值Vbackgroud设定为整个调查数据区的平均值。在另一个实施例中,Vbackgroud设定为2,000。如果P的新值大于4,000,则使用4,000。如果它小于1,则使用1。先前为0的P值替换为2,000。该步之后,调查数据任何大区域内的平均值大约为2,000。换言之,整个背景接近一个具有任意背景值Vbackground的平面。图32显示了如背景照明均匀化之前软件所显示的调查数据,而图33显示了如背景照明均匀化之后软件所显示的调查数据。已成像的调查数据572,其是以图象格式提供的调查数据,表示了调查目标的数据。在图32和图33中,已成像的调查数据572标记或相应于被捕获的细胞。不使用滑动平均值,另一个实施例在背景照明均匀化步骤中使用傅立叶变换(FT)。在执行傅立叶变换时,调查数据首先转化为频域。然后,除去频域内的一部分频谱。这除去了一部分由光生物盘以及电路中的电噪音和其它不规则产生的背景噪音。在一个实施例中,除去了波长非常短或非常长的频谱(频率除以波长为1)。这些被除去波长的阈值经验地加以确定。最后,执行逆变换从而将数据还原成空间域。
步骤4:标准化
图34详细示出了图2的步骤250中所包含的步骤。如本文上面所概述的,需要进行标准化使标准偏差为600左右的数值,而调查数据的平均值为2,000左右的数值。标准化还使软件对硬件修改和调谐较不敏感。例如,探测电路,例如图10A的顶部探测器158,中的信号增益可以改变而不会显著影响结果细胞计数。
为了实现它,在一个实施例中,处理在步骤390开始,为来自上一步骤(背景照明均匀化)的结果调查数据计算平均值A和标准偏差S。忽略数值为0的点。在步骤392,图34,如下地对调查数据中的每一个点执行标准化。在一个实施例中,对于每个点P,P的值v替代为2,000+(v-A)*600/(S)。成分(v-A)使每个点集中,同时成分(600/S)调节振幅。结果与背景值2,000相加。数值600能够被调整以便获得期望的振幅范围。图36的数值曲线400显示了标准化之后曲线的一个实例。注意,数值曲线在2,000附近波动,波动振幅保持在600左右。
一旦P值被标准化,便如下地执行截取(步骤394):1)如果P的新值超过4,000,则数值截取为4,000;2)如果P的新值小于1,则数值截取为1。
在一个实施例中,软件的绘图用户界面显示调查数据中所有点的直方图,从而使用户看到标准化过程。图35显示了在逐步标准化期间软件所显示的示例调查数据的一部分。已成像调查数据572表示调查目标的数据。在图35中,已成像调查数据572表示或相应于被捕获细胞。输入框(input box)要求标准化的范围。如图所示,使用的数值为0-4,000。应当意识到,能够根据期望地使用任何范围的数值。
图36显示了标准化步骤之后的软件显示。顶部窗口396显示了一部分示例调查数据的近视图(close-up view)。已成像调查数据572表示被捕获细胞。底部曲线400显示了窗口396中水平虚线398所经过的点的相应数值。底部曲线400显示出,调查数据的背景区域被标准化(也就是具有小噪音的稳定值),其中含有细胞的区域在曲线上具有明显的“尖峰”。例如,在底部曲线400中,尖峰信号402相应于特定位置处的特定示例细胞404。尖峰信号402与背景噪音406显著不同。这简化了细胞识别的处理。
步骤5:过滤
图37详细示出了图22的步骤250中所包含的步骤。在图37的步骤410中,选择相邻矩形的尺寸。这些相邻矩形与背景照明均匀化步骤(图22的步骤248)中所使用的在概念上相似。然而,这一步骤中使用的相邻矩形大约为细胞的尺寸,且小于在背景照明均匀化中使用的相邻矩形。在步骤412,对调查数据中的每个点P,计算数值“足够区分”于Vbackground的P的附近点的数目。在一个实施例中,Vbackground的数值设定为2,000。通过阈值数目确定多大的差异(所调查点与Vbackground之间)构成“足够区分”。在一个实施例中,阈值数目通过检查调查数据的信号模式指出(note)背景数值与细胞(或者样品中目标的其它物质)的数值之间的差异加以确定。在其它化验中,调查特征除细胞之外,还能够是珠子(bead)(基于珠子的化验)、凝集物质、沉淀(酶反应)或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。在其它实施例中,阈值数目通过校准机制产生,该校准机制基于可变条件确定背景噪音和背景值。这些条件包括光生物盘的反射率,光生物盘驱动内光生物盘的不平衡,光生物盘的颤动(rattle)、震动(vibration)或不稳定,电噪音,光生物盘的金属化,所涉及样品的类型(白细胞或其它),和任何其它需要通过校准调节加以补偿或修正的条件。
在步骤414,检测“足够区分”点的数目,看是否该数目大于预先确定的过滤标准。如果是,则P值保持与步骤416一样。否则变为步骤418的Vbackground(或2,000)。这一步骤的期望效果是除去噪音,从而只有细胞保留在调查数据中并且背景均匀地等于Vbackground。图38显示了过滤步骤之后软件所显示的示例调查数据。已成像调查数据572标记或表示被捕获的细胞。注意,背景与调查数据中的细胞具有良好的对比度。图39提供了图39调查数据一部分的近视图。已成像调查数据572标记并相应于被捕获细胞。如图39的底部曲线420所示,现在背景具有平线值(flat line value),且含有细胞的区域被尖峰422良好地限定。对比度的改良有助于下一步骤的细胞计数。
步骤5a:除去不合需要的成分
图40详细显示了图22的步骤254所涉及的步骤。这是一个设定为消除不合需要的成分,例如可能干扰细胞计数的空气泡、灰尘和裂缝,的可选择步骤。这里使用的处理与结合图25说明的发现窗口中使用的处理相似。在步骤428,图40,选择阈值T。在一个实施例中,T设定为在背景照明均匀化步骤中发现的Vbackground。然后在步骤430,执行二值化。类似于在窗口发现(图25的步骤334)中使用的二值化步骤,在步骤430中,数值超过阈值T的点显示黑色,而其余的显示白色。从而现在,调查数据能够认为是黑与白的图象,并能够用1比特表示。在步骤432,以和图25的步骤336中相同的方法对调查数据进行规则化。
在执行规则化之后,结果调查数据传递到步骤434用于提取关联成分。在这一步骤中,扫描调查数据从而定义关联成分。对于任何给定的一对调查数据中的黑点,如果两个点能够用它们之间一连串的黑点连接,则该对定义为相同的成分。本步骤的主要目的是将调查数据分解成用白色空间分隔的关联黑色成分的集合。这里采用的关联成分提取处理可以与图27中详细显示的处理相同。
本去除处理的下一步是步骤436,其致力于除去不规则尺寸的成分。用户选择的尺寸阈值应用于所有的关联成分。如果成分小于或大于尺寸阈值,则从调查数据中除去整个成分。该方法是有效的,因为不规则成分(例如气泡,裂缝)的尺寸通常比典型的细胞尺寸大得多。同样,用户能够根据计数细胞的类型选择阈值。在一个实施例中,通过用具有恒定值2,000(背景值)的成分替代所有的点实现。优选地,在实现图40的所有步骤时,调查数据应当返回步骤248-步骤252(图22)进行重新处理。
图41显示了除去裂缝之前调查数据的实例。已成像调查数据572相应于被捕获细胞。如图所示,裂缝574遍布该区域。图42显示了除去裂缝后图41的相同调查数据。只保留标记被捕获细胞的已成像调查数据572。
步骤6:通过明亮中心计数细胞
现在参考图43,其详细显示了图22的步骤256中所涉及的步骤。在步骤440,图43,对调查数据进行卷积。在卷积期间,形成代表卷积图象的辅助阵列。卷积图象的每个点P是调查数据在P的圆形邻域中过滤后积分(integration)的结果。如本领域技术人员将意识到的,图象处理中使用的通用卷积方法包括两个函数。两个函数f和g:R2->R的卷积是函数:
F ( x , y ) = fog ( x , y ) = ∫ ∫ R 2 f ( x + u ) , y + v ) g ( u , v ) dudv .
更精确地,在本发明的一个实施例中,被积分的函数f是函数:
f(x,y)=h(x,y)-2,000  如果h(x,y)>2,000  或
      =0               如果h(x,y)<=2,000
其中h(x,y)是描述先前步骤中x,y点的数值的函数。一旦确定了f(x,y),便用圆形邻域(neighborhood)指示函数g进行卷积,其中:
g(u,v)=1                如果u2+v2≤r2,  或者
       =0                其它
其中r是细胞预期的半径。卷积积分是
F ( x , y ) = ∫ ∫ ( u - x ) 2 + ( v - y ) 2 ≤ r 2 f ( x + u , y + v ) dudv .
用所有格点(u,v)内f的总和值代替积分,格点处于如下限定的圆形邻域内:
(u-x)2+(v-y)2<r2
卷积之后,对步骤442,图43,中的卷积图象进行局部最大值的寻找。卷积步骤使调查数据中的明亮中心突出作为局部最大值,从而更容易识别。因为使用整数值,所以舍入(rounding)能够产生多余的局部最大值。为了修正它,在步骤444中除去位于同样紧邻的邻域内的冗余局部最大值。然后在步骤446中,所有保留的局部最大值显示为细胞中心。在其它化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物质、沉淀(酶反应),或其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。局部最大值用于发现这些目标或者这些化验中靶定的调查特征的中心。在本发明其它实施例中,计数方法考虑细胞凝集的效果。彼此接近的最大值不会被自动地忽略。在一个实施例中,局部峰在一个附近的下陷(nearby dip)之后被断言为细胞中心。参数能够加以调节,从而如果凝集细胞显示在调查数据上,限定局部冗余最大值的距离阈值能够被调节得更小。类似地,能够为正在计数的细胞类型调节距离阈值。例如,因为红细胞具有更一致的尺寸,所以假定的细胞尺寸能够用作距离阈值。
在另一个实施例中,能够对细胞的分布进行统计分析。从而每个区域上细胞的平均数目能够用于估计位于由于可见性低或者细胞凝集使它们不能计数的区域中的细胞数目。在其它化验中,细胞能够替换为珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。另一个实施例允许用户以更高的分辨率重新采样调查数据区域以执行更加完全而精确的计数。
图44显示了充满用明亮中心法计数的细胞的示例调查数据。已成像调查数据572标记被捕获细胞。在明亮中心法中执行的步骤有助于突出细胞,从而它们相对于背景的黑暗对比显得明亮。如软件所示,它们被单个地标记和计数。图45显示了图44所示一部分调查数据的近视图和数值迹线曲线。已成像调查数据572相应于被捕获细胞。
步骤7和8:细胞标记和通过黑暗边缘附加计数细胞
步骤7和8,图22,是能够执行用于提高细胞计数精度的可选择步骤。在图2中,它们参考步骤258和260。在硬件构型可能会剩下没有明亮中心的细胞的情况下,这两个可选择步骤可以用于修正细胞的欠计数。在其它化验中,本方法不是应用于细胞,而是能够应用于珠子(基于珠子的化验)、凝集物、沉淀(酶反应)或者其它尺寸能够通过本发明光学系统的入射波束检测到的生物指示标记。选择地,可以使用这两个步骤代替通过识别明亮中心计数细胞。
如果一些细胞已经通过识别明亮中心法加以计数了,那么执行步骤7标记这些已计数细胞并将它们从调查数据中除去。然后通过识别具有黑暗边缘的细胞进行计数。图46A详细示出了图22的步骤260(主要步骤8)(principal step 8)中所包含的步骤。在步骤450,对调查数据进行反转(inversion)。每个点P的数值v替换为2,000-v。如果结果数值为负,则用0代替。以等式形式表达反转,我们进行如下:
假定h(x,y)为以前的调查数据,我们引入f(x,y)使
f(x,y)=2,000-h(x,y)    如果h(x,y)<2,000
       =0                         其它
这保证在我们执行卷积时,具有低数据值的黑暗边缘具有较高的值。在步骤452,执行使用移动环(shifted ring)的卷积。本领域技术人员会意识到,在图象处理中使用的通用卷积方法包括两个函数。两个函数,f和g:R2->R,的卷积是函数
F ( x , y ) = fog ( x , y ) = ∫ ∫ R 2 ( x + u , y + v ) g ( u , v ) dudv
以等式形式表达卷积有:
F ( x , y ) = ∫ ∫ g ( u , v ) f ( x + u , y + v ) dudv
用圆形邻域指示函数(circular neighborhood indicator function)g执行卷积,其中
g(u,v)=1                如果u2+v2≤r2
       =0                其它
其中r是细胞的预期半径。在本实施例中,g是具有内径r1和外径r2的环的指示函数,其中r1和r2限定细胞预期半径r的界限。这产生:
F ( x , y ) = ∫ ∫ r 1 2 ≤ ( u - x ) 2 + ( v - y ) 2 ≤ r 2 2 f ( x + u , y + v ) dudv
用环内所有格点(u,v)的总和值f代替积分。进行四次卷积,我们得到四个函数:f1(x,y)=f(x+hx,y);
            f2(x,y)=f(x-hx,y);
            f3(x,y)=f(x,y+hy);和
            f4(x,y)=f(x,y-hy),
其中hx和hy是沿x和y方向的特定移动。它们等于细胞估计尺寸的一半。四个函数意味着用内径为r1,外径为r2的环的指示函数进行四次卷积。数值r1和r2分别限定细胞预期半径r的最小值和最大值。用沿左、右、上和下方向移动距离r的环进行四次卷积。图46B显示了这样的一个示例。首先产生环458,限定细胞黑暗边缘的边界。四次移动的卷积产生四个环。在步骤454,图46A,总和四次移动的结果。再次返回图46B,我们看到总和的环在点457产生局部最大值。然后断言点457为这一细胞的局部最大值并加以计数。注意,图46B是应用到细胞的一个示例。在卷积环不限定细胞黑暗边缘边界的位置,不存在最大值。这样,卷积步骤通过突出细胞的黑暗边缘发现潜在的细胞。在步骤456,图46A,在卷积计数局部最大值之后,计数步骤检查所有(go through)调查数据。
图47显示了调查数据的图象,在调查数据中通过该方法发现(计数)的细胞用十字标记。已成像调查数据572相应于被捕获的细胞。十字580标记已计数细胞。
选择地,能够根据如下等式形式执行卷积步骤:
F(x,y)=f1(x,y)+f2(x,y)+f3(x,y)+f4(x,y),如果f1、f2、f3或f4中至少有三个(或者选择地,两个)大于0,和
F(x,y)=0  其它。
注意,如果函数中有两个或三个大于0,则该卷积无移动环地加以执行。卷积步骤的另一个选择是使用众所周知的光滑函数来卷积函数f。在另外一个能够用于本卷积步骤和在识别明亮中心中使用的卷积步骤的选择中,采用不同的指示函数g。在其一个特殊实施例中,g能够是高斯形式:
g ( u , v ) = e - ( u 2 a 2 + v 2 b 2 )
或者其它用于执行卷积的合适函数,其中执行卷积的目的在于突出细胞的特征。
步骤9:数据输出
在步骤9,数据输出到合适的显示装置上。软件的一个实施例具有用户界面,以便为由估计矩形限定边界的调查数据区域显示细胞计数的结果。另一个实施例显示了调查数据区域的图象,该区域的每个细胞用十字标记,如图47所示。
B.红细胞实例
本领域技术人员会意识到,数据分析的各种步骤和方法能够以不同的方式组合用于分析各种类型的调查数据。图48提供了一个流程图,其示出了计数调查数据中红细胞的实例。在步骤460,选择阈值并进行二值化,其中数值超过阈值的点显示为黑色,其余的显示白色。二值化步骤将高数值的点与低数值的点分离开,高数值的点通常表示细胞,而低数值的点通常表示背景或背景噪音。
图49显示了执行步骤460之前调查数据的图示。已成像调查数据572相应于被捕获细胞。图50显示了二值化结果(步骤460)。二值化成像数据576标记或指示调查目标的数据。在这一例子中,二值化成像数据576指示细胞。非细胞标记二值化成像数据578指出了不代表细胞的数据。
然后在步骤462,图48,执行构成规则化的两个部分,侵蚀和扩展,以填充细胞边界的缺失部分。此处的目的是获得清晰划分出单个细胞的细胞边界。图51显示了侵蚀和扩展(步骤462)的结果。二值化成像数据576标记调查目标的数据。在这个实例中,二值化成像数据576标记细胞。非细胞标记的二值化成像数据578指出了不代表细胞的数据。
在图48的步骤464中,为每个细胞提取一个像素宽的细胞边界。图52显示了提取一个像素宽度边界(步骤464,图48)的结果。二值化成像数据576标记调查目标的数据。在这一实例中,二值化成像数据576标记细胞。非细胞标记的二值化成像数据578指出了不代表细胞的数据。在提取一个像素宽度的边界时,选择出所有同时具有黑色和白色邻域的黑色点使其保持黑色。相邻点不同时具有这两种颜色的黑色点被转变为白色。这产生了粗边界,其位于数个像素宽度的位置内。然后应用细化处理(thinning process)消除边界内多余的点,直到仅剩下一个像素显示调查数据中的轮廓。细化处理通过从边界上除去冗余的黑色点开始,直到只有一个像素宽度的线标记每个黑色区域的边界为止。如图52所示,在图51中存在的非细胞标记的二值化成像图象578现在不存在了。
细化处理首先以粗边界开始。粗边界包括所有同时具有位于细胞内部和外部的相邻点的像素。在粗边界提取之后,数据包括三类:
(1)细胞内的像素,
(2)标记细胞边界的像素,和
(3)细胞外的像素。
这三类根据三个条件相关。
(A)所有三类都相连,
(B)(1)通过(2)与(3)不相连,和最后
(C)(2)中的每个点都具有位于(1)或(2)内的相邻点,但不同时具有位于(1)和(2)内的相邻点。
然后细化处理逐个检查(2)内的像素。如果(2)内的像素P具有(1)或(3)内的相邻点,那么进行检查,看重新涂色P(例如由黑色变为白色)使之从(2)变到(1)或(3)是否仍保持(preserve)条件(A)、(B)和(C)。如果是,则进行重新涂色。执行该涂色直到为每个细胞获得一个像素宽的边界。
在提取一个像素宽的边界之后,被一个像素宽度的边界限定的区域用黑色点填充(步骤466,图48)。图53显示了步骤466的结果。二值化成像数据576标记细胞。非细胞标记的二值化成像数据578指出不代表细胞的数据,例如图51所示。使用这一系列黑色和白色点作为标记,填充原始数据点以代替黑色点。从而,细胞区域被隔离并且现在容易加以分析。
图54显示了填充原始数据点的结果。已成像调查数据572标记被捕获细胞。该方法的优点在于能够精确地提取细胞。该提取使用户能够测量细胞的直径,检查细胞内的其它特征,如被染色细胞核的形态。关于这类应用进一步的细节在共同转让(commonly assigned)美国专利申请序列No.10/xxx,xxx中有所讨论,其题目为“基于使用光生物盘系统鉴定和量化白细胞类型的细胞核形态学”,(NuclearMorphology Based Identification and Quantification of White BloodCell Types Using Optical Bio-Disc Systems),于2002年9月6日提出申请,在这里引用作为参考。
除了检查已成像细胞以外,用户能够通过采用本发明的用户特征(user feature)标记和计数细胞。图55是图49-54的调查数据原始图示的近视图,该调查数据具有用十字标记以显示它们已被计数的红细胞。已成像调查数据572相应于被捕获细胞。十字580标记已计数的细胞。
C.可选择算法
本发明包括大量用于处理在细胞计数操作期间可能出现的特殊情况的可选择算法。
图56A-64显示了本发明的实施例,其处理计数没有明显的明亮中心或黑暗边缘的细胞的情况。该方法称作“绝对值计数”,主要处理细胞看上去没有明显的明亮中心或黑暗边缘的情况。明亮中心法依赖于隔离调查数据的高数值区域(亮点)以计数细胞。黑暗边缘法依赖于隔离调查数据的低数值区域(暗点)以计数细胞。相反,该绝对值计数法隔离如下区域,该区域可以不具有由先前的一种或两种方法可探测的明显的高或低的数值,但是仍然含有能够区分于背景噪音的数值模式。
图56A是在它们根据绝对值计数法用十字标记之前,分散细胞的图示屏幕拍摄(pictorial screen shot)。已成像调查数据572代表被捕获细胞。图56B是描述绝对值计数所含步骤的流程图。在步骤480,对调查数据进行标准化和过滤。图57是最初在图56A中显示的分散细胞在标准化和过滤步骤(步骤480,图56B)之后的图示屏幕拍摄。已成像调查数572相应于被捕获细胞。标准化和过滤处理与本文上面所说明的处理相同。
在标准化和过滤之后,下一步骤包括背景去除和二值化(步骤482,图56B)。图58是最初在图56A中显示的分散细胞在背景去除和二值化步骤(步骤482,图56B)之后的图示屏幕拍摄。在该实例中,二值化成像数据576表示细胞。非细胞标记二值化成像图象578指出不代表细胞的数据。除去背景以隔离细胞定位的位置。然后对调查数据进行二值化从而在调查数据中产生黑色和白色的点。二值化处理的进行如下。首先检查每个点的数值。当点值与背景值之差大于预先确定的阈值时,点显示为黑色。否则点显示为白色。在一个实施例中,选择阈值使与背景值差异小的点(背景或背景噪音)成为白色,而与背景值差异大的点(细胞的黑暗或明亮区域)成为黑色。在其它实施例中,阈值通过校准机制产生,该校准机制根据可变条件确定背景噪音和背景值。这些条件包括光生物盘的反射率,光生物盘驱动内光生物盘的不平衡,光生物盘的颤动(rattle)、震动(vibration)或不稳定,电噪音,光生物盘的金属化,所涉及样品的类型(白细胞或其它),和任何其它需要补偿或修正的条件类型。
规则化(步骤484,图56B)是处理的下一步。图59是最初在图56A中显示的分散细胞在规则化步骤(步骤482,图56B)之后的图示屏幕拍摄。二值化成像数据576标记调查目标的数据。在该实例中,二值化成像数据576表示细胞。非细胞标记二值化成像图象578指出不代表细胞的数据。规则化包括侵蚀和扩展。侵蚀的执行如下。对于图象P,构建相应的图象P’。如果(1)P中相应的点X为白色,或者(2)X的任何邻域为白色,则P’中的点X’显示为白色。如果哪个条件都不满足,则X’显示为黑色。P’是侵蚀的结果图象。扩展操作是相反的形式。对于图象R,构建相应的图象R’。如果(1)R中相应的点Y为黑色,或者(2)Y的任何邻域为黑色,则R’中的点Y’显示为黑色。如果哪个条件都不满足,则Y’显示为白色。R’是扩展的结果图象。多个侵蚀和扩展的组合使得二值化图象更加规则(单一黑色和单一白色点消失)。
下一步骤是提取一个像素宽度的边界。该提取步骤参考图56B的步骤486。图60是一些最初在图56A中显示的分散细胞在应用一个像素宽度边界提取的步骤486之后的图示屏幕拍摄。二值化成像数据576标记调查目标的数据。在该实例中,二值化成像数据576相应于细胞。非细胞标记二值化成像图象578指出不代表细胞的数据。在提取一个宽度的边界中,选择所有同时具有黑色和白色邻域的黑色点保持其黑色。相邻点不同时具有两种颜色的黑色点转换为白色。这产生了粗边界,其位于数个像素宽度的位置内。然后应用细化处理(thinningprocess)消除边界内冗余的点,直到仅剩下一个像素显示调查数据中的轮廓。细化处理通过从边界上除去冗余的黑色点开始,直到只有一个像素宽度的线标记每个黑色区域的边界为止。
细化处理首先以粗边界开始。粗边界包括所有同时具有位于细胞内部和外部的相邻点的像素。在粗边界提取之后,数据包括三类:
(1)细胞内的像素,
(2)标记细胞边界的像素,和
(3)细胞外的像素。
这三类根据三个条件相关。
(A)所有三类都相连,
(B)(1)通过(2)与(3)不相连,和最后
(C)(2)中的每个点都具有位于(1)或(2)内的相邻点,但不同时具有位于(1)和(2)内的相邻点。
然后细化处理逐个检查(2)内的像素。如果(2)内的像素P具有(1)或(3)内的相邻点,那么进行检查,看从(2)到(1)或(3)的重新涂色P(例如由黑色变为白色)是否仍保留(preserve)条件(A)、(B)和(C)。如果是,则进行重新涂色。执行该涂色直到为每个细胞获得一个像素宽的边界。
在提取一个像素宽的边界之后,被一个像素宽度的边界限定的区域用黑色点填充(步骤488,图56B)。图61是一些最初在图56A中显示的分散细胞在执行根据本发明的填充成分步骤488之后的图示屏幕拍摄。二值化成像数据576标记调查目标的数据。在该实例中,二值化成像数据576标记细胞。非细胞标记的二值化成像数据578指出不代表细胞的数据。使用这一系列黑色和白色点作为标记,填充原始数据点以代替黑色点(步骤490,图56B)。从而,细胞区域被隔离并且能够加以计数。在一个实施例中,卷积应用于调查数据中被隔离的区域,并标记局部最大值以鉴定细胞。如先前所显示的实施例,能够应用具有圆形邻域的卷积。因为在卷积中使用的圆形邻域的尺寸大约为细胞的尺寸,所以局部最大值能够确定为细胞中心。图62是最初在图56A中显示的分散细胞在调查数据填充步骤之后的图示屏幕拍摄。成像调查数据572标记被捕获细胞。最后,图63显示了根据图56B的步骤492加以计数并用十字标记的细胞。成像调查数据572相应于被捕获细胞。十字580标记已计数的细胞。然而,图63只显示了凝集的细胞,图64显示了应用于稀疏群集(packed)样品的绝对值计数法,该样品具有单个的细胞和多个凝集细胞区域。类似地,成像调查数据572相应于被捕获细胞,且十字580标记已计数细胞。
本发明的另一个实施例改良了标记用于细胞计数的估计矩形的屏上显示。在细胞计数期间,有时需要向用户显示估计矩形的图象。除了提供样品的可视表示以外,高品质的图象能够帮助用户确定应当使用什么方法计数和分析细胞。例如,清晰的图象可以提醒用户存在许多没有明亮中心的细胞。这样,用户可以选择也通过黑暗边缘方法计数细胞。本实施例通过快速傅立叶变换法改良了图象的品质。如本公开所给出领域的技术人员会意识到的,快速傅立叶变换(FFT)是傅立叶变换(FT)的变型。傅立叶变换的任何变型都能够在这里应用。傅立叶变换在较早前结合背景照明均匀化时有所讨论。图65提供了示出本发明该实施例的流程图。在步骤520,对调查数据执行快速傅立叶变换。调查数据转换为频域。然后在步骤522,除去频域内的一部分频谱。最后,在步骤524,执行逆变换。图66显示了快速傅立叶变换之前的示例调查数据。已成像调查数据572表示被捕获细胞。图67显示了快速傅立叶变换之后的同一调查数据。已成像调查数据572表示被捕获细胞,且十字580标记已计数细胞。从而改良了屏上显示。
本发明的另一个实施例处理窗口区域的屏上显示。在具有窗口的光盘上进行细胞计数期间,有时需要向用户显示窗口区域。有时窗口区域的图象是歪斜的(skewed),如图68所示。为了适当地显示窗口区域,需要修正歪斜。修正方法的第一步发现歪斜的方向。该发现歪斜的步骤利用如下事实,即窗口是黑暗背景中的明亮区域——这方便地显示为白色。为了澄清这些项目,窗口是在其顶部和底部具有半圆形接触的矩形,且该矩形的宽度在下文称作窗口宽度。歪斜发现步骤进行如下。首先,数字地差分图象每条线的点。这意味着,对于每个点(x,y),从(x,y)到(x+dx,y)间隔内的平均值中减去(x-dx,y)到(x,y)间隔内所有点的平均值。这里dx是选择用于消除噪音的特殊间隔长度。
对于与窗口一致的图象的线,该减法的结果在窗口的左边界采用最大值,在右边界采用最小值。这是因为明亮窗口内的平均值大大高于黑暗背景的平均值。在位于窗口外的图象的线中,这最大值和最小值可以发生在别的任意位置,因为平均值来自于可能具有一些噪音的黑暗背景。利用这一性质,在下一步骤中,为每条线计算最大值和最小值之间的距离D。接着,该处理选择具有选择作为标准窗口宽度的D值的线。之后,在这些被选线中标记最大值点。最后,直线拟合通过这些最大值点。该直线的方向是歪斜的方向。图69显示了具有已发现方向线的结果。最后,使用该方向线作为导向并且据此除去所有相应点以修正歪斜。图象适当地加以对准并显示给用户,如图70所示。
本发明的一个实施例包括用于从调查数据中除去气泡痕迹的方法。有时空气泡可以截留在光盘上的沟道内。该空气泡可能通过样品,并且在其通过捕获区时除去一些沿着其路径的细胞。该细胞去除可能导致细胞分布的不均匀。因为报告的最终结果为每平方毫米区域(mm2)的细胞数目,所以不规则细胞必须加以修正。修正的执行如下。图71A显示了处理的流程图。在步骤540,如前地执行细胞计数。然后在步骤542分析细胞的分布。在步骤544忽略细胞局部浓度太小的区域,因为这些区域内的细胞可能已经被气泡擦除。
在一个实施例中,正在计数的整个区域划分成盒形的格子。图71B显示了这种实例。已成像调查数据572标记被捕获细胞(用圆表示)。气泡痕迹548通过由方形550、552、554和556限定边界的区域。忽略由这些方形限定边界的区域。一旦忽略这些区域,在图71A的步骤546中重新计算细胞计数。图71C和71D显示了气泡痕迹的实例。图71C显示了通过样品的气泡痕迹(包括痕迹548和558),如5X显微镜下所见。图71D显示了另一个实例,其中示出了40X显微镜下通过样品的气泡痕迹548。
III.白细胞计数方法
图72提供了如何利用本发明的方法执行用于快速确定CD4+和CD8+T-淋巴细胞种群绝对数值的类均匀固相细胞捕获化验(generichomogeneous solid phase cell capture assay)。该检测在整合于光生物盘上的小流动沟道内执行,确定7-15μl从全血中分离的单核细胞(MNC)中被捕获区内特异抗体捕获的CD4+、CD8+CD2+、CD3+、CD19+和CD45+细胞的数目。该检测根据在光盘的局部位置上进行特定细胞捕获的原理。在光盘上通过基于单克隆或多克隆特殊血细胞表面抗原的抗体局部应用捕获化学剂产生多个特殊的细胞捕获区。在充满25-100μl具有MNC血液(10,000-30,000细胞/μl)的室时,表达CD4、CD8、CD2、CD3、CD19和CD45抗原的细胞被捕获在光盘的捕获区内。同时在捕获区内限定阴性和阳性控制区域(negative andpositive control areas)。
在图72的步骤1中,血液(4-8ml)和抗凝血剂如EDTA、ACD或肝素直接收集到4或8ml的Becton Dickinson CPT VacutainerTM中。在本发明另一实施例的等价步骤中,3ml含有抗凝血剂的血液覆盖含有分离梯度(separation gradient)176如Histopaque 1077的管172。在任何情况下,血液样品174优选地用在2小时以内的收集中。血液样品所覆盖的含有分离梯度176的管172在具有水平马达的生物危害物质离心机(biohazard centrifuge)中在1,500-1,800RCF(2,800rpm)下离心。离心之后,除去血浆层178(步骤2),在单核细胞(MNC)部分180上方留下大约2mm的血浆。收集并用磷酸盐缓冲液(PBS)冲洗MNC层180。细胞在室温下以300RCF(1200rpm)离心10分钟成为小球(pelleted),从而除去任何残留的血小板。除去上清液,并通过轻轻敲打管使MNC小球180在PBS中重新悬浮。根据光生物盘110流动沟道130的高度,将最终的小球180(步骤3)重新悬浮为10,000-30,000细胞/μl的细胞数。
光生物盘110的流动沟道130充满7μl的MNC悬液,且室的入口122和出口124(图3和5)用密封标签密封(步骤4)。光生物盘110在室温下培育15分钟,然后用光驱动器112内的780nm激光器扫描以成像捕获区域(步骤5)。应当理解,如果使用透射光生物盘110,光驱动器112选择地包括顶部探测器158(图10A)以成像捕获区域。软件优选地编码在光盘上指导驱动自动地执行如下操作:(a)离心光盘从而在一个或多个阶段中旋出多余的未键合细胞,(b)在显示器114上成像特殊捕获窗口,和(c)处理数据。数据处理包括,但不仅限于,计数每个捕获区内特异捕获的细胞并得出CD4+/CD8+的比率或任何其它可以相应编程的期望计数或定量。
如图72进一步示出的,本发明致力于用光盘和光驱动器执行聚类指示计数(cluster designation count)的方法。该方法包括如下步骤:在含有分离梯度的第一管内提供血液样品,在某时旋转第一管且速度足以使血液样品分离成层,重新悬浮含有T-细胞的MNC层形成MNC悬液,在包括至少一个含有至少一种捕获剂的捕获区的光盘表面提供MNC悬液样品,将光盘装载到光学阅读器内,旋转光盘,将电磁辐射的入射波束导向到捕获区,探测在与光盘捕获区相互作用之后形成的电磁辐射波束,将探测到的波束转换为输出信号,和分析输出信号以提取有关捕获区内捕获的细胞数目的信息。在本发明的一个实施例中,光盘用反射层加以构建,从而反射被导向到捕获区的并与细胞相互作用的光线。在本发明的另一个实施例中,构建光盘使得被导向到捕获区的并与细胞相互租用的光线透过光盘。
在分析/处理步骤期间,软件通读每个捕获区图象并在其遇到细胞图象时标记细胞图象。例如,在估计CD4+和CD8+细胞之后,软件计算CD4+/CD8+细胞的比率,并同时显示CD4+、CD8+、CD3+和CD45+捕获区内每毫升全血的细胞绝对数目和CD4+/CD8+的比率。整个处理从将光盘插入光驱动器内到获得数目和比率花费大约12分钟。
在一个实施例中,光盘是前摆动系列FDL21:13707或FDL21:1270 CD-R光盘,其涂覆有300nm的金作为编码信息层。在反射光盘上,通过众所周知的平版印刷技术从反射层腐蚀掉尺寸为2×1mm的观察窗口。在透射光盘的一些设计中,不腐蚀单独的观察窗口,且整个光盘都可以使用。附加层是Fraylock粘合剂DBL 201 Rev C3M94661。顶盖是具有48个直径为0.040英寸的样品入口的干净光盘,样品入口等距离地分布在26mm的半径上。使用CD4+/CD8+计数软件以4X速度和2.67MHz的取样速度扫描并阅读数据光盘。
IV.结论
这样,结合一个或多个特殊的实施例说明了用于成像实验室样品中的细胞和分析该图象的方法和装置。尽管参考某些优选实施例详细地说明了本发明,但是应当意识到,本发明并不仅限于那些确定的实施例。而是,由于说明了用于实践本发明的当前最佳模式的本公开,许多修改和变化将自身提供给本领域的技术人员而不背离本发明的范围和精神。因此本发明的范围由下文的权利要求指明而不是前述的说明。在权利要求等价物的意思和范围内出现的所有改变、修改和变化均认为包含在其范围之内。

Claims (56)

1.一种对细胞进行计数的方法,该方法包括如下步骤:
获得含细胞样品的调查数据;
在所述调查数据中选择估计矩形;
在所述估计矩形内强化所述调查数据;和
对所述估计矩形内的细胞进行计数。
2.权利要求1的方法,其中所述选择步骤进一步包括对所述估计矩形选择自定义尺寸的步骤。
3.权利要求1的方法,其中所述选择步骤选择多个估计矩形。
4.权利要求1的方法,其中所述强化所述调查数据区域的步骤进一步包括如下步骤:
对所述调查数据执行背景照明均匀化;
对所述调查数据执行标准化;和
过滤所述调查数据。
5.权利要求4的方法,其中所述执行背景照明均匀化的步骤进一步包括如下步骤:
为相邻矩形选择尺寸;
在所述调查数据中选取一个点;
执行水平扫描从而为将中心定于所述点的所述相邻矩形内的所有相邻点计算第一滑动平均值;
执行垂直扫描从而为将中心定于所述点的所述相邻矩形内的所有相邻点计算第二滑动平均值;
组合所述第一滑动平均值与第二滑动平均值产生整体平均值;
将所述点的原始值重新分配给一个结果值,该结果值通过获得所述整体平均值与所述原始值之间的差值并将所述差值加到背景值而计算出来;和
对所述调查数据中的所有点重复所述执行水平扫描、执行垂直扫描、组合两个所述平均值和重新分配原始值的步骤。
6.权利要求4的方法,其中所述执行背景照明均匀化的步骤进一步包括如下步骤:
对所述调查数据执行傅立叶变换以产生频域函数;
从所述频域函数中除去低波长函数;
从所述频域函数中除去高波长函数;和
对所述频域函数执行逆变换以获得所述调查数据的修改版本。
7.权利要求4的方法,其中所述执行标准化的步骤进一步包括如下步骤:
计算所述调查数据中所有点数值的平均值和标准偏差;
使用所述平均值和标准偏差使所述调查数据中的所有点的数值标准化;和
如果需要则截去一些点的所述数值。
8.权利要求4的方法,其中所述过滤步骤进一步包括如下步骤:
为相邻矩形选择一个尺寸;
在所述调查数据中选取一个点;
发现位于将中心定于所述点的所述相邻矩形内的所有足够区别的点;
如果所述足够区别的点的数目大于预先确定的过滤标准则重新分配所述点的数值;和
重复所述发现所有足够区别的点和为所述调查数据中的所有点重新分配数值的步骤。
9.权利要求4的方法,进一步包括步骤:
在所述过滤步骤之后从所述调查数据中除去不合需要的成分;和
重复所述执行背景照明的步骤、所述执行标准化的步骤和所述过滤步骤。
10.权利要求9的方法,其中所述除去不合需要成分的步骤进一步包括如下步骤:
选择阈值;
使用该阈值对所述调查数据执行二值化;
对所述调查数据执行规则化;
提取关联成分;
选择尺寸阈值;和
除去不满足所述尺寸阈值的成分。
11.权利要求10的方法,其中所述执行规则化的步骤进一步包括执行多个侵蚀和扩展的步骤。
12.权利要求10的方法,其中所述提取关联成分的步骤进一步包括如下步骤:
对所述调查数据上的所有黑色点分配初始成分数目;
选取开始点;
设定扫描方向;
扫描所述调查数据的所有点以重新分配每个所述黑色点的成分数目从而与相邻黑色点的成分数目相匹配;
根据一组预先确定的规则改变所述扫描方向;和
重复所述扫描和改变的步骤从而使关联的黑色点的所述成分数目相同。
13.权利要求1的方法,其中对所述显示于所述估计矩形内的细胞进行计数的步骤进一步包括如下步骤:
对所述调查数据执行卷积;
寻找所述调查数据的多个局部最大值;
从所述多个局部最大值中除去冗余的局部最大值;
断言剩余的最大值为细胞的明亮中心;和
通过识别所述细胞的明亮中心对细胞进行计数。
14.权利要求13的方法,其中所述执行卷积的步骤使用定义了一个圆形邻域的指示函数,其中所述圆形邻域限定了细胞预期尺寸的边界。
15.权利要求13的方法,其中所述执行卷积的步骤使用高斯指示函数。
16.权利要求13的方法,其中所述除去冗余局部最大值的步骤进一步包括如下步骤:
选择距离阈值;和
使用所述距离阈值以确定一个局部最大值是否是冗余的。
17.权利要求13的方法,进一步包括执行统计分析的步骤,其包括如下步骤:
根据已计数细胞获得细胞的分布;和
估计细胞凝集或可见性低的区域内的细胞数量。
18.权利要求13的方法,进一步包括如下步骤:
以更高的分辨率重新取样所述调查数据;和
重复所述执行卷积、寻找多个局部最大值、除去冗余局部最大值、断言剩余的最大值为细胞的明亮中心和通过识别所述细胞的明亮中心对细胞进行计数的步骤。
19.权利要求13的方法,进一步包括如下步骤:
从所述调查数据中除去所述通过明亮中心计数的细胞;
通过识别黑暗边缘对细胞进行计数;和
将从通过识别明亮中心对细胞进行计数的步骤获得的总数和从通过识别黑暗边缘而计数的步骤获得的总数相加。
20.权利要求19的方法,其中所述通过识别黑暗边缘对细胞进行计数的步骤进一步包括如下步骤:
对所述调查数据执行反转;
执行多次使用移动环的卷积;
总和来自所述多次卷积的结果;
发现局部最大值;
断言最大值为细胞的中心;和
对所述细胞的中心进行计数。
21.权利要求20的方法,其中所述执行多次卷积的步骤无移动环地执行卷积。
22.权利要求20的方法,其中所述执行卷积的步骤使用高斯指示函数。
23.权利要求20的方法,其中所述执行卷积的步骤使用光滑函数。
24.权利要求1的方法,其中所述对显示于所述估计矩形内的细胞进行计数的步骤进一步包括如下步骤:
对所述调查数据执行反转;
执行多次使用移动环的卷积;
总和来自所述多次卷积的结果;
发现局部最大值;
断言最大值为细胞的中心;和
对所述细胞的中心进行计数。
25.权利要求24的方法,其中所述执行多次卷积的步骤无移动环地执行卷积。
26.权利要求1的方法,其中所述强化步骤进一步包括如下步骤:
对所述调查数据执行标准化;
对所述调查步骤执行过滤;
选择阈值数目;
通过确定所述调查数据是否与一组背景值相差大于所述阈值数目的数值而对所述调查数据执行二值化;
对所述调查数据执行规则化;
在所述调查数据中提取一个像素宽度的边界;
用调查数据填充由所述一个像素宽度的边界限定的区域;和
在所述区域填充中应用卷积。
27.权利要求1的方法,进一步包括在计算机监视器上显示表示所述调查数据的图象的步骤。
28.权利要求27的方法,其中所述显示步骤进一步包括如下步骤:
对所述调查数据进行快速傅立叶变换以在频域中产生调查数据;
除去频域内的一部分频谱;和
对频域内的所述调查数据执行逆变换以强化用于显示的所述调查数据。
29.权利要求1的方法,其中所述获得含有细胞样品的调查数据的步骤包括如下步骤:
在光盘表面上提供血液样品,所述表面包括一个或多个具有一种或多种捕获剂的捕获区;
将所述光盘装载到光学阅读器内;
旋转所述光盘;
将电磁辐射的入射波束从光源导向到所述捕获区的其中一个;
用探测器探测电磁辐射的结果波束,该结果波束在所述入射波束与光盘在所述捕获区相互作用之后形成;
将探测到的波束转换成模拟输出信号;和
将所述模拟输出信号转换成含有在所述捕获区捕获的细胞的数字数据。
30.权利要求29的方法,其中所述将所述模拟输出转换成所述数字数据的步骤进一步包括如下步骤:
以固定的间隔取样所述模拟信号的幅度;
将所述取样振幅记录在一维阵列内;
使用所述取样和记录步骤产生多个一维阵列;和
组合所述多个一维阵列产生含有所述样品的数字数据的二维阵列。
31.根据权利要求29的方法,其中所述光盘用反射层构建,从而导向到所述捕获区的光线被反射到所述探测器。
32.权利要求31的方法,其中所述探测器是底部探测器。
33.根据权利要求29的方法,其中光盘构建成使导向到所述捕获区的光线透射通过所述光盘,所述光盘在所述光源与所述探测器之间。
34.权利要求33的方法,其中所述探测器是顶部探测器。
35.权利要求33的方法,其中所述探测器是分波探测器。
36.权利要求29的方法,其中所述一个或多个捕获区定位在所述光盘内的一个或多个室内。
37.权利要求29的方法,其中所述光盘包括多个相应于所述捕获区的窗口。
38.权利要求37的方法,其中所述选择估计矩形的步骤进一步包括如下步骤:
在所述调查数据中发现所述多个窗口的其中一个;和
在所述窗口内修剪标准尺寸的估计矩形。
39.权利要求37的方法,其中发现所述多个窗口的其中一个的步骤进一步包括如下步骤:
对所述调查数据执行压缩;
对所述调查数据执行阈值估计;
对所述调查数据执行二值化;
对所述调查数据执行规则化;
从所述调查数据中提取关联成分;和
从所述关联成分中发现对应于一个窗口的成分。
40.权利要求39的方法,其中所述提取关联成分的步骤进一步包括如下步骤:
为所述调查数据上的所有黑色点指定初始成分数目;
选取开始点;
设定扫描方向;
扫描所述调查数据的所有点以重新分配每个所述黑色点的成分数目以便与相邻黑色点的成分数目相匹配;
根据一组预先确定的规则改变所述扫描方向;和
重复所述扫描和改变的步骤从而使关联的黑色点的所述成分数目相同。
41.权利要求29的方法,其中所述光盘的表面含有标记所述捕获区位置的暗点。
42.权利要求41的方法,其中所述选择估计矩形的步骤进一步包括:
在所述调查数据中发现所述暗点的其中一个;和
产生标准尺寸的估计矩形,中心位于通过从所述暗点移动预先确定的距离发现的点处。
43.权利要求42的方法,其中所述发现所述暗点的其中一个的步骤进一步包括如下步骤:
对所述调查数据进行压缩;
对所述调查数据进行阈值估计;
对所述调查数据进行二值化;
对所述调查数据进行规则化;
从所述调查数据提取关联成分;和
从所述关联成分中发现相应于暗点的成分。
44.权利要求43的方法,其中所述提取关联成分的步骤进一步包括如下步骤:
为所述调查数据上的所有黑色点指定初始成分数目;
选取开始点;
设定扫描方向;
扫描所述调查数据的所有点以重新分配每个所述黑色点的成分数目以便与相邻黑色点的成分数目相匹配;
根据一组预先确定的规则改变所述扫描方向;和
重复所述扫描和改变的步骤从而使关联的黑色点的所述成分数目相同。
45.权利要求42的方法,其中所述发现所述暗点其中一个的步骤进一步包括读取位置信息的步骤。
46.权利要求29的方法,其中所述光盘含有用于定位所述捕获区的计算机可读的位置信息。
47.权利要求37的方法,进一步包括在计算机监视器上显示所述窗口的图象的步骤。
48.权利要求47的方法,其中所述显示所述窗口的所述图象的步骤进一步包括:
确定所述图象是否歪斜;
发现歪斜的方向;和
修正所述图象的歪斜。
49.权利要求1的方法,其中所述获得所述含有细胞样品的调查数据的步骤包括从存档中检索先前存储的样品的调查数据。
50.权利要求49的方法,其中所述存档根据患者的特征对所述存储的调查数据进行分类。
51.权利要求50的方法,其中所述检索先前存储的样品调查数据的步骤进一步包括如下步骤:选择与多个从所述患者的特征中选则出的标准相匹配的样品从而进行人口健康趋势分析。
52.权利要求1的方法,进一步包括输出来自所述对细胞进行计数步骤的结果的步骤。
53.权利要求52的方法,其中所述细胞是白细胞。
54.权利要求53的方法,其中所述结果包括CD4+细胞和CD8+细胞的计数,和CD4+与CD8+细胞的比率。
55.权利要求54的方法,其中所述结果进一步包括CD3+细胞和CD45+细胞的计数。
56.权利要求1的方法,其中所述对细胞进行计数的步骤进一步包括如下步骤:
对气泡痕迹分析细胞分布;
忽略具有过小的局部细胞浓度的区域;和
重新计算细胞的计数。
CN02820990.7A 2001-09-12 2002-09-11 用于差分细胞计数的方法以及用于执行该方法的相关装置和软件 Pending CN1575475A (zh)

Applications Claiming Priority (16)

Application Number Priority Date Filing Date Title
US32286301P 2001-09-12 2001-09-12
US60/322,863 2001-09-12
US35330002P 2002-01-31 2002-01-31
US35392102P 2002-01-31 2002-01-31
US60/353,300 2002-01-31
US60/353,921 2002-01-31
US35564402P 2002-02-05 2002-02-05
US60/355,644 2002-02-05
US35530402P 2002-02-08 2002-02-08
US60/355,304 2002-02-08
US35847902P 2002-02-19 2002-02-19
US60/358,479 2002-02-19
US36394902P 2002-03-12 2002-03-12
US60/363,949 2002-03-12
US40492102P 2002-08-21 2002-08-21
US60/404,921 2002-08-21

Publications (1)

Publication Number Publication Date
CN1575475A true CN1575475A (zh) 2005-02-02

Family

ID=27575386

Family Applications (1)

Application Number Title Priority Date Filing Date
CN02820990.7A Pending CN1575475A (zh) 2001-09-12 2002-09-11 用于差分细胞计数的方法以及用于执行该方法的相关装置和软件

Country Status (6)

Country Link
US (1) US20030096324A1 (zh)
EP (1) EP1425696A2 (zh)
JP (1) JP2005502369A (zh)
CN (1) CN1575475A (zh)
AU (1) AU2002333589A1 (zh)
WO (1) WO2003023571A2 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102472703A (zh) * 2009-10-30 2012-05-23 西门子公司 体液分析系统以及用于分析体液的图像处理设备和方法
CN104508481A (zh) * 2012-06-01 2015-04-08 原子能和替代能源委员会 用于对诸如血液颗粒之类的包含在流体中的颗粒的移动速率进行表征的方法和系统
CN106332551A (zh) * 2014-03-31 2017-01-11 Jvc 建伍株式会社 分析装置和分析方法
CN106644897A (zh) * 2016-10-14 2017-05-10 北京海岸鸿蒙标准物质技术有限责任公司 一种用于颗粒计数标准物质的计数装置
CN108169081A (zh) * 2017-12-14 2018-06-15 四川大学华西医院 血细胞分析的差值核查模型及其应用方法
CN109085332A (zh) * 2018-09-29 2018-12-25 厦门大学 一种动态数据采集装置
CN111539354A (zh) * 2020-04-27 2020-08-14 易普森智慧健康科技(深圳)有限公司 一种液基细胞玻片扫描区域识别方法
CN111832389A (zh) * 2020-05-25 2020-10-27 中国人民解放军陆军军医大学第二附属医院 一种骨髓细胞形态学自动检测系统的计数及分析方法

Families Citing this family (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2518677A1 (en) * 2003-03-03 2004-12-09 Nagaoka & Co., Ltd. Methods and apparatus for use in detection and quantitation of various cell types and use of optical bio-disc for performing same
GB2429058B (en) * 2004-03-06 2008-12-03 Michael Trainer Method and apparatus for determining the size and shape of particles
US9297737B2 (en) 2004-03-06 2016-03-29 Michael Trainer Methods and apparatus for determining characteristics of particles
WO2006017642A2 (en) * 2004-08-04 2006-02-16 Agilent Technologies, Inc. Filtering of signals during array scanning
EP1825317B1 (en) * 2004-11-24 2013-04-17 Battelle Memorial Institute Optical system for cell imaging
JP2007080327A (ja) * 2005-09-13 2007-03-29 Hitachi Ltd 記録媒体欠陥検査装置及び欠陥検査方法
US7518723B2 (en) * 2005-09-19 2009-04-14 Jmar Technologies, Inc. Systems and methods for detecting radiation, biotoxin, chemical, and biological warfare agents using a multiple angle light scattering (MALS) instrument
US7504629B2 (en) * 2007-03-07 2009-03-17 Canon Kabushiki Kaisha Image forming apparatus
WO2009102159A2 (ko) * 2008-02-14 2009-08-20 Jae Chern Yoo 바이오 디스크 판독 장치 및 이를 이용한 분석 방법
EP2593771B1 (en) * 2010-07-16 2019-09-04 Luminex Corporation Methods, storage mediums, and systems for analyzing particle quantity and distribution within an imaging region of an assay analysis system and for evaluating the performance of a focusing routing performed on an assay analysis system
WO2012043499A1 (ja) 2010-09-30 2012-04-05 日本電気株式会社 情報処理装置、情報処理システム、情報処理方法、プログラム及び記録媒体
US10114020B2 (en) 2010-10-11 2018-10-30 Mbio Diagnostics, Inc. System and device for analyzing a fluidic sample
JP5776433B2 (ja) * 2011-08-11 2015-09-09 株式会社Jvcケンウッド 光学的分析装置、光学的分析方法
US9506935B2 (en) 2012-06-01 2016-11-29 Commissariat à l'énergie atomique et aux énergies alternatives Method and system for estimating the quantity of an analyte contained in a liquid
WO2014070235A1 (en) 2012-10-29 2014-05-08 Mbio Diagnostics, Inc. Biological particle identification system, cartridge and associated methods
US9176112B2 (en) 2012-12-31 2015-11-03 Beckman Coulter, Inc. Systems and methods for platelet count with clump adjustment
JP6449339B2 (ja) * 2015-01-29 2019-01-09 オリンパス株式会社 細胞解析装置および方法
WO2016133915A1 (en) * 2015-02-17 2016-08-25 Siemens Healthcare Diagnostics Inc. Classification of barcode tag conditions from top view sample tube images for laboratory automation
DE112015006266T5 (de) 2015-04-07 2017-11-30 Olympus Corporation Zellanalysegerät und Zellanalyseverfahren
CN104951782A (zh) * 2015-07-01 2015-09-30 携程计算机技术(上海)有限公司 图像识别的背景过滤方法及系统
CN106085829B (zh) * 2016-07-06 2018-02-16 上海海洋大学 一种卵母细胞集聚器
US10809178B2 (en) * 2017-02-22 2020-10-20 JVC Kenwood Corporation Analysis device and analysis method
JP6958166B2 (ja) * 2017-02-22 2021-11-02 株式会社Jvcケンウッド 分析装置及び分析方法
CN107123102A (zh) * 2017-05-24 2017-09-01 天津工业大学 一种贴壁细胞生长融合度自动分析方法
JPWO2021100191A1 (zh) 2019-11-22 2021-05-27
CN113554690B (zh) * 2021-07-20 2024-04-05 北京邮电大学 一种基于归一化卷积的细胞膜等效半径识别方法
JP7760882B2 (ja) * 2021-09-30 2025-10-28 株式会社プロテリアル 窒化珪素基板の評価方法、窒化珪素基板の評価装置、窒化珪素基板の評価システム、窒化珪素基板の製造方法、及びパワーモジュールの製造方法

Family Cites Families (48)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3736432A (en) * 1971-03-22 1973-05-29 Varian Associates Bacterial colony counting method and apparatus
US3710323A (en) * 1971-12-13 1973-01-09 Ibm Pattern-size normalizing for recognition apparatus
CA923621A (en) * 1972-10-16 1973-03-27 Burroughs Corporation Multi-stage queuer system
US3966322A (en) * 1973-11-08 1976-06-29 Vickers Limited Device for use in producing a scanning beam of radiation and apparatus for use in investigating specimens
US4199748A (en) * 1976-11-01 1980-04-22 Rush-Presbyterian-St. Luke's Medical Center Automated method and apparatus for classification of cells with application to the diagnosis of anemia
US4209548A (en) * 1976-11-01 1980-06-24 Rush-Presbyterian-St. Luke's Medical Center Method for the preparation of blood samples for automated analysis
JPS5952359A (ja) * 1982-09-02 1984-03-26 Hitachi Medical Corp 画像間演算時の画像歪み自動補正装置
FR2535058B1 (fr) * 1982-10-21 1987-08-21 Materiel Biomedical Dispositif pour la detection et la quantification d'agglutinats
FR2545610B1 (fr) * 1983-05-02 1989-04-21 Materiel Biomedical Procede et dispositif pour la detection et la quantification d'agglutinats
US5112134A (en) * 1984-03-01 1992-05-12 Molecular Devices Corporation Single source multi-site photometric measurement system
US4790640A (en) * 1985-10-11 1988-12-13 Nason Frederic L Laboratory slide
US5087556A (en) * 1989-05-17 1992-02-11 Actimed Laboratories, Inc. Method for quantitative analysis of body fluid constituents
US5143854A (en) * 1989-06-07 1992-09-01 Affymax Technologies N.V. Large scale photolithographic solid phase synthesis of polypeptides and receptor binding screening thereof
US5107422A (en) * 1989-10-02 1992-04-21 Kamentsky Louis A Method and apparatus for measuring multiple optical properties of biological specimens
JP3036049B2 (ja) * 1990-10-31 2000-04-24 スズキ株式会社 粒子凝集パターン判定方法
US5548661A (en) * 1991-07-12 1996-08-20 Price; Jeffrey H. Operator independent image cytometer
US6192320B1 (en) * 1991-07-30 2001-02-20 The University Of Virginia Patent Foundation Interactive remote sample analysis system
US5366896A (en) * 1991-07-30 1994-11-22 University Of Virginia Alumni Patents Foundation Robotically operated laboratory system
USH1183H (en) * 1991-12-19 1993-05-04 Laser scanner method for determining number and size of particles
JPH05215750A (ja) * 1992-02-06 1993-08-24 Idemitsu Petrochem Co Ltd 免疫学的分析方法
US5329461A (en) * 1992-07-23 1994-07-12 Acrogen, Inc. Digital analyte detection system
AU2593192A (en) * 1992-09-14 1994-04-12 Oystein Fodstad Detection of specific target cells in specialized or mixed cell population and solutions containing mixed cell populations
US5747265A (en) * 1992-10-30 1998-05-05 T Cell Diagnostics, Inc. Method for measuring the amount of a cell-associated molecule
US5478750A (en) * 1993-03-31 1995-12-26 Abaxis, Inc. Methods for photometric analysis
US5587833A (en) * 1993-07-09 1996-12-24 Compucyte Corporation Computerized microscope specimen encoder
US5793969A (en) * 1993-07-09 1998-08-11 Neopath, Inc. Network review and analysis of computer encoded slides
US5537669A (en) * 1993-09-30 1996-07-16 Kla Instruments Corporation Inspection method and apparatus for the inspection of either random or repeating patterns
EP0688113A2 (en) * 1994-06-13 1995-12-20 Sony Corporation Method and apparatus for encoding and decoding digital audio signals and apparatus for recording digital audio
US6143510A (en) * 1994-07-29 2000-11-07 Iatron Laboratories Inc. Measuring method using whole blood sample
US5978497A (en) * 1994-09-20 1999-11-02 Neopath, Inc. Apparatus for the identification of free-lying cells
US6327031B1 (en) * 1998-09-18 2001-12-04 Burstein Technologies, Inc. Apparatus and semi-reflective optical system for carrying out analysis of samples
GB9418981D0 (en) * 1994-09-21 1994-11-09 Univ Glasgow Apparatus and method for carrying out analysis of samples
US5795716A (en) * 1994-10-21 1998-08-18 Chee; Mark S. Computer-aided visualization and analysis system for sequence evaluation
US5585069A (en) * 1994-11-10 1996-12-17 David Sarnoff Research Center, Inc. Partitioned microelectronic and fluidic device array for clinical diagnostics and chemical synthesis
US6143247A (en) * 1996-12-20 2000-11-07 Gamera Bioscience Inc. Affinity binding-based system for detecting particulates in a fluid
US5871697A (en) * 1995-10-24 1999-02-16 Curagen Corporation Method and apparatus for identifying, classifying, or quantifying DNA sequences in a sample without sequencing
US5741213A (en) * 1995-10-25 1998-04-21 Toa Medical Electronics Co., Ltd. Apparatus for analyzing blood
JPH09173300A (ja) * 1995-10-25 1997-07-08 Toa Medical Electronics Co Ltd 血液分析装置
US5798994A (en) * 1996-01-16 1998-08-25 Kamatani; Yasuo Multi-layered optical disk reading method using liquid crystal diffraction device
JP3996644B2 (ja) * 1996-06-07 2007-10-24 イムニベスト・コーポレイション 外部及び内部勾配を有する磁気分離
JP3679512B2 (ja) * 1996-07-05 2005-08-03 キヤノン株式会社 画像抽出装置および方法
CN1249816A (zh) * 1997-02-28 2000-04-05 伯斯坦恩实验室股份有限公司 盘片中的实验室
JPH10275150A (ja) * 1997-03-28 1998-10-13 Toa Medical Electronics Co Ltd 画像ファイリングシステム
JP3735190B2 (ja) * 1997-10-28 2006-01-18 オリンパス株式会社 走査型サイトメータ
US6466695B1 (en) * 1999-08-04 2002-10-15 Eyematic Interfaces, Inc. Procedure for automatic analysis of images and image sequences based on two-dimensional shape primitives
US6203992B1 (en) * 1999-10-15 2001-03-20 Abbott Laboratories Nucleic acid primers and probes for detecting tumor cells
US6734401B2 (en) * 2000-06-28 2004-05-11 3M Innovative Properties Company Enhanced sample processing devices, systems and methods
US7087203B2 (en) * 2000-11-17 2006-08-08 Nagaoka & Co., Ltd. Methods and apparatus for blood typing with optical bio-disc

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102472703A (zh) * 2009-10-30 2012-05-23 西门子公司 体液分析系统以及用于分析体液的图像处理设备和方法
CN104508481A (zh) * 2012-06-01 2015-04-08 原子能和替代能源委员会 用于对诸如血液颗粒之类的包含在流体中的颗粒的移动速率进行表征的方法和系统
CN104508481B (zh) * 2012-06-01 2017-06-16 原子能和替代能源委员会 用于对诸如血液颗粒之类的包含在流体中的颗粒的移动速率进行表征的方法和系统
CN106332551A (zh) * 2014-03-31 2017-01-11 Jvc 建伍株式会社 分析装置和分析方法
CN106332551B (zh) * 2014-03-31 2019-10-15 Jvc 建伍株式会社 分析装置和分析方法
CN106644897A (zh) * 2016-10-14 2017-05-10 北京海岸鸿蒙标准物质技术有限责任公司 一种用于颗粒计数标准物质的计数装置
CN108169081B (zh) * 2017-12-14 2021-06-04 四川大学华西医院 血细胞分析的差值核查方法及其应用方法
CN108169081A (zh) * 2017-12-14 2018-06-15 四川大学华西医院 血细胞分析的差值核查模型及其应用方法
CN109085332A (zh) * 2018-09-29 2018-12-25 厦门大学 一种动态数据采集装置
CN111539354B (zh) * 2020-04-27 2020-12-15 易普森智慧健康科技(深圳)有限公司 一种液基细胞玻片扫描区域识别方法
CN111539354A (zh) * 2020-04-27 2020-08-14 易普森智慧健康科技(深圳)有限公司 一种液基细胞玻片扫描区域识别方法
CN111832389A (zh) * 2020-05-25 2020-10-27 中国人民解放军陆军军医大学第二附属医院 一种骨髓细胞形态学自动检测系统的计数及分析方法
CN111832389B (zh) * 2020-05-25 2022-12-09 中国人民解放军陆军军医大学第二附属医院 一种骨髓细胞形态学自动检测系统的计数及分析方法

Also Published As

Publication number Publication date
US20030096324A1 (en) 2003-05-22
JP2005502369A (ja) 2005-01-27
WO2003023571A3 (en) 2004-03-04
WO2003023571A2 (en) 2003-03-20
EP1425696A2 (en) 2004-06-09
AU2002333589A1 (en) 2003-03-24

Similar Documents

Publication Publication Date Title
CN1575475A (zh) 用于差分细胞计数的方法以及用于执行该方法的相关装置和软件
CN1592684A (zh) 细胞的定性和定量分析方法及相关的光学生物盘系统
CN1659439A (zh) 利用光学生物盘系统基于细胞核形态学的白血细胞类型的识别和定量
US7336812B2 (en) System for microvolume laser scanning cytometry
US6159740A (en) Method and apparatus for screening obscured or partially obscured cells
CN101044404A (zh) 用于不同类型细胞检测和定量的方法和装置以及生物光盘在该检测和定量中的应用
AU2009212193A1 (en) Method and system for analysis of flow cytometry data using support vector machines
WO2002003052A2 (en) Laser scanning cytology with digital image capture
CN1062212A (zh) 光学筛选显微细胞的方法和仪器
CN1759321A (zh) 包含等径和/或螺旋分析区域的光盘以及相关的光盘驱动系统与方法
US5260192A (en) Method and apparatus for screening cells or formed bodies with populations expressing selected characteristics utilizing at least one sensing parameter
EP1203339A1 (en) System for microvolume laser scanning cytometry
CN1636141A (zh) 用于生物驱动器的分段区域检测器及其相关方法
Yoon et al. Analytical performance of the digital morphology analyzer Sysmex DI-60 for body fluid cell differential counts
CN1050771A (zh) 用至少一个检测参数筛选具有特定特征值群体的细胞或生物体的方法和装置
CN1034242C (zh) 筛选细胞或成形体的方法和设备
EP0381669A1 (en) AUTOMATIC ANALYZER AND METHOD FOR SELECTING CELLS OF BODIES FORMED FOR THE ENUMERATION OF POPULATIONS EXPRESSING SELECTED CHARACTERISTICS.
Tzanoudaki et al. Basic Principles of Flow Cytometry
EP4462101A1 (en) Methods, system and non-transitory computer-readable storage medium for classifying analyte data into clusters
CN1829914A (zh) 用于分析细胞的光学生物盘和射流回路及其相关方法
Sulçe et al. Applicability of flow cytometry in identifying and staging lymphoma, leukemia and mast cell tumors in dogs: an overview
CN1636140A (zh) 包括生物和医学成像有关方法的光盘分析系统
CN1894586A (zh) 使用光学生物盘的血型测定的方法和装置
WO2005085842A2 (en) Method for the simultaneous detection of populations of several different biological entities using flow cytometry, device and computer program therefor
EP4713685A1 (en) Apparatus and methods for microscopic analysis of a biological sample

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
AD01 Patent right deemed abandoned

Effective date of abandoning: 20050202

C20 Patent right or utility model deemed to be abandoned or is abandoned