一种基于显微镜图像的微囊藻类细胞统计方法与流程

1.本发明涉及水生态环境监测技术领域，尤其涉及一种基于显微镜图像的微囊藻类细胞统计方法。


背景技术：

2.藻类细胞统计是先利用显微镜和高清工业相机来采集藻类图像，然后通过深度学习检测模型识别出微囊藻类和它的像素坐标，需要设计一种图像模式识别方法来统计图像中微囊藻类的细胞个数。
3.微囊藻类细胞经常会呈现比较密集的分布或者是相互叠加在一起的现象，现有的方法直接采用深度学习模型来统计微囊藻类的细胞个数，会产生较严重的漏检问题，因此，本发明提出一种基于显微镜图像的微囊藻类细胞统计方法，以解决现有技术中的不足之处。


技术实现要素：

4.针对上述问题，本发明提出了一种基于显微镜图像的微囊藻类细胞统计方法，该方法利用微囊藻类细胞区域均值和方差特征，有效地计算出了藻类细胞密度分布图，并对密度分布较稀疏的图像采用斑点检测的方法修正了微囊藻类细胞的像素尺寸，不仅降低了细胞间重叠引起的统计误差，也提高了细胞个数统计的精度；另外，本方法不需要在深度学习模型中直接检测微囊藻类的细胞，减少了数据标注的工作，提高了模型训练和模型优化的效率。
5.本发明提出一种基于显微镜图像的微囊藻类细胞统计方法，包括以下步骤：
6.步骤一、获取彩色图像并利用高斯模糊化对图像进行去噪处理，然后将rgb图像转换成hsv图像，并在hsv图像的v通道上进行对比度增强处理，再转换为rgb图像，获得对比度增强后的rgb图像；
7.步骤二、在对比度增强后的v通道图像上进行otsu(大律)二值化处理，并对处理后的图像进行形态学闭(close)的操作，获得处理后的二值图像；
8.步骤三、在对比度增强后的rgb图像上，分析二值图像白色区域(像素值为255)的均值与方差，并通过计算得出微囊藻类细胞密度分布图；
9.步骤四、计算图像中藻类细胞占比和分析藻类细胞区域内部密度情况，由此判断微囊藻类细胞分布的稀疏程度，并对分布稀疏的微囊藻类细胞像素尺寸进行修正；
10.步骤五、根据微囊藻类细胞密度分布图计算密度加权系数，再对统计的微囊细胞个数进行修正，获得最终微囊藻类细胞统计数据。
11.进一步改进在于：所述步骤一中获取图像去噪处理后，将rgb图像转化为hsv图像，并分离出h、s、v通道得到himg、simg和vimg三个单通道图像，再对vimg图像采用对比度受限的自适应直方图均衡化算法进行对比度拉伸操作，得到增强后的图像vimg1，然后将himg、simg和vimg1三个单通道的图像合并成新的三通道hsv图像，最后将新的三通道hsv图像转
化为对比度增强后的rgb图像。
12.进一步改进在于：所述步骤二中在图像vimg1上利用otsu方法进行二值化分割，得到二值图像并进行形态学close操作，得到处理后的二值图像，其中二值图像中白色为微囊藻类细胞区域。
13.进一步改进在于：所述步骤三具体包括
14.s1、在对比度增强后的rgb图像上，计算二值图像中像素值为255对应的像素点的r、g和b三个通道的均值和平均方差；
15.s2、设置比例阈值，分别将它乘以三个通道的平均方差，得到三个通道的偏离阈值；
16.s3、调查二值图像中像素值为255的坐标点，得到该坐标点对应彩色图像的像素值r、g和b，再分别计算三个通道当前像素值与均值的偏差；
17.s4、生成新的二值图像，大小与原二值图像保持一致，其中默认初始值为0，得到的图像为细胞密度分布图，当计算的偏差绝对值大于对应的偏离阈值，则将细胞密度分布图上对应的像素赋值为255，即白色。
18.进一步改进在于：所述步骤四中修正细胞尺寸具体包括
19.e1、统计二值图像中白色区域的像素面积，将其除以图像的像素面积，获得图像中微囊藻类细胞的占比数据；
20.e2、采用自适应阈值分割法计算藻类细胞区域内部的细胞占比数据；
21.e3、将图像中微囊藻类细胞的占比数据和藻类细胞区域内部的细胞占比数据与阈值进行对比判断，并根据判断结果对微囊藻类细胞尺寸进行修正。
22.进一步改进在于：所述e3中细胞尺寸修正判断时当微囊藻类细胞的占比数据或藻类细胞区域内部的细胞占比数据中的任意一个数据值小于阈值时，就会对微囊藻类细胞尺寸进行修正，反之不进行修正。
23.进一步改进在于：所述微囊藻类细胞尺寸进行修正时，先采用斑点检测算法在灰度图像上检测出微囊藻类细胞，然后找出细胞的中位值像素面积，通过中位值像素面积验证所有检测到的微囊藻类细胞，将像素面积偏差大的数据删除，最后计算剩余微囊藻类细胞的平均像素面积尺寸，以平均像素面积尺寸来修正微囊藻类单位细胞像素尺寸。
24.进一步改进在于：所述步骤五中具体包括
25.g1、计算藻类细胞密度分布图上所有像素的密度平均值，并根据它的数值计算出密度加权系数；
26.g2、统计藻类细胞密度分布图中白色区域的像素面积，然后计算该面积与藻类细胞区域面积的比例；
27.g3、计算出细胞未叠加状态下微囊藻类细胞的个数；
28.g4、最后根据藻类细胞密度分布情况对统计的细胞个数进行加权修正，得到微囊藻类细胞个数。
29.本发明的有益效果为：本发明利用微囊藻类细胞区域均值和方差特征，有效地计算出藻类细胞密度分布图，并对密度分布较稀疏的图像采用斑点检测的方法修正了微囊藻类细胞的像素尺寸，不仅降低了细胞间重叠引起的统计误差，也提高了细胞个数统计的精度；
30.另外，本方法不需要在深度学习模型中直接检测微囊藻类的细胞，减少了数据标注的工作，提高了模型训练和模型优化的效率；
31.且本发明方法为独立的模块，便于算法后期的优化和维护。
附图说明
32.图1为本发明实施例流程图。
33.图2为本发明实施例图像预处理流程图。
34.图3为本发明实施例图像二值化流程图。
35.图4为本发明实施例计算细胞密度分布图流程图。
36.图5为本发明实施例细胞密度分布效果图。
37.图6为本发明实施例修正细胞尺寸流程图。
38.图7为本发明实施例计算细胞区域内部细胞占比流程图。
39.图8为本发明实施例删除误检测微囊藻类细胞流程图。
40.图9为本发明实施例统计微囊藻类的细胞个数流程图。
41.图10为本发明实施例按细胞密度对统计结果进行修正流程图。
42.图11为本发明实施例微囊藻类细胞个数统计效果图。
具体实施方式
43.为了使发明实现的技术手段、达成目的与功效易于明白了解，下面结合具体实施方式，进一步阐述本发明。
44.根据图1-11所示，本实施例提出一种基于显微镜图像的微囊藻类细胞统计方法，包括以下步骤：
45.步骤一、图像预处理：增强图像的对比度；
46.步骤二、图像二值化：在增强后的图像上进行二值化处理；
47.步骤三、计算细胞密度分布图：计算图像上微囊藻类细胞密度分布图；
48.步骤四、修正细胞尺寸：修正图像中微囊藻类的单位细胞像素尺寸；
49.步骤五、统计微囊藻类细胞个数：根据微囊藻类的像素面积和密度分布图，统计出它的细胞个数。
50.根据说明书附图2所示，所述图像预处理模块包括以下内容:
51.(1)首先对彩色图像进行缩放处理，提高整体运行效率；
52.(2)对缩放后的彩色图像进行高斯模糊化去噪；
53.(3)将rgb图像转换为hsv图像；
54.(4)从hsv三通道图像中分离出h、s、v通道形成三个单通道图像himg、simg和vimg；
55.(5)对vimg图像采用clahe(对比度受限的自适应直方图均衡化)算法进行对比度拉伸操作，得到增强后的图像vimg1；
56.(6)将himg、simg和vimg1三个单通道的图像合并成新的三通道图像hsv1；
57.(7)将hsv1图像转换为rgb图像rgbimg。
58.根据说明书附图3所示，所述图像二值化包括以下内容:
59.(1)在图像vimg1上进行otsu(大律)方法进行二值化分割，得到二值图像
binaryimg,在二值图像中微囊藻类细胞为白色，对应的像素值为255，其它地方为黑色，像素值为0；
60.(2)对二值图像binaryimg进行形态学闭(close)操作。
61.根据说明书附图4和说明书附图5所示，所述计算细胞密度分布图包括以下内容:
62.(1)在三通道彩色图像rgbimg上分别计算r、g和b三个通道的均值和平均方差，它们分别表示为meanr、meang、meanb和stdr、stdg、stdb；
63.(2)在rgbimg上计算均值和方差时，只统计二值图像binaryimg中白色(像素值255)对应的像素；
64.(3)设置一个比例阈值(0.7)，分别将它乘以stdr、stdg和stdb，得到三个通道的最小偏离阈值，将它们表示为thresr、thresg和thresb；
65.(4)新建一个二值图像binaryimg1，图像大小与binaryimg一致，图像初始值为0。找到binaryimg中像素值为255的坐标点，得到该坐标点对应的rgbimg的像素值r、g和b。分别计算三个通道当前像素值与均值的偏差(r-meanr、g-meang和b-meanb)，如果三个通道偏差的绝对值都大于对应的偏离阈值(abs(r-meanr)》thresr&&abs(g-meang)》thresg&&abs(b-meanb)》thresb)，那么将binaryimg1图像上该像素赋值为255，标记为细胞叠加点；
66.(5)对二值图像binaryimg1进行形态学开(open)操作,降低孤立点对密度分布计算造成的干扰；
67.(6)通过观察微囊藻类图像发现：当微囊藻类细胞叠加在一起时，图像上该区域的颜色会变深，密度越大，颜色越深。因此，根据这个特性，先通过二值图像找到微囊藻类细胞区域(二值图像上白色区域)，然后计算细胞区域内部的均值和方差。如果细胞未叠加在一起，那么它的方差会很小，计算出的细胞密度分布图上白色覆盖的区域也会很少；如果细胞叠加得较多，那么它的方差会较大，计算出的细胞密度分布图上白色覆盖的区域也会较多。
68.根据说明书附图6所示，所述修正细胞尺寸包括以下内容:
69.(1)统计二值图像binaryimg上微囊藻类细胞(白色区域)的像素面积areat，用该面积除以图像的像素面积，计算出图像中藻类细胞占比；
70.(2)如果图像中藻类细胞占比较小，则说明图像上微囊藻类细胞整体分布较稀疏，可以对藻类细胞尺寸进行修正；
71.(3)如果计算的图像中藻类细胞占比大于设定的阈值，将进一步分析藻类细胞区域内部的细胞占比；
72.(4)如果藻类细胞区域内部的细胞占比小于设定的阈值，则说明藻类细胞区域内部的细胞分布较稀疏，那么，也可以对藻类细胞尺寸进行修正；
73.(5)如果图像中藻类细胞占比和藻类细胞区域内部的细胞占比都大于设定的阈值，那么将不对微囊藻类细胞尺寸进行修正；
74.(6)在vimg1图像上，通过斑点(simpleblobdetector)检测算法来检测藻类细胞；
75.(7)删除斑点检测中的误检测数据；
76.(8)如果图像上检测到多个大小相似的藻类细胞，则计算它们的平均像素尺寸，并对微囊藻类单位细胞尺寸进行修正。
77.根据说明书附图7所示，所述计算细胞区域内部细胞占比包括以下内容:
78.(1)在图像vimg1上采用adaptivethreshold(自适应阈值分割)方法进行二值化分
割，得到二值图像binaryimg2；
79.(2)对binaryimg和binaryimg2进行and(与)操作，更新二值图像binaryimg2；
80.(3)统计二值图像binaryimg2的像素面积，并将该面积除以areat，计算出微囊藻类细胞区域内部的细胞占比。
81.根据说明书附图8所示，所述删除误检测微囊藻类细胞包括以下内容:
82.(1)分析检测的藻类细胞外接矩形框，合并重叠的矩形框；
83.(2)如果检测到矩形框内部存在另一个矩形框，则将内部的矩形框删除掉；
84.(3)按矩形框的像素面积进行排序，取中间值作为细胞像素面积参考值；
85.(4)遍历所有的矩形框，计算它们的像素面积和参考像素面积之间的偏差，删除偏差较大的矩形框。
86.根据说明书附图9所示，所述统计微囊藻类细胞个数包括以下内容:
87.(1)根据二值图像binaryimg，计算微囊藻类细胞的像素总面积areat；
88.(2)如果当前图像修正了藻类细胞尺寸，那么直接将areat除以单位细胞尺寸，统计出平面上微囊藻类的细胞个数count；
89.(3)对(2)统计的细胞个数，按藻类细胞密度分布图进行加权，统计出最终结果；
90.(4)如果当前图像未修正藻类细胞尺寸，那么直接以图像的分辨率作为像素总面积来除以默认的单位细胞尺寸，计算出平面上微囊藻类细胞个数的最大值；
91.(5)结合图像中藻类细胞像素面积占比来对(4)统计的结果进行修正；
92.(6)将(5)统计的平面上藻类细胞个数，再按照藻类细胞密度分布图进行加权，计算出最终的微囊藻类细胞个数。
93.根据说明书附图10所示，所述按细胞密度对统计结果进行修正包括以下内容:
94.(1)根据微囊藻类细胞密度分布图binaryimg1，统计每个白色像素点周围9*9邻域范围内白色像素点的个数num，然后将num除以邻域(9*9)面积，计算每个像素点的密度值；
95.(2)如果(1)计算的密度值低于设定的阈值，那么将binaryimg1图像上该像素点的值设置为0。通过这个方法，可以对藻类细胞密度分布图binaryimg1进行修正，提高平均密度加权系数计算的精度；
96.(3)统计微囊藻类细胞密度分布图binaryimg1中白色区域的像素面积aread，然后将aread除以areat，计算藻类细胞叠加区域相对于藻类细胞区域的面积占比rd；
97.(4)将binaryimg1中白色区域对应的密度值进行累加，然后除以aread，得到平均密度值，根据该值给定一个细胞密度加权系数r1；
98.(5)将前面计算的平面上微囊藻类细胞个数乘以藻类细胞叠加区域相对于藻类细胞区域的面积占比(count*rd)计算出叠加的细胞个数count1；
99.(6)将叠加的细胞个数乘以密度加权系数(count1*r1)，计算出最终叠加的细胞个数count2；
100.(7)将平面上统计的微囊藻类细胞的个数和最终叠加的微囊藻类细胞个数相加(count+count2)，计算出图像上微囊藻类细胞的个数。
101.根据说明书附图11所示，微囊藻类图像左上角深色数字为本发明统计的细胞个数，下方亮色的数字为人工统计的微囊藻类细胞个数
102.以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该
了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。