CN106127264A - 一种控制水体系中微生物生长的方法 - Google Patents

Abstract

一种控制水体系中微生物生长的方法，包括以下步骤：(1)采用细胞识别模块对确定微生物种类；(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。本发明的有益效果为：能够有效控制微生物生长。

Description

一种控制水体系中微生物生长的方法

技术领域

本发明涉及微生物领域，具体涉及一种控制水体系中微生物生长的方法。

背景技术

为了控制工业制备用水中微生物引起的问题，众多抗微生物剂已经被用于消除、抑制或减少微生物生长。抗微生物剂被单独或组合使用以防止或控制由微生物生长引起的问题。通常将抗微生物剂直接加入制备用水流中或制备中使用的材料中。当用于防止生物膜形成时，典型的加入方法使得抗微生物剂分布遍及整个制备体系。这样，可以控制浮游微生物以及与制备用水接触的表面上的生物膜内的微生物。

发明内容

为解决上述问题，本发明旨在提供一种控制水体系中微生物生长的方法。

本发明的目的采用以下技术方案来实现：

一种控制水体系中微生物生长的方法，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

本发明的有益效果为：能够有效控制微生物生长。

附图说明

利用附图对本发明作进一步说明，但附图中的实施例不构成对本发明的任何限制，对于本领域的普通技术人员，在不付出创造性劳动的前提下，还可以根据以下附图获得其它的附图。

图1是本发明流程示意图；

图2是细胞识别模块的结构示意图。

附图标记：

细胞识别模块1、细胞图像分割单元11、特征提取单元12、分类识别单元13。

具体实施方式

结合以下应用场景对本发明作进一步描述。

应用场景1

参见图1、图2，本应用场景的一个实施例的一种控制水体系中微生物生长的方法，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

优选地，通过将铵或胺源与卤化氧化剂接触或者通过在卤素源存在下将该铵或胺源与氧化剂接触来制备一卤代胺。

本优选实施例制备方法简单，便于工业上实现。

优选地，所述一卤代胺包括一氯代胺。

本优选实施例节约成本，效果好。

优选的，所述细胞识别模块1包括细胞图像分割单元11、特征提取单元12、分类识别单元13；所述细胞图像分割单元11用于区分由细胞图像采集模块采集的细胞图像中的背景、细胞核和细胞质；所述特征提取单元12用于对细胞图像的纹理特征进行提取；所述分类识别单元13用于根据纹理特征利用分类器实现对细胞图像分类识别。

本优选实施例构建了细胞识别模块1的单元架构。

优选的，所述细胞图像分割单元11包括图像转换子单元、噪声去除子单元、粗分割子单元、细胞核中心标定子单元、精确分割子单元，具体为：

(1)图像转换子单元，用于将采集的细胞图像转化为灰度图像；

(2)噪声去除子单元，用于对灰度图像进行去噪处理，包括：

对于像素点(x,y)，选取其3×3的邻域S_x,y和(2N+1)×(2N+1)的邻域L_x,y，N为大于等于2的整数；

首先对像素点是否为边界点进行判断，设定阈值T，T∈[13,26]，计算像素点(x,y)与其邻域S_x,y中每个像素点的灰度差值，并与阈值T进行比较，若灰度差值大于阈值T的个数大于等于6，则像素点(x,y)为边界点，否则，像素点(x,y)为非边界点；

若(x,y)为边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{q(i,j)\∈\[q(x,y)-1.5\mathrm{\σ},q(x,y)+1.5\mathrm{\σ}\]}q(i,j)}{k}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(x,y)为降噪前像素点(x,y)的灰度值，σ为像素点(x,y)邻域L_x,y内灰度值标差,q(i,j)∈[q(x,y)-1.5σ,q(x,y)+1.5σ]表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点，k表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点的数量；

若(x,y)为非边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)q(i,j)}{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(i,j)代表图像中点(i,j)处的灰度值，w(i,j)为邻域L_x,y内点(i,j)对应的高斯权重；

(3)粗分割子单元，用于对去噪后的细胞图像中的背景、细胞质、细胞核进行粗划分，具体为：

将每个像素(x,y)用一个四维特征向量表示：

$\stackrel{\→}{u}(x,y)=\[h(x,y),h_{ave}(x,y),h_{med}(x,y),h_{sta}(x,y)\]$

式中，h(x,y)代表(x,y)的灰度值，h_ave(x,y)代表其邻域S_x,y灰度均值，h_med(x,y)代表其邻域S_x,y灰度中值，h_sta(x,y)代表其邻域S_x,y灰度方差；

采用K-means聚类法将其划分为背景、细胞质、细胞核三类；

(4)细胞核中心标定子单元，用于对细胞核中心进行标定：

由粗分割子单元得到细胞核大致区域，设细胞核区域包含n个点：(x₁,y₁),…,(x_n,y_n)，对该区域进行灰度加权标定和几何中心标定，取其平均值作为细胞核中心(x_z,y_z)：

$x_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{x}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{x}_i}{n}\right)$

$y_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{y}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{y}_i}{n}\right)$

(5)精确分割子单元，用于对细胞核、细胞质进行精确分割；

构建从细胞核中心(x_z,y_z)到细胞核和细胞质边界点(x_p,y_p)的有向线段的距离表示向下取整；

沿线段以单位长度进行采样可以得到dis_p个点若采样点的坐标不是整数，其灰度值通过周围像素线性插值得到；

点(x_i,y_i)处沿线段方向的灰度差：

hd(x_i,y_i)＝h(x_i-1,y_i-1)-h(x_i,y_i)

定义灰度差抑制函数：

$Y\left(x\right)=\left\{\begin{array}{ccc}x& if& x\≤0\\ 0.5x& if& x>0\end{array}\right.$

点(x_i,y_i)处沿线段方向的梯度gra(x_i,y_i)：

$gra(x_i,y_i)=\frac{\left|Y\left(hd\right(x_i,y_i\left)\right)\right|+\left|Y\left(hd\right(x_{i+1},y_{i+!}\left)\right)\right|}{2}$

选取梯度最大的值点作为细胞核和细胞质的精确边缘。

本优选实施例设置噪声去除子单元，有效融合了中心像素与邻域像素的空间临近性和灰度相似性来进行降噪处理，在图像中的平坦区域，邻域内像素灰度值相差不大，采用高斯滤波器对灰度值进行加权滤波，在变化剧烈的边界区域，行边界保持滤波，有利于图像边缘的保持；采用K均值聚类提取细胞核和细胞质粗轮廓，可有效去除噪声的干扰；设置细胞核中心标定子单元，便于后续对细胞核和细胞质轮廓进行精确定位；精确分割子单元充分利用了方向信息，克服了炎症细胞对边缘图的干扰，能够准确提取出细胞核和细胞质边缘。

优选的，所述对细胞图像的纹理特征进行提取，包括：

(1)基于改进的灰度共生矩阵法求取细胞图像的综合灰度共生矩阵，所述综合灰度共生矩阵体现了细胞在不同方向上的纹理特征：

设在0°、45°、90°、135°四个方向上的灰度共生矩阵分别为h(x,y,d,0°)、h(x,y,d,45°)、h(x,y,d,90°)、h(x,y,d,135°)，所对应的矩阵元素项目为X₁、X₂、X₃、X₄，则综合灰度共生矩阵的计算公式为：

H(x,y,d)＝w₁h(x,y,d,0°)+w₂h(x,y,d,45°)+w₃h(x,y,d,90°)+w₄h(x,y,d,135°)

综合灰度共生矩阵元素数目为：

$X=\underset{i=1}{\overset{4}{\Σ}}{w}_i{X}_i$

式中，d表示距离，d的取值范围为[2,4]，w_i为加权系数,i＝1,2,3,4，其由四个方向中的每个方向上的灰度共生矩阵对应的对比度参数计算，设四个方向上的灰度共生矩阵对应的对比度参数分别为D_i，均值为i＝1,2,3,4，则加权系数w_i的计算公式为：

$w_i=\frac{1}{|D_i-\stackrel{\‾}{D}|+1}/\underset{i=1}{\overset{4}{\Σ}}\frac{1}{|D_i-\stackrel{\‾}{D}|+1}$

(2)利用所述综合灰度共生矩阵和矩阵元素项目获取所需的四个纹理特征参数：对比度、方差和、能量和均值；

(3)对所述四个纹理特征参数进行归一化处理，最终获得归一化的纹理特征值。

本优选实施例基于改进的灰度共生矩阵法，采用设置加权系数的方式求取细胞图像的综合灰度共生矩阵，进而提取细胞在指定四个方向上的纹理特征，解决了由于外部干扰(如细胞图像采集时光照角度造成的影响、气体的流动干扰等)造成的细胞的纹理特征参数值在不同方向上有较大差别的问题，提高了细胞图像纹理特征提取的精度；选定对比度、方差和、能量和均值四个纹理特征，去掉了冗余和重复的特征参数；对所述四个纹理特征参数进行归一化处理，方便了后续的细胞图像的分类识别处理。

在此应用场景中，设定阈值T＝13，d＝2，图像去噪效果相对提高了5％，细胞图像特征的提取精度提高了8％。

应用场景2

参见图1、图2，本应用场景的一个实施例的一种控制水体系中微生物生长的方法，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

优选地，通过将铵或胺源与卤化氧化剂接触或者通过在卤素源存在下将该铵或胺源与氧化剂接触来制备一卤代胺。

本优选实施例制备方法简单，便于工业上实现。

优选地，所述一卤代胺包括一氯代胺。

本优选实施例节约成本，效果好。

优选的，所述细胞识别模块1包括细胞图像分割单元11、特征提取单元12、分类识别单元13；所述细胞图像分割单元11用于区分由细胞图像采集模块采集的细胞图像中的背景、细胞核和细胞质；所述特征提取单元12用于对细胞图像的纹理特征进行提取；所述分类识别单元13用于根据纹理特征利用分类器实现对细胞图像分类识别。

本优选实施例构建了细胞识别模块1的单元架构。

优选的，所述细胞图像分割单元11包括图像转换子单元、噪声去除子单元、粗分割子单元、细胞核中心标定子单元、精确分割子单元，具体为：

(1)图像转换子单元，用于将采集的细胞图像转化为灰度图像；

(2)噪声去除子单元，用于对灰度图像进行去噪处理，包括：

对于像素点(x,y)，选取其3×3的邻域S_x,y和(2N+1)×(2N+1)的邻域L_x,y，N为大于等于2的整数；

首先对像素点是否为边界点进行判断，设定阈值T，T∈[13,26]，计算像素点(x,y)与其邻域S_x,y中每个像素点的灰度差值，并与阈值T进行比较，若灰度差值大于阈值T的个数大于等于6，则像素点(x,y)为边界点，否则，像素点(x,y)为非边界点；

若(x,y)为边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{q(i,j)\∈\[q(x,y)-1.5\mathrm{\σ},q(x,y)+1.5\mathrm{\σ}\]}q(i,j)}{k}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(x,y)为降噪前像素点(x,y)的灰度值，σ为像素点(x,y)邻域L_x,y内灰度值标差,q(i,j)∈[q(x,y)-1.5σ,q(x,y)+1.5σ]表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点，k表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点的数量；

