首页 / 专利库 / 图形技术 / 颜色数据 / 一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法

一种基于图像处理模式识别技术的稻瘟病孢子显微图像识别方法

阅读:801发布:2022-09-04

专利汇可以提供一种基于图像处理模式识别技术的稻瘟病孢子显微图像识别方法专利检索,专利查询,专利分析的服务。并且一种基于 图像处理 和 模式识别 技术的稻瘟病孢子显微 图像识别 方法,本 发明 涉及稻瘟病孢子显微图像识别领域。本发明是要解决大田检测通常根据肉眼观察判别不准确、目前实验对样本数需求巨大、传统 显微镜 孢子计数需要大量时间以及增大对初期灾情的发现难度的问题,而提出的一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法。该方法是通过步骤一成原始图像转灰度图;步骤二将灰度图进行图像增强处理;步骤三得到二值化效果图;步骤四得到去噪效果图;步骤五得到疑似稻瘟病孢子的图形轮廓;步骤六识别出稻瘟病孢子显微图像,最终完成稻瘟病孢子的识别并统计稻瘟病孢子数量等步骤实现的。本发明应用于稻瘟病孢子显微图像领域。,下面是一种基于图像处理模式识别技术的稻瘟病孢子显微图像识别方法专利的具体信息内容。

1.一种基于图像处理模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法具体是按照以下步骤进行的:
步骤一、成像系统获取图像即原始图像转灰度图;
步骤二、将灰度图进行图像增强处理,得到的直方图均衡化效果图;
步骤三、将直方图均衡化效果图进行局部自适应阈值分割,得到二值化效果图;
步骤四、将二值化效果图进行形态学变换的去噪处理,得到去噪效果图;
步骤五、将去噪效果图通过Canny边缘检测,得到只含有边缘信息的前景图,将只含边缘信息的前景图进行轮廓提取得到轮廓的个数,并将各轮廓以轮廓树的形式保存起来;按照轮廓树逐一绘制出前景图的各轮廓,即得到疑似稻瘟病孢子的图形轮廓;
步骤六、采集提取出的疑似稻瘟病孢子的图形轮廓特征值,利用图形轮廓特征值对疑似稻瘟病孢子的图形轮廓进行分类统计出稻瘟病孢子的图形轮廓和杂质的图形轮廓;从而识别出稻瘟病孢子显微图像;最终完成稻瘟病孢子的识别并统计稻瘟病孢子数量;其中,采集提取出的疑似稻瘟病孢子的图形轮廓特征值包括形态特征和光密度特征两种;形态特征包括:宽度、高度、面积、周长、偏心率、圆形度和矩形度特征;光密度特征包括综合光密度、平均灰度、光密度方差以及特征灰度;即完成了一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法。
2.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤一中成像系统获取图像即原始图像转灰度图的具体过程为:
将三通道R、G和B图像转换为单通道的灰度图;即R、G和B三通道颜色数据转换为灰度数据Gray的换算公式如下:
Gray=(R+G+B)/3 (1)
其中,R为图像中的红色通道,G为图像中的绿色通道,B为图像中的蓝色通道。
3.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤二中将灰度图进行图像增强处理具体过程为:
采用空间域的图像增强方法;
(1)设变量r代表灰度图中像素灰度级;灰度图是根据式(1)由RGB空间转换而来的,r的值在下述范围之内:
0≤r≤255
在灰度级中,r=0代表黑,r=255代表白;r为灰度图的像素的强度是一个随机变量;
(2)利用直方图均衡化拉伸灰度图的像素的强度r分布来增强灰度图对比度
将一幅图像中灰度级rk出现的概率pr(rk)近似为:
其中,n是灰度图中像素的总数,nk是灰度级为rk的像素个数,l为灰度图中灰度级的总数;在直坐标系中作出rk与pr(rk)的关系图形,这个图形称为直方图;
灰度图的累积分布函数即输出图像灰度级sk为:
通过式(2),将输人灰度图中灰度级为rk的各像素映射到输出图像灰度级为sk的对应像素得到输出图像;得到的直方图均衡化效果图。
4.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤三中将直方图均衡化效果图进行局部自适应阈值分割,得到二值化效果图具体过程为:
根据直方图均衡化效果图灰度级像素的邻域的像素值分布来确定直方图均衡化效果图灰度级像素位置上的二值化阈值。
5.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤四中将二值化效果图进行形态学变换的去噪处理,得到去噪效果图的具体过程为:
(1)利用腐蚀方法将二值化效果图上的点与核进行卷积,并把卷积结果放到与二值化效果图上的点坐标相同的新的图像中的点上;从而计算核区域像素的最小值;并得到一个腐蚀后的图像;其中,核是一个3*3像素的中间带有参考点的实心圆盘;
(2)膨胀是腐蚀的反操作,即核与腐蚀后的图像卷积后,计算核覆盖的区域的像素点最大值,并把这个最大值赋值给与腐蚀后的图像上的点坐标相同的新的图像中的点上;则新的图像为去噪后效果图。
6.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤五中将去噪效果图通过Canny边缘检测,得到只含有边缘信息的前景图具体过程:
(1)设二维高斯函数G(x,y)
在某一方向n上的一阶方向导数为:其中x,y为去噪效果图横纵坐标,σ为方差;
(2)根据Canny边缘检测算子的定义,中心边缘点为算子Gn:
其中, 是方向矢量, 是梯度矢量与图像f(x,y)的卷积在边缘
梯度方向上的区域中的最大值即为边缘点;θ为方向矢量的方向角;
(3)Canny算法采用双阈值技术,应用边缘提取时,将去噪效果图通过高斯卷积Gn*f(x,y)进行平滑,得到平滑后的图像;
(4)将图像f(x,y)与Gn作卷积,同时改变n的方向,Gn*f(x,y)取得最大值时,求得正交于检测边缘的方向n:

因此,对应于极值的方向n有
(5)在方向n上Gn*f(x,y)有最大输出响应|Gn*I|,此时
二维最优阶跃边缘算子是以卷积Gn*f(x,y)为基础的,去噪效果图的边缘强度由决定,而去噪效果图的边缘方向为 确
定去噪效果图的每一个点的梯度方向强度为最大值的点为疑似边缘点;疑似边缘点的梯度数量图像呈现屋脊状;
(6)算子Gn沿着梯度数量图像中屋脊的最大值开始进行边缘跟踪,并将不在屋脊最大值的梯度数量图像像素设为0值,从而输出一条很细的梯度数量图像疑似边缘线,即非最大抑制;其中,边缘追踪的过程采用了滞后策略,由T1和T2两个阈值(T1>T2)控制,从屋脊的最大值大于T1的点开始追踪,随后沿着屋脊的两个方向继续进行追踪,直到屋脊的两个方向中去噪效果图的边缘强度小于T2停止;高阈值T1和低阈值T2是通过对图像直方图的统计确定的具体做法是:
首先去除去噪效果图像素为零的点,然后求取去噪效果图中剩余像素的均值μ和标准差σ,高阈值T1和低阈值T2的选取如下式:
(7)如果疑似边缘线中点的边缘强度大于高阈值T1,那么疑似边缘线中的点一定是边缘点;如果疑似边缘线中点的边缘强度低于低阈值,那么疑似边缘线中的点一定不是边缘点;如果疑似边缘线中点的边缘强度在低阈值和高阈值之间,判断疑似边缘线中的点中
3×3邻域中的疑似边缘线中的点的8个邻接像素点有没有大于高阈值T1的边缘点,如果有大于高阈值T1的边缘点则疑似边缘线中的点为边缘点,所有边缘点得到只含有边缘信息的前景图。
7.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤五中将只含边缘信息的前景图进行轮廓提取得到轮廓的个数,并将各轮廓以轮廓树的形式保存起来;按照轮廓树逐一绘制出前景图的各轮廓,即得到疑似稻瘟病孢子的图形轮廓的具体为:
提取边缘信息的前景图的轮廓信息;在OpenCV里提取目标轮廓的函数是
findContours,识别前景图的轮廓;其中findContours的输入图像是一幅二值图像即边缘信息的前景图,findContours输出的是边缘信息的前景图中每一个连通区域的轮廓点的集合;轮廓点的集合包括轮廓个数和各轮廓上点的个数信息,并可根据轮廓点的集合绘制出疑似稻瘟病孢子的图形轮廓。
8.根据权利要求1所述一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,其特征在于:步骤六中采集提取出的疑似稻瘟病孢子的图形轮廓特征值,利用图形轮廓特征值对疑似稻瘟病孢子的图形轮廓进行分类统计出稻瘟病孢子的图形轮廓和杂质的图形轮廓具体过程为:
1)采集提取出疑似稻瘟病孢子的图形轮廓特征值;
2)利用疑似稻瘟病孢子的图形轮廓特征值创建样本,利用样本训练决策树、贝叶斯、支持向量机和K近邻四种分类器模型;
3)用四个基于特征值的分类器即决策树、贝叶斯、支持向量机和K近邻四种分类器模型对疑似稻瘟病孢子的图形轮廓进行判别,检测出稻瘟病孢子;其中,对疑似稻瘟病孢子的图形轮廓进行判别具体过程为:
(1)当决策树、贝叶斯、支持向量机和K近邻四种分类器中至少有三个分类器对疑似稻瘟病孢子的图形轮廓进行判别的结果一致时,则一致的结果作为判别的结果即检测出稻瘟病孢子;
(2)当决策树、贝叶斯、支持向量机和K近邻四种分类器中对疑似稻瘟病孢子的图形轮廓进行判别结果两两一致时,采用与贝叶斯分类器结果不同的的那一个判别结果即检测出稻瘟病孢子。

