CN103180858A - 移动微实体流中的运动相关误差的补偿 - Google Patents

移动微实体流中的运动相关误差的补偿 Download PDF

Info

Publication number
CN103180858A
CN103180858A CN2012800034774A CN201280003477A CN103180858A CN 103180858 A CN103180858 A CN 103180858A CN 2012800034774 A CN2012800034774 A CN 2012800034774A CN 201280003477 A CN201280003477 A CN 201280003477A CN 103180858 A CN103180858 A CN 103180858A
Authority
CN
China
Prior art keywords
micro
entity
signal
stream
value
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
CN2012800034774A
Other languages
English (en)
Other versions
CN103180858B (zh
Inventor
G·都拉克
S·D·弗雷施尔
J·哈特彻尔
D·罗伯特
M·卒尔丹
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.)
Sony Corp
Sony Corp of America
Sony Pictures Entertainment Inc
Original Assignee
Sony Corp
Sony Pictures Entertainment 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 Sony Corp, Sony Pictures Entertainment Inc filed Critical Sony Corp
Publication of CN103180858A publication Critical patent/CN103180858A/zh
Application granted granted Critical
Publication of CN103180858B publication Critical patent/CN103180858B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10—Investigating individual particles
    • G01N15/14—Optical investigation techniques, e.g. flow cytometry
    • G01N15/1404—Handling flow, e.g. hydrodynamic focusing
    • 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/1425—Optical investigation techniques, e.g. flow cytometry using an analyser being characterised by its control arrangement
    • G01N15/1427—Optical investigation techniques, e.g. flow cytometry using an analyser being characterised by its control arrangement with the synchronisation of components, a time gate for operation of components, or suppression of particle coincidences
    • 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/1456—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals
    • G01N15/1459—Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals 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
    • 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/149—Optical investigation techniques, e.g. flow cytometry specially adapted for sorting particles, e.g. by their size or optical properties
    • 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
    • G01N2015/1027—Determining speed or velocity of a particle
    • 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/1404—Handling flow, e.g. hydrodynamic focusing
    • G01N2015/1415—Control of particle position
    • 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/1438—Using two lasers in succession
    • 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
    • G01N2015/1477—Multiparameters

Landscapes

  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明描述了用于检测、表征和补偿移动微实体的运动相关误差的设备和方法。运动相关误差可出现在移动微实体的流中，并且可表示期望到达时间的偏差或流中微实体位置的不确定性。微实体的运动相关误差（例如脉冲抖动）在流式细胞计中观测到，并且被发现具有与系统的参数的函数相关性。根据一个实施例，可通过调节数据获取观测窗来补偿脉冲抖动。对于流式细胞计，脉冲抖动的减小可以改善测量准确性、双液滴的分辨率、系统处理量以及能增大用于探测微实体的询问区域。

Description