若(x,y)为非边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)q(i,j)}{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(i,j)代表图像中点(i,j)处的灰度值，w(i,j)为邻域L_x,y内点(i,j)对应的高斯权重；

(3)粗分割子单元，用于对去噪后的细胞图像中的背景、细胞质、细胞核进行粗划分，具体为：

将每个像素(x,y)用一个四维特征向量表示：

$\stackrel{\→}{u}(x,y)=\[h(x,y),h_{ave}(x,y),h_{med}(x,y),h_{sta}(x,y)\]$

式中，h(x,y)代表(x,y)的灰度值，h_ave(x,y)代表其邻域S_x,y灰度均值，h_med(x,y)代表其邻域S_x,y灰度中值，h_sta(x,y)代表其邻域S_x,y灰度方差；

采用K-means聚类法将其划分为背景、细胞质、细胞核三类；

(4)细胞核中心标定子单元，用于对细胞核中心进行标定：

由粗分割子单元得到细胞核大致区域，设细胞核区域包含n个点：(x₁,y₁),…,(x_n,y_n)，对该区域进行灰度加权标定和几何中心标定，取其平均值作为细胞核中心(x_z,y_z)：

$x_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{x}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{x}_i}{n}\right)$

$y_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{y}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{y}_i}{n}\right)$

(5)精确分割子单元，用于对细胞核、细胞质进行精确分割；

构建从细胞核中心(x_z,y_z)到细胞核和细胞质边界点(x_p,y_p)的有向线段的距离表示向下取整；

沿线段以单位长度进行采样可以得到dis_p个点若采样点的坐标不是整数，其灰度值通过周围像素线性插值得到；

点(x_i,y_i)处沿线段方向的灰度差：

hd(x_i,y_i)＝h(x_i-1,y_i-1)-h(x_i,y_i)

定义灰度差抑制函数：

