CN109477784A - 具有增强图像分辨率的荧光成像流式细胞术 - Google Patents

具有增强图像分辨率的荧光成像流式细胞术 Download PDF

Info

Publication number
CN109477784A
CN109477784A CN201780025417.5A CN201780025417A CN109477784A CN 109477784 A CN109477784 A CN 109477784A CN 201780025417 A CN201780025417 A CN 201780025417A CN 109477784 A CN109477784 A CN 109477784A
Authority
CN
China
Prior art keywords
sample
image
signal
frequency
radiation
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.)
Granted
Application number
CN201780025417.5A
Other languages
English (en)
Other versions
CN109477784B (zh
Inventor
埃里克·迪博尔德
基根·奥斯利
乔纳森·林
马修·巴尔
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.)
BD Bioscience
Original Assignee
BD Bioscience
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 BD Bioscience filed Critical BD Bioscience
Publication of CN109477784A publication Critical patent/CN109477784A/zh
Application granted granted Critical
Publication of CN109477784B publication Critical patent/CN109477784B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/50—Depth or shape recovery
    • G06T7/521—Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
    • 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
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429—Signal processing
    • G01N15/1433—Signal processing using image recognition
    • 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
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429—Signal processing
    • 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
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1468—Optical investigation techniques, e.g. flow cytometry with spatial resolution of the texture or inner structure of the particle
    • G01N15/147—Optical investigation techniques, e.g. flow cytometry with spatial resolution of the texture or inner structure of the particle the analysis being performed on a sample stream
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/645—Specially adapted constructive features of fluorimeters
    • G01N21/6456—Spatial resolved fluorescence measurements; Imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • 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/01—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials specially adapted for biological cells, e.g. blood cells
    • 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
    • G01N2015/1006—Investigating individual particles for cytology
    • 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
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434—Optical arrangements
    • G01N2015/144—Imaging characterised by its optical setup
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/06—Illumination; Optics
    • G01N2201/061—Sources
    • G01N2201/06113—Coherent sources; lasers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/10—Scanning
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing
    • G01N2201/129—Using chemometrical methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30024—Cell structures in vitro; Tissue sections in vitro

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Signal Processing (AREA)
  • Optics & Photonics (AREA)
  • Multimedia (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

在一个方面，公开了用于进行流式细胞术的***，其包括用于产生用于照射样品的激光辐射的激光器，至少一个用于检测响应于所述照射而从样品发出的辐射的至少一部分以便产生对应于所述检测到的辐射的时间信号的检测器，以及用于接收所述时间信号并基于正演模型对所述信号进行统计分析以重建所述样品的图像的分析模块。

Description

具有增强图像分辨率的荧光成像流式细胞术
相关申请
本申请要求于2016年5月12日提交的标题为“Fluorescence Imaging FlowCytometry With Enhanced Image Resolution”的临时专利申请第62/335,359号的优先权，其通过引用整体并入本文。本申请还涉及2017年3月17日提交的标题为“Cell SortingUsing A High Throughput Fluorescence Flow Cytometer”的实用新型专利申请第15/461,124号，所述实用新型专利申请通过引用整体并入本文。
背景技术
本发明总体上涉及用于样品的荧光分析的装置和方法，更具体地涉及用于基于荧光的流式细胞术的装置和方法。
荧光成像具有多种生物医学应用，例如，在获得关于生物样品的分子组成的信息。在生物医学流式细胞术中，收集并分析由外源和/或内源细胞荧光团发射的荧光辐射，以导出关于细胞的化学和/或物理性质的信息。
然而，在常规的基于荧光的流式细胞术中，获取高速流动的细胞或其它快速现象(诸如亚毫秒生化动力学)的无模糊图像可具有挑战性。特别地，许多荧光团的弱光学发射与在千赫兹范围内的成像帧速率下的样品的短曝光相结合使得难以获得无模糊图像。此外，许多常规***仅在一种成像模式下操作，因此不能为样品分析提供足够的灵活性。
因此，需要用于荧光分析的增强方法和***，特别是需要用于进行基于荧光的流式细胞术的增强方法和***。
发明内容
在一个方面，公开了用于进行流式细胞术的***，其包括用于产生用于照射样品(例如，细胞)的激光辐射的激光器，至少一个被布置来检测响应于所述照射而从样品发出的辐射的至少一部分并产生对应于所述检测到的辐射的时间信号的检测器，以及被连接到接收所述时间信号的检测器并基于诸如正演模型(forward model)的模型执行所述信号的统计分析以重建所述样品的图像的分析模块。在一些实施方案中，正演模型是线性模型，而在其它实施例中，正演模型是非线性模型。
在一些实施方案中，从样品发出的辐射可以是荧光辐射。从而时间信号为荧光信号，分析模块对所述荧光信号进行操作以形成样品的荧光图像。
在一些实施方案中，从样品发出的辐射可以是由样品散射的激光辐射的一部分，并且分析模块可对对应于散射的激光辐射的时间信号进行操作以生成样品的暗场图像。
在一些实施方案中，从样品发出的辐射可以是透射通过样品的照射激光辐射的一部分，并且分析模块可对对应于透射的激光辐射的时间信号进行操作以生成样品的明场图像。
在一些实施方案中，统计分析基于生成所述模型的一个或多个参数的最大似然估值。举例来说，统计分析可采用最小二乘法方法通过最小化对应于检测到的时间信号与从正演模型推断的相应时间信号之间的差的残差平方和来重建样品的图像。
适用于本发明实践的统计分析的一些实例包括但不限于最小二乘优化、梯度下降优化、粒子群优化、遗传算法优化和贝叶斯频谱估计。
在一些实施方案中，分析模块形成样品的初始图像，并将所述图像用作所述统计分析的初始(种子)图像。举例来说，分析模块可采用傅里叶变换和锁定检测中的任一种来形成初始图像。
在一些实施方案中，分析模块将检测到的时间信号的模型用作多个时间区段。举例来说，可将每个时间区段建模为具有正弦和余弦项。在此类实施方案中，统计分析可采用最小二乘回归法，以便通过最小化对应于建模的时间区段与对应的测量的时间区段之间的差的残差平方和来获得建模的时间信号的参数的值，例如正弦项和余弦项的幅度和相位。
在一些实施方案中，统计分析采用梯度下降优化方法。在此类实施方案中，梯度下降优化方法计算指示基于正演模型的预期时间信号与测量的时间信号之间的距离的误差梯度，并且将先前样品图像的估计沿着误差梯度向下步进以获得改善的图象。可以迭代地重复该过程，直至获得最佳拟合图像。梯度下降优化方法可利用基于测量的时间信号计算的图像的初始估计来开始迭代过程。举例来说，分析模块可通过对测量的时间信号应用傅里叶变换(例如，FFT)来计算初始图像的估计。
在一些实施方案中，分析模块采用关于与贝叶斯频谱估计结合的测量的时间信号的先验信息以重建样品的图像。举例来说，在一些实施方案中，此类先验信息可以指示时间信号由许多具有未知频率、幅度和相位的正弦曲波组成。另外，在一些情况下，可以已知时间信号包括噪声分量，例如具有泊松分布的噪声分量。在这种情况下，贝叶斯频谱估计可提供正弦波的未知频率、幅度和相位的估计以及与建模的噪声相关的参数。
在一些实施方案中，激光辐射包括彼此偏移一个射频的至少两个光学频率。举例来说，射频偏移可以在约10MHz至约250MHz的范围内。在一些此类实施方案中，激光辐射可呈多个角度分离的激光束的形式，每个激光束相对于另一个具有射频偏移。在一些此类实施方案中，时间信号(例如，荧光信号)可以是表现出与所述至少两个光学频率之间的差相对应的至少一个拍频的时频波形。
在一些实施方案中，研究中的样品可包括任何细胞(例如，哺乳动物细胞)、微囊泡、细胞碎片、脂质体、珠粒或小生物体。
如下面更详细讨论的，分析模块可以在硬件、固件和/或软件中实现。在一些实施方案中，此类分析模块，除了其它组件以外，还可包括可以经由一个或多个通信总线与一个或多个存储器模块通信的处理器。处理器可以执行实现正演模型的程序以及根据本教导的统计分析，用于在流式细胞仪中重建处于询问中的样品的图像。
在相关方面，公开了一种用于执行流式细胞术的方法，其包括用激光辐射照射样品(例如，细胞)，部署检测器以检测响应于所述照射而从样品发射的辐射的至少一部分并且以产生对应于所述检测到的辐射的时间信号，以及使用计算机处理器基于模型(例如，正演模型)执行时间信号的统计分析，以便重建所述样品的图像。举例来说，重建图像可以是荧光、暗场和明场图像中的任一种。在一些实施方案中，激光辐射可包括彼此偏移一个射频的至少两个光学频率，其中选择光学频率以引发来自样品的荧光辐射。在此类实施方案中，时间信号可以包括与射频移位的光频率相关联的一个或多个拍频。可以处理此类时间信号以生成种子荧光图像，在该种子荧光图像上可进行统计分析以重建样品的增强图像。
举例来说，统计分析可采用最小二乘法方法通过最小化对应于检测到的时间信号与从正演模型推断的相应时间信号之间的差的残差平方和来重建样品的图像。时间信号可以是例如荧光信号、散射信号或传输信号中的任一种。
在一些实施方案中，时间信号可以被建模为多个时间区段。在一些此类实施方案中，每个时间区段可以被建模为具有一个正弦和一个余弦项。在此类实施方案中，统计分析可以采用最小二乘回归方法，通过最小化对应于建模的与相应的测量的时间区段之间的差的残差平方和来获得建模的时间信号的参数的值。
在一些实施方案中，统计分析采用梯度下降优化方法来重建图像。梯度下降优化方法可以计算指示基于正演模型计算的预期时间信号与测量的时间信号之间的距离的误差梯度，并且可以通过将先前图像沿着误差梯度向下步进来迭代地计算更新的图像。迭代过程可以从基于测量的时间信号计算的样品的初始估计图像开始，例如，通过对测量的时间信号应用傅里叶变换。
在上述方法的一些实施方案中，统计分析可以采用关于与贝叶斯频谱估计结合的测量的时间信号的先验信息以重建样品的图像。举例来说，在一些实施方案中，此类先验信息可以指示时间信号由许多具有未知频率、幅度和相位的正弦曲波组成。贝叶斯频谱估计可提供频率、幅度和相位的估计。另外，在某些情况下，先验信息可以指示时间信号可包括具有泊松频率分布的噪声分量。贝叶斯估计可以确定噪声分量的参数。
在上述方法的其它实施方案中，统计分析可以采用粒子群或遗传优化方法中的任一种。
上述方法可用于获得各种不同样品诸如含有细胞(例如，哺乳动物细胞)、微囊泡、细胞碎片、脂质体、珠粒或小生物的样品的图像。
通过参考以下结合相关附图的详细描述，可以获得对本发明各个方面的进一步理解，所述附图在下面简要描述。
附图说明
图1示意性地描绘了根据本发明的实施方案的***，
图2A是垂直于光束传播方向的平面中的高斯光束的示意性示例性轮廓，
图2B是通过使图2A中所示的高斯光束通过顶帽式光束整形器并聚焦光束整形器的输出光束获得的示意性顶帽式光束轮廓。
图3示意性地描绘了示例性顶帽式光束整形器的组件，
图4示意性地描绘了多个RF梳状光束的横截面光束轮廓，
图5示意性地描绘了图4中描绘的RF梳状光束和具有顶帽式光束轮廓的LO光束的叠加，
图6示意性地描绘了图5中所示的照射分析中的样品的组合光束，
图7示意性地描绘了假想荧光团的示例性能级，
图8示意性地描绘了对应于图7的假想荧光团的吸收曲线，
图9A示意性地描绘了根据本教导的实施方案的检测***，其包括用于传输荧光辐射的光纤
图9B示意性地描绘了根据本教导的实施方案的另一检测***，其中荧光辐射通过自由空间传播以到达多个光电探测器，
图9C示意性地描绘了用于本教导的一些实施方案的明场和暗场图像生成臂，
图9D示意性地描绘了用于本教导的一些实施方案的检测***，其包括用于生成明场图像的检测臂以及将检测从样品散射的激发辐射以及由样品发射的荧光辐射的能力集成的检测臂，
图10示意性地描绘了由根据本发明的***的实施方案中的光电探测器产生的荧光信号可通过放大器放大，并且可通过分析模块分析放大的信号以构建分析中的样品的荧光图像，
图11A和11B描绘了用于分析通过用由多个RF梳状形光束和顶帽式轮廓LO光束组成的组合光束照射样品而获得的荧光信号的根据本发明实施方案的方法中的各个步骤，
图12示意性地描绘了根据本发明实施方案的分析模块的示例性硬件实现的所选组件，
图13A和13B描绘了用于分析通过用由多个RF梳状形光束和顶帽式轮廓LO光束组成的组合光束照射样品而获得的荧光信号的根据本发明实施方案的另一方法中的各个步骤，
图14A和14B描绘了用于分析通过用由多个RF梳状形光束和顶帽式轮廓LO光束组成的组合光束照射样品而获得的荧光信号的根据本发明实施方案的另一方法中的各个步骤，
图15A示意性地描绘了在单个激发频率下由顶帽式轮廓光束照射样品，
图15B是允许荧光寿命测量和荧光寿命成像的根据本教导的实施方案的***的示意图，
图16A是通过根据本教导的实施方案，使用细胞计数***获得的用8个离散水平的荧光染料染色的聚苯乙烯珠的暗场强度与明场强度的散点图，
图16B显示由多个聚苯乙烯珠粒发射的红色荧光(PE)对比绿色荧光(FITC)的散点图，其中通过使用图16A中所示的图的矩形部分作为门控来获得数据，
图16C和16D是对应于图16B中所示的数据的直方图，
图17A显示根据本教导使用细胞计数***获得的用抗CD45-FITC和碘化丙锭染色的含有固定的外周血白细胞的样品的明场、暗场和荧光图像，其中样品还含有一部分的用钙黄绿素-AM染色的活的HeLa细胞，
图17B是散点图，其中群体A代表白细胞，群体B代表图17A为其提供图像的样品中的HeLa细胞，
图18A显示根据本教导的一个实施方案，使用细胞计数***获得的用抗CD45-FITC和碘化丙啶染色的含有固定的外周血白细胞的样品的明场、暗场和荧光图像，其中样品掺有一小部分用抗EpCAM-FITC和碘化丙啶染色的固定的MCF-7细胞，
图18B是散点图，其中群体A代表白细胞，群体B代表图18A为其提供图像的样品中的MCF-7细胞，
图19A是描绘根据本教导的方法的实施方案中的各个步骤的流程图，所述方法用于进行流式细胞术，其中统计分析用于重建处于询问中的样品的图像，
图19B是描绘根据本教导的方法的实施方案中的各个步骤的流程图，所述方法用于进行流式细胞术，其中细胞计数信号被建模为多个时间区段，并且最小二乘回归分析用于获得模型参数的估计值，
图19C是描绘根据本教导的方法的实施方案中的各个步骤的流程图，所述方法用于从获自响应于其照射的样品的细胞计数信号获得样品的图像，
图19D是描绘根据本教导的方法的实施例中的各个步骤的流程图，所述方法用于进行流式细胞术，其中梯度下降优化方法用于迭代地改进样品的图像估计以便获得样品的最佳拟合图像，
图20示意性地描绘了根据本教导的实施方案的细胞计数***的各种组件的建模，
图21A描绘了细长的眼虫(Euglena)细胞的消光和散射的测量，
图21B显示通过提取图21A中描绘的信号的傅里叶分量生成的眼虫细胞的图像。
图21C描绘了根据本教导的实施方案使用***的激发以及电子和相位响应的模型生成的眼虫细胞的改进图像，
图22A描绘了图21B和图21C中描绘的眼虫细胞的脂质含量的FIRE(射频标记发射(radio-frequency-tagged-emission))测量，
图22B是通过提取图22A中描绘的信号的傅里叶系数而获得的图像，
图22C是通过合并关于收集电子器件的频率响应、驱动***的基础FIRE信号以及***的相位响应的信息而获得的图22B中所示图像的改进版本，以及
图22D是通过采用基础信号的模型以及光子散粒噪声的影响而获得的图22B所示图像的改进形式。
具体实施方式
本教导总体上涉及用于进行流式细胞术，例如样品的荧光分析的方法和***。在一些实施方案中，统计技术可用于从时间数据(例如，以时频波形的形式存在)重建荧光、明场和暗场图像中的任一种图像，所述时间数据获自响应于询问辐射(例如，激光辐射)对样品的照射的样品。可以使用各种不同的***获得时间数据。在下面公开的一些实施方案中，可通过用具有多个射频移位子光束的光束照射样品来生成时间数据。下文中描述这些实施方案的实例，但是本发明的教导不限于这些实施方案。
在许多实施方案中，此类统计技术的使用可有利地改善重建图像的空间分辨率、强度分辨率和/或信噪比。在一些实施方案中，可使用例如信号的FFT分析来重建初始(种子)荧光、明场或暗场图像。随后，可以将下面讨论的至少一种统计技术应用于图像以生成增强图像。在一些实施方案中，统计技术可用于例如，以下面讨论的方式将模型的参数直接与由一个或多个检测器检测的荧光、明场和/或暗场信号拟合。
除非另有说明，否则以下用于描述本教导的各种术语具有其在本领域中的普通含义。例如，术语“荧光团”在本文中与其在本领域的常规含义一致，是指可响应激发辐射的照射而发射辐射的荧光化合物。
术语“细胞计数”和“流式细胞术”也与它们在本领域的常规含义一致。特别地，术语“细胞计数”可以指用于鉴定和/或分选或以其它方式分析细胞的技术。术语“流式细胞术”可以指细胞计数技术，其中存在于流体流中的细胞可以被鉴定和/或分类，或以其它方式被分析，例如，通过用荧光标记物标记它们并通过辐射激发来检测荧光标记物。如本文中所用，术语“约”和“基本上”表示相对于包括数值的性质的最大10％或5％的变化。
术语“正演模型”是指用于在数学上预测***的输出或用于给定输入(例如，一个或多个输入参数)的过程的***或过程的数学模型。举例来说，在以下实施方案中的一些实施方案中，正演模型可以是流式细胞仪中处于询问中的样品的像素化图像的数学模型，其可以用于预测从在流式细胞仪中响应于其照射的样品获得的时间信号。然后可以使用统计分析来细化模型的参数，例如，通过最小化预测的时间信号和测量的时间信号(例如，通过检测响应于询问辐射(诸如RF调制的照射辐射)照射而从样品发出的辐射而测量的信号)之间的差异。
图1示意性地描绘了根据本教导的实施方案的用于进行细胞计数的***10，其可以以三种操作模式操作。如下面更详细讨论的，在一种操作模式中，研究中的样品可以用多个激发频率同时照射，每个激发频率可以例如通过移动激光束的中心频率来获得。更具体地，多个样品位置可以由激光束同时照射，该激光束通过将参考激光束(本文中也称为本地振荡器光束)与多个射频移位的激光束混合而产生，使得每个样品位置由参考光束和射频移位光束之一照射，以在该位置激发目标荧光团(如果存在的话)。在一些实施方案中，参考光束本身可以通过激光束的射频移位产生。因此，样品的每个空间位置可以用不同的拍频“标记”，该拍频对应于参考光束的频率与射频移位光束之一的频率之间的差。换句话说，荧光团发射的荧光辐射将在空间上编码拍频。可检测荧光发射并且可分析其频率分量以构建样品的荧光图像。
在另一种操作模式中，可以通过多个激发频率的激光束在一段时间间隔内连续照射样品。在一些此类实施立案中，可以通过将时变驱动信号施加到接收激光束的声光偏转器(AOD)来获得激发频率。在许多实施方案中，激光束具有在数百太赫兹(THz)范围内的频率，例如，在约300THz至约1000THz的范围内。施加到AOD的驱动信号通常在射频范围内，例如，在大约10MHz至约250MHz的范围内。激光束通过AOD而产生多个衍射光束，每个衍射光束对应于不同的衍射级。然而零级衍射光束相对于输入激光束的频率没有表现出频移，高阶衍射光束表现出相对于输入激光束的频率的频移(对应于驱动信号的频率或其倍数)。在一些实施方案中，具有与偏移了驱动信号的输入激光束的频率相对应的频率的一阶衍射光束被用作激发电影票荧光团(如果存在于分析中的样品中的疾)的激发光束。随着驱动信号随时间变化，一阶衍射光束的频率和角度偏移也变化，从而允许在不同位置处以不同激发频率照射样品。可收集并分析来自每个照射位置的荧光发射(如果有的话)以构建样品的荧光图像。
在另一种操作模式中，可操作***10以通过单个激发频率同时照射样品的多个位置，该激发频率可以例如通过使激光束的中心频率偏移一个射频来产生。例如，可通过激光束以单个激发频率照射样品的水平范围。检测到的荧光辐射可用于分析样品(例如，细胞/颗粒)的荧光含量。
因此，除了下面讨论的其它有利方面外，***10的一个有利方面在于其提供了在不同模式下获得荧光发射数据的显著灵活性，而无需使用不同的仪器或在不同操作模式之间切换时对***进行任何机械修改。
在某些实施方案中，***包括一个或多个光源。在一些情况下，光源是窄带光源，包括但不限于窄波长LED、激光器或耦合到一个或多个光学带通滤波器、衍射光栅、单色器或其任何组合的宽带光源(其组合在一起产生一束窄频带照明光)。在某些情况下，光源是单波长激光器，诸如单波长二极管激光器(例如，488nm激光器)。在一些实施方案中，主题***包括单个光源(例如，激光器)。在其它实施方案中，主题***包括两个或更多个不同的光源，诸如3个或更多个不同的光源，诸如4个或更多个不同的光源，并且包括5个或更多个不同的光源。例如，***可包括输出第一波长的第一光源(例如，激光器)和输出第二波长的第二光源。在其它实施方案中，***包括输出第一波长的第一光源、输出第二波长的第二光源和输出第三波长的第三光源。
每个光源的波长可以为300nm至1000nm，诸如350nm至950nm，例如400nm至900nm，并且包括450nm至850nm。在某些实施方案中，光源具有对应于一种或多种荧光团的最大吸收的波长(如下所述)。例如，光源可以输出波长在280-310nm、305-325nm、320-350nm、370-425nm、400-450nm、440-500nm、475-550nm、525-625nm、625-675nm和650-750nm中的一个或多个的范围内的光。在某些实施方案中，每个光源输出具有选自348nm、355nm、405nm、407nm、445nm、488nm、640nm和652nm的波长的光。
***10包括产生激光束14的激光辐射源12。举例来说，激光束的频率可以在约1000THz至约300THz的范围内，对应于在约300nm至约1000nm范围内的真空波长。激光束的光束直径(例如，当采用高斯激光束时的束腰)可以在例如约0.1mm至约10mm的范围内。在不失一般性的情况下，在该实施方案中，激光器12发射波长为488nm的光束，光束直径约为1mm。
可以基于***所针对的一个或多个特定应用来选择激光束的频率。具体地，如下面更详细讨论的，激光频率可以适合于激发目标荧光团的电子跃迁(例如，通过吸收辐射)，从而使荧光团以较低的频率发射荧光辐射。可以使用各种激光源。此类激光源的一些实例包括但不限于Sapphire488-SF(由Coherent,Inc.of Santa Clara,CA U.S.A.销售)、Genesis MX-488-1000-STM(Coherent,Inc.)、OBIS 405-LX(Coherent,Inc.)、Stadus405-250(由Vortran Laser Technology,Inc.of Sacramento,CA U.SA.销售)和NewportCorporation of Irvine,CA U.S.A.的LQC-660-110。在不失一般性的情况下，在本实施方案中，假设激光束在垂直于其传播方向的平面中具有高斯强度分布。
反射镜16接收激光辐射束14并通过反射将激光束引导至声光偏转器(AOD)18。在该实施方案中，AOD 18被安装在可调节的柱式支架底座(A)上，其允许AOD围绕垂直于光束14的传播方向的轴旋转。在控制器21的控制下操作的直接数字合成器(DDS)20可将一个或多个驱动信号施加到AOD 18。举例来说，在一些实施方案中，这些驱动信号可以跨越约50MHz至约250MHz的频率范围。在一些实施方案中，驱动信号可以通过约0.1MHz至约4MHz范围内的频率彼此分离。在该实施方案中，电子功率放大器21'放大由DDS 20产生的射频信号以施加到AOD 18。
在用多个激发频率同时照射样品的操作模式中，RF梳状发生器20同时将多个RF驱动信号施加到AOD 18。举例来说，同时施加的RF驱动信号的数量可以在约20至约200的范围内。激光束与驱动信号的相互作用导致产生多个角度分离的激光束，每个激光束具有对应于驱动信号之一的相对于激光器12产生的激光束的频率的频移。不受任何特定理论的限制，在AOD中，压电换能器可以在晶体(例如，石英晶体)中产生射频声子，并且这种射频声子对激光束的光学光子的散射可以导致产生频移激光束。这些频移光束22中的一个在本文中被称为“本地振荡器”(LO)光束，而频移光束24的其余部分在本文中被称为“RF梳状光束”。频移光束的角度间隔可以例如在约1毫弧度至约100毫弧度的范围内。
LO和RF梳状光束穿过透镜26，在该实施方案中，透镜26是焦距为约50mm的正透镜。在通过透镜26之后，LO激光束被反射镜28截取，反射镜28使LO光束在不同方向上重定向(在该实施方案中，在与LO光束的原始传播方向基本正交的方向上)。将反射镜28相对于RF梳状光束放置，使得这些光束错过反射镜28并传播到透镜30(在该实施方案中，透镜30的焦距为200mm)。以这种方式，LO光束和RF梳状光束沿不同的传播方向被引导。以上面公开的方式使用拾取镜28允许利用单个AOD产生LO光束和RF梳状光束，并以下面讨论的方式将它们组合以产生用于照射样品的激发光束。使用单个AOD而不是多个AOD(例如，两个AOD，一个用于产生LO光束，另一个用于产生RF梳光束)，简化了***的设计，并且还允许以如在下文中更详细讨论的多个不同的操作模式高效使用***。
在一些实施方案中，在与RF光束梳重组之前修改LO光束的光束轮廓。例如，可以在空间维度、光束形状、强度、光束的空间分布或其任何组合方面调整(增大或减小)LO光束的光束轮廓。在某些实施方案中，修改LO光束的光束轮廓的空间尺寸。例如，可调节光束轮廓以在一个或多个维度上(诸如沿着与流动流的纵向轴正交的长轴)延长光束轮廓。在根据这些实施方案的一个实例中，光束轮廓的空间维度(例如，在一个或多个维度中)可以增加1％或更多，诸如2％或更多，诸如3％或更多，诸如5％或更多，诸如10％或更多，诸如25％或更多，诸如50％或更多，诸如75％或更多，诸如25％或更多，诸如50％或更多，诸如75％或更多，诸如90％或更多，诸如1.5倍或更多倍，诸如2倍或更多倍，诸如3倍或更多倍，包括5倍或更多倍。在根据这些实施方案的一个实例中，光束轮廓的空间维度(例如，在一个或多个维度中)可以增加1％或更多，诸如2％或更多，诸如3％或更多，诸如5％或更多，诸如10％或更多，诸如25％或更多，诸如50％或更多，诸如75％或更多，诸如90％或更多，诸如1.5倍或更多倍，诸如2倍或更多倍，诸如3倍或更多倍，包括5倍或更多倍。
在其它实施方案中，修改LO光束的光束形状。例如，可以修改光束形状以在一个或多个维度上拉长光束轮廓。在某些情况下，在垂直于LO光束传播的平面中拉长LO光束的光束形状。在某些实施方案中，LO光束轮廓的形状从圆形光束轮廓变为椭圆形光束轮廓，该椭圆形光束轮廓在垂直于流动流的纵向轴线的轴线上拉长。在某些实施方案中，LO光束轮廓的形状从圆形光束轮廓变为椭圆形光束轮廓，该椭圆形光束轮廓在垂直于流动流的纵向轴线的轴线上拉长。在其它实施方案中，修改LO光束的强度。例如，可使LO光束的强度增加诸如1％或更多，诸如2％或更多，诸如3％或更多，诸如5％或更多，诸如10％或更多，诸如25％或更多，诸如50％或更多，诸如75％或更多，诸如90％或更多，诸如1.5倍或更多，诸如2倍或更多倍，诸如3倍或更多倍，包括5倍或更多倍。在某些实施例中，修改LO光束的强度以匹配RF梳状光束的强度。例如，LO光束的强度可以与RF梳头光束的强度相差10％或更少，诸如9％或更小，诸如8％或更小，诸如3％或更小，诸如2％或更少，诸如1％或更少，诸如0.01％或更少，并且包括LO光束的强度与RF梳状光束相差0.001％或更少的情况。在某些情况下，LO光束和RF梳状光束的强度是相同的。
在其它实施方案中，还可以修改光束轮廓的空间分布。例如，可以修改LO光束，使得LO光束的强度在一个或多个维度上不再是高斯分布。例如，可以修改LO光束以沿着与流动流的纵向轴线平行的第一轴线具有高斯分布，并且沿着与流动流的纵向轴线正交的第二轴线具有非高斯分布。
可采用任何光束整形方案来修改LO光束的光束轮廓，包括但不限于折射和衍射光束整形方案。在一些实施方案中，LO光束由顶帽式光束整形器修改。
在该实施方案中，LO光束传播到另一个正透镜32(在该实施方案中，其具有约200mm的焦距)。透镜26和透镜32的组合放大并准直LO光束，以便适当地填充顶帽式光束整形器34的后孔。更具体地，LO光束22穿过透镜32并被反射镜33和35反射到顶帽式式光束面形器34。
顶帽式光束整形器34对高斯LO光束的相前(phase front)进行整形，以形成顶帽式强度分布。更具体地，离开顶帽式光束整形器的LO激光束22'被分束器44反射，并被透镜46(在该实施方案中具有100mm的焦距)聚焦到中间像平面48上。中间像平面48上的激光束在垂直于光束传播方向的平面中沿水平方向具有顶帽式强度分布。类似于AOD 18，在该实施方案中，将分束器44安装在可调节的柱式支架底座(B)上。在该实施方案中，顶帽式光束整形器产生顶帽式光束轮廓，其中辐射的偏振沿着光束的顶帽方向(在该实施方案中沿着水平方向)基本上是均匀的。
作为说明，图2A示意性地描绘了LO激光束进入顶帽式光束整形器时的高斯强度分布。如图2B中示意性所示，在中间像平面48上，LO激光束呈现出在水平方向上(在该说明图中垂直于页面的方向上)拉伸的光束轮廓，并且沿着延伸通过轮廓的每条水平线(例如，水平线A)基本恒定，但按照高斯轮廓垂直变化。
可以使用各种顶帽式光束整形器。举例来说，具有非球面表面或衍射光学元件的折射光学元件可用于产生具有适当空间相前的光束，其在通过透镜聚焦后将在透镜的焦平面处生成顶帽式轮廓图案。对于此类顶帽式光束整形器存在多种形状因子，并且该方法的各种实施方式可用于在本教导的各种实施方案中在样品处产生适当的LO光束形状。例如，标题为“Refractive optical system that converts a laser beam to a collimatedflat-top beam”的美国专利第6,295,168号和标题为“Rectangular flat-top beamshaper”的美国专利第7,400,457号(这两个专利均通过引用整体并入本文)公开了可在根据本教导的一些实施例的***中用作平顶式光束整形器的光束整形***。作为说明，图3是美国专利第7,400,457号的图1的再现(具有不同的参考标号)，所述图1示意性地描绘了用于提供正方形或矩形光束的光束整形***300，其包括两个正交设置的非圆柱透镜302和304。第一非圆柱透镜302用于沿X轴对入射光束A进行成形，第二非圆柱透镜304用于沿Y轴对入射光束A进行成形。两个交叉的非圆柱透镜适于提供具有沿X轴的平顶轮廓的所得矩形激光束B。非圆柱透镜302的输入表面302a是凸面的非圆柱形表面，其具有可变曲率半径，该曲率半径在表面的中心处较小并且朝向透镜的两个X-末端平稳地增加。第二非圆柱透镜304类似于第一非圆柱透镜，但相对于透镜302正交设置，以便沿Y轴对光束成形。可独立地选择透镜302和304的输入表面302a/304a和输出表面302b/304b的轮廓作为入射光束A的X和Y轮廓以及所得矩形光束B的所需强度分布(参见例如，专利的第5和第6列)的函数。
可使用的商购可得的顶帽式光束整形器的实例包括例如由Osela,Inc.ofLachine,Canada销售的DTH-1D-0.46deg-4mm。
如下面更详细讨论的，使用光束整形器沿水平方向拉伸LO光束提供了许多有利方面。例如，其可确保LO光束和RF梳状光束的组合以基本相似的照明强度照射多个样品位置，以便在整个样品位置上匹配LO和RF梳状光束的强度，从而产生具有高调制深度的荧光信号的强度幅度调制。在此类强度匹配不存在的情况下，成像***可具有小视图并且可以不利用驱动AOD的所有频率(像素)。由于荧光信号的调制深度在***重建样品的荧光图像的能力中起重要作用，因此在所有像素处的激发拍频的均匀地高的调制深度对于***的操作特别有利。另外，可通过控制直接数字合成器的输出(例如，通过使用控制器21)来调整施加到AOD以产生RF梳状光束的电子信号的幅度，以便均衡RF梳状光束，使得它们的强度在RF梳状光束和LO光束重叠的所有空间位置上等于LO光束的强度。该特征提供的有利方面在于其确保了荧光辐射的强度幅度调制的高调制深度。
再次参考图1，RF梳状光束24通过透镜26和30的组合成像到中间像平面38上。更具体地，RF梳状光束24穿过透镜26并且错过反射镜28以到达透镜30，所述透镜30将RF梳状光束经由反射镜40和42引导到中间像平面38。
图4示意性地描绘了中间图像平面38中的示例性数量的RF梳状束的分布(不失一般性，出于说明目的，RF梳状光束的数量被选择为6(标记为RF1，...，RF6)，虽然也可以使用其它数字)。如图4所示，在中间像平面38中，RF梳状束24沿水平方向在空间上彼此分离。在其它实施方案中，两个或更多个RF梳状束24可以部分重叠。因此，透镜26和30的组合将角度分离的RF梳状光束变换成跨越水平范围的一组空间分离的光束。
再次参考图1，如上所述，分束器44接收离开顶帽式光束整形器34的激光束22'并将该光束反射到透镜46，所述透镜46又将光束聚焦在其中LO光束展现出顶帽式光束轮廓的中间像平面48上。分束器还从中间像平面38接收RF梳状光束24，并允许RF梳状光束通过其中。透镜46将RF梳状光束24聚焦到中间像平面48上，以与具有顶帽式光束轮廓的LO光束组合，以产生组合光束49。
作为说明，图5示意性地描绘了组合光速49在垂直于其传播轴的平面中的一个示例性轮廓。组合光束的强度分布被生成为顶帽式LO光束的强度分布(由正方形示意性地显示)与RF梳状光束24的强度分布(每个由圆圈之一示意性地显示)的叠加。如下面更详细讨论的，LO光束和RF梳状光束的这种叠加沿水平范围提供多个拍频，每个拍频对应于沿该水平范围的一个空间位置。在照射样品的水平范围时，从样品的位置发射的荧光辐射经由幅度调制编码与照射该位置的辐射相关联的拍频。
再次参考图1，正透镜50(在该实施方案中为200mm透镜)和物镜52，在本实施方案中安装在可调节的柱式支架底座C上，形成望远镜，用于将中间平面48处的图像中继到流过流动池54的样品上。在该实施方案中，反射镜56将组合光束49反射到透镜50，并且分色镜58在组合光束通过透镜50后将其朝向物镜52反射。
如图6中示意性所示，组合光束49同时照射样品62流动通过流动池54的多个空间位置60。因此，每个位置60由RF梳形光束之一与帽式LO激光束之一的重叠照射。在这些空间位置，辐射将激发样品中的目标荧光团(如果存在的话)。更具体地，在该实施方案中，LO光束和RF梳状光束同时激发荧光团(例如，通过在多个样品位置60处使其电子跃迁到激发电子态)。
在一些实施方案中，样品可包括流动的流体，其中夹带有多个细胞。在一些情况下，细胞可以用一种或多种荧光标记物(荧光团)标记。荧光标记物的一些实例包括但不限于荧光蛋白(例如，GFP、YFP、RFP)、用荧光团标记的抗体(例如，异硫氰酸荧光素(FITC)、藻红蛋白(PE)、别藻蓝蛋白(APC))、核酸染色剂(例如，4',6-二脒基-2-苯基吲哚(DAPI)、SYTO16、碘化丙啶(PI))、细胞膜染色剂(例如，FMI-43)和细胞功能染料(例如，Fluo-4、Indo-1)。在其它情况下，细胞中存在的内源荧光团可用于引发细胞的荧光辐射。如下面更详细讨论的，此类种外源或内源荧光团响应于照射辐射而经历电子激发并发射荧光辐射(通常以比激发频率低的频率)，收集并分析所述荧光辐射。
作为说明而不限于任何特定理论，图7显示对应于电子基态A的假设能级以及荧光团的两个电子激发电子态B和C。荧光团可以通过吸收辐射能从其电子基态(A)激发到电子激发态(B)。然后荧光团可以弛豫到低激发态B，例如，通过由荧光团的振动模式介导的无辐射跃迁。荧光团可通过光学跃迁进一步从低电子态C弛豫到基态，从而以小于激发频率的频率发射荧光辐射。应当理解，提供该假想实例仅用于说明目的，而不是指示可籍此发射荧光辐射的唯一机制。
在许多情况下，荧光团可在一定频率范围内吸收电磁辐射，以从基态激发到激发电子态。作为说明，图8显示了结合图7讨论的假想荧光团的吸收曲线。在根据本教导的实施方案的一个实现中，可选择LO频率以与对应于目标荧光团的峰值吸收的频率一致。射频移位光束可具有通过它们各自的拍频与峰值吸收分开的频率。通常，这些频率分离与荧光团的吸收带宽相比较小，以避免激发频率的任何降低。作为示例并且仅作为说明，虚线A和B示意性地描绘了LO光束和RF梳状光束之一的频率(为了便于描述，附图未按比例绘制)。通过LO激光束和所描绘的RF梳状光束之一同时照射样品的空间位置导致荧光辐射在对应于LO与RF梳状光束频率之间的差的拍频处表现出幅度调制。
再次通过说明而不限于任何特定理论，通过其LO光束和RF梳状光束之一的同时照射施加到荧光团的电场可在数学上定义如下：
其中，
Ecom表示组合光束的电场，
ERF表示与一个RF梳状光束相关的电场幅度，
ELO表示与LO光束相关的电场的幅度，
表示激光器12产生的激光束的频率，
表示与RF梳状光束相关的频移，以及
表示与LO光束相关的频移。
响应LO和RF梳状光束的电场叠加而发射的荧光辐射的强度将表现出在对应于(ωRF-ωLO)的拍频处的调制。因此，从被LO光束和RF梳状光束之一的叠加照射的样品的每个空间位置发出的荧光辐射表现出在对应于与LO光束相关的射频偏移与与照射该空间位置的RF梳状光束的射频偏移之间的差的拍频处的调制。
由于荧光发射过程需要有限量的时间(对于普通有机荧光团通常为1-10纳秒)，如果激发拍频过高，则发射的荧光将不会表现出高调制深度。因此，在许多实施方案中，激励拍频被选择为远小于1/τf，其中τf是荧光团的特征荧光寿命。在某些情况下，激发拍频可比1/τf小1％或更多，诸如2％或更多，诸如3％或更多，诸如5％或更多，诸如10％或更多，诸如25％或更多，诸如50％或更多，诸如75％或更多，诸如90％或更多，诸如1.5倍或更多，诸如2倍或更多，诸如3倍或更多，包括5倍或更多。例如，激励拍频可比1/τf小0.01MHz或更多，诸如0.05MHz或更多，诸如0.1MHz或更多，诸如0.5MHz或更多，诸如1MHz或更多，诸如5MHz或更多，诸如10MHz或更多，诸如25MHz或更多，诸如50MHz或更多，诸如100MHz或更多，诸如250MHz或更多，诸如500MHz或更多，包括750MHz或更多。在一些实施方案中，光电探测器被配置来从被照射的样品中检测光(例如，诸如荧光的发光)。在一些实施方案中，光电探测器可包括一个或多个检测器诸如2个或更多个检测器，诸如3个或更多个检测器，诸如4个或更多个检测器，诸如5个或更多个检测器，诸如6个或更多个检测器，诸如7个或更多个检测器并包括8个或更多个检测器。可采用任何光检测方案，包括但不限于有源像素传感器(APS)、象限光电二极管、图像传感器、电荷耦合器件(CCD)、强化电荷耦合器件(ICCD)、发光二极管、光子计数器、辐射热测量计、热电探测器、光敏电阻器、光伏电池、光电二极管、光电倍增管、光电晶体管、量子点光电导体或光电二极管及其组合，以及其它光电探测器。在一些实施方案中，目标光电探测器被配置来检测在350nm至1200nm，诸如450nm至1150nm，诸如500nm至1100nm，诸如550nm至1050nm，诸如从500nm至1000nm并且包括400nm至800nm的范围内的光。在某些实施方案中，光电探测器被配置来检测发光的最大发射处，诸如在395nm、421nm、445nm、448nm、452nm、478nm、480nm、485nm、491nm、496nm、500nm、510nm、515nm、519nm、520nm、563nm、570nm、578nm、602nm、612nm、650nm、661nm、667nm、668nm、678nm、695nm、702nm、711nm、719nm、737nm、785nm、786nm或805nm处的发射的光。
在一些实施方案中，样品发射的荧光辐射可以以各种不同的方式(例如，沿着垂直于激发光束的传播方向的光学路径)收集。在其它实施方案中，在外延方向上检测荧光辐射。
再次参考图1，在该实施方案中，由照射样品中存在的一种或多种荧光团发射的荧光辐射穿过物镜52并透过分色镜58到达光电探测器64。更具体地，在该实施方案中，透镜65将透过分色镜58的荧光辐射聚焦到狭缝孔66上。透过狭缝的荧光辐射通过荧光发射滤波器68到达光电探测器64。设置在光检测器前面的狭缝孔66(或下面讨论的其它实施方案中的滤波器)基本上允许从样品的特定平面发射的荧光辐射通过，同时拒绝面外荧光发射。另外，荧光发射滤波器68(例如，通带滤波器)允许荧光辐射通过光电探测器64，同时基本上阻挡其它频率的辐射通过。
光电探测器64具有足够的RF带宽，以检测和发送来自整个拍频范围的信号。合适的光电探测器的一些实例包括但不限于光电倍增管、雪崩光电二极管、PIN光电二极管和混合光电探测器等。举例来说，在一些实施方案中，可以使用由Hamamatsu Corporation销售的光电倍增管(例如，R3896、R10699、H11462)。光电探测器响应于接收到的荧光辐射的检测，在该实施方案中产生信号，诸如模拟信号。
作为另一个实例并参考图9A，由响应于LO光束和空间上分离的RF梳状光束同时照射的样品发出的荧光辐射穿过物镜52和分色镜58，通过透镜100耦合到多模光纤102上，所述多模光纤从近端102a延伸到远端102b。更具体地，光纤102的近端102a被置于透镜100的焦平面附近，以便接收荧光辐射。耦合到光纤的远端102b的输出耦合透镜104使离开光纤的辐射准直。
在许多情况下，照射样品的激发辐射激发多个荧光团(例如，有机荧光团)，所述荧光团可具有足够宽的辐射吸收光谱，使得激发频率落在样品中多个荧光团的吸收光谱内。然后，每个荧光团将以不同的频率发射荧光辐射。在不失一般性且为了说明的目的的情况下，在该实施方案中，检测***包括四个光电倍增管106、108、110和112，其每一个接收对应于由四个荧光团之一发射的荧光辐射的准直辐射的一部分，所述荧光团被经照射的样品中的激发辐射激发。更具体地，分色镜114将由荧光团之一以第一频率发射的荧光辐射反射到光电倍增管106，同时允许其它频率的荧光辐射通过。另一个分色镜116将由不同荧光团以不同的第二频率发射的荧光辐射反射到光电倍增管108，同时允许包含由另一种荧光团以第三频率发射的荧光辐射的其余辐射到达第三分色镜118，其将该荧光辐射反射到光电倍增管110。分色镜118允许包括其余辐射(包括由第四荧光团以第四辐射频率发射的荧光辐射)穿过以到达光电倍增管112。
每个以四个荧光频率之一为中心的多个带通滤波器120、122、124和126分别放置在光电倍增管106、108、110和112的前面。以下面讨论的方式分析由每个光电倍增管检测到的信号，以相应的荧光频率生成荧光图像。在一些实施方案中，不是使用多个光电探测器，而是可使用单个光电探测器，例如单个光电倍增管来检测荧光辐射，例如对应于来自单个荧光团的发射的荧光频率。
在一些实施方案中，当样品流过流动池时，照射样品的不同水平行，并且通过一个或多个光电探测器(诸如光电倍增器106、108、110和112)检测与每个水平行相关的荧光辐射。
图9B示意性地描绘了检测***，所述检测***除在该***中不使用光纤外与上文结合图9A所讨论的检测***类似，包含通过分色镜58的来自多个荧光团的荧光发射的荧光辐射在自由空间中传播以到达光电倍增管106、108和112。更具体地，透镜100将荧光辐射聚焦到设置在透镜100与104之间的孔126上，其中所述孔可以拒绝离焦辐射。透镜104准直通过孔径的辐射，其中经准直的辐射以上面结合图9A所讨论的方式分布在光电倍增管之间。
在一些实施方案中，***10可被配置来使用激发辐射提供(没有样品的流动池的)样品的暗场图像和明场图像。举例来说，图9C示意性地描绘了***10的实施方案，所述***包括两个检测臂200和202，用于分别检测样品的暗场图像和明场图像。
更具体地，检测臂200垂直于激发辐射的传播放置，以便接收一部分由流过流动池的样品散射的激发辐射。检测臂200包括两个透镜204和206，它们共同将由样品散射到由透镜204对向的立体角中的激发辐射的至少一部分引导到光电倍增管208上。更具体地，透镜204准直接收的散射辐射，并且透镜206将准直的散射辐射聚焦到光电倍增管208上。在该实施方案中，适当的带通滤波器210设置在光电倍增管208的前面，以允许具有所需频率的辐射通过光电倍增管208，同时阻挡不需要的频率的辐射。光电倍增管208的输出可以以本领域已知的方式处理，例如通过诸如下面讨论的分析模块来处理以生成暗场图像。
检测臂202又包括两个透镜212和214，其中透镜212以向前方向(基本上沿着进入流动池54的激发辐射的传播方向)准直离开流动池的激发辐射，并且透镜214将准直的辐射聚焦到光电探测器216上。设置在光电探测器前面的适当滤波器218(例如带通滤波器)允许将激发频率传输到光电探测器216，同时基本上阻挡其它辐射频率。可以以本领域已知的方式处理光电探测器216的输出，以生成流动池的明场图像。
因此，检测臂200检测由流过所述池的流体散射的激发辐射，并且检测臂202检测传输通过流动池的激发辐射。当没有流体流过流动池时，由光电倍增管208检测到的信号是低的，并且由光电探测器216检测到的信号是高的，因为通过流动池的激发辐射几乎没有散射，因此很大百分比的激发辐射，在某些情况下所有激发辐射，都通过流动池传输。相反，流体样品通过流动池的流动可以使由光电倍增管208产生的信号由于样品的一部分激发辐射的散射而增加，并且由光电探测器216产生的信号随着通过流动池传输的激发辐射的水平降低而减弱。
通过进一步的实例并参考图9D，在根据本教导的***的一个实施方案中，以基本垂直于激发辐射的传播方向的方向相对于流动池54放置的检测臂220a，包括用于检测由样品中的多个荧光团发射的荧光辐射以及由样品散射的激发辐射的光电探测器。更具体地，检测臂220包括透镜222和224，其将荧光辐射以及散射的激发辐射引导到孔226上，孔226拒绝未聚焦的辐射。透镜228准直穿过孔的辐射。分色镜230将激发频率处的辐射的部分反射到光电倍增管232上，用于检测暗场图像，同时允许荧光辐射通过。设置在光电倍增管232前面的适当滤波器232a，例如带通滤波器，允许激发频率的辐射通过至光电倍增管232，同时阻挡不需要的辐射频率。另一个分色镜234将荧光团以第一频率发射的荧光辐射反射到光电倍增管236上，同时允许其它荧光团以其它频率发射的荧光辐射通过。另一个分色镜238将由另一荧光团以第二频率发射的荧光辐射反射到光电倍增管240上，同时允许另一荧光团以第三频率发射的荧光辐射通过，其中光电倍增管242检测到所述荧光辐射。类似于前面的实施方案，多个滤波器236a、240a和242a分别设置在光电倍增管236、240和242的前面，以允许以所需频率传输辐射，同时基本上阻挡不需要的辐射频率。
继续参考图9D，根据本教导的***的该实现还包括另一个检测臂220b，用于例如以结合图9C所讨论的方式生成明场图像。更具体地，检测臂202包括两个透镜212和214，其将光聚焦到光电探测器216上，以生成激发辐射的明场图像。滤波器218(例如，带通滤波器)放置在光电探测器216的前面，以允许激发辐射通过探测器，同时拒绝不需要的辐射频率。
再次参考图1以及图10，在该实施方案中，可将跨阻抗放大器70耦合到光电探测器64(或者结合图9A-9D讨论的光电探测器中的每一个)的输出以放大由光电探测器产生的信号。数据分析单元72(在本文中也称为分析模块或分析仪)接收放大的信号并分析信号以生成样品的荧光图像。数据分析单元72可以用硬件、固件和/或软件实现。举例来说，用于分析检测到的荧光数据的方法可以存储在分析模块的只读存储器(ROM)单元中，以在处理器的控制下访问以分析接收的荧光信号。
如下面更详细讨论的，分析方法确定时变光电探测器输出的频率分量，并基于那些频率分量构建样品的荧光图像。可采用用于测定光电探测器输出的频率含量的各种方法。此类合适方法的一些实例包括但不限于傅里叶变换、锁定检测、滤波、I/Q解调、零差检测和外差检测。
举例来说，图图11A和11B显示可由分析模块72执行以生成样品的荧光图像的示例性分析步骤。在步骤(1)中，将模拟放大信号数字化以产生数字化荧光数据。在步骤(2)中，选择数字化数据的适当部分(长度)用于分析。例如，可以选择对应于照射的样品行(本文中也称为框架)的荧光数据用于分析。或者，可选择数据帧的一部分用于分析。
在步骤(3)中，执行所选数据的傅里叶变换。举例来说，在一些实施方案中，执行数据的快速傅里叶变换(FFT)以测定数据的频率分量。在一些此类实施方案中，FFT的仓(bin)可对应于为数据采集选择的频率。例如，对于256MHz采样率，256个样品可以产生例如从DC至128MHz的彼此分开1MHz的频率仓(frequency bin)。FFT分析提供对应于发射的荧光发射呈现幅度调制时所处的拍频的频率。
继续参考图图11A和11B，在该实施方案中，在步骤(4)中，通过获得该频率分量的实部和虚部分量的平方和的平方根来计算FFT数据中存在的每个频率分量的幅度的度量。由于每个频率分量对应于用于从样品的特定位置引发荧光辐射的拍频之一，因此频率分量的幅度的测量可提供与沿水平样品行的该频率分量相关联的位置的像素值。以这种方式，可以测定水平样品行的图像的像素值。当样品在垂直方向上流过流动池时，对于针对每个水平样品行获得的荧光数据可以重复上述步骤。像素值可用于构建荧光图像(步骤5)。
如上所述，分析模块72可以使用本领域已知的技术并根据本教导以硬件、固件和/或软件实现。举例来说，图12示意性地描绘了分析仪72的示例性实现，所述分析仪包括模数转换器74，用于接收来自放大器70的放大的荧光信号并将该信号数字化以生成数字化的荧光数据。分析模块还包括中央处理单元(CPU)76，其用于控制分析模块的操作，包括执行计算和逻辑操作。分析模块还包括ROM(只读存储器)78、RAM(随机存取存储器)80和永久存储器82。通信总线84便于分析模块的各种组件之间的通信，包括CPU 76与其它组件之间的通信。存储器模块可用于存储用于分析荧光数据和分析结果的指令。举例来说，在一些实施方案中，用于数据分析的指令，例如，用于执行结合图11A和11B所讨论的上述步骤的指令，可以存储在ROM 78中。CPU可采用存储在ROM 78中的指令来操作存储在RAM 80中的数字化荧光数据，以生成样品的荧光图像(例如，一维或二维图像)。CPU可以将荧光图像存储在永久存储器82中，例如，存储在数据库中。如图12中示意性所示，分析模块可以任选地包括图形处理单元(GPU)76'，其用于从接收的数据(例如，荧光数据)执行像素强度和其它量的计算。
在一些实施方案中，可以使用锁定检测技术来实现由光电探测器产生的输出信号的频率解调。举例来说，参考图13A和13B，在一个此类实施方案中，将放大的荧光信号数字化(步骤1)并产生数字化荧光信号的若干个拷贝(步骤2)，其中数字化拷贝的数量(N)对应于与RF梳状光束相关的频率的数量。将信号的每个数字化拷贝与正弦和余弦波相乘以产生多个中间信号(步骤2)，所述正弦和余弦波的频率对应于等于RF梳状光束与LO光束之一的频率之间的差的拍频。使每个中间信号通过低通滤波器(步骤3)，其带宽等于RF梳状光束频率之间的频率间隔的一半。
对于对应于RF梳状频率之一的每个拍频(换句话说，对于对应于被照射样品的空间位置的每个频率)，获得对应于该频率的两个滤波后的中间信号的平方和的平方根作为对应于由LO光束和具有该频率的RF梳状光束照射的样品位置的图像像素的幅度的度量(步骤4)。在一些实施方案中，可以以上面讨论的方式处理对应于相同拍频(即，对应于相同样品位置)的多个荧光数据信号，并且可以对像素值进行平均以便获得平均像素值。
当样品在垂直方向上流过流动池时，对于针对每个水平样品行获得的荧光数据可以重复上述步骤。像素值可用于构建荧光图像(步骤5)。
上述锁定检测方法可用软件、固件和/或硬件实现。举例来说，在一个实施方案中，上述锁定检测方法可使用现场可编程门阵列(FPGA)来实现，特别是如果用多于6个频率时。在一些实施方案中，可以采用多频锁定放大器，诸如由瑞士苏黎世的苏黎世仪器公司(Zurich Instruments of Zurich，Switzerland)销售的HF2L-MF多频放大器。
通过进一步的实例，在一些实施方案中，可通过采用基于带通滤波器的图像解调技术来实现检测到的荧光信号的频率解调。参考图如图14A和14B，在此类频率解调方法的一个实施方案中，由光电探测器64和放大器70提供的荧光信号被数字化(步骤1)并且产生数字化信号的若干拷贝(步骤2)，其中数字化拷贝的数量(N)对应于与RF梳状光束相关的频率数。通过使数字化荧光信号通过以与RF梳状光束之一相关的拍频(即，与样品的特定位置相关联的拍频)为中心的带通滤波器来过滤所述数字化荧光信号的每个拷贝(步骤3)。更具体地，每个带通滤波器以N拍频之一为中心并且具有等于相邻拍频之间的频率间隔的一半的带宽。
每个拍频处的包络检测器用于为每条水平线估计对应于该频率的每个像素的幅度(步骤4)。在一些情况下，将通过处理多个对应于与该像素相关的样品位置的荧光信号而获得的多个对应于像素的像素值进行平均以获得平均像素值。当样品在垂直方向上流过流动池时，对于针对每个水平样品行获得的荧光数据可以重复上述步骤。像素值可用于构建样品的一维或二维荧光图像(步骤5)。
分析模块还可被配置来接收和处理明场和暗场图像数据。例如，参考图9C和图10，分析模块72还可以被配置来从光电探测器208和218接收暗场和明场图像数据，以生成暗场和明场图像。例如，参考图12所示，例如可将以本领域已知的方式生成暗场和明场图像的指令存储在永久存储器82中。处理器76可采用这些指令来处理接收的暗场和明场图像数据以生成图像。分析模块还可被配置为通过覆盖例如荧光图像以及明场和暗场图像中的一种或两者来生成合成图像。
荧光图像以及由根据本教导的***(诸如上述***10)生成的明场和暗场图像可以用于各种不同的方式。例如，可对荧光图像进行积分以生成与通过常规流式细胞仪生成的数据相当的值。还可分析荧光图像以确定生成该图像的荧光探针的位置(例如，可确定探针是在细胞核、细胞质、定位于细胞器中还是存在于细胞膜外部)。另外，在一些应用中，通过检测不同荧光带获得的多个荧光图像(所述图像全部取自同一细胞)，可用于测定细胞内多个荧光探针的共定位程度。另外，使用多色荧光、明场和暗场图像，可能分析细胞形态、细胞信号传导、内化、细胞间相互作用、细胞死亡、细胞周期和斑点计数(例如，FISH)等等。
如上所述，***10可以以至少三种不同的模式操作。在上面讨论的一种模式中，LO光束和多个RF梳状光束同时照射样品的一部分(例如，沿水平范围设置的位置)，并且检测和分析从照射位置发射的荧光辐射，以便构建样品的荧光图像。在另一种操作模式中，不是将多个RF驱动信号同时施加到AOD，而是将包含驱动信号的频率斜坡施加到AOD，使得激光束的频率随着时间从起始频率改变(f1)至一个结束频率(f2)。对于频率斜坡中的每个驱动频率，激光束的频率偏移了该驱动频率，并且通过频移激光束照射样品以引发来自样品的荧光辐射。换句话说，在该模式中，操作所述***以通过用从中心激光频率偏移的多个频率在时间间隔内连续照射样品来获得来自样品的荧光辐射。由AOD产生的频移伴随着角度偏转，使得使用相同的光路，光束以高速扫描过样品。
更具体地，在该操作模式中，RF频率合成仪10用于将驱动信号以从起始频率(f1)至结束频率(f2)的梯度施加到AOD 18。举例来说，驱动信号斜坡的频率范围可以为约50MHz至约250MHz。在一些实施方案中，驱动信号从约100MHz斜坡上升至约150MHz。在该实施例中，驱动频率随时间连续变化，例如，以实现高速。在其它实施方案中，可以以从起始频率(f1)至结束频率(f2)的离散的步进改变驱动频率。
选择驱动频率，使得频移光束将错过反射镜28并沿着由透镜26、透镜30、反射镜40/42、分束器44、透镜46、反射镜56、透镜50、反射镜58和物镜52限定的光路传播，以照射流过样品架的样品的一部分。斜坡率优选足够快，以便改善并优选地防止当样品流过光束时基于发射的荧光辐射生成的荧光图像的垂直方向上的任何模糊。例如，这可以通过将斜坡率与样品的流速相匹配来实现。样品处的激光光斑尺寸可用于估计适当的速率。举例来说，对于1微米的激光光斑尺寸，对于0.1米/秒的样品流速，1行的扫描时间应为10微秒或更短，以避免图像模糊。
以上述方式被收集和检测响应于激发辐射的照射而从样品发射的荧光辐射。具体而言，参考图10，荧光辐射由光电探测器64检测。检测到的荧光由放大器70放大，并且放大的信号由分析模块72分析以重建样品的荧光图像。通过向从开始频率(f1)至结束频率(f2)的扫描周期内的特定时间分配水平像素位置来进行图像的重建。与在上述操作模式中分析频率分量的幅度以获得像素值相反，在该操作模式中使用的解调方法仅使用检测到的荧光信号的时域值来将值分配给图像的像素。当样品在垂直方向上流动时，可以重复该过程，以获得样品的二维荧光图像。
由样品发射的荧光辐射(如果有的话)由光电探测器64收集。参考图10，检测到的荧光辐射由放大器70放大。分析模块72接收放大的信号。在该操作模式中，分析模块分析荧光信号以测定样品(例如细胞/颗粒)的荧光含量。由于在该操作模式中仅有一个光束激发样品，因此响应于激发样品而不产生拍频。因此，在荧光信号的频域中没有图像信息。相反，检测到的荧光信号具有在时域中编码的图像信息。在该操作模式中，可使用检测到的荧光信号的时间值作为水平像素坐标，并且将荧光信号的数字化电压值作为像素值(亮度)来数字地重建图像。应用于AOD的驱动频率的每次扫描生成图像的一条水平线(行)。当样品流过照明区域(点)时，通过连续扫描实现图像重建。
在另一种操作模式中，可操作***10以通过单个激发频率同时照射样品的多个位置，该激发频率可以例如通过使激光束的中心频率偏移一个射频来产生。更具体地，再次参见图1，在此类操作模式中，可以将单个驱动射频施加到AOD 18，以产生具有相对于进入AOD18的激光束偏移的频率的激光束。另外，频移激光束相对于进入AOD的激光束呈现角度偏移，使得射频激光束被截取并被反射镜28反射，经由透头32以及反射镜33和35反射到顶帽式式光束整形器34。离开顶帽式光束整形器的光束被分束器44反射，并被透镜46聚焦到中间像平面48上。在该平面中，如图15A中示意性所示，激光束1000显示沿水平方向的拉伸轮廓。
水平拉伸的激光束被反射镜56反射到正透镜50。在通过透镜50之后，激光束被反射镜58反射到物镜52。如上所述，正透镜50和物镜52形成望远镜，用于将来自中间像平面48的顶帽形轮廓激光束中继到流过流动池54的样品上。
水平拉伸的激光束照射样品的水平范围，以沿着该水平范围激发目标荧光团(如果存在于样品中的话)。因此，在该操作模式中，与以不同激发频率照射样品的多个水平位置的第一操作模式不同，以相同的激发频率照射样品的多个水平位置。该操作模式不能使用户获得流过的细胞或颗粒的图像。然而，在该操作模式中，通常可向样品施加比在其它两种操作模式中更高的光功率，如果不需要图像，这对于获得更高的信噪比数据是有用的。仅通过改变驱动声光偏转器的电子信号就可以实现该操作模式，而无需对***进行任何机械改变。
因此，***10可以以三种不同的操作模式操作以引发来自样品的荧光辐射。
在一些实施方案中，可以例如通过将每个射频移位和本地振荡器光束的节拍的相位与检测到的荧光信号中的相应射频分量的相位进行比较，来在样品上的每个空间位置处进行荧光寿命测量。举例来说，图15B显示了***10'，即上面讨论的***10的改进形式，其允许这种荧光寿命测量(为简洁起见，在该图中未示出图1中所示的某些组件)。具体地，入射在分束器44上的RF梳状光束的一部分被分束器反射到会聚透镜400上(作为说明，在该实施方案中，透镜400具有200mm的焦距，但是也可以使用其它焦距)。透镜400将RF梳状光束的该部分聚焦到光电二极管402上，光电二极管402检测激发光束。光电二极管402的输出可由分析模块72(参见，图10)接收。分析模块可例如使用上面讨论的解调技术之一来提供激发光束的频率解复用，并测定激发光束中的每个射频分量的相位。这可为检测到的荧光信号中的每个射频分量提供参考相位，利用该参考相位可以比较该射频分量的相位。例如，可采用激励信号的FFT的实部和虚部分量或锁定型解调的I和Q分量。或者，检测样品/流动细胞的明场图像的检测器的输出可用于获得参考相位，利用该参考相位可以比较荧光拍频的相位。
更具体地，分析模块72可以例如以上面讨论的方式提供检测到的荧光信号的频率解复用。如本领域技术人员将理解的，对于荧光信号中的每个拍频，可以将射频分量的相位与激发光束的相应参考相位进行比较，以获得空间分辨的荧光寿命测量和荧光寿命图像。
在某些实施方案中，主题***包括采用上述光学配置的流式细胞术***，用于检测流动流中样品发出的光。在某些实施方案中，所述主题***是流式细胞术***，其包括以下美国专利号(其公开内容通过引用整体并入本文)中描述的流式细胞仪的一种或多种组件：3,960,449；4,347,935；4,667,830；4,704,891；4,770,992；5,030,002；5,040,890；5,047,321；5,245,318；5,317,162；5,464,581；5,483,469；5,602,039；5,620,842；5,627,040；5,643,796；5,700,692；6,372,506；6,809,804；6,813,017；6,821,740；7,129,505；7,201,875；7,544,326；8,140,300；8,233,146；8,753,573；8,975,595；9,092,034；9,095,494和9,097,640。
如上所述，在一些实施方案中，主题***被配置用于对作为流动流(诸如在流式细胞仪的流动流中)流动的样品中的颗粒(例如，细胞)成像。流动流中的颗粒流速可以是0.00001m/s或更大，诸如0.00005m/s或更大，诸如0.0001m/s或更大，诸如0.0005m/s或更大，诸如0.001m/s或更大，诸如0.005m/s或更大，诸如0.01m/s或更大，诸如0.05m/s或更大，诸如0.1m/s或更大，诸如0.5m/s或更大，诸如1m/s或更大，诸如2m/s或更大，诸如3m/s或更大，诸如4m/s或更大，诸如5m/s或更大，诸如6m/s或更大或0μL/分钟或或更大，诸如7m/s或或更大，诸如8m/s或或更大，诸如9m/s或或更大，诸如10m/s或或更大，诸如15m/s或或更大，包括25m/s或更大。例如，取决于流动流的尺寸(例如，流动喷嘴孔)，流动流可以在主题***中具有0.001μL/分钟或更大，诸如0.005μL/分钟或更大，诸如0.01μL/分钟或更大，诸如0.05μL/分钟或更大，诸如0.1μL/分钟或更大，诸如0.5μL/分钟或更大，诸如1μL/分钟或更大，诸如5μL/分钟或更大，诸如10μL/分钟或更大，诸如25μL/分钟或更大，诸如50μL/分钟或更大，诸如100μL/分钟或更大，诸如250μL/分钟或更大，并且包括500μL/分钟或更大的流速。
如上所述，可采用许多技术从通过询问样品获得的时间数据(信号)诸如荧光数据重建图像(例如，荧光图像)。在一些实施方案中，采用正演模型，其基于底层图像提供预期信号(例如，荧光信号)以分析测量信号，以便获得最可能已经生成该信号的底层图像。在一些实施方案中，可反转正演模型以产生反向模型，用于对测量信号进行操作以获得图像。例如，正演模型可以将底层图像处理为频谱图(spectrogram)的幅度。在此类情况下，可通过对图像应用逆傅里叶变换来获得生成图像的信号。然后，反向模型需要将傅里叶变换应用于测量信号以获得图像。
在一些实施方案中，构建指示根据本教导的细胞计数***如何将由彼此相隔一个或多个RF频率的多个光束照射的样品的底层“真实图像”转换为待测量的时间信号的模型。通过反转该模型并结合***的已知或测量的参数，可以将测量的信号转换成更接近地匹配底层真实图像的估计。通过进一步说明，图20示意性地显示如何构建此类模型。生成样品的“真实图像”(表示为二维强度图(I(x,y)))，及其至测量的信号(y(t))的转换受细胞计数***的部件诸如照明***、光学收集***和检测臂中使用的电子***的多个特性影响。举例来说，照明***可包括，除其它以外，激光束、接收激光束的声光偏转器以及用于将RF信号施加到偏转器以产生频率上相隔一个或多个RF频率的多个光束的RF发生器。关于与这些***相关的参数(由椭圆标示的)的知识提高了估计样品图像的能力。例如，在一些实施方案中，可以从不包含图像信号的电子部件的测量的输出响应的傅里叶谱的部分测量电子部件的响应。可以假设响应滚降的模型，并且可将模型与测量的响应拟合。
样品上的照射点可以通过具有略微不同的频率的两个光束的光学干涉产生。光束的路径长度的差异可以导致样品处的拍频的相移。在许多实施方案中，即使由于***的机械振动引起的(例如，由冷却风扇引起的)微小路径长度变化也会影响样品处的光束的相位。为了部署锁定检测，必须考虑这些相位变化。在一些实施方案中，可使用明场信号的相含量(phase content)来测量每个事件之前和之后的相移，该相含量通常具有高信噪比(SNR)并且否则不应该改变相位。所得的相移可被建模为信号的傅里叶分量的旋转，将其在频域中分出来。
在一些实施方案中，与电子***相关的散粒噪声本身可在图像中表现为成行的高斯白噪声，其与该行中存在的信号量成比例。这与常规成像中的大多数噪声模型不同，所述常规成像假设独立高斯噪声或逐像素泊松噪声(每个像素中的噪声与该像素的强度成比例)。在一些实施方案中，假设图像被具有该特性的噪声破坏，则可通过找到底层图像的最大似然性来改善图像质量。一种技术应用了自然发生的图像与像素没有相关性的白噪声不同的事实。此类技术可以假设真实图像包含许多具有相对平坦的强度的区域。这表示估计是一个优化问题；需要找到尽可能接近测量的信号的图像的估计值，但最小化估计值内的变化量。在一些实施方案中，这可通过以下成本函数来实现：
||G-y||+λ(r)*TV(G)，
其中||G-y||表示图像的估计(G)与测量的图像之间的距离加上猜测G的总变化(由惩罚项λ加权而来，其是行依赖性的)。可最小化该成本函数以获得样品的最佳拟合图像。任何优化方法都可用于进行此最小化(包括粒子群优化、遗传算法或梯度下降)。
如上所述，在某些情况下，正演模型还可并入对信号的各种噪声贡献的模型，诸如散粒噪声模型和来自光电倍增管和/或电子放大器的噪声贡献(诸如乘法噪声、热噪声、谐波失真和非线性失真)的模型。
在一些实施方案中，可采用统计技术将模型参数与测量数据进行拟合。在一些此类实施方案中，不能直接反转正演模型。在一些情况下，此类统计技术作为最大似然估计器操作，以生成最有可能生成测量的数据的模型参数。另外，在一些实施方案中，使用反向模型(例如，上文讨论的FFT或基于锁定的处理)重建的图像可用作用于此类统计技术的初始(种子)图像，以进一步增强图像的信噪比和/或分辨率。
举例来说，参考图19A的流程图，在根据本教导的用于进行流式细胞术的一个示例性方法中，用激光辐射照射样品，并且检测响应于照射而从样品发出的至少一部分辐射以产生对应于检测到的辐射的时间信号。进行基于正演模型的信号的统计分析，以便重建样品的图像。图像可以是样品的荧光、暗场和明场图像中的任一种。
通过进一步说明，下面提供了此类统计技术的一些实例。
此类统计技术的一个实例是最小二乘优化，其可将***/过程建模为多个具有正态分布误差的线性分量的总和。举例来说，参考图19B，细胞计数信号(例如，荧光、明场和/或暗场)可通过将信号分解成多个短时间片段(例如，持续约1-100微秒的片段)来建模(步骤1)。在一些实施例中，每个片段可被建模为具有一个正弦项和一个余弦项(对于每个图像像素)以及零阶、一阶和二阶项，导致(2xN+3)xM参数拟合NxM像素图像。另外，还可对每个片段进行建模以包括噪声，例如，正态分布的噪声。然后可采用最小二乘回归分析(步骤2)以通过最小化对应于由模型预测的时间片段与与测量的细胞计数信号相关的相应时间片段之间的差的残差平方和，来获得模型的参数的估计值。用于重建细胞计数图像的最小二乘优化方法可以是快速的并且可以例如使用矩阵运算来实现。
作为另一个实例，参考图19C，可用激光辐射照射样品(步骤1)，并且可检测响应于其激光辐射照射而从样品发出的辐射(例如，在该实施方案中为荧光辐射)，以产生对应于检测到的辐射的荧光信号(步骤2)。例如，上述***可用于获取射频标记的荧光数据。可采用最小二乘分析，通过最小化对应于从模型推断的荧光信号与测量的荧光信号之间的差的残差平方和来获得与样品的像素化图像(例如，在该实施例中为样品的荧光图像)的模型相关的参数的最佳估计(步骤3)。
在一些实施方案中，梯度下降优化方法可用于重建细胞计数图像，例如荧光、明场和/或暗场图像。该方法允许采用非线性模型，从而通常以计算时间为代价减少了自由参数的数量。利用较少的自由参数，模型方差可以减小图像的信噪比(SNR)，从而改善其信噪比。梯度下降是迭代优化方法，其可始于图像(例如，荧光图像)的初始估计。更具体地，参考图19D的流程图，在此类实施方案中，可用激光辐射照射样品(步骤1)，并且可检测响应于照射(例如，荧光辐射)而从样品发出的辐射以生成时间数据(步骤2)。可处理时间数据以生成样品图像的初始估计。例如，在一些实施方案中，可通过将FFT或上述最小二乘法应用于从样品获得的时间数据(例如，荧光数据)来获得此类估计。迭代地改进图像的初始估计以获得样品的最终图像。在每次迭代时，计算指示预期信号与测量的信号之间的距离的误差函数和相关的误差梯度，并且使图像逐步降低误差梯度。重复迭代过程直到误差最小化，从而生成最佳拟合图像(步骤3)。可以例如基于计算时间与精度之间的折衷来选择用于使图像沿着误差梯度向下步进的步长。
在一些情况下，梯度下降优化方法可能无法达到最佳全局图像，因为优化过程可能会陷入参数空间的局部最小值。因此，在一些实施方案中，粒子群优化可用于重建细胞计数图像。类似于梯度下降优化方法，粒子群优化可始于图像的估计(例如，通过对检测到的荧光数据应用FFT或最小二乘法获得的荧光图像的估计)。然后，粒子群优化采用多个(“群体”)模型，所述模型中的每个模型被称为粒子，其以偏向其负误差梯度的方式对参数空间进行随机地采样。通过控制随机采样的程度，群体可以以较慢的收敛为代价来逃避局部最小值。
在一些实施方案中，将遗传算法优化用于图像的重建。在此类方法中，采用具有随机变化参数的模型群。例如，基于预测的图像数据与测量的图像数据之间的误差程度来评估每个模型“适合度”。然后选择具有大于阈值的“适合度”的模型，并且随机组合它们的参数以产生新一代模型。这些后代经历了相同的过程，直到最适合模型仍然存在。特别地，图像中的像素是遗传算法中使用的参数。***的正演模型将该组参数转换为预期的时间信号，然后使用预期信号和测量的信号之间的差来评估模型的适合度。产生与测量的信号非常不同的时间信号的图像被认为不太合适。
在一些实施方案中，可将关于图像信号(例如，荧光、明场或暗场信号)的先验信息与贝叶斯频谱估计一起使用以获得高分辨率图像。举例来说，在一些实施方案中，可以知道图像信号由许多具有未知频率和幅度的正弦波组成。在一些此类实施例中，贝叶斯频谱估计方法可用于以比使用常规功率谱方法可获得的更高的准确度和精度来估计那些正弦波的参数。这可反过来导致比常规可能的更高的频率分辨率，从而允许超分辨率成像。换句话说，在此类实施方案中，时域信号可以被处理为具有有限数量的参数的多个独立测量，并且可获得参数的最大似然估值。
用于图像重建的上述方法可用硬件、固件和/或软件实现。举例来说，上述数据分析模块72(参见例如图10)可被配置来使用统计分析技术进行上述图像重建方法。举例来说，包含用于进行上述方法的必要步骤的程序可以存储在分析模块的存储器单元(例如ROM)中。然后，模块的处理器可以访问程序以将图像重建方法的步骤应用于从研究中的样品接收的时间数据，例如荧光数据。
提供以下实施例仅用于进一步阐明本教导的各个方面，并且不意图必须说明实施本发明的教导的最佳方式或可获得的最佳结果。
实施例1
将***用于测量来自利用8个离散水平的荧光染料染色的聚苯乙烯珠粒的荧光辐射，所述***类似于上面结合图1描述的具有类似于结合图9A描述的检测***的检测***的***，所述荧光染料由Spherotech Inc.of Lake Forest,IL以商品名RCP-30-5A销售。该***还用于以上述方式生成明场和暗场图像。
图16A是暗场强度对比明场强度的散点图。将该图的矩形部分用作门以生成图16B、16C和16D中的数据，其包含所测量的所有事件的约32％(检测到总共50,000个事件)。图16B显示由每个颗粒发射的红色荧光(PI)对绿色荧光(FITC)的散点图。该图清楚地显示了具有不同亮度水平的8个群体。图16C和图16D是相同数据的直方图。
实施例2
使用***获得用CD45-FITC和碘化丙啶染色的固定外周血白细胞的FIRE、明场和暗视野图像，所述***类似于上面结合图1描述的具有类似于上面结合图9A讨论的检测***的检测***的***。样品还含有用钙黄绿素-AM染色的活HeLa细胞的一部分。在数据采集期间，将细胞以0.5米/秒的速率流过流动池。
图17A中所示的图像从上到下为：明场、CD45-FITC和碘化丙锭荧光通道的叠加、明场、暗场、CD45-FITC和PI通道荧光。未应用任何补偿，并且所有图像都自动缩放以供查看。编号为2、3、5、8、9、12、17、20和21的细胞为HeLa细胞(B群)，其它为白细胞(A群)。
图17B是散点图，其中群体A代表白细胞，群体B代表HeLa细胞。
实施例3
使用***获得用CD45-FITC和碘化丙啶染色的固定外周血白细胞的FIRE、明场和暗视野图像，所述***类似于上面结合图1描述的具有类似于上面结合图9A讨论的检测***的检测***的***。将样品掺入固定的用抗-EpCAM-FITC和碘化丙啶(PI)染色的MCF-7细胞的一小部分。在数据采集期间，将细胞以0.5米/秒的速率流过流动池。
图18A中所示的图像从上到下为：明场、FITC和PI通道荧光的叠加、明场、暗场、FITC和PI荧光。在白细胞群中，绿色荧光是来自PI染色的荧光漏入的伪影。所有图像都在亮度上自动缩放以供观察，因此，白细胞似乎表现出FITC荧光，但这是来自PI的少量荧光漏入的信号。编号为1、2、4、5、10、13、15和16的细胞是MCF-7细胞。
图18B是散点图，其中群体A代表白细胞，群体B代表MCF-7细胞。
实施例4
图21A描绘了细长的眼虫(藻类)细胞的消光和散射的高信噪比(SNR)FIRE测量，使得细胞的存在导致信号减少。图21B显示通过提取图21A中描绘的信号的傅里叶分量生成的眼虫细胞的图像。这是细胞的可接受图像。然而，通过以上面讨论的方式结合关于FIRE激励信号、电子响应和***的相位响应的信息，可以改善该图像，以生成图21C中描绘的图像。
实施例5
图22A描绘了通过荧光显示的实施例4的眼虫细胞的脂质含量的低信噪比FIRE测量。虽然存在荧光靶标产生的清晰脉冲，但大量的噪声破坏信号。图22B显示提取信号的傅里叶系数的结果；细胞无法识别噪声。结合关于收集电子器件的频率响应、驱动***的基础FIRE信号以及***的相位响应的信息生成图22C中所示的改进的但有噪声的图像。以上面讨论的方式，更复杂的基础信号模型和光子散粒噪声的影响生成图22D中所示的图像。
本领域普通技术人员将理解，在不脱离本教导的范围的情况下，可以进行各种改变。特别地，可将上面讨论的实施方案的各种特性、结构或特征以合适的方式组合。例如，结合一个实施方案讨论的检测***可以用于另一个实施方案。

Claims (49)

1.一种用于进行流式细胞术的***，所述***包括：
激光器，其用于产生用于照射样品的激光辐射，
至少一个检测器，其被布置来检测响应所述照射而从所述样品发出的辐射的至少一部分并且产生时间信号，
分析模块，其被连接到接收所述时间信号的检测器并基于正演模型对所述信号进行统计分析以重建所述样品的图像。
2.如权利要求1所述的***，其中所述激光辐射包括彼此偏移一个射频的至少两个光学频率，并且所述时间信号是表现出对应于所述至少两个光学频率之间的差的至少一个拍频的时频波形。
3.如权利要求2所述的***，其中所述时间信号是时间荧光信号。
4.如权利要求1所述的***，其中从所述样品发出的所述辐射是荧光辐射，并且所述时间信号是所述分析模块所操作以形成所述样品的荧光图像的荧光信号。
5.如权利要求1所述的***，其中从所述样品发出的所述辐射是散射激光辐射，并且所述分析模块对所述时间信号进行操作以生成所述样品的暗场图像。
6.如权利要求1所述的***，其中从样品发出的所述辐射是透过所述样品的所述照射激光辐射的一部分，并且所述分析模块对所述时间信号进行操作以形成所述样品的明场图像。
7.如权利要求1所述的***，其中所述统计分析基于产生所述模型的一个或多个参数的最大似然估值。
8.如权利要求1所述的***，其中所述统计分析采用最小二乘优化、梯度下降优化、粒子群优化、遗传算法优化、参数估计和贝叶斯频谱估计中的任一种。
9.如权利要求1所述的***，其中所述分析模块形成所述样品的初始图像，并将所述图像用作用于所述统计分析的种子图像，以生成相对于所述种子图像具有改善的信噪比或改善的分辨率的所述重建图像。
10.如权利要求9所述的***，其中所述分析模块采用傅里叶变换和锁定检测中的任一种来形成所述初始图像。
11.如权利要求1所述的***，其中所述分析模块采用所述正演模型来将所述检测的时间信号的模型生成为多个时间片段。
12.如权利要求11所述的***，其中所述统计分析采用最小二乘回归分析，以便通过最小化对应于所述建模的时间片段与相应的测量的时间片段之间的差的残差平方和来获得与所述建模时间片段相关的参数的值。
13.如权利要求12所述的***，其中每个所述时间片段被建模为包括一个正弦项和一个余弦项。
14.如权利要求1所述的***，其中所述分析模块采用最小二乘分析，以通过最小化对应于所述检测到的时间信号与从所述正演模型推断的时间信号之间的差的残差平方和来获得与所述照射的样品的图像的模型相关的一个或多个参数的一个或多个估值。
15.如权利要求1所述的***，其中所述正演模型包括非线性模型。
16.如权利要求15所述的***，其中所述统计分析包括梯度下降优化方法。
17.如权利要求16所述的***，其中所述梯度下降优化方法计算指示基于所述正演模型的预期时间信号与所述测量的时间信号之间的距离的误差梯度，并通过将先前的图像沿着所述误差梯度向下步进来迭代地计算所述样品的更新图像。
18.如权利要求17所述的***，其中所述梯度下降优化方法利用基于所述测量的时间信号计算的初始估计图像开始所述迭代计算。
19.如权利要求18所述的***，其中所述分析模块通过对所述测量的时间信号应用傅里叶变换和锁定检测中的任一种来计算所述初始图像。
20.如权利要求1所述的***，其中所述统计分析采用关于与贝叶斯频谱估计结合的所述测量的时间信号的先验信息以重建所述样品的图像。
21.如权利要求1的***，其中所述先验信息表示所述时间信号由多个具有未知频率和幅度的正弦波组成。
22.如权利要求21所述的***，其中所述贝叶斯频谱估计提供所述正弦曲线的所述未知频率和幅度的估值。
23.如权利要求1所述的***，其中所述统计分析采用粒子群优化方法。
24.如权利要求1所述的***，其中所述统计分析采用遗传算法。
25.如权利要求2所述的***，其中所述射频在约10MHz至约250MHz的范围内。
26.如权利要求1所述的***，其中所述样品包括细胞、微囊泡、细胞碎片、脂质体、珠粒和小生物体中的任一种。
27.如权利要求1所述的***，其中所述激光辐射的频率范围为约300THz至约1000THz。
28.一种用于进行流式细胞术的***，所述***包括：
用激光辐射照射样品，
部署检测器以检测响应于所述照射而从所述样品发出的辐射的至少一部分并且产生对应于所述检测到的辐射的时间信号，以及
利用计算机处理器基于正演模型对所述信号进行统计分析，以便重建所述样品的图像。
29.如权利要求28所述的方法，其中所述图像是荧光、暗场和明场图像中的任一种。
30.如权利要求28所述的方法，其中所述激光辐射包括彼此偏移一个射频的至少两个光学频率，以引发来自所述样品的荧光辐射，其中所述时间信号包含与所述射频分隔的光学频率相关的一个或多个拍频。
31.如权利要求30所述的方法，所述方法还包括处理所述时间信号以生成荧光图像，并且使用所述荧光图像作为种子图像以用于进行所述统计分析。
32.如权利要求31所述的方法，其中所述统计分析采用最小二乘法，以通过最小化对应于所述检测到的时间信号与从所述正演模型推断的相应时间信号之间的差的残差平方和来获得所述重建图像。
33.如权利要求32所述的方法，其中所述时间信号是荧光信号。
34.如权利要求32所述的方法，其中所述时间信号是散射信号。
35.如权利要求30所述的方法，所述方法还包括将所述检测到的时间信号建模为多个时间片段，每个时间片段具有一个正弦项和一个余弦项。
36.如权利要求35所述的方法，其中所述统计分析采用最小二乘回归分析，以便通过最小化对应于所述建模的时间片段与相应的测量的时间片段之间的差的残差平方和来获得所述时间信号模型的参数值。
37.如权利要求28所述的***，其中所述正演模型包括非线性模型。
38.如权利要求37所述的***，其中所述统计分析包括梯度下降优化方法。
39.如权利要求38所述的***，其中所述梯度下降优化方法计算指示基于所述正演模型的预期时间信号与所述测量的时间信号之间的距离的误差梯度，并通过将先前的图像沿着所述误差梯度向下步进来迭代地计算更新图像。
40.如权利要求39所述的***，其中所述梯度下降优化方法利用基于所述测量的时间信号计算的初始估计图像开始所述迭代计算。
41.如权利要求40所述的方法，所述方法还包括使用所述处理器通过对所述测量的时间信号应用快速傅里叶变换(FFT)来计算所述初始估计的图像。
42.如权利要求28所述的***，其中所述统计分析采用关于与贝叶斯频谱估计结合的所述测量的时间信号的先验信息以重建所述样品的图像。
43.如权利要求42的***，其中所述先验信息表示所述时间信号由多个具有未知频率和幅度的正弦波组成。
44.如权利要求43所述的***，其中所述贝叶斯频谱估计提供所述正弦曲线的所述未知频率和幅度的估值。
45.如权利要求28所述的***，其中所述统计分析采用粒子群优化方法。
46.如权利要求28所述的***，其中所述统计分析采用遗传算法。
47.如权利要求30所述的***，其中所述射频在约10MHz至约250MHz的范围内。
48.如权利要求28所述的***，其中所述样品包括细胞、微囊泡、细胞碎片、脂质体、珠粒和小生物体中的任一种。
49.如权利要求28所述的***，其中所述激光辐射的频率范围为约300THz至约1000THz。
CN201780025417.5A 2016-05-12 2017-05-12 具有增强图像分辨率的荧光成像流式细胞术 Active CN109477784B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201662335359P 2016-05-12 2016-05-12
US62/335,359 2016-05-12
PCT/US2017/032417 WO2017197271A1 (en) 2016-05-12 2017-05-12 Fluorescence imaging flow cytometry with enhanced image resolution

Publications (2)

Publication Number Publication Date
CN109477784A true CN109477784A (zh) 2019-03-15
CN109477784B CN109477784B (zh) 2022-04-01

Family

ID=58765980

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201780025417.5A Active CN109477784B (zh) 2016-05-12 2017-05-12 具有增强图像分辨率的荧光成像流式细胞术

Country Status (6)

Country Link
US (3) US10935485B2 (zh)
EP (1) EP3455608A1 (zh)
JP (1) JP7023244B2 (zh)
KR (1) KR102381176B1 (zh)
CN (1) CN109477784B (zh)
WO (1) WO2017197271A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112037217A (zh) * 2020-09-09 2020-12-04 南京诺源医疗器械有限公司 一种基于荧光成像的术中血流成像方法
CN114088606A (zh) * 2021-10-23 2022-02-25 广州市艾贝泰生物科技有限公司 细胞分析装置

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4242631A3 (en) * 2016-09-13 2024-03-06 Becton, Dickinson and Company Flow cytometer with optical equalization
JP6781987B2 (ja) * 2017-02-17 2020-11-11 国立大学法人大阪大学 電磁波検出装置、フローサイトメーター、電磁波検出方法及び電磁波検出プログラム
US10613413B1 (en) * 2017-05-31 2020-04-07 Facebook Technologies, Llc Ultra-wide field-of-view scanning devices for depth sensing
US10181200B1 (en) 2017-06-28 2019-01-15 Facebook Technologies, Llc Circularly polarized illumination and detection for depth sensing
US10574973B2 (en) * 2017-09-06 2020-02-25 Facebook Technologies, Llc Non-mechanical beam steering for depth sensing
WO2019245709A1 (en) * 2018-06-19 2019-12-26 Becton, Dickinson And Company Variable multiplexing switches for detector arrays, systems and methods of use thereof
US10977815B1 (en) * 2018-07-02 2021-04-13 Facebook Technologies, Llc Structured light eye-tracking
CN110888290B (zh) * 2018-09-07 2022-03-04 深圳光峰科技股份有限公司 光源***及投影***
US20220034785A1 (en) * 2018-09-14 2022-02-03 The Regents Of The University Of California Cell sorting device and method
WO2020197787A1 (en) * 2019-03-22 2020-10-01 Becton, Dickinson And Company Spectral unmixing of fluorescence imaging using radiofrequency-multiplexed excitation data
WO2020241652A1 (ja) * 2019-05-30 2020-12-03 ソニー株式会社 光学測定装置及び光学測定システム
KR20220003129A (ko) 2019-05-30 2022-01-07 벡톤 디킨슨 앤드 컴퍼니 무선 주파수 다중화 신호의 위상 보정
US11531112B2 (en) 2019-06-20 2022-12-20 Cilag Gmbh International Offset illumination of a scene using multiple emitters in a hyperspectral, fluorescence, and laser mapping imaging system
US11903563B2 (en) * 2019-06-20 2024-02-20 Cilag Gmbh International Offset illumination of a scene using multiple emitters in a fluorescence imaging system
US20220247962A1 (en) 2019-06-26 2022-08-04 Sony Group Corporation Optical measurement device and optical measurement system
JP7420551B2 (ja) * 2019-12-27 2024-01-23 リオン株式会社 粒子測定装置
US20210278333A1 (en) 2020-01-31 2021-09-09 Becton, Dickinson And Company Methods and systems for adjusting a training gate to accommodate flow cytometer data
US11879827B2 (en) 2020-04-29 2024-01-23 Becton, Dickinson And Company Methods for modulation and synchronous detection in a flow cytometer and systems for same
EP4147029A4 (en) 2020-05-06 2023-11-01 Becton, Dickinson and Company METHODS AND SYSTEMS FOR CHARACTERIZING OVERFLOW SCATTERING IN FLOW CYTOMETER DATA
AU2021275676A1 (en) 2020-05-19 2022-12-08 Becton, Dickinson And Company Methods for modulating an intensity profile of a laser beam and systems for same
CN115916409A (zh) 2020-06-24 2023-04-04 贝克顿·迪金森公司 流式细胞术液滴分配***及使用该***的方法
WO2021262285A1 (en) 2020-06-26 2021-12-30 Becton, Dickinson And Company Dual excitation beams for irradiating a sample in a flow stream and methods for using same
KR20220107378A (ko) 2021-01-25 2022-08-02 금오공과대학교 산학협력단 Led 광원을 이용한 세포분석장치
US20230053122A1 (en) 2021-08-10 2023-02-16 Becton, Dickinson And Company Clamps for operably coupling an optical component to a mounting block, and methods and systems for using the same
US20230046207A1 (en) 2021-08-10 2023-02-16 Becton, Dickinson And Company Outlet fittings for reducing bubbles at the interface with a flow cell, and flow cytometers and methods using the same
KR20230117000A (ko) 2022-01-28 2023-08-07 금오공과대학교 산학협력단 3d 프린터를 이용한 이동성 형광 암세포 분석기

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030031352A1 (en) * 2001-08-10 2003-02-13 Nelson Alan C. Optical projection imaging system and method for automatically detecting cells with molecular marker compartmentalization associated with malignancy and disease
US20050207633A1 (en) * 2003-04-02 2005-09-22 Nick Arini Method of, and computer software for, classification of cells into subpopulations
US20110275558A1 (en) * 2010-05-04 2011-11-10 Virginia Tech Intellectual Properties, Inc. Lanthionine synthetase component c-like proteins as molecular targets for preventing and treating diseases and disorders
US20120128264A1 (en) * 2010-11-23 2012-05-24 General Electric Company Methods and systems of optical imaging for target detection in a scattering medium
CN102650527A (zh) * 2012-05-25 2012-08-29 北京航空航天大学 一种基于时间序列分析消噪的光纤陀螺温度补偿方法
CN102792146A (zh) * 2010-02-10 2012-11-21 赫拉巴Abx公司 用于流体中微颗粒的多参数测量的装置和方法
US20130078625A1 (en) * 2011-09-25 2013-03-28 Theranos, Inc., a Delaware Corporation Fluid handling apparatus and configurations
US20130323825A1 (en) * 2012-05-29 2013-12-05 Sony Corporation Information processing apparatus, information processing method, and program
US20150182136A1 (en) * 2013-12-26 2015-07-02 Fundacio Institut De Ciencies Fotoniques Speckle contrast optical tomography
US20160003741A1 (en) * 2013-01-09 2016-01-07 The Regents Of The University Of California Apparatus and methods for fluorescence imaging using radiofrequency-multiplexed excitation
US20160054293A1 (en) * 2013-03-22 2016-02-25 Map Diagnostics Limited Prenatal Screening
WO2016054293A1 (en) * 2014-09-30 2016-04-07 The Regents Of The University Of California Imaging flow cytometer using spatial-temporal transformation

Family Cites Families (116)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US504732A (en) 1893-09-12 And louis
US3960449A (en) 1975-06-05 1976-06-01 The Board Of Trustees Of Leland Stanford Junior University Measurement of angular dependence of scattered light in a flowing stream
US4347935A (en) 1979-05-16 1982-09-07 The United States Of America As Represented By The United States Department Of Energy Method and apparatus for electrostatically sorting biological cells
US4667830A (en) 1981-06-15 1987-05-26 The Board Of Trustees Of The Leland Stanford Junior University Method and means for sorting individual particles into containers for culturing, cloning, analysis, or the like
US4545677A (en) 1984-03-05 1985-10-08 Becton, Dickinson And Company Prismatic beam expander for light beam shaping in a flow cytometry apparatus
US4770992A (en) 1985-11-27 1988-09-13 Den Engh Gerrit J Van Detection of specific DNA sequences by flow cytometry
DE3610396A1 (de) 1986-03-27 1987-10-01 Wella Ag Mittel und verfahren zur oxidativen faerbung von haaren
US4704891A (en) 1986-08-29 1987-11-10 Becton, Dickinson And Company Method and materials for calibrating flow cytometers and other analysis instruments
US5245318A (en) 1987-07-24 1993-09-14 Canon Kabushiki Kaisha Particle analyzing apparatus having pressure control system
US5040890A (en) 1987-11-25 1991-08-20 Becton, Dickinson And Company Sheathed particle flow controlled by differential pressure
US5030002A (en) 1989-08-11 1991-07-09 Becton, Dickinson And Company Method and apparatus for sorting particles with a moving catcher tube
JPH0810306B2 (ja) 1990-01-26 1996-01-31 大日本スクリーン製造株式会社 光ビームの偏向制御方法
DE69132273D1 (de) 1990-10-10 2000-08-03 Univ Maryland Verfahren und vorrichtung zur bestimmung von lebensdauer von fluoreszenz in der durchflusscytometrie
EP0568596A1 (en) 1991-01-24 1993-11-10 The University Of Maryland Method and apparatus for multi-dimensional phase fluorescence lifetime imaging
EP0515211A3 (en) 1991-05-23 1993-04-07 Becton Dickinson And Company Apparatus and method for phase resolved fluorescence lifetimes of independent and varying amplitude pulses
WO1993005478A1 (en) 1991-08-28 1993-03-18 Becton, Dickinson & Company Gravitational attractor engine for adaptively autoclustering n-dimensional data streams
US5309912A (en) 1991-11-08 1994-05-10 The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services Multidimensional imaging using a single point detector for a phase encoded modulated optical carrier
US5192870A (en) 1992-01-14 1993-03-09 International Business Machines Corporation Optical submicron aerosol particle detector
US5255257A (en) 1992-03-04 1993-10-19 Lasertape Systems, Inc. Frequency, phase and amplitude control apparatus and method for acousto-optic deflector optimization
US6057814A (en) 1993-05-24 2000-05-02 Display Science, Inc. Electrostatic video display drive circuitry and displays incorporating same
US5296911A (en) 1992-07-16 1994-03-22 Schiapparelli Biosystems, Inc. Optical test system for a chemical analyzer
US5270548A (en) 1992-07-31 1993-12-14 The United States Of America As Represented By The United States Department Of Energy Phase-sensitive flow cytometer
US5293213A (en) 1992-08-12 1994-03-08 Klein Uwe K A Utilization of a modulated laser beam in heterodyne interferometry
US5464581A (en) 1993-08-02 1995-11-07 The Regents Of The University Of California Flow cytometer
US5483469A (en) 1993-08-02 1996-01-09 The Regents Of The University Of California Multiple sort flow cytometer
US5489977A (en) 1993-08-11 1996-02-06 Texaco Inc. Photomeric means for monitoring solids and fluorescent material in waste water using a falling stream water sampler
US5700692A (en) 1994-09-27 1997-12-23 Becton Dickinson And Company Flow sorter with video-regulated droplet spacing
US5602039A (en) 1994-10-14 1997-02-11 The University Of Washington Flow cytometer jet monitor system
US5643796A (en) 1994-10-14 1997-07-01 University Of Washington System for sensing droplet formation time delay in a flow cytometer
US5620842A (en) 1995-03-29 1997-04-15 Becton Dickinson And Company Determination of the number of fluorescent molecules on calibration beads for flow cytometry
JPH0961858A (ja) 1995-08-28 1997-03-07 Minolta Co Ltd 音響光学光走査装置
US5968738A (en) 1995-12-06 1999-10-19 The Board Of Trustees Of The Leland Stanford Junior University Two-reporter FACS analysis of mammalian cells using green fluorescent proteins
JPH10148778A (ja) 1996-11-18 1998-06-02 Minolta Co Ltd マルチビーム発生装置
DE19721882C2 (de) 1997-05-26 1999-04-29 Bosch Gmbh Robert Interferometrische Meßvorrichtung
DE19721881C2 (de) 1997-05-26 1999-05-20 Bosch Gmbh Robert Interferometrische Meßvorrichtung
US6016196A (en) 1997-06-17 2000-01-18 Massachusetts Institute Of Technology Multiple beam pair optical imaging
US6236454B1 (en) 1997-12-15 2001-05-22 Applied Materials, Inc. Multiple beam scanner for an inspection system
US6821740B2 (en) 1998-02-25 2004-11-23 Becton, Dickinson And Company Flow cytometric methods for the concurrent detection of discrete functional conformations of PRB in single cells
US6642018B1 (en) 1998-03-27 2003-11-04 Oncosis Llc Method for inducing a response in one or more targeted cells
AU3897999A (en) 1998-05-14 1999-11-29 Luminex Corporation Multi-analyte diagnostic system and computer implemented process for same
US6271924B1 (en) 1998-12-29 2001-08-07 Bryan Kok Ann Ngoi Noncontact acoustic optic scanning laser vibrometer for determining the difference between an object and a reference surface
US6396069B1 (en) 1999-06-25 2002-05-28 Macpherson David C. Topographer for real time ablation feedback having synthetic wavelength generators
US6372506B1 (en) 1999-07-02 2002-04-16 Becton, Dickinson And Company Apparatus and method for verifying drop delay in a flow cytometer
US6813017B1 (en) 1999-10-20 2004-11-02 Becton, Dickinson And Company Apparatus and method employing incoherent light emitting semiconductor devices as particle detection light sources in a flow cytometer
US6295168B1 (en) 1999-12-15 2001-09-25 International Business Machines Corporation Refractive optical system that converts a laser beam to a collimated flat-top beam
US6809804B1 (en) 2000-05-11 2004-10-26 Becton, Dickinson And Company System and method for providing improved event reading and data processing capabilities in a flow cytometer
US7630063B2 (en) 2000-08-02 2009-12-08 Honeywell International Inc. Miniaturized cytometer for detecting multiple species in a sample
JP2002296178A (ja) 2001-03-30 2002-10-09 Shimadzu Corp フローセル検出器
US6683314B2 (en) 2001-08-28 2004-01-27 Becton, Dickinson And Company Fluorescence detection instrument with reflective transfer legs for color decimation
FR2830629B1 (fr) 2001-10-04 2004-01-09 Commissariat Energie Atomique Procede et dispositif de prelevement d'une partie d'un faisceau lumineux, notamment pour appareil d'analyse de fluorescence
DE20216583U1 (de) 2001-12-20 2003-01-23 Leica Microsystems Heidelberg Gmbh, 68165 Mannheim Mikroskop und Durchflusszytometer
US20060014212A1 (en) 2002-05-10 2006-01-19 Epitome Biosystems, Inc. Proteome epitope tags and methods of use thereof in protein modification analysis
DE10225838B4 (de) 2002-06-11 2021-10-21 Leica Microsystems Cms Gmbh Verfahren zur Scanmikroskopie, Scanmikroskop und Vorrichtung zum Codieren eines Beleuchtungslichtstrahles
US7201875B2 (en) 2002-09-27 2007-04-10 Becton Dickinson And Company Fixed mounted sorting cuvette with user replaceable nozzle
US7745221B2 (en) * 2003-08-28 2010-06-29 Celula, Inc. Methods and apparatus for sorting cells using an optical switch in a microfluidic channel network
US7803624B2 (en) * 2003-09-30 2010-09-28 Cytyc Corporation Automated cytological sample classification
US7590991B2 (en) 2003-10-09 2009-09-15 Terayon Communication Systems, Inc. Method and apparatus for determining channel to which a TV or VCR is tuned
DE10357584B4 (de) 2003-12-08 2006-06-14 Leica Microsystems Cms Gmbh Verfahren zum Trennen unterschiedlicher Emissionswellenlängen in einem Scanmikroskop
US7133187B2 (en) 2004-06-07 2006-11-07 Electro Scientific Industries, Inc. AOM modulation techniques employing plurality of transducers to improve laser system performance
US20070121196A1 (en) 2005-09-29 2007-05-31 The General Hospital Corporation Method and apparatus for method for viewing and analyzing of one or more biological samples with progressively increasing resolutions
US7889348B2 (en) 2005-10-14 2011-02-15 The General Hospital Corporation Arrangements and methods for facilitating photoluminescence imaging
GB0525072D0 (en) 2005-12-09 2006-01-18 Enigma Diagnostics Ltd Fluorescence-based detection methods and apparatus
EP1991887B1 (en) 2006-02-17 2018-10-17 Regents of the University of Minnesota High field magnetic resonance
US20090323061A1 (en) 2006-02-28 2009-12-31 Lukas Novotny Multi-color hetereodyne interferometric apparatus and method for sizing nanoparticles
JP2007285999A (ja) 2006-04-20 2007-11-01 Furukawa Electric Co Ltd:The 光測定装置
ATE440301T1 (de) 2006-05-29 2009-09-15 Olympus Corp Laserscan-mikroskop und mikroskopisches überwachungsverfahren
JP5133600B2 (ja) 2006-05-29 2013-01-30 オリンパス株式会社 レーザ走査型顕微鏡および顕微鏡観察方法
US7697885B2 (en) 2006-09-15 2010-04-13 Aeroflex High Speed Test Solutions, Inc. Multi-band jammer
US7804594B2 (en) 2006-12-29 2010-09-28 Abbott Laboratories, Inc. Method and apparatus for rapidly counting and identifying biological particles in a flow stream
US7400457B1 (en) 2007-01-04 2008-07-15 Stockeryale Canada Inc. Rectangular flat-top beam shaper
US7738094B2 (en) 2007-01-26 2010-06-15 Becton, Dickinson And Company Method, system, and compositions for cell counting and analysis
US8101426B2 (en) 2007-03-02 2012-01-24 Icyt Mission Technology, Inc. System and method for the measurement of multiple fluorescence emissions in a flow cytometry system
JP5619344B2 (ja) 2007-05-04 2014-11-05 マツクス−プランク−ゲゼルシャフト ツール フエルデルング デル ヴイツセンシャフテン エー フアウ モノリシックマイクロ共振器を使用した光周波数コム発生のための装置と方法
US7982944B2 (en) * 2007-05-04 2011-07-19 Max-Planck-Gesellschaft Zur Forderung Der Wissenschaften E.V. Method and apparatus for optical frequency comb generation using a monolithic micro-resonator
US20100233676A1 (en) * 2007-06-14 2010-09-16 Kimberly Kelly High affinity fluorochrome binding peptides
US8848199B2 (en) 2007-07-10 2014-09-30 Massachusetts Institute Of Technology Tomographic phase microscopy
JP4509163B2 (ja) 2007-10-26 2010-07-21 ソニー株式会社 微小粒子の測定方法
US8159670B2 (en) 2007-11-05 2012-04-17 Abbott Laboratories Method and apparatus for rapidly counting and identifying biological particles in a flow stream
GB0800333D0 (en) 2008-01-09 2008-02-20 Ucl Business Plc Beam deflection apparatus and methods
JP4365439B2 (ja) 2008-03-04 2009-11-18 三井造船株式会社 蛍光検出方法及び蛍光検出装置
US8565861B2 (en) 2008-05-02 2013-10-22 Olympus Corporation Optical inspection device, electromagnetic wave detection method, electromagnetic wave detection device, organism observation method, microscope, endoscope, and optical tomographic image generation device
US8140300B2 (en) 2008-05-15 2012-03-20 Becton, Dickinson And Company High throughput flow cytometer operation with data quality assessment and control
US8184279B2 (en) 2008-06-16 2012-05-22 The Regents Of The University Of Colorado, A Body Corporate Fourier domain sensing
EP2333526A4 (en) 2008-09-19 2012-04-04 Mitsui Shipbuilding Eng FLUORESCENCE DETECTION DEVICE USING INTENSITY MODULATED LASER LIGHT AND FLUORESCENCE DETECTION METHOD
US8440952B2 (en) * 2008-11-18 2013-05-14 The Regents Of The University Of California Methods for optical amplified imaging using a two-dimensional spectral brush
US8233146B2 (en) 2009-01-13 2012-07-31 Becton, Dickinson And Company Cuvette for flow-type particle analyzer
WO2010084720A1 (ja) * 2009-01-23 2010-07-29 三井造船株式会社 蛍光検出装置及び蛍光検出方法
TWI594828B (zh) 2009-05-28 2017-08-11 伊雷克托科學工業股份有限公司 應用於雷射處理工件中的特徵的聲光偏轉器及相關雷射處理方法
EP2267430A1 (de) * 2009-06-24 2010-12-29 Masterrind GmbH Vorrichtung und Verfahren zur Selektion von Partikeln
WO2011023593A1 (en) 2009-08-24 2011-03-03 INSERM (Institut National de la Santé et de la Recherche Médicale) Method of and apparatus for imaging a cellular sample
GB0919854D0 (en) 2009-11-12 2009-12-30 Stfc Science & Technology Detecting species in a dilute medium
JP4918178B2 (ja) * 2010-01-21 2012-04-18 三井造船株式会社 蛍光検出方法
JP5649828B2 (ja) 2010-02-03 2015-01-07 オリンパス株式会社 レーザ顕微鏡装置
JP2011163787A (ja) 2010-02-05 2011-08-25 Sony Corp 微小粒子分析装置及び微小粒子分析方法
WO2011097032A1 (en) 2010-02-05 2011-08-11 Cytonome/St, Llc Multiple flow channel particle analysis system
JP5551477B2 (ja) 2010-03-15 2014-07-16 オリンパス株式会社 光源装置およびレーザ走査型顕微鏡装置
JP2012008055A (ja) * 2010-06-25 2012-01-12 Olympus Corp 画像解析方法および画像解析装置
US8384045B2 (en) 2010-07-01 2013-02-26 Sony Corporation Minute particle analyzing device and method
US8528427B2 (en) 2010-10-29 2013-09-10 Becton, Dickinson And Company Dual feedback vacuum fluidics for a flow-type particle analyzer
DE102010044013A1 (de) 2010-11-16 2012-05-16 Carl Zeiss Microimaging Gmbh Tiefenauflösungsgesteigerte Mikroskopie
JP5930220B2 (ja) 2011-03-18 2016-06-08 国立研究開発法人理化学研究所 非線形光学顕微鏡および非線形光学顕微鏡法
CN103649294B (zh) 2011-04-29 2016-05-18 贝克顿·迪金森公司 多路分类***和方法
WO2012158727A2 (en) 2011-05-16 2012-11-22 Oewaves, Inc. Generation of single optical tone, rf oscillation signal and optical comb in a triple-oscillator device based on nonlinear optical resonator
JP6035329B2 (ja) * 2011-05-26 2016-11-30 ザ ジェネラル ホスピタル コーポレイション 光学的トロンボエラストグラフィーシステム及び血液凝固基準の評価方法
US9095494B2 (en) 2011-09-30 2015-08-04 Becton, Dickinson And Company Fluid exchange methods and devices
KR101260051B1 (ko) 2012-01-30 2013-05-06 한국화학연구원 살아있는 세포에 대한 명시야 이미징 및 형광 이미징의 동시 수행이 가능한 세포 이미징 장치 및 방법
US9395425B2 (en) * 2012-08-24 2016-07-19 The Trustees Of Dartmouth College Method and apparatus for magnetic susceptibility tomography, magnetoencephalography, and taggant or contrast agent detection
WO2014152048A2 (en) 2013-03-14 2014-09-25 Cytonome/St, Llc Assemblies and methods for reducing optical crosstalk in particle processing systems
WO2014169231A1 (en) 2013-04-12 2014-10-16 Norton Pierce O Automated set-up for cell sorting
CN103762487B (zh) 2014-01-04 2016-05-25 天津奇谱光电技术有限公司 一种具有双输出光束的可调谐激光器
EP3120130B1 (en) 2014-03-18 2023-07-26 The Regents of the University of California Parallel flow cytometer using radiofrequency mulitplexing, and method
JP6652559B2 (ja) 2014-11-12 2020-02-26 オーボテック リミテッド 光学装置及び光学方法
CN108474728B (zh) 2015-10-13 2022-04-26 贝克顿·迪金森公司 多模态荧光成像流式细胞术***
EP3411713B1 (en) * 2016-02-05 2021-06-30 NanoView Biosciences, Inc. Detection of exosomes having surface markers
KR102641469B1 (ko) 2016-03-17 2024-02-28 벡톤 디킨슨 앤드 컴퍼니 고효율 형광 유세포 분석기를 사용하는 세포 선별
EP4242631A3 (en) 2016-09-13 2024-03-06 Becton, Dickinson and Company Flow cytometer with optical equalization

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030031352A1 (en) * 2001-08-10 2003-02-13 Nelson Alan C. Optical projection imaging system and method for automatically detecting cells with molecular marker compartmentalization associated with malignancy and disease
US20050207633A1 (en) * 2003-04-02 2005-09-22 Nick Arini Method of, and computer software for, classification of cells into subpopulations
CN102792146A (zh) * 2010-02-10 2012-11-21 赫拉巴Abx公司 用于流体中微颗粒的多参数测量的装置和方法
US20110275558A1 (en) * 2010-05-04 2011-11-10 Virginia Tech Intellectual Properties, Inc. Lanthionine synthetase component c-like proteins as molecular targets for preventing and treating diseases and disorders
US20120128264A1 (en) * 2010-11-23 2012-05-24 General Electric Company Methods and systems of optical imaging for target detection in a scattering medium
US20130078625A1 (en) * 2011-09-25 2013-03-28 Theranos, Inc., a Delaware Corporation Fluid handling apparatus and configurations
CN102650527A (zh) * 2012-05-25 2012-08-29 北京航空航天大学 一种基于时间序列分析消噪的光纤陀螺温度补偿方法
US20130323825A1 (en) * 2012-05-29 2013-12-05 Sony Corporation Information processing apparatus, information processing method, and program
US20160003741A1 (en) * 2013-01-09 2016-01-07 The Regents Of The University Of California Apparatus and methods for fluorescence imaging using radiofrequency-multiplexed excitation
JP2016504598A (ja) * 2013-01-09 2016-02-12 ザ リージェンツ オブ ザ ユニバーシティ オブ カリフォルニアＴｈｅ Ｒｅｇｅｎｔｓ Ｏｆ Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ Ｏｆ Ｃａｌｉｆｏｒｎｉａ 無線周波数多重励起を用いた蛍光イメージングのための装置及び方法
US20160054293A1 (en) * 2013-03-22 2016-02-25 Map Diagnostics Limited Prenatal Screening
US20150182136A1 (en) * 2013-12-26 2015-07-02 Fundacio Institut De Ciencies Fotoniques Speckle contrast optical tomography
WO2016054293A1 (en) * 2014-09-30 2016-04-07 The Regents Of The University Of California Imaging flow cytometer using spatial-temporal transformation

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112037217A (zh) * 2020-09-09 2020-12-04 南京诺源医疗器械有限公司 一种基于荧光成像的术中血流成像方法
CN114088606A (zh) * 2021-10-23 2022-02-25 广州市艾贝泰生物科技有限公司 细胞分析装置

Also Published As

Publication number Publication date
US20210131945A1 (en) 2021-05-06
JP2019521316A (ja) 2019-07-25
JP7023244B2 (ja) 2022-02-21
WO2017197271A1 (en) 2017-11-16
KR20190006164A (ko) 2019-01-17
US20170328826A1 (en) 2017-11-16
US10935485B2 (en) 2021-03-02
CN109477784B (zh) 2022-04-01
US20230349812A1 (en) 2023-11-02
KR102381176B1 (ko) 2022-04-01
EP3455608A1 (en) 2019-03-20

Similar Documents

Publication Publication Date Title
CN109477784A (zh) 具有增强图像分辨率的荧光成像流式细胞术
US10935482B2 (en) Multi-modal fluorescence imaging flow cytometry system
US20210356381A1 (en) Cell Sorting Using A High Throughput Fluorescence Flow Cytometer
US20230228668A1 (en) Spectral Unmixing of Fluorescence Imaging Using Radiofrequency-Multiplexed Excitation Data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant