基于HSV颜色模型的植保雾滴检测方法

技术领域

本发明涉及农药喷洒效果检测领域，具体公开了一种基于HSV颜色模型的植保雾滴检测方法。

背景技术

在植物生长过程中，需要喷洒农药以保护植物免受虫害，此时，需要判断农药喷洒效果。常用的判断方法是先采集雾滴，然后对雾滴进行检测，计算出雾滴沉积密度、雾滴沉积覆盖率和雾滴粒径谱，再进行判断。现有雾滴检测方法主要有激光法、机械法、图像法，激光法操作较为便捷，但测量仪器成本价格昂贵，并且对检测的条件较为严格，具有一定的局限性；对于机械法来说，如油盘法，需要操作者进行大量手工操作，获取采集雾滴后，还需通过显微镜进行测定、计算、和矫正，操作过程复杂且精度很低；图像法又可分为传统图像法和基于深度学习的图像法，传统图像法在对雾滴进行图像提取时不能很好的去除背景信息，容易引入背景造成检测误差，基于深度学习的方法在实现方法上较为繁琐，需要预先训练网络模型，然后再通过模型进行检测，在检测过程中速度较慢且容易出现误检的情况。上述的激光法、机械法和图像法在某些特点场景下不能准确快速的对雾滴进行分析。

发明内容

有鉴于此，本发明的目的在于提供一种基于HSV颜色模型的植保雾滴检测方法。

为达到上述目的，本发明提供如下技术方案：

一种基于HSV颜色模型的植保雾滴检测方法，包括以下步骤：

步骤S1、使用农药喷洒设备将测试液喷洒到水敏纸上，在水敏纸上形成雾滴区域；

步骤S2、获取水敏纸的图像信息作为雾滴原图，并对雾滴原图进行RGB三通道拆分，得到R、G、B三个通道的图像信息；

步骤S3、提取出R通道的图像信息，使用R通道的图像信息代替G通道和B通道的图像信息，并将三个通道的图像信息合并得到R三通道模型；

步骤S4、将R三通道模型转化为HSV颜色模型，通过HSV颜色模型得到雾滴图像各部分对应的色相分量H、饱和度分量S和明度分量V，调节明度分量V的最大值，生成去除背景区域并保留雾滴区域的掩膜图；

步骤S5、将雾滴原图转化为BGR颜色模式并与掩膜图进行叠加，从而去除BGR颜色模式图像的背景区域，得到蒙版图；

步骤S6、对蒙版图依次进行灰度化操作和高斯滤波，得到高斯滤波图；

步骤S7、对高斯滤波图依次进行二值化操作、形态学操作和开运算去噪，得到雾滴去噪图；

步骤S8、对雾滴去噪图进行轮廓检测，提取出雾滴轮廓，得到雾滴轮廓图；

步骤S9、计算雾滴轮廓最小包围圆形的圆心及半径，并根据雾滴轮廓半径及圆心信息计算出雾滴沉积密度、雾滴沉积覆盖率和雾滴粒径谱。

进一步的，在所述步骤S4中，生成掩膜图时，先在水敏纸背景部分的明度分量值和雾滴部分的明度分量值之间选择一个值作为明度阈值V

进一步的，在所述步骤S6中，进行灰度化操作时，先定义BGR模式的图像中单个像素的灰度值为Y，然后通过下式计算每个像素的灰度值Y：

Y＝0.299R+0.587G+0.114B

其中，R表示单个像素的红色分量的值；G表示单个像素的绿色分量的值；B表示单个像素的蓝色分量的值。

进一步的，在所述步骤S7中，二值化操作的方法为：预设一个灰度阈值T

进一步的，在所述步骤S8中，进行轮廓检测并提取雾滴轮廓的方法为：压缩雾滴图像二值图的水平方向、垂直方向和对角线方向的元素，只保留该方向的终点坐标，得到雾滴的外侧轮廓信息，并通过轮廓信息绘制出雾滴轮廓。

进一步的，在执行完步骤S8后，先执行以下步骤：

步骤S801、预设一个最小面积阈值T

进一步的，雾滴沉积密度的计算公式如下：

上式中，K为雾滴沉积密度；N为过滤后雾滴的数量；A

进一步的，雾滴沉积覆盖率的计算公式如下：

上式中，C为雾滴沉积覆盖率；

本方案的工作原理及有益效果在于：本发明在传统图像法的基础上进行了改进，通过对R三通道进行合并生成掩膜图，采用蒙版的方法对雾滴进行检测，能够在雾滴数量密集的水敏纸上较为准确的统计出雾滴数量中径、体积中径、沉积密度、沉积覆盖率等衡量农药喷洒效果的重要技术指标，检测过程简单。

附图说明

图1为本发明基于HSV颜色模型的植保雾滴检测方法的一个优选实施例的流程图；

图2为步骤S2中所述的雾滴原图的示意图；

图3为步骤S2中所述的R、G、B三通道对比图的示意图；

图4为步骤S3中所述的R三通道图的示意图；

图5为步骤S4中所述的掩膜图的示意图；

图6为步骤S5中所述的蒙版图的示意图；

图7为步骤S6中所述的灰度图的示意图；

图8为步骤S6中所述的高斯滤波图的示意图；

图9为步骤S7中所述的雾滴二值图的示意图；

图10为步骤S7中所述的雾滴去噪图的示意图；

图11为步骤S8中所述的雾滴轮廓图的示意图；

图12为步骤S9中所述的检测结果图的示意图。

具体实施方式

下面通过具体实施方式进一步详细说明：

实施例

本发明基于HSV颜色模型的植保雾滴检测方法的一个优选实施例包括以下步骤：

步骤S1、使用农药喷洒设备将测试液喷洒到水敏纸上，在水敏纸上形成雾滴区域。

步骤S2、获取水敏纸的图像信息作为雾滴原图，如图2所示；对雾滴原图进行RGB三通道拆分，得到R、G、B三个通道的图像信息，如图3所示。

步骤S3、通过对比可知R通道图像能够更好地突出雾滴原图的背景与雾滴的差异，更有利于雾滴轮廓的提取，因此，提取出R通道的图像信息，并使用R通道的图像信息代替G通道和B通道的图像信息，并将三个通道的图像信息合并得到R三通道模型，生成雾滴的R三通道图片，如图4所示。

步骤S4、将R三通道模型转化为HSV颜色模型，通过HSV颜色模型得到雾滴图像各部分对应的H色相分量、S饱和度分量、V明度分量，调节明度分量V的最大值，生成去除背景区域并保留雾滴区域的掩膜图，如图5所示。具体方法为，在水敏纸背景部分的明度分量值和雾滴部分的明度分量值之间选择一个值作为明度阈值V

通过上式可将雾滴部分的颜色变为白色，把背景区域变为黑色，从而生成掩膜图。通过掩膜图可知其只能显示黑白两色，即目标区域为白色，背景为黑色，当然，也可将白色变成所需要的其他色彩，只需对图像进行合并叠加即可。

步骤S5、将雾滴原图转化为BGR颜色模式并与掩膜图进行叠加，从而去除BGR颜色模式图像的背景区域，得到蒙版图，如图6所示。具体方法为，利用掩膜屏蔽作用对图像上某些区域作屏蔽,使其不参加处理或不参加处理参数的计算。当两张图片进行叠加时，会在掩膜图的非屏蔽区域将原始图片保留并显示，在掩膜图的屏蔽区域显示遮罩部分。利用掩膜图与原图像进行“与”操作，即掩膜图像白色区域是对需要处理图像像素的保留，黑色区域是对需要处理图像像素的剔除，通过二者叠加能够有效的过滤水敏纸背景信息并最大限度的保持雾滴信息。

步骤S6、对蒙版图依次进行灰度化操作和高斯滤波。

去掉背景后需要对图像进一步进行灰度化，通过灰度化完成颜色空间从BGR模式转化为灰度模式，减小图像原始数据量，降低计算量加速后续处理，同时避免条带失真，保持图片信息的完整性。具体方法为：定义BGR模式的图像中单个像素的灰度值为Y，通过下式计算每个像素的灰度值Y：