$Y\left(x\right)=\left\{\begin{array}{ccc}x& if& x\≤0\\ 0.5x& if& x>0\end{array}\right.$

点(x_i,y_i)处沿线段方向的梯度gra(x_i,y_i)：

$gra(x_i,y_i)=\frac{\left|Y\left(hd\right(x_i,y_i\left)\right)\right|+\left|Y\left(hd\right(x_{i+1},y_{i+!}\left)\right)\right|}{2}$

选取梯度最大的值点作为细胞核和细胞质的精确边缘。

本优选实施例设置噪声去除子单元，有效融合了中心像素与邻域像素的空间临近性和灰度相似性来进行降噪处理，在图像中的平坦区域，邻域内像素灰度值相差不大，采用高斯滤波器对灰度值进行加权滤波，在变化剧烈的边界区域，行边界保持滤波，有利于图像边缘的保持；采用K均值聚类提取细胞核和细胞质粗轮廓，可有效去除噪声的干扰；设置细胞核中心标定子单元，便于后续对细胞核和细胞质轮廓进行精确定位；精确分割子单元充分利用了方向信息，克服了炎症细胞对边缘图的干扰，能够准确提取出细胞核和细胞质边缘。

优选的，所述对细胞图像的纹理特征进行提取，包括：

(1)基于改进的灰度共生矩阵法求取细胞图像的综合灰度共生矩阵，所述综合灰度共生矩阵体现了细胞在不同方向上的纹理特征：

设在0°、45°、90°、135°四个方向上的灰度共生矩阵分别为h(x,y,d,0°)、h(x,y,d,45°)、h(x,y,d,90°)、h(x,y,d,135°)，所对应的矩阵元素项目为X₁、X₂、X₃、X₄，则综合灰度共生矩阵的计算公式为：

H(x,y,d)＝w₁h(x,y,d,0°)+w₃h(x,y,d,45°)+w₃h(x,y,d,90°)+w₄h(x,y,d,135°)

综合灰度共生矩阵元素数目为：

$X=\underset{i=1}{\overset{4}{\Σ}}{w}_i{X}_i$

式中，d表示距离，d的取值范围为[2,4]，w_i为加权系数,i＝1,2,3,4，其由四个方向中的每个方向上的灰度共生矩阵对应的对比度参数计算，设四个方向上的灰度共生矩阵对应的对比度参数分别为D_i，均值为i＝1,2,3,4，则加权系数w_i的计算公式为：

$w_i=\frac{1}{|D_i-\stackrel{\‾}{D}|+1}/\underset{i=1}{\overset{4}{\Σ}}\frac{1}{|D_i-\stackrel{\‾}{D}|+1}$

(2)利用所述综合灰度共生矩阵和矩阵元素项目获取所需的四个纹理特征参数：对比度、方差和、能量和均值；

(3)对所述四个纹理特征参数进行归一化处理，最终获得归一化的纹理特征值。

本优选实施例基于改进的灰度共生矩阵法，采用设置加权系数的方式求取细胞图像的综合灰度共生矩阵，进而提取细胞在指定四个方向上的纹理特征，解决了由于外部干扰(如细胞图像采集时光照角度造成的影响、气体的流动干扰等)造成的细胞的纹理特征参数值在不同方向上有较大差别的问题，提高了细胞图像纹理特征提取的精度；选定对比度、方差和、能量和均值四个纹理特征，去掉了冗余和重复的特征参数；对所述四个纹理特征参数进行归一化处理，方便了后续的细胞图像的分类识别处理。

在此应用场景中，设定阈值T＝15，d＝2，图像去噪效果相对提高了6％，细胞图像特征的提取精度提高了8％。

应用场景3

参见图1、图2，本应用场景的一个实施例的一种控制水体系中微生物生长的方法，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

优选地，通过将铵或胺源与卤化氧化剂接触或者通过在卤素源存在下将该铵或胺源与氧化剂接触来制备一卤代胺。

本优选实施例制备方法简单，便于工业上实现。

优选地，所述一卤代胺包括一氯代胺。

本优选实施例节约成本，效果好。

优选的，所述细胞识别模块1包括细胞图像分割单元11、特征提取单元12、分类识别单元13；所述细胞图像分割单元11用于区分由细胞图像采集模块采集的细胞图像中的背景、细胞核和细胞质；所述特征提取单元12用于对细胞图像的纹理特征进行提取；所述分类识别单元13用于根据纹理特征利用分类器实现对细胞图像分类识别。

本优选实施例构建了细胞识别模块1的单元架构。

优选的，所述细胞图像分割单元11包括图像转换子单元、噪声去除子单元、粗分割子单元、细胞核中心标定子单元、精确分割子单元，具体为：

(1)图像转换子单元，用于将采集的细胞图像转化为灰度图像；

(2)噪声去除子单元，用于对灰度图像进行去噪处理，包括：

对于像素点(x,y)，选取其3×3的邻域S_x,y和(2N+1)×(2N+1)的邻域L_x,y，N为大于等于2的整数；

首先对像素点是否为边界点进行判断，设定阈值T，T∈[13,26]，计算像素点(x,y)与其邻域S_x,y中每个像素点的灰度差值，并与阈值T进行比较，若灰度差值大于阈值T的个数大于等于6，则像素点(x,y)为边界点，否则，像素点(x,y)为非边界点；

若(x,y)为边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{q(i,j)\∈\[q(x,y)-1.5\mathrm{\σ},q(x,y)+1.5\mathrm{\σ}\]}q(i,j)}{k}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(x,y)为降噪前像素点(x,y)的灰度值，σ为像素点(x,y)邻域L_x,y内灰度值标差,q(i,j)∈[q(x,y)-1.5σ,q(x,y)+1.5σ]表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点，k表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点的数量；

若(x,y)为非边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)q(i,j)}{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(i,j)代表图像中点(i,j)处的灰度值，w(i,j)为邻域L_x,y内点(i,j)对应的高斯权重；

(3)粗分割子单元，用于对去噪后的细胞图像中的背景、细胞质、细胞核进行粗划分，具体为：

将每个像素(x,y)用一个四维特征向量表示：

$\stackrel{\→}{u}(x,y)=\[h(x,y),h_{ave}(x,y),h_{med}(x,y),h_{sta}(x,y)\]$

式中，h(x,y)代表(x，y)的灰度值，h_ave(x,y)代表其邻域S_x,y灰度均值，h_med(x,y)代表其邻域S_x,y灰度中值，h_sta(x,y)代表其邻域S_x,y灰度方差；

采用K-means聚类法将其划分为背景、细胞质、细胞核三类；

(4)细胞核中心标定子单元，用于对细胞核中心进行标定：

由粗分割子单元得到细胞核大致区域，设细胞核区域包含n个点：(x₁,y₁),…,(x_n,y_n)，对该区域进行灰度加权标定和几何中心标定，取其平均值作为细胞核中心(x_z,y_z)：

$x_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{x}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{x}_i}{n}\right)$

$y_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{y}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{y}_i}{n}\right)$

(5)精确分割子单元，用于对细胞核、细胞质进行精确分割；

构建从细胞核中心(x_z,y_z)到细胞核和细胞质边界点(x_p,y_p)的有向线段的距离表示向下取整；

沿线段以单位长度进行采样可以得到dis_p个点若采样点的坐标不是整数，其灰度值通过周围像素线性插值得到；

点(x_i,y_i)处沿线段方向的灰度差：

hd(x_i,y_i)＝h(x_i-1,y_i-1)-h(x_i,y_i)

定义灰度差抑制函数：

$Y\left(x\right)=\left\{\begin{array}{ccc}x& if& x\≤0\\ 0.5x& if& x>0\end{array}\right.$

点(x_i,y_i)处沿线段方向的梯度gra(x_i,y_i)：

$gra(x_i,y_i)=\frac{\left|Y\left(hd\right(x_i,y_i\left)\right)\right|+\left|Y\left(hd\right(x_{i+1},y_{i+!}\left)\right)\right|}{2}$

选取梯度最大的值点作为细胞核和细胞质的精确边缘。

本优选实施例设置噪声去除子单元，有效融合了中心像素与邻域像素的空间临近性和灰度相似性来进行降噪处理，在图像中的平坦区域，邻域内像素灰度值相差不大，采用高斯滤波器对灰度值进行加权滤波，在变化剧烈的边界区域，行边界保持滤波，有利于图像边缘的保持；采用K均值聚类提取细胞核和细胞质粗轮廓，可有效去除噪声的干扰；设置细胞核中心标定子单元，便于后续对细胞核和细胞质轮廓进行精确定位；精确分割子单元充分利用了方向信息，克服了炎症细胞对边缘图的干扰，能够准确提取出细胞核和细胞质边缘。

优选的，所述对细胞图像的纹理特征进行提取，包括：

(1)基于改进的灰度共生矩阵法求取细胞图像的综合灰度共生矩阵，所述综合灰度共生矩阵体现了细胞在不同方向上的纹理特征：

设在0°、45°、90°、135°四个方向上的灰度共生矩阵分别为h(x,y,d,0°)、h(x,y,d,45°)、h(x,y,d,90°)、h(x,y,d,135°)，所对应的矩阵元素项目为X₁、X₂、X₃、X₄，则综合灰度共生矩阵的计算公式为：

H(x,y,d)＝w₁h(x,y,d,0°)+w₂h(x,y,d,45°)+w₃h(x,y,d,90°)+w₄h(x,y,d,135°)

综合灰度共生矩阵元素数目为：

$X=\underset{i=1}{\overset{4}{\Σ}}{w}_i{X}_i$

式中，d表示距离，d的取值范围为[2,4]，w_i为加权系数,i＝1,2,3,4，其由四个方向中的每个方向上的灰度共生矩阵对应的对比度参数计算，设四个方向上的灰度共生矩阵对应的对比度参数分别为D_i，均值为i＝1,2,3,4，则加权系数w_i的计算公式为：

$w_i=\frac{1}{|D_i-\stackrel{\‾}{D}|+1}/\underset{i=1}{\overset{4}{\Σ}}\frac{1}{|D_i-\stackrel{\‾}{D}|+1}$

(2)利用所述综合灰度共生矩阵和矩阵元素项目获取所需的四个纹理特征参数：对比度、方差和、能量和均值；

(3)对所述四个纹理特征参数进行归一化处理，最终获得归一化的纹理特征值。

本优选实施例基于改进的灰度共生矩阵法，采用设置加权系数的方式求取细胞图像的综合灰度共生矩阵，进而提取细胞在指定四个方向上的纹理特征，解决了由于外部干扰(如细胞图像采集时光照角度造成的影响、气体的流动干扰等)造成的细胞的纹理特征参数值在不同方向上有较大差别的问题，提高了细胞图像纹理特征提取的精度；选定对比度、方差和、能量和均值四个纹理特征，去掉了冗余和重复的特征参数；对所述四个纹理特征参数进行归一化处理，方便了后续的细胞图像的分类识别处理。

在此应用场景中，设定阈值T＝18，d＝3，图像去噪效果相对提高了7％，细胞图像特征的提取精度提高了7％。

应用场景4

参见图1、图2，本应用场景的一个实施例的一种控制水体系中微生物生长的方法，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

优选地，通过将铵或胺源与卤化氧化剂接触或者通过在卤素源存在下将该铵或胺源与氧化剂接触来制备一卤代胺。

本优选实施例制备方法简单，便于工业上实现。

优选地，所述一卤代胺包括一氯代胺。

本优选实施例节约成本，效果好。

优选的，所述细胞识别模块1包括细胞图像分割单元11、特征提取单元12、分类识别单元13；所述细胞图像分割单元11用于区分由细胞图像采集模块采集的细胞图像中的背景、细胞核和细胞质；所述特征提取单元12用于对细胞图像的纹理特征进行提取；所述分类识别单元13用于根据纹理特征利用分类器实现对细胞图像分类识别。

本优选实施例构建了细胞识别模块1的单元架构。

优选的，所述细胞图像分割单元11包括图像转换子单元、噪声去除子单元、粗分割子单元、细胞核中心标定子单元、精确分割子单元，具体为：

(1)图像转换子单元，用于将采集的细胞图像转化为灰度图像；

(2)噪声去除子单元，用于对灰度图像进行去噪处理，包括：

对于像素点(x,y)，选取其3×3的邻域S_x,y和(2N+1)×(2N+1)的邻域L_x,y，N为大于等于2的整数；

首先对像素点是否为边界点进行判断，设定阈值T，T∈[13,26]，计算像素点(x,y)与其邻域S_x,y中每个像素点的灰度差值，并与阈值T进行比较，若灰度差值大于阈值T的个数大于等于6，则像素点(x,y)为边界点，否则，像素点(x,y)为非边界点；

若(x,y)为边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{q(i,j)\∈\[q(x,y)-1.5\mathrm{\σ},q(x,y)+1.5\mathrm{\σ}\]}q(i,j)}{k}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(x,y)为降噪前像素点(x,y)的灰度值，σ为像素点(x,y)邻域L_x,y内灰度值标差,q(i,j)∈[q(x,y)-1.5σ,q(x,y)+1.5σ]表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点，k表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点的数量；

