CN107356515B - 一种荧光图像的荧光强度确定方法和系统 - Google Patents
一种荧光图像的荧光强度确定方法和系统 Download PDFInfo
- Publication number
- CN107356515B CN107356515B CN201710596449.2A CN201710596449A CN107356515B CN 107356515 B CN107356515 B CN 107356515B CN 201710596449 A CN201710596449 A CN 201710596449A CN 107356515 B CN107356515 B CN 107356515B
- Authority
- CN
- China
- Prior art keywords
- fluorescence
- value
- target
- image
- fluorescence image
- 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
- 238000002073 fluorescence micrograph Methods 0.000 title claims abstract description 157
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000001514 detection method Methods 0.000 claims abstract description 106
- 238000000799 fluorescence microscopy Methods 0.000 claims abstract description 35
- 238000007621 cluster analysis Methods 0.000 claims abstract description 34
- 238000000605 extraction Methods 0.000 claims abstract description 24
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 230000011218 segmentation Effects 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 11
- 238000000926 separation method Methods 0.000 claims description 11
- 238000003708 edge detection Methods 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 8
- 238000012632 fluorescent imaging Methods 0.000 claims description 2
- 238000005206 flow analysis Methods 0.000 abstract description 9
- 238000003384 imaging method Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 10
- 238000001917 fluorescence detection Methods 0.000 description 5
- 230000001186 cumulative effect Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000000684 flow cytometry Methods 0.000 description 2
- 238000000265 homogenisation Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 238000003703 image analysis method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
-
- 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—Electro-optical investigation, e.g. flow cytometers
- G01N15/1429—Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its signal processing
- G01N15/1431—Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its signal processing the electronics being integrated with the analyser, e.g. hand-held devices for on-site investigation
-
- G01N15/1433—
-
- 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—Electro-optical investigation, e.g. flow cytometers
- G01N15/1434—Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its optical arrangement
-
- 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—Electro-optical investigation, e.g. flow cytometers
- G01N15/1434—Electro-optical investigation, e.g. flow cytometers using an analyser being characterised by its optical arrangement
- G01N2015/144—Imaging characterised by its optical setup
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Electro-optical investigation, e.g. flow cytometers
- G01N2015/1493—Particle size
Abstract
一种荧光图像的荧光强度确定方法和系统，方法包括：采用显微荧光成像系统对目标样品进行明场成像和荧光成像得到明场图像和荧光图像；对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，得到所述荧光图像中每个检测目标对应的荧光图像区域；计算得到每个检测目标的荧光图像区域的累计灰度值、最大灰度值、平均灰度值、每个检测目标的明场图像区域的直径值，基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析。可以开展或输入到流式聚类分析软件进行分析，并获得与流式分析术对应和类似的结果。
Description
技术领域
本发明涉及细胞检测技术领域，具体涉及一种荧光图像的荧光强度确定方法和系统。
背景技术
流式分析技术是一种重要的细胞分析技术，其核心的原理就是统计检测目标(细胞或其他检测颗粒)的荧光强度进行聚类分析。流式分析技术是通过记录荧光信号的电压脉冲值来表达荧光强度的。现在，随着显微荧光成像技术和智能识别技术的发展，通过对大量检测目标进行荧光成像，并通过图像识别方法进行分析，可以实现基于荧光成像技术的聚类分析，这种以显微成像分析为基础融合流式聚类分析的方法(图像类流式分析方法)，成为一种新的分析方法。
但是，基于显微荧光成像和智能识别的图像分析是以像素和灰度为基础来表达荧光强度的，与流式分析以电压脉冲来表达荧光强度的方式完全不同。如何建立者两者之间的对应关系，让图像表达的荧光强度与流式表达的荧光强度具有可比性，并保证进行类流式分析结果的准确性，成为图像类流式分析方法中亟待解决的关键点和难点。
发明内容
有鉴于此，本发明实施例提供一种荧光图像的荧光强度确定方法和系统，以使得图像表达的荧光强度与流式表达的荧光强度之间具有可比性。
为实现上述目的，本发明实施例提供如下技术方案：
对应权利要求
基于上述技术方案，本发明实施例提供的荧光图像的荧光强度确定方法和系统，通过采用显微荧光成像系统对目标样品进行荧光成像得到荧光图像；对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，得到所述荧光图像中每个检测目标对应的荧光图像区域；计算得到每个检测目标的荧光图像区域的累计灰度值、最大灰度值、平均灰度值、每个检测目标的明场图像区域的直径值作，基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析。可以开展或输入到流式聚类分析软件进行分析，并获得与流式分析术对应和类似的结果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本申请实施例公开的一种荧光图像的荧光强度确定方法的流程示意图；
图2为本申请实施例公开的一种荧光图像的荧光强度确定方法中智能识别处理的流程示意图；
图3为本申请另一实施例公开的一种荧光图像的荧光强度确定方法中智能识别处理的流程示意图；
图4为本申请再一实施例公开的一种荧光图像的荧光强度确定方法中智能识别处理的流程示意图；
图5为本申请又一实施例公开的一种荧光图像的荧光强度确定方法中智能识别处理的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在实际的应用中，流式分析术提供了多种指标供客户选择，让客户可以从不同侧面去对目标细胞进行分析，主要指标有：1、以细胞通过时采集到的对应侧向荧光通道的累计电压脉冲值作为累计荧光强度。2、以细胞通过时采集到的对应侧向荧光通道的最强电压脉冲作为最高荧光强度。3、以细胞通过时采集到的对应侧向荧光通道的平均电压脉冲值作为平均荧光强度。4、细胞通过时的正向散射通道采集到的电压脉冲值作为目标的大小值。。
在细胞荧光图像分析中，对目标细胞的荧光图像进行边缘提取和分割后，可以分析统计出该细胞的荧光像素面积、荧光图像直径、荧光图像长轴、荧光图像短轴、荧光图像周长、荧光图像圆度、荧光图像最大灰度值、最小灰度值、平均灰度值、累计灰度值、中值灰度值、灰度值落差值(最大灰度值减最小灰度值)……等几十种特征指标。
但是，细胞荧光图像分析虽然提取了非常多的图像特征值，哪种特征值可以与上述流式细胞术中的累计电压脉冲值、高度和宽度具有对应关系，其统计分析数据可以得到类比于流式分析的结果和精度，是现有图像分析方法没有解决的问题。这阻碍了图像分析方法与流式分析方法结果的可比性。
针对于此，申请人通过进行大量实验验证得到4个可以与流式细胞术荧光强度对接的指标。分别是：累计荧光灰度值与流式对应荧光通道的累计电压脉冲值对应；单细胞最高荧光灰度值与流式对应荧光通道的最强电压脉冲值对应；单细胞平均荧光灰度值与流式对应荧光通道的平均电压脉冲值对应；单细胞明场直径与流式正向散射通道采集到的电压脉冲值对应。
针对于现有技术中图像表达的荧光强度与流式表达的荧光强度不具有可比性，导致用户无法更好地判断检测结果的准确性的问题，本申请公开了一种基于上述四种对应关系的荧光图像的荧光强度确定方法，参见图1，该方法包括：
步骤S101：采用显微荧光成像系统对目标样品进行荧光成像得到荧光图像；
在本步骤中，首先采用显微荧光成像系统对目标样品进行图像类流式荧光检测，获得目标样品的荧光图像；
步骤S102：确定荧光图像区域；
在得到目标样品的荧光图像之后，还可以采用智能识别系统对所述荧光图像进行图像处理与分析，得到所述荧光图像中所包含的每个细胞(检测目标)的荧光图像区域，具体的，在图像处理与分析过程中通过对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，即可得到所述荧光图像中每个检测目标对应的荧光图像区域；
步骤S103：计算得到每个检测目标的荧光图像区域的累计灰度值；
在本步骤中，依据得到的单个检测目标的荧光图像区域的区域范围，通过图像识别技术，计算得到每个检测目标的荧光图像区域的累计灰度值，将所述累计灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的累计电压脉冲值表达的累计荧光强度值；所述累计荧光灰度值为单个检测目标的荧光图像区域内所有的像素灰度值的总和；
步骤S104：计算得到每个检测目标的荧光图像区域的最大灰度值；
在本步骤中，依据得到的单个检测目标的荧光图像区域的区域范围，通过图像识别技术，计算得到每个检测目标的荧光图像区域的最大灰度值，将所述最大灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的最大电压脉冲值表达的最大荧光强度值；所述最大灰度值为单个检测目标的荧光图像区域内灰度值最高的像素点的灰度值；
步骤S105：计算得到每个检测目标的荧光图像区域的平均灰度值；
在本步骤中，依据得到的单个检测目标的荧光图像区域的区域范围，通过图像识别技术，计算得到每个检测目标的荧光图像区域的平均灰度值；将所述平均灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的平均电压脉冲值表达的平均荧光强度值；所述平均灰度值为单个检测目标的荧光图像区域内所有像素点的平均值；
步骤S106：计算每个检测目标的明场图像区域的直径值作为目标大小值；
在本步骤中，获取每个检测目标的明场图像，对所述明场图像进行数据分析，计算得到每个检测目标的明场图像区域的直径值作为目标大小值，并确定该值作为开展流式聚类分析中正向散射通道以光电倍增管采集到的电压脉冲值表达的检测目标的大小值的对应替代指标；
上述方法中，步骤S103-106之间的顺序可以依据用户需求自行调整；
步骤S107：基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析；
在本步骤中，将所述累计灰度值作为累计荧光强度值、将所述最大灰度值作为最大荧光强度值、将所述平均灰度值作为平均荧光强度值；将所述明场图像区域的直径值作为检测目标的大小值，采用累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析，生成与累计灰度值、最大灰度值、平均灰度值和明场直径值中的一项或多项对应的荧光检测波形图，将生成的荧光检测波形图可直接与采用常规流式分析术得到的荧光检测波形图进行直接对比，依据对比结果，调整流式聚类分析中检测目标的大小值，采用累计灰度值、最大灰度值、平均灰度值和明场直径值的权重值，从而使得其借测结果更接近于现有的采用常规流式分析术得到的荧光检测波形，从而提高检测结果的准确性。
具体的，在实际使用过程中，在本申请上述实施例公开的方法中无需通过光电倍增管获得所述检测目标对应的累计电压脉冲值、最大电压买充值、平均电压买充值和检测目标的大小值，即可对检测目标进行流式分析，且具有较高的精确度。
在本申请上述实施例公开的技术方案中，基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析时，可采用所述累计灰度值、最大灰度值、平均灰度值和明场直径值直接进行流式聚类分析，当然，在采用累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析之前，需要对流式聚类分析中的公式进行训练，调整公式中的权重值，以提高计算结果的精准度。当然在本步骤中也可以依据预设的映射表查找得到与所述累计灰度值、最大灰度值、平均灰度值和明场直径值一一对应的累计荧光强度值、最大荧光强度值、平均荧光强度值和荧光面积值，所述预设映射表中存储有预设的累计灰度值与累计荧光强度值之间的映射关系、最大灰度值与最大荧光强度值之间的映射关系、平均灰度值与平均荧光强度值之间的映射关系、明场直径值与荧光面积值之间的映射关系。
在本申请上述实施例公开的技术方案中，可以采用智能识别系统对目标样品进行荧光成像得到的荧光图像进行智能识别，通过对荧光预想中每个检测目标的荧光图像区域进行边缘提取和分割的方式得到每个检测目标的荧光图像区域，参见图2，其过程具体可以包括：
S1021：对得到的荧光图像进行高斯平滑滤波；
在本步骤中，首先对荧光图像进行高斯平滑滤波，滤除荧光图像中的孤立噪声点，平滑图像边缘，以方便后续的边界的提取和确定；
S1022：对平滑滤波后的荧光图像进行梯度边缘检测；
本步骤中，对滤波后的荧光图像进行梯度边缘的检测，主要是利用滤波后的荧光图像对应的灰度图像，计算出荧光图像的梯度特征，按照梯度边缘的强弱，从背景目标中分离出来待识别的前景目标区域；
S1023：依据梯度边缘检测的检测结果对所述荧光图像进行荧光目标前景提取；
本步骤中，对步骤S202的前景检测区域进行目标进行提取，利用边缘的跟踪，把待识别的前景目标区域按照连通区域的属性逐一标记出来，记录边缘的位置；
S1024：依据提取到的目标前景对所述荧光图像进行单个荧光区域的分割，得到每个检测目标的荧光图像区域；
步骤S1023处理完成以后，必然会发现存在大量的交叠的检测目标，因此还需要对交叠的检测目标进行分离，本步骤可通过利用距离变换等一系列算法，对边缘进行标记，最终对单个检测目标进行分离得到每个检测目标的荧光图像区域。
在本申请另一实施例公开的技术方案中除了采用荧光成像系统采集到目标样品的荧光图像之外，还需要采集得到其明场图像，因此，上述步骤中，所述采用显微荧光成像系统对目标样品进行荧光成像得到荧光图像，具体为：采用显微荧光成像系统采集得到目标样品的荧光图像和明场图像，此时，在对荧光图像进行高斯平滑滤波之前还可以包括：
步骤S1011：通过明场识别技术对明场图像中的明场目标进行识别，并且对识别出的明场目标中的团块区域进行分离，得到明场图像中每个检测目标的位置信息以及尺寸信息；
具体的，本步骤中，针对明场拍摄的明场图像，从背景目标中分离出前景，具体的方法是对明场图像进行灰度化、亮度均一化、图像增强、二值化操作，从背景目标中提取出前景的检测目标图像；上述提取到的前景的检测目标图像会在大量的结团等情况，使得多个检测目标图像重叠在一起，必须对检测目标图像进行分离，获得单个的检测目标图像，求取单个检测目标图像的边缘，并且记录边缘信息；
步骤S1012：将所述明场识别技术获取到的检测目标的位置信息以及尺寸信息与荧光图像中对应的区域进行对比，对荧光图像进行识别；判断所述荧光图像中对应明场图像检测到的各个检测目标的位置信息、尺寸信息的区域是否为空，如果否，执行步骤S1013，如果是，执行步骤S1014；
步骤S1013：依据采用得到的明场图像识别技术检测到的每个检测目标的位置、尺寸信息在所述荧光图像中相对应的区域，由所述荧光图像中确定对应的每个检测目标的荧光图像区域；
步骤S1014：如果所述述荧光图像中对应明场图像中检测目标的位置、尺寸信息的区域为空，对所述荧光图像进行高斯平滑滤波并执行后续操作；
所述对荧光图像进行识别，具体为：针对拍摄的荧光图像，从背景目标中分离出荧光目标所在的区域，具体的方法是图像进行灰度化、亮度均一化、边缘增强以及图像的二值化操作，然后进行图像边缘的提取，把荧光目标所在的区域都标记存储下来；
所述明场识别具体为：是以明场图像中各个检测目标分割的区域(单个检测目标的边缘，并且记录边缘信息)为标记，利用明场图像技术的分割结果对荧光图像标进行分割，提取得到单个检测目标在荧光图像中的区域，也就是该处存在明场目标，那么在荧光目标中也标记出这个位置，使用明场目标一一对应进行荧光目标分离。
上述方案中，采用显微荧光成像系统采集得到目标样品的荧光图像和明场图像，具体过程为：
采用显微荧光成像系统采集得到目标样品的荧光图像，判定得到的目标样品的荧光图像区域是否满足预设条件，如果否，采用显微荧光成像系统采集得到目标样品的明场图像。
如果该处荧光目标用图2中公开的方法明显分割不合理，比如边缘分割碎裂化，或者边缘弥散开来导致荧光的区域明显不正确，系统会自动采用图3中方法进行校正识别，也就是用明场目标的识别结果为依据来操作，对此在本实施例中，还公开了一种将上述两种实施例公开的智能识别处理方法相结合的处理方式，该方式，以图2中公开的处理方式为主，以图3公开的处理方式为辅，具体的，参见图4，
得到每个检测目标的荧光图像区域之后，还包括：
步骤S401：判定得到的检测目标的荧光图像区域是否满足预设条件，如果否执行步骤S402；
步骤S402：采用显微荧光成像系统采集得到目标样品的明场图像，执行步骤S1011。
与上述方法相对应，本申请还公开了一种荧光图像的荧光强度确定系统，所述方法和系统可以相互借鉴，参见图5，该系统可以包括：
显微荧光成像系统100，用于对目标样品进行荧光成像；；
智能识别系统200，用于对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，得到所述荧光图像中每个检测目标对应的荧光图像区域；
数据处理系统300，所述数据处理系统300包括：第一计算单元310、第二计算单元320、第三计算单元330、第四计算单元以及340以及流式聚类分析单元350；
所述第一计算单元310用于计算得到每个检测目标的荧光图像区域的累计灰度值，将所述累计灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的累计电压脉冲值表达的累计荧光强度值；
所述第二计算单元320用于计算得到每个检测目标的荧光图像区域的最大灰度值，将所述最大灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的最大电压脉冲值表达的最大荧光强度值；
所述第三计算单元330用于计算得到每个检测目标的荧光图像区域的平均灰度值，将所述平均灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的平均电压脉冲值表达的平均荧光强度值；
所述第四计算单元340用于计算每个检测目标的明场图像区域的直径值作为目标大小值，并确定该值作为开展流式聚类分析中正向散射通道以光电倍增管采集到的电压脉冲值表达的检测目标的大小值的对应替代指标
所述流式聚类分析单元350，用于用于基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析。
与上述方法相对应，本系统中所述第一计算单元310、第二计算单元320、第三计算单元330、第四计算单元340可以择一存在，即，本系统中可以选择性的有第一计算单元310、第二计算单元320、第三计算单元330和/或第四计算单元340。
与上述方法相对应，所述智能识别系统200被配置为：对所述荧光图像进行智能识别处理，提取单个检测目标荧光图像的边缘，边缘区域内的所有像素即为单个检测目标的荧光图像区域。
与上述方法相对应，所述智能识别系统200，包括至少荧光目标提取单元，所述荧光目标提取单元包括：
滤波单元，用于对荧光图像进行高斯平滑滤波；
梯度检测单元，用于对平滑滤波后的荧光图像进行梯度边缘检测；
前景提取单元，用于依据梯度边缘检测的检测结果对所述荧光图像进行荧光目标前景提取；
分割单元，用于依据提取到的目标前景对所述荧光图像进行单个荧光区域的分割，得到每个检测目标的荧光图像区域。
与上述方法相对应，所述显微荧光成像系统100还用于采集得到目标样品的明场图像；
所述智能识别系统200，还可以包括：明场细胞分离系统，所述明场细胞分离系统可以包括：
明场识别单元，用于对明场图像中的明场目标进行识别，并且对识别出的明场目标中的团块区域进行分离，得到明场图像中每个检测目标的位置、尺寸信息；
明场细胞分离单元，用于判断所述述荧光图像中对应明场图像中各个检测目标的位置、尺寸信息的区域是否为空，如果否，采用得到的明场图像中每个检测目标的位置、尺寸信息相对应的区域，由所述荧光图像中确定对应的每个检测目标的荧光图像区域，如果是，向荧光目标输出用于触发荧光目标动作的触发信号。
与上述方法相对应，在本申请另一实施例公开的技术方案中，所述显微荧光成像系统还用于采集得到目标样品的明场图像；
所述智能识别系统300，包括：荧光目标提取单元和明场细胞分离系统，本实施例虽然也包含荧光目标和明场细胞分离系统，但是两者的工作时序与上述实施例中的荧光目标和明场细胞分离系统的工作时序不同，具体的：
在本实施例中，所述荧光目标，还包括：判断单元，用于判定得到的细胞的荧光图像区域是否满足预设条件，如果否，控制显微荧光成像系统采集得到目标样品的明场图像，向所述明场细胞分离系统输出触发信号。
为了描述的方便，描述以上系统时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统或系统实施例而言，由于其基本相似于方法实施例，所以描述得比较简单，相关之处参见方法实施例的部分说明即可。以上所描述的系统及系统实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.一种荧光图像的荧光强度确定方法，其特征在于，包括：
采用显微荧光成像系统对目标样品进行荧光成像得到荧光图像；
对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，得到所述荧光图像中每个检测目标对应的荧光图像区域；
计算得到每个检测目标的荧光图像区域的累计灰度值，将所述累计灰度值作为流式聚类分析中对应荧光通道中以光电倍增管采集到的累计电压脉冲值表达的累计荧光强度值；
计算得到每个检测目标的荧光图像区域的最大灰度值，将所述最大灰度值作为流式聚类分析中对应荧光通道中以光电倍增管采集到的最大电压脉冲值表达的最大荧光强度值；
计算得到每个检测目标的荧光图像区域的平均灰度值，将所述平均灰度值作为流式聚类分析中对应荧光通道中以光电倍增管采集到的平均电压脉冲值表达的平均荧光强度值；
计算每个检测目标的明场图像区域的直径值作为目标大小值，并确定该值作为开展流式聚类分析中正向散射通道中以光电倍增管采集到的电压脉冲值表达的检测目标的大小值的对应替代指标
基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析。
2.根据权利要求1所述的荧光图像的荧光强度确定方法，其特征在于，基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析，包括：
依据预设映射表查找得到分别与所述累计灰度值、最大灰度值、平均灰度值和明场直径值一一对应的累计荧光强度值、最大荧光强度值、平均荧光强度值和目标大小值，基于查表得到的累计荧光强度值、最大荧光强度值、平均荧光强度值和目标大小值进行流式聚类分析。
3.根据权利要求1所述的荧光图像的荧光强度确定方法，其特征在于，对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，包括：
对荧光图像进行高斯平滑滤波；
对平滑滤波后的荧光图像进行梯度边缘检测；
依据梯度边缘检测的检测结果对所述荧光图像进行荧光目标前景提取；
依据提取到的目标前景对所述荧光图像进行单个荧光区域的分割，得到每个检测目标的荧光图像区域。
4.根据权利要求3所述的荧光图像的荧光强度确定方法，其特征在于，
所述采用显微荧光成像系统对目标样品进行荧光成像得到荧光图像，具体为：采用显微荧光成像系统采集得到目标样品的荧光图像和明场图像；所述对荧光图像进行高斯平滑滤波之前还包括：
对明场图像中的明场目标进行识别，并且对识别出的明场目标中的团块区域进行分离，得到明场图像中每个检测目标的位置、尺寸信息；
对荧光图像进行识别，判断所述述荧光图像中对应明场图像中各个检测目标的位置、尺寸信息的区域是否为空，如果否，采用得到的明场图像中每个检测目标的位置、尺寸信息相对应的区域，由所述荧光图像中确定对应的每个检测目标的荧光图像区域；
如果所述述荧光图像中对应明场图像中检测目标的位置、尺寸信息的区域为空，对荧光图像进行高斯平滑滤波并执行后续操作。
5.根据权利要求4所述的荧光图像的荧光强度确定方法，其特征在于，采用显微荧光成像系统采集得到目标样品的荧光图像和明场图像，具体为：
采用显微荧光成像系统采集得到目标样品的荧光图像，判定得到的目标样品的荧光图像区域是否满足预设条件，如果否，采用显微荧光成像系统采集得到目标样品的明场图像。
6.一种荧光图像的荧光强度确定系统，其特征在于，包括：
显微荧光成像系统，用于对目标样品进行荧光成像；
智能识别系统，用于对荧光图像中每个检测目标的荧光图像区域进行边缘提取和分割，得到所述荧光图像中每个检测目标对应的荧光图像区域；
数据处理系统，所述数据处理系统包括：第一计算单元、第二计算单元、第三计算单元、第四计算单元和流式聚类分析单元；
所述第一计算单元用于计算得到每个检测目标的荧光图像区域的累计灰度值，将所述累计灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的累计电压脉冲值表达的累计荧光强度值；
所述第二计算单元用于计算得到每个检测目标的荧光图像区域的最大灰度值，将所述最大灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的最大电压脉冲值表达的最大荧光强度值；
所述第三计算单元用于计算得到每个检测目标的荧光图像区域的平均灰度值，将所述平均灰度值作为流式聚类分析中对应荧光通道以光电倍增管采集到的平均电压脉冲值表达的平均荧光强度值；
所述第四计算单元用于计算每个检测目标的明场图像区域的直径值作为目标大小值，并确定该值作为开展流式聚类分析中正向散射通道以光电倍增管采集到的电压脉冲值表达的检测目标的大小值的对应替代指标；
所述流式聚类分析单元，用于基于所述累计灰度值、最大灰度值、平均灰度值和明场直径值进行流式聚类分析。
7.根据权利要求6所述的荧光图像的荧光强度确定系统，其特征在于，所述流式聚类分析单元具体用于：依据预设映射表查找得到分别与所述累计灰度值、最大灰度值、平均灰度值和明场大小值一一对应的累计荧光强度值、最大荧光强度值、平均荧光强度值和目标大小值，基于查表得到的累计荧光强度值、最大荧光强度值、平均荧光强度值和目标大小值进行流式聚类分析。
8.根据权利要求6所述的荧光图像的荧光强度确定系统，其特征在于，智能识别系统，包括至少荧光目标提取单元，所述荧光目标提取单元包括：
滤波单元，用于对荧光图像进行高斯平滑滤波；
梯度检测单元，用于对平滑滤波后的荧光图像进行梯度边缘检测；
前景提取单元，用于依据梯度边缘检测的检测结果对所述荧光图像进行荧光目标前景提取；
分割单元，用于依据提取到的目标前景对所述荧光图像进行单个荧光区域的分割，得到每个检测目标的荧光图像区域。
9.根据权利要求8所述的荧光图像的荧光强度确定系统，其特征在于，所述显微荧光成像系统还用于采集得到目标样品的明场图像；
所述智能识别系统，还包括：明场细胞分离系统，所述明场细胞分离系统包括：
明场识别单元，用于对明场图像中的明场目标进行识别，并且对识别出的明场目标中的团块区域进行分离，得到明场图像中每个检测目标的位置、尺寸信息；
明场细胞分离单元，用于判断所述述荧光图像中对应明场图像中各个检测目标的位置、尺寸信息的区域是否为空，如果否，采用得到的明场图像中每个检测目标的位置、尺寸信息相对应的区域，由所述荧光图像中确定对应的每个检测目标的荧光图像区域，如果是，向荧光目标提取单元输出用于触发荧光目标动作的触发信号。
10.根据权利要求9所述的荧光图像的荧光强度确定系统，其特征在于，
所述荧光目标提取单元，还包括：判断单元，用于判定得到的细胞的荧光图像区域是否满足预设条件，如果否，控制显微荧光成像系统采集得到目标样品的明场图像，向所述明场细胞分离系统输出触发信号。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710596449.2A CN107356515B (zh) | 2017-07-20 | 2017-07-20 | 一种荧光图像的荧光强度确定方法和系统 |
| EP18835191.0A EP3657154A4 (en) | 2017-07-20 | 2018-07-11 | METHOD AND SYSTEM FOR DETERMINING THE FLUORESCENT INTENSITY OF A FLUORESCENT IMAGE |
| JP2020523477A JP6851670B2 (ja) | 2017-07-20 | 2018-07-11 | 蛍光画像の蛍光強度確定方法及びシステム |
| US16/629,400 US11249005B2 (en) | 2017-07-20 | 2018-07-11 | Method and system for determining fluorescence intensity of fluorescence image |
| PCT/CN2018/095239 WO2019015514A1 (zh) | 2017-07-20 | 2018-07-11 | 一种荧光图像的荧光强度确定方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710596449.2A CN107356515B (zh) | 2017-07-20 | 2017-07-20 | 一种荧光图像的荧光强度确定方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107356515A CN107356515A (zh) | 2017-11-17 |
| CN107356515B true CN107356515B (zh) | 2020-04-03 |
Family
ID=60284537
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710596449.2A Active CN107356515B (zh) | 2017-07-20 | 2017-07-20 | 一种荧光图像的荧光强度确定方法和系统 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US11249005B2 (zh) |
| EP (1) | EP3657154A4 (zh) |
| JP (1) | JP6851670B2 (zh) |
| CN (1) | CN107356515B (zh) |
| WO (1) | WO2019015514A1 (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107356515B (zh) * | 2017-07-20 | 2020-04-03 | 上海睿钰生物科技有限公司 | 一种荧光图像的荧光强度确定方法和系统 |
| CN109800631B (zh) * | 2018-12-07 | 2023-10-24 | 天津大学 | 基于掩膜区域卷积神经网络的荧光编码微球图像检测方法 |
| CN109815870B (zh) * | 2019-01-17 | 2021-02-05 | 华中科技大学 | 细胞表型图像定量分析的高通量功能基因筛选方法及系统 |
| CN110009619A (zh) * | 2019-04-02 | 2019-07-12 | 清华大学深圳研究生院 | 一种基于荧光编码的液相生物芯片的图像分析方法 |
| CN111289484A (zh) * | 2020-03-11 | 2020-06-16 | 哈尔滨工业大学（威海） | 一种基于罗丹明b荧光特性的冷表皮检测方法 |
| CN113237818A (zh) * | 2021-05-28 | 2021-08-10 | 上海睿钰生物科技有限公司 | 一种细胞分析方法和系统 |
| JP2024506817A (ja) * | 2021-01-25 | 2024-02-15 | ベクトン・ディキンソン・アンド・カンパニー | 空間的に分離されたレーザにわたる散乱信号を使用して液滴遅延を判定するための方法及びシステム |
| CN112819795B (zh) * | 2021-02-05 | 2023-03-31 | 复旦大学附属中山医院 | 基于多色荧光图片的细胞数量和面积的统计方法 |
| CN112907582B (zh) * | 2021-03-24 | 2023-09-29 | 中国矿业大学 | 面向矿井的图像显著性提取去雾方法、装置及人脸检测 |
| CN113450279B (zh) * | 2021-07-01 | 2023-03-14 | 维柯基科技(上海)有限公司 | 一种多孔荧光微阵列图像的荧光强度检测方法、装置、计算机设备及计算机可读存储介质 |
| CN116519650B (zh) * | 2023-04-17 | 2024-03-12 | 北京大学口腔医学院 | 一种骨组织结构状态检测方法 |
| CN116124754B (zh) * | 2023-04-17 | 2023-07-25 | 山东省中源联科生物工程集团有限公司 | 通过图像分析检测透明质酸钠凝胶中交联剂残留量的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104266956A (zh) * | 2014-09-02 | 2015-01-07 | 上海凯度机电科技有限公司 | 高内涵图像流式生物显微分析方法 |
| CN105092550A (zh) * | 2015-07-01 | 2015-11-25 | 上海睿钰生物科技有限公司 | 一种荧光显微成像方法和装置 |
| CN105928921A (zh) * | 2016-06-07 | 2016-09-07 | 河南科技大学 | 一种通过荧光染色图片数字化进行细胞pH值计算的方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3411115B2 (ja) * | 1994-12-26 | 2003-05-26 | シスメックス株式会社 | フローサイトメータ |
| DE10218706B4 (de) * | 2002-04-26 | 2008-02-28 | Leica Microsystems Cms Gmbh | Verfahren zur zeitoptimierten Erfassung von speziellen Spektren mit einem Scanmikroskop |
| JP2005315862A (ja) * | 2004-03-30 | 2005-11-10 | Sysmex Corp | 子宮頸癌のスクリーニング方法及び子宮頸癌診断薬 |
| US7264794B2 (en) * | 2004-12-17 | 2007-09-04 | The General Hospital | Methods of in vivo cytometry |
| JP4982385B2 (ja) * | 2005-02-01 | 2012-07-25 | アムニス コーポレイション | イメージングフローサイトメータを使用した血液及び細胞の分析 |
| JPWO2011099565A1 (ja) * | 2010-02-10 | 2013-06-17 | シスメックス株式会社 | 血液中癌細胞の検出方法及びそれに用いるプログラム |
| US9250229B2 (en) * | 2011-09-25 | 2016-02-02 | Theranos, Inc. | Systems and methods for multi-analysis |
| US11249799B2 (en) * | 2011-09-26 | 2022-02-15 | Labrador Diagnostics Llc | Methods, systems, and devices for real time execution and optimization of concurrent test protocols on a single device |
| US20150024385A1 (en) * | 2013-07-22 | 2015-01-22 | Wisconsin Alumni Research Foundation | Prediction of fertility in males |
| EP3054279A1 (en) | 2015-02-06 | 2016-08-10 | St. Anna Kinderkrebsforschung e.V. | Methods for classification and visualization of cellular populations on a single cell level based on microscopy images |
| CN107356515B (zh) | 2017-07-20 | 2020-04-03 | 上海睿钰生物科技有限公司 | 一种荧光图像的荧光强度确定方法和系统 |
-
2017
- 2017-07-20 CN CN201710596449.2A patent/CN107356515B/zh active Active
-
2018
- 2018-07-11 WO PCT/CN2018/095239 patent/WO2019015514A1/zh unknown
- 2018-07-11 EP EP18835191.0A patent/EP3657154A4/en active Pending
- 2018-07-11 US US16/629,400 patent/US11249005B2/en active Active
- 2018-07-11 JP JP2020523477A patent/JP6851670B2/ja active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104266956A (zh) * | 2014-09-02 | 2015-01-07 | 上海凯度机电科技有限公司 | 高内涵图像流式生物显微分析方法 |
| CN105092550A (zh) * | 2015-07-01 | 2015-11-25 | 上海睿钰生物科技有限公司 | 一种荧光显微成像方法和装置 |
| CN105928921A (zh) * | 2016-06-07 | 2016-09-07 | 河南科技大学 | 一种通过荧光染色图片数字化进行细胞pH值计算的方法 |
Non-Patent Citations (1)
| Title |
|---|
| RNA干扰Apollon基因逆转白血病K562细胞多药耐药;贾秀红;《中国癌症杂志》;20131231;第713-720 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US11249005B2 (en) | 2022-02-15 |
| JP2020527274A (ja) | 2020-09-03 |
| WO2019015514A1 (zh) | 2019-01-24 |
| CN107356515A (zh) | 2017-11-17 |
| JP6851670B2 (ja) | 2021-03-31 |
| US20200141856A1 (en) | 2020-05-07 |
| EP3657154A1 (en) | 2020-05-27 |
| EP3657154A4 (en) | 2021-04-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107356515B (zh) | 一种荧光图像的荧光强度确定方法和系统 | |
| KR102469620B1 (ko) | 생물학적 입자의 분류 시스템 및 방법 | |
| US10445619B2 (en) | Automatic glandular and tubule detection in histological grading of breast cancer | |
| US8379961B2 (en) | Mitotic figure detector and counter system and method for detecting and counting mitotic figures | |
| JP5376024B1 (ja) | 画像処理装置、プログラム及び画像処理システム | |
| EP2377073B1 (en) | Multi-nucleated cell classification and micronuclei scoring | |
| Chang et al. | Nuclear segmentation in H&E sections via multi-reference graph cut (MRGC) | |
| US20230177341A1 (en) | Using a First Stain to Train a Model to Predict the Region Stained by a Second Stain | |
| JP5413501B1 (ja) | 画像処理装置、画像処理システム及びプログラム | |
| US9443129B2 (en) | Methods and apparatus for image analysis using threshold compactness features | |
| JP2016507229A5 (zh) | ||
| KR102393268B1 (ko) | 도로 모양 분석을 통한 교통사고 발생 확률 예측 장치 및 방법 | |
| Mabaso et al. | Spot detection methods in fluorescence microscopy imaging: a review | |
| Habibzadeh et al. | Application of pattern recognition techniques for the analysis of thin blood smear images | |
| CN114638832A (zh) | 一种基于分水岭算法的dpcr液滴荧光检测方法 | |
| CN107367456B (zh) | 一种免洗涤图像类流式荧光检测方法和系统 | |
| JP2020518045A (ja) | 医療画像検出 | |
| EP2776974B1 (en) | Methods and apparatus for image analysis using threshold compactness features | |
| Jeffrey et al. | High definition licence plate detection algorithm | |
| CN113837255B (zh) | 预测以细胞为基质的抗体核型类别的方法、设备和介质 | |
| TWI637146B (zh) | Cell classification method | |
| Restif | Towards safer, faster prenatal genetic tests: Novel unsupervised, automatic and robust methods of segmentation of nuclei and probes | |
| CN117808836A (zh) | 基于直线边缘检测的谷糙分界线检测方法、设备及介质 | |
| KR20140147454A (ko) | 차량패턴 히스토그램 생성방법 및 이를 이용한 차량패턴 검출방법 | |
| CN117711010A (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 | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A method and system for determining fluorescence intensity of fluorescence image Effective date of registration: 20210521 Granted publication date: 20200403 Pledgee: The Bank of Shanghai branch Caohejing Limited by Share Ltd. Pledgor: SHANGHAI RUIYU BIOTECH Co.,Ltd. Registration number: Y2021310000038 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20231103 Granted publication date: 20200403 Pledgee: The Bank of Shanghai branch Caohejing Limited by Share Ltd. Pledgor: SHANGHAI RUIYU BIOTECH Co.,Ltd. Registration number: Y2021310000038 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right |