CN107532989B - 用于分析颗粒的方法 - Google Patents
用于分析颗粒的方法 Download PDFInfo
- Publication number
- CN107532989B CN107532989B CN201680022598.1A CN201680022598A CN107532989B CN 107532989 B CN107532989 B CN 107532989B CN 201680022598 A CN201680022598 A CN 201680022598A CN 107532989 B CN107532989 B CN 107532989B
- Authority
- CN
- China
- Prior art keywords
- particles
- sample
- matrix array
- image
- photodetector
- 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.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 151
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims abstract description 65
- 239000002131 composite material Substances 0.000 claims description 77
- 210000003743 erythrocyte Anatomy 0.000 claims description 38
- 238000001514 detection method Methods 0.000 claims description 35
- 210000000265 leukocyte Anatomy 0.000 claims description 24
- 230000014509 gene expression Effects 0.000 claims description 23
- 210000001772 blood platelet Anatomy 0.000 claims description 22
- 230000003287 optical effect Effects 0.000 claims description 16
- 210000004027 cell Anatomy 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 5
- 230000001427 coherent effect Effects 0.000 claims description 4
- 230000001902 propagating effect Effects 0.000 claims description 4
- 210000000601 blood cell Anatomy 0.000 claims description 3
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 239000000523 sample Substances 0.000 description 78
- 230000000295 complement effect Effects 0.000 description 24
- 230000006870 function Effects 0.000 description 21
- 210000004369 blood Anatomy 0.000 description 14
- 239000008280 blood Substances 0.000 description 14
- 239000013598 vector Substances 0.000 description 14
- 210000002381 plasma Anatomy 0.000 description 8
- 210000001124 body fluid Anatomy 0.000 description 7
- 239000010839 body fluid Substances 0.000 description 7
- LOKCTEFSRHRXRJ-UHFFFAOYSA-I dipotassium trisodium dihydrogen phosphate hydrogen phosphate dichloride Chemical compound P(=O)(O)(O)[O-].[K+].P(=O)(O)([O-])[O-].[Na+].[Na+].[Cl-].[K+].[Cl-].[Na+] LOKCTEFSRHRXRJ-UHFFFAOYSA-I 0.000 description 6
- 239000002953 phosphate buffered saline Substances 0.000 description 6
- 239000000872 buffer Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 239000011325 microbead Substances 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000011144 upstream manufacturing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- 239000007864 aqueous solution Substances 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 229920001817 Agar Polymers 0.000 description 2
- 239000008272 agar Substances 0.000 description 2
- 239000013060 biological fluid Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 229920002307 Dextran Polymers 0.000 description 1
- 240000004808 Saccharomyces cerevisiae Species 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000010170 biological method Methods 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 238000004820 blood count Methods 0.000 description 1
- 238000009534 blood test Methods 0.000 description 1
- 210000001185 bone marrow Anatomy 0.000 description 1
- 239000005018 casein Substances 0.000 description 1
- BECPQYXYKAMYBN-UHFFFAOYSA-N casein, tech. Chemical compound NCCCCC(C(O)=O)N=C(O)C(CC(O)=O)N=C(O)C(CCC(O)=N)N=C(O)C(CC(C)C)N=C(O)C(CCC(O)=O)N=C(O)C(CC(O)=O)N=C(O)C(CCC(O)=O)N=C(O)C(C(C)O)N=C(O)C(CCC(O)=N)N=C(O)C(CCC(O)=N)N=C(O)C(CCC(O)=N)N=C(O)C(CCC(O)=O)N=C(O)C(CCC(O)=O)N=C(O)C(COP(O)(O)=O)N=C(O)C(CCC(O)=N)N=C(O)C(N)CC1=CC=CC=C1 BECPQYXYKAMYBN-UHFFFAOYSA-N 0.000 description 1
- 235000021240 caseins Nutrition 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 210000001808 exosome Anatomy 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 210000004623 platelet-rich plasma Anatomy 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000001376 precipitating effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 208000010110 spontaneous platelet aggregation Diseases 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000003053 toxin Substances 0.000 description 1
- 231100000765 toxin Toxicity 0.000 description 1
- 108700012359 toxins Proteins 0.000 description 1
- 210000002700 urine Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/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/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
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
- G01N21/453—Holographic interferometry
-
- 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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4788—Diffraction
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
-
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
- G06V20/698—Matching; Classification
-
- 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/1454—Optical arrangements using phase shift or interference, e.g. for improving contrast
-
- 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/1488—Methods for deciding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/49—Blood
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
- G03H2001/0447—In-line recording arrangement
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
- G03H2001/0883—Reconstruction aspect, e.g. numerical focusing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Signal Processing (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Multimedia (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明是一种用于识别包含在样本中的颗粒的方法,其包括使用光源照射样本,所述光源产生向样本传播的入射光波,然后使用矩阵阵列光电检测器获取样本的图像,样本被放置在所述光源和光电检测器之间,使得矩阵阵列光电检测器暴露于由入射光波和由每个颗粒产生的衍射波之间的干涉引起的光波。所述方法的特征在于,其包括对由光电检测器获取的图像应用数字重建算法,以便在距检测器多个距离处估计到达检测器的光波的特征量。随距离变化的特征量允许识别颗粒。
Description
技术领域
本发明涉及对存在于液体中,特别是体液、例如血液中的颗粒进行计数和识别的领域。
背景技术
体液、特别是血液中可能含有颗粒、例如细胞,知道其类型和数量是有用的。
例如,对于血液,全血细胞计数或全血检查通常是在医学实验室进行的检验。这种类型的检验允许识别血液的主要成分(特别是红细胞、白细胞或血小板)及其数量。这些检查通常使用高性能自动化装置进行,但是正在研究比较便宜同时允许获得相当性能水平的更简单的方法。
研究所寻求的途径之一就是使用简单的光学方法,例如无透镜成像。2000年底以来,通过无透镜成像观察生物颗粒已经有一定的发展。该技术在于将样本放置在光源和矩阵阵列光电检测器或图像传感器之间。由光电检测器捕获的图像通过由光源产生的入射波和由构成样本的颗粒衍射的波之间的干涉形成。该图像通常称为“全息图”。因此,对于每个颗粒,可以在传感器上记录其特定的衍射图案。适用于生物样本的这种技术已经在文献WO2008090330中描述。然后可以通过将其产生的衍射图案和事先建立的衍射图案相比较并且与已知颗粒相对应来对每个颗粒进行简单分析。然而,随着颗粒浓度增加,这种方法可能达到极限。
具体地,仅基于由图像传感器检测到的衍射图案对颗粒进行计数和识别随着样本中颗粒浓度的增加而达到一定的极限。尤其是,当样本是血液且颗粒是每μl超过100000个颗粒的红细胞时,计数不再可靠,如在2010年6月1日,Anal Chem,Seo Sungjyu的“High-throughput lensfree blood analysis on a chip”的出版物中所报道的。可以应用称为数字全息重建技术的数学技术以构建样本中存在的每个颗粒的所谓复合图像。这种类型的技术在于将光波反向传播到颗粒位于其中的物体平面,所述物体平面位于距成像器的已知距离处。上述出版物表明,这样的全息重建允许对样本中的高浓度红细胞进行计数。该出版物报道,在重建的复合图像中,事先经过标记的白细胞具有与红细胞不同的特征。
文献US2014/0327944还描述了基于全息图、通过对由图像传感器获得的全息图与模拟的全息图的库进行比较来对颗粒、例如血液颗粒进行分类的方法。然而,这种方法也有相同的限制,即当粒子密度高时变得难以实施。
在文献US2012/0148141和WO2014/012031中也描述了允许重建细胞(在该情况下是***)的复合图像的方法。然而,这些方法允许从重建的复合图像估计所述细胞的特性及其路径。文献US2009/0290156也是如此,其描述了基于所述颗粒的复合图像来分类颗粒以及跟踪所述颗粒的路径。样本的复合图像可能不足以识别颗粒。
寻求可以应用于颗粒浓度高的样本的、用于识别颗粒、特别是血细胞的方法。此外,该方法还具有大的观察范围并且易于实施,同时特别是避免需要事先标记颗粒。此外,该方法必须允许容易在体液中发现的颗粒被可靠地区分,所述颗粒特别是红细胞、白细胞和血小板。
此外,寻求一种不需要精确了解颗粒和光电检测器之间的距离的方法。
发明内容
为应对该问题,本发明提供了一种用于识别样本中存在的颗粒的方法,所述样本是例如生物液体、比如血液的样本,该方法包括以下步骤:
-使用光源照射所述样本,所述光源产生沿着传播轴线向样本传播的入射光波;
-使用矩阵阵列光电检测器获取样本的图像,所述样本置于所述光源和所述光电检测器之间使得矩阵阵列光电检测器暴露于光波,所述光波包含入射光波与由每个颗粒产生的衍射波之间的干涉;
所述方法特征在于,其还包括以下步骤:
-确定所述颗粒在与矩阵阵列光电检测器所在平面平行的平面中的位置;
-对所述获取的图像应用数字重建算法,以便在距离所述矩阵阵列光电检测器的多个重建距离处估计矩阵阵列光电检测器所暴露于的所述光波的至少一个特征量;
-确定表示所述特征量沿着与所述传播轴线平行并通过所述位置的轴线随所述重建距离的变化的曲线;
-根据所述曲线识别颗粒。
通过应用数字重建算法意味着通常以卷积乘积的形式将传播算子应用于图像。
特征量可以通过在每个重建距离处估计矩阵阵列光电检测器所暴露于的光波的复数表达式来获得。
特征量可以由所述复数表达式的模来确定,在这种情况下,其表示检测器暴露于的所述光波的幅度。
特征量可以根据所述复数表达式的幅角来确定,在这种情况下,其表示矩阵阵列光电检测器所暴露于的所述光波的相位。
根据一个实施方案,所述方法包括:
-通过对由矩阵阵列光电检测器获取的图像应用数字重建算法来确定称为参考复合图像的复合图像;
-在距离所述矩阵阵列光电检测器的多个重建距离处,根据所述参考复合图像估计矩阵阵列光电检测器所暴露于的光波的至少一个特征量。
然后,所述方法可以包括:
-将传播算子应用于参考复合图像,以便针对距重建平面或矩阵阵列光电探测器所在的平面的多个距离,计算所谓二次复合图像;
-从每个二次复合图像确定在所述距离的每一处的特征量。
所述参考复合图像可以是在远离样本平面的重建平面中形成的复合图像。它也可以是在检测平面中形成的复合图像。
可以通过将所述特征量的变化与在学习阶段确定的参考曲线进行比较来实现所述识别。
可以使用由光电检测器获取的图像或使用光电检测器所暴露于的光波的复数表达式来确定每个颗粒在平行于矩阵阵列光电检测器平面的平面中的位置。
光源优选为空间相干光源,例如发光二极管,在这种情况下,优选在光源和样本之间放置空间滤波器。所述光源可以是时间相干光源,例如激光二极管。
矩阵阵列光电检测器包括能够收集光电检测器暴露于的波的像素的矩阵阵列。像素和样本之间的距离可以为50μm至2cm,优选为100μm至5mm。优选地,不将样本放置成与光电检测器的像素直接接触。
优选地,在样本和矩阵阵列光电检测器之间不放置放大光学元件。
特别地,样本可以包括血细胞。在这种情况下,可以在白细胞、红细胞或血小板的细胞系中识别颗粒。
本发明的另一主题是一种用于识别颗粒的装置,所述颗粒包含在样本中,所述装置包括:
-光源,其被布置成在所述样本的方向沿着传播轴线产生入射光波;
-固定器,用于将样本固定在所述光源和矩阵阵列光电检测器之间;
-矩阵阵列光电检测器,其被布置成获取样本的图像,能够暴露于由所述入射光波与由所述颗粒形成的衍射波之间的干涉引起的光波;
其特征在于,该装置包括被配置为实施以下操作的处理器,例如微处理器或电子计算机:
-确定所述颗粒在与矩阵阵列光电检测器的平面平行的平面中的位置;
-对所述获取的图像应用数字重建算法,以便在距离所述矩阵阵列光电检测器的多个重建距离处估计矩阵阵列光电检测器所暴露于的所述光波的至少一个特征量;
-确定表示所述特征量沿着与所述传播轴线平行并通过所述位置的轴线随所述重建距离的变化的曲线;
-根据所述曲线识别颗粒。
优选地,所述装置在矩阵阵列光电检测器和被分析的样本之间不包括放大光学元件。
处理器可以包括或连接到可编程的存储器,其包括允许执行上述步骤的一系列指令。
特别地,其能够:
-在每个重建距离处确定检测器所暴露于的光辐射的复数表达式;
-通过确定所述复数幅值的模或幅角,估计在每个重建距离处的所述特征量。
附图说明
图1示出根据本发明的一个实施方案的装置。
图2A示出由矩阵阵列光电检测器获取的图像。
图2B示出根据第一实施例,对于各种类型的颗粒,光电检测器所暴露于的光波的特征量、称为互补幅度随相对于光电检测器的距离变化的曲线。
图3A示出根据该第一实施例,对于各种白细胞,光电检测器所暴露于的光波的特征量、称为互补幅度随相对于光电检测器的距离变化的曲线。
图3B示出根据该第一实施例,对于各种红细胞,光电检测器所暴露于的光波的特征量、称为互补幅度随相对于光电检测器的距离变化的曲线。
图4A示出根据该第一实施例,由光电检测器获取的图像的感兴趣区域,其以血小板聚集群为中心。
图4B示出对于形成图4A示出的集群的部分的各种血小板,光电检测器所暴露于的光波的特征量、称为互补幅度的曲线。
图5A示出根据第二实施例,对于各种类型的颗粒,光电检测器所暴露于的光波的特征量、称为互补幅度随相对于光电检测器的距离变化的曲线。
图5B示出根据第二实施例,对于各种类型的颗粒,光电检测器所暴露于的光波的相位随相对于光电检测器的距离z变化的曲线。
图6A示出根据第三实施例,对于各种类型的颗粒,光电检测器所暴露于的光波的特征量、称为互补幅度随相对于光电检测器的距离变化的曲线。
图6B示出根据第三实施例,对于各种类型的颗粒,光电检测器所暴露于的光波的相位随相对于光电检测器的距离变化的曲线。
图7A、7B和7C分别示出根据第一、第二和第三实施例,光电检测器所暴露于的光波的复合特征量随距离变化的曲线。
图8A、8B、8C和8D分别示出:
-允许在重建平面中计算的样本的复合图像、称为参考复合图像的方法;
-由矩阵阵列光电检测器获取的全息图;
-在图8A所示的方法的多次迭代之后重建的图像、称为参考复合图像的表示;
-基于由参考复合图像形成的二次复合图像获得的曲线。
图9A是由图像传感器获取的全息图,样本包括分散在水溶液中的红细胞。图9B和9C分别示出称为参考图像的复合图像的模和相位,该复合图像在重建平面中形成。图9D和9E是分别示出图像传感器所暴露于的光波的模和相位沿着穿过红细胞的传播轴线的变化的曲线。
具体实施方式
图1示出作为本发明的一个主题的装置的实施例。光源11能够沿着传播轴线Z在样本14的方向上产生被称为入射光波的光波12。样本14包括介质10,例如生物液体,其包括期望在预设类型的颗粒之间识别的颗粒1、2、3、4、5、...、9。
颗粒可以是细胞。特别地,当介质10是血液或包含血液的溶液时,颗粒可以是红细胞、白细胞或血小板。
颗粒也可以是有机或无机微珠,例如金属微珠或聚合物或玻璃的微珠,这种类型的微珠通常在实施生物方法时被使用。颗粒也可以是浸入介质10中的液滴,例如脂滴。其也可以是微生物、例如细菌或酵母,或外来体。通常,颗粒的尺寸有利地小于1mm,甚至小于500μm,优选为0.5μm至500μm的尺寸。因此,术语颗粒是指最初存在于被检测样本中的内源性颗粒和在分析之前添加到该样本中的外源性颗粒。
介质10最通常地是液体介质,特别是体液,但其也可以是琼脂、或空气、或液体的干残留物。
作为本发明的一个主题的方法允许识别每个观察到的颗粒。识别是将颗粒分类成预定类型的颗粒。其可以是确定颗粒在预设性质中的性质、或者确定颗粒在预设性质中的大小的问题。
光源和样本之间的距离Δ优选大于1cm。其优选为2至30cm。优选地,从样本的角度看,光源可以认为是点状的。这意味着,其直径(或其对角线)优选为小于样本和光源之间的距离的十分之一,更好的是小于样本和光源之间的距离的百分之一。因此,光以平面波的形式、或可以认为是平面波的波的形式到达样本。
光源11可以是点光源,或与隔板或空间滤波器(图1中未示出)相关联以便呈现点状。隔板的孔径通常为5μm至1mm,优选为50μm至500μm。
隔板可以由光纤代替,光纤的第一端面向光源放置,其第二端面向样本放置。在这种情况下,所述第二端可以比作点光源11。
样本14由包括基座15和盖13的室所限定。该室的侧壁没有示出。在所考虑的实施例中,所述室是Neubauer C-芯片流体室。基座15和盖13之间的距离为100μm。通常,室的沿着传播轴线Z的厚度小于几厘米,例如小于1cm,甚至小于1mm,例如为50μm至500μm。
光源11可以是时间相干的,但这不是必需的。
在该第一实施例中,光源是以450nm的波长发光的激光二极管。其位于距样本15cm的距离处。
样本14放置在光源11和矩阵阵列光电检测器或图像传感器16之间。后者优选平行地或基本上平行地延伸到包围样本的室的基座15。
表述“基本平行地”意味着两个元件可以不是严格平行的,几度、小于10°的角度公差是可接受的。
优选地,光源具有小的光谱宽度,例如光谱宽度小于100nm或甚至小于20nm,甚至优选小于5nm。表述“光谱宽度”表示在光源发射峰的半最大值处的全宽。
光电检测器16可以是矩阵阵列光电检测器,其包括CCD或CMOS像素的矩阵阵列。CMOS光电检测器是优选的,因为像素的尺寸较小,这允许更有利于获取的图像的空间分辨率。在该实施例中,检测器是参考MT9P031的12位APTINA传感器,其是RGB CMOS传感器,其像素间间距为2.2μm。光电检测器的有用面积是5.7×4.3mm2。
光电检测器位于优选垂直于入射光波12的传播轴线Z的检测平面P中。
优选地,光电检测器包括像素的矩阵阵列,其上方放置有透明保护窗。像素的矩阵阵列与保护窗之间的距离通常为几十μm至150或200μm。为了提高图像的空间分辨率,光电检测器的像素间间距优选小于3μm。
在该实施例中,颗粒1、2、...9与光电检测器16的像素的矩阵阵列之间的距离d等于1.5mm。然而,其可能根据所使用的流体室的厚度而波动。通常,无论哪种实施方案,颗粒与光电检测器的像素之间的距离d优选为50μm至2cm,优选为100μm至2mm。
注意到矩阵阵列光电检测器16和样本14之间不存在放大光学元件。这不妨碍聚焦微透镜可选地与光电检测器16的每个像素保持同一水平。
在该第一实施例中,样本是根据常规方法获得的白细胞富血浆,在6%右旋糖苷(Sigma Aldrich参考D4876)的存在下在Alsever溶液中沉淀红细胞后,收集血小板富血浆和白细胞富血浆。然后将所获得的血浆在生理pH的磷酸盐缓冲盐水(PBS)缓冲液中稀释。红细胞的消耗是不完全的,所获得的富血浆含有残留的红细胞。
颗粒可以分为多种类型的颗粒,特别是红细胞、白细胞或血小板。优选地,颗粒预先没有经过标记。
图2A示出由光电检测器16所获得的图像。该图像示出总衍射图案,其中可以看到基本衍射图案,每个基本衍射图案与相应的颗粒相关联。每个基本衍射图案包括中心的圆盘形区域,深和浅的同心圆围绕其交替延伸。这种基本图案允许选择待识别的颗粒,以及确定在检测平面P中所述颗粒的径向坐标(x,y)。这些坐标例如是与所述颗粒对应的基本衍射图案的中心。
每个基本衍射图案是由样本上游的源11产生的入射光波12与由颗粒的入射波的衍射所产生的波之间的干涉形成的。因此,光电检测器16暴露于由以下叠加形成的光波22:
-由样本14上游的源11发射的光波12;
-由样本14中存在的每个颗粒或其它衍射元素衍射的光波。
处理器20(例如微处理器)接收矩阵阵列光电检测器16的图像,并沿着传播轴线Z执行矩阵阵列光电检测器所暴露于的光波22的特征量的重建。微处理器20连接到能够存储用于实施本申请中所描述的计算步骤的指令的存储器23。其可以连接到屏幕25。特别地,在矩阵阵列光电检测器和观察到的样本之间形成重建。
处理器20能够执行存储在存储器中的一系列指令,以便实施识别方法的步骤。处理器可以是微处理器或能够处理矩阵阵列光电检测器传送的图像的任何其它电子计算机,以便执行本说明书中描述的一个或更多个步骤。
由矩阵阵列光电检测器获取的图像I,其图像如图2A所示,表示光波22的强度I(x,y)的空间分布,x和y是光电检测器的平面P中的坐标。
根据众所周知的数字全息重建原理,其在Ryle等人的“Digital in-lineholography of biological specimens”,Proc.of SPIE第6311卷(2006)中描述,可以通过确定由光电检测器测量的强度I(x,y)和传播算子h(x,y,z)的卷积乘积来重建在空间坐标(x,y,z)任何点处的、特别是在位于距光电检测器距离|z|的平面中的光波22的复数表达式U(x,y,z)。
-u(x,y,z)=abs[U(x,y,z)];
运算符abs和arg分别返回模和幅角。
特别地,传播算子的应用允许在光电检测器上游、距光电检测器距离|z|处估计复数表达式。由此,在光波22到达所述检测器之前,光波22的复合值被重建。然后,说到反向传播。如果坐标z=0被认为是检测平面P,则该反向传播通过应用传播算子h(x,y,-|z|)实施。术语“上游”和“下游”理解为相对于入射波12的传播方向。
如果I(x,y)=I(x,y,z=0)对应于由光电检测器测量的信号强度,则在检测平面P中,所测量的强度I(x,y)和的光波的复数表达式U(x,y)之间的关系表示为:I(x,y)=|U(x,y)|2。
在坐标(x,y,z)处,光波(22)的复数表达式表示为:
其中:
-在由检测平面P所描绘的且包括样本14的半空间中,z<0;
-在由检测平面P所描绘且不包括样本14的半空间中,z>0。
在由检测平面P所描绘的且包括样本14的半空间中,光波的复数表达式也可以写为:
优选地,在全息重建之前,将数学预处理事先应用于所测量的强度I(x,y)。这允许尤其是通过减少当应用传播算子时所产生的伪影的数量来改善结果的质量。
其中
-I(x,y)=在坐标(x,y)处通过光电检测器测量的强度;
-Average(I)=在图像I的感兴趣区域中测量的强度的平均值,所述区域包括所述坐标(x,y)。该感兴趣区域可以对应于通过光电检测器形成的整个图像。
该预处理将入射光波(12)的强度等同于所测量的强度的归一化,入射光波(12)的强度通过算子Average(I)估计。
数字重建可以特别地基于菲涅耳衍射模型。在该实施例中,传播算子是菲涅耳-亥姆霍兹函数,使得:
其中λ是波长。
其中
-x′和y′是光电检测器平面中的坐标;
-x和y是重建平面中的坐标,重建平面位于距光电检测器的距离|z|处;
-z是重建图像沿着入射光波(12)的传播轴线Z的坐标。
由复数表达式U(x,y,z)的值,可以提取源11发射的入射光波12的、由颗粒(1、2、...、9)引起的衍射产生的光波22的特征量。如上文所描述,可以求幅度u(x,y,z)或相位的值,但也可以求幅度或相位的任何函数的值。
从每个重建复数表达式U(x,y,z),可以形成:
-在平行于检测器的平面、距离检测器的平面|z|处的平面中,波22的幅度的图像uz,其中uz(x,y)=abs[U(x,y,z)];
在该第一实施例中,
-颗粒1至4:白细胞,由首字母缩写词WBC标示;
-颗粒5和6:红细胞,由首字母缩写词RBC标示;
-颗粒7至9:血小板,由字母PLT标示。
重建距离|z|在zmin=1000和zmax=1500μm之间变化。
与这些操作类似地,在显微镜下观察每个颗粒(1、…、9),在显微镜下观察所得作为参考测量,允许获得无可争议的识别。
在该实施例中所研究的样本中,可以看出:
与现有技术相比,不是通过从距样本预定距离处进行全息重建形成颗粒的复合图像,而是沿着入射波的传播方向在距光电检测器多个距离处重建由具有入射波12的颗粒衍射所产生的波22的特征。所获得的信息更加丰富,并允许在各种类型的颗粒之间的清晰分类。
图3A和3B示出检测器所暴露于的波22的互补幅度曲线所述曲线分别是针对50个白细胞WBC和240个红细胞RBC获取的。曲线的重复性足以允许基于颗粒的鲁棒分类。这些曲线在与上述实施例类似的实验条件下获得的。
用于获得图3A所示测量值的样本是类似于参照图2A和2B所描述的样本的富血浆。
用于获得图3B所示测量值的样本包括在如上所述的磷酸盐缓冲盐水(PBS)缓冲液中以1/400的稀释因子稀释的全部血液。
图4A示出通过光电检测器16获取的图像I的感兴趣区域。其示出待识别聚集群的每个血小板的坐标(x101,y101)、(x102,y102)、(x103,y103)、(x104,y104)。
在与第一实施例类似的实验条件下获得这些曲线。可以看出,无论血小板聚集或不聚集,曲线是相似的,并且限于基线BL,在间隔BL±ε内。因此,利用作为本发明一个主题的识别方法,即使血小板聚集,其也能被正确地识别。
在第二实施例中,光源11是耦接到以波长λ=485nm为中心的485-DF-22Omega光学滤波器的白色发光二极管,并且其半极大处全宽度为22nm。光源和检测器之间的距离Δ等于8cm。样本是如上所述的富血浆。
与前述实施例一样,通过在显微镜下观察来证实观察到的颗粒的性质。
因此,对于每个检测到的颗粒n,其在平行于检测器平面的平面中的位置是(xn,yn),可以建立曲线其表示在多个重建高度z处,检测器暴露于的波22的如上所述的互补幅度并且可以使用该曲线来分类所述颗粒是红细胞RBC、白细胞WBC或血小板PLT。
因此用除了激光源之外的光源进行分类是可能的。
图5B示出:
因此,对于每个检测到的颗粒n,其在与检测器平面平行的平面中的位置是(xn,yn),可以建立表示在多个重建高度z处、检测器所暴露于的辐射的相位变化的曲线并且可以使用该曲线来分类颗粒是是红细胞、白细胞或血小板。
在第三实施例中,光源11是耦接到以波长λ=610nm为中心的610-DF-20Omega光学滤波器的白色发光二极管,且其半极大处全宽度为20nm,所述光源置于距样本距离Δ=8cm处。操作模式和样本类似于前述实施例中的那些。
因此,如上所述,对于每个检测到的颗粒n,其在平行于检测器平面的平面中的位置是(xn,yn),可以建立曲线其表示在多个重建高度z处,检测器暴露于的辐射的互补幅度变化,并且可以使用该曲线来分类所述颗粒是红细胞RBC、白细胞WBC或血小板PLT。
图6B示出:
也可以使用所谓的复合光学参数、符号k,以比例的形式组合互补幅度和相位,例如:
图7A、7B和7C示出在第一实施例(405nm激光源)、第二实施例(白色LED光源与以λ=485nm为中心的滤波器的组合)和第三实施例(白色LED光源与以λ=610nm为中心的滤波器的组合)的配置中,沿着传播轴线Z的所述复合光学参数的变化。在每个配置中,观察到的样本是如上所述的富血浆。
在每个图中,经分析的颗粒是3个白细胞(WBC)和1个红细胞(RBC)。
可以看出,无论那种源,白细胞的曲线k(z)的波动大于红细胞的曲线k(z)的波动。
特别地,可以确定第一复合阈值kthreshold 1和第二复合阈值kthreshold 2,使得当曲线k(z)保持在低于第一复合阈值kthreshold 1且高于第二复合阈值kthreshold 2时,经分析的颗粒是红细胞。当曲线与这些阈值中的一个交叉,所检查的颗粒被识别为白细胞。
对于由矩阵阵列光电检测器16所获取的图像I或全息图,应用数字传播算子h可能具有一定的限制,因为所获取的图像不包括与相位相关联的信息。因此,在建立曲线前,优选地获得与光电检测器16所暴露于的光波22的相位相关联的信息。该与相位相关联的信息可以使用现有技术所描述的方法通过重建样本14的复合图像Uz获得,以便获得矩阵阵列光电检测器16的平面P中的光波22的幅度和相位的估计,或者获得位于距矩阵阵列光电检测器16为距离|z|的重建平面Pz中的光波22的幅度和相位的估计。本发明人已开发了基于参考复合图像的计算的方法,参考图8A描述了该方法。该方法包括以下步骤:
-用矩阵阵列光电检测器16获取样本14的图像I,该图像形成全息图(步骤100)。
-在重建平面Pz或检测平面P中计算称为样本14的参考图像Uref的复合图像,该参考复合图像包括关于矩阵阵列光电检测器16所暴露于的光波22的相位和幅度的信息;该步骤通过将如上所述的传播算子h应用于所获取的图像I(步骤110至170)来执行。该复合图像称为参考图像,因为基于其形成表征颗粒的曲线。
-使用参考复合图像Uref或由光电检测器16所获取的图像I来选择颗粒在检测平面或平行于检测平面的平面中的径向位置(x,y)(步骤180)。
-沿着传播轴线Z将传播算子h应用于参考复合图像Uref以便计算称为二级图像的复合图像Uref,z(步骤185)。
-由每个二级复合图像Uref,z,在预先选择的颗粒的径向位置(x,y)处以及在距重建平面Pz(或检测平面P)的多个距离处估计光波22的特征量,然后形成表示所述特征量沿着传播轴线Z变化的曲线(步骤190)。
-根据所述曲线表征颗粒。如上所述,可以通过将所获得的曲线与在校准阶段中使用标准样本所获得的标准曲线进行比较来实现该表征(步骤200)。
图8A中给出的算法详述如下,在某些步骤中获得的结果在图8B至8D中示出。步骤110至170是获得表示为Uref的参考复合图像的优选方式,该图像表示波22在重建平面Pz中的复数表达式的空间分布。本领域技术人员将理解,其他算法允许重建这样的复合图像,例如也可以想到使用现有技术提及的算法。
步骤100:图像获取。
在该步骤中,图像传感器16获取样本14的图像I,更准确地,由样本14发射的、图像传感器暴露于的光波22。这种图像或全息图示于图8B中。
该图像使用浸入盐水缓冲液中的包含红细胞的样本14产生,该样本包含在距离CMOS传感器1500μm的距离d的100μm厚的流动室中,如上文所述的装置中。
步骤110:初始化。
在该步骤中,由图像传感器16所获取的图像I定义样本14的初始图像该步骤是下面参照步骤120至180所描述的迭代算法的初始化,指数k表示每次迭代的阶。初始图像的模可以通过将平方根算子应用于由图像传感器获得的图像I来获得,在这种情况下在该实施例中,图像I通过表示入射到样本14上的光波12的强度的项来归一化。所述项例如可以是图像I的平均值的平方根,在这种情况下所获取的图像的每个像素I(x,y)被所述平均值除,使得
初始图像的相位或认为在每个像素(x,y)中是0,或预设为任意的值。特别地,初始图像直接来自于由矩阵阵列光电检测器16所获得的图像I。然后,后者不包括与通过样本14传输的光波22的相位相关联的信息,图像传感器16仅对该光波的强度敏感。
步骤120:传播。
符号*表示卷积运算符。下标–z表示传播在与传播轴线Z的方向相反的方向上进行,称为反向传播。
重建平面Pz是远离检测平面P的平面,优选地平行于检测平面。优选地,重建平面Pz是样本14位于的平面P14。特别地,在该平面中重建的图像允许获得通常高的空间分辨率。重建平面也可以是另一个平面,其位于距检测平面的非零距离处,优选平行于检测平面,例如位于矩阵阵列光电检测器16和样本14之间的平面。
步骤130:计算多个像素中的指标。
在该步骤中,计算与复合图像的多个像素(x,y)的每个像素相关联的量∈k(x,y),优选在这些像素中的每个中。在所计算的像素(x,y)中,该量取决于图像的值或其模。其也可以取决于该像素中的图像的维度导数,例如该图像的维度导数的模。
在步骤130中,在复合图像的多个像素中计算了量εk(x,y)。这些量可以形成向量Ek,其项是与每个像素(x,y)相关联的量εk(x,y)。在该步骤中,根据向量Ek的范数计算称为噪声指标的指标。通常,阶与范数有关,使得坐标(x1,x2,....xn,)的维数n的向量x的阶p的范数‖x‖p为:其中p≤0。
在这种情况下,使用1阶的范数,换句话说p=1。具体地,发明人估计1阶的或者低于或等于1阶的范数特别适用于如下所述的这种样本。
因此,εk=∑(x,y)εk(x,y)。
因为使用了1阶的或者低于或等于1阶的范数,噪声指标εk的值随着复合图像变得越来越类似样本而降低。具体地,在第一迭代中,粗略估计图像的每个像素(x,y)中的相位的值。然后,如现有技术中所提到的,将样本从检测平面P传播到重建平面Pz伴随着大量的重建噪声。该重建噪声采用重建图像中的波动的形式。由于这些波动,随着重建图像中重建噪声的贡献增加,如上定义的噪声指标εk的值增加。具体地,重建噪声引起的波动会增加该指标的值。
在第一迭代中,如上所述,仅关于光波22的强度的相关信息可得到,而关于其相位的相关信息不可得到。因此,由于没有关于检测平面P中的光波22的相位的相关信息,在重建平面Pz中重建的第一图像受到大量重建噪声的影响。因此,指标εk=1是高的。在随后的迭代中,算法在检测平面P中逐渐调整相位以便使指标εk逐渐最小化。
从光波22的强度及其相位的角度来看,检测平面中的图像表示在检测平面P中的光波22,步骤120至160旨在迭代地建立图像的每个像素的相位的值,其最小化指标εk,指标εk是通过将图像传播到重建平面Pz而获得的图像所获得的。
最小化算法可以是梯度下降算法或共轭梯度下降算法,后者在下文中描述。
步骤150:调整检测平面中的相位值。
步骤150旨在确定复合图像的每个像素的相位的值,以便在下一个迭代k+1中最小化指标εk+1,其由复合图像传播到重建平面Pz而产生。为此,建立相位向量其每一项是复合图像的像素(x,y)的相位该向量的维数是(Npix,1),其中Npix是所讨论的像素数。在每个迭代中,使用以下更新表达式来更新该向量:
-αk是整数,称为“步长”,表示距离;
该表达式可以表达为如下的向量形式:
其可以示为:
其中:
-pk-1是在上一次迭代中建立的方向向量;
-βk是应用于方向向量pk-1的比例因子。
其中Im是返回操作数虚部的算子,r′是检测平面中的坐标(x,y)。
比例因子βk可以表示为:
步长αk可以根据迭代变化,例如第一次迭代中的0.03至最后一次迭代的0.0005。
更新等式允许调整待获得的向量这导致复合图像的每个像素中的相位的迭代更新。然后用与每个像素相关联的相位的这些新值更新检测平面中的该复合图像将注意到的是,不修改复合图像的模,复合图像是由矩阵阵列光电检测器16所获取的图像确定的,使得
步骤160:重复或退出算法。
假设没有达到收敛标准,则步骤160在步骤150中更新的复合图像的基础上重复步骤120至160的新迭代以重复算法。收敛标准可以是迭代的预定数K、或者指标的梯度的最小值、或两个连续相位向量之间的被认为可忽略的差值。当达到收敛标准时,该估计被认为是在检测平面P或重建平面Pz中的样本的复合图像的正确估计。
步骤170:获得参考复合图像。
在最后一次迭代结束时,该方法可以包括将从最后一次迭代得到的复合图像传播到重建平面Pz,以获得参考复合图像替代地,参考复合图像Uref是由检测平面P中的最后一次迭代产生的复数图像当颗粒密度高时,该替代不太有利,因为检测平面P中的空间分辨率比在重建平面Pz中低,特别是当重建平面Pz对应于样本14所在的平面P14时。
步骤180:选择颗粒径向坐标。
在该步骤中,从参考图像例如从其模或从其相位的图像中选择颗粒的径向坐标(x,y)。如上所述,表述径向坐标表示在检测平面或重建平面中的坐标。也可以设想基于全息图I0或者基于最后一次迭代后的检测平面中获得的复合图像来进行选择。然而,当颗粒数增加时,优选地在重建平面中形成的图像上进行选择,这是由于其更好的空间分辨率,特别是当重建平面Pz对应于样本P14的平面时。在图8C中,示出颗粒的选择,其被虚线围绕。
步骤185:应用传播算子。
在该步骤185中,使用如上所述的传播算子h将参考复合图像Uref传播到多个重建距离,以便获得在距检测平面P或重建平面Pz不同距离处的多个所谓的二级复合图像Uref,z。因此,该步骤包括确定多个复合图像Uref,z,使得:
Uref,z=Uref*hz其中zmin≤z≤zmax。
值zmin和zmax是沿着轴线Z传播参考复合图像的最小和最大坐标。优选地,在样本14和图像传感器16之间的多个坐标z处重建复合图像。可以在样本14的任一侧上形成复合图像。
通过将全息重建算子h应用于参考图像Uref建立这些二次复合图像。参考图像Uref是在步骤120至160的迭代之后,正确描述图像传感器所暴露于的光波22的、特别是其相位的复合图像。因此,二次图像Uref,z形成沿着传播轴线Z传播光波22的良好的描述符。
步骤190:形成曲线
在该步骤中,从每个二级复合图像Uref,z确定如上所述的光波22的特征量,以便限定表示所述特征量沿着传播轴线Z的变化的曲线。特征量可以例如是模或相位,或其组合。图8D示出光波22的相位沿着传播轴线Z的变化。
步骤200:表征
然后可以由上一步骤中形成的曲线来表征颗粒。优选地,存在使用已知的标准样本在学习阶段形成的标准曲线的数据库。在标准曲线的基础上,通过比较或分类所形成的曲线来进行表征。
该实施方案已经在包括红细胞的样本上进行了测试。另一实施例如图9A至9E所示。在这些实施例中,样本包括稀释在包含稀释至1/400的磷酸盐缓冲盐水(PBS)缓冲液的水溶液中的红细胞。将样本14放置在100μm厚的流体室15中,该室位于与上述发光二极管的距离为8cm处,所述发光二极管的光谱带以450nm为中心。将样本放置在距上述CMOS图像传感器1.5mm的距离处。空间滤波器18的孔径尺寸为150μm。
图像形成如上所述应用传播算子h的参考图像Uref,以便获得沿着传播轴线Z的多个二级复合图像Aref,z。此外,在参考图像的模图像或相位图像中,识别出红细胞,红细胞在这些图像的每一个中被虚线围绕。提取该红细胞的径向坐标(x,y)。由二级复合图像Aref,z,形成表示模的曲线u(z)和表示到达图像传感器16的光波22的相位的曲线通过确定所述径向坐标处的相应二级图像的模和相位来获得曲线的每个点的值。图9D和9E分别示出如此选择的红细胞的模和相位的曲线。在坐标zmin=1000μm与zmax=2000μm之间确定曲线,z方向步长为5μm。重建平面位于距检测平面1380μm处,对应于图9D和9E中的横坐标76。
所描述的方法不限于血液,并且可以应用于其它体液,例如尿液、脑脊髓液、骨髓等。此外,该方法可以应用于非体液,特别是用于分析在水或任何其他水溶液中的污染物或毒素。
该方法还适用于放置在非液体介质中的颗粒的检测和识别,例如琼脂或体液的干残余物,例如引起在载玻片上大量沉积干血的血液涂片。在后者的情况下,颗粒通过干残留物或空气彼此分离。
此外,如上所述,颗粒可以是内源性的(例如血液颗粒)或外源性的(微珠、液滴)。
上述实施例提供了基于随重建距离变化的特征量的曲线的变化以及基于与预定阈值比较的简单识别标准。标准的有效性与在其中放置颗粒的介质以及样本的制备方法有关。其他标准可以适用于经历不同制备方法的颗粒。因此,对于给定类型的样本,可以在包括已知颗粒的标准样本上进行的学习阶段来定义识别标准。
此外,在不脱离本发明范围的情况下,可以实施更复杂和更鲁棒的其他分类方法。
Claims (15)
1.一种用于识别包含在样本(14)中的颗粒(1、2、…9、101…104)的方法,所述方法包括以下步骤:
a)使用光源(11)照射所述样本,所述光源产生沿着传播轴线(Z)向所述样本(14)传播的入射光波(12);
b)使用矩阵阵列光电检测器(16)获取所述样本的图像,所述样本放置于所述光源和所述矩阵阵列光电检测器之间,以使得所述矩阵阵列光电检测器暴露于光波(22),所述光波(22)包含入射光波(12)和由所述颗粒产生的衍射波之间的干涉;
所述方法的特征在于其还包括以下步骤:
c)确定所述颗粒在与所述矩阵阵列光电检测器(16)所在的平面(P)平行的平面中的位置(x,y);
f)通过将所述特征量的变化曲线与在学习阶段中确定的参考曲线进行比较来识别所述颗粒;
其中步骤f)包括确定各种类型的颗粒中的颗粒的类型。
2.根据权利要求1所述的方法,其中,通过在每个重建距离处估计所述矩阵阵列光电探测器所暴露于的光波(22)的复数表达式(U(x,y,z))来获得所述特征量。
3.根据权利要求2所述的方法,其中,所述特征量由所述复数表达式(U(x,y,z))的模或幅角确定。
4.根据权利要求1或2所述的方法,其中,使用由所述光电检测器获取的图像或使用所述光电检测器暴露于的光波(22)的复数表达式(U(x,y,z)) 来确定每个颗粒在与所述矩阵阵列光电检测器的平面平行的平面中的位置。
6.根据权利要求5所述的方法,包括:
-将传播算子(h)应用于所述参考复合图像(Uref),以便针对距所述重建平面或距所述矩阵阵列光电检测器所在平面的多个距离计算二次复合图像(Uref,z);
-根据每个二次 复合图像(Uref,z)确定在所述距离的每一距离处的特征量。
7.根据权利要求1或2所述的方法,其中所述光源是空间相干源。
8.根据权利要求1或2所述的方法,其中所述光源是发光二极管或激光二极管。
9.根据权利要求1或2所述的方法,其中在所述样本和所述矩阵阵列光电检测器之间不放置放大光学元件。
10.根据权利要求1或2所述的方法,其中所述样本包括血细胞。
11.根据权利要求1或2所述的方法,其中在白细胞或红细胞或血小板的细胞系中识别所述颗粒。
12.一种用于识别颗粒的装置,所述颗粒包含在样本(14)中,所述装置包括:
-光源(11),其被布置为在所述样本(14)的方向沿着传播轴线(Z) 产生入射光波(12);
-固定器,其用于将所述样本(14)固定在所述光源(11)和矩阵阵列光电检测器(16)之间;
-矩阵阵列光电检测器(16),其被布置为获取所述样本的图像,且被暴露于由所述入射光波(12)与由所述颗粒形成的衍射波之间的干涉引起的光波(22);
所述装置的特征在于其包括被配置为实施以下步骤的处理器(20):
-确定所述颗粒在平行于所述矩阵阵列光电检测器(16)所在平面(P)的平面中的位置(x,y);
-通过将所述特征量的变化曲线与在学习阶段中确定的参考曲线进行比较来识别颗粒;
其中,识别颗粒包括确定各种类型的颗粒中的颗粒的类型。
13.根据权利要求12所述的装置,其中所述装置在所述矩阵阵列光电检测器和所分析的样本之间不包括放大光学元件。
14.根据权利要求12或13所述的装置,其中,所述处理器被配置为:
-在每个重建距离处确定所述检测器所暴露于的光辐射的复数表达式(U(x,y,z));
-通过确定所述复数幅值(U(x,y,z))的模或幅角,估计在每个重建距离处的所述特征量。
15.根据权利要求12所述的装置,所述处理器是微处理器。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1552443A FR3034196B1 (fr) | 2015-03-24 | 2015-03-24 | Procede d'analyse de particules |
FR1552443 | 2015-03-24 | ||
PCT/FR2016/050643 WO2016151248A1 (fr) | 2015-03-24 | 2016-03-23 | Procédé d'analyse de particules |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107532989A CN107532989A (zh) | 2018-01-02 |
CN107532989B true CN107532989B (zh) | 2021-08-17 |
Family
ID=53366075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680022598.1A Active CN107532989B (zh) | 2015-03-24 | 2016-03-23 | 用于分析颗粒的方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US10379027B2 (zh) |
EP (1) | EP3274689B1 (zh) |
JP (1) | JP6727229B2 (zh) |
KR (1) | KR102479862B1 (zh) |
CN (1) | CN107532989B (zh) |
FR (1) | FR3034196B1 (zh) |
WO (1) | WO2016151248A1 (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3049347B1 (fr) * | 2016-03-23 | 2018-04-27 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede d’observation d’un echantillon par calcul d’une image complexe |
FR3049348B1 (fr) * | 2016-03-23 | 2023-08-11 | Commissariat Energie Atomique | Procede de caracterisation d’une particule dans un echantillon |
WO2017171762A1 (en) * | 2016-03-30 | 2017-10-05 | Siemens Healthcare Diagnostics Inc. | Systems, methods, and apparatus for processing platelet cell data |
DE102016215419A1 (de) * | 2016-08-17 | 2018-02-22 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Messanordnung und Verfahren zum Lenken und Detektieren von Partikeln |
FR3056749B1 (fr) * | 2016-09-28 | 2018-11-23 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede de numeration de leucocytes dans un echantillon |
FR3060746B1 (fr) * | 2016-12-21 | 2019-05-24 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede de numeration de particules dans un echantillon par imagerie sans lentille |
FR3063085B1 (fr) | 2017-02-17 | 2022-06-17 | Commissariat Energie Atomique | Procede optique de suivi de l'amplification in-vitro d'une sequence de nucleotides |
FR3066503B1 (fr) | 2017-05-22 | 2021-05-07 | Commissariat Energie Atomique | Procede d'analyse de microorganismes |
FR3071609B1 (fr) * | 2017-09-27 | 2019-10-04 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede de detection de microorganismes dans un echantillon |
FR3073047B1 (fr) | 2017-11-02 | 2021-01-29 | Commissariat Energie Atomique | Procede optique d'estimation d'un volume representatif de particules presentes dans un echantillon |
FR3076617B1 (fr) | 2018-01-08 | 2020-02-07 | Horiba Abx Sas | Systeme d'imagerie holographique et procede d'analyse par imagerie holographique avec detection de defauts dans la chambre d'observation |
FR3081997B1 (fr) | 2018-05-31 | 2021-04-16 | Commissariat Energie Atomique | Dispositif et procede d'observation de particules, en particulier des particules submicroniques |
FR3082944A1 (fr) | 2018-06-20 | 2019-12-27 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede d'observation d'un echantillon par imagerie sans lentille, avec prise en compte d'une dispersion spatiale dans l'echantillon |
FR3082943A1 (fr) | 2018-06-20 | 2019-12-27 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede de comptage de particules de petite taille dans un echantillon |
JP2020041928A (ja) * | 2018-09-11 | 2020-03-19 | 株式会社東芝 | セルフチェックシステム |
FR3086758B1 (fr) | 2018-09-28 | 2020-10-02 | Commissariat Energie Atomique | Procede et dispositif d'observation d'un echantillon sous lumiere ambiante |
FR3087009B1 (fr) | 2018-10-09 | 2020-10-09 | Commissariat Energie Atomique | Procede de determination de parametres d'une particule |
EP3839479B1 (en) * | 2019-12-20 | 2024-04-03 | IMEC vzw | A device for detecting particles in air |
FR3118169A1 (fr) | 2020-12-22 | 2022-06-24 | Commissariat à l'Energie Atomique et aux Energies Alternatives | Procédé de caractérisation de spermatozoïdes |
FR3138522A1 (fr) | 2022-07-29 | 2024-02-02 | Horiba Abx Sas | Dispositif de détection de particules en imagerie sans lentille |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090073521A1 (en) * | 2005-03-16 | 2009-03-19 | Consiglio Nazionale Delle Richerche | Holographic Method With Numerical Reconstruction for Obtaining an Image of a Three-Dimensional Object Which Even Points out of the Depth of Field Are in Focus, and Holographic Apparatus Using Such a Method |
US20130280752A1 (en) * | 2011-01-06 | 2013-10-24 | The Regents Of The University Of California | Lens-free tomographic imaging devices and methods |
US20140139625A1 (en) * | 2011-07-19 | 2014-05-22 | Ovizio Imaging Systems NV/SA | Method and system for detecting and/or classifying cancerous cells in a cell sample |
WO2015024020A1 (en) * | 2013-08-16 | 2015-02-19 | The General Hospital Corporation | Portable diffraction-based imaging and diagnostic systems and methods |
Family Cites Families (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1692296A (zh) * | 2002-09-30 | 2005-11-02 | 独立行政法人科学技术振兴机构 | 共焦点显微镜及采用它的荧光测量方法和偏振光测量方法 |
US20080208511A1 (en) * | 2005-03-07 | 2008-08-28 | Michael Trainer | Methods and apparatus for determining characteristics of particles |
SE530750C2 (sv) * | 2006-07-19 | 2008-09-02 | Hemocue Ab | En mätapparat, en metod och ett datorprogram |
GB0701201D0 (en) | 2007-01-22 | 2007-02-28 | Cancer Rec Tech Ltd | Cell mapping and tracking |
CN101842751B (zh) * | 2007-10-30 | 2014-10-22 | 纽约大学 | 用全息视频显微术来跟踪和表征颗粒 |
US8184298B2 (en) * | 2008-05-21 | 2012-05-22 | The Board Of Trustees Of The University Of Illinois | Spatial light interference microscopy and fourier transform light scattering for cell and tissue characterization |
KR20100098107A (ko) * | 2009-02-27 | 2010-09-06 | 휴먼전자 주식회사 | 휴대용 세포 분석장치 및 분석방법 |
EP3671176B1 (en) * | 2009-10-20 | 2022-04-13 | The Regents of the University of California | Incoherent lensfree cell holography and microscopy on a chip |
US9176152B2 (en) * | 2010-05-25 | 2015-11-03 | Arryx, Inc | Methods and apparatuses for detection of positional freedom of particles in biological and chemical analyses and applications in immunodiagnostics |
US9569664B2 (en) | 2010-10-26 | 2017-02-14 | California Institute Of Technology | Methods for rapid distinction between debris and growing cells |
BR112013010262B1 (pt) * | 2010-11-12 | 2022-01-25 | Universite Libre De Bruxelles | Processo para caracterização de objetos transparentes em um meio transparente |
WO2012082776A2 (en) | 2010-12-14 | 2012-06-21 | The Regents Of The University Of California | Method and device for holographic opto-fluidic microscopy |
WO2012112114A1 (en) * | 2011-02-16 | 2012-08-23 | Phase Holographic Imaging Phi Ab | Determination of physical parameters associated with erythrocytes |
BR112014013351A2 (pt) * | 2011-12-02 | 2017-06-13 | Csir | sistema, método e dispositivo de análise de material |
CN102660457A (zh) * | 2012-04-17 | 2012-09-12 | 南昌航空大学 | 无透镜全息衍射成像血细胞分析计数的装置及方法 |
EP2667178A3 (en) * | 2012-05-21 | 2018-01-03 | IMEC vzw | Holographic imaging for analyzing molecules |
US9588037B2 (en) | 2012-07-13 | 2017-03-07 | The Regents Of The University Of California | High throughput lens-free three-dimensional tracking of sperm |
FR2993372B1 (fr) * | 2012-07-13 | 2015-04-10 | Commissariat Energie Atomique | Procede et systeme de reconstruction de proprietes optiques d'objets diffractants baignant dans un milieu liquide |
FR3002634B1 (fr) * | 2013-02-28 | 2015-04-10 | Commissariat Energie Atomique | Procede d'observation d'au moins un objet, tel qu'une entite biologique, et systeme d'imagerie associe |
FR3020682B1 (fr) * | 2014-04-30 | 2016-05-27 | Commissariat Energie Atomique | Procede et systeme de detection d'au moins une particule dans un liquide corporel, et procede associe de diagnostic de la meningite |
FR3030749B1 (fr) * | 2014-12-19 | 2020-01-03 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Methode d'identification de particules biologiques par piles d'images holographiques defocalisees |
-
2015
- 2015-03-24 FR FR1552443A patent/FR3034196B1/fr not_active Expired - Fee Related
-
2016
- 2016-03-23 CN CN201680022598.1A patent/CN107532989B/zh active Active
- 2016-03-23 KR KR1020177030554A patent/KR102479862B1/ko active IP Right Grant
- 2016-03-23 WO PCT/FR2016/050643 patent/WO2016151248A1/fr active Application Filing
- 2016-03-23 US US15/560,763 patent/US10379027B2/en active Active
- 2016-03-23 EP EP16717422.6A patent/EP3274689B1/fr active Active
- 2016-03-23 JP JP2017550232A patent/JP6727229B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090073521A1 (en) * | 2005-03-16 | 2009-03-19 | Consiglio Nazionale Delle Richerche | Holographic Method With Numerical Reconstruction for Obtaining an Image of a Three-Dimensional Object Which Even Points out of the Depth of Field Are in Focus, and Holographic Apparatus Using Such a Method |
US20130280752A1 (en) * | 2011-01-06 | 2013-10-24 | The Regents Of The University Of California | Lens-free tomographic imaging devices and methods |
US20140139625A1 (en) * | 2011-07-19 | 2014-05-22 | Ovizio Imaging Systems NV/SA | Method and system for detecting and/or classifying cancerous cells in a cell sample |
WO2015024020A1 (en) * | 2013-08-16 | 2015-02-19 | The General Hospital Corporation | Portable diffraction-based imaging and diagnostic systems and methods |
Also Published As
Publication number | Publication date |
---|---|
FR3034196A1 (fr) | 2016-09-30 |
JP2018514759A (ja) | 2018-06-07 |
US10379027B2 (en) | 2019-08-13 |
WO2016151248A1 (fr) | 2016-09-29 |
KR102479862B1 (ko) | 2022-12-21 |
EP3274689B1 (fr) | 2022-02-16 |
US20180080760A1 (en) | 2018-03-22 |
JP6727229B2 (ja) | 2020-07-22 |
EP3274689A1 (fr) | 2018-01-31 |
FR3034196B1 (fr) | 2019-05-31 |
CN107532989A (zh) | 2018-01-02 |
KR20180011762A (ko) | 2018-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107532989B (zh) | 用于分析颗粒的方法 | |
US10481076B2 (en) | Method for determining the state of a cell | |
CN109074025B (zh) | 用于表征样本中的颗粒的全息方法 | |
US20180052425A1 (en) | Device and method for nanoparticle sizing based on time-resolved on-chip microscopy | |
JP2022008632A (ja) | 分析方法 | |
JP2019512697A (ja) | デジタルホログラフィ顕微鏡検査および無傷の(untouched)末梢血白血球を用いる高精度の5部鑑別(5−part Differential) | |
Walcutt et al. | Assessment of holographic microscopy for quantifying marine particle size and concentration | |
JP6986083B2 (ja) | レンズレス撮像によるサンプル内の粒子の計数方法 | |
US10908060B2 (en) | Method for determining parameters of a particle | |
US9581429B2 (en) | Method and system for reconstructing optical properties of diffracting objects immersed in a liquid medium | |
Kim et al. | Volumetric monitoring of airborne particulate matter concentration using smartphone-based digital holographic microscopy and deep learning | |
Go et al. | Three-dimensional volumetric monitoring of settling particulate matters on a leaf using digital in-line holographic microscopy | |
Pégard et al. | Three-dimensional deconvolution microfluidic microscopy using a tilted channel | |
EP3729053B1 (en) | Fast and robust fourier domain-based cell differentiation | |
van Grootheest et al. | Large volume holographic imaging for biological sample analysis | |
Santoro et al. | A Novel Method for Evaluation of Ore Minerals Based on Optical Microscopy and Image Analysis: Preliminary Results | |
US20240020949A1 (en) | Method for classifying an input image representing a particle in a sample | |
Cornelis et al. | Fast and robust Fourier domain-based classification for on-chip lens-free flow cytometry | |
US20230386176A1 (en) | Method for classifying an input image representing a particle in a sample | |
Ronen et al. | Plankton reconstruction through robust statistical optical tomography | |
CN116888644A (zh) | 用于对代表随时间推移的样本中粒子的输入图像序列进行分类的方法 | |
Dyomin et al. | Investigation of particles located in the water by digital holography | |
US20230386232A1 (en) | Method for classifying an input image containing a particle in a sample | |
Robison | Imaging White Blood Cells using a Snapshot Hyper-Spectral Imaging System | |
Röhrl et al. | Autoencoder Features for Differentiation of Leukocytes Based on Digital Holographic Microscopy (DHM) |
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 |