若(x,y)为非边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)q(i,j)}{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(i,j)代表图像中点(i,j)处的灰度值，w(i,j)为邻域L_x,y内点(i,j)对应的高斯权重；

(3)粗分割子单元，用于对去噪后的细胞图像中的背景、细胞质、细胞核进行粗划分，具体为：

将每个像素(x,y)用一个四维特征向量表示：

$\stackrel{\→}{u}(x,y)=\[h(x,y),h_{ave}(x,y),h_{med}(x,y),h_{sta}(x,y)\]$

式中，h(x,y)代表(x,y)的灰度值，h_ave(x,y)代表其邻域S_x,y灰度均值，h_med(x,y)代表其邻域S_x,y灰度中值，h_sta(x,y)代表其邻域S_x,y灰度方差；

采用K-means聚类法将其划分为背景、细胞质、细胞核三类；

(4)细胞核中心标定子单元，用于对细胞核中心进行标定：

由粗分割子单元得到细胞核大致区域，设细胞核区域包含n个点：(x₁,y₁),…,(x_n,y_n)，对该区域进行灰度加权标定和几何中心标定，取其平均值作为细胞核中心(x_z,y_z)：

$x_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{x}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{x}_i}{n}\right)$

$y_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{y}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{y}_i}{n}\right)$

(5)精确分割子单元，用于对细胞核、细胞质进行精确分割；

构建从细胞核中心(x_z,y_z)到细胞核和细胞质边界点(x_p,y_p)的有向线段的距离表示向下取整；

沿线段以单位长度进行采样可以得到dis_p个点若采样点的坐标不是整数，其灰度值通过周围像素线性插值得到；

点(x_i,y_i)处沿线段方向的灰度差：

hd(x_i,y_i)＝h(x_i-1,y_i-1)-h(x_i,y_i)

定义灰度差抑制函数：

