CN103620509A - 测定与红细胞相关的物理参数 - Google Patents

测定与红细胞相关的物理参数 Download PDF

Info

Publication number
CN103620509A
CN103620509A CN201280007685.1A CN201280007685A CN103620509A CN 103620509 A CN103620509 A CN 103620509A CN 201280007685 A CN201280007685 A CN 201280007685A CN 103620509 A CN103620509 A CN 103620509A
Authority
CN
China
Prior art keywords
lambda
erythrocytic
mean corpuscular
epsiv
interference pattern
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
CN201280007685.1A
Other languages
English (en)
Other versions
CN103620509B (zh
Inventor
彼得·埃格尔伯格
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.)
Phase Holographic Imaging PHI AB
Original Assignee
Phase Holographic Imaging PHI AB
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 Phase Holographic Imaging PHI AB filed Critical Phase Holographic Imaging PHI AB
Publication of CN103620509A publication Critical patent/CN103620509A/zh
Application granted granted Critical
Publication of CN103620509B publication Critical patent/CN103620509B/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
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1429Signal processing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1468Optical investigation techniques, e.g. flow cytometry with spatial resolution of the texture or inner structure of the particle
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
    • G01N21/453Holographic interferometry
    • 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/08Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • 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/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434Optical arrangements
    • G01N2015/1454Optical arrangements 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

Landscapes

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

Abstract

提供了一种基于数字全息术来测定至少一个红细胞的参数的方法,所述方法包括以下步骤:构建包括表示所述至少一个红细胞的波前的干涉图样,所述干涉图样由所述至少一个红细胞和电磁辐射的相互作用产生,根据所述干涉图样对表示所述至少一个红细胞波前的振幅信息和相位信息进行重构,以及通过对表示所述至少一个红细胞的相位信息和振幅信息进行组合,来测定所述至少一个红细胞的平均红细胞体积和/或平均红细胞血红蛋白浓度和/或氧饱和度和/或平均红细胞血红蛋白。

Description