Y＝0.299R+0.587G+0.114B

其中，R表示单个像素的红色分量的值；G表示单个像素的绿色分量的值；B表示单个像素的蓝色分量的值。通过灰度化操作，得到灰度图，如图7所示。

在数字图像后期应用中噪声会造成很大的误差，而误差在不同的处理操作中会累积传递，为了能够得到较好的图像，需要采用高斯滤波器除去高斯噪声，让图像更加平滑。高斯滤波器能够根据高斯函数的形状来选择权值的线性平滑滤波器，高斯滤波器对于服从正太分布的噪声非常有效，其计算方法为，建立二维的UV坐标系，定义输入图像为S(u，v)，通过二维高斯函数处理后得到G(u,v)，用一个模板扫描灰度图中的每一个像素，确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值，公式为：

上式中，σ为正态分布的标准偏差。在高斯滤波后的得到高斯滤波图，如图8所示。

步骤S7、对高斯滤波图依次进行二值化操作、形态学操作和开运算去噪。

进行二值化操作的具体方法为：预设一个灰度阈值T

形态学处理的核心就是定义结构元素，包括结构元素的大小及形状。对雾滴二值图进行形态学操作时，先定义一个3×3的十字形的结构元素

步骤S8、对雾滴去噪图进行轮廓检测，提取出雾滴轮廓；具体方法为：压缩雾滴图像二值图的水平方向、垂直方向和对角线方向的元素，只保留该方向的终点坐标，得到雾滴的外侧轮廓信息，并通过轮廓信息绘制出雾滴轮廓，生成雾滴轮廓图，如图11所示。

步骤S801、预设一个最小面积阈值T

步骤S9、通过过滤后的轮廓信息检测出雾滴轮廓最小包围圆形的圆心及半径信息，并在雾滴图片上圈出，得到检测结果图，如图12所示。计算雾滴轮廓最小包围圆形的圆心及半径，并根据雾滴轮廓半径及圆心信息计算出雾滴沉积密度、雾滴沉积覆盖率和雾滴粒径谱。

1、雾滴沉积密度

雾滴沉积密度K是影响雾滴防治效果的重要参数，是指雾滴收集装置单位面积上所承接的雾滴个数，单位为：个/cm

上式中，N为过滤后雾滴的数量；A

2、雾滴沉积覆盖率

雾滴沉积覆盖率C表示雾滴覆盖区域面积占统计总面积的百分比，其计算公式如下：

上式中，

3、雾滴粒径谱

雾滴粒径谱主要包括数量中径和体积中径；当然，雾滴粒径谱还包括其他参数，由于雾滴粒径谱的其他参数在判断农药喷洒效果时可以忽略，因此本实施例中只计算了数量中径和体积中径。

在一次喷雾样本中，将所有雾滴按照雾滴直径从小到大顺序累积，当累积到雾滴数量等于雾滴总数的50％时，所对应的雾滴直径即为此次雾滴样本中雾滴的数量中径，其单位为μm(微米)；数量中径NMD的计算公式表示为：

上式中，r

在一次喷雾样本中，将所有雾滴按照体积由小到大的顺序累积，当累积叠加到的体积值等于雾滴体积总和的50％时，所对应的雾滴直径即为此次雾滴样本中雾滴的体积中径，其单位为μm(微米)；体积中径VMD的计算公式表示为：

上式中，R表示体积中径对应的像素直径，其计算公式为：

上式中，V表示累积叠加的体积值；Vm表示所有雾滴体积之和的一半，其计算公式为：

之后，即可根据雾滴沉积密度、雾滴沉积覆盖率和雾滴粒径谱判断农药的喷洒效果。

以上所述的仅是本发明的实施例，方案中公知的具体结构及特性等常识在此未作过多描述。应当指出，对于本领域的技术人员来说，在不脱离本发明结构的前提下，还可以作出若干变形和改进，这些也应该视为本发明的保护范围，这些都不会影响本发明实施的效果和本发明的实用性。