$Y\left(x\right)=\left\{\begin{array}{ccc}x& if& x\≤0\\ 0.5x& if& x>0\end{array}\right.$

点(x_i,y_i)处沿线段方向的梯度gra(x_i,y_i)：

$gra(x_i,y_i)=\frac{\left|Y\left(hd\right(x_i,y_i\left)\right)\right|+\left|Y\left(hd\right(x_{i+1},y_{i+!}\left)\right)\right|}{2}$

选取梯度最大的值点作为细胞核和细胞质的精确边缘。

本优选实施例设置噪声去除子单元，有效融合了中心像素与邻域像素的空间临近性和灰度相似性来进行降噪处理，在图像中的平坦区域，邻域内像素灰度值相差不大，采用高斯滤波器对灰度值进行加权滤波，在变化剧烈的边界区域，行边界保持滤波，有利于图像边缘的保持；采用K均值聚类提取细胞核和细胞质粗轮廓，可有效去除噪声的干扰；设置细胞核中心标定子单元，便于后续对细胞核和细胞质轮廓进行精确定位；精确分割子单元充分利用了方向信息，克服了炎症细胞对边缘图的干扰，能够准确提取出细胞核和细胞质边缘。

优选的，所述对细胞图像的纹理特征进行提取，包括：

(1)基于改进的灰度共生矩阵法求取细胞图像的综合灰度共生矩阵，所述综合灰度共生矩阵体现了细胞在不同方向上的纹理特征：

设在0°、45°、90°、135°四个方向上的灰度共生矩阵分别为h(x,y,d,0°)、h(x,y,d,45°)、h(x,y,d,90°)、h(x,y,d,135°)，所对应的矩阵元素项目为X₁、X₂、X₃、X₄，则综合灰度共生矩阵的计算公式为：

H(x,y,d)＝w₁h(x,y,d,0°)+w₂h(x,y,d,45°)+w₃h(x,y,d,90°)+w₄h(x,y,d,135°)

综合灰度共生矩阵元素数目为：

$X=\underset{i=1}{\overset{4}{\Σ}}{w}_i{X}_i$

式中，d表示距离，d的取值范围为[2,4]，w_i为加权系数,i＝1,2,3,4，其由四个方向中的每个方向上的灰度共生矩阵对应的对比度参数计算，设四个方向上的灰度共生矩阵对应的对比度参数分别为D_i，均值为i＝1,2,3,4，则加权系数w_i的计算公式为：

$w_i=\frac{1}{|D_i-\stackrel{\‾}{D}|+1}/\underset{i=1}{\overset{4}{\Σ}}\frac{1}{|D_i-\stackrel{\‾}{D}|+1}$

(2)利用所述综合灰度共生矩阵和矩阵元素项目获取所需的四个纹理特征参数：对比度、方差和、能量和均值；

(3)对所述四个纹理特征参数进行归一化处理，最终获得归一化的纹理特征值。

本优选实施例基于改进的灰度共生矩阵法，采用设置加权系数的方式求取细胞图像的综合灰度共生矩阵，进而提取细胞在指定四个方向上的纹理特征，解决了由于外部干扰(如细胞图像采集时光照角度造成的影响、气体的流动干扰等)造成的细胞的纹理特征参数值在不同方向上有较大差别的问题，提高了细胞图像纹理特征提取的精度；选定对比度、方差和、能量和均值四个纹理特征，去掉了冗余和重复的特征参数；对所述四个纹理特征参数进行归一化处理，方便了后续的细胞图像的分类识别处理。

在此应用场景中，设定阈值T＝20，d＝4，图像去噪效果相对提高了8％，细胞图像特征的提取精度提高了6％。

应用场景5

参见图1、图2，本应用场景的一个实施例的一种控制水体系中微生物生长的方法，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

优选地，通过将铵或胺源与卤化氧化剂接触或者通过在卤素源存在下将该铵或胺源与氧化剂接触来制备一卤代胺。

本优选实施例制备方法简单，便于工业上实现。

优选地，所述一卤代胺包括一氯代胺。

本优选实施例节约成本，效果好。

优选的，所述细胞识别模块1包括细胞图像分割单元11、特征提取单元12、分类识别单元13；所述细胞图像分割单元11用于区分由细胞图像采集模块采集的细胞图像中的背景、细胞核和细胞质；所述特征提取单元12用于对细胞图像的纹理特征进行提取；所述分类识别单元13用于根据纹理特征利用分类器实现对细胞图像分类识别。

本优选实施例构建了细胞识别模块1的单元架构。

优选的，所述细胞图像分割单元11包括图像转换子单元、噪声去除子单元、粗分割子单元、细胞核中心标定子单元、精确分割子单元，具体为：

(1)图像转换子单元，用于将采集的细胞图像转化为灰度图像；

(2)噪声去除子单元，用于对灰度图像进行去噪处理，包括：

对于像素点(x,y)，选取其3×3的邻域S_x,y和(2N+1)×(2N+1)的邻域L_x,y，N为大于等于2的整数；

首先对像素点是否为边界点进行判断，设定阈值T，T∈[13,26]，计算像素点(x,y)与其邻域S_x,y中每个像素点的灰度差值，并与阈值T进行比较，若灰度差值大于阈值T的个数大于等于6，则像素点(x,y)为边界点，否则，像素点(x,y)为非边界点；

若(x,y)为边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{q(i,j)\∈\[q(x,y)-1.5\mathrm{\σ},q(x,y)+1.5\mathrm{\σ}\]}q(i,j)}{k}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(x,y)为降噪前像素点(x,y)的灰度值，σ为像素点(x,y)邻域L_x,y内灰度值标差,q(i,j)∈[q(x,y)-1.5σ,q(x,y)+1.5σ]表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点，k表示邻域L_x,y内灰度值落于区间[q(x,y)-1.5σ,q(x,y)+1.5σ]的点的数量；

若(x,y)为非边界点，则进行如下降噪处理：

$h(x,y)=\frac{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)q(i,j)}{{\mathrm{\Σ}}_{(i,j)\∈L_{x,y}}w(i,j)}$

式中，h(x,y)为降噪后像素点(x,y)的灰度值，q(i,j)代表图像中点(i,j)处的灰度值，w(i,j)为邻域L_x,y内点(i,j)对应的高斯权重；

(3)粗分割子单元，用于对去噪后的细胞图像中的背景、细胞质、细胞核进行粗划分，具体为：

将每个像素(x,y)用一个四维特征向量表示：

$\stackrel{\→}{u}(x,y)=\[h(x,y),h_{ave}(x,y),h_{med}(x,y),h_{sta}(x,y)\]$

式中，h(x,y)代表(x,y)的灰度值，h_ave(x,y)代表其邻域S_x,y灰度均值，h_med(x,y)代表其邻域S_x,y灰度中值，h_sta(x,y)代表其邻域S_x,y灰度方差；

采用K-means聚类法将其划分为背景、细胞质、细胞核三类；

(4)细胞核中心标定子单元，用于对细胞核中心进行标定：

由粗分割子单元得到细胞核大致区域，设细胞核区域包含n个点：(x₁,y₁),…,(x_n,y_n)，对该区域进行灰度加权标定和几何中心标定，取其平均值作为细胞核中心(x_z,y_z)：

$x_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{x}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{x}_i}{n}\right)$

$y_z=\frac{1}{2}\left(\frac{{\Σ}_{i=1}^n{y}_{i}h(x_i,y_i)}{{\Σ}_{i=1}^{n}h(x_i,y_i)}+\frac{{\Σ}_{i=1}^n{y}_i}{n}\right)$