移动微实体流中的运动相关误差的补偿
相关申请的交叉引用
本申请根据35U.S.C.119(e)要求于2011年8月25日提交的名为“System and Method for Correction of Pulse Jitter in FlowCytometers”的美国临时申请No.61/527,340的优先权，其整体在此引入作为参考。
技术领域
所述方法和设备涉及对移动微实体流中的运动相关误差进行补偿。对更高阶运动特征的评估和补偿能够改善对移动微实体的测量和/或操作。应用领域包括流式细胞术、微流体（microfluidics）及纳米制造。
背景技术
近年来，对微实体的分析、开发和制造领域已经变得在各种技术中广泛应用。这里使用的微实体可包括细胞、微生物、微颗粒或纳米颗粒、液滴、分子、蛋白质、肽、校准颗粒、以及微制造结构，该微制造结构诸如为微机电（MEM）结构、微电子芯片及微传感器。微实体可包括最大尺寸小于约1毫米（mm）的任何实体（人造或自然生成的）。在一些研究、医学或制造应用中，大数目的微实体可以流的方式移动。所述流可以是控制对进入流的微实体的移动和计量的系统的一部分。一个示例应用为对流式细胞计中的细胞进行表征和分选。另一个实例可以为对微流体通道中的微实体进行移动、观察和分选。
在一些应用中，可能期望在一个或多个选择的时间精确获知移动流中的所选择的微实体的位置，或者在其它实施例中，精确获知所选择的微实体将到达一个或多个位置的时间。例如，参考图1，微实体125可在微实体流120中移动。微实体的流120可包括能够输送微实体的流体流（例如气体、液体或颗粒流）或流动路径。微实体可以悬浮在流中。在一些情况中，微实体的流120可包括以类似流的方式移动的微实体的稠密集合。微实体的流120可通过机械、电学和磁方式的任何组合传送。
可在第一时间t1于第一位置P1检测到微实体125，并且可能期望知道微实体125到达沿流的预先选择的第二位置P2或多个其它位置的时间t2。相反地，可在第一时间t1于第一位置P1检测到微实体125，并且可能期望知道微实体125在预先选择的经过时间t3或其它多个时间之后将位于哪个位置P3或其它多个位置。在前者的情况中，预测到达时间t2有利于于确定何时可以在第二位置P2处对微实体进行测量或操作（例如，分选操作、成像操作、充电操作、转化操作、照射实体、信号检测等）。在后者的情况中，预测微实体125的位置P3有利于精确跟踪流中的实体的移动和/或演变，例如，当微实体沿流移动时获取并重叠微实体的多个图像。
一种用于预测移动流中的微实体的到达时间t2或位置P3的方法为基于平均流速度vavg计算t2或P3的期望值。平均流速度可以以任何合适的方法测量或确定。如果求出或已知vavg，则可以通过关系式d=vavg×t确定t2或P3，其中d为微实体125在经过时间t中行进的距离。一种用于预测在预定位置P2的到达时间t2的可选方法为，在P1处观测到实体之后将平均传输时间或延迟时间Δtavg加到t1的值。时间Δtavg可以以任何合适的方法测量或确定。在一些情况中，这些预测移动微实体的到达时间或位置的方法已经足够。
流式细胞计是使用流以传输微实体用于对诸如生物细胞的微实体进行表征和分选的系统实例。流式细胞计广泛用于快速分析异质细胞悬浮物，以识别构成子种群。其中使用流式细胞计细胞分选的许多应用实例包括：隔离稀少种群的免疫系统细胞用于艾滋病（AIDS）研究，隔离遗传非典型细胞用于癌症研究，隔离特定染色体用于基因研究，以及隔离各种微生物品种用于环境研究。例如，荧光标记的单克隆抗体被经常用作“标记物”以识别诸如T淋巴细胞和B淋巴细胞的免疫细胞，临床实验室使用该技术以对HIV感染病人的“CD4阳性”T淋巴细胞数目进行计数，他们还使用该技术以识别与各种白血病和淋巴癌相关的细胞。
近来，相比于严格地研究应用，两种关注领域为面向临床的移动细胞分选，病患护理应用。首先从化学制药开发转移到生物制药开发。例如，许多新的癌症治疗使用生物材料。这些包括一类基于抗体的癌症疗法。细胞分选器在识别、开发、提纯和最终制造这些产品中起到重要的作用。
与此相关的是转移到使用细胞置换疗法用于病患护理。目前在干细胞上的很多关注围绕经常称为再生疗法或再生医学的新医学领域。这些疗法可能经常要求将相对稀少的细胞从病人的组织分离。例如，可从骨髓分离出成人干细胞并最终作为再注入物的一部分回到其被移出的病人。流式细胞术和细胞分选是使得可以进行这些治疗的重要组织处理工具。
发明内容
描述了用于对在流中移动的微实体的运动相关误差进行补偿的设备和方法。移动微实体的运动相关误差可包括任何类型和形式的误差，该误差造成在已知在第一位置或第一时间检测到微实体的情况下相同微实体在第二位置或第二时间的期望到达时间或期望位置的偏移。移动微实体的运动相关误差可包括实际物理移动和/或微实体的似动（apparent movement）。可以理解，运动相关误差可与沿流的各个点的多个到达时间和位置相关。一个或多个期望到达时间或一个或多个期望位置可基于流中的运动的基本假设，例如，恒定速度和/或恒定加速度。在一些实例中，对运动相关误差的表征可获得说明各种类型的运动相关误差的微实体运动的更精确的模型。对微实体的运动相关误差的表征和补偿有利于诊断和改善被设计用于对流中移动的微实体进行测量和/或操作的设备的性能。
在一个实例中，开发技术和设备以表征和补偿流式细胞计的运动相关误差。发现在流式细胞计流中移动的微实体的询问点到达时间的偏移是稳定的并且函数地依赖于关于系统参数测量的值。可使用偏移的函数相关性的模型以减小测量误差并改善流式细胞计的操作（例如，细胞处理能力和/或分选纯度）。
根据一些实施例，一种用于补偿与在第一位置与至少第二位置之间移动的微实体相关联的运动相关误差的系统，该系统包括：检测装置，被配置为当微实体穿过第一位置时产生第一信号，以及当该微实体穿过至少第二位置时产生至少第二信号。该第一位置可以在流中位于第二位置之前，或者可以在流中位于第二位置之后。该系统可以包括一个或多个处理器，该一个或多个处理器被配置为至少接收所述第一信号，以及根据所述第一信号相对于系统的参数确定一个值。可以针对穿过第一位置的每个微实体确定该值。作为示例，该值可以是相对于系统度量（例如系统时钟或周期信号）测量的微实体到达第一位置的时间偏移或相位偏移。系统的一个或多个处理器还被配置为，按照补偿与微实体相关联的运动相关误差的校正量来调节在第二位置处对微实体进行观测或操作的观测时间或操作时间。可以根据误差模型确定校正量，该误差模型根据相对于系统参数确定的值来预测与微实体相关联的运动相关误差。可以确定校正量并使用该校正量来调节在沿着微实体移动通过的路径上的其它位置的测量或观测。
在一些实施例中，用于补偿与移动微实体相关联的运动相关误差的系统包括周期性能量源（例如换能器），该周期性能量源被配置为把周期性能量耦合到移动微实体的流中。相对于系统的参数确定的值可以包括相对于提供给周期性能量源的驱动信号的特征测量的时间偏移或相位偏移。该系统还包括激光或聚焦辐射源，用于在沿着流的一个或多个位置处光学地探测移动微实体。系统还可以包括光学检测器，被配置为检测来自被光学探测的微实体的辐射。
还提出了用于补偿在流中移动的微实体的运动相关误差的相关方法。例如，根据一个实施例，一种补偿与在系统中的第一位置和至少第二位置之间移动的微实体相关联的运动相关误差的方法可以包括：利用检测装置检测表示第一微实体出现在第一位置处的信号。该方法还可以包括：利用至少一个处理器，根据针对所述第一微实体的所述信号，相对于系统的参数确定一个值。另外，该方法可以包括：利用所述至少一个处理器，按照补偿与所述第一微实体相关联的运动相关误差的校正量来调节用于在第二位置处对所述第一微实体进行观测或操作的观测时间或操作时间。根据误差模型确定所述校正量，所述误差模型根据相对于系统的参数确定的值来预测与第一微实体相关联的运动相关误差。
在本发明范围内还提出了一种有形的制造的计算机可读介质和/或存储装置，其包括可机读指令，该指令当由至少一个处理器执行时调节至少一个处理器以执行用于对在流中移动的微实体的运动相关误差进行表征和/或补偿的方法的一个或多个步骤。计算机可读介质和/或存储装置可用于调节被设计为在流中移动的微实体上测量和/或操作的仪器。
通过下面结合附图的描述，可以更全面地理解本发明前述和其它方面、实施例和特征。
附图说明
本领域技术人员可以理解，这里所述附图仅为了说明的目的。可以理解，在一些情况中，为了方便理解本发明，可能扩大或放大地示出本发明各方面。在附图中，在各个图中，相同的标号通常表示相同的特征、功能相似和/或结构相似的元件。附图不一定成比例，重点是为了说明本发明的原理。附图不会以任何方式限制本发明范围。
图1示出移动微实体的流；
图2示意示出根据一个实施例的流式细胞计的选定部件；
图3A示意示出用于说明液滴动力学的校准的流式细胞计的选定部件；
图3B-3G示出根据一些实施例的可用于指示微实体在相应询问点的存在或表征微实体的信号图；
图4A-4F示出根据一些实施例的可在流式细胞计中使用或观察到的各种系统信号和数据；
图5为在流式细胞计中观察到的脉冲抖动的一个实例；
图6A示出根据一个实施例的当关于流式细胞计中的参考信号测量时在第一位置检测到的来自微实体的触发信号的延迟；
图6B示出图6A中观察的微实体在第二位置处的到达时间的对应变化；
图7A-7C为对于流式细胞计测量的抖动图的实例；
图8A和8B示出了根据一个实施例的脉冲抖动补偿；
图9示出用于对运动相关误差进行表征和补偿的方法的一个实施例；
图10A示出根据一个实施例的对记录的抖动图的正弦拟合（黑线）；
图10B表示根据正弦拟合补偿的图10A的抖动数据；
图11A是图10A的抖动数据的直方图；
图11B是图10B的经过补偿的抖动数据的直方图。
通过下面结合附图的详细描述，本发明的特征和优点将变得更加清楚。
具体实施方式
如上所述，在一些涉及移动微实体的应用中，可能期望预测或获知微实体到达预选位置的到达时间，或相反地预测或获知微实体在选定时间的位置。对于低精度和/或低速系统，基于基本假设（例如，流中的实体的平均速度vavg，或平均传输时间Δtavg、恒定速度、恒定加速度）的预测对于获得微实体到达时间或其位置的合理精确的预测以对微实体进行测量或操作已经足够。然而，对于高精度和/或高速系统，系统中的微实体的更高阶运动特征可能导致流中的实体的到达时间或位置的不能接受的大的不确定性或误差。大的不确定性可能限制对移动微实体的操作和/或采集代表移动微实体的数据的能力。从而，时间或位置的大的不确定性可能对速度、精度、及系统的性能施加限制。
根据一些实施例，可以对于在系统中在第一位置和至少第二位置之间移动的微实体的运动相关误差进行表征和补偿。在一些系统中，运动相关误差可表现出与对移动微实体关于系统的参数测量的值的函数的和/或稳定的相关性。可进行对移动微实体的测量以表征运动相关误差并揭示误差与系统参数的相关性。通过这样的表征测量，可构建运动相关误差的模型并使用该模型来预测和/或补偿在系统中移动的微实体的运动相关误差。
图2示意示出流式细胞计200的某些元件。流式细胞计是包括移动微实体225（通常为在单个纵列中移动的细胞）的流的系统的一个实例。尽管将主要结合流式细胞计描述本发明的实施例，相同的技术和设备可用于使用移动微实体的流的其它系统或设备中。例如，本发明实施例可实施于微流体系统中和微米或纳米制造系统中，其中在流体流中制造微米或纳米实体。各个系统可适于对移动微实体进行一个或多个操作。这些系统中的操作类型可包括下列操作的任何组合：测量微实体的特征、成像、照射以检测荧光性、分选、转换（例如，引入化学或电学电荷、照射以交联聚合物、照射以转换微实体的成分）、尺寸分析、以及重量分析。
I.流式细胞计系统概述
30多年前首次将基于流式细胞术的细胞分选引入研究群体。这是一种已经广泛应用于生命科学研究的许多领域的技术，对于在诸如遗传学、免疫学、分子生物学和环境科学的领域中工作的技术人员是极为有用的工具。与诸如免疫淘选或磁柱分离的批量细胞分离技术不同，基于流式细胞术的细胞分选设备以每秒几千个细胞或更高的速度串行地测量、归类、然后分选各个微实体。该对单个细胞的快速“逐个”处理使得流式细胞术成为用于从来自另外异质细胞悬浮物的细胞提取高纯度子群的有价值的工具。
通常通过荧光材料以一些方式标记将分选的细胞。束缚到细胞的荧光探针在该细胞通过紧密聚焦、高强度光束（通常为激光束）时发出荧光。计算机记录每个细胞的发射强度。这些数据然后被用于对每个细胞分类，以用于具体分选操作。基于流式细胞术的细胞分选已经成功应用于多种细胞类型、细胞组分和微生物，以及多种尺寸相当的无机颗粒。
在当前的使用中有两种基本类型的细胞分选器。它们是“液滴细胞分选器”和“流体切换细胞分选器”。液滴细胞分选器使用微液滴作为将选定细胞传输或传送到采集器皿的容器。通过将超声能量耦合到喷射流而形成微液滴。包含选择的细胞的液滴然后被静电导引到希望的位置以对来自较大群的选择细胞进行分选。这可以是非常有效的过程，可以从单个流每秒分选多达90,000个细胞。
第二种基于流式细胞术的细胞分选器是流体切换细胞分选器。一些流体切换细胞分选器使用压电装置以驱动机械系统，该机械系统将流动样本流的一段偏转到采集器皿中。相比于液滴细胞分选器，流体切换细胞分选器由于用于偏转样本流的机械系统的循环时间可能具有较低的最大细胞分选速率。该循环时间，即初始样本偏转与恢复稳定非分选流动的时间之间的时间，通常显著大于在液滴细胞分选器上的液滴生成器的周期。该较长的循环时间将一些流体切换细胞分选器限定为几百个细胞每秒的处理速率。由于相同的原因，通过流体细胞分选器切换的一段流通常是来自液滴生成器的单个微液滴的体积的至少十倍。这导致相比于液滴分选器的采集器皿，在流体切换分选器的采集器皿中的细胞浓度相应较低。
近来，微流体分选器已经发展为根据需要使用MEMS致动器或电场以偏转微实体。这些分选器可以在诸如微流体芯片的完全封装的芯片中分选。使用微流体系统使得可以从相同的输入通道进行多个并行分选，从而增大了微流体分选器中的净处理量。
参考图2，流式细胞计200可包括喷嘴202，从喷嘴流出流体流220。喷嘴202可包括容纳鞘液（sheath fluid）的第一室205。喷嘴还可以包括导管203，该导管203将样本流体传送到流220的芯部。微实体225可悬浮在样本流体中。鞘液和样本流体可朝向喷嘴孔208流动，以流220的形式从喷嘴202流出。鞘液可形成流外围的鞘层。
喷嘴202还可以包括换能器212，该换能器212被配置为将能量耦合到至少鞘液中。换能器212可以是声学换能器，其基于规则的周期将声能耦合到鞘液中。将周期性能量耦合到鞘液和/或样本流体中可以高一致性调节从喷嘴202流出的流220中的液滴228的形成。由耦合到喷嘴内的流体的声能（例如，来自喷嘴中的压电换能器的周期性振动）导致流220中的扰动。由于如通过瑞利预测的表面张力，在声学频率下形成液滴228。在可接受的参数范围内（频率、鞘层压力、换能器驱动幅度等），液滴在距离喷嘴一致的距离处分离出来。在一些实施例中，细胞分选器可以实现一种反馈系统，该反馈系统把分离位置226保持在距离喷嘴孔208近似恒定的距离处。还可以与液滴228一起形成卫星液滴227。
如图2所示，从液滴细胞分选器的喷嘴流出的射流沿着流220向下承载微实体225行进了几个毫米，直到微实体255被封入从流同步分离出来的微液滴228之一中。在一些实施例中，微实体可以沿流以近似恒定的速度移动。在一些液滴细胞分选器中，从测量点（例如在束B1的第一询问点）到细胞通过分离点226（流分离为自由液滴的位置）的位置的飞行时间（tf）恒定在液滴生成时段的约1％内。
液滴228的形成可受到施加到换能器212的信号的水滴驱动频率fd和振幅Ad的影响。在一些实施例中，可通过水滴时钟265调节水滴驱动频率fd，水滴时钟265又可通过系统的处理器250控制。可通过被配置为向换能器212提供驱动信号的驱动器260调节水滴驱动振幅Ad。在一些实施例中，驱动器260可从处理器250接收振幅控制信号。根据一些实施例，增加水滴驱动振幅Ad将使液滴分离点226向喷嘴移动。改变驱动频率fd改变形成液滴的频率。在一些实施例中，驱动频率可选为具有在约10kHz和约250kHz之间的值。
处理器250可以是任何类型和形式的数据处理装置，例如，微处理器、微控制器、连接到系统的计算机、或现场可编程门阵列（FPGA）。在一些实施例中在系统中可有多于一个的处理器。在一些情况中，可以存在多种处理器类型的组合，例如微处理器和一个或多个FPGA。
系统200还可以包括用于在每个液滴228上放置电荷的设备。对液滴充电可用于各个液滴的分选操作，从而分选液滴中包含的微实体。可通过处理器250控制电荷源270，以在选定的时间量中对与鞘液和/或样本流体电接触的导电探针207施加选定的电势值。例如，当在流220中形成最后附着的液滴228-0时，可对导电探针207施加第一电势。由于流体是导电的，流体的表面带电，并且最后附着的液滴获取电荷。当最后附着的液滴从流分离以形成液滴时，例如移动到前一个液滴228-1的位置，新形成的液滴保存的电荷与在最后附着的液滴228-0从流220分离时施加到流体和流的电势相关。在一些情况中，在形成的液滴上的电荷与施加到流体的电势的方根成比例。在最后附着的液滴分离后，可以对导电探针207施加第二电势，从而形成为最后附着的液滴的下一个液滴将获得与其前任不同或相同的电荷。在一些情况中，可以不对液滴施加电荷（例如，接地电势或中性电荷）。
系统200还可以包括用于分选液滴的设备，从而对通过液滴传送的任何微实体进行分选。分选设备可包括静电偏转板290，对该静电偏转板施加电势以在板之间建立电场。当带电的液滴横穿过电场时，该液滴可以被横向偏转并被采集到样本器皿（未示出）中。被提供中性电荷的液滴可以在偏转板290之间通过而不被偏转板偏转。以该方式或类似方式，可以对微实体进行分选。
在各个实施例中，系统200可包括用于检测操作的部件，例如检测可用于表征移动微实体的一个或多个信号。该一个或多个信号可用于确定在含微实体的流分段上进行的后续操作（例如，分选操作）。对于图2所示的实施例，检测操作可包括通过来自一个或多个辐射源S1、S2、...S5的光学辐射照射微实体。辐射源可发出不同波长的辐射束B1、B2、...B5。在一些实施例中，一个或多个源可发出相同波长的辐射。在一些实施例中的辐射源S1、S2、...S5可以是激光源、在一些实施例中可以是高强度弧光灯源、在一些实施例中可以是高强度发光二极管、或者任何合适的辐射源。辐射束可以是分离的，并被设置为在喷嘴202附近与流220相交。在一些实施例中，可将辐射束聚焦到流220上。辐射束B1、B2、...B5可与微实体相互作用以产生散射的辐射或激发荧光辐射。例如可通过检测器D1检测散射或荧光辐射，检测器D1可生成代表检测到的辐射的电信号。可通过例如一个或多个检测器D2、D3、...D5检测代表荧光辐射的一个或多个信号。可通过光学元件215采集来自微实体的荧光辐射和/或散射的辐射，并将其引导到检测器。
在一些流式细胞计和其它系统中，例如微流体系统，可以使用其它类型的颗粒检测和感测。在一些实施例中，可以对沿流移动的微实体进行直接电感测或库尔特型（Coulter-type）测量。在一些情况中，可进行库尔特体积测量以表征微实体。在一些实施中，可使用其它类型的感测，例如磁的、磁光的、电光的、热学的感测。可使用利用这些检测和感测技术中的任意技术的信号或测量结果以表征和/或补偿微实体的运动相关误差。例如，可使用利用这些检测和感测技术中的一个或多个获得的信号或测量结果来开发用于微实体的运动相关误差的模型，并且可以使用随后的信号或测量结果来确定补偿运动相关误差的校正量。
液滴细胞分选器可被配置为作为空气中感测（sense-in-air，SIA）设备操作。在这样的系统中，至少一个激励激光器束被聚焦到流220上，该流220包括流动鞘层流和传送微实体（以单个纵列对准）的内部芯流。该流被定位为，使得该相交在一个或多个光采集系统的焦点处进行。共同的焦点经常被称为询问点。如图2所示，沿流220可以空间分离地具有多个激光系统和多个询问点。这些询问点中的每一个可具有一个或多个分离的激光束和信号检测设备。可将询问点空间限制在喷嘴202附近的询问区域222中。可使用空间过滤以减小与单个询问点相关的观测区域（还称为视场（FOV）），并减小附近询问点的交叉检测。用于单个询问点的FOV可以在50微米的量级。从而，要求流、激光器焦点、和光学系统的焦点的严格和精确的对准。在一些实施例中，与单个询问点相关的观测区域可限于几十微米。当微实体沿着流行进时，其顺序通过每个询问点。可以从每个询问点测量和记录来自微实体的辐射。这些测量结果在从包含微实体的分段的流分离前是相关的。带电的液滴然后可以通过空气行进到静电板290，其在该静电板290中将被偏转以用于分选。
尽管图2中示出了五个探测束和五个检测器，在一些系统中可使用更少或更多的探测束和更少或更多的检测器。在一些实施例中，可以设置多于一个的检测器用于探测束，并且光学部件可以将光发射分割并引导到多于一个的检测器。在一些情况中，可在检测光学器件中使用光学和/或空间过滤，以分离荧光信号和/或阻挡来自辐射源的辐射。使用五个光源以探测微实体的流式细胞计的一个实例为sy3200TM流式细胞计（可从iCyt Mission Technology of Champaign, Illinois获得）。
例如可以通过数据获取硬件230数字采样地采集来自检测器D1、D2、D3...D5的时间变化信号，并将其存储到存储器240中用于后续通过处理器250的分析，所述数据获取硬件230可通过模拟到数字转换时钟245来定时。在一些实施例中，可将数据直接提供到处理器250用于分析。在一些实施例中，将数据临时存储到存储器240中，例如缓冲器或适当的随机存取装置，并然后通过处理器250读取用于分析。
如从上述对流式细胞计的描述可以理解，来自检测器的数据信号可以是代表横穿辐射束B1、B2、...B5的微实体225的脉冲。每个脉冲可具有峰值高度、宽度和面积。所述脉冲高度、宽度和面积可通过处理器250评估并用于表征每个微实体225。然后可使用评估的特征确定对每个微实体225或包含一个或多个微实体的液滴的后续操作。
尽管可以以周期间隔形成液滴228，微实体225可以以随机间隔分散在流220中。例如，微实体可以在供给到喷嘴202的样本流体中处于稀释悬浮状态。微实体然后可以在随机时间到达喷嘴的开口208。从而，如图2所示，一些液滴可能不包含微实体，一些液滴可包含多于一个的微实体，例如两个、三个或更多。
对细胞分选器的正确操作要求将从每个询问点记录的测量数据与生成数据的每个微实体正确地相关。为正确表征流220中的选定微实体，需要获知来自检测器D1、D2、D3、...D5的什么信号对应于选定的微实体。为了在封入液滴时在选定微实体上正确进行操作（例如分选），需要进一步知道微实体将被封入哪个液滴228中。例如，需要将每个检测信号正确分配给从其产生检测信号的同一微实体（以及将包含微实体的液滴），从而可以确定对该微实体的精确表征。另外，为了流式细胞计实例中的分选目的，需要知道微实体何时到达最后附着的液滴228-0，以在正确的时间对导电探针207施加正确的电势，以便给予包含该微实体的液滴228-0希望的电荷。当适当地带电时，最后附着的液滴将随后分离为自由的液滴并根据其电荷被分选。
作为理解实施为流式细胞计的本发明内容的辅助，现在将描述用于空气中分选流式细胞计的数据处理、操作定时、以及校准方法的细节。然而，所述技术和设备可应用于除流式细胞计以外的系统，在这些系统中，微实体从至少第一位置移动到一个或多个第二位置，并且其中可在移动微实体上进行操作。
II校准
可利用校准颗粒在校准试验中建立系统的数据处理和系统操作的定时，在所述系统中，微实体从第一位置移动到一个或多个其它位置。作为一个实例，可校准液滴细胞分选器使得数据获取系统可将对每个细胞或微实体采集的数据与将承载其（并且可能是其它微实体）离开流并进入采集器皿的具体液滴相关。可在制造设备中执行校准试验或运行，以初始设置系统，以及/或者其可以日常地运行（例如在每次实验之前），以验证系统的操作。可利用高荧光的微实体通过经验观测完成该校准。可通过将液滴分选到用于人工观察的显微镜载片上或通过使用自动检测系统以记录来自液滴的发射来观测分选结果。根据一些实施例，用于校准目的的微实体可包括校准珠，例如荧光标记的聚苯乙烯微珠。在一些实施例中，可使用荧光细胞用于校准目的。可以将校准珠稀释到样本流体中至这样的程度，即，例如，如图3A所示，在任意给定时刻从喷嘴202流出的流220中有两个或更少的珠。
当微实体225或校准珠沿流220行进时，其通过位于流中的询问束B1、B2、...B5，并之后由从流分离的液滴228-1封装。对于多询问点系统，当单个微实体顺序通过每个询问点时，需要把对该单个微实体的多个测量结果（例如如图3B-3F所示的信号）进行相关。相关的数据用于对每个微实体进行分选决定，并影响对包含微实体的液滴的分选。
可在校准进程中把来自多个询问点的信号正确识别为与单个微实体和液滴相关。通过使校准微实体运行通过系统，数据获取系统可被配置为使得，单个微实体通过选定“触发”询问点（例如，具有辐射束B1的第一询问点）触发一个定时序列，该定时序列在期望微实体到达每个其它询问点时精确地“打开”或选择窄数据获取窗口（例如，从每个询问点选择数据流的分段）。为了改善系统处理量（对于细胞分选器较重要），应将这些定时窗口的持续时间保持为尽可能地短。
在现代数字数据获取系统中，通常对来自全部询问点的测量数据缓存模拟到数字转换（ADC）样本。在一些情况中，将数据缓存足够长，以允许从下游询问点发生触发，而不是仅根据来自沿流的第一询问点的触发。换句话说，当缓存数据时，触发询问点可以是这样的点，在该点，在已经在其它询问点生成针对微实体的一个或多个信号后生成触发信号。随后，选定的“触发”询问点可以是系统中的询问点中的任意一个。
根据一些实施例，如果将获知流的速度和到流中的每个询问点的距离以及到断开点226的距离，可以从基于第一辐射束B1的通过事件的计算估计辐射束通过时间和在最后附着液滴的到达。然而，在一些系统中，距离随每天的操作或每次运行的操作而变化，并且没有方便的方法以直接测量距离或流体速度。
在一些实施例中，可通过校准进程确定束横穿的有效传输时间和断开点。作为实例，如图3B所示，横穿第一束B1的荧光校准珠可在时间t0生成荧光信号。在一些实施例中，该信号可用作“触发”信号，其向系统表示将分析并随后处理（例如，分选）的新的微实体的到达。在其它实施例中，并且如上所述，可使用另一个询问点来提供触发信号而不是第一询问点。与该触发信号一起并且对应于触发信号，如图3C-3G所示，可通过导致触发信号的相同微实体生成在近似时间t1、t2、t3、t4、tc发生的一组信号或事件。
通过将多个校准珠运行通过系统，可确定在各个信号和/或事件之间发生的平均时间偏移或传输时间dt1、dt2、dt3、dt4、dtc。当这些时间已知时，然后可以确定哪个信号与启动触发事件的微实体相关。例如，在与时间t0的触发事件隔开的时间dt1、dt2、dt3、dt4的其它检测器数据流中出现正确信号。在处理信号之后，还可以确定在将包含该微实体的液滴上将施加什么电荷值。在一些实施例中，可以处理来自在承载所关注的微实体的流的分段中或附近的其它微实体的信号，以确定对从流分段形成的液滴将施加什么电荷。当已知时间偏移时，系统使用被分析和处理的有效样本来运行，而不是校准珠。
在一些实施例中，当同一微实体通过其它询问束B2、B3、B4、B5时，可能生成或不生成其它信号，例如在时间t1、t2、t3、t4的脉冲。如果荧光微实体对询问点处的激发辐射没有响应，例如如图3E所示，在该询问点可以不生成信号。当微实体在时间tc到达最后附着液滴228-0时可以生成信号。根据一些实施例，当微实体到达最后附着液滴228-0以及对最后附着液滴充电时，可以通过频闪观测仪激发来对该微实体成像。
由于可以数字采样和记录来自询问点的信号，可通过数字数据流（例如来自模拟到数字转换器的数据流）中的一组数字样本表示来自询问点的信号。因此，校准过程可识别，在来自询问点的检测设备的数据流中的哪个数据流分段或观测窗对应于在穿过触发询问点时启动触发事件的微实体。校准过程可配置数据获取系统，以从与相同微实体相关的每个询问点识别数据流中的数据分段的定时或位置。在一些实施例中，可以使用与来自每个检测器的生成数据相关的时间戳，以辅助识别属于每个微实体的脉冲。
从上述讨论可以理解，在出现触发脉冲时系统200要解决两个问题。第一个问题涉及识别来自其它检测器数据流的哪些信号是由生成触发脉冲的相同微实体生成的。如上所述，可通过从校准进程导出的时间偏移dt1、dt2、dt3、dt4、或使用其它合适的内部相关性识别在其它检测器数据流中的正确信号或数据流分段。当正确地识别信号或数据流分段时，可以对诸如分选的随后的操作适当表征微实体。
第二个问题涉及水滴从属关系或确定微实体最后将进入哪个液滴，从而可以对该液滴施加正确的电荷。第二个问题有些针对液滴细胞分选器，但是可以与形成传送微实体的规则隔开的封装的任何系统相关。例如，第二个问题可以涉及这样的微流体系统，其中将交替类型的不融和的流体引入通道，或者其中，周期地将空气泡引入通道，以形成分隔的流体封装。在这样的系统中，可以周期性地形成液滴或封装，并且微实体可以随机地进入流中。
在关于形成液滴和确定水滴从属关系的更详细内容中，可使用如图4A所示的水滴时钟信号401，以生成施加到流式细胞计喷嘴202中的换能器212的周期性换能器信号403（图4B）。如上所述，换能器212对喷嘴和流220中的流体施加声能。另外，可将充电脉冲（如图4C所示）施加到导电探针207，以对最后附着液滴228-0充电。换能器信号可与水滴时钟信号同步或相位锁定。可调节换能器信号的幅度以将断开点226移动到距喷嘴202的选定位置或距离处。充电脉冲可以与水滴时钟相位锁定，并且可以与水滴时钟和/或换能器信号同相或不同相。根据一些实施例，在系统200的实验或校准进程中，没有用于调节充电脉冲与换能器信号和/或水滴时钟之间的相位的可变相位调节。
根据一些实施例，当对最后附着液滴施加不同的分选充电脉冲时，调节到换能器的驱动幅度，以产生纯净（clean）的分选流。通过设置换能器信号幅度以形成纯净的分选流并且将断开点226保持在到喷嘴202近似固定的距离，可以使荧光校准珠或微实体行进通过系统以确定水滴从属关系。在一些校准进程中，可以以低速将校准实体供应通过系统，从而没有数据信号的混乱。通过检测断开区域中的校准实体，可确定对于特定微实体的来自询问点的哪些信号对应于断开区域中的液滴形成和/或液滴充电时间。在一些实施中，该确定可获得dtc的值。一旦进行确定，可以分配微实体的信号到特定液滴的从属关系。
当确定了水滴从属关系，可以对液滴中的一个或多个微实体进行分选决定。可以基于将处于液滴中的一个或多个微实体的信号的相关性进行分选决定。可通过将液滴电荷值分配到液滴实施分选的物理操作。例如，可以将液滴电荷值412分配到液滴。在一些实施例中，施加到液滴的电荷值可以是准离散的。例如，可通过对应于离散分选流的第一量来使各个电荷值分开，而可以调节最终电荷量以考虑更高阶效应，例如与相邻液滴的电容耦合。
在一些实施例中，可按照系统时钟周期数，例如，按照水滴时钟周期数（图4A）或模拟至数字（ADC）时钟周期数（图4D）或两种时钟周期数的组合，计算在校准进程中发现的平均时间偏移dt1、dt2、dt3、dt4、dtc。可通过使用以高于其它时钟的频率运行的系统时钟获得高分辨率的偏移。在流式细胞计中，ADC时钟可以以几十或几百MHz运行，而水滴时钟可以以比这些值慢几个数量级的频率运行。按照ADC时钟周期数计算的延迟值将提供更高的分辨率。在一些实施中，可按照更慢的时钟周期数（例如水滴时钟周期数）来计算长延迟。在一些情况中，可按照较慢和较快的时钟周期数的混合来计算长延迟。例如，如果水滴时钟的周期为Td，ADC时钟的周期为TADC，则可以将延迟计算为
dt＝NTd+aTd     (1）
其中N为整数，a表示小数：0<a<1。可将值a计算为
a＝MTADC/Td         （2）
其中M为整数。组合公式1-2获得dt=NT_d+MT_ADC
III.数据获取
如上所述，时间偏移或延迟dt1、dt2、dt3、dt4可用于选择用于观测定时窗310或数据流分段的适当时间，其中期望来自微实体的信号在每个询问点出现并且对于已经启动触发事件（例如，穿过B1处的触发询问点）的相同微实体出现。在观测窗内出现的数据将归属于该微实体，并被处理以确定对该微实体的操作（例如，分选操作）。如上所述，由于可以缓存来自询问点的数据，任意询问点可用作触发点，并且用于单个微实体的观测窗可在触发信号之前和之后出现。
在一些实施例中，观测窗310识别来自检测器D1、D2、...D5的数据流中的时间分段，在所述时间分段，信号将被归属于特定微实体。例如，系统中的检测器D1、D2、...D5的每一个可输出包括通过穿过相应询问束B1、B2、...B5的微实体生成的信号的样本的数据流。例如，如图4D-4F所示，可以以根据模拟到数字转换器时钟的高数据速率获取样本。穿过将用作触发器的选定束（在该实例中为束B1）的微实体，将生成触发信号（例如脉冲），其可被检测并用于识别微实体到达触发器询问点的到达时间。在下面的示出实例中，用于标记到达时间或参考时间t01的触发信号的特征被取为触发脉冲的峰，在其它实施例中可以使用其它值，例如由系统用户选择的上升沿阈值水平、脉冲积分的阈值。触发事件的参考时间t01可与时钟周期数相关，如图4D所示。
通过对微实体建立参考时间t01，并从校准进程获知延迟（例如延迟dt2），可以按照时钟周期数或数据样本确定在数据流中的何处查找移动通过系统的相同微实体的后续事件（例如穿过束B3）。例如，并参考图4D-4F，第一微实体可横穿触发束（在该实例中取为束B1），并在时间t01生成触发脉冲410。触发脉冲410可在来自检测器D1的数据流中被检测。数据流可包括其它微实体的后续触发脉冲412、414。通过获知延迟时间dt2，可以在来自检测器D3的数据流中找到当第一微实体横穿束B3时由其生成的后续信号。后续信号将出现在横跨微实体在束B3处的期望到达时间t21的观测窗310中。观测窗将从触发事件t01偏移对应于延迟时间dt2的多个时钟周期数（在该实例中为ADC时钟周期数）。在该观测窗310中，脉冲422表现为到达得晚，如图4F所示。对应于其它微实体的其它脉冲也出现在图4F的轨迹中。
从检测器D1、D2、...D5接收的数据可被存储在临时数据存储装置中，例如FIFO缓存器中。在一些实施例中，可将数据存储在长期数据存储装置中。由于数据被存储，在操作和处理数据时存在更大的灵活性，因为不需要实时处理。例如，通过数据存储，可以从选定事件在时间上向后和向前查找。根据一些实施例，根据时钟周期数确定观测窗的位置可以变换为根据存储器地址或存储器位置确定数据信号的位置（例如，假设将数据顺序存储在地址空间中或者与观测时间相关联地存储）。时间的延迟值，例如dt2，可被变换为地址空间中的偏移。
在这样的实施方式中，流220中的任何束横穿事件可以用作特定微实体的参考触发事件。例如，针对流中横穿第二束B2的微实体检测的信号可用作触发事件。通过以对应于适当的延迟时间的地址偏移在地址空间中向前和向后查找，可以找到对于相同微实体的其它束横穿的对应数据。
确定数据获取观测窗310的定时或位置确保对于穿过系统的每个微实体正确记录信号（可包括没有信号）。如图3E所示，一些微实体在询问点可以不生成信号，例如，不提供超过背景水平的可测量信号。例如，标记样本中的一些微实体的荧光标签可以对一些询问束敏感而对其它询问束不敏感。在背景噪声水平之上的可识别信号的缺失对于分类或表征微实体可能是重要的。从而，对观测窗310的正确确定可以帮助确定在该询问点或横穿束的微实体的信号的缺失。
IV.运动相关误差
在具有含移动微实体的流的系统中，如流式细胞计，微实体225到达沿流的询问点可能与期望到达时间存在偏差。例如，即使流220的流速可以基本恒定，并且通过校准仔细确定时间偏移dt1、dt2、dt3、dt4，然而对于由微实体启动的触发事件，仍可能存在微实体到达第二询问点和各个询问点或沿流220的位置的不确定性。微实体的到达时间或位置的不确定性可能来源于一种或多种运动相关误差。
移动微实体的运动相关误差可包括任何类型和形式的误差，该误差造成，在给定在第一位置或第一时间检测到微实体的情况下，相同微实体在第二位置或第二时间（或其它位置或其它时间）与期望到达时间（或多个时间）或期望位置（或多个位置）的偏离。移动微实体的运动相关误差可由实际和/或表观的原因导致。实际运动相关误差可包括流中的微实体从期望位置的物理移动，例如，耦合到流的能量在微实体沿流移动时可以扰动或抖动流中的微实体的位置，即使在一些实施例中，流可以没有湍流地流动。表观运动相关误差可包括观测影响，其在没有相应的物理运动时给出微实体的移动。表观运动相关误差的实例可包括但不限于，电子检测影响（例如，可变电子延迟，数字化误差）和光学影响（例如，由于波状流体流导致的可变透镜化影响）。
在流式细胞计领域，与微实体在询问点的到达时间的不确定性相关联的运动相关误差可称为“脉冲抖动”。脉冲抖动提供了商业系统中的移动微实体的运动相关误差的一个实例，并且当在示波器上观察来自微实体的信号脉冲时可以看到。脉冲抖动可由实际和/或表观原因导致。如图4F所示，相对于期望峰值到达时间t21，通过脉冲422的迟到示出脉冲抖动的实例。本发明人已经发现，在不通过换能器212对流220施加周期性能量的情况下下可以观测到脉冲抖动，并且在对流施加周期性能量的情况下，该脉冲抖动明显变得更加显著。在目前商业使用的全部多激光器、SIA细胞分选器中都在一定程度上观测到脉冲抖动。
图5示出脉冲抖动的实例。在图5中，在图中重叠了在相同的下游位置观测的并从上游位置的标称相同的条件触发的来自多个微实体的多个脉冲。图上的脉冲位置的变化反映了在下游位置的到达时间的变化，并表示系统中的运动相关误差。在一些系统中，到达时间的变化可大于几个脉冲长度。图上的脉冲位置的变化反映获知微实体在选定时间的位置的不确定性。
在图3C-3G中，由垂直虚线表示的观测窗310的尺寸被设置为反映在每个询问束位置处的脉冲抖动的量。如上所述，观测窗是其中可以采集数据并将数据归于启动触发脉冲的相同微实体的时间。观测窗随着到喷嘴的距离而尺寸增加，这反映了脉冲抖动随着到喷嘴的距离而增加的观测。更大量的脉冲抖动需要观测窗310的更大的时间跨度。
在各个实施例中，观测窗310的尺寸被设置为，使得由微实体在相应询问点生成的大多数脉冲在观测窗310中被捕获。在一些实施例中，观测窗310的尺寸被设置为，使得由微实体在相应询问点生成的脉冲的期望部分在观测窗310中被捕获。在各个实施例中，期望部分可以为在大约0.95和大约1之间、在大约0.90和大约0.95之间、在大约0.8和大约0.90之间、在大约0.7和大约0.8之间、或大约0.6和大约0.7之间。在一些实施中，如果脉冲到达一个位置的时间分布由标准偏差σ表征，则观测窗310的尺寸可以设置为与脉冲的平均宽带wavg成比例的第一值加上与标准偏差成比例的第二值。例如，观测窗的宽度W可通过下面的公式表示：W=cwavg+fσ，其中c是大于0的实数，且f是大于0的实数。如果一项显著大于另一项，实施例可具有通过cwavg或fσ确定的窗口宽度。
可以理解，运动相关误差（例如，流式细胞计中的脉冲抖动）可以对系统操作和处理量施加限制。例如，为了将来自询问点的数据正确归属于正确的微实体，每个观测窗优选不与数据流中的前一个窗口或连续的观测窗显著重叠。在一些情况中，少量的重叠可以接受。在一些实施方式中，观测窗可以没有重叠，并且在各个观测窗之间可能存在缓存或死区时间。保持小重叠或防止观测窗的重叠可对微实体处理速率和/或询问区域222的尺寸施加限制。例如，流220中的微实体必须被足够远地隔开，从而观测结果可以被正确地去卷积并以期望的概率或精度被分配到关注的微实体。在一些实施中，用于连续微实体的窗最小地重叠或者不重叠。可选地，或者另外，询问点必须被限制在喷嘴开口208附近的区域，在该区域中脉冲抖动的量在可接受的限制内。
发现当测量脉冲来自最接近流出开口208的询问点时脉冲抖动最小，并且相应的数据获取窗对于开口附近的询问点可具有最短的持续时间，同时保持期望水平的测量精度。本发明人已经观测到，抖动随着到喷嘴202的距离略微增加。该抖动的增加可对询问区域222的尺寸施加限制。例如，询问区域可被限制为这样的尺寸，使得对于给定处理量，抖动不导致用于给定处理量的观测窗的显著重叠。在一些系统中，脉冲抖动可将总询问区域限制为靠近喷嘴开口208的流的小区域（约1mm或更小）。对询问区域的尺寸的限制可限制可用于系统中的询问束的数目，并且对询问束的数目的限制限制了对微实体的分析。
另外，本发明人已经发现，抖动随着施加到换能器212的驱动信号的幅度增加。该抖动的增加可对驱动换能器的强度施加限制。在一些应用中可能期望以较大的幅度驱动换能器212，因为这通常引起液滴形成的更好的均匀性和稳定性。然而，相关联的脉冲抖动的增加可能阻止以高处理量、更高幅度驱动换能器。
V.表征运动相关误差
进行了一组实验以表征液滴细胞分选器流式细胞计中的运动相关误差，并导出运动相关误差的模型。由于液滴生成基本是周期性的和正弦的（例如，由于施加到换能器212的驱动信号的形式），假设流式细胞计的运动相关误差可以与施加到流220的周期扰动相关，并且包含周期函数的模型可代表系统中的运动相关误差。因此，假设在沿流220的任何位置处的流式细胞计中的脉冲抖动J可以表现出与时间偏移Δtdc或相位
Figure BDA00003107272400232
的近似正弦相关性，所述时间偏移Δtdc或相位
Figure BDA00003107272400233
是从水滴时钟周期中的参考或特征点（例如前缘）到检测的“触发”脉冲上的参考点（例如，在脉冲的前沿或后沿上的脉冲峰值或选定阈值）测量的。
为了理解并不限于任何特定理论，该相关性可如下数学表示：
J～Asin(ΔtdcF+θ)    （3）
其中A是将确定的幅度，F表示液滴生成频率或水滴时钟的频率，以及θ表示将确定的相位。根据一些实施例，A和θ可在校准进程中对于系统实验地确定。在其它实施例中，A和θ可以在分析和处理有用样本期间通过系统累积数据时现场确定。对于本领域技术人员显而易见的，可以使用其它模型（例如，简单或更高阶的Bessel函数，更高阶三角函数、高斯函数、相同或不同函数的叠加），以表示具有移动微实体的系统中的运动相关误差。还可以认识到，可以存在流式细胞计中的运动相关误差的其它成分，诸如由于移动和波状流导致的在光学检测设备中引入的误差、以及在用于从询问点采集数据的数据获取系统的电子元件中引入的误差。误差的一些其它成分可以或不以相同的模型表示。一些可能需要不同的模型。
根据一个实施例，可在具有移动微实体流的系统上将用于流式细胞计的抖动J测量为从观测窗中的期望到达时间到实际到达时间的时间偏移。可通过如上所述的校准过程建立到流220中的任意位置的期望到达时间（例如，时间t1、t2、t3、t4，参考图3C-3F）和观测窗310。用于选定位置的观测窗可被设置为横跨期望到达时间。实际到达时间可被测量为检测脉冲上的选定参考点的时间。在观测窗中何处出现期望到达时间（例如，开始、中间、结尾）以及在脉冲上使用哪个作为参考点（例如，前沿阈值、峰值、后沿阈值、积分阈值）的选择可以取决于特定实施方式，并且对于运动相关误差的模型开发不重要，只要选定的到达时间和参考点被一直使用。
为了理解应用于流式细胞计的本发明的各个方面，图6A-6B示意示出抖动J和与水滴时钟的时间偏移Δtdc。图6A可表示从触发询问点接收的信号，以及图6B可表示从沿流220的另一个询问点接收的信号。在一些实施例中，两个值都可以按照时钟循环被记录并通过将时钟周期数乘以时钟周期而转换为实际时间。例如，如上文结合公式（1）-（2）所述，可以按照ADC时钟周期数、水滴时钟周期数、或其组合测量时间偏移Δtdc和抖动J。在该实例中，从水滴时钟的上升沿到检测的触发脉冲610的峰值来测量时间偏移Δtdc。从观测窗310的中心的期望到达时间ta到第二检测脉冲620的峰值测量抖动J。在一些实施例中，可由微实体在横穿第一束B3时生成触发脉冲610，且第二检测脉冲620可以是从流220中的另一个询问点检测的脉冲。在一些实施例中，可使用脉冲上的其它参考点、水滴时钟以及观测窗。如果在脉冲抖动与液滴生成之间的假设关系正确，则观测到的多个微实体的抖动J将根据构建的误差模型（例如公式3的模型）而与测量的时间偏移Δtdc相关。
图7A示出对于流式细胞计中的多个微实体记录的作为时间偏移Δtdc或相位
Figure BDA00003107272400241
的函数的脉冲抖动幅度J的测量结果。图7A所示的抖动数据的图可以称为“抖动图”。利用校准珠在上述sy3200TM系统上进行测量。喷嘴开口直径为70微米，并且流体上的压力为50psi。液滴生成频率为84.9kHz，并且利用离喷嘴最远（离喷嘴大约2mm）的询问束进行抖动测量。获得的抖动图清楚示出在抖动振幅J与触发脉冲定时的偏移Δtdc之间的近似正弦的相关性。
图7B示出在与图7A相同的条件下的抖动测量的结果，但是其中利用距离喷嘴212第二远的询问束进行测量。抖动图示出抖动振幅的变化（在更靠近喷嘴处减小），但抖动图的相位θ没有变化。对于图7A和图7B中的图，已经减去与两个询问点之间的沿流的延迟相关的相位延迟。结果表示，无论是由实际还是表观原因导致，运动相关误差以相同的相位沿流随着微实体行进并增大振幅。
图7C中示出抖动图的相位θ的变化。对于该试验，液滴频率被改变至77.9kHz，并且所有其它条件相同。频率变化导致抖动图相位的变化。当喷嘴212的位置相对用作触发器的询问束的位置变化时，也观测到抖动图相位的变化。当换能器驱动信号的振幅变化时，未观测到在任意询问点的抖动图相位的变化。
当换能器驱动信号的振幅变化时，观测到抖动图的振幅的变化。另外，如上所述，抖动图的振幅的变化与询问束到喷嘴的距离相关。（比较图7A和7B）在距离喷嘴较远的束测量的抖动图的振幅大于在距离喷嘴较近的束测量的抖动图的振幅。
图7A-7C所示的抖动图结果示出，在流式细胞计系统或类似配置的任何系统中存在可以校准或补偿的微实体的运动相关误差的预测或确定分量。例如，如果对触发脉冲测量到距离水滴时钟信号上的特征点的时间偏移Δtdc，则可以从公式3和/或从抖动图的数据计算用于调节微实体到达流220中的另一个位置的到达时间的校正量。所计算的校正量可用于例如当微实体穿过系统时调节用于微实体的每个询问点的观测窗310的定时。这样的调节可用于重新使各个观测定时窗以各个脉冲的期望出现处为中心。关于从询问点传出的数据流，可以使用计算的校正量以调节数据流中的数据分段的位置或存储对应于相同微实体的信号的数据的存储器地址。因为存在脉冲抖动的预测分量，可在系统的数据获取设备中补偿该预测分量。运动相关误差自身不能被物理校正。可以理解，可以在其它系统中进行用于运动相关误差的补偿，在这些其它系统中，存在可以建模的运动相关误差的预测或确定分量，例如，由诸如公式3的数学模型表示。
VI．抖动补偿
如上面针对流式细胞计所描述的，可以为系统表征和补偿运动相关误差的与能量到流220的周期耦合相关联的分量。类似的技术可以应用到其它系统，例如微流体系统中。图8A-8B示出了根据在流式细胞计中的一个实施例的抖动补偿结果。图8A示出了在按照没有抖动补偿的传统方式操作的流式细胞计中的一个询问点测量的脉冲抖动。在各个脉冲的重叠中可以看到，微实体按照一定时间分布到达该询问点。对于没有补偿的系统，图8A的脉冲抖动大约为1微秒并且大于所记录的脉冲的半高宽值（full-width-half-maximum value）。
图8B示出了针对实施了抖动补偿的相同系统记录的类似数据。通过使用误差模型（例如根据公式3的模型）实现抖动的补偿，并且调节用来检测询问点处的脉冲的数据获取观测窗的定时，从而使残差最小化。下面描述在流式细胞计系统中实现运动相关误差补偿的进一步细节。对于图8B的抖动补偿数据，明显的抖动减小到脉冲宽度的一部分。图8B的结果表明可以在数据获取系统中对脉冲抖动的主要由通过换能器212引入的周期性扰动导致的分量进行补偿。脉冲抖动可能存在不能通过上述技术补偿的残留分量（例如背景抖动），这可以通过公式3的替换公式来补偿，或者通过等价的技术来减小。
应当理解，可以在系统的数据获取装置中执行对运动相关误差的补偿。在物理系统中仍然会出现微实体到达时间的变化。在一些实施例中，误差补偿包括改变观测脉冲的方法，例如调节观测窗的定时或者调节对微实体执行操作的定时。
可对系统进行进一步改善。例如当在系统中校正运动相关误差的主要来源时，更高阶运动误差可以变得可观测。按照这里所述的处理，还可以表征、建模和补偿更高阶误差。
从上述讨论可以理解，对移动流中的微实体的运动相关误差的测量和表征可用于诊断和改善系统性能。在由共有专利实体同时提交名为“Characterization of Motion-Related Error in a Stream of MovingMicro-Entities”的专利申请中描述了用于测量和表征运动相关误差的技术，其整体在此引入作为参考。
图9的流程图示出了用于表征和补偿移动微实体的流中的运动相关误差的方法1100的一个实施例。图9所示的步骤不必全部包括在用于表征运动相关误差的方法的各个实施例中。一些实施例可以仅包括所示步骤的子集。一些实施例可包括图中未示出而在这里描述的其它步骤。尽管图9中的步骤涉及表征和补偿流式细胞计中的脉冲抖动，可以理解，该方法包括用于表征其它系统中的运动相关误差的等同实施例，所述其它系统诸如为微流体系统和微加工系统。
根据一个实施例，用于表征和补偿运动相关误差的方法1100可包括启动1105系统，在该系统中微实体从第一位置移动到第二位置或一个或多个其它位置。系统可以是流式细胞计、微流体系统或被配置为对移动微实体进行测量和/或操作的任何系统。根据一些实施例，系统可以是其中在流体流中形成微结构的系统，例如，在Doyle等人的美国专利申请公开No.2010/0172898中描述的系统。该方法可包括初始化1110系统的步骤。初始化动作可包括验证和/或设置系统的运行参数，例如，设置移动微实体的速度。在流式细胞计中，初始化步骤1110的动作可包括对上述信号相关性和水滴从属关系进行校准。
方法1100还可以包括测量1120运动相关误差，诸如脉冲抖动数据。测量的运动相关误差可表示实际到达时间相对于期望到达时间的偏差或变化。微实体到达选定位置的期望到达时间可基于由微实体在触发位置导致的单独的检测事件（例如，触发事件）。可基于其它因素，例如微实体的平均速度或在“触发”位置与测量位置之间的平均传输时间，相对于期望到达时间评估实际到达时间的偏差。除了测量实际到达时间的偏差之外，测量1120运动相关误差的动作可包括测量触发事件的方面（例如，关于系统度量，诸如系统参考信号，或系统时钟，脉冲分布、距离或波长，测量触发事件的发生）。在液滴细胞分选流式细胞计中，触发事件的方面可以是其相对于用于对液滴的形成定时的水滴时钟的参考点（例如上升沿）的发生时间。在其它实施例中，可相对于将能量耦合到传送微实体的流中的定时测量触发事件的定时。在一些实施例中，系统度量可包括感测的值，例如流的温度、流体的压力、流的瞬时速度、外部声学水平、对系统构件测量的加速度值。
继续流式细胞计的实例，测量1120运动相关误差可包括测量每个脉冲或微实体的下列值：当脉冲到达选定触发通道时的水滴时钟相位（或距离水滴时钟上的参考点的时间偏移）、以及对于选定询问点或与询问点相关的检测通道在期望到达时间和实际到达时间之间的时间偏移。可以存储采集的数据以用于在方法1100中的下一阶段的处理。
在描述可以对多个微实体的每一个进行测量时，应当理解，可以不必对流中每个微实体进行测量。该多个微实体可以是全部微实体的一个子集，并且可以指为了表征运动相关误差而获得的正确测量的一部分微实体。
可以存在几个可能的实施例，其中在选定触发通道测量脉冲到达时间。一种方法为测量每个脉冲的峰值位置。另一个实施例为测量脉冲的上升边沿通过选定阈值的点。另一个实施例包括测量脉冲的积分（即脉冲下的面积）通过选定阈值的点。
在一些实施例中，可在误差校准和表征进程中获取运动相关误差数据（诸如抖动图数据）。在一些实施方式中，可在液滴延迟校准进行之后执行误差校准和表征。在其它实施例中，在其中正分析微实体的实验中，可以现场获取运动相关误差数据。当现场采集运动相关误差数据时，误差表征和补偿可以是动态进程，其中基于观测数据周期地更新误差模型。因此，误差表征和补偿可以随时间变化以反映采集的误差数据中的任何变化。在一些实施例中，使用贝叶斯先验对于构建初始误差模型并在给定系统中的运动相关误差的进一步观测值时支持误差模型的更新是有用的。
在一些实施例中，对于流220中的任意询问点表征运动相关误差包括确定1130在选定询问点的运动相关误差E与可决定误差的可测量系统质量之间的关系。一旦测量了误差数据，可处理该数据以确定1130运动相关误差E与如上所述的用于流式细胞计的可测量系统变量（例如相对于系统参考信号的触发事件的偏移（Δtdc或
Figure BDA00003107272400281
））之间的关系。在一些实施例中，可通过误差校准和表征测量，例如如上文结合图7A-7C所述的测量，经验地确定E与Δtdc或
Figure BDA00003107272400282
之间的关系。校准珠或荧光细胞可运行通过系统以对于每个询问点获得如同图7A-7C的抖动图。在一些实施中，该关系可通过函数建模，例如识别适合测量的数据的函数。
可以对移动微实体流中的每个关注点确定1130运动相关误差E与可测量系统变量之间的关系。关注点可以是流中对微实体进行测量或对微实体进行操作的位置。在一些实施例中，可以关注点或关注位置进行多于一个的测量。例如，在流式细胞计中，对于单个询问点可以有几个检测通道。由于在一些情况中系统中的每个检测通道可表现出不同的误差（例如，不同量的脉冲抖动），可以对每个通道确定运动相关误差E与可测系统变量之间的关系。
用于表征和补偿运动相关误差的方法1100可包括确定1140是否已经对于每个关注点评估运动相关误差E与可测量系统变量（例如，Δtdc）之间的关系的动作，所述关注点例如为流中的各个询问点和与各个询问点相关的各个检测通道。在一些情况中，可能每个询问点具有多于一个的检测通道，并且方法1100可包括确定1140是否已经对每个检测通道评估运动相关误差E与可测量系统变量之间的关系的动作。如果未对流中的每个关注点和/或检测通道评估运动相关误差E与可测量系统变量之间的关系，则方法流程可返回到对流中的下一个关注点和/或检测通道测量1120运动相关误差的步骤。对于其中内插一个或多个关注点的数据的实施例，方法流程可返回到如图9中虚线箭头所示确定1130运动相关误差E与可测量系统变量之间的关系的步骤。可通过测量端点之间的内插技术确定误差补偿值。如果已经对流中的每个关注点和/或检测通道评估运动相关误差E与可测量系统变量之间的关系，则方法流程可前进到计算1150可用于表征和补偿运动相关误差的校正值的动作。
根据一些实施例，计算1150校正值包括进行对测量的误差数据的曲线拟合。可以在或不在一个或多个图上绘制测量的误差数据。可将例如公式3中的周期函数的一个或多个函数拟合到每个图或其相应数据，可在每个询问点使用所获得的拟合函数，以补偿在随后的实验进程中的运动相关误差。通过确定拟合函数的参数（例如，振幅、相位、延迟、系数），可以对每个接收的脉冲根据分析函数计算表征和补偿值。
可选地，可以将每个询问点的运动相关误差数据编制并存储到查找表（LUT）中，该查找表可用于每个询问点在后续实验进程中补偿脉冲抖动。在一些实施例中，计算1150校正值可包括，生成各个值的查找表，其中LUT的每个元素对应于触发相位或延迟分段（bin）,并且从用于相位或延迟分段的运动相关误差数据确定每个元素的值。用于确定LUT中的每个元素的值的可能方法包括但不限于，对该分段计算运动相关误差数据的平均值、中值及模式。对该阶段的可选补充为跟踪残差，并使其趋于零平均值。在一些实施例中，这可以通过重新计算标称偏移以使得LUT中的值为零平均值而完成。这可以通过计算测量数据的平均值、将其从数据减去，以及将其增加到标称偏移而完成。在其它实施例中，增加标称偏移的一部分，例如所述偏移的1/2，以获得到零平均值的收敛。
在一些实施例中，所计算的LUT值的数目可受到限制，其中每个LUT值可跨越一定范围的运动相关误差条件。继续流式细胞计的上述实例，计算的LUT值可限于离散的触发时间或相位偏移分段，例如，如图10A所示，36个分段每个分段跨越约10度的触发相位偏移值。从而，任何落入10度分段中的测量的触发偏移相位将接收相同的平均值。可以理解，在一些实施例中可以用时间表达触发偏移，例如，相对于水滴时钟信号上的参考点的时间延迟（Δtdc），或者可以用相位
Figure BDA00003107272400301
表达，例如水滴时钟循环的相位。
在一些实施例中，可以在正常系统操作期间，例如，在细胞分选运行期间，现场执行计算1150校正值的动作，或者在其它实施例中，其可以在正常操作之前（例如，在校准运行期间）执行。作为一个实例，可使用测量的触发偏移以从现场的函数拟合计算校正值。在上述流式细胞计的情况中，计算1150校正值的动作可包括计算可用于偏移时间观测窗310的偏移时间，所述时间观测窗310用于从每个询问点获取脉冲数据。计算的时间偏移可基于运动相关误差测量数据，并且可以使得计算的偏移时间补偿在选定位置测量的运动相关误差的平均量。计算的校正值可用于补偿测量的误差并然后弃用，或可以被存储到查找表中并后续根据测量的触发偏移进行访问以补偿运动相关误差。
在一些实施例中，用于补偿运动相关误差的方法1100包括当各个微实体穿过第一位置（例如触发询问点）时接收1160来自各个微实体的信号。各个接收信号可以表示各个微实体出现在第一位置，例如到达时间。该信号可以是由系统的至少一个处理器接收的数字信号或模拟信号。
用于补偿运动相关误差的方法还包括针对每个第一信号确定1170相对于系统的参数的各个微实体的值。根据一个实施例，该值可以是根据信号确定的各个微实体到达第一位置处的到达时间的时间偏移或相位偏移，其中相对于系统的度量（例如系统时钟、系统提供的参考信号、系统提供的周期信号）来测量时间偏移或相位偏移。在一些实施例中，系统的度量可以是感测的值，例如流的温度、流体的压力、流的即时速度、外部声学水平、针对系统部件测量的加速度值。在一些实施例中，该值可以与根据运动相关特征测量而得到的误差模型的分量相关联，该误差模型预测与各个微实体在沿着系统的流的一个或多个其它位置处的运动相关误差。
用于补偿运动相关误差的方法1100还包括应用1180所计算的校正值以基于所确定的值来补偿运动相关误差。可以基于误差模型即时计算（例如根据对运动相关误差的特征数据拟合的函数来计算）校正值，或者可以从LUT中选择校正值（例如使用所确定的值和/或流的其它位置作为LUT中的索引）。如上所述，在一些实施例中，可以应用校正值来偏移数据获取观测窗310。例如，系统处理器可以使用校正值从来自询问点的数据流中识别包括穿过该询问点的微实体的信号的一段数据。在一些实施方式中，可以应用校正值来更准确地确定微实体在流220中的位置。在一些实施例中，可以应用校正值来调节对微实体执行操作的时间，所述操作例如分选操作、成像操作、探测操作、转变操作等。
在一些实施例中，应用1180所计算的校正值还包括采用插值技术连同LUT一起来补偿运动相关误差。例如，如上所述，可以进行测量以针对10度分段的增量确定LUT值。然后，所确定的LUT值可以存储在LUT数据存储器中供以后使用。在接收1160到微实体的信号并确定了微实体的值以后，可以使用所确定的值作为表格中的索引来访问LUT。如果该值对应于LUT的已有的时间偏移或相位偏移，则返回对应的校正量。如果所确定的值不对应于已有的时间偏移或相位偏移，则可以使用插值处理来计算校正量。插值处理可以使用来自LUT的两个校正量，这两个校正量对应于距离根据接收信号确定的值最近的索引。可以通过处理器250基于线性或更高阶插值来计算插入的值。在一些实施例中，可以在现场可编程门阵列中执行插值。然后，可以使用插入的值（而不是特定的LUT值）来补偿微实体的运动相关误差。
根据一些实施例，可以高重复率针对穿过系统的各个微实体执行接收1160、确定1170以及应用1180的动作。在一些系统中，重复率可以高达每秒执行50000次，在一些实施例中大于约每秒执行50000次，在一些实施例中大于约每秒执行100000次，在一些实施例中大于约每秒执行150000次。在一些实施方式中，可以在沿着流的一个或多个其它位置针对单个微实体重复应用动作。例如，在确定1170动作中确定的值可以用于针对一个或多个其它位置中的每一个确定一个或多个校正值。
在一些实施例中，针对系统中的每个询问点和数据获取通道，可以并行地进行以下动作：测量1120运动相关误差、确定1130误差E和可测量的系统变量之间的关系、以及计算1150校正值。在一些实施例中，可以不需要确定1140是否已经对每个关注点评估了误差E和可测量的系统变量之间的关系的动作。
在一些实施中，可以实时地，例如，在其中分析微实体的正常分选进程期间，进行测量1120运动相关误差数据和确定1130运动相关误差E与可测量系统变量之间的关系的动作。实时获取的运动相关误差数据可被处理并在系统运行时用于计算误差补偿值。在一些实施例中，计算的值可以被随后用于补偿运动相关误差。在这样的实施方式中，用于测量和表征运动相关误差的单独校准进程不是必需的，可以在正常进程期间实时更新先前计算的任何运动相关误差补偿数据。在一些实施例中，使用贝叶斯统计可提供用于实时更新用于系统的误差模型的有效架构。
在一些实施例中，可以对每个关注点或关注位置和/或每个检测通道执行以下动作：测量1120运动相关误差数据、确定1130误差E与可测量系统变量之间的关系、以及计算1150校正值。可以单独地处理测量的数据。例如，可以将来自每个测量的数据拟合为函数形式、或编制到专用于对应于测量的数据的位置和/或通道的查找表中。
如上所述，当缓存来自询问点的数据时，流中的询问点的任一个询问点可以用作触发询问点。因此，系统可被配置为基于一个触发询问点（例如，自流式细胞计喷嘴沿流遇到的第一个询问点）确定校正值，以及随后改变操作以使用另一个询问点（例如，第二个询问点）作为触发询问点。在一些实施例中，可以不需要在系统上执行重新校准过程以确定新的触发询问点的校正值。而是，可以数学变换存储的校正值、函数拟合、或LUT以反映触发询问点中的变化。根据流式细胞计的一个实施例，如果校正值被实现为由在颗粒横穿触发点时确定的相位或时间偏移指示的LUT时，并且触发询问点从沿流的第一个询问点变化到沿流的第二个点，则可以应用下面的数学变换：首先，从基于每个分段的询问点LUT的每一个减去用于第二询问点的初始LUT值。对于第二询问点，这从其自身减去初始LUT值，而对于新的询问点正确地剩下LUT中的零值。第二，用于其它询问点的新LUT的每个必须被循环地（circularly）偏移多个分段，所偏移的多个分段等于第一和第二询问点之间的相位差（其对于当前水滴时钟频率和流速率被确定）。
由于可以对选定询问点使用多个数据获取通道，并且全部这些通道具有基本相等的运动相关误差特征，可以有利于使横穿这些共用通道的脉冲抖动测量数据一致。一个实施例包括选择位于询问点的共用通道中的一个（例如，对于检测的微实体具有最佳信噪比的通道），以用于确定运动相关误差值。根据另一个实施例，在计算校正值之前可组合对于位于询问点的全部通道测量的全部运动相关误差数据。另一个方法包括，对每个通道计算指示校正值适合测量数据的程度的量度，并从具有最佳量度的通道选择校正值。例如，量度可使用均方误差（MSE）或均方根误差（RMSE）算法，或者可以计算总残差。
在一些实施中，确定1130运动相关误差E与可测量系统变量之间的关系和/或计算1150校正值可包括测量第一询问点和第二询问点处的运动相关误差数据。在第一和第二询问点之间可存在一个或多个询问点。根据一些实施例，可从对于第一和第二询问点测量的数据内插用于第一与第二询问点之间的一个或多个询问点的数据。内插可符合线性模型、非线性模型、或选定函数形式。用于每个点的数据，无论测量的或内插的，可被拟合为函数形式，或可被编制到专用于每个询问点的查找表中。
在系统硬件和/或软件方面，参考图2，可通过组合一个或多个处理器250、数据获取装置230和存储器240实施运动相关误差表征和补偿。根据一个实施例，可通过调节数据获取窗的定时实施运动相关误差表征和补偿，所述数据获取窗在检测和处理触发事件之后由数据获取装置230打开用于检测脉冲。处理器250可识别触发事件，例如来自检测器D1的脉冲，并将触发脉冲上的参考点（例如前沿阈值）与水滴时钟上的参考点（例如上升沿）比较。基于这两个事件之间的时间差，处理器250可调节通过数据获取装置230打开的观测窗的定时，以接收来自其它询问点检测器D2、D3、D4、D5的脉冲数据。观测窗可以在时间上提前或延迟。
根据另一个实施例，从询问点检测器接收的全部数据可以临时存储到存储器240中，例如FIFO或环形缓存器中。处理器250可被配置为基于触发事件相对于水滴时钟的定时，例如基于时间偏移Δtdc或水滴时钟相位
Figure BDA00003107272400341
的测量值，从每个检测器数据流中选择数据分段。
如上所述，对于流式细胞计，抖动图存在依赖于一定系统条件的方面。抖动图的振幅可取决于从液滴换能器212耦合到流220中的能量的振幅或量和从喷嘴或换能器到询问点的距离。抖动图的相位可取决于水滴驱动信号的频率以及触发询问点到喷嘴的距离。从而，如果水滴驱动信号的振幅、水滴驱动信号的频率、以及/或者询问点到喷嘴的距离改变，则系统可能需要被重新校准以获得用于运动相关误差的更新的模型。在配置用于现场误差数据采集、表征和补偿的系统中，可连续更新表征和补偿值，从而系统跟踪水滴驱动信号的振幅、水滴驱动信号的频率、以及/或者询问点到喷嘴的距离的变化，从而可以不需要单独的误差重新校准进程。
表征和补偿运动相关误差具有用于具有移动微实体流的系统的几个有利方面。表征运动相关误差可通过使用补偿技术对误差建模并减小误差。在数据分析中补偿运动相关误差可减小位于任何询问点的观测窗310的宽度。减小观测窗可以允许对流中具有更近间隔的微实体进行检测、测量和分选，从而允许系统的更高处理量。
与在流式细胞计中表征和补偿运动相关误差（例如脉冲抖动）相关的一个益处涉及用于检测微实体的检测设备。因为可以补偿流式细胞计中的脉冲抖动，在系统中可以容许更大量的脉冲抖动。因此，询问点可以在流中隔开较远。在一些实施例中，询问区域可以在2毫米上延伸。增大的询问点间隔减少了来自每个询问点的信号的串扰，从而改善了每个询问点的信噪比。
表征和补偿运动相关误差还可以改善在对于每个微实体确定水滴从属关系中的精确度，因为更好地确定了在流220的将形成液滴的分段中微实体的位置距离中心位置的偏差。较窄的观测窗可以允许用于在流的将形成液滴的分段内确定微实体的位置的更精细的分辨率。以更精细的分辨率获知分段内的微实体的位置改善了对水滴从属关系的确定。对水滴从属关系的更精确的确定可以改善分选质量，例如，分选纯度。
还可以通过两个或多个的更好的分辨率改善分选纯度。一些将形成液滴的流分段可以承载一个以上的微实体。较窄的观测窗可允许更清晰地分辨两个或多个，例如，指示液滴中存在多于一个的微实体，而之前的系统仅可以指示一个实体。较窄观测窗还可以允许更好地确定，如果微实体在将形成的两个液滴之间的流分段边界附近，微实体最终进入哪个液滴。
如上所述，运动相关误差表征和补偿可以允许更大的询问区域222。例如，不能容许在距离喷嘴较大距离（例如大于约1mm，大于约1.5mm）处出现大量脉冲抖动的流式细胞计可以使用抖动表征和补偿技术以减小脉冲抖动的不希望的影响。更大的观测区域可允许使用额外的询问束以用于对微实体的更复杂的分析。
VII.对移动微实体流的调制
在上述流式细胞计实施例中，将声能量耦合到传送微实体的流以形成液滴，从而可以将微实体包含在可根据施加的电荷分选的液滴中。根据规则周期性循环将声能耦合到流220，但是微实体在随机时间到达流220中。将周期性声能施加到流有助于控制随机到达流中的各个微实体，例如有助于识别并分选每个微实体。尽管耦合声能对于流式细胞计以脉冲抖动的形式将一定量的运动相关误差引入系统中，该误差具有确定性质，例如如图7A-7C所示，这允许补偿误差以允许对微实体的更精确的分析和更高的系统处理量。
可以存在其它技术领域，其中对移动微实体流施加声能可以有利于微实体分析和控制。例如，在可以用于分析或制造微实体的微流体系统中，可将周期性声能耦合到微实体移动通过的通道中。尽管微实体可以随机到达通道中，上述周期性声能和误差补偿技术可以用于更好地确定通道中的各个微实体的位置、或者在对微实体进行操作的选定位置的时间。更精确地获知微实体在通道中的位置或者微实体到达选定位置的到达时间可以改善微流体系统的分析和/或分选操作。在一些实施例中，可通过位于微流体芯片上的集成声学换能器将声能耦合到通道中。
在一些实施例中，耦合到移动微实体流的声能可改善对微实体的成像。例如，系统可被配置为利用时间延迟集成（TDI）多像素成像器件成像在流中流动的微实体。在D.Basiji等人的美国专利No.6,608,680中描述了TDI成像的一个实例，其在此引入作为参考。更精确地获知在流中移动的微实体的连续位置允许更好地对准由移动微实体形成的连续图像。更好地对准连续图像可用于减小图像噪声并改善图像分辨率。
VIII.实例
通过说明，并且不限于任何特定实施方式，在下面的实例中描述了检测流式细胞计中的补偿脉冲抖动的实例。
实例1
在该实例中，液滴细胞分选器由称为反射采集光学元件-5（RCO-5）的光采集光学元件构成，其允许探测沿流220的五个询问点并从其采集光。（可从以下地址获得RCO-5， iCyt MissionTechnology，2100South Oak Street,Champaign,IL61820,USA）。RCO-5允许各询问点之间的各种间隔距离。在该实例中，细胞将从顶部询问点（＃1）行进到底部询问点（＃5）的路径长度（垂直距离）为2mm。通过具有70、100和126微米的开口直径的液滴细胞分选器喷嘴和从10m/s到40m/s的各个流速度测试RCO-5。通过使用来自iCyt反射系统的数据获取系统（产品BW-5）（可从以下地址获得，iCytMission Technology, 2100South Oak Street,Champaign,IL61820,USA）测量脉冲抖动。该数据获取系统以105MHZ数字地采样通过细胞分选器的检测器（Hamamatsu光电倍增管（PMT），具有45MHz的带宽）产生的模拟信号。通过对于在询问点＃2产生的信号和在询问点＃5的来自相同微实体的信号测量10％峰值高度的脉冲前沿的交叉点之间的Δt评估脉冲抖动。全部5个询问点在喷嘴开口208的大约2mm内。液滴直到距离喷嘴开口大约6mm处从流分离。
当液滴生成器未启动时（即，流未分离为液滴而进入废液器皿），观测到抖动小于100ns。当启动液滴生成器时，对于每个操作点记录各种抖动量（鞘层压力、液滴频率、液滴换能器驱动强度的组合）。观察到的脉冲抖动的最高值为3微妙的总变化（或脉冲的平均宽带的大约三倍）。
还进行了以下观测：（1）没有液滴生成时（即，只有流动流），任意两个点的Δt几乎恒定；（2）当开始液滴生成时，在第一与最后（间隔最远的）询问点之间观测到显著的Δt的变化，其在一些情况中可以是4个脉冲长度；（3）脉冲抖动的量或剧烈程度随着液滴换能器驱动振幅的增大而增大；以及（4）当液滴断开点移为更接近喷嘴的出口时，脉冲抖动的量或剧烈程度增加。
对在实验室实验中采集的数据进行分析，以确定：（a）正弦函数与数据的拟合程度，以及（b）是否能使用拟合的正弦函数充分地校正所观测到的抖动。图10A-10B示出了这些分析的结果。在图10A中，示出用于从距离喷嘴最远的询问点（询问点＃5）获得的测量的抖动图。记录数据把测量抖动（纳秒）表示为相对于水滴时钟信号上的参考点的触发脉冲偏移（以相位表示）的函数。图中所示为适配数据的最小二乘正弦曲线（黑线）和来自查找表计算的分段平均值（白色圆形）。在该实例中，通过各个离散相位角分段中的抖动数据的平均值计算LUT值。每个相位角度分段跨越大约10度。尽管拟合函数和LUT值都良好地表示测量的数据，LUT值可以更精确地表征抖动。例如，LUT值反映了抖动数据中的更高阶的影响，其在抖动图的最大和最小值附近更清楚可视。
在图10B中，使用对数据的正弦拟合校正了图10A的测量的抖动数据。结果是把脉冲抖动从大约1000纳秒有效地减小或补偿到了大约200纳秒。图11A和图11B中绘出的直方图对应于图10A和图10B的相应数据。
图11A的直方图表示系统中没有抖动补偿的脉冲分布在宽度大于约1000纳秒的窗口上。当对系统应用抖动补偿时，脉冲的FWHM分布小于约100纳秒，如图11B所示。
鉴于这些结果，数据获取窗口的宽度需要至少比脉冲宽度多1.5微秒（持续时间），从而捕获具有如10A和图11A中所示的抖动量的脉冲。然而，通过如图10B和11B所示补偿或校正抖动（例如通过调节数据获取窗的定时），窗口宽度只需要比脉冲宽度大100纳秒。
IX.结论
在本申请中引用的全部文献和类似材料，包括但不限于，专利、专利申请、文章、书、专著和网页，无关于这些文献和类似材料的格式，其整体在此详细引入作为参考。在一个或多个引入的文献和类似材料与本申请不同或矛盾的情况中，包括但不限于定义的术语、术语使用、描述的技术等，以本申请为主。
这里所用的章节标题仅为了有条理性，而不应理解为以任何方式限制所述主题。
尽管结合了各个实施例和实例描述了本发明，本发明不限于这些实施例或实例。相反，如本领域技术人员可以理解，本发明包括各种替换、修改和等同物。
尽管在本文已经描述和示出了本发明各个实施例，本领域普通技术人员容易想到用于进行这里所述的功能和/获得结果和/或一个或多个优点的多种其它方法和/或结构，并且认为该变化和/或修改的每个落入这里所述的本发明实施例的范围内。更通常地，本领域技术人员容易理解，这里所述的全部参数、尺寸、材料和配置都是示例的，实际参数、尺寸、测量和/或配置将取决于具体应用或者使用本发明技术的应用。本领域技术人员可以仅利用常规实验而了解或者能够确定这里所述的具体本发明实施例的多个等同物。从而，可以理解，前述实施例仅以示例的方式示出，在所附权利要求及其等同物的范围内，本发明实施例可以以不同于具体描述和要求的方式实施。本发明的新颖实施例涉及这里所述的每个单独特征、系统、物件、材料、套装、系统升级和/或方法。另外，如果所述特征、系统、物件、材料、套装、系统升级和/或方法不是相互一致的，两个或多个所述特征、系统、物件、材料、套装和/或方法中的任意组合被包括在本发明范围内。
本发明的上述实施例可以以多种方法中的任一方法实现。例如，一些实施例可以利用硬件、软件或其组合实现。当实施例的任何方面至少部分地以软件实现时，可以在任意合适的处理器或处理器的集合上执行软件代码，所述处理器或者被设置在单个计算机中或者分布在多个计算机中。
在这方面，本发明各方面，例如，分析检测脉冲、校准和补偿脉冲抖动，可以指示部分地实施和/或实现为以一个或多个程序编码的计算机可读存储介质（或多个计算机可读存储介质）（例如，计算机存储器，一个或多个软盘、压缩盘、光盘、磁带、闪存、现场可编程门阵列（FPGA）中的电路配置或其它半导体器件、或其它有形计算机存储介质或非瞬时介质），所述程序当在一个或多个计算机或其它处理器上执行时进行实现上述技术的各个实施例的方法。计算机可读介质可以是可输运的，从而可以将其上存储的程序或多个程序装载到一个或多个不同的计算机或其它处理器上，以实现上述本发明的各个方面。作为一个实例，在一些实施例中，可以在FPGA中完整或至少部分地实现对数据的处理和系统操作。
这里所用的术语“程序”或“软件”概括地表示任何类型的计算机代码或计算机可执行指令集，其可被用于编程计算机或其它处理器以如上所述实现本发明各个方面。另外，可以理解，根据该实施例的一个方面，在执行时进行本发明方法的一个或多个计算机程序不需要位于单个计算机或处理器上，而是可以以模块方式分布于多个不同计算机或处理器中以实施本发明的各个方面。
计算机可执行指令可以为多种形式，诸如通过一个或多个计算机或其它装置执行的程序模块。通常，程序模块包括进行特定任何或实现特定抽象数据类型的路线、程序、对象、组成、数据结构等。通常，程序模块的功能性可根据需要组合或分布于各个实施例中。
另外，这里所述的技术可实现为已经提供其至少一个实例的方法。作为该方法的部分进行的步骤可以以任何方式排序。因此，实施例可以以不同于所示顺序的顺序进行的步骤构成，其中可包括同时进行即使在所述实施例中示出为顺序步骤的一些步骤。
这里定义和使用的全部定义应理解为在词汇定义上控制引入作为参考的文献中的定义和/或定义术语的普通意义。
说明书和权利要求中使用的不定冠词“一”和“一个”除非明确地相反的指出，应理解为表示“至少一个”。
说明书和权利要求中所用的短语“和/或”应理解为表示结合的元件的“任一个或同时二者”，即在一些情况中结合出现、在其它情况中不一起出现的元件。以“和/或”列出的多个元件应以相同的方式构成，即“一个或多个”这样结合的元件。除了由“和/或”短句具体识别的元件之外的其它元件可以可选地出现，无论其与这些具体识别的元件是否相关。从而，作阿非限制性实例，结合开放性表述（诸如“包括”）的表示“A和/或B”在一个实施例中可以仅指A（可选地包括除B之外的元件）；在另一个实施例中，仅指B（可选地包括除A以外的元件）；在另一个实施例中，指A和B（可选地包括其它元件）；等。
在说明书和权利要求中，所用的“或者”应理解为具有与上述“和/或”相同的意思。例如，当在列表中分开项目，“或者”或“和/或”应理解为包括的，即，在多个元件或其列表中，包括至少一个，还包括一个以上，并可选地包括其它未列出的项目。仅清楚指出相反意思的术语，诸如“仅一个”或“恰好一个”，或者，在权利要求中使用的，“由...构成”将表示包括多个元件或元件列表中的仅仅一个元件。通常，这里所用的术语“或者”在由诸如“任一个”、“之一”、“仅一个”或者“恰好一个”的排他性术语在前限定时仅应解释为指示排他的可选项（即，“一个或另一个，但不是同时二者”）。当在权利要求中使用时，“基本由...构成”应具有如在专利法的领域中使用的普通意思。
在说明书和权利要求中，在指示一个或多个元件的列表中使用的短语“至少一个”应理解为表示选自元件列表中的任何一个或多个元件的至少一个元件，而不一定包括在元件列表中具体列出的每个元件的至少一个，并且不排除元件列表中的元件的任意组合。该定义还允许，还可以可选地存在除了在短语“至少一个”表示的元件列表中具体识别的元件之外的元件，无论其与所述具体识别的元件是否相关。从而，作为非限制性实例，“A和B中的至少一个”（或，等同地，“A或B的至少一个”，或，等同地，“A和/或B的至少一个”）在一个实施例中可以指至少一个，可选地包括多于一个的A，而不存在B（并可选地包括除B之外的元件）；在另一个实施例中，可以指至少一个，可选地包括多于一个的A，和至少一个，可选地包括多于一个的B（并可选地包括其它元件），等。
在权利要求以及上述说明书中，所有的过渡短语，诸如“包括”、“包含”、“承载”、“具有”、“含有”、“涉及”、“保存”、“由...组成”等应理解为开放性的，即表示包括但不限于。如美国专利局专利审查手册第2111.03部分中所述，仅仅过渡短语“由...构成”和“主要由...构成”应分别为闭合或半闭合型的过渡短语。
权利要求不应理解为限于所述顺序或要素，除非陈述该影响。应理解，在不偏离所附权利要求的精神和范围下，本领域技术人员可以进行形式和内容上的各种变化。在此要求落入下面的权利要求及其等同物的精神和范围内的全部实施例。

Claims (24)

1.一种用于补偿与在第一位置与至少第二位置之间移动的微实体相关联的运动相关误差的系统，该系统包括：
检测装置，被配置为当微实体穿过第一位置时产生第一信号，以及当该微实体穿过至少第二位置时产生至少第二信号；
一个或多个处理器，被配置为：
至少接收所述第一信号；
根据所述第一信号相对于系统的参数确定一个值；以及
按照补偿与微实体相关联的运动相关误差的校正量来调节在第二位置处对微实体进行观测或操作的观测时间或操作时间，其中根据误差模型来确定所述校正量，所述误差模型根据所确定的值来预测与所述微实体相关联的运动相关误差。
2.根据权利要求1所述的系统，其中所述一个或多个处理器还被配置为：按照补偿与所述微实体相关联的运动相关误差的至少一个其它校正量来调节在相应的一个或多个其它位置处对所述微实体进行观测或操作的相应的一个或多个观测时间或操作时间，其中使用所确定的值根据所述误差模型来确定所述至少一个其它校正量。
3.根据权利要求1所述的系统，其中所述微实体在流式细胞计或微流体设备的流体流中移动。
4.根据权利要求1所述的系统，其中相对于所述系统的参数确定的值是测量的在所述微实体到达第一位置的到达时间与系统的参考信号之间的时间偏移或相位偏移。
5.根据权利要求4所述的系统，其中使用由系统提供的时钟信号的时钟周期数来测量所述时间偏移或相位偏移。
6.根据权利要求4所述的系统，其中所述参考信号包括由系统产生的周期信号。
7.根据权利要求6所述的系统，其中使用所述周期信号把能量耦合到在第一位置和第二位置之间传送所述微实体的流中。
8.根据权利要求7所述的系统，其中根据用于对流式细胞计中的液滴的形成进行计时的水滴时钟来得到所述周期信号。
9.根据权利要求1所述的系统，还包括换能器，该换能器被配置为周期性地把能量耦合到在第一和第二位置之间传送所述微实体的流中，其中相对于驱动所述换能器的信号的特征来确定所述值。
10.根据权利要求1所述的系统，其中第一和至少第二位置间隔超过大约1毫米。
11.根据权利要求10所述的系统，还包括一个或多个辐射源，所述一个或多个辐射源被配置为在所述第一和至少第二位置探测所述微实体。
12.一种补偿与在系统中的第一位置和至少第二位置之间移动的微实体相关联的运动相关误差的方法，该方法包括：
利用检测装置检测表示第一微实体出现在第一位置处的信号；
利用至少一个处理器，根据针对所述第一微实体的所述信号，确定相对于系统的参数的值；以及
利用所述至少一个处理器，按照补偿与所述第一微实体相关联的运动相关误差的校正量来调节用于在第二位置处对所述第一微实体进行观测或操作的观测时间或操作时间，其中根据误差模型确定所述校正量，所述误差模型根据所确定的值来预测与第一微实体相关联的运动相关误差。
13.根据权利要求12所述的方法，还包括：
针对所述第一微实体移动到的其它位置重复调节动作，其中针对每个其它位置相应地确定不同的校正量。
14.根据权利要求13所述的方法，其中以大于每秒50000次的速率执行检测、确定和调节动作。
15.根据权利要求12所述的方法，还包括：
把能量周期性地耦合到传送所述微实体的流中；以及
相对于驱动能量的周期性耦合的信号的特征而确定所述值。
16.根据权利要求12所述的方法，其中调节观测时间或操作时间包括：
从查找表中选择校正量，所述查找表与所确定的值相关联地存储校正量；以及
利用所述校正量补正所述观测时间或操作时间。
17.根据权利要求16所述的方法，其中补正所述观测时间包括：
根据调节后的观测时间设置数据获取窗口的定时，从而从所接收的来自检测装置的数据流中选择针对第二位置的数据；以及
将从所述数据流中选择的数据识别为表示第二位置处的第一微实体的特征信号。
18.根据权利要求17所述的方法，还包括通过至少一个微处理器，基于所述特征信号来决定对所述第一微实体的操作。
19.根据权利要求18所述的方法，其中所述决定包括流式细胞计中的分选决定。
20.根据权利要求17所述的方法，还包括基于所述特征信号对所述第一微实体进行分选或分类。
21.根据权利要求12所述的方法，其中所述观测时间是拍摄所述微实体的图像的时间。
22.根据权利要求12所述的方法，其中所述校正量与具有包括所确定的值的自变量的函数成比例。
23.根据权利要求22所述的方法，其中所述函数是正弦的。
24.一种实现机器可读指令的制造的数据存储设备，当由至少一个处理器执行所述机器可读指令时，使所述至少一个处理器被配置为：
接收表示第一微实体出现在系统中的第一位置处的信号，其中微实体在系统中的第一位置和至少第二位置之间移动；
根据第一信号相对于所述系统的参数确定一个值；以及
按照补偿与所述微实体相关联的运动相关误差的校正量来调节在第二位置处对所述微实体进行观测或操作的观测时间或操作时间，其中根据误差模型确定所述校正量，所述误差模型根据所确定的值来预测与所述微实体相关联的运动相关误差。
CN201280003477.4A 2011-08-25 2012-08-24 移动微实体流中的运动相关误差的补偿 Active CN103180858B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161527340P 2011-08-25 2011-08-25
US61/527,340 2011-08-25
PCT/US2012/052204 WO2013028948A1 (en) 2011-08-25 2012-08-24 Compensation of motion-related error in a stream of moving mirco-entities

Publications (2)

Publication Number Publication Date
CN103180858A true CN103180858A (zh) 2013-06-26
CN103180858B CN103180858B (zh) 2016-01-20

Family

ID=47746889

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201280003477.4A Active CN103180858B (zh) 2011-08-25 2012-08-24 移动微实体流中的运动相关误差的补偿
CN201280003432.7A Active CN103180713B (zh) 2011-08-25 2012-08-24 移动微实体流中的运动相关误差的表征

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201280003432.7A Active CN103180713B (zh) 2011-08-25 2012-08-24 移动微实体流中的运动相关误差的表征

Country Status (4)

Country Link
US (2) US20140220621A1 (zh)
EP (2) EP2732395B1 (zh)
CN (2) CN103180858B (zh)
WO (2) WO2013028947A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015154288A1 (zh) * 2014-04-10 2015-10-15 深圳迈瑞生物医疗电子股份有限公司 一种进行自动补偿的方法、装置及相应的流式细胞分析仪
CN108548817A (zh) * 2018-03-26 2018-09-18 清华大学 多相浮射流实验发生装置及油滴气泡阴影图像处理方法
CN116957992A (zh) * 2023-09-20 2023-10-27 南京木木西里科技有限公司 基于特征追踪的实时显微图像防抖方法

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201019521D0 (en) * 2010-11-18 2010-12-29 Durham Scient Crystals Ltd Radiation detection
EP2732395B1 (en) 2011-08-25 2018-04-04 Sony Corporation Compensation of motion-related error in a stream of moving mirco-entities
US10190959B2 (en) * 2012-02-09 2019-01-29 Beckman Coulter, Inc. Sorting flow cytometer
JP5924077B2 (ja) 2012-03-30 2016-05-25 ソニー株式会社 微小粒子分取装置及び微小粒子分取装置における軌道方向判定方法
US9915935B2 (en) * 2012-03-30 2018-03-13 Sony Corporation Microchip-type optical measuring apparatus and optical position adjusting method thereof
JP5601424B2 (ja) 2012-03-30 2014-10-08 ソニー株式会社 微小粒子分取装置及び該装置における流体ストリーム最適化方法
EP3910318B1 (en) 2013-01-28 2024-06-12 Sony Group Corporation Microparticle sorting device, and method and program for sorting microparticles
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
US10126225B2 (en) 2013-08-16 2018-11-13 Bio-Rad Laboratories, Inc. Timing and/or phase adjustment of the separation and/or charging of drops from a fluid stream in a flow cytometer
JP6447506B2 (ja) 2013-10-16 2019-01-09 ソニー株式会社 粒子分取装置及び粒子分取方法
WO2015122071A1 (ja) 2014-02-13 2015-08-20 ソニー株式会社 粒子分取装置、粒子分取方法、プログラム及び粒子分取システム
JP6102783B2 (ja) 2014-02-14 2017-03-29 ソニー株式会社 粒子分取装置、粒子分取方法及びプログラム
JP2017511884A (ja) * 2014-03-06 2017-04-27 ライフ テクノロジーズ コーポレーション 流体装置システムを診断するため及びフローサイトメータ用のデータ処理設定を決定するためのシステム及び方法
JP6657625B2 (ja) 2014-09-05 2020-03-04 ソニー株式会社 液滴分取装置、液滴分取方法及びプログラム
JP2016080563A (ja) * 2014-10-20 2016-05-16 日本光電工業株式会社 分析システム、分析装置
JP6282969B2 (ja) * 2014-10-31 2018-02-21 日本光電工業株式会社 フロー解析装置、フローサイトメータ、及びフロー解析方法
EP3215855A4 (en) * 2014-11-06 2018-03-28 Bio-Rad Laboratories, Inc. Droplet velocity detection
CN108139312B (zh) 2015-10-19 2021-02-05 索尼公司 图像处理设备、微粒分选设备和图像处理方法
US10732094B2 (en) * 2016-10-03 2020-08-04 Becton, Dickinson And Company Methods and systems for determining a drop delay of a flow stream in a flow cytometer
JP2017122734A (ja) * 2017-03-02 2017-07-13 ソニー株式会社 粒子分取装置、粒子分取方法及びプログラム
WO2018235383A1 (ja) * 2017-06-21 2018-12-27 ソニー株式会社 サンプル送液装置、フローサイトメータ、およびサンプル送液方法
JP2020051936A (ja) * 2018-09-27 2020-04-02 ソニー株式会社 微小粒子測定装置及び微小粒子測定方法
CN109799247B (zh) * 2019-02-15 2023-12-22 河北大学 基于微波传输时间的两相流相含率检测装置及方法
US11169899B2 (en) * 2019-04-15 2021-11-09 Toyota Motor Engineering & Manufacturing North America, Inc. Mitigating data offsets for machine learning
WO2020241226A1 (ja) * 2019-05-30 2020-12-03 ソニーセミコンダクタソリューションズ株式会社 光学測定装置及び光学測定システム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5150313A (en) * 1990-04-12 1992-09-22 Regents Of The University Of California Parallel pulse processing and data acquisition for high speed, low error flow cytometry
US6248590B1 (en) * 1998-02-27 2001-06-19 Cytomation, Inc. Method and apparatus for flow cytometry
US7653509B2 (en) * 2007-08-29 2010-01-26 Verity Software House Probability state models
CN102077174A (zh) * 2008-05-02 2011-05-25 英格朗公司 流式细胞计数仪远程监视系统

Family Cites Families (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3826364A (en) * 1972-05-22 1974-07-30 Univ Leland Stanford Junior Particle sorting method and apparatus
US4070617A (en) * 1974-05-08 1978-01-24 Max-Planck-Gesellschaft Zur Forderung Der Wissenschaften E.V. Device for controlling the particle flow in an apparatus for measuring the properties of particles suspended in liquid
US3963606A (en) * 1974-06-03 1976-06-15 Coulter Electronics, Inc. Semi-automatic adjusting delay for an electronic particle separator
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
JPS5319891A (en) * 1976-06-10 1978-02-23 Coulter Electronics Method and apparatus for folling drop formation and separation
US4184766A (en) * 1976-10-28 1980-01-22 Coulter Electronics, Inc. Method and apparatus for correlating measurements of tandem sensing zones
US4318480A (en) * 1979-08-20 1982-03-09 Ortho Diagnostics, Inc. Method and apparatus for positioning the point of droplet formation in the jetting fluid of an electrostatic sorting device
US4318482A (en) * 1979-08-20 1982-03-09 Ortho Diagnostics, Inc. Method for measuring the velocity of a perturbed jetting fluid in an electrostatic particle sorting system
US4318483A (en) 1979-08-20 1982-03-09 Ortho Diagnostics, Inc. Automatic relative droplet charging time delay system for an electrostatic particle sorting system using a relatively moveable stream surface sensing system
US4325483A (en) * 1979-08-20 1982-04-20 Ortho Diagnostics, Inc. Method for detecting and controlling flow rates of the droplet forming stream of an electrostatic particle sorting apparatus
US4538733A (en) * 1983-10-14 1985-09-03 Becton, Dickinson And Company Particle sorter with neutralized collection wells and method of using same
JPS63307332A (ja) * 1987-06-08 1988-12-15 Hitachi Ltd 光計測式生体細胞分析装置
US5030002A (en) 1989-08-11 1991-07-09 Becton, Dickinson And Company Method and apparatus for sorting particles with a moving catcher tube
EP0422616B1 (en) * 1989-10-11 1996-02-07 Canon Kabushiki Kaisha Apparatus for and method of fractionating particle in particle-suspended liquid in conformity with the properties thereof
US5548395A (en) * 1991-09-20 1996-08-20 Toa Medical Electronics Co., Ltd. Particle analyzer
US5483469A (en) * 1993-08-02 1996-01-09 The Regents Of The University Of California Multiple sort flow cytometer
US5700692A (en) * 1994-09-27 1997-12-23 Becton Dickinson And Company Flow sorter with video-regulated droplet spacing
US5602349A (en) * 1994-10-14 1997-02-11 The University Of Washington Sample introduction system for a flow cytometer
US5602039A (en) * 1994-10-14 1997-02-11 The University Of Washington Flow cytometer jet monitor system
DE19549015C1 (de) * 1995-12-28 1997-04-03 Siemens Ag Verfahren und Anordnung zur Überwachung eines abreißenden Flüssigkeitstrahls
US6079836A (en) * 1998-07-20 2000-06-27 Coulter International Corp. Flow cytometer droplet break-off location adjustment mechanism
US6372506B1 (en) * 1999-07-02 2002-04-16 Becton, Dickinson And Company Apparatus and method for verifying drop delay in a flow cytometer
WO2001063270A1 (en) * 2000-02-23 2001-08-30 Caliper Technologies, Inc. Multi-reservoir pressure control system
US6608680B2 (en) 2000-08-25 2003-08-19 Amnis Corporation TDI imaging system for kinetic studies
CN102830053A (zh) 2003-03-28 2012-12-19 英格朗公司 用于分拣颗粒和提供性别分拣的动物精子的设备、方法和程序
WO2004104178A2 (en) * 2003-05-15 2004-12-02 Xy, Inc. Efficient haploid cell sorting for flow cytometer systems
KR101149880B1 (ko) 2003-08-13 2012-05-24 루미넥스 코포레이션 유세포 분석기식 측정 시스템의 하나 이상의 파라미터의 제어 방법
US7274316B2 (en) 2004-11-17 2007-09-25 Luminex Corporation System and method for managing data from a flow analyzer
US7709544B2 (en) 2005-10-25 2010-05-04 Massachusetts Institute Of Technology Microstructure synthesis by flow lithography and polymerization
JP5032792B2 (ja) * 2006-05-22 2012-09-26 浜松ホトニクス株式会社 細胞選別装置
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
US7691636B2 (en) * 2007-05-23 2010-04-06 Beckman Coulter, Inc. Method and apparatus for compensating for variations in particle trajectories in electrostatic sorter for flowcell cytometer
CA2699319A1 (en) * 2007-09-17 2009-03-26 Luminex Corporation Systems, storage mediums, and methods for identifying particles in flow
JP4990746B2 (ja) * 2007-12-14 2012-08-01 ベイバイオサイエンス株式会社 液体フローに含まれる生物学的粒子を分別する装置ならびにその方法
JP5446563B2 (ja) * 2009-08-06 2014-03-19 ソニー株式会社 微小粒子分取装置、および該微小粒子分取装置を用いたフローサイトメーター
CN102792147B (zh) * 2010-03-10 2015-05-13 贝克曼考尔特公司 在颗粒分析器中产生脉冲参数
EP2732395B1 (en) 2011-08-25 2018-04-04 Sony Corporation Compensation of motion-related error in a stream of moving mirco-entities

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5150313A (en) * 1990-04-12 1992-09-22 Regents Of The University Of California Parallel pulse processing and data acquisition for high speed, low error flow cytometry
US6248590B1 (en) * 1998-02-27 2001-06-19 Cytomation, Inc. Method and apparatus for flow cytometry
US7653509B2 (en) * 2007-08-29 2010-01-26 Verity Software House Probability state models
CN102077174A (zh) * 2008-05-02 2011-05-25 英格朗公司 流式细胞计数仪远程监视系统

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015154288A1 (zh) * 2014-04-10 2015-10-15 深圳迈瑞生物医疗电子股份有限公司 一种进行自动补偿的方法、装置及相应的流式细胞分析仪
US9970857B2 (en) 2014-04-10 2018-05-15 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Automatic compensation method, device, and corresponding flow cytometer
CN108548817A (zh) * 2018-03-26 2018-09-18 清华大学 多相浮射流实验发生装置及油滴气泡阴影图像处理方法
CN108548817B (zh) * 2018-03-26 2021-08-03 清华大学 多相浮射流实验发生装置及油滴气泡阴影图像处理方法
CN116957992A (zh) * 2023-09-20 2023-10-27 南京木木西里科技有限公司 基于特征追踪的实时显微图像防抖方法
CN116957992B (zh) * 2023-09-20 2024-01-05 南京木木西里科技有限公司 基于特征追踪的实时显微图像防抖方法

Also Published As

Publication number Publication date
US9897530B2 (en) 2018-02-20
CN103180858B (zh) 2016-01-20
WO2013028947A1 (en) 2013-02-28
US20140220621A1 (en) 2014-08-07
US20140212917A1 (en) 2014-07-31
EP2732395B1 (en) 2018-04-04
WO2013028948A1 (en) 2013-02-28
EP2724144A4 (en) 2015-03-11
EP2732395A4 (en) 2015-03-11
EP2724144B1 (en) 2018-04-11
CN103180713B (zh) 2015-07-15
EP2732395A1 (en) 2014-05-21
CN103180713A (zh) 2013-06-26
EP2724144A1 (en) 2014-04-30

Similar Documents

Publication Publication Date Title
CN103180858B (zh) 移动微实体流中的运动相关误差的补偿
CN103460018B (zh) 颗粒分选设备和方法
EP2545357B1 (en) Generating pulse parameters in a particle analyzer
US7569789B2 (en) Cantilevered coaxial flow injector apparatus and method for sorting particles
US8263407B2 (en) Method for non-contact particle manipulation and control of particle spacing along an axis
EP2950079A1 (en) Fine particle fractionation device, fine particle fractionation method and program
CN104096687A (zh) 微粒分选装置和微粒分选方法
CN113316713A (zh) 颗粒分析仪的自适应分选
Poole et al. Second order macroscopic traffic flow model validation using automatic differentiation with resilient backpropagation and particle swarm optimisation algorithms
EP4024027A1 (en) Sampling control device, particle sampling device and particle sampling system in which sampling control device is used, method for controlling sampling, and control program
Sano et al. Microfluidic mechanotyping of a single cell with two consecutive constrictions of different sizes and an electrical detection system
Van Den Engh High‐speed cell sorting
Shellard et al. Frictiotaxis underlies adhesion-independent durotaxis
Welf et al. A unified role for membrane-cortex detachment during cell protrusion initiation
CN108474727A (zh) 用于分析在待检测液体中的颗粒的流动室
CN112903569A (zh) 用于计算液滴延迟时间的系统和方法以及分选装置
Bause et al. Reweighting non-equilibrium steady-state dynamics along collective variables
US20240042491A1 (en) Systems and methods for particle sorting with automated adjustment of operational parameters
EP4130710A2 (en) Label free cell sorting
US10690583B2 (en) Microparticle measuring device and microparticle measuring method
Spagnolo et al. Trajectory Analysis in Single Particle Tracking: From Mean Square Displacement to Machine Learning Approaches
Perez‐Gonzalez et al. Evaluation of Sort Recovery via Rmax
Kimmerling A toolset for linking phenotype and gene expression at the single-cell level
Beer et al. Passive chip-based droplet sorting

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant