CN108496070B - 通过无透镜成像观察样本的方法 - Google Patents

通过无透镜成像观察样本的方法 Download PDF

Info

Publication number
CN108496070B
CN108496070B CN201680079063.8A CN201680079063A CN108496070B CN 108496070 B CN108496070 B CN 108496070B CN 201680079063 A CN201680079063 A CN 201680079063A CN 108496070 B CN108496070 B CN 108496070B
Authority
CN
China
Prior art keywords
intensity
image
sample
laser diode
image sensor
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
Application number
CN201680079063.8A
Other languages
English (en)
Other versions
CN108496070A (zh
Inventor
皮埃尔·布朗丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
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 Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Publication of CN108496070A publication Critical patent/CN108496070A/zh
Application granted granted Critical
Publication of CN108496070B publication Critical patent/CN108496070B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Electro-optical investigation, e.g. flow cytometers
    • G01N15/1434Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its optical arrangement
    • G01N15/1433
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/0443Digital holography, i.e. recording holograms with digital recording means
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/0465Particular recording light; Beam shape or geometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N2015/0007Investigating dispersion of gas
    • G01N2015/0011Investigating dispersion of gas in liquids, e.g. bubbles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N2015/0042Investigating dispersion of solids
    • G01N2015/0053Investigating dispersion of solids in liquids, e.g. trouble
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N2015/1006Investigating individual particles for cytology
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Electro-optical investigation, e.g. flow cytometers
    • G01N15/1434Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its optical arrangement
    • G01N2015/1454Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its optical arrangement using phase shift or interference, e.g. for improving contrast
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/0443Digital holography, i.e. recording holograms with digital recording means
    • G03H2001/0447In-line recording arrangement
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • G03H1/04Processes or apparatus for producing holograms
    • G03H1/0443Digital holography, i.e. recording holograms with digital recording means
    • G03H2001/0452Digital holography, i.e. recording holograms with digital recording means arranged to record an image of the object
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H2222/00Light sources or light beam properties
    • G03H2222/10Spectral composition
    • G03H2222/12Single or narrow bandwidth source, e.g. laser, light emitting diode [LED]
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H2240/00Hologram nature or properties
    • G03H2240/50Parameters or numerical values associated with holography, e.g. peel strength
    • G03H2240/51Intensity, power or luminance
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03HHOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H2240/00Hologram nature or properties
    • G03H2240/50Parameters or numerical values associated with holography, e.g. peel strength
    • G03H2240/52Exposure parameters, e.g. time, intensity

Abstract

一种用于通过无透镜成像观察样本(10)的方法,其中样本(10)设置在激光二极管(11)和图像传感器(20)之间,所述激光二极管(11)被供应有强度(i)小于或等于临界强度(ic)的电源电流。在预备步骤期间确定该临界强度(ic),在此期间强度(i)初始大于称为二极管(11)的激光阈值的阈值(is)。通过观察在图像传感器(20)处形成的图像(25),强度(i)降低直到在所形成的图像(25)上观察到干涉图像的衰减,对应于该衰减处的强度的临界强度(ic)是最佳的。

Description

通过无透镜成像观察样本的方法
技术领域
本发明的技术领域涉及通过无透镜成像观察样本,特别是生物样本。
背景技术
通过无透镜成像观察样本、特别是生物样本,在近十年中得到了显著的进步。该技术允许通过将样本设置在光源和图像传感器之间来观察样本,而无需在样本和图像传感器之间设置光学放大镜。图像传感器因此收集经样本透射的光波的图像。
该图像是由光波和衍射波之间的干涉图案形成的,其中光波是由光源发射并经样本透射的,衍射波是由光源发射的光波经样本的衍射而产生的。有时用术语“衍射图案”来指代这些干涉图案。在图像传感器上形成的图像可以通过数字传播算法来进行处理,以估计样本的光学特性。这种算法在全息重建领域是众所周知的。为此,在样本与图像传感器之间的距离已知的情况下,应用考虑了该距离的全息重建算法。
Garcia-Sucerquia J.出版的“数字同轴全息显微术”(Digital in-lineholographic microscopy),应用光学,卷45(5),2006年2月10日描述了使用激光束观察颗粒(例如生物颗粒),以及将重建算法应用于形成在CCD传感器上的图像。
文献WO 2008090330已经示出,通过用发光二极管代替激光光源,利用空间滤波,可以通过无透镜成像获得生物样本(在这种情况下是细胞)的可用图像。该文献中描述的装置使得可以将每个细胞与干涉图案相关联,其中干涉图案的形态使得可以识别细胞的类型。接下来其他出版物证实了这种技术的益处,例如申请US 2012/0218379。
目前,无透镜成像似乎是传统显微镜的简单且廉价的替代品。而且,其观察范围显著地大于显微镜的可能的观察范围。可以理解,与该技术相关的应用前景非常广阔。
发明人提出了一种与上述方法相比得到改进的用于观察样本的方法,使得可以获得与现有技术相比其信噪比得到改善的图像。
发明内容
本发明的第一主题是一种用于观察样本的方法,其包括以下步骤:
a)使用激光二极管照射样本,激光二极管能够发射入射光波,所述激光二极管被供应有电源电流,使得当所述电源电流超过称为阈值强度的强度时获得激光效应;
b)使用图像传感器获取样本的图像,样本位于激光二极管和图像传感器之间,每个图像表示在所述照射的作用下经样本透射的光波;
该方法的特征在于,激光二极管的电源电流的强度小于或等于所谓的临界强度,在步骤a)之前根据以下预备步骤来确定所述临界强度:
i)将所述激光二极管定位成使其发射的入射光波传播至图像传感器;
ii)向所述激光二极管供应所谓的初始电源电流,所述初始电源电流的强度大于或等于所述阈值强度,并且使用所述图像传感器获取称为初始图像的图像,所述初始图像包括延伸穿过图像的横向干涉条纹;
iii)改变所述电源电流的强度并且使用图像传感器获取图像,重复该步骤直到获取具有与所述初始图像相比衰减的横向干涉条纹的图像,电流的强度则对应于所述临界强度。
衰减的干涉条纹被理解为是指这些干涉条纹是不可见的,或者当它们可见时,与初始图像的干涉条纹相比是衰减的。在初始图像上的两个相邻的横向干涉条纹的特征可以由表示所述干涉条纹之间的对比度的对比度指标来表征。衰减对应于所述对比度指标降低至少25%,或甚至更多。
该方法可以包括以下特征中的任何一个,单独地或以技术上可行的组合:
-所述临界强度小于所述阈值强度。临界强度可以在阈值强度的0.5倍和阈值强度之间,并且优选地在阈值强度的0.7倍或0.8倍和阈值强度之间。
-在预备步骤i)至iii)中,表示所观察的样本的参考样本位于图像传感器和激光二极管之间。作为替代方案,在预备步骤i)至iii)中,没有样本位于图像传感器和激光二极管之间。
-样本包括衍射元素,衍射元素能够在所述光波照射的作用下形成衍射波,使得图像传感器暴露于称为曝光波的波,所述称为曝光波的波包括经样本透射的入射光波以及衍射波。样本可以特别地包括液体,例如其中悬浮有颗粒的生物液体,样本的厚度例如在50μm和2cm之间。样本也可以是厚度小于50μm的组织切片。
-在步骤b)中,电源电流的强度在阈值强度的0.5倍和所述阈值强度之间,或者甚至在阈值强度的0.5倍和阈值强度的0.9倍之间。优选地,该强度在阈值强度的0.7和0.9倍至阈值强度的0.95或1倍之间。
-激光二极管在可见光谱范围内发射。
本发明的另一主题是一种用于观察样本的方法,其包括以下步骤:
a)使用能够产生朝向样本传播的光波(称为入射光波)的激光二极管照射所述样本;
b)使用图像传感器获取样本在检测平面中形成的图像,样本位于激光二极管和图像传感器之间,每个图像表示在所述照射的作用下经样本透射的光波;
该方法的特征在于,所述激光二极管具有与其相关联的阈值强度,该阈值强度对应于二极管的电源电流的强度,超过该强度则获得激光效应,并且当样本被照射时,向激光二极管供应强度严格小于阈值强度(例如在阈值强度的0.5倍和所述阈值强度之间)的电源电流。
优选地,该强度在阈值强度的0.7倍和阈值强度的0.9到0.95倍之间。
从以下通过非限制性示例给出且在附图中示出的本发明的特定实施例的描述中,其他优点和特征将变得更加显而易见,其中:
附图说明
图1示出了用于观察样本的示例性装置。
图2A示出了由激光二极管发射的光功率随所述激光二极管的电源电流的强度的变化。图2B示出了激光二极管的光谱随激光二极管的电源电流的强度的变化。图2C示出了由所述激光二极管发射的光波的相干长度随激光二极管的电源电流的强度的变化。
图3A、3B、3C、3D、3E和3F示出了使用参照图1所描述的装置获得的样本的图像,考虑了用于照射样本的激光二极管的电源电流的各种强度。
图4A示出了图3A、3B、3C、3D、3E和3F中存在的衍射图案的径向轮廓随激光二极管的电源电流的强度的变化。
图4B示出了图4A所示的轮廓的标准偏差随激光二极管的电源电流的强度的变化。
图5A和5B示出了使用图1所示的装置获得的图像,其中没有在激光二极管和样本之间插入样本,激光二极管的电源电流的强度分别为两个不同的值,分别低于和高于所谓的阈值强度。
图5C示出了表示在填充有液体的流体室插入在激光二极管和样本之间的情况下,使用图1所示的装置获得的图像的背景的变异系数的变化,该变异系数随激光二极管的电源电流的强度变化。
图6A、6B、6C、6D、6E和6F示出了随重建距离变化的、沿着重建图像中的线延伸的像素的强度,所述图像是通过分别对由图像传感器获取的两个珠的图像应用全息重建算法来获得的,其中所述两个珠的图像是利用激光二极管的电源电流的各种强度来获得的。
图7A示出了使用现有技术装置获得的样本图像,其中光源是发光二极管。图7B示出了使用根据本发明的装置产生的样本图像,其中适当地调整激光二极管的强度。
具体实施方式
图1示出了本发明的一个主题的示例性装置。激光二极管11能够沿着传播轴Z在样本10的方向上产生称为入射光波的光波12。在该示例中,激光二极管是以Civillaser销售的二极管,其以接近405nm的波长发射。激光二极管因此在可见光谱范围内发射。
样本10可以是期望表征的生物样本。样本10还可以包括含有待表征的颗粒13的固体或液体介质14。这可能涉及例如培养基或体液中的生物颗粒。生物颗粒被理解为是指细胞、细菌或其他微生物,真菌、孢子等。术语颗粒还可以指生物学方案中常用的微珠,例如金属微珠、玻璃微珠或有机微珠。这也可能涉及悬浮在液体介质中的不溶性液滴,例如水包油乳液中的脂滴。因此,术语颗粒既指最初存在于被检样本中的内源性颗粒,也指分析前添加到该样本中的外源性颗粒。通常,颗粒的尺寸有利地小于1mm,或者甚至小于500μm,并且优选地在0.5μm和500μm之间。样本10的厚度ε允许入射光12中的一部分透射到图像传感器20,形成所谓的透射光波22。当样本是含有颗粒的液体时,取决于样本对入射光的衰减,厚度ε可以在几十微米到几毫米或甚至几厘米之间。
根据另一示例,样本可以是薄组织切片或解剖病理切片,包括沉积在透明载玻片上的厚度薄的组织。薄的厚度理解为指优选小于100μm,并且优选小于10μm,通常为几微米的厚度。
激光二极管与样本之间的距离Δ优选大于1cm。所述距离优选在2cm到30cm之间。优选地,从样本处看,激光二极管被认为是点状的。这意味着激光二极管的直径(或其对角线)优选地小于样本和激光二极管之间的距离的十分之一,并且最好为百分之一。因此,光优选地以平面波的形式或者以可以被认为是平面波的波的形式到达样本。
装置包括能够在检测平面P20中形成图像的图像传感器20。在所示的示例中,这是一个包括CCD型或CMOS的像素矩阵阵列的矩阵阵列光电检测器。CMOS是优选的图像传感器,因为像素的尺寸较小,从而允许获得空间分辨率更好的图像。图像传感器包括像素矩阵阵列17,在该矩阵阵列上方设置有透明保护窗18。在该示例中,透明保护窗18设置在距离像素矩阵阵列17 125μm处,并且其厚度达到400μm。像素矩阵阵列17与保护窗18之间的距离通常在几十μm到150或200μm之间。在该示例中,检测器是一个8位CMOS传感器,包括3884x 2764像素,像素间距为1.67μm。检测平面P20优选垂直于入射光波12的传播轴Z延伸。像素间距小于3μm的图像传感器是优选的,以改善图像的空间分辨率。
在该示例中,样本10和图像传感器20的像素矩阵阵列17之间的距离d小于2cm,并且优选小于1cm。例如大约为5mm。样本10和图像传感器20的接近使得可以获得宽的观察范围。可替代地,可以移动样本使其更靠近激光二极管,在这种情况下观察范围会减小。然而,这样的配置使得可以获得关于样本的更详细的信息。
要注意的是,没有成像光学器件,例如图像传感器20和样本10之间的放大光学器件。这并不排除在图像传感器20的每个像素处存在聚焦微透镜的可能性,这些微透镜不具有图像放大功能。
在入射光波12的作用下,样本10可以产生在检测平面P20中易于产生干涉(特别是与经样本透射的光波22的干涉)的衍射波23。该干涉在由图像传感器20获取的图像25上产生多个基本衍射图案26,每个基本衍射图案包括中心区域和多个同心衍射环,每个环构成圆形干涉条纹。每个基本衍射图案都是由样本中的衍射对象13引起的,在这种情况下是颗粒。图像传感器20所暴露于的光波24包括:
-由入射光波12经样本的衍射产生的分量23;
-由入射光波12经样本的透射产生的分量22。
这两个分量在由图像传感器获取的图像25上产生干涉,该干涉为基本衍射图案26的形式。稍后将参照图3A至3F或7A和7B来示出基本衍射图案的示例。通常,基本衍射图案26可以与样本中的每个衍射元素相关联。样本10越靠近激光二极管11,基本衍射图案26就越详细,图像传感器上可识别的环的数量就越多。样本10越靠近图像传感器20,观察范围就越大,导致在同一个图像25上观察到更多数量的基本衍射图案26。
处理器30,例如微处理器,能够处理由图像传感器20获取的每个图像25。特别地,处理器是连接到可编程存储器32的微处理器,在可编程存储器32中存储有用于执行本说明书中描述的图像处理和计算操作的一系列指令。处理器30也可以连接到屏幕34。
图像处理的一个示例是全息重建算法的应用,其使得可以根据由图像传感器获取的图像25来重建称为重建图像的图像28。这种类型的算法是常规的,并且稍后将在说明书中参照图6A至6F来进行讨论。
激光二极管11连接到电路11c,从而允许向激光二极管供应电流。激光二极管11包括由半导体材料制成的有源区。超过一定的强度值(称为激光阈值强度is),可以观察到激光发射。该激光阈值强度is对应于足以在激光二极管的有源区中获得激光效应的电子注入。
旨在向所使用的激光二极管供应额定强度in等于150mA的电源电流。发明人使用光功率计测量随电源电流的强度i变化的由激光二极管传输的光功率。通过在额定强度两侧在15mA和170mA之间改变电源电流的强度i,在三个不同的时刻1、2和3进行三次测量。图2A中示出了所获得的曲线。观察到的是激光二极管的光功率特征的变化。每条曲线具有:
-第一部分,其中电荷到有源区的注入是不充足的。在第一操作状态I中,激光二极管的光发射被认为是自发的,并且发射的光功率与注入到有源区的电荷载流子的密度成比例,也就是说与电源电流的强度i成比例;
-第二部分,其中电荷的注入足以显著地增加激光二极管的输出,反映在所谓的受激光发射的现象中。在第二操作状态II,或激光操作状态中,获得了激光效应,并且光功率随着强度i,随着反映在曲线梯度的增加中的放大效应而线性增加。激光二极管的额定强度in等于150mA,对应于该操作状态。
第一和第二操作状态之间的转换发生在125mA和135mA之间,但是不能够建立精确的值。具体地,激光二极管的性能发生波动,特别是随着温度波动。此外,在上面概述的两种操作状态之间存在所谓的瞬态操作状态,其中二极管的发射是受激的从而不会补偿二极管的某些损耗。由于所述损耗,光功率随着强度的变化不再是线性的。在所示的曲线中,该瞬态状态在120和135mA之间的电源电流强度i之间延伸。然后可以认为对应于二极管的电源电流的、被称为激光阈值电流的阈值强度值is为140mA。
在本说明书中,阈值强度is表示这样的强度,在超过该强度时,会发生受激光的放大,这对应于上述激光操作状态(状态Ⅱ)的最小强度。如在图2A中可以看出的,该阈值强度is是这样的强度,当超过该强度时,由激光效应引起的光功率随着激光二极管的电源电流的强度线性增加。
在强度i的每个值处,还使用光谱仪确定由激光二极管11产生的光波12的发射光谱。在图2B中描绘了各种光谱。每个光谱具有以401nm的中心波长λc为中心的发射峰。如所预期的那样,只要电源电流的强度i低于阈值is,就会观察到发射峰的放大。
可以使用以下表达式来确定由激光二极管11产生的光波12的相干长度Lc
Figure BDA0001731529190000081
其中,λ表示波长,并且Δλ表示发射峰的宽度。
该表达式适用于发射光谱的峰值类似于高斯分布的条件,这是本示例中所采用的假设。
相干长度Lc反映了这样的事实:如果作为辅助光源的两个对象(例如反射或衍射对象)间隔开比相干长度短的距离,则由这些对象发射的光波会形成干涉。
图2C示出了由激光二极管发射的光波12的相干长度Lc随电源电流的强度i的变化。当电源电流的强度i小于阈值强度is时,该相干长度很短,大约几十微米,并且小于50微米或甚至20微米。超过强度阈值is时,相干长度达到或甚至超过100μm,并且由于所使用的光谱仪的光谱分辨率而无法被高精度地确定。
相干长度Lc是确定观察样本形成能够被检测到的干涉的能力的参数。在图1所示的示例中,样本包括透明液体14或能够被认为是其中悬浮有颗粒13的液体。这些颗粒是衍射的,并且引起在由图像传感器20形成的图像25上的衍射图案26的出现。
发明人研究了改变激光二极管11的电源电流的强度i对由图像传感器20获取的图像的影响。图3A、3B、3C、3D、3E和3F分别表示样本的图像,其中激光二极管的电源电流的强度i分别等于41mA、67mA、92mA、110mA、126mA和135mA。在这些图中出现基本衍射图案26,这些图案的大部分能够与存在于样本10中的颗粒13相关联。在该示例中,样本由直径为3μm和6μm的乳胶珠混合物组成,所述乳胶珠位于PBS(phosphate buffered saline(磷酸盐缓冲盐水)的首字母缩写)类型的盐水缓冲液中。在这些图像中的每一个上,调整曝光时间使得每个图像的平均灰度是相似的,使得不管照射功率如何变化,这些图像是相当的。可以看出,当电源电流的强度增加时,形成每个衍射图案的环的数量趋于增加,这是所预期的,因为这些环对应于圆形干涉条纹。因此,图像中存在的信息更加丰富。此外,如参考现有技术所表明的,这样的衍射图案可以经历全息重建,并且公认的是,当衍射环的数量和对比度增加时,这种重建的结果会得到改善。
另一个观察结果是横向干涉条纹29的出现,当电源电流的强度大于或接近阈值强度is时,横向干涉条纹29基本上是直线的并且延伸穿过图像。这些横向干涉条纹的出现显著地损害了图像的信噪比。与位于由图像传感器20形成的图像25的一些部分中的形成基本衍射图案26的同心干涉条纹相比,横向干涉条纹29延伸穿过整个图像。它们构成了既影响具有基本衍射图案26的图像区域,又影响形成图像背景的、不具有基本衍射图案的图像区域的噪声。这些横向条纹可以是直线的、椭圆的或圆形的。术语“横向”表示它们影响整个图像25,从该图像的一个边缘延伸到另一个边缘。
在每个图像25上已经选择了一个基本衍射图案26,对于该图案已经建立了平均径向轮廓。平均径向轮廓被理解为是指通过确定各种基本轮廓的平均值而获得的幅度轮廓,每个基本轮廓沿着在所选的衍射图案的中心和所述图的端部之间延伸的轴线27来截取。对于同一个衍射图案26,如图3A所示可以考虑几个轴线27,其中轴线27以白点的形式示出。每个图像上所选的基本图案26都由一个白点框构成。图4A示出了在图像3A至3F中的每一个上获得的平均径向轮廓。可以看出,i强度增加地越多,环的数量就越多。
已经确定了图4A中所示的每个径向轮廓的标准偏差。图4B示出了该标准偏差随激光二极管11的电源电流的强度i的变化。该指标表示所选的衍射图案的对比度,考虑了所有可识别的衍射环。其在i=121mA时达到最大值,然后在强度超过该值时减小。衍射图案的对比度越高,该图案越能够被识别,并且基于该衍射图案重建的图像的质量将越好。图3A至3F以及图4A和4B示出了二极管的电源电流的最佳强度范围为60mA至125mA之间,在该范围中衍射图案26的对比度较高。它们还示出了在超过130mA时,出现了干涉条纹,影响了整个图像25,大大增加了图像的噪声水平。而且,这种噪声不是固定的,使得难以消除。因此,重要的是通过调节其电源电流来使用激光二极管,使得其强度i足够高以使衍射图案具有一定对比度并且具有足够数量的衍射环,同时其强度低于称为临界强度ic的强度,超过该临界强度会出现延伸穿过由图像传感器20获取的图像25的干涉条纹。可以看出,在这个示例中,电源电流的强度的最佳范围在60mA和135mA之间,即在阈值强度is的大约0.5倍和所述阈值强度is之间,并且似乎特别优选在80mA和120mA之间,即在所述阈值强度is的0.6倍和0.9倍之间。
图5A和5B示出了在图像传感器20和激光二极管11之间没有插入样本10的情况下由图像传感器20获得的图像,电源电流的强度i分别等于120mA(即低于阈值强度is)和140mA(即大于或等于阈值强度)。当强度超过临界值ic时可以观察到干涉条纹的出现,临界值小于或等于阈值强度is。看起来过长的相干长度对于传感器获取的图像的质量是有害的。因此,优选调整激光二极管的电源电流的强度i,使得由激光二极管发射的光波12的相干长度小于100μm,或甚至小于80μm或50μm。发明人将这种干涉条纹的出现降低至在像素矩阵阵列17和检测器的距离所述像素矩阵阵列约100μm的保护窗口18之间发生的杂散反射。当样本具有相隔这个量级的两个平行界面时,这些干涉条纹也可能出现。通常,所观察的样本的厚度ε在20μm和200μm之间,或者甚至500μm。
发明人通过在70mA和150mA之间改变激光二极管11的电源电流的强度i来评估由图像传感器25获取的图像的变异系数CV的变化,样本10包括厚度ε=100μm,填充有其中没有衍射颗粒的盐水缓冲液的流体室。图像的变异系数CV被理解为是指图像的标准偏差除以所述图像的平均值。调整每个图像的曝光时间,使得每个图像的平均强度是恒定的,以补偿二极管的照射功率的变化。通过排除具有杂散衍射图案的某些图像部分来计算该变异系数。因此,其表示了每个图像的“背景”,即表示了不具有可能的衍射图案的区域。图5C示出了该变异系数CV随激光二极管的电源电流的强度i的变化。观察到在70mA和120mA之间是稳定的,随后显著增加,该显著增加反映了在图像中出现图案,即以上参照图3F和图5B概述的干涉条纹。
发明人观察到,通过将样本设置在图像传感器和激光二极管之间,临界强度变化了。因此看起来干涉条纹的出现还取决于位于激光二极管11和图像传感器20之间的样本10。
由此,激光二极管11的电源电流的强度i优选地通过执行以下步骤(称为预备步骤)来根据具体情况进行调整:
-在激光二极管的额定使用工况下获取初始图像,即通过采用大于阈值强度is的电源电流的强度i,使得通过激光效应获得由激光二极管发射的光波12。该图像具有横向干涉条纹,使得其不适于全息重建;
-改变激光二极管的电源电流的强度i,特别是降低该强度,并使用图像传感器获取图像,重复该步骤,直到观察到横向干涉条纹的显著衰减。当观察到这种衰减时,则认为激光二极管的电源电流的强度对应于所述临界强度ic
然后,根据该临界强度ic确定激光二极管11的电源电流的强度i,该强度i有助于通过无透镜成像观察样本。它可以等于该临界强度,或者稍微低于该临界强度,例如在0.8ic和ic之间。具体地,认为该强度i应该优选地接近临界强度。在临界强度以上,出现干涉条纹。在临界强度以下,光功率降低。通过低于临界强度,例如低于0.9ic或0.95ic,克服了临界强度的可能的波动,例如在激光二极管的漂移(例如热漂移)的影响下。这样的“安全边际”可能被证明是有用的。
上述步骤被认为是预备步骤,因为它们是在观察样本之前实施的。它们优选通过定位所谓的参考样本10ref来执行,参考样本10ref表示下文中要观察的样本。它们也可以在没有样本的情况下实施,图像传感器20在没有插入样本的情况下直接面向激光二极管11定位。
两个连续的干涉条纹可以由对比度指标ind来表征,该对比度指标ind代表表示亮条纹的最大强度水平Smax和表示邻近所述亮条纹的暗条纹的最小强度水平Smin之间的差。对比度指标可以通过简单的减法Smax-Smin来获得,或者通过加法Smax+Smin标准化该差值来获得。因此,对比度指标可以如下表示:
ind=Smax-Smin(2)或
Figure BDA0001731529190000111
表达式(3)对应于对比度的常规定义。
当干涉条纹不再显著地出现在由图像传感器获取的图像上,或者当它们充分显著地衰减时,则达到临界强度ic。因此,所述显著的衰减对应于初始图像和在临界强度下获得的图像之间的对比度指标ind相对于横向干涉条纹降低25%,或者甚至50%,或者甚至80%或更多。
如上所述,由图像传感器获取的每个图像25可以经历全息重建,旨在使用传播算子h(x,y,z)来对所述图像进行卷积运算。传播算子h(x,y,z)的作用是描述图像传感器20和坐标点(x,y,z)之间的光传播。例如在Ryle等人的出版物“生物标本的数字同轴全息(Digital in-line holography of biological specimens)”,国际光学工程学会会议(SPIE),第6311卷(2006)中描述了全息重建的原理。全息重建使得可以确定光波24的复数表达式U(x,y,z),其中检测器在空间中(特别是在距离图像传感器的距离|z|处的平面中)的任何坐标点(x,y,z)处暴露于光波24。然后可以使用以下表达式来确定该光波的幅度u(x,y,z)和相位
Figure BDA0001731529190000121
u(x,y,z)=abs[U(x,y,z)](4),或者
Figure BDA0001731529190000122
Figure BDA0001731529190000123
运算abs和arg分别返回模和幅角。
数字重建特别可以基于菲涅耳衍射模型。在该示例中,传播算子是菲涅耳-亥姆霍兹方程,如下所示:
Figure BDA0001731529190000124
其中λ表示波长。
由此,
Figure BDA0001731529190000125
其中
-x’和y’表示图像传感器的平面中的坐标;
-S(x’,y’)表示图像传感器在坐标x’和y’处获取的图像的强度;
-x和y表示重建平面中的坐标,该重建平面位于距图像传感器的距离|z|处;
-z表示重建图像沿着入射光波的传播轴Z的坐标。
发明人使用图1所示的装置,通过将两个边到边直径为6μm的乳胶珠定位在厚度为100μm的、填充有上述定义的PBS盐水缓冲液的流体室的底部来获得图像25。使用图像传感器20通过在41mA和135mA之间调制激光二极管11的电源电流的强度i来获取图像25。
根据上述原理,考虑等于图像传感器和所述珠之间的距离的z值,该图像经历全息重建。图像6A、6B、6C、6D、6E和6F示出了随重建距离|z|变化的重建图像28的像素的强度。更准确地说,每个图像表示随着重建距离|z|变化的位于穿过两个珠的线上的像素的强度分布,激光二极管11的电源电流的强度i分别达到41mA、67mA、92mA、110mA、126mA和135mA。虚线对应于等于检测器和珠之间的距离的重建距离。可以看到对应于i=126mA的图6E是最佳的,并且看到图像6D(i=110mA)或6C(i=92mA)也能够清楚地呈现出重建图像28上的珠的分离。图6F(i=135mA)实现了满意的分离,但是由于存在上述横向干涉条纹,在图像的其余部分中受到高背景噪声的影响。
图7A和图7B示出了当光源分别是蓝色发光二极管和激光二极管时使用图像传感器获取的图像,其中光源的强度等于临界强度ic,在这种情况下为125mA。图7B中形成的衍射图案包括比图7A中形成的更多的可识别的衍射环。
在调节其电源电流的强度的条件下使用激光二极管使得可以获得样本的更精确的观察。
可以在体液分析中实施本发明,以实现检测和/或识别形成所述体液的颗粒的目的。本发明也可以允许表征生物组织薄片。除生物学以外,本发明可以应用于检测和/或识别工业流体中、农产品中或固体基质中的颗粒。

Claims (10)

1.一种观察样本(10)的方法,包括以下步骤:
a)使用激光二极管(11)照射所述样本,所述激光二极管被配置为发射朝向所述样本传播的入射光波(12),所述激光二极管被供应有电源电流,使得当所述电源电流的强度(i)超过阈值强度(is)时获得激光效应;
b)使用图像传感器(20)获取所述样本的图像(25),所述样本位于所述激光二极管(11)和所述图像传感器(20)之间,每个图像表示由所述入射光波(12)经所述样本的透射产生的透射光波(22);
其中,所述激光二极管的电源电流的强度(i)小于或等于临界强度(ic),在步骤a)之前根据以下预备步骤来确定所述临界强度:
i)将所述激光二极管(11)定位成使其发射的所述入射光波(12)向所述图像传感器(20)传播;
ii)向所述激光二极管(11)供应初始电源电流,所述初始电源电流的强度大于或等于所述阈值强度(is),并且使用所述图像传感器获取初始图像,所述初始图像具有横向干涉条纹;
iii)改变所述电源电流的强度并且使用所述图像传感器(20)获取图像,重复该步骤直到获取具有与所述初始图像相比衰减的横向干涉条纹的图像,所述电源电流的强度则对应于所述临界强度(ic)。
2.根据权利要求1所述的观察样本(10)的方法,其中,所述临界强度(ic)小于所述阈值强度(is)。
3.根据权利要求1所述的观察样本(10)的方法,其中,所述临界强度在所述阈值强度(is)的0.7倍和所述阈值强度(is)之间,或者在所述阈值强度(is)的0.8倍和所述阈值强度(is)之间。
4.根据权利要求1所述的观察样本的方法,其中,在所述初始图像上,两个相邻的横向干涉条纹由对比度指标(ind)来表征,所述对比度指标表示所述干涉条纹之间的对比度,所述临界强度对应于所述对比度指标降低超过25%。
5.根据权利要求1所述的观察样本(10)的方法,其中,在预备步骤i)至iii)中,表示所观察的样本的参考样本位于所述图像传感器(20)和所述激光二极管(11)之间。
6.根据权利要求1所述的观察样本(10)的方法,其中,在预备步骤i)至iii)中,没有样本位于所述图像传感器(20)和所述激光二极管(11)之间。
7.根据权利要求1所述的观察样本(10)的方法,其中,所述样本(10)包括衍射元素(13),所述衍射元素在所述入射光波(12)照射的作用下形成衍射波(23),使得所述图像传感器暴露于曝光波(24),所述曝光波(24)包括所述透射光波(22)和所述衍射波(23)。
8.根据权利要求1所述的观察样本(10)的方法,其中,在步骤b)中,所述电源电流的强度(i)在所述阈值强度的0.5倍和所述阈值强度(is)之间,或者在所述阈值强度的0.7倍和所述阈值强度之间。
9.根据权利要求1所述的观察样本(10)的方法,其中,在步骤b)中,所述电源电流的强度(i)在所述阈值强度(is)的0.5倍和所述阈值强度的0.9倍之间,或者在所述阈值强度的0.7倍和所述阈值强度的0.9倍之间。
10.根据权利要求1所述的观察样本(10)的方法,其中,没有放大光学器件在所述样本和所述图像传感器(20)之间延伸。
CN201680079063.8A 2015-12-24 2016-12-22 通过无透镜成像观察样本的方法 Active CN108496070B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1563320A FR3046238B1 (fr) 2015-12-24 2015-12-24 Procede d’observation d’un echantillon par imagerie sans lentille
FR1563320 2015-12-24
PCT/FR2016/053648 WO2017109428A1 (fr) 2015-12-24 2016-12-22 Procédé d'observation d'un échantillon par imagerie sans lentille