(5)精确分割子单元，用于对细胞核、细胞质进行精确分割；

构建从细胞核中心(x_z,y_z)到细胞核和细胞质边界点(x_p,y_p)的有向线段的距离表示向下取整；

沿线段以单位长度进行采样可以得到dis_p个点若采样点的坐标不是整数，其灰度值通过周围像素线性插值得到；

点(x_i,y_i)处沿线段方向的灰度差：

hd(x_i,y_i)＝h(x_i-1,y_i-1)-h(x_i,y_i)

定义灰度差抑制函数：

$Y\left(x\right)=\left\{\begin{array}{ccc}x& if& x\≤0\\ 0.5x& if& x>0\end{array}\right.$

点(x_i,y_i)处沿线段方向的梯度gra(x_i,y_i)：

$gra(x_i,y_i)=\frac{\left|Y\left(hd\right(x_i,y_i\left)\right)\right|+\left|Y\left(hd\right(x_{i+1},y_{i+!}\left)\right)\right|}{2}$

选取梯度最大的值点作为细胞核和细胞质的精确边缘。

本优选实施例设置噪声去除子单元，有效融合了中心像素与邻域像素的空间临近性和灰度相似性来进行降噪处理，在图像中的平坦区域，邻域内像素灰度值相差不大，采用高斯滤波器对灰度值进行加权滤波，在变化剧烈的边界区域，行边界保持滤波，有利于图像边缘的保持；采用K均值聚类提取细胞核和细胞质粗轮廓，可有效去除噪声的干扰；设置细胞核中心标定子单元，便于后续对细胞核和细胞质轮廓进行精确定位；精确分割子单元充分利用了方向信息，克服了炎症细胞对边缘图的干扰，能够准确提取出细胞核和细胞质边缘。

优选的，所述对细胞图像的纹理特征进行提取，包括：

(1)基于改进的灰度共生矩阵法求取细胞图像的综合灰度共生矩阵，所述综合灰度共生矩阵体现了细胞在不同方向上的纹理特征：

设在0°、45°、90°、135°四个方向上的灰度共生矩阵分别为h(x,y,d,0°)、h(x,y,d,45°)、h(x,y,d,90°)、h(x,y,d,135°)，所对应的矩阵元素项目为X₁、X₂、X₃、X₄，则综合灰度共生矩阵的计算公式为：

H(x,y,d)＝w₁h(x,y,d,0°)+w₂h(x,y,d,45°)+w₃h(x,y,d,90°)+w₄h(x,y,d,135°)

综合灰度共生矩阵元素数目为：

$X=\underset{i=1}{\overset{4}{\Σ}}{w}_i{X}_i$

式中，d表示距离，d的取值范围为[2,4]，w_i为加权系数,i＝1,2,3,4，其由四个方向中的每个方向上的灰度共生矩阵对应的对比度参数计算，设四个方向上的灰度共生矩阵对应的对比度参数分别为D_i，均值为i＝1,2,3,4，则加权系数w_i的计算公式为：

$w_i=\frac{1}{|D_i-\stackrel{\‾}{D}|+1}/\underset{i=1}{\overset{4}{\Σ}}\frac{1}{|D_i-\stackrel{\‾}{D}|+1}$

(2)利用所述综合灰度共生矩阵和矩阵元素项目获取所需的四个纹理特征参数：对比度、方差和、能量和均值；

(3)对所述四个纹理特征参数进行归一化处理，最终获得归一化的纹理特征值。

本优选实施例基于改进的灰度共生矩阵法，采用设置加权系数的方式求取细胞图像的综合灰度共生矩阵，进而提取细胞在指定四个方向上的纹理特征，解决了由于外部干扰(如细胞图像采集时光照角度造成的影响、气体的流动干扰等)造成的细胞的纹理特征参数值在不同方向上有较大差别的问题，提高了细胞图像纹理特征提取的精度；选定对比度、方差和、能量和均值四个纹理特征，去掉了冗余和重复的特征参数；对所述四个纹理特征参数进行归一化处理，方便了后续的细胞图像的分类识别处理。

在此应用场景中，设定阈值T＝26，d＝2，图像去噪效果相对提高了7.5％，细胞图像特征的提取精度提高了8％。

最后应当说明的是，以上实施例仅用以说明本发明的技术方案，而非对本发明保护范围的限制，尽管参照较佳实施例对本发明作了详细地说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和范围。

Claims (3)

1.一种控制水体系中微生物生长的方法，其特征是，包括以下步骤：

(1)采用细胞识别模块对确定微生物种类；

(2)向水体系中加入有效量的一卤代胺和二卤代胺的混合物，其中选择一卤代胺对二卤代胺的比率以得到体系的协同指数小于1。

2.根据权利要求1所述的一种控制水体系中微生物生长的方法，其特征是，通过将铵或胺源与卤化氧化剂接触或者通过在卤素源存在下将该铵或胺源与氧化剂接触来制备一卤代胺。

3.根据权利要求2所述的一种控制水体系中微生物生长的方法，其特征是，所述一卤代胺包括一氯代胺。