一种活性污泥相差显微图像絮体和丝状菌分割方法

1.本发明涉及污水处理智能检测方法，尤其涉及一种活性污泥相差显微图像絮体和丝状菌分割方法。


背景技术：

2.活性污泥显微图像由絮体和丝状菌组成，其正确的分割是提取活性污泥形态学参数的关键，但由于疏密度不同的絮体附着在丝状菌骨架上，存在非均质、结构多样等复杂特点，在相差显微晕圈和遮光效应影响下，难以区分絮体、丝状菌和背景。另外，不同批次污泥样品采集时光源强度变化、絮体和丝状菌分布不均衡且丝状菌目标小的问题，也增加分割难度。
3.活性污泥显微图像分割主要包括基于数字图像技术的分割和基于深度学习的分割。其中数字图像方法包括基于阈值的分割方法、基于边缘的分割方法以及基于区域的分割方法，这些方法均利用单一的方法对目标进行处理，容易出现欠分割、过分割以及无法分割的情况。深度学习的方法包括基于deeplab的分割方法、基于u-net的分割方法、基于pspnet分割方法，其分割方法发展速度很快，产生了各式各样的模型变体，但采集数据后需大量人力标注样本，所耗时间较多，成本较大，不适用于实际工程。


技术实现要素：

4.本发明公开了一种活性污泥相差显微图像絮体和丝状菌分割方法，提出根据絮体和丝状菌的特点采用不同的方法进行分割，无需标定样本，在保证模型精度的前提下，用于数字图像处理时容易出现的欠分割、过分割以及不能分割的情况，实现目标准确快速的分割，并提取絮体和丝状菌的形态学特征，满足线下部署以及应用在工程实际上的需求。其显微图像分割和定量信息提取技术在污泥膨胀的预报、絮状物和丝状菌的个体浓度测量生物传感器方面具有潜在价值。
5.本发明为实现上述目标，采用的技术方案是：一种活性污泥相差显微图像絮体和丝状菌分割方法，其特征在于，包括下述步骤：s1，图像采集：在污水处理厂的曝气池出口采集活性污泥样品，采用光学显微镜、工业数码相机和图像采集系统捕获活性污泥相差显微图像；s2，絮体分割：采用otsu阈值法对经过线性变换、图像反转和开操作后的图像进行阈值化检测，得到絮体的大致形状，并用边缘检测的结果确定最终的絮体；s3，丝状菌分割：采用多尺度retinex(msr)对图像进行增强后，通过顶帽底帽滤波、边缘检测、约束模糊聚类、回旋半径特征方法实现活性污泥丝状菌分割；s4，特征提取：从分割好的絮体和丝状菌图像中，分别提取出絮体和丝状菌的尺寸和形状等形态学特征参数。
6.进一步可选地，步骤s1包括如下步骤：s11，从污水处理厂曝气池出口取得的活性污泥样品；s12，采用光学显微镜、工业
数码相机和图像采集系统，从左上到右下以z字路径扫描载玻片上的污泥样品，污泥样品拍摄张rgb彩色图像，创建活性污泥显微图像数据集。
7.进一步可选地，步骤s2包括如下步骤：s21，利用灰度变换、图像反转以及开操作对图像进行预处理；s22，选择otsu阈值法对经过预处理的图像进行阈值化检测，得到絮体的大致结果；s23，通过形态学闭操作对分割出来的絮体去除断边，并适当填充孔洞部分，随后利用边缘检测的结果确定最终的絮体边缘。
8.进一步可选地，步骤s3包括如下步骤：s31，图像预处理，包括多尺度retinex、滤波和边缘检测；s32，利用约束模糊聚类对预处理后的图像进行分割；s33，对分割后的图像进行后处理操作，利用减小回转半径去碎片，同时判定老化絮体等，以便更准确的将丝状菌分割出来。
9.进一步可选地，减少回转半径公式如下：其中，和为二阶矩阵，定义如下：为二阶矩阵，定义如下：，为一阶矩阵，即图像质心的坐标：。
10.进一步可选地，步骤s4包括提取污泥絮体和丝状菌形态学特征。对图像中絮体和丝状菌的每个目标依次提取特征，以下为每个目标提取的形态学特征参数：
本发明有益效果：1.采用无监督模式分割活性污泥显微图像絮体和丝状菌，无需标定样本，节省人力；2.组合多尺度retinex频率域处理和顶帽底帽滤波、形态学运算空间域处理，抑制相差显微图像固有的伪影现象，消除背景影响，增强絮体和丝状菌目标；3.基于约束模糊聚类和形态学特征进一步提升微小目标丝状菌分割准确性和鲁棒性。
附图说明
11.附图作为本发明的一部分，用来提供对本发明的进一步的理解，本发明的示意性实施例及其说明用于解释本发明，但不构成对本发明的不当限定。显然，下面描述中的附图仅仅是一些实施例，对于本领域普通技术人员来说，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。在附图中：图1为基于活性污泥相差显微图像分割方法的流程图；图2为活性污泥相差显微图像采集流程图；图3为本发明在污泥显微图像测试集上的分割结果示例图；图4本方法在污泥显微图像数据集上的性能。
具体实施方式
12.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。如附图1所示，为本发明的一种活性污泥相差显微图像絮体和丝状菌分割方法流程图。本发明是一种活性污泥相差显微图像絮体和丝状菌分割方法，包括下述步骤：s1，图像采集。其中，s1包括如下步骤：s11，活性污泥相差显微图像的采集流程如图2所示，污泥样本取自污水处理厂曝气池出口，每次污泥样品取500 ml，样品采集、图像采集之间的间隔时间不超过3小时；s12，使用一个带有截面尖端枪头的校准移液枪，取10μl样品于载玻片上，并用24mm
×
24mm的正方形盖玻片覆盖，采用光学显微镜(nikon eclipse ts100)、工业数码相机(touptek toupcamucoms03100kpa)和图像采集软件(toupview)，从左上到右下以z字路径扫描载玻片上的污泥样本，污泥样品拍摄张rgb彩色图像，创建活性污泥显微图像数据集。显微镜相差放大倍数为100x，以确定絮凝物以及丝状菌的大小、形状和结构参数。由于相差显微镜对相位变化过于敏感，并且在成像过程中存在一些缺陷，使得所呈现的图像中产生一些伪影，如光晕，光斑和阴影。
13.s2，絮体分割。其中，s2包括如下步骤：s21，对灰度化的图像进行灰度变换，调整图像灰度值。原图像灰度范围为0-255，经过反复试验，设置调整范围为[0.3,0.8]，将小于255
×
0.3的灰度值设置为0，将大于255
×
0.8的灰度值设置为255，超出此区间的灰度保持不变。将经过灰度变换的图像再分别进行图像反转和开操作得到图像，利用灰度图像减去经过处理的图像，能极大地抑制光晕和阴影现象，同时使絮体对比度增强；s22，选择otsu阈值法对经过预处理的图像进行阈值化检测，得到絮体的大致结果；s23，由于原始灰度图像中灰度分布不均匀以及亮度明暗等原因，经过otsu算法分割出来的絮体效果并不完美，絮体内部有空隙，同时絮体和絮体的连接处出现断点等现象，因此通过形态学闭操作对分割出来的絮体去除断边，并适当填充孔洞部分，随后利用边缘检测的结果确定最终的絮体边缘。
[0014]
s3，丝状菌分割。其中，s3包括如下步骤：s31，首先，将rgb图像(红、绿、蓝)经过多尺度retinex进行图像增强，经过多尺度retinex增强，克服了不同样本因光源明暗所带来的纹理结构影响，提高了整个数据集的对
比度。对灰度图像进行顶帽和底帽滤波，将顶帽滤波图像添加到灰度图像中，使得絮体区域变得更亮和更大，随后减去底帽滤波图像使得细丝更暗，整个图像絮体和丝状菌轮廓清晰，便于下一步进行边缘提取。对滤波后的图像进行sobel梯度边缘检测，边缘检测后的图像由于噪声的干扰，其边界特征很少能够被完整地描述，在亮度不一致的地方会中断，采用形态学闭操作将整个边界进行连接。随后将得到的二值图和原灰度图进行相乘操作，去除杂乱背景；s32，利用约束模糊聚类的方法对去除杂乱背景后的图像进行分割，约束模糊聚类数目三类分别为整体、絮体以及丝状菌；经过聚类，每一类的像素点被归类到相应的标签，即标签1、2、3，提取丝状菌的标签3，得到的丝状菌的分割图像；s33，对分割后的图像进行后处理操作。经过s32分割后的图像由于一些目标被误检测，因此还需要进一步去碎片。第一阶段是将面积小于300的区域去除，第二阶段是利用减小的回转半径来区分碎片和丝状菌，rg大于1的物体是丝状菌，得到所需的丝状菌图像，第三阶段是去除老化絮体,采集的污泥图像中往往会碰到老化的絮体，老化的絮体呈黑褐色，同丝状菌的颜色相似，在聚类分割时，容易将老化絮体识别成丝状菌，这一阶段是判断和絮体模型分割出的絮体图像的像素，如果像素点，则认为此像素点絮体的像素点，减去属于絮体的部分，即，得到最终的丝状菌。
[0015]
s4，特征提取。其中，s4包括提取污泥絮体和丝状菌形态学特征，经过步骤（s2-s3），污泥图像中的絮体和丝状菌与背景分割开，分别计算出絮体和丝状菌的尺寸和形状等形态学特征参数，包括全局特征和局部特征。全局特征包括絮状物和丝状菌的数量及其总面积，局部特征包括尺寸度量和形状描述符。其中尺寸度量如面积、长度、宽度和周长等，形状描述符如圆度、形状系数、减小的回转半径。
[0016]
其中欧拉数，可作为识别物体的特征，欧拉数定义为连通成份数减去空洞数，这种拓扑特征具有平稳、旋转和比例不变特性；形状因子()主要用来于比较目标与圆之间的偏差，形状因子用来判断边界是否光滑，其取值范围为0到1，其中圆的形状因子为1；圆度，用r表示，圆度由其目标的伸长率决定，其取值范围为0到1，其中圆的圆度是1；回旋半径，用rg表示，回旋半径也由其目标的伸长率决定，越长絮体的旋转半径越大，圆的旋转半径是；纵横比，该特征表示物体形状的对称性，它用ar表示。其中每张图像中可能会有许多目标，对多个目标的图像计算出的上述参数（除总面积外）都取均值。采用22个参数定量表征絮体形态学特征、13个参数定量表征丝状菌形态学特征，即每个样本提取了35个形态学参数，这部分研究的目的是寻求图像分析信息和svi之间的相关性。
[0017]
本发明包括活性污泥显微图像采集、絮体分割、丝状菌分割、特征提取四个步骤，实施例具体实施进行了从某污水厂采集的污泥显微图像的实施分割，结果如图3所示。图3中，每行所示均为一张实验图像过程，在每列中，上图为原始采集污泥显微图像，中图为labelme标注的标签真值图（其中灰色为絮体，白色为丝状菌），下图为采用本发明处理获得
的分割结果图。
[0018]
本实施例最后使用333张污泥显微图像作为分割数据集，分割目标为絮体和丝状菌。使用图像分割的标准评判指标分割精度(accurary、precision、recall、dice)进行评测，图4给出了本方法与其他方法从污泥显微图像分割的测试结果，值越大，说明性能越好。
[0019]
图4本方法在污泥显微图像数据集上的性能，从图4可见，本方法在絮体和丝状菌的图像分割指标上均比其他方法高。
[0020]
以上所描述的装置实施例仅仅是示意性的，通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件的方式来实现，在不付出创造性的劳动的情况下，即可以理解并实施。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行各个实施例或者实施例的某些部分所述的方法。
[0021]
以上所述仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专利的技术人员在不脱离本发明技术方案范围内，当可利用上述提示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明方案的范围内。