说明书全文

一种基于图像处理模式识别技术的稻瘟病孢子显微图像

识别方法

技术领域

背景技术

[0002] 稻稻瘟病病害的早期检测和病害程度判别是对稻瘟病进行预测预报和化学防治基础和关键。由于病害初期阶段症状表现不明显,农业生产者往往缺乏作物病害诊断的专业知识,使得受害作物得不到有效控制,作物受害面积迅速扩大,病害加重。稻瘟病病害识别与分级检测,主要分为大田检测和实验室检测两方面进行。大田检测通常根据农学专家或者农业生产者的肉眼观察来完成,这种诊断方式常常会受到各种主观因素的影响,导致判别不准确。而实验室检测是指采集染病样本,在实验室环境下由专业的技术检测人员完成对稻瘟病孢子的识别与数目统计。由于样本数巨大,并且孢子个体很小,传统显微镜孢子计数需要大量时间,而长时间观察易造成眼睛疲劳。此外,如果样本中的病原孢子中混有其他种类的病原孢子,将对传统显微镜孢子计数造成很大的干扰,大大增加了对初期灾情的发现难度。
[0003] 图像识别是计算机应用领域中利用计算机对图像进行处理、分析和理解,以识别各种不同模式的目标和对象的技术。基于图像识别技术对稻瘟病孢子进行自动识别与统计具有快速、低成本、智能化等优点,但目前尚无其他针对稻瘟病孢子显微图像的计算机自动识别技术。

发明内容

[0004] 本发明的目的是为了解决大田检测通常根据肉眼观察判别不准确、目前实验对样本数需求巨大,并且孢子个体很小,传统显微镜孢子计数需要大量时间,而长时间观察易造成眼睛疲劳增加了对初期灾情的发现难度的问题,而提出的一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法。
[0005] 上述的发明目的是通过以下技术方案实现的:
[0006] 步骤一、成像系统获取图像即原始图像转灰度图;
[0007] 步骤二、将灰度图进行图像增强处理,得到的直方图均衡化效果图;
[0008] 步骤三、将直方图均衡化效果图进行局部自适应阈值分割,得到二值化效果图;
[0009] 步骤四、将二值化效果图进行形态学变换的去噪处理,得到去噪效果图;
[0010] 步骤五、将去噪效果图通过Canny边缘检测,得到只含有边缘信息的前景图,将只含边缘信息的前景图进行轮廓提取得到轮廓的个数,并将各轮廓以轮廓树的形式保存起来;按照轮廓树逐一绘制出前景图的各轮廓,即得到疑似稻瘟病孢子的图形轮廓;
[0011] 步骤六、采集提取出的疑似稻瘟病孢子的图形轮廓特征值,利用图形轮廓特征值对疑似稻瘟病孢子的图形轮廓进行分类统计出稻瘟病孢子的图形轮廓和杂质的图形轮廓;从而识别出稻瘟病孢子显微图像;最终完成稻瘟病孢子的识别并统计稻瘟病孢子数量;其中,采集提取出的疑似稻瘟病孢子的图形轮廓特征值包括形态特征和光密度特征两种;形态特征包括:宽度、高度、面积、周长、偏心率、圆形度和矩形度特征;光密度特征包括综合光密度、平均灰度、光密度方差以及特征灰度;即完成了一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法。
[0012] 发明效果
[0013] 本发明是一种智能化的稻瘟病孢子显微图像识别处理方法。其步骤为:图像预处理、图像分割、孢子特征提取、模式识别。优点:引入自适应阈值分割与自适应边缘检测算法,使轮廓提取的信息损失降到最小;分类器设计采用投票的方法,避免单个分类器分类精度较低的情况,提高整体的分类精度。当各分类器的分类结果差别较大时,采用对原图像的分类结果作为参考和降低贝叶斯分类器权值的方法进行判别,并返回提示信息,可进行人为更改分类结果。图像质量的影响通过以上研究,已完成了对稻瘟病孢子显微图像的预处理、边缘检测、轮廓提取工作,但考虑到实际拍摄图像的质量差异问题,在此对于质量较差的图像处理进行如下分析。
[0014] 如图7所示为一质量较差样本,图像中存在着大量的噪声即干扰信息,对该样本进行直方图均衡化处理后的效果如图8。
[0015] 可知,稻瘟病孢子的轮廓与背景的对比度得到了加强,但部分干扰色仍存在较强干扰,按照本设计的处理方法继续进行Canny边缘检测。
[0016] 如图9(a)~(d)可发现,不同的阈值对于检测的效果影响很大。对于图9(a)~(d),选择10的低阈值和350的高阈值比较合理,而选取其他阈值则将无法识别出所需信息。选取该阈值进行轮廓识别可以得到如图10(a)~(b)的效果。
[0017] 可见,阈值的选择对图像的分割提取至关重要,只有当阈值选择合适时才能完成对稻瘟病孢子的识别。而人为修改阈值费时费,丧失了计算机图像识别的优越性,因此对原算法进行了改良,采用自适应阈值的二值化与边缘检测算法。附图说明
[0018] 图1具体实施方式一提出的稻瘟病孢子图像识别过程框图
[0019] 图2具体实施方式二提出的稻瘟病孢子显微原图;
[0020] 图3具体实施方式二提出的灰度化效果图;
[0021] 图4具体实施方式一提出的直方图均衡化效果图;
[0022] 图5具体实施方式六提出的Canny边缘检测效果图;
[0023] 图6具体实施方式七提出的灰度图轮廓提取效果图;
[0024] 图7具体实施方式一提出的成像系统获取图像即原始图像;
[0025] 图8具体实施方式一提出的直方图均衡化效果图;
[0026] 图9(a)具体实施方式一提出的阈值为10100的Canny边缘检测效果图;
[0027] 图9(b)具体实施方式一提出的阈值为10350Canny边缘检测效果图;
[0028] 图9(c)具体实施方式一提出的阈值为80100Canny边缘检测效果图;
[0029] 图9(d)具体实施方式一提出的阈值为100350Canny边缘检测效果图;
[0030] 图10(a)具体实施方式一提出的轮廓提取效果图;
[0031] 图10(b)具体实施方式一提出的形状提取效果图;
[0032] 图11具体实施方式四提出的自适应二值化示意图;
[0033] 图12具体实施方式五提出的膨胀腐蚀去噪处理后的图像示意图;
[0034] 图13具体实施方式六提出的自适应Canny边缘检测示意图;
[0035] 图14具体实施方式七提出的原始图像轮廓提取示意图;
[0036] 图15具体实施方式三提出的原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布的示意图;
[0037] 图16具体实施方式二提出的RGB模型示意图。