Publications (2)

Publication Number Publication Date
CN108496070A CN108496070A (zh) 2018-09-04
CN108496070B true CN108496070B (zh) 2021-06-22

Family

ID=56084122

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680079063.8A Active CN108496070B (zh) 2015-12-24 2016-12-22 通过无透镜成像观察样本的方法

Country Status (7)

Country Link
US (1) US11221290B2 (zh)
EP (1) EP3394591B1 (zh)
JP (1) JP6832355B2 (zh)
KR (1) KR102562719B1 (zh)
CN (1) CN108496070B (zh)
FR (1) FR3046238B1 (zh)
WO (1) WO2017109428A1 (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN109001900A (zh) * 2018-09-05 2018-12-14 南京大学 一种明场和荧光双模态的显微成像系统及方法
FR3090107B1 (fr) 2018-12-18 2020-12-25 Commissariat Energie Atomique Procédé de caractérisation d'une particule à partir d'un hologramme.
KR102111124B1 (ko) * 2019-05-28 2020-05-15 (주)위아프렌즈 무렌즈 단층 촬영 회절영상의 재구성방법 및 이 방법이 사용되는 휴대용 수질입자 분석장치
CN110836867A (zh) * 2019-10-18 2020-02-25 南京大学 一种基于卷积神经网络的无透镜全息显微微粒表征方法
KR102497160B1 (ko) 2020-11-08 2023-02-06 단국대학교 산학협력단 무렌즈 3차원 측정장치 및 이의 운용방법

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101471332A (zh) * 2007-12-28 2009-07-01 友嘉科技股份有限公司 含有激光二极管与发光二极管的光学装置
CN102077077A (zh) * 2008-06-25 2011-05-25 赫拉巴Abx公司 分类和计数微观元素的电光测量装置和电光测量方法
US8293524B2 (en) * 2006-03-31 2012-10-23 Fluxion Biosciences Inc. Methods and apparatus for the manipulation of particle suspensions and testing thereof
CN104204767A (zh) * 2012-02-13 2014-12-10 奥维奇奥成像系统股份有限公司 有数字全息显微镜的流式细胞仪

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3213333B2 (ja) * 1991-05-14 2001-10-02 シスメックス株式会社 尿中の細胞分析装置および方法。
JP3267685B2 (ja) * 1992-07-28 2002-03-18 セイコーインスツルメンツ株式会社 微細表面観察装置
US6784981B1 (en) * 2000-06-02 2004-08-31 Idexx Laboratories, Inc. Flow cytometry-based hematology system
US7515620B2 (en) * 2002-10-29 2009-04-07 Finisar Corporation Method for optimizing laser diode operating current
DE102005009317B4 (de) * 2005-02-18 2007-05-31 Atmel Germany Gmbh Verfahren und Vorrichtung zur Bestimmung einer Laserschwelle einer Laserdiode
US20060279819A1 (en) * 2005-05-26 2006-12-14 Inphase Technologies, Inc. Laser mode stabilization using an etalon
US8368897B2 (en) * 2008-08-22 2013-02-05 Ciencia, Inc. Versatile surface plasmon resonance analyzer with an integral surface plasmon resonance enhanced fluorescence mode
KR101030303B1 (ko) * 2009-03-17 2011-04-19 연세대학교 산학협력단 이미지 센서를 이용한 배열 광원의 특성 평가 장치 및 방법
JP2011038829A (ja) * 2009-08-07 2011-02-24 Topcon Corp 干渉顕微鏡及び測定装置
EP2596320A4 (en) * 2010-07-24 2015-08-26 Focused Innovation Inc METHOD AND APPARATUS FOR IMAGING
JP2014507645A (ja) * 2011-01-06 2014-03-27 ザ・リージェンツ・オブ・ザ・ユニバーシティー・オブ・カリフォルニア 無レンズ断層撮影装置及び方法
US8866063B2 (en) * 2011-03-31 2014-10-21 The Regents Of The University Of California Lens-free wide-field super-resolution imaging device
FR2991457B1 (fr) * 2012-06-01 2014-07-18 Commissariat Energie Atomique Procede et systeme de caracterisation de la vitesse de deplacement de particules contenues dans un liquide, telles que des particules sanguines
CN104755908B (zh) * 2012-07-27 2017-12-12 统雷有限公司 敏捷成像系统
FR3009084B1 (fr) * 2013-07-23 2015-08-07 Commissariat Energie Atomique Procede pour trier des cellules et dispositif associe.
EP2887008A1 (en) * 2013-12-17 2015-06-24 Universitat Politècnica de Catalunya A method of measuring a change in an optical path using differential laser self-mixing interferometry and a differential laser self-mixing interferometry measuring system
GB2544507B (en) * 2015-11-18 2017-11-22 Servomex Group Ltd Method and system for reduction of influence of baseline distortion in absorption spectroscopy measurements

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8293524B2 (en) * 2006-03-31 2012-10-23 Fluxion Biosciences Inc. Methods and apparatus for the manipulation of particle suspensions and testing thereof
CN101471332A (zh) * 2007-12-28 2009-07-01 友嘉科技股份有限公司 含有激光二极管与发光二极管的光学装置
CN102077077A (zh) * 2008-06-25 2011-05-25 赫拉巴Abx公司 分类和计数微观元素的电光测量装置和电光测量方法
CN104204767A (zh) * 2012-02-13 2014-12-10 奥维奇奥成像系统股份有限公司 有数字全息显微镜的流式细胞仪

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
determination of the size of irregular particles using interferometric out-of-focus imaging;Marc Brunel;《international journal of optics》;20140109;第1-8页 *
透射电镜多图像拼接法测量纳米氧化镍的颗粒粒径;韩冬;《电子显微学报》;20090630;第200-203页 *

Also Published As

Publication number Publication date
FR3046238B1 (fr) 2018-01-26
JP6832355B2 (ja) 2021-02-24
KR20180097679A (ko) 2018-08-31
JP2019509466A (ja) 2019-04-04
EP3394591A1 (fr) 2018-10-31
US11221290B2 (en) 2022-01-11
FR3046238A1 (fr) 2017-06-30
WO2017109428A1 (fr) 2017-06-29
EP3394591B1 (fr) 2019-11-06
CN108496070A (zh) 2018-09-04
KR102562719B1 (ko) 2023-08-03
US20210208053A1 (en) 2021-07-08

Similar Documents

Publication Publication Date Title
CN108496070B (zh) 通过无透镜成像观察样本的方法
US10816454B2 (en) Method for observing a sample, by calculation of a complex image
Zeng et al. 3D resolved mapping of optical aberrations in thick tissues
CN113960908B (zh) 用于表征样本中的颗粒的全息方法
EP2144034A1 (en) Multilayer structure measuring method and multilayer structure measuring apparatus
KR20180011762A (ko) 입자 분석 방법
US20180113064A1 (en) Method for determining the state of a cell
EP3467595A1 (en) Holographic observation method and device
US9404857B2 (en) White light diffraction tomography of unlabeled live cells
JP6667539B2 (ja) 光遠隔測定装置
Weigel et al. Resolution in the ApoTome and the confocal laser scanning microscope: comparison
JP2020514704A (ja) レンズレス撮像によるサンプル内の粒子の計数方法
Liu et al. Visualizing the ultra-structure of microorganisms using table-top extreme ultraviolet imaging
US9581429B2 (en) Method and system for reconstructing optical properties of diffracting objects immersed in a liquid medium
León-Rodríguez et al. High topographical accuracy by optical shot noise reduction in digital holographic microscopy
US10724937B2 (en) Device and method for bimodal observation of an object
Owen et al. Super-resolution imaging by localization microscopy
US10852696B2 (en) Device and method allowing observation of an object with a large field of observation without use of magnifying optics between a light source and the object
EP3230713B1 (fr) Procédé d'obtention d'une image d'un échantillon, et système associé d'imagerie sans lentille
JP2020190459A5 (zh)
EP3853651B1 (en) Confocal laser scanning microscope configured for generating line foci
FR3082943A1 (fr) Procede de comptage de particules de petite taille dans un echantillon
Patel et al. Identification of malaria infected red blood samples by digital holographic quantitative phase microscope
McFarlane et al. Characterisation of a deep-ultraviolet light-emitting diode emission pattern via fluorescence
CN114624193A (zh) 一种基于宽光谱偏焦相位恢复的生物成像装置及方法

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