测定与红细胞相关的物理参数
发明的技术领域
本发明涉及一种基于数字全息术来测定与红细胞相关的物理参数的方法和设备。
背景技术
出于诊断的目的,通常使用的工具为全血细胞计数(CBC)。CBC可以分为用于分析的三个重要的细胞类型,即,白细胞(白血球)、红细胞(红血球)和血小板(凝血细胞)。通过测定红细胞的重要特征(例如平均红细胞体积(MCV)-红细胞的平均体积、红细胞分布宽度(RDW)-红细胞群变异的测量、平均红细胞血红蛋白(MCH)-每个红细胞的平均血红蛋白量、氧饱和度(SO2)-红细胞中氧所占据的血红蛋白结合位点的比例测量和平均红细胞血红蛋白浓度(MCHC)-红细胞中血红蛋白的平均浓度),可以诊断多种疾病。例如基于MCV的值是否高于或低于预期的正常范围,贫血被分为小红血球性贫血或巨红血球性贫血。目前,这些参数可以通过复杂且昂贵的系统来测量,例如流式细胞仪和共聚焦激光扫描显微镜,或仅能够测量一个参数的系统(例如阻抗体积分析仪)。
Rappaz等人的文章“Comparative Study of Human Erythrocytes byDigital Holographic Microscopy,Confocal Microscopy,and ImpedanceVolume Analyzer”(Journal of the International Society for Advancement ofCytometry,73A:895-903,2008)中描述的一项技术,说明了数字全息显微镜技术结合去耦过程对于折射率和完整的单个红血球体积的测量的适用性。Rappaz等人的文章“Measure ment of the integral refractive index anddynamic cell morphometry of living cells with digital holographic microscopy”(OPTICS EXPRESS,Vol.13,No.23,936-9373,2005)进一步对去耦过程做了更详细地解释,其中包含名为“去耦过程”的实验方案,其目的在于根据活细胞的定量相位图像分别测量整体折射率和细胞厚度。
因此,为了计算红细胞参数(例如MCV和/或MCHC和/或RDW和/或SO2和/或MCH),必须配制红细胞溶液,并将该红细胞溶液与第二溶液进行交换,其中这两种溶液的折射率必须是已知的。这种方法费力且耗时,对实验技能(例如红细胞不能受到污染)要求高。此外,交换溶液的过程是在测量中通过用户交互增加意外风险的另一步骤。事实上,交换溶液时的另一缺点是,存在于第一溶液样品中的所有红细胞在冲洗掉第一溶液并用第二溶液替换后将可能不再存在,这意味着将会有样品的丢失或污染。溶液交换过程中红细胞的减少,归因于在进行首次测量前红细胞需要沉淀并粘附在观察容器上。如果红细胞不能适当地粘附,它们将松开并随第一溶液冲洗掉和/或移动到新的位置并粘附在观察容器上。这将由于引入伪像,因此影响红细胞的重构波前的质量。此外,流体的交换和存在可能以不希望的方式影响红细胞的形态和/或预期的行为。观察容器也必须只能使用一次,以便不会对测量造成污染(例如通过稀释和/或混合溶液从而影响折射率使其成为未知而非假设的已知因数)。此外,使用所需精度在纳米范围内的一次性消耗品生产具有特定高度的观察容器是非常困难的。
因此,一种测定上述参数的改进方法是有利的,并且特别是更快、更容易和更可靠的方法是所期望的。
发明内容
因此,本发明寻求单独地或以任何组合的方式来减轻、缓和或消除一个或多个上述缺点。特别地,本发明的一个目的可以为在不需要交换溶液和/或存在溶液的情况下,测定红细胞的平均红细胞体积、红细胞分布宽度、平均红细胞血红蛋白、氧饱和度和平均红细胞血红蛋白浓度。
根据本发明的第一方面,上述及其他目的通过一种基于数字全息术来测定至少一个红细胞的参数的方法来解决,该方法包括以下步骤:
构建包括表示所述至少一个红细胞的波前的干涉图样,所述干涉图样由所述至少一个红细胞和电磁辐射的相互作用产生;
根据所述干涉图样对表示所述至少一个红细胞波前的振幅信息和相位信息进行重构;以及
通过对表示所述至少一个红细胞的相位信息和振幅信息进行组合,来测定所述至少一个红细胞的平均红细胞体积和/或平均红细胞血红蛋白浓度和/或氧饱和度和/或平均红细胞血红蛋白。
文章“Comparative Study of Human Erythrocytes by Digital HolographicMicroscopy,Confocal Microscopy,and Impedance Volume Analyzer”中公开了一种不涉及或使用振幅信息的方法。Rappaz等人的文章中所建议的方法涉及使用在使用不同的样品液体时获得的两个不同的图像,并利用这些图像的差异能够测定平均红细胞体积和平均红细胞血红蛋白。然而,本发明涉及一种将振幅信息和相位信息结合起来使用的方法,其创建了仅从一个图像(即不使用例如不同的细胞液体)测定参数的基础。本发明的此优点及其他优点在下面进一步公开。
通过本发明的方法,通过使用两个或多个波长的重构的振幅信息和/或相位信息来测定氧饱和度。
本发明的方法在不需要交换溶液和/或存在溶液的情况下,能够测定红细胞的平均红细胞体积、红细胞分布宽度、平均红细胞血红蛋白、氧饱和度和平均红细胞血红蛋白浓度。Rappaz等人的方法并非如此。另外,通过去除交换溶液的过程,显著降低了对用户所需的实验技能的要求,且由于用户处理不当而造成样品污染的风险也降低了。此外,不需要交换溶液,减少了测定所述参数的时间,原因在于根据本发明的方法仅需对干涉图样进行一次检测并对表示至少一个红细胞的振幅信息和相位信息进行重构来测定所述参数。本发明的方法从而可以以非常快的速率来研究现象,其中该速率主要取决于检测干涉图样的数字检测器的获取速率。
根据本发明,使用了振幅信息和相位信息二者。因此可以直接根据一个图像来确定每个像素的厚度进而折射率。由于折射率应仅取决于血红蛋白,根据本发明可以从血液中测定所需的参数,例如红细胞的平均红细胞体积、平均红细胞血红蛋白浓度、氧饱和度或平均红细胞血红蛋白。因此,根据本发明的一个特定实施例,对表示所述至少一个红细胞的相位信息和振幅信息进行组合,涉及测定包含至少一个红细胞的样品的折射率。根据又一实施例,所述方法通过仅创建一个用于测定所述参数的图像来实施。下面对此进行进一步讨论。
另外,本发明的方法去除了对一次性观察容器的需要和将红细胞粘附到所述观测容器上的需要。而且,在根据本发明的方法中,由于单次检测的要求,去除了在交换溶液时由于红细胞的可能移动或从观察容器中冲洗掉而引起的伪像。根据本发明的方法,由于使用了振幅信息和相位信息,不需要具有明确高度的观察容器和具有已知折射率的溶液。该方法由于能够促使流体结合至少一个红细胞,还允许至少一个红细胞存活于类似体内的环境中。
根据第二方面,提供了一种用于基于数字全息术来分析包含至少一个红细胞的样品的装置。该装置包括:至少一个光源,其被布置为创建至少一个物体光束和至少一个基准光束,其中所述至少一个物体光束和所述至少一个基准光束是相互相干的;用于将所述样品暴露于所述至少一个物体光束的工具;用于将通过了所述样品的所述至少一个物体光束和所述至少一个基准光束进行叠加从而产生干涉图样的装置;被布置为检测所述干涉图样的传感器;及处理单元,其被布置为根据所述干涉图样对表示所述至少一个红细胞的物体波前的相位信息和/或振幅信息进行重构,以及根据所述相位信息和/或振幅信息对平均红细胞体积和/或平均红细胞血红蛋白浓度和/或平均红细胞血红蛋白和/或氧饱和度进行测定
第一方面的特征和优点一般适用于第二方面。下面更详细地公开了本发明的其他特征和优点。
附图说明
图1是根据本发明的方法的流程图。
图2是在透明生物体通过前后的物体光束的波前的示意图。
图3是根据本发明的红细胞模型的一个实施例和方法的概念性示意图。
图4是根据本发明的一个实施例的装置的示意图。
本发明的优选实施方式的详细描述
如图1中的流程图所示,提供了一种基于数字全息术来测定至少一个红细胞的参数的方法,所述方法包含以下步骤:
·构建干涉图样,所述干涉图样包含表示所述至少一个红细胞的波前,所述干涉图样由所述至少一个红细胞和电磁辐射的相互作用产生;
·根据所述干涉图样对表示所述至少一个红细胞波前的振幅信息和相位信息进行重构;以及
·通过对表示所述至少一个红细胞的相位信息和振幅信息进行组合,来测定所述至少一个红细胞的平均红细胞体积和/或平均红细胞血红蛋白浓度和/或氧饱和度和/或平均红细胞血红蛋白。
根据本发明的方法的一个特定实施例,电磁辐射具有在可见光谱范围内的波长。根据本发明的方法的其他实施例,可以使用例如在红外和/或X射线和/或紫外光谱范围内的电磁辐射。电磁辐射由能够辐射至少一个波长的辐射源产生。例如,产生在光谱范围内的辐射的电磁辐射源可以是例如任何类型的激光源,例如二极管或He-Ne激光器。优选使用二极管激光器。
在根据本发明的方法的一个实施例中,相互相干的至少一个物体光束和至少一个基准光束,是通过将从相干光源发出的光束分成两个光束(例如通过分束器)生成的。物体光束和基准光束是相互相干的,这意味着它们在这一过程内具有相同的频率且具有恒定的相位关系。
物体光束穿过包含在观察容器中和/或上的至少一个红细胞,该观察容器适于包含所述至少一个红细胞且能够允许电磁辐射穿过所述观察容器和包含在观察容器中和/或上的所述至少一个红细胞。观察容器优选能够允许电磁辐射穿过观察容器且对辐射没有或具有很小的影响或具有已知的辐射干扰效应。观察容器还可以包含能够突出显示至少一个红细胞的所需物理特性(例如形态变化和/或折射率变化)的介质。观察容器还可以包含使至少一个红细胞保持在所需状态(例如在磷酸盐缓冲盐水和/或三羟甲基氨基甲烷缓冲盐水和/或Hank’s缓冲盐溶液和/或细胞生长培养基(例如在有或无补充剂的RPMI1640)的溶液中保持死的和/或活的和/或生长的和/或萎缩的状态)的介质。
基准光束不受包含在观察容器中和/或上的至少一个红细胞的影响,原因在于基准光束是沿与物体光束不同的光路引导的(例如通过反射镜或光纤)。
如图2所示,物体光束50在穿过包含在观察容器58中和/上的至少一个红细胞55之前具有已知的波前51。当物体光束50穿过至少一个红细胞55时,至少一个红细胞55将吸收物体光束的至少一部分强度,并且物体光束50由于其穿过至少一个红细胞55,与周围光束相比将经历一个光程长度差。由已知的物体光束波前51和至少一个红细胞55之间的相互作用产生的物体波前52将因此发生相移,且其振幅将受到来自至少一个红细胞55的吸收的影响。类似地,基准光束(未示出)也具有已知的波前,该波前因此不受至少一个红细胞的影响。光路被定义为物理/几何厚度乘以折射率。吸收被定义为通过红细胞由于波传播转换到其他类型能量的电磁辐射能。
在根据本发明的方法的一个实施例中,已穿过包含在观察容器中和/或上的至少一个红细胞的至少一个物体光束与至少一个基准光束的叠加,是通过将两个光束合并在一起(例如通过另一个分束器)来实现的。至少一个红细胞和电磁辐射的相互作用和叠加产生干涉图样,该干涉图样中包含关于表示至少一个红细胞的物体波前的信息。
在根据本发明的方法的一个实施例中,干涉图样通过数字传感器(例如电荷耦合器件(CCD)或互补金属氧化物半导体(CMOS))进行检测。所检测的干涉图样也称为全息图。
例如可以使用夫琅和费(Fraunhofer)或傅立叶(Fourier)或菲涅耳(Fresnel)装置,来叠加已穿过至少一个红细胞的至少一个物体光束和至少一个基准光束,从而生成干涉图样,并检测该干涉图样。优选使用菲涅耳装置。所述三种装置之间的差别在于光学配置满足允许表示所述至少一个红细胞的波前的重构近似值的某些条件。在菲涅耳装置中需满足的允许重构波前的近似值的条件是传感器与至少一个红细胞之间的距离与至少一个红细胞的尺寸和传感器的尺寸相比较大。这可通过使用显微镜物镜来实现,该显微镜物镜收集来自所述至少一个红细胞的散射光并将其以几乎平行的光束引导至所述数字传感器,从而形成一个光学结构,其中至少一个红细胞和数字传感器之间的距离与至少一个红细胞的尺寸和传感器的尺寸相比较大。
根据所检测的干涉图样,对物体波前的振幅信息和相位信息进行重构。重构通过任何常见的数值重构过程来实现,例如傅立叶变换重构和/或菲涅耳变换重构和/或角频谱重构和/或卷积重构和/或小波重构。表示所述至少一个红细胞的重构的振幅信息和相位信息可以在例如每个域、振幅或相位中进行表示,在至少一个红细胞信息的二维空间内表示为图像和/或图谱,和/或在三维空间内表示为3D表象。重构的信息优选表示为图像。
重构的信息还用于由计算机执行的图像处理,来测定有关至少一个红细胞的参数,例如平均红细胞体积、红细胞分布宽度、平均红细胞血红蛋白、氧饱和度和平均红细胞血红蛋白浓度。若干获得的图像也可以进一步用于图像处理以测定有关所研究的至少一个红细胞的其他有用信息。所述若干图像可以在不同的时间点获得并存储,以进行之后的图像处理和重构,这意味着根据本发明的方法与时间无关。这样的图像处理可能涉及到其中至少一个红细胞被定位的分割算法。这些被定位的至少一个红细胞可以进行进一步的形态分析,例如关于它们的宽度、形状等。可以使用的分割算法是分水岭分割和/或梯度分析,但可以使用用于检测和/或在图像中描绘至少一个红细胞的轮廓的任何类型的算法。优选使用分水岭分割。
或者,为了生成一个基准光束和一个物体光束,可以使用其他类型的数字全息术装置来产生干涉图样,例如在同轴数字全息术中利用了生成和使用一个光束。在研究本说明书时,如何修改本发明的方法以便将其用于例如同轴数字全息术中,对本领域技术人员来说是明显的。
当测定红细胞的参数(例如平均红细胞体积和/或红细胞分布宽度和/或和/或平均红细胞血红蛋白和/或氧饱和度和/或平均红细胞血红蛋白浓度)时,红细胞根据光学视图通常被模拟为,主要包含血红蛋白、盐和包含在厚度可忽略的透明细胞膜中的其他有机化合物的均匀的水溶液。水和血红蛋白分子占据了大于95%的红细胞体积。因此,对于使用前向光散射的技术(例如流式细胞仪、相差显微镜和数字全息术),红细胞的特征可表示为体积V和复合折射率nc=nr-ni。由于红细胞的内部几乎完全由水和血红蛋白占据,红细胞之间的折射率nc的变化可以只归因于血红蛋白浓度的变化。
通过应用上述模型,米式散射理论是适用的,并且包含在由电磁辐射和至少一个红细胞之间的相互作用产生的干涉图样中的物体波前可以通过数字传感器进行检测。这使得在重构表示至少一个红细胞的所述物体波前时,可以测定与红细胞相关的所述参数。此种常用的红细胞模型和米式散射理论使用中的问题在于,重构的波前取决于两个未知参数,体积V和复合折射率nc和其他已知参数。这意味着这个或一个重构的波前包含以非明显的方式混合的两个未知参数。
图3中示出了根据本发明的方法的概念性示意图。为了简单起见,下文将关于一个红细胞进行描述,但根据本发明的方法并不限于一个红细胞,而是优选至少一个红细胞。波前40在其穿过至少一个红细胞45之前是已知的,亦如图2(51,55)所示。在已知波前41(物体波前)的至少一部分之后已穿过所述至少一个红细胞45后,物体波前由于光程差(nery·l,与未穿过所述至少一个红细胞45的波前40的至少一部分的光程,nm·l相比)将发生相移,物体波前41与已知的波前40相比,也会受到来自至少一个红细胞的吸收Io的影响。
下面,将描述一种用于测定本发明所关注的参数的可行方法。
在根据本发明的方法的一个实施例中,红细胞的折射率被建模为:
nery=n0+αch  (1)
其中,n0是红细胞在缺乏血红蛋白时的折射率,α是在特定波长处血红蛋白的折射率增量,ch是血红蛋白浓度。此外,比尔-朗伯定律表明,吸收液体和非吸收液体的混合物的电磁辐射传输(A),与吸收液体的浓度(c)和所述电磁辐射穿过所述吸收液体的距离(l)成正比:
A = - log 10 I I 0 = log 10 I 0 - log 10 I = ∈ cl - - - ( 2 )
其中∈是所述吸收液体的摩尔吸收系数,I是从所述混合物透过的光强度,I0是所述混合物上的入射光强度。
假设红细胞内的血红蛋白浓度是恒定的和/或至少在通过数字传感器检测干涉图样的过程中是恒定的,该干涉图样包含表示至少一个红细胞的波前,并且由于仅由血红蛋白引起吸收,红细胞的总吸收(Ω)包含在表示所述红细胞的波前中,进一步包含在干涉图样中并通过数字传感器进行检测并重构到例如以振幅图像表示的振幅信息中,该红细胞的总吸收(Ω)表示为:
Ω = Σ i N p A i = ∈ h c h Σ i N p l i - - - ( 3 )
其中Np是表示所述重构的振幅图像中的所述红细胞的像素总数,∈h是在波长λ处的血红蛋白的摩尔吸光系数,ch是血红蛋白浓度,li是在像素i处所述红细胞的厚度。
如上所述,该模型取决于波长,并且可以在等吸光的波长处(例如在波长525-575nm处)使用,其中血红蛋白的摩尔吸光系数(∈h)不取决于血红蛋白的饱和度。可选地或另外地,在所使用的波长处的摩尔吸光系数∈h取决于血红蛋白的饱和度和进而红细胞时,总吸收被建模为:
Ω = S ∈ λ 0 c h Σ i N p l i + ( 1 - S ) ∈ λ d c h Σ i N p l i - - - ( 4 )
其中是在波长λ处氧合血红蛋白的摩尔吸光系数,S是血红蛋白氧饱和度,及
Figure BDA00003619277700096
是在波长λ处脱氧血红蛋白的摩尔吸光系数。
通过使用等式(1),在像素i处的相对相移
Figure BDA00003619277700097
包含在干涉图样中并通过数字传感器检测并重构到以相位图像表示的相位信息中,该相对相移
Figure BDA00003619277700098
表示为:
Figure BDA00003619277700094
其中nm是在像素i处非红细胞周围的流体的折射率,红细胞(φ)的总相移表示为:
Figure BDA00003619277700101
其中a是通过像素i成像的区域,
Figure BDA00003619277700102
是红细胞体积。通过使用根据等式(3)得出的血红蛋白的总吸收和根据等式(5)得出的总相移,至少一个红细胞的总相移还可表示为:
φ = 2 π ( n 0 - n m ) aλ V + 2 πα λ ∈ h Ω - - - ( 7 )
通过重新安排各项,至少一个红细胞的体积因此可表示为:
V = aλ 2 π ( n 0 - n m ) φ - aα ( n 0 - n m ) ∈ h Ω - - - ( 8 )
应认识到的是,根据本发明的上述方法,使用至少一个红细胞的相移和吸收的组合,红细胞平均体积MCV可由下式测定:
MCV = a N rbc ( n 0 - n m ) [ λ 2 π φ - α ∈ h Ω ] - - - ( 9 )
其中Nrbc是红细胞总数。
此外,像素中的相移通过结合等式(2)和(5)可以表示为:
Figure BDA00003619277700106
求解该式来计算像素厚度,得出:
Figure BDA00003619277700107
将其代入等式(2)中,求出像素中的血红蛋白浓度,得出:
Figure BDA00003619277700108
因此,平均红细胞血红蛋白浓度(MCHC)可以计算为:
Figure BDA00003619277700109
其中Np是红细胞像素的总数。
此外,通过使用等式3和4中的总吸收和等式6中对红细胞体积的定义,红细胞的氧饱和度被表示为:
S = [ ( Ω λ Ω λi ) ϵ λi - ϵ λ d ] / ( ϵ λ o - ϵ λ d ) - - - ( 14 )
(13)其中λi是等吸收光波长。这意味着氧饱和度可以通过对两个波长处(即等吸收光波长和非等吸收光波长)的振幅信息和/或相位信息进行重构来测定。或者,氧饱和度可以通过使用等式4在两个非等吸收光波长处进行测定。因此,氧饱和度被表示为:
S = ϵ λ 1 d - ϵ λ 2 d ( Ω λ 1 / Ω λ 2 ) ( ϵ λ 2 o - ϵ λ 2 d ) ( Ω λ 1 / Ω λ 2 ) - ( ϵ λ 1 o / ϵ λ 1 d ) - - - ( 15 )
其中
Figure BDA00003619277700113
Figure BDA00003619277700114
是在第一波长处脱氧血红蛋白和氧合血红蛋白的摩尔吸光系数,
Figure BDA00003619277700115
是在第二波长处脱氧血红蛋白和氧合血红蛋白的摩尔吸光系数。
另外地或可选地,平均红细胞血红蛋白MCH通过以下等式测定:
MCH=MCV·MCHC  (16)
另外地或可选地,红细胞分布宽度RDW可以通过任何常用的边缘检测算法来测定,该边缘检测算法被应用于重构的振幅图像和/或相位图像来测定红细胞的宽度并因此测定红细胞宽度的变化。
在本发明中应当注意的是,用于测定例如MCH和/或MCHC的上述方法仅构成根据本发明的一个实施例。存在可以执行来获得根据本发明的所关注参数的可选的衍生实施例。例如,通过首先计算每个血细胞的平均浓度,然后使用所有推导的平均浓度的平均值,可以得到MCHC。同样,这也仅是根据本发明的一个可能的实施例,然而也可以根据本发明使用其他方法、推导和/或公式。
现在将参考图4对根据本发明的一个实施例的装置进行描述。图4是用于分析样品10的装置100的示意图。装置100包含至少一个光源300、传感器500、分束器7a和7b、反射表面9a和9b,以及处理单元13。分束器7a和7b与反射表面9a和9b一起用作将样本暴露于至少一个物体光束21的工具和将至少一个物体光束21与至少一个基准光束23进行叠加的工具。
光源300被布置成产生至少一个相干光束19。该相干光束19可以是例如由任何种类的激光源(如发射波长为635nm的光的二极管激光器)发射的激光束。这里仅示出一个光源。然而,一般来说,装置100可以包含若干个可以同时使用的光源300。如果使用多个光源300,这些光源优选产生不同波长的光。这可能是有利的,原因在于样品10对不同波长的光的反应不同。因此,通过使用具有不同波长的多个光源300,可以获得样品的更多信息。光源300发射的相干光束19被引导至分束器7a。该分束器7a将至少一个相干光束19分为至少一个物体光束21和至少一个基准光束23。
采用这样的结构,至少一个物体光束21和至少一个基准光束23变得相互相干,这意味着它们在这一过程期间具有相同的频率且具有恒定的相位关系。如果存在多于一个光源300,分束器7a可以将由光源300发射的每个光束19分成相互相干的物体光束21和基准光束23。装置100在使用时,将包含至少一个红细胞的样品10布置在至少一个物体光束21的光路中。例如,反射表面9a可以用于重定向至少一个物体束21,使得样品10被暴露于至少一个物体光束21。
当物体光束21入射到样品10时,物体光束21将穿过样品10,并且特别是穿过至少一个红细胞11。由于与未穿过样本的光束(例如基准光束23)相比,物体光束21穿过样品10,由于折射率和吸收的差异,因而光程和强度存在差异。这又将导致物体光束21和基准光束23之间的相移差和振幅差。光路长度被定义为物理/几何厚度与折射率的乘积,吸收被定义为通过波传播转换到其他类型能量的电磁辐射能。反射表面9b可以被布置为反射并由此重定向至少一个基准光束23。具体地,可以对反射表面9a和9b进行布置,使得至少一个物体光束21和至少一个基准光束23被引导至用于叠加的装置,在此是分束器7b的形式。
分束器7b对至少一个物体光束21和至少一个基准光束23进行叠加,并将叠加的光束25引导至传感器500。传感器500被布置成检测从物体光束21和基准光束23产生的干涉图样。由于物体光束21和基准光束23是相互相干的,它们将在传感器500处产生干涉图样。特别地,由于至少一个物体光束21穿过样品10,至少一个物体光束21和至少一个基准光束23传播不同的光程并具有不同的强度,则干涉图样表示物体光束21和基准光束23之间的相移和强度。
传感器500例如可以是数字传感器(例如CCD(电荷耦合器件))或CMOS(互补金属氧化物半导体)图像传感器。传感器500还可以可操作地连接到处理单元13,处理单元13可以包含执行任何常见的重构过程的软件和/或硬件。重构过程可以根据通过传感器500检测的干涉图样对相位信息和/或振幅信息进行重构。重构的信息例如可以用于获得所研究的至少一个红细胞11的至少一个图像。例如,该信息也可以用于测定至少一个红细胞11的形态参数。具体地,处理单元13被布置为根据相位信息和/或振幅信息来测定参数(例如MCH、MCV、SO2、MCHC和RDW)。换言之,处理单元13可以被布置为执行根据本发明的实施例的方法的任何数据处理步骤。
此外,装置100可以包含存储介质或存储器(未示出)。该存储介质或存储器可以可操作地连接到处理单元13。具体地,该存储介质或存储器可以被布置为存储分析图像和与在分析图像中检测到变化的时间点有关的时间数据。装置100可以适于在环境舱15的内部操作。在这样的布置中,样品10不必从保存该样品的环境舱中取出来进行分析。如图4所示,光源300、传感器500和光学部件(例如分束器7a和7b以及反射表面9a和9b)包含在环境舱15的内部。
处理单元13可以位于环境舱15的内部,也可以位于其外部。如果处理单元位于环境舱15的外部,处理单元可以优选与传感器500进行无线通信。例如,该装置可以包含发送器(未示出),该发送器被布置为从传感器500发送与干涉图样有关的信息到与处理单元13可操作地连接的接收器上。发送器可以构成收发器的一部分,该收发器被布置为发送信号到处理单元13和接收来自处理单元13的信号。可选地,装置100可以包含环境舱15。在这种情况下,环境舱15优选与装置100形成一体。这样的设置提供了一种非常紧凑和灵活的方案。例如,包含一体形成的环境舱的装置100可以放置在工作台上,还可以在需要时方便地移动到另一个位置。有利的是,通过使环境舱15与装置100形成一体,并非装置的所有部件都必须被布置在环境舱的内部。例如,一些对环境舱的内部条件敏感的光学部件(例如分束器7a和7b以及反射表面9a和9b)可以被布置在环境舱的外部。
在一个实施例中,仅盛放样品10的样品架位于环境舱15的内部。装置100还可以包含控制器(未示出),用于控制环境舱15内部的湿度、温度和二氧化碳百分比中的至少一项。这样,可以实现针对至少一个红细胞的最佳条件。红细胞需要某些严格控制的参数。首先,介质需要保持pH为7.4。由于介质通常使用碳酸盐进行缓冲,这通常通过供应5%-10%的CO2到空气中来实现,尽管可以使用例如HEPES的有机缓冲剂来代替。为了阻止介质蒸发,空气应该是水饱和的。例如,相对湿度可以约为90%-95%或以上。由于细胞物种的依赖性(例如大多数哺乳动物细胞的适宜温度为37℃,昆虫细胞为20℃,而鸟类细胞为40℃~42℃),红细胞优选保持在37℃下。此外,对装置100,尤其是装置100的光学部件(例如分束器7a和7b以及反射表面9a和9b)可以使用不同的方案,以应对环境舱15的内部条件。
根据一个方案,装置100被布置在环境舱15的内部,而培养箱内部温度基本上等于室温。然后,环境舱的内部温度低速增加。这样,可以避免在装置100的光学部件上的冷凝,因为在环境舱15内部的光学部件与空气之间不会产生温度差。根据另一方案,光学部件(例如分束器7a和7b和反射表面9a和9b)可以被布置和/或处理,以能够应对环境舱15内部的条件。具体地,可以对光学部件进行处理以避免冷凝。这样,在环境舱15内部的温度并非基本上等于室温时,装置100可以被布置在环境舱15的内部。因此,可以避免上述进行缓慢升温的过程。
虽然已经结合特定的实施例对本发明进行了描述,但其目的并非限定这里所阐述的特定形式。相反,本发明的范围仅由所附权利要求书进行限定。在权利要求中,术语“包含”不排除存在其他元件或步骤。另外,虽然单个特征可以包含在不同的权利要求中,这并不意味着这些特征的组合是不可行的和/或不利的。另外,单数引用并不排除多个。这样,对“一”、“第一”、“第二”的引用并不排除多个。

Claims (17)

1.一种基于数字全息术来测定至少一个红细胞的参数的方法,所述方法包括以下步骤:
·构建包含表示所述至少一个红细胞的波前的干涉图样,所述干涉图样由所述至少一个红细胞和电磁辐射的相互作用产生;
·根据所述干涉图样对表示所述至少一个红细胞波前的振幅信息和相位信息进行重构;以及
·通过对表示所述至少一个红细胞的相位信息和振幅信息进行组合来测定至少一个红细胞的平均红细胞体积和/或平均红细胞血红蛋白浓度和/或氧饱和度和/或平均红细胞血红蛋白。
2.根据权利要求1所述的方法,其中对表示所述至少一个红细胞的相位信息和振幅信息进行组合涉及测定包含所述至少一个红细胞的样品的折射率。
3.根据权利要求1或2所述的方法,其中所述方法通过仅创建一个图像来实施,所述一个图像用于通过对表示所述至少一个红细胞的相位信息和振幅信息进行组合来测定至少一个红细胞的平均红细胞体积和/或平均红细胞血红蛋白浓度和/或氧饱和度和/或平均红细胞血红蛋白。
4.根据权利要求1至3中任一项所述的方法,其中所述平均红细胞体积通过以下等式来确定:
MCV = a N rbc ( n 0 - n m ) [ λ 2 π φ - α ∈ h Ω ] .
5.根据权利要求1至4中任一项所述的方法,其中所述平均红细胞血红蛋白浓度通过以下等式来确定:
Figure FDA00003619277600012
6.根据权利要求1至5中任一项所述的方法,其中所述平均红细胞血红蛋白通过以下等式来确定:
MCH=MCV·MCHC。
7.根据权利要求1至6中任一项所述的方法,其中氧饱和度通过以下等式来确定:
S = [ ( Ω λ Ω λi ) ϵ λi - ϵ λ d ] / ( ϵ λ o - ϵ λ d ) .
8.根据权利要求1至7中任一项所述的方法,其中来自至少两个波长的表示至少一个红细胞的重构的振幅信息和/或相位信息用于测定氧饱和度。
9.根据权利要求1至8中任一项所述的方法,其中氧饱和度(S)通过以下等式来确定:
S = ϵ λ 1 d - ϵ λ 2 d ( Ω λ 1 / Ω λ 2 ) ( ϵ λ 2 o - ϵ λ 2 d ) ( Ω λ 1 / Ω λ 2 ) - ( ϵ λ 1 o / ϵ λ 1 d ) .
10.根据权利要求1至9中任一项所述的方法,其中,所述电磁辐射具有在可见光谱范围内的波长。
11.一种用于通过数字全息术分析包含至少一个红细胞的样品的装置,所述装置包括:
至少一个光源,所述至少一个光源被布置为产生至少一个物体光束和至少一个基准光束,其中所述至少一个物体光束和所述至少一个基准光束是相互相干的;
·用于将所述样品暴露于所述至少一个物体光束的工具;
·用于对通过了所述样品的所述至少一个物体光束和所述至少一个基准光束进行叠加从而产生干涉图样的工具;
·被布置为检测所述干涉图样的传感器;及
·处理单元,所述处理单元被布置为根据所述干涉图样对表示所述至少一个红细胞的物体波前的相位信息和/或振幅信息进行重构,以及根据所述相位信息和/或振幅信息对平均红细胞体积和/或平均红细胞血红蛋白和/或氧饱和度和/或平均红细胞血红蛋白浓度进行测定。
12.根据权利要求11所述的装置,其中所述处理单元包括对红细胞分布的测定。
13.根据权利要求11或12所述的装置,其中所述处理单元包括通过使用至少两个波长的、表示至少一个红细胞的重构的相位信息和/或振幅信息对氧饱和度的测定。
14.根据权利要求13所述的装置,其中使用了至少两个波长。
15.根据权利要求11至14中任一项所述的装置,还包括环境舱,其中所述至少一个光源、所述用于暴露的工具、所述用于叠加的工具以及所述传感器位于所述培养箱内部。
16.根据权利要求11至15中任一项所述的装置,其中所述处理单元位于所述环境舱外部,并且所述装置还包括被布置为将所述干涉图样从所述传感器传送到所述处理单元的发送器。
17.根据权利要求11至16中任一项所述的装置,还包括控制器,所述控制器用于控制所述环境舱内部的湿度、温度和二氧化碳百分比中的至少一项。
CN201280007685.1A 2011-02-16 2012-02-15 测定与红细胞相关的物理参数 Active CN103620509B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201161443487P 2011-02-16 2011-02-16
SE1150125 2011-02-16
US61/443,487 2011-02-16
SE1150125-1 2011-02-16
PCT/SE2012/050161 WO2012112114A1 (en) 2011-02-16 2012-02-15 Determination of physical parameters associated with erythrocytes

Publications (2)

Publication Number Publication Date
CN103620509A true CN103620509A (zh) 2014-03-05
CN103620509B CN103620509B (zh) 2016-08-10

Family

ID=46672836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280007685.1A Active CN103620509B (zh) 2011-02-16 2012-02-15 测定与红细胞相关的物理参数

Country Status (5)

Country Link
US (1) US9423342B2 (zh)
EP (1) EP2676173B1 (zh)
JP (1) JP6001566B2 (zh)
CN (1) CN103620509B (zh)
WO (1) WO2012112114A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108426994A (zh) * 2014-06-16 2018-08-21 西门子医疗保健诊断公司 分析数字全息显微术数据以用于血液学应用
CN109564209A (zh) * 2016-05-11 2019-04-02 思迪赛特诊断有限公司 对样品实施的光学测量
CN110646377A (zh) * 2019-09-27 2020-01-03 京东方科技集团股份有限公司 血液检测系统及其检测方法
CN110799829A (zh) * 2017-07-06 2020-02-14 特拉维夫大学拉莫特有限公司 用于在环境室中的生物细胞样品的三维无标记光学成像的系统和方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10614287B2 (en) 2014-06-16 2020-04-07 Siemens Healthcare Diagnostics Inc. Virtual staining of cells in digital holographic microscopy images using general adversarial networks
FR3034196B1 (fr) * 2015-03-24 2019-05-31 Commissariat A L'energie Atomique Et Aux Energies Alternatives Procede d'analyse de particules
WO2017046806A1 (en) * 2015-09-17 2017-03-23 Technion Research & Development Foundation Limited Reflectance confocal microscopy of blood cells
JP6951738B2 (ja) * 2017-07-10 2021-10-20 国立大学法人千葉大学 ホログラムデータ作成プログラム及びホログラムデータ作成方法
JP7194030B2 (ja) * 2019-01-30 2022-12-21 株式会社エビデント 赤血球モニタリング装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6262818B1 (en) * 1998-10-07 2001-07-17 Institute Of Applied Optics, Swiss Federal Institute Of Technology Method for simultaneous amplitude and quantitative phase contrast imaging by numerical reconstruction of digital holograms
US20030033102A1 (en) * 2001-08-09 2003-02-13 Thomas Dietiker System and method for a self-calibrating non-invasive sensor
WO2006090320A1 (en) * 2005-02-23 2006-08-31 Lyncee Tec S.A. Wave front sensing method and apparatus
CN101346673A (zh) * 2005-12-22 2009-01-14 相位全息成像Phi有限公司 用于分析细胞的样本的方法和装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005114678A (ja) * 2003-10-10 2005-04-28 Hamamatsu Photonics Kk 散乱吸収体計測装置及び計測方法
JP4639321B2 (ja) * 2005-11-14 2011-02-23 コニカミノルタセンシング株式会社 生体情報測定装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6262818B1 (en) * 1998-10-07 2001-07-17 Institute Of Applied Optics, Swiss Federal Institute Of Technology Method for simultaneous amplitude and quantitative phase contrast imaging by numerical reconstruction of digital holograms
US20030033102A1 (en) * 2001-08-09 2003-02-13 Thomas Dietiker System and method for a self-calibrating non-invasive sensor
WO2006090320A1 (en) * 2005-02-23 2006-08-31 Lyncee Tec S.A. Wave front sensing method and apparatus
CN101346673A (zh) * 2005-12-22 2009-01-14 相位全息成像Phi有限公司 用于分析细胞的样本的方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
RAPPAZ ET AL.: "Comparative Study of Human Erythrocytes by Digital Holographic Microscopy,Confocal Microscopy,and Impedance Volume Analyzer", 《CYTOMETRY.PART A:JOURNAL OF THE INTERNATIONAL SOCIETY FOR ADVANCEMENT OF CYTOMETRY》, vol. 73, 31 October 2008 (2008-10-31), pages 895 - 903, XP055032058, DOI: doi:10.1002/cyto.a.20605 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108426994A (zh) * 2014-06-16 2018-08-21 西门子医疗保健诊断公司 分析数字全息显微术数据以用于血液学应用
CN109564209A (zh) * 2016-05-11 2019-04-02 思迪赛特诊断有限公司 对样品实施的光学测量
CN110799829A (zh) * 2017-07-06 2020-02-14 特拉维夫大学拉莫特有限公司 用于在环境室中的生物细胞样品的三维无标记光学成像的系统和方法
CN110646377A (zh) * 2019-09-27 2020-01-03 京东方科技集团股份有限公司 血液检测系统及其检测方法
CN110646377B (zh) * 2019-09-27 2022-07-05 京东方科技集团股份有限公司 血液检测系统及其检测方法

Also Published As

Publication number Publication date
CN103620509B (zh) 2016-08-10
EP2676173B1 (en) 2023-03-29
EP2676173A4 (en) 2017-11-08
US20140113323A1 (en) 2014-04-24
JP6001566B2 (ja) 2016-10-05
US9423342B2 (en) 2016-08-23
WO2012112114A9 (en) 2012-11-22
JP2014512514A (ja) 2014-05-22
EP2676173A1 (en) 2013-12-25
WO2012112114A1 (en) 2012-08-23

Similar Documents

Publication Publication Date Title
CN103620509A (zh) 测定与红细胞相关的物理参数
KR102479862B1 (ko) 입자 분석 방법
Liebi et al. Nanostructure surveys of macroscopic specimens by small-angle scattering tensor tomography
JP6559555B2 (ja) 光計測方法および装置
Kirshner et al. 3‐D PSF fitting for fluorescence microscopy: implementation and localization application
Baczewska et al. Refractive index changes of cells and cellular compartments upon paraformaldehyde fixation acquired by tomographic phase microscopy
US20200300768A1 (en) Determination device, determination method, and determination program
US9581429B2 (en) Method and system for reconstructing optical properties of diffracting objects immersed in a liquid medium
Carpignano et al. Low-coherence reflectometry for refractive index measurements of cells in micro-capillaries
Liao et al. Polarization measurements and evaluation based on multidimensional polarization indices applied in analyzing atmospheric particulates
Brenner et al. Optical coherence tomography images simulated with an analytical solution of Maxwell’s equations for cylinder scattering
KR101563810B1 (ko) 공통 광 경로 3차원 회절 광 현미경 및 그의 시편 변형성 측정 방법
Schrangl et al. Kinetic analysis of single molecule FRET transitions without trajectories
EP1729261A1 (en) Method for tomographic reconstruction
JP7228189B2 (ja) 細胞毒性を評価する方法及び装置
Roselli et al. New data‐driven method from 3D confocal microscopy for calculating phytoplankton cell biovolume
Wu et al. Comparison of high-speed polarization imaging methods for biological tissues
Aharon et al. Second harmonic generation for moisture monitoring in dimethoxyethane at a gold-solvent interface using plasmonic structures
Merolle et al. Impact of Sample Preparation Methods on Single-Cell X-ray Microscopy and Light Elemental Analysis Evaluated by Combined Low Energy X-ray Fluorescence, STXM and AFM
US20240102912A1 (en) Plaque counting assay method
Zhang et al. Surface plasmon resonance microscopy based on total internal reflection
Hevisov et al. Confocal Laser Scanning Microscope Imaging of Custom-Made Multi-Cylinder Phantoms: Theory and Experiment
東京大学大学院理学系研究科物理学専攻 Mid-infrared photothermal quantitative phase imaging
EP4293342A1 (en) Automatic analyzer and method for optically analyzing a biological sample
KR20240107127A (ko) 샘플 중의 나노입자 특성들을 결정하는 방법 및 장치

Legal Events

Date Code Title Description
PB01 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