具体实施方式

[0038] 具体实施方式一:本实施方式的一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法,具体是按照以下步骤制备的:
[0039] 步骤一、成像系统获取图像即原始图像转灰度图;
[0040] 步骤二、将灰度图进行图像增强处理,得到的直方图均衡化效果图如图4;
[0041] 步骤三、将直方图均衡化效果图进行局部自适应阈值分割,得到二值化效果图;
[0042] 步骤四、将二值化效果图进行形态学变换的去噪处理,得到去噪效果图;
[0043] 步骤五、将去噪效果图通过Canny边缘检测,得到只含有边缘信息的前景图,将只含边缘信息的前景图进行轮廓提取得到轮廓的个数,并将各轮廓以轮廓树的形式保存起来;按照轮廓树逐一绘制出前景图的各轮廓,即得到疑似稻瘟病孢子的图形轮廓;
[0044] 步骤六、采集提取出的疑似稻瘟病孢子的图形轮廓特征值,利用图形轮廓特征值对疑似稻瘟病孢子的图形轮廓进行分类统计出稻瘟病孢子的图形轮廓和杂质的图形轮廓;从而识别出稻瘟病孢子显微图像;最终完成稻瘟病孢子的识别并统计稻瘟病孢子数量;其中,采集提取出的疑似稻瘟病孢子的图形轮廓特征值包括形态特征和光密度特征两种;形态特征包括:宽度、高度、面积、周长、偏心率、圆形度和矩形度特征;光密度特征包括综合光密度、平均灰度、光密度方差以及特征灰度;如图1即完成了一种基于图像处理和模式识别技术的稻瘟病孢子显微图像识别方法。
[0045] 本实施方式效果:
[0046] 本实施方式是一种智能化的稻瘟病孢子显微图像识别处理方法。其步骤为:图像预处理、图像分割、孢子特征提取、模式识别。优点:引入自适应阈值分割与自适应边缘检测算法,使轮廓提取的信息损失降到最小;分类器设计采用投票的方法,避免单个分类器分类精度较低的情况,提高整体的分类精度。当各分类器的分类结果差别较大时,采用对原图像的分类结果作为参考和降低贝叶斯分类器权值的方法进行判别,并返回提示信息,可进行人为更改分类结果。图像质量的影响通过以上研究,已完成了对稻瘟病孢子显微图像的预处理、边缘检测、轮廓提取工作,但考虑到实际拍摄图像的质量差异问题,在此对于质量较差的图像处理进行如下分析。
[0047] 如图7所示为一质量较差样本,图像中存在着大量的噪声即干扰信息,对该样本进行直方图均衡化处理后的效果如图8。
[0048] 可知,稻瘟病孢子的轮廓与背景的对比度得到了加强,但部分干扰色块仍存在较强干扰,按照本设计的处理方法继续进行Canny边缘检测。
[0049] 如图9(a)~(d)可发现,不同的阈值对于检测的效果影响很大。对于图9(a)~(d),选择10的低阈值和350的高阈值比较合理,而选取其他阈值则将无法识别出所需信息。选取该阈值进行轮廓识别可以得到如图10(a)~(b)的效果。
[0050] 可见,阈值的选择对图像的分割提取至关重要,只有当阈值选择合适时才能完成对稻瘟病孢子的识别。而人为修改阈值费时费力,丧失了计算机图像识别的优越性,因此对原算法进行了改良,采用自适应阈值的二值化与边缘检测算法。
[0051] 具体实施方式二:本实施方式与具体实施方式一不同的是:步骤一中成像系统获取图像即原始图像转灰度图的具体过程为:
[0052] 由于显微图像的色彩识别度不高,对孢子的识别主要依据光强的明暗;而单通道的图像数据更有利于后续的图像处理,可以缩短处理时间;
[0053] 在RGB模型如图16中,每种颜色出现在红、绿和蓝的原色光谱分量中,这个模型基于笛卡儿坐标系统,所考虑的彩色子空间是图16所示的立方体;图中,R,G,B(红、绿和蓝)位于立方体的3个上;青、深红和黄位于另外3个角上,黑色在原点处,白色位于离原点最远的角上;在该模型中,灰度等级沿着这两点的连线分布;在RGB模型中,不同的颜色处在立方体上或其内部,并用从原点分布的向量来定义,每个通道的取值范围为(0,255);
[0054] 将三通道R、G和B图像转换为单通道的灰度图;即R、G和B三通道颜色数据转换为灰度数据Gray的换算公式如下:
[0055] Gray=(R+G+B)/3 (1)
[0056] 其中,R为图像中的红色通道,G为图像中的绿色通道,B为图像中的蓝色通道;如图2和图3可以看到灰度图较好地保留了显微原始图像的识别信息。其它步骤及参数与具体实施方式一相同。
[0057] 具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤二中将灰度图进行图像增强处理具体过程为:
[0058] 一般情况下,成像系统获取图像即原始图像由于受到种种条件限制和随机干扰,往往不能在视觉系统中直接使用,必须利用视觉信息处理的早期阶段对原始图像进行灰度校正、噪声过滤等图像预处理.对机器视觉系统来说,所用的图像预处理方法并不考虑图像降质原因,只将图像中感兴趣的特征有选择地突出,衰减其不需要的特征,故预处理后的输出图像并不需要去逼近原图像.这类图像预处理方法统称为图像增强.图像增强技术主要有两种方法:空间域法和频率域法.空间域方法主要是在空间域内对图像像素直接运算处理.本发明采用空间域的图像增强方法;
[0059] (1)图像增强是根据一定的要求突出图像中的某些信息,去除或削弱某些不需要的信息的方法,其目的是突出图像中的“有用”信息,扩大图像中不同物体特征之间的差别;本发明选用直方图均衡化的方法进行图像增强;灰度级的直方图就是反映一幅图像中的灰度级与出现这种灰度的概率之间的关系的图形;
[0060] 设变量r代表灰度图中像素灰度级;灰度图是根据式(1)由RGB空间转换而来的,r的值在下述范围之内:
[0061] 0≤r≤255
[0062] 在灰度级中,r=0代表黑,r=255代表白;对于灰度图来说,每一个灰度图像素取得[0,255]区间的灰度级是随机的,也就是说r为灰度图的像素的强度是一个随机变量;
[0063] (2)利用直方图均衡化拉伸灰度图的像素的强度r分布来增强灰度图对比度;通过增强灰度图对比度,亮度可以更好地在直方图上分布;直方图均衡化处理的“中心思想”是把原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布如图15;
[0064] 由于我们要处理的是数字图像,其像素为离散值,将一幅图像中灰度级rk出现的概率pr(rk)近似为:
[0065]
[0066] 其中,n是灰度图中像素的总数,nk是灰度级(灰度图像素的强度)为rk的像素个数,l为灰度图中灰度级的总数;在直角坐标系中作出rk与pr(rk)的关系图形,这个图形称为直方图;
[0067] 灰度图的累积分布函数即输出图像灰度级sk为:
[0068]
[0069] 比较概率分布函数和累积分布函数,概率分布函数的二维图像是参差不齐的,累积分布函数是单调递增的;因此,通过式(2),将输人灰度图中灰度级为rk的各像素映射到输出图像灰度级为sk的对应像素得到输出图像,式(2)给出的变换(映射)称做直方图均衡化或直方图线性化(空间域图像增强方法);得到的直方图均衡化效果图如图4所示;图4为直方图均衡化效果图,可以看到,灰度图的对比度得到了明显增强,稻瘟病孢子的特征变得更加明显。其它步骤及参数与具体实施方式一或二相同。
[0070] 具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:步骤三中将直方图均衡化效果图进行局部自适应阈值分割,得到二值化效果图具体过程为:
[0071] 根据直方图均衡化效果图灰度级像素的邻域块的像素值分布来确定直方图均衡化效果图灰度级像素位置上的二值化阈值;这样做的好处在于每个像素位置处的二值化阈值不是固定不变的,而是由其周围邻域像素的分布来决定的;不同亮度、对比度和纹理的局部图像区域拥有相对应的局部二值化阈值如图11。其它步骤及参数与具体实施方式一至三之一相同。
[0072] 具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:步骤四中将二值化效果图进行形态学变换的去噪处理,得到去噪效果图的具体过程为:
[0073] 对图像进行形态学变换,采用先腐蚀后膨胀算法;
[0074] (1)利用腐蚀方法将二值化效果图上的点与核进行卷积,并把卷积结果放到与二值化效果图上的点坐标相同的新的图像中的点上;从而计算核区域像素的最小值;并得到一个腐蚀后的图像;其中,核可以是任何的形状或大小,它拥有一个单独定义出来的参考点(anchor point);本发明中,核是一个3*3像素的中间带有参考点的实心圆盘;
[0075] (2)膨胀是腐蚀的反操作,即核与腐蚀后的图像卷积后,计算核覆盖的区域的像素点最大值,并把这个最大值赋值给与腐蚀后的图像上的点坐标相同的新的图像中的点上;这样就会使图像中的高亮区域逐渐增长,这样的增长就是“膨胀操作”;
[0076] 先腐蚀后膨胀的算法称为开运算,其效果是消除了高于其邻近点的孤立点,即达到了去噪的目的;则新的图像为去噪后效果图如图12。其它步骤及参数与具体实施方式一至四之一相同。
[0077] 具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:步骤五中将去噪效果图通过Canny边缘检测,得到只含有边缘信息的前景图具体过程:
[0078] 图像的边缘是指图像局部区域亮度变化显著的部分,该区域的灰度剖面一般可以看作是一个阶跃,既从一个灰度值在很小的缓冲区域内急剧变化到另一个灰度相差较大的灰度值;图像边缘是图像分割所依赖的重要特征,通过Canny边缘检测算法完成去噪效果图的灰度变化的度量、灰度变化的检测和灰度变化的定位工作:
[0079] (1)根据对信噪比与定位乘积进行测度,得到最优化逼近算子即Canny边缘检测算子;对于阶跃的边缘,Canny推出的最优边缘检测器的形状与高斯函数的一阶导数类似,而根据二维高斯函数的圆对称性和可分解性,可以很容易地计算高斯函数在任意方向上的方向导数与图像的卷积;设二维高斯函数G(x,y)
[0080]
[0081] 在某一方向n上的一阶方向导数为:其中x,y为去噪效果图横纵坐标,σ为方差;
[0082] (2)根据Canny边缘检测算子的定义,中心边缘点为算子Gn:
[0083]
[0084] 其中, 是方向矢量, 是梯度矢量与图像f(x,y)的卷积在边缘梯度方向上的区域中的最大值即为边缘点;这样,就可以在每一点的梯度方向上判断此点强度是否为其最大值来确定该点是否为边缘点;θ为方向矢量的方向角;
[0085] (3)Canny算法采用双阈值技术,应用边缘提取时,将去噪效果图通过高斯卷积Gn*f(x,y)进行平滑,得到平滑后的图像;
[0086] (4)将图像f(x,y)与Gn作卷积,同时改变n的方向,Gn*f(x,y)取得最大值时,求得正交于检测边缘的方向n:
[0087]
[0088] 得
[0089]
[0090] 因此,对应于极值的方向n有
[0091]
[0092] (5)在方向n上Gn*f(x,y)有最大输出响应|Gn*I|,此时
[0093]
[0094] 二维最优阶跃边缘算子是以卷积Gn*f(x,y)为基础的,去噪效果图的边缘强度由决定,而去噪效果图的边缘方向为 确定去噪效果图的每一个点的梯度方向强度为最大值的点为疑似边缘点;疑似边缘点的梯度数量图像呈现屋脊状;
[0095] (6)算子Gn沿着梯度数量图像中屋脊的最大值开始进行边缘跟踪,并将不在屋脊最大值的梯度数量图像像素设为0值,从而输出一条很细的梯度数量图像疑似边缘线,即非最大抑制;其中,边缘追踪的过程采用了滞后策略,由T1和T2两个阈值(T1>T2)控制,从屋脊的最大值大于T1的点开始追踪,随后沿着屋脊的两个方向继续进行追踪,直到屋脊的两个方向中去噪效果图的边缘强度小于T2停止;高阈值T1和低阈值T2是通过对图像直方图的统计确定的具体做法是:
[0096] 首先去除去噪效果图像素为零的点,然后求取去噪效果图中剩余像素的均值μ和标准差σ,高阈值T1和低阈值T2的选取如下式如图13:
[0097]
[0098] (7)如果疑似边缘线中点的边缘强度大于高阈值T1,那么疑似边缘线中的点一定是边缘点;如果疑似边缘线中点的边缘强度低于低阈值,那么疑似边缘线中的点一定不是边缘点;如果疑似边缘线中点的边缘强度在低阈值和高阈值之间,判断疑似边缘线中的点中3×3邻域中的疑似边缘线中的点的8个邻接像素点有没有大于高阈值T1的边缘点,如果有大于高阈值T1的边缘点则疑似边缘线中的点为边缘点,所有边缘点得到只含有边缘信息的前景图如图5为Canny边缘检测效果图,可以看到,轮廓清晰的孢子信息被保留了下来,而多余的背景信息被过滤掉了。其它步骤及参数与具体实施方式一至五之一相同。
[0099] 具体实施方式七:本实施方式与具体实施方式一至六之一不同的是:步骤五中将只含边缘信息的前景图进行轮廓提取得到轮廓的个数,并将各轮廓以轮廓树的形式保存起来;按照轮廓树逐一绘制出前景图的各轮廓,即得到疑似稻瘟病孢子的图形轮廓的具体为:
[0100] 提取边缘信息的前景图的轮廓信息;在OpenCV(Open Source Computer Vision Library是一个基于(开源)发行的跨平台计算机视觉库)里提取目标轮廓的函数是findContours,识别前景图的轮廓;其中findContours的输入图像是一幅二值图像即边缘信息的前景图,findContours输出的是边缘信息的前景图中每一个连通区域的轮廓点的集合;轮廓点的集合包括轮廓个数和各轮廓上点的个数等信息,并可根据轮廓点的集合绘制出疑似稻瘟病孢子的图形轮廓;如图6和图14为轮廓提取效果图。其它步骤及参数与具体实施方式一至六之一相同。
[0101] 具体实施方式八:本实施方式与具体实施方式一至七之一不同的是:步骤六中采集提取出的疑似稻瘟病孢子的图形轮廓特征值,利用图形轮廓特征值对疑似稻瘟病孢子的图形轮廓进行分类统计出稻瘟病孢子的图形轮廓和杂质的图形轮廓具体过程为:
[0102] 1)采集提取出疑似稻瘟病孢子的图形轮廓特征值;
[0103] 2)利用疑似稻瘟病孢子的图形轮廓特征值创建样本,利用样本训练决策树、贝叶斯、支持向量机和K近邻四种分类器模型;
[0104] 3)用四个基于特征值的分类器即决策树、贝叶斯、支持向量机和K近邻四种分类器模型对疑似稻瘟病孢子的图形轮廓进行判别,检测出稻瘟病孢子;其中,对疑似稻瘟病孢子的图形轮廓进行判别具体过程为:
[0105] (1)当决策树、贝叶斯、支持向量机和K近邻四种分类器中至少有三个分类器对疑似稻瘟病孢子的图形轮廓进行判别的结果一致时,则一致的结果作为判别的结果即检测出稻瘟病孢子;
[0106] (2)当决策树、贝叶斯、支持向量机和K近邻四种分类器中对疑似稻瘟病孢子的图形轮廓进行判别结果两两一致时,采用与贝叶斯分类器结果不同的的那一个判别结果即检测出稻瘟病孢子。其它步骤及参数与具体实施方式一至七之一相同。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