CN108052886A - 一种小麦条锈病菌夏孢子自动统计计数方法 - Google Patents
一种小麦条锈病菌夏孢子自动统计计数方法 Download PDFInfo
- Publication number
- CN108052886A CN108052886A CN201711270040.8A CN201711270040A CN108052886A CN 108052886 A CN108052886 A CN 108052886A CN 201711270040 A CN201711270040 A CN 201711270040A CN 108052886 A CN108052886 A CN 108052886A
- Authority
- CN
- China
- Prior art keywords
- spore
- spores
- contour
- summer
- counting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
- G06V20/695—Preprocessing, e.g. image segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
Abstract
本发明提供了一种小麦条锈病菌夏孢子自动统计计数方法，其可用于野外田间空气中小麦条锈病菌夏孢子的定量监测和准确计数。该方法包括以下步骤：通过TPBZ3型孢子捕捉仪、研究级显微镜将小麦条锈病菌夏孢子的显微图像采集到计算机；然后进行K‑means聚类算法、形态学处理、基于形状因子和面积的粘连孢子判别、基于凹度的粘连轮廓分割和粘连孢子轮廓段融合与计数，其计数准确率可达98.6％，实现了夏孢子的自动统计与计数。本发明的方法计数精度高，实用性强，解决了野外田间空气中小麦条锈病菌夏孢子计数这一亟待解决的难题。为在线式小麦条锈病夏孢子监测装备的开发提供技术支持。
Description
技术领域
本发明属于农业图像处理技术领域，特别涉及一种小麦条锈病菌夏孢子自动统计计数方法。
背景技术
由条形柄锈菌(Puccinia striiformis f.sp.tritici)引起的小麦条锈病是我国西北、西南、华北和淮北等冬麦区和西北春麦区小麦的重要病害，一般流行年份可致小麦减产10％～20％，特大流行年份可致小麦减产50％～60％，对我国粮食安全生产具有重大威胁。小麦条锈病是一种气传病害，病原菌主要以夏孢子随气流传播到小麦上完成周年侵染循环，因此，田间空气中的夏孢子数量与小麦条锈病的流行程度紧密相关。对田间空气中条锈病菌夏孢子的定量监测和准确计数，并利用计数结果进行病害的早期预测，是及早制定正确的防控策略和采取防治措施的一个重要依据。
目前，小麦条锈病菌夏孢子数量多采用孢子捕捉器取样和监测，然后将黏附有孢子的载玻片拿回实验室在显微镜下人工计数或分子生物学方法计数。传统的人工计数方法是在光学显微镜下通过肉眼观测载玻片确定孢子个数，或者是将黏着的孢子洗刷下来，配成悬浮液，制成临时玻片进行镜检计数，均存在工作量大、效率低、且随工作时间准确性降低等缺点。分子生物学方法需要根据病原孢子DNA的量实现计数，过程复杂，对技术、试验条件和仪器要求高、花费高，很难大面积推广应用。因此，急需一种简便、快捷、准确、高效的小麦条锈病菌夏孢子计数方法。
近年来，基于显微图像处理技术被逐步应用到气传植物病原真菌孢子检测的相关研究上，并逐步应用到小麦条锈病菌夏孢子的相关研究中。齐龙等提出了基于距离变换和高斯滤波的改进分水岭算法的稻瘟病菌显微图像孢子自动检测和计数方法，但当粘连孢子的接触线长度超过单个孢子的宽度时，容易造成局部极小值点间的距离小于结构元素的长度，引起漏分割现象。李小龙等基于K-means聚类和分水岭分割算法，对小麦条锈病菌夏孢子显微图像进行处理，实现了对夏孢子的自动计数。然而图像中的噪声、结构特征和孢子表面细微的灰度变化等因素，易导致出现极小值点过多、分割位置不准确，从而产生过度分割现象。张荣标等提出了图像处理和支持向量机的圆褐固氮菌浓度快速检测方法，但方法仅适用于圆褐固氮菌相互粘连较少的简单情况。此外，上述方法中，大多数方法对未粘连的孢子能很好的计数，但对经常发生的多孢子粘连分割计数问题尚未解决。
发明内容
为了克服上述现有技术的缺点，本发明的目的在于提供一种小麦条锈病菌夏孢子自动统计计数方法，利用一系列的显微图像处理方法实现夏孢子的自动统计与计数，以期提高小麦条锈病菌夏孢子计数精度，解决多孢子粘连分割计数的难题，并为在线式小麦条锈病夏孢子监测装备的开发提供技术支持。
为了实现上述目的，本发明采用的技术方案是：
一种小麦条锈病菌夏孢子自动统计计数方法，其特征在于包括以下步骤：
步骤1，夏孢子显微图像采集
用TPBZ3型孢子捕捉仪模拟捕捉田间空气中的小麦条锈病菌夏孢子。在野外小麦田间，抽出捕捉仪的载波器并将表面均匀涂抹一薄层凡士林的载玻片放其上(涂面朝上)，捕捉一定时间后取出载玻片。为了使玻片上黏附的孢子密度不同，按照上述方法制备凡士林玻片30片，获取捕捉时间分别为120、180、240min的条锈病菌夏孢子载玻片，重复10次，共获不同孢子密度的载玻片30片；
用BX52型倒置显微镜对孢子临时玻片进行观察、拍照，放大倍数为10×20。在显微镜下，分别对30片载玻片进行显微图像采集，每片随机选取5个视野拍照，共获得150幅供试夏孢子图像，图像尺寸为4140×3096Pixel，图像分辨率为72dpi，存贮为bmp格式；
步骤2，显微图像选取
由于显微图像尺寸过大，为减少图像处理运算量，提高检测和计数速度，将所述步骤S1的原始夏孢子显微图像用双线性插值方法缩小为911×682Pixel。然后随机选取30幅用于算法训练，其余120幅图像用于测试；
S3、K-means聚类分割
将所述步骤S2的夏孢子显微图像用K-means聚类算法实现孢子和背景的分割，即将其分为孢子目标和背景2大类，然后采用全局阈值算法进行二值化处理得到二值图像；
S4、形态学处理
二值图像中有小面积区域等噪声，且夏孢子目标边缘有凸刺或内部含有孔洞，故首先用填充操作对所述步骤S3中的二值图像进行孔洞填充，再用完整孢子中最小面积值为阈值，去除小面积区域噪声，此时当边界上孢子所占面积小于完整孢子最小面积的均从图像中移除，计数时作为其他视野中的孢子。最后，用圆盘结构元素对二值图像进行形态学开运算，以消除夏孢子边界小的凸刺；
S5、基于形状因子和面积的粘连孢子判别
所述步骤S4中的二值图像中既有单个孢子，也有2个或多个孢子粘连在一起的孢子群。通过观察可知，粘连孢子群的区域轮廓要比单个孢子区域的轮廓复杂，故选择描述目标边界复杂程度的形状因子作为粘连孢子判别的依据。形状因子公式为：
SF＝4πS/L2
式中SF为形状因子；S为一个连通区域的面积像素值；L为连通区域的周长像素值；
经粘连孢子判别后，对判定为单个孢子的区域直接采用最小二乘椭圆拟合算法进行拟合并记录椭圆个数Num供后续统计计数。对粘连孢子区域将每个孢子一一分割出来并自动计数是本发明下面拟解决的问题；
S6、基于凹度的粘连孢子轮廓分割
对所述步骤S5二值图像中的粘连孢子区域通过移除内部像素点操作提取出边缘轮廓，然后遍历图像跟踪边缘轮廓点，并将轮廓点坐标保存在一个有序表中。pt(xt,yt)是边缘轮廓上任意一点，两向量ptpt-k和ptpt+k之间的夹角称为pt的凹度，pt-k和pt+k表示点pt的相邻轮廓点。凹度concavity的计算公式为：
其中凹点是满足如下两个条件的轮廓点：
(1)凹度concavity(pt)在角angle(δ1,δ2)范围之内；
(2)直线不包含在粘连孢子内；
粘连孢子边缘轮廓由轮廓上的所有凹点和由凹点分割成的多个轮廓段组成，如下式：
C＝CS1+...+CSi+...+CSN+cp1+...+cpj+...+cpK
式中C为粘连孢子轮廓；CSi为第i个轮廓段；N为轮廓段总数；cpj为第j个凹点；K为凹点总数；
删除凹点后，将粘连孢子边缘轮廓分割成多个轮廓段，通过连通区域标记，将轮廓段保存在一个有序的cell结构轮廓段集中，用于后续孢子轮廓段融合处理；
S7、粘连孢子轮廓段融合
经所述步骤S6的轮廓分割后同一孢子的轮廓可能被分割成多个轮廓段，所以，识别出同一孢子的轮廓段并融合，以获得各个孢子尽可能完整的椭圆拟合数据是本发明的关键。故本发明以轮廓段距离测量方法为准则确定候选轮廓段，对候选轮廓段进行最小二乘椭圆拟合并用偏移误差方法进行评价，将符合条件的轮廓段融合成新的轮廓段，对新的轮廓段进行椭圆拟合即为最后的正确椭圆；
S8、最后统计所有椭圆的总数即为夏孢子的个数。
所述步骤S4中的圆盘结构元素，大小为3×3。
所述步骤S5中的基于形状因子和面积的粘连孢子判别，形状因子阈值SF0＝0.8，面积阈值Smax＝560，依次提取单个连通区域i，若i满足公式：
SFi>SF0&Si<Smax
则判定区域i为未粘连的单个夏孢子，否则为粘连孢子。
所述步骤S6中的凹度的角δ1和δ2分别设定为50°和150°。
所述步骤S7的轮廓段距离方法，第i和第j个轮廓段之间的距离DM定义为公式：
式中d(pi1,pj1)、分别为轮廓段CSi和CSj的2端点位置的欧氏距离；为轮廓段CSi和CSj的中点位置的欧式距离。
所述步骤S7的最小二乘椭圆拟合算法的具体步骤为：
在二维平面坐标系中，一般的椭圆曲线方程可以用2个向量相乘的隐式方程来表示：
F(α,X)＝X·α＝ax2+bxy+cy2+dx+ey+f＝0
式中，α＝[a b c d e f]T为椭圆方程的系数；x和y分别是曲线上点横、纵坐标，X＝[x2 xy y2 x y 1]；
加约束条件：αTCα＝1，以保证拟合结果为椭圆，
其中C为一个6×6的矩阵，
方程F(α,Xi)为二值图像中边缘点(xi,yi)到给定椭圆方程的代数距离。根据最小二乘原理，椭圆拟合问题通过将代数距离平方之和
最小化来解决。方程改写为向量形式：E＝||Dα||2
其中D为一个m×6的矩阵，
引入拉格朗日系数并微分得：
式中λ为特征值；W为一个6×6的散射矩阵，W＝DTD；
根据广义特征值求解方法可得到：
式中λi和ui分别为特征值和特征向量。由此解得至多6个实数解(λi,αi)；
推导出数据点到椭圆的代数距离平方和公式为：
E＝||Dα||2＝αTDTDα＝αTWα＝λαTCα＝λ
由上式可知所需要的是最小的正特征值λi所对应的特征向量αi。得到解αi之后，即可实现对椭圆的拟合。
所述步骤S7的偏移误差公式为：
DEM(CS#,CE)＝E/M#＝λM#
式中E为由数据点到椭圆的代数距离平方和公式求出的给定轮廓段CS#到椭圆CE的最小二乘代数距离；M#为轮廓段CS#上点的总数。
若偏移误差小于阈值，则该轮廓段属于同一孢子的轮廓段，否则不是。
所述步骤S7的轮廓段融合步骤具体为：
(1)从轮廓段集中选取轮廓段中最长轮廓段CS1；
(2)将孢子轮廓段中剩余的轮廓段CSi分别与CS1作距离测量，设定距离测量阈值ωDM，若满足DM(CS1,CSi)<ωDM条件，则有
CS*＝{CSi|DM(CS1,CSi)<ωDM,i＝1,2,...,k}；
(3)将CS1依次与CS*中每一个轮廓段CSi作椭圆拟合，同时求出相应的候选椭圆和偏移误差DEM；
(4)若偏移误差DEM小于预设的阈值σDEM，则融合成新的轮廓段CSnew：
CSnew＝{CS1+CSi|DEM(CS1,CSi)<σDEM,i＝1,2,...,k}；
(5)对CSnew作椭圆拟合，作为一个孢子的正确拟合椭圆，孢子数量计数Num+1，并且删除轮廓段集中的CSnew；
(6)最后，判断轮廓段集中的轮廓段是否全部删除。若是，则输出所有正确拟合椭圆并且结束程序，否则返回到步骤1。
所述步骤S7中的距离测量阈值ωDM和偏移误差阈值σDEM由30幅夏孢子图像样本的训练确定，分别为40和95。
与现有技术相比，本发明通过K-means聚类分割及形态学处理，能够准确地将夏孢子从背景中分割出来，且保留了夏孢子原有形态。提出基于夏孢子形状因子和面积特征的粘连孢子群判别方法，可正确地判别出未粘连的单个孢子和粘连孢子。
同时，本发明小麦条锈病菌夏孢子自动统计计数方法，基于凹度的轮廓分割和轮廓段融合方法对粘连孢子进行处理，可准确识别出粘连孢子中同一孢子的轮廓段，通过最小二乘椭圆拟合算法对夏孢子进行拟合，可有效对粘连孢子进行分割并计数，提高了孢子计数结果的准确性。经过试验结果表明，本发明最低计数准确率为92.7％，最高计数准确率为100％，总平均计数准确率为98.6％，具有较高的孢子计数精度。为在线式小麦条锈病夏孢子监测装备的开发提供技术支持。
附图说明
图1为本发明的自动统计计数方法处理流程图。
图2为条锈病菌夏孢子原始图像。
图3为K-means聚类后的夏孢子图像。
图4为夏孢子二值图像。
图5为形态学处理后的夏孢子二值图像。
图6为单个孢子二值图像。
图7为粘连孢子二值图像。
图8为本发明粘连夏孢子椭圆拟合和计数过程实例：a.预处理后的二值图像；b.粘连孢子的轮廓；c.轮廓上十字状黑色凹点；d.轮廓段；e.候选拟合椭圆；f.孢子正确拟合椭圆；g.最终分割结果。
图9为本发明分割及计数效果图。
具体实施方式
下面结合附图并通过实施例对本发明作进一步的详细说明，以下实施例是对本发明的解释而本发明并不局限于此。
本发明的小麦条锈病菌夏孢子自动统计计数方法处理流程如图1所示，具体包含以下步骤：
1、供试材料与设备
供试小麦条锈病菌夏孢子在西北农林科技大学植物保护学院小麦抗病种质资源创制与利用研究中心的东南窑低温温室内繁育，其他材料和仪器有凡士林、载玻片、洗耳球、便携式定量风流孢子捕捉仪(型号：TPBZ3，浙江托普云农公司)、研究级倒置显微镜(型号：BX52，Olympus，Japan)、数码成像系统(型号：DP72，Olympus，Japan)、CellSensstandard显微图像采集软件(Olympus，Japan)和台式计算机。
2、夏孢子显微图像采集
用TPBZ3型孢子捕捉仪模拟捕捉田间空气中的小麦条锈病菌夏孢子。在野外小麦田间，抽出捕捉仪的载波器并将表面均匀涂抹一薄层凡士林的载玻片放其上(涂面朝上)，捕捉一定时间后取出载玻片。为了使玻片上黏附的孢子密度不同，按照上述方法制备凡士林玻片30片，获取捕捉时间分别为60、120、180、240min的条锈病菌夏孢子载玻片，重复10次，共获不同孢子密度的载玻片30片。
用BX52型倒置显微镜对孢子临时玻片进行观察、拍照，放大倍数为10×20。在显微镜下，分别对30片载玻片进行显微图像采集，每片随机选取5个视野拍照，共获得150幅供试夏孢子图像，图像尺寸为4140×3096Pixel，图像分辨率为72dpi，存贮为bmp格式。
3、显微图像选取
上述二个步骤之后，由于显微图像尺寸过大，为减少图像处理运算量，提高检测和计数速度，将原始图像用双线性插值方法缩小为911×682Pixel，然后随机选取30幅用于算法训练，其余120幅图像用于测试。
4K-means聚类分割
为了对夏孢子进行计数，首先需要将夏孢子从背景中分割出来。夏孢子目标分割可看成聚类问题，即图像中像素点类别未知的前提下，根据像素点的特征值，将图像划分为若干个区域。考虑到K-means聚类算法是一种准确、高效的目标分割算法，故本发明用K-means聚类算法实现对孢子图像的分割，即将其分为孢子目标和背景2大类，图2所示的夏孢子图像经K-means聚类后的分割结果如图3所示，可见，该算法能够较准确地将夏孢子分割出来。对图3采用全局阈值算法进行二值化处理后的结果如图4所示。
5形态学处理
由图4可知，二值图像中有小面积区域等噪声，且夏孢子目标边缘有凸刺或内部含有孔洞，故首先用填充操作对图像进行孔洞填充，再用完整孢子中最小面积值为阈值，去除小面积区域噪声，此时当边界上孢子所占面积小于完整孢子最小面积的均从图像中移除，计数时作为其他视野中的孢子。最后，为消除夏孢子边界小的凸刺，经试验用3×3的圆盘结构元素对二值图像进行形态学开运算，处理结果如图5所示，可见，经上述处理后，更好地保留了图像中孢子目标的原有形态，有利于进行后续处理。
6基于形状因子和面积的粘连孢子判别
从图5可看出，二值图像中既有单个孢子，也有2个或多个孢子粘连在一起的孢子。通过观察可知，粘连孢子的区域轮廓要比单个孢子区域的轮廓复杂，故选择描述目标边界复杂程度的形状因子和面积作为粘连孢子判别的依据。形状因子公式为：
SF＝4πS/L2 (1)
式中SF为形状因子；S为一个连通区域的面积像素值；L为连通区域的周长像素值。
当多个孢子相互粘连时，粘连孢子边界会出现凹陷而变得复杂，形状因子会相应变小，故SF对孢子是否粘连有很好的区分度。对30幅夏孢子图像中单个孢子和粘连孢子SF和S的统计可知，单个孢子SF范围为0.9080～1.0912，S为200～523Pixel；粘连孢子的SF在0.2625～0.7606，S为600～2301Pixel。故设置形状因子阈值SF0＝0.8、面积阈值Smax＝560，作为判定单个孢子和粘连孢子的依据。依次提取单个连通区域i，若i满足
SFi>SF0&Si<Smax (2)
则判定区域i为未粘连的单个夏孢子，否则为粘连孢子。对图5进行粘连孢子判别后，单个夏孢子的二值图像如图6所示，对判定为单个孢子的区域直接采用椭圆拟合算法进行拟合并记录椭圆个数Num供后续统计计数。粘连夏孢子的二值图像如图7所示。可见，单个夏孢子和粘连孢子被正确判别出来。
7基于凹度的粘连孢子轮廓分割
由图7可知，图像中孢子目标粘连或部分重叠，故轮廓往往是多个目标轮廓的混合集，由于孢子的外形似椭圆，椭圆边界上所有点的凹度变化是连续的，当椭圆间出现粘连或重叠时凹度会变小，基于此思想，本发明依次提取粘连孢子边界轮廓，并基于凹度提取出轮廓上的凹点，由凹点将轮廓分割成多个轮廓段。
图8为本发明粘连夏孢子椭圆拟合和计数过程，图7中矩形框内的粘连孢子放大后如图8a所示。对二值图像中的粘连孢子通过移除内部像素点操作提取出边缘轮廓(图8b)，然后遍历图像跟踪边缘轮廓点，并将轮廓点坐标保存在一个有序表中。pt(xt,yt)是边缘轮廓上任意一点，两向量ptpt-k和ptpt+k之间的夹角称为pt的凹度，pt-k和pt+k表示点pt的相邻轮廓点，根据预试验本文设定k值为3。凹度concavity的计算公式为：
其中凹点是满足如下两个条件的轮廓点：
(1)凹度concavity(pt)在角angle(δ1,δ2)范围之内；
(2)直线不包含在粘连孢子内；
本文根据30幅夏孢子图像样本的训练统计得出凹度的角δ1和δ2分别设定为50°和150°。
粘连孢子边缘轮廓由轮廓上的所有凹点和由凹点分割成的多个轮廓段组成，如式(4)：
C＝CS1+...+CSi+...+CSN+cp1+...+cpj+...+cpK (4)
式中C为粘连孢子轮廓；CSi为第i个轮廓段；N为轮廓段总数；cpj为第j个凹点；K为凹点总数。
基于凹度选取的轮廓凹点如图8c所示，删除凹点后，将粘连孢子边缘轮廓分割成多个轮廓段，通过连通区域标记，将轮廓段保存在一个有序的cell结构轮廓段集中，用于后续孢子轮廓段融合处理(图8d)。
8粘连孢子轮廓段融合
轮廓分割后同一孢子的轮廓可能被分割成多个轮廓段，所以，识别出同一孢子的轮廓段并融合，以获得各个孢子尽可能完整的椭圆拟合数据是本发明的关键。故本发明以轮廓段距离测量方法为准则确定候选轮廓段，对候选轮廓段进行最小二乘椭圆拟合并用偏移误差方法评价，将符合条件的轮廓融合成新的轮廓段，对新的轮廓段进行椭圆拟合即为最后的正确椭圆。
8.1轮廓段的距离测量方法
若2个轮廓段之间的距离很大，则属于同一孢子的可能性就很小，反之亦然。基于这一思想，本发明提出轮廓段的距离测量方法，判断轮廓段之间的位置关系，以排除属于同一孢子可能性小的轮廓段，确定候选轮廓段组合。第i和第j个轮廓段之间的距离DM定义为：
式中d(pi1,pj1)、分别为轮廓段CSi和CSj的2端点位置的欧氏距离；为轮廓段CSi和CSj的中点位置的欧式距离。
8.2最小二乘椭圆拟合算法
小麦条锈病菌夏孢子形状类似椭圆，故本发明用最小二乘椭圆拟合算法进行夏孢子轮廓拟合。
在二维平面坐标系中，一般的椭圆曲线方程可以用2个向量相乘的隐式方程来表示：
F(α,X)＝X·α＝ax2+bxy+cy2+dx+ey+f＝0 (6)
式中α＝[a b c d e f]T为椭圆方程的系数；x和y分别是曲线上点横、纵坐标，X＝[x2 xy y2 x y 1]。
为保证拟合结果为椭圆，需要加约束条件：
αTCα＝1 (7)
其中C为一个6×6的矩阵，
方程F(α,Xi)为二值图像中边缘点(xi,yi)到给定椭圆方程的代数距离。根据最小二乘原理，椭圆拟合问题通过将代数距离平方之和
最小化来解决。方程(9)改写为向量形式：
E＝||Dα||2 (10)
其中D为一个m×6的矩阵，
引入拉格朗日系数并微分可得：
式中λ为方程(12)的特征值；W为一个6×6的散射矩阵，
W＝DTD (13)
根据广义特征值求解方法可得到：
式中λi和ui分别为方程的特征值和特征向量。由此解得至多6个实数解(λi,αi)。
推导出数据点到椭圆的代数距离平方和公式：
E＝||Dα||2＝αTDTDα＝αTWα＝λαTCα＝λ (15)
由式(15)可知所需要的是最小的正特征值λi所对应的特征向量αi。得到解αi之后，即可实现对椭圆的拟合。
8.3轮廓段到拟合椭圆的偏移误差方法
通过候选轮廓段进行椭圆拟合时，并未考虑轮廓段与椭圆的拟合程度，导致得到的候选椭圆与实际有较大偏差，故需要对候选轮廓段进行评价筛选，本文以偏移误差作为去除错误候选轮廓段的评价条件，偏移误差公式为：
DEM(CS#,CE)＝E/M#＝λ/M# (16)
式中E为由式(15)求出的给定轮廓段CS#到椭圆CE的最小二乘代数距离；M#为轮廓段CS#上点的总数。
若偏移误差小于阈值，则该轮廓段属于同一孢子的轮廓段，否则不是。
8.4轮廓段融合步骤
以图8d为例对融合算法进一步说明，具体步骤为：
(1)首先选取图8d中最长的轮廓段CS8；
(2)依次计算CS8与剩余的CSi的距离DM，通过计算，CS8与CS4、CS5、CS6和CS7的距离测量值小于预设的阈值ωDM，故一个新的轮廓段集CS*包括CS4、CS5、CS6和CS7；
(3)然后依次作轮廓段(CS8和CS4、CS8和CS5、CS8和CS6、CS8和CS7)的候选拟合椭圆(如图8e所示)，同时求出各轮廓段的偏移误差值DEM；
(4)无轮廓段满足式DM(CS1,CSi)<ωDM，则只有CS8融合成一个新的轮廓段；
(5)用新的轮廓段拟合椭圆(如图8f所示)，同时孢子数量Num+1，并且删除轮廓段集中的CS8；
(6)重复上述步骤，直到所有轮廓段都拟合出正确的椭圆(如图8g所示)。
所述步骤8.4中的距离测量阈值ωDM和偏移误差阈值σDEM由30幅夏孢子图像样本的训练确定，分别为40和95。
9最后统计所有椭圆的总数即为夏孢子的个数。本发明对图2的分割计数结果如图9所示，从图9可以看出，粘连孢子被很好地分割出来。本发明的小麦条锈病菌夏孢子计数准确率结果如表1所示。应用本发明方法，最低计数准确率为92.7％，最高计数准确率为100％，总平均计数准确率为98.6％，具有较高的孢子计数精度。为在线式小麦条锈病夏孢子监测装备的开发提供技术支持。
表1本发明的小麦条锈病菌夏孢子计数准确率结果
本说明书中所描述的以上内容仅仅是对本发明所作的举例说明。本发明所述技术领域的技术人员可以对所描述的具体实施例做各种修改或补充或采用类似的方式替代，只要不偏离本发明说明说的内容或者超越本权利要求书所定义的范围，均应属于本发明的保护范围。
Claims (10)
1.一种小麦条锈病菌夏孢子自动统计计数方法，其特征在于，包括以下步骤：
S1、夏孢子显微图像采集
用孢子捕捉仪捕捉田间空气中的小麦条锈病菌夏孢子，并在显微镜下进行显微图像采集，获得供试夏孢子显微图像；
S2、显微图像选取
将所述步骤S1的原始夏孢子显微图像用双线性插值方法缩小尺寸，随机选取部分供试夏孢子图像用于算法训练，其余图像用于测试；
S3、K-means聚类分割
将所述步骤S2的夏孢子显微图像用K-means聚类算法实现孢子和背景的分割，即将其分为孢子目标和背景两大类，然后采用全局阈值算法进行二值化处理得到二值图像；
S4、形态学处理
首先用填充操作对所得二值图像进行孔洞填充，再用完整孢子中最小面积值为阈值，去除小面积区域噪声，所占面积小于该阈值的边界上孢子，均从图像中移除，计数时作为其他视野中的孢子，最后，用圆盘结构元素对二值图像进行形态学开运算，以消除夏孢子边界小的凸刺；
S5、基于形状因子和面积的粘连孢子判别
选择描述目标边界复杂程度的形状因子和面积作为粘连孢子判别的依据，形状因子公式为：
SF＝4πS/L2
式中SF为形状因子；S为一个连通区域的面积像素值；L为连通区域的周长像素值；
经粘连孢子判别后，对判定为单个孢子的区域直接采用最小二乘椭圆拟合算法进行拟合并记录椭圆个数Num供后续统计计数，对粘连孢子区域将每个孢子一一分割出来并自动计数；
S6、基于凹度的粘连孢子轮廓分割
对所述步骤S5二值图像中的粘连孢子区域通过移除内部像素点操作提取出边缘轮廓，然后遍历图像跟踪边缘轮廓点，并将轮廓点坐标保存在一个有序表中，pt(xt,yt)是边缘轮廓上任意一点，两向量ptpt-k和ptpt+k之间的夹角称为pt的凹度，pt-k和pt+k表示点pt的相邻轮廓点，凹度concavity的公式为：
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>n</mi>
<mi>c</mi>
<mi>a</mi>
<mi>v</mi>
<mi>i</mi>
<mi>t</mi>
<mi>y</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>a</mi>
<mi>n</mi>
<mi>g</mi>
<mi>l</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>arccos</mi>
<mfrac>
<mrow>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>&CenterDot;</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>p</mi>
<mi>t</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
其中凹点是满足如下两个条件的轮廓点：
(1)凹度concavity(pt)在角angle(δ1,δ2)范围之内；
(2)直线不包含在粘连孢子内；
粘连孢子边缘轮廓由轮廓上的所有凹点和由凹点分割成的多个轮廓段组成，如下式：
C＝CS1+...+CSi+...+CSN+cp1+...+cpj+...+cpK
式中C为粘连孢子边缘轮廓；CSi为第i个轮廓段；N为轮廓段总数；cpj为第j个凹点；K为凹点总数；
删除凹点后，将粘连孢子边缘轮廓分割成多个轮廓段，通过连通区域标记，将轮廓段保存在一个有序的cell结构轮廓段集中，用于后续孢子轮廓段融合处理；
S7、粘连孢子轮廓段融合
以轮廓段距离测量方法为准则确定候选轮廓段，对候选轮廓段进行最小二乘椭圆拟合并用偏移误差方法进行评价，将符合条件的轮廓段融合成新的轮廓段，对新的轮廓段进行椭圆拟合即为最后的椭圆；
S8、最后统计所有椭圆的总数即为夏孢子的个数。
2.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S1中，用TPBZ3型孢子捕捉仪捕捉田间空气中的小麦条锈病菌夏孢子，是在野外小麦田间，抽出捕捉仪的载波器并将表面均匀涂抹一薄层凡士林的载玻片放其上，涂面朝上，将载波器放入捕捉仪仓口内，为了使玻片上黏附的孢子密度不同，获取捕捉时间分别为120、180、240min的条锈病菌夏孢子载玻片，重复10次，共获不同孢子密度的载玻片30片；
用BX52型倒置显微镜，放大倍数为10×20，分别对30片载玻片进行显微图像采集，每片随机选取5个视野拍照，共获得150幅供试夏孢子图像，图像尺寸为4140×3096Pixel，图像分辨率为72dpi，存贮为bmp格式；
所述步骤S2中，将所述供试夏孢子显微图像用双线性插值方法缩小为911×682Pixel，然后随机选取30幅用于算法训练，其余120幅图像用于测试。
3.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S4中的圆盘结构元素，大小为3×3。
4.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S5中的基于形状因子和面积的粘连孢子判别，形状因子阈值SF0＝0.8，面积阈值Smax＝560，依次提取单个连通区域i，若i满足公式：
SFi>SF0&Si<Smax
则判定区域i为未粘连的单个夏孢子，否则为粘连孢子。
5.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S6中的凹度的角δ1和δ2分别设定为50°和150°。
6.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S7的轮廓段距离方法，第i和第j个轮廓段之间的距离DM定义为公式：
<mrow>
<mi>D</mi>
<mi>M</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>CS</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>CS</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>5</mn>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>d</mi>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mi>d</mi>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>jM</mi>
<mi>j</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mi>d</mi>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>im</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>jm</mi>
<mi>j</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>+</mo>
<mi>d</mi>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>iM</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mi>d</mi>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>iM</mi>
<mi>i</mi>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<msub>
<mi>jM</mi>
<mi>j</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中d(pi1,pj1)、分别为轮廓段CSi和CSj的2端点位置的欧氏距离；为轮廓段CSi和CSj的中点位置的欧式距离。
7.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S7的最小二乘椭圆拟合算法的具体步骤为：
在二维平面坐标系中，椭圆曲线方程用2个向量相乘的隐式方程来表示：
F(α,X)＝X·α＝ax2+bxy+cy2+dx+ey+f＝0
式中α＝[a b c d e f]T为椭圆方程的系数；x和y分别是曲线上点横、纵坐标，X＝[x2xy y2 x y 1]；
加约束条件：αTCα＝1，以保证拟合结果为椭圆，
其中C为一个6×6的矩阵，
方程F(α,Xi)为二值图像中边缘点(xi,yi)到给定椭圆方程的代数距离。根据最小二乘原理，椭圆拟合问题通过将代数距离平方之和
<mrow>
<mi>E</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mi>F</mi>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
最小化来解决，方程改写为向量形式：E＝||Dα||2，
其中D为一个m×6的矩阵，
引入拉格朗日系数并微分得：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>W</mi>
<mi>&alpha;</mi>
<mo>=</mo>
<mi>&lambda;</mi>
<mi>C</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>&alpha;</mi>
<mi>T</mi>
</msup>
<mi>C</mi>
<mi>&alpha;</mi>
<mo>=</mo>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
式中λ为方程特征值；W为一个6×6的散射矩阵，W＝DTD；
根据广义特征值求解方法得到：
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mfrac>
<mn>1</mn>
<mrow>
<msubsup>
<mi>u</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>Cu</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</msqrt>
<mo>)</mo>
</mrow>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mrow>
<msubsup>
<mi>u</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msub>
<mi>Su</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</msqrt>
<mo>)</mo>
</mrow>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
</mrow>
式中λi和ui分别为方程的特征值和特征向量，由此解得至多6个实数解(λi,αi)；
推导出数据点到椭圆的代数距离平方和公式为：
E＝||Dα||2＝αTDTDα＝αTWα＝λαTCα＝λ
由上式可知所需要的是最小的正特征值λi所对应的特征向量αi，得到解αi之后，即可实现对椭圆的拟合。
8.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S7的偏移误差公式为：
DEM(CS#,CE)＝E/M#＝λ/M#
式中E为由数据点到椭圆的代数距离平方和公式求出的给定轮廓段CS#到椭圆CE的最小二乘代数距离；M#为轮廓段CS#上点的总数；
若偏移误差小于阈值，则该轮廓段属于同一孢子的轮廓段，否则不是。
9.根据权利要求1所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S7的轮廓段融合步骤具体为：
(1)从轮廓段集中选取轮廓段中最长轮廓段CS1；
(2)将孢子轮廓段中剩余的轮廓段CSi分别与CS1作距离测量，设定距离测量阈值ωDM，若满足DM(CS1,CSi)<ωDM条件，则有
CS*＝{CSi|DM(CS1,CSi)<ωDM,i＝1,2,...,k}；
(3)将CS1依次与CS*中每一个轮廓段CSi作椭圆拟合，同时求出相应的候选椭圆和偏移误差DEM；
(4)若偏移误差DEM小于预设的阈值σDEM，则融合成新的轮廓段CSnew：
CSnew＝{CS1+CSi|DEM(CS1,CSi)<σDEM,i＝1,2,...,k}；
(5)对CSnew作椭圆拟合，作为一个孢子的正确拟合椭圆，孢子数量计数Num+1，并且删除轮廓段集中的CSnew；
(6)最后，判断轮廓段集中的轮廓段是否全部删除，若是，则输出所有正确拟合椭圆并且结束程序，否则返回到步骤1。
10.根据权利要求8所述的小麦条锈病菌夏孢子自动统计计数方法，其特征在于，所述步骤S7中的距离测量阈值ωDM和偏移误差阈值σDEM由30幅夏孢子图像样本的训练确定，分别为40和95。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711270040.8A CN108052886B (zh) | 2017-12-05 | 2017-12-05 | 一种小麦条锈病菌夏孢子自动统计计数方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711270040.8A CN108052886B (zh) | 2017-12-05 | 2017-12-05 | 一种小麦条锈病菌夏孢子自动统计计数方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108052886A true CN108052886A (zh) | 2018-05-18 |
| CN108052886B CN108052886B (zh) | 2018-12-18 |
Family
ID=62122070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711270040.8A Active CN108052886B (zh) | 2017-12-05 | 2017-12-05 | 一种小麦条锈病菌夏孢子自动统计计数方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108052886B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109344874A (zh) * | 2018-08-31 | 2019-02-15 | 华侨大学 | 一种基于深度学习的染色体自动分析方法及系统 |
| CN110208190A (zh) * | 2019-07-17 | 2019-09-06 | 四川沃文特生物技术有限公司 | 一种用于显微镜镜检追踪的方法及系统 |
| CN110458808A (zh) * | 2019-07-10 | 2019-11-15 | 山东仕达思生物产业有限公司 | 基于形态学与yolo算法的女性生殖道病原体识别方法 |
| CN111724379A (zh) * | 2020-06-24 | 2020-09-29 | 武汉互创联合科技有限公司 | 基于组合视图的显微图像细胞计数与姿态识别方法及系统 |
| CN111724378A (zh) * | 2020-06-24 | 2020-09-29 | 武汉互创联合科技有限公司 | 一种显微图像细胞计数与姿态识别方法及系统 |
| CN113567392A (zh) * | 2021-07-20 | 2021-10-29 | 西北农林科技大学 | 一种基于近红外光谱的小麦气传病原菌孢子快速无损识别方法 |
| CN113673396A (zh) * | 2021-08-10 | 2021-11-19 | 佛山市海天(高明)调味食品有限公司 | 一种孢子发芽率计算方法、装置及存储介质 |
| CN113869246A (zh) * | 2021-09-30 | 2021-12-31 | 安徽大学 | 基于改进CenterNet技术的小麦条锈病菌夏孢子显微图像检测方法 |
| CN118552539A (zh) * | 2024-07-29 | 2024-08-27 | 三峡金沙江云川水电开发有限公司 | 一种基于图像识别技术的座环法兰面垫块识别方法及系统 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103745231A (zh) * | 2014-01-16 | 2014-04-23 | 华南农业大学 | 小麦矮腥黑穗病tck及其近似种tct的冬孢子图像鉴定方法 |
| CN104091179A (zh) * | 2014-07-01 | 2014-10-08 | 北京工业大学 | 小麦白粉病菌孢子图片的智能识别方法 |
| CN104651462A (zh) * | 2015-01-29 | 2015-05-27 | 华南农业大学 | 基于显微图像分析的稻瘟病菌孢子检测方法 |
| CN105160675A (zh) * | 2015-08-31 | 2015-12-16 | 北京工业大学 | 白粉病菌孢子图像的自动分割方法 |
| CN105447859A (zh) * | 2015-11-18 | 2016-03-30 | 扬州大学 | 一种田间麦蚜虫计数方法 |
| CN106056118A (zh) * | 2016-06-12 | 2016-10-26 | 合肥工业大学 | 一种用于细胞的识别计数方法 |
| CN106295789A (zh) * | 2015-06-10 | 2017-01-04 | 浙江托普云农科技股份有限公司 | 一种基于图像处理的作物种子计数方法 |
| CN106855951A (zh) * | 2016-12-16 | 2017-06-16 | 陕西科技大学 | 一种基于计算机视觉的粮种品质检测方法 |
-
2017
- 2017-12-05 CN CN201711270040.8A patent/CN108052886B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103745231A (zh) * | 2014-01-16 | 2014-04-23 | 华南农业大学 | 小麦矮腥黑穗病tck及其近似种tct的冬孢子图像鉴定方法 |
| CN104091179A (zh) * | 2014-07-01 | 2014-10-08 | 北京工业大学 | 小麦白粉病菌孢子图片的智能识别方法 |
| CN104651462A (zh) * | 2015-01-29 | 2015-05-27 | 华南农业大学 | 基于显微图像分析的稻瘟病菌孢子检测方法 |
| CN106295789A (zh) * | 2015-06-10 | 2017-01-04 | 浙江托普云农科技股份有限公司 | 一种基于图像处理的作物种子计数方法 |
| CN105160675A (zh) * | 2015-08-31 | 2015-12-16 | 北京工业大学 | 白粉病菌孢子图像的自动分割方法 |
| CN105447859A (zh) * | 2015-11-18 | 2016-03-30 | 扬州大学 | 一种田间麦蚜虫计数方法 |
| CN106056118A (zh) * | 2016-06-12 | 2016-10-26 | 合肥工业大学 | 一种用于细胞的识别计数方法 |
| CN106855951A (zh) * | 2016-12-16 | 2017-06-16 | 陕西科技大学 | 一种基于计算机视觉的粮种品质检测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 焦林 等: "小麦条锈病孢子图像实时传输及处理系统设计", 《西北农林科技大学学报（自然科学版）》 * |
| 陶明超 等: "田间孢子捕捉及自动计数设计", 《内蒙古农业大学学报（自然科学版）》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109344874A (zh) * | 2018-08-31 | 2019-02-15 | 华侨大学 | 一种基于深度学习的染色体自动分析方法及系统 |
| CN109344874B (zh) * | 2018-08-31 | 2020-09-15 | 华侨大学 | 一种基于深度学习的染色体自动分析方法及系统 |
| CN110458808A (zh) * | 2019-07-10 | 2019-11-15 | 山东仕达思生物产业有限公司 | 基于形态学与yolo算法的女性生殖道病原体识别方法 |
| CN110458808B (zh) * | 2019-07-10 | 2021-06-15 | 山东仕达思生物产业有限公司 | 基于形态学与yolo算法的女性生殖道病原体识别方法 |
| CN110208190A (zh) * | 2019-07-17 | 2019-09-06 | 四川沃文特生物技术有限公司 | 一种用于显微镜镜检追踪的方法及系统 |
| CN111724379A (zh) * | 2020-06-24 | 2020-09-29 | 武汉互创联合科技有限公司 | 基于组合视图的显微图像细胞计数与姿态识别方法及系统 |
| CN111724378A (zh) * | 2020-06-24 | 2020-09-29 | 武汉互创联合科技有限公司 | 一种显微图像细胞计数与姿态识别方法及系统 |
| CN111724379B (zh) * | 2020-06-24 | 2024-05-24 | 武汉互创联合科技有限公司 | 基于组合视图的显微图像细胞计数与姿态识别方法及系统 |
| CN113567392A (zh) * | 2021-07-20 | 2021-10-29 | 西北农林科技大学 | 一种基于近红外光谱的小麦气传病原菌孢子快速无损识别方法 |
| CN113673396A (zh) * | 2021-08-10 | 2021-11-19 | 佛山市海天(高明)调味食品有限公司 | 一种孢子发芽率计算方法、装置及存储介质 |
| CN113869246A (zh) * | 2021-09-30 | 2021-12-31 | 安徽大学 | 基于改进CenterNet技术的小麦条锈病菌夏孢子显微图像检测方法 |
| CN118552539A (zh) * | 2024-07-29 | 2024-08-27 | 三峡金沙江云川水电开发有限公司 | 一种基于图像识别技术的座环法兰面垫块识别方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108052886B (zh) | 2018-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108052886A (zh) | 一种小麦条锈病菌夏孢子自动统计计数方法 | |
| US8605981B2 (en) | Centromere detector and method for determining radiation exposure from chromosome abnormalities | |
| CN103400151B (zh) | 一体化的光学遥感影像与gis自动配准与水体提取方法 | |
| CN109190567A (zh) | 基于深度卷积神经网络的异常宫颈细胞自动检测方法 | |
| CN112396619B (zh) | 一种基于语义分割的内部复杂组成的小型颗粒分割方法 | |
| CN107742099A (zh) | 一种基于全卷积网络的人群密度估计、人数统计的方法 | |
| CN111860459B (zh) | 一种基于显微图像的禾本科植物叶片气孔指数测量方法 | |
| CN102324032B (zh) | 一种基于极坐标系下灰度共生矩阵的纹理特征提取方法 | |
| CN108537751B (zh) | 一种基于径向基神经网络的甲状腺超声图像自动分割方法 | |
| CN105825169B (zh) | 一种基于道路影像的路面裂缝识别方法 | |
| CN105844228A (zh) | 一种基于卷积神经网络的遥感图像云检测方法 | |
| CN105930852B (zh) | 一种气泡图像识别方法 | |
| CN111191628B (zh) | 基于决策树与特征优化的遥感影像震害建筑物识别方法 | |
| CN108596165A (zh) | 基于无人机低空航拍图像的道路交通标线检测方法及系统 | |
| CN109684938A (zh) | 一种基于作物冠层航拍俯视图的甘蔗株数自动识别方法 | |
| CN105046259B (zh) | 基于多特征融合的日冕物质抛射检测方法 | |
| CN102346854A (zh) | 前景物体检测方法和设备 | |
| CN115115627B (zh) | 基于数据处理的土壤盐碱地监测方法 | |
| CN113435460A (zh) | 一种亮晶颗粒灰岩图像的识别方法 | |
| CN108629297A (zh) | 一种基于空域自然场景统计的遥感图像云检测方法 | |
| CN105719294A (zh) | 一种乳腺癌病理学图像有丝分裂核自动分割方法 | |
| CN111882573B (zh) | 一种基于高分辨率影像数据的耕地地块提取方法及系统 | |
| CN110197113A (zh) | 一种高精度锚点匹配策略的人脸检测方法 | |
| CN112183237A (zh) | 基于颜色空间自适应阈值分割的白细胞自动分类方法 | |
| CN117853722A (zh) | 一种融合超像素信息的钢材金相组织分割方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |