CN109145848B - 一种麦穗计数方法 - Google Patents
一种麦穗计数方法 Download PDFInfo
- Publication number
- CN109145848B CN109145848B CN201811002276.8A CN201811002276A CN109145848B CN 109145848 B CN109145848 B CN 109145848B CN 201811002276 A CN201811002276 A CN 201811002276A CN 109145848 B CN109145848 B CN 109145848B
- Authority
- CN
- China
- Prior art keywords
- image
- ear
- wheat
- area
- color space
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/50—Context or environment of the image
- G06V20/52—Surveillance or monitoring of activities, e.g. for recognising suspicious objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- 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
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/56—Extraction of image or video features relating to colour
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Analysis (AREA)
Abstract
本发明提一种麦穗计数方法,包括以下步骤:(1)采集麦穗图像;(2)对采集的麦穗图像进行颜色空间转换;(3)对颜色空间转换后的图像进行图像分割,使得图像内的内容被分割为麦穗、耕地、麦叶和茎秆三部分;(4)对分割处理后的图像进行二值化处理,以使麦穗单独从图像中分离出来;(5)对二值化图像进行滤波,滤除掉面积过大或过小的二值化连通区域,然后根据过滤后的二值化图像来统计麦穗数。本发明方法通过颜色特征,建立从图像低层特征到图像中包含麦穗数的一个直接映射关系,从而不需要再对目标进行分割或检测。这类方法更加适应于复杂的光照和密集麦穗的场景,相比基于图像分割的方法具有更高的计算效率和准确率。
Description
技术领域
本发明涉及作物生长信息测量技术领域,尤其涉及一种麦穗计数方法。
背景技术
单位面积麦穗数是小麦产量和种植密度构成的重要因素之一,小麦种植密度不仅影响小麦生长对水、肥料的需求,而且影响小麦群体内的通风和透光,进而影响到小麦个体的生长发育,从而影响小麦的产量。所以,单位面积的麦穗数是实现局部范围小麦产量预测和种植密度评估的重要参数之一。
目前国内外已有研究学者将图像处理和视频监控技术应用于作物生长监控中,并且对小麦产量预测做了一些有意义的研究工作,尤其是在基于图像处理技术的小麦麦穗自动计数方面做了一些有益的探索。国内,一些学者提出利用图像的颜色特征和纹理特征对大田环境下麦穗图像进行分割,该方法在麦穗图像没有粘连时计数精度较高,当麦穗粘连时该方法精度较差。国内另外一些学者提出提取麦穗的颜色、纹理特征,采用SVM学习的方法,精确提取小麦麦穗轮廓,同时构建麦穗特征数据库,对麦穗的二值图像细化得到麦穗骨架,最后通过计算麦穗骨架的数量以及麦穗骨架有效交点的数量,即可得到图像中麦穗的数量。国外一些学者提出基于RGB图像,利用边缘检测、膨胀腐蚀、孔洞区域填充等方法对麦穗进行计数,并且对在不同生长期条件下生长的五种无芒小麦品种进行检测试验,取得了不错的效果。另有一些国外学者提出利用拉普拉斯频域滤波器去除图像中高频和低频部分,接着用中值滤波器去除麦穗图像的噪声,最后找到图像的局部峰值的个数,从而实现麦穗的计数。这些方法都存在效率低,计数准确率不高等问题。
发明内容
本发明的目的在于解决上述现有技术存在的缺陷,针对传统的基于图像颜色特征、纹理特征、Haar特征的麦穗图像分割计数方法,本发明提出一种基于密度聚类和面积回归的麦穗计数方法。
一种麦穗计数方法,包括以下步骤:
(1)采集麦穗图像;
(2)对采集的麦穗图像进行颜色空间转换;
(3)对颜色空间转换后的图像进行图像分割,使得图像内的内容被分割为麦穗、耕地、麦叶和茎秆三部分;
(4)对分割处理后的图像进行二值化处理,以使麦穗单独从图像中分离出来;
(5)对二值化图像进行滤波,滤除掉面积过大或过小的二值化连通区域,然后根据过滤后的二值化图像来统计麦穗数。
进一步地,如上所述的麦穗计数方法,步骤(2)包括对采集的麦穗图像采用RGB颜色空间转换为超绿色波段EG空间
EG=2g-r-b (1)
式中,r、g、b是RGB图像的红色、绿色和蓝色分量。
进一步地,如上所述的麦穗计数方法,步骤(3)中采用密度聚类算法对颜色空间转换后的图像进行图像分割,所述密度聚类算法包括以下步骤:
步骤1:按照公式(2)计算数据xi与数据集中其它数据点xi之间的距离dij;
步骤3:按照公式(3)计算第j个数据点的局部密度ρj;
步骤4:计算第j个数据点的最近邻类距离δj,所述ε-邻域最近邻类距离δj为xj距最近的局部密度比它小的邻域距离;
步骤5:计算tj=ρjδj,根据tj值的大小,选取合适的聚类中心;
步骤6:根据选取的聚类中心,按照距离最近的原则,将数据分配到每个类别。
进一步地,如上所述的麦穗计数方法,所述步骤(5)中滤波采用面积滤波法来滤除掉面积过大或过小的二值化连通区域,该面积滤波法包括:
根据莱特准则,设计面积滤波器,当第i个连通区域ui的面积si满足
则认为ui是一个麦穗区域,于是麦穗的计数公式为
有益效果:
本发明方法通过颜色特征,建立从图像低层特征到图像中包含麦穗数的一个直接映射关系,从而不需要再对目标进行分割或检测。这类方法更加适应于复杂的光照和密集麦穗的场景,相比基于图像分割的方法具有更高的计算效率和准确率。
附图说明
图1(a)为本发明用于验证密度聚类算法性能的数据分布图;
图1(b)为本发明用于验证密度聚类算法性能的(ρj,δj)分布图;
图2(a)为3幅不同的用于麦穗计数的原始图;
图2(b)为与图2(a)对应的3幅不同的麦穗反色图;
图3(a)为与图2(b)对应的3幅不同的基于密度聚类算法对麦穗图像分割结果图;
图3(b)为与图3(a)对应的3幅不同的二值化麦穗图;
图4为利用本发明方法标记的与图2(a)对应的3幅不同的麦穗计数图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种麦穗的统计方法,具体包括以下步骤:
1.颜色空间的选择:颜色信息是小麦的重要特征,对于成长期的小麦,其颜色与背景有显著的差异。按照颜色的显著不同,麦穗图像分为麦穗、耕地、麦叶和茎秆三部分。RGB颜色空间是最常用的一种,但由于在颜色空间上对图像进行处理和分割时,RGB空间存在不直观和不均匀的缺点,因此可以考虑用绿色波段G和超绿色波段EG更加显著的颜色空间来进行麦穗分割。超绿色波段EG可以用式(1)表示。
EG=2g-r-b (1)
式中,r、g、b是RGB图像的红色、绿色和蓝色分量。
2.密度聚类算法:本发明提出基于密度的聚类算法,首先假定类别可以通过样本分布的紧密程度决定,即密度大小决定。同一类别的样本,他们之间是紧密相连的,也就是说,这个类别围绕一个聚类中心,其它样本紧紧围绕它。通过将紧密相连的样本划为一类,这样就得到了一个聚类类别。通过将所有各组紧密相连的样本划为各个不同的类别,则我们就得到了最终的所有聚类类别结果。为了刻画这个算法的思想,首先定义几个概念,即局部密度和最近邻类距离。
假设待聚类的数据集Is={i=1,2,…,N}为相应的指标集,该方法是基于一组邻域来描述样本集的紧密程度的,参数(ε,MinPts)用来描述邻域的样本分布紧密程度,其中ε描述了某一样本的邻域距离阈值,MinPts描述了某一样本的距离为ε的邻域中样本个数的阈值。
(1)ε-邻域:对于xj∈D,其ε-邻域包含样本集D中与xj的距离不大于ε的子样本集,即Nε(xj)={xi∈D|distance(xi,xj)≤ε}这个子样本集的个数记为Nε(xj)。
(2)核心对象:对于任一样本xj∈D,如果其ε-邻域样本集Nε(xj)至少包含MinPts个样本,即如果Nε(xj)≥MinPts,则xj是核心对象。
(3)局部密度:对于ε-邻域,其核心样本集Nε(xj)以xj为中心,其它样本xi与xj之间的距离dij用欧式距离表示:
于是定义ε-邻域的局部密度ρj用式(3)表示,
(4)最近邻类距离:分别以xj∈D为中心构建ε-邻域,该邻域的局部密度为ρj,按照局部密度大小由大到小排序,定义ε-邻域最近邻类距离δj为xj距最近的局部密度比它小的邻域距离。
根据以上定义,对于每个样本点xj,相应的可以求出(ρj,δj),设tj=ρjδj,选tj值明显大的xj值为聚类中心。具体基于密度聚类算法的步骤描述如下:
a)按照公式(2)计算数据xi与数据集中其它数据点xi之间的距离dij;
c)按照公式(3)计算第j个数据点的局部密度ρj;
d)计算第j个数据点的最近邻类距离δj
e)计算tj=ρjδj,根据tj值的大小,选取合适的聚类中心;
f)根据选取的聚类中心,按照距离最近的原则,将数据分配到每个类别。
本发明提出的密度聚类算法不需要预先设定聚类数和初始聚类中心,该算法能自动得到聚类中心和聚类数。图1(a)是一组数据的分布图,用本发明密度聚类算法对数据集D进行聚类。首先求出数据集D中每个数据点的(ρj,δj),其分布图1(b)所示。由图1(b)可以看出,第8个数据点、第16个数据点、第12个数据点所对应的(ρj,δj)值明显大于其他数据点所对应的值,所以选取数据点8、16、12做为聚类中心。与图1(a)数据分布图对比分析,本发明计算结果与原图分布完全一样。
3.基于密度聚类的麦穗图像分割:将获取的麦穗图像,从RGB颜色空间转换到超绿颜色空间(转换公式见式(1))。由麦穗原始图像和它的反色图像可以看出(见图2(a)、图2(b)),麦穗图像分为麦穗、耕地、叶子和茎秆三部分。用提出的密度聚类算法对麦穗图像进行分类,分类结果如图3(a)所示,由图3(a)可以看出,利用本发明提出的密度聚类算法能将麦穗、耕地、麦叶和茎秆很好的区分出来,麦穗的图像用绿色表示、耕地的颜色用蓝色表示、叶子和茎秆用红色表示。这样麦穗图像被划分成三个部分,将其中麦穗部分单独提取出来并进行二值化,如图3(b)所示,图3(b)中二值化图像大部分区域剩余麦穗。
4.麦穗计数:由图3(b)可以看出,图中绝大部分二值化区域是麦穗图像,少部分过小或过大的二值化区域是由叶子影响造成,为了准确进行麦穗计数,需要过滤掉这些过小或过大的二值化区域。为此本发明提出用面积滤波来进行麦穗的计数。
根据莱特准则,设计面积滤波器,当第i个连通区域ui的面积si满足
就认为ui是一个麦穗区域,于是麦穗的计数公式为,
图4是最终的麦穗标记计数结果图。
实施例:
首先,图像采集。用手机或数码相机正对着麦田,距离地面80cm,调整好焦距进行拍照,图像要求清晰、没有出现抖动模糊;
其次,颜色空间转换。按照公式(1)将麦穗图像的RGB颜色空间转换为超绿色波段EG空间,便于进行麦穗分割;
再其次,图像分割和二值化。用本发明提出的基于密度聚类方法对图像进行分割,并进行二值化,使得麦穗单独从图像中分离出来。
再其次,面积滤波。用本发明提出的面积滤波器对二值化图像进行滤波,滤除掉面积过大或过小的二值化连通区域。
最后,麦穗计数。运用公式(7),完成对麦穗的计数。
为了验证本专利方法的有效性,将本发明提出的麦穗算法计数与人工计数精度进行比较。运用统计误差、统计精度进行综合比较分析,见表1。由表1可以统计出,统计的平均测量精度是94.69%,统计误差在2.17%~7.41%之间,平均为5.31%。以上比较分析表明,说明利用本发明方法可以有效地进行麦穗自动计数,且精度高。
表1人工计数与算法计数比较
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种麦穗计数方法,其特征在于,包括以下步骤:
(1)采集麦穗图像;
(2)对采集的麦穗图像进行颜色空间转换;
(3)对颜色空间转换后的图像进行图像分割,使得图像内的内容被分割为麦穗、耕地、麦叶和茎秆三部分,步骤(3)中采用密度聚类算法对颜色空间转换后的图像进行图像分割,所述密度聚类算法包括以下步骤:
步骤1:按照公式(2)计算数据xi与数据集中其它数据点xi之间的距离dij;
步骤3:按照公式(3)计算第j个数据点的局部密度ρj;
步骤4:计算第j个数据点的最近邻类距离δj,ε-邻域最近邻类距离δj为比xj距最近的局部密度小的邻域距离;
步骤5:计算tj=ρjδj,根据tj值的大小,选取合适的聚类中心;
步骤6:根据选取的聚类中心,按照距离最近的原则,将数据分配到每个类别;
(4)对分割处理后的图像进行二值化处理,以使麦穗单独从图像中分离出来;
(5)对二值化图像进行滤波,滤除掉面积过大或过小的二值化连通区域,然后根据过滤后的二值化图像来统计麦穗数。
2.根据权利要求1所述的麦穗计数方法,其特征在于,步骤(2)包括对采集的麦穗图像采用RGB颜色空间转换为超绿色波段EG空间
EG=2g-r-b (1)
式中,r、g、b是RGB图像的红色、绿色和蓝色分量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811002276.8A CN109145848B (zh) | 2018-08-30 | 2018-08-30 | 一种麦穗计数方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811002276.8A CN109145848B (zh) | 2018-08-30 | 2018-08-30 | 一种麦穗计数方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109145848A CN109145848A (zh) | 2019-01-04 |
CN109145848B true CN109145848B (zh) | 2022-01-28 |
Family
ID=64829360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811002276.8A Active CN109145848B (zh) | 2018-08-30 | 2018-08-30 | 一种麦穗计数方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109145848B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766690B (zh) * | 2019-11-07 | 2020-08-14 | 四川农业大学 | 基于深度学习点监督思想的麦穗检测和计数方法 |
CN113011220A (zh) * | 2019-12-19 | 2021-06-22 | 广州极飞科技股份有限公司 | 穗数识别方法、装置、存储介质及处理器 |
CN111259925B (zh) * | 2020-01-08 | 2023-05-02 | 安徽大学 | 基于k均值聚类和宽度突变算法的田间麦穗计数方法 |
CN111369494B (zh) * | 2020-02-07 | 2023-05-02 | 中国农业科学院农业环境与可持续发展研究所 | 冬小麦穗密度检测方法及装置 |
CN113628186A (zh) * | 2021-08-09 | 2021-11-09 | 海南青峰生物科技有限公司 | 一种基于5g通讯及图像识别的测量水稻穗数的*** |
CN115330786B (zh) * | 2022-10-13 | 2023-04-25 | 南京邮电大学 | 基于ceci算法的水稻植株深度学习计数数据集的创建方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2949984B1 (fr) * | 2009-09-15 | 2023-09-29 | Francois Paradis | Dispositif d'ordonnance de grains pour en realiser une caracterisation par imagerie et traitement d'image |
CN104392240A (zh) * | 2014-10-28 | 2015-03-04 | 中国疾病预防控制中心寄生虫病预防控制所 | 一种基于多特征融合的寄生虫虫卵识别方法 |
CN106373136A (zh) * | 2015-07-22 | 2017-02-01 | 崔胡晋 | 一种基于颜色和特征的害虫识别计数方法 |
CN106296695B (zh) * | 2016-08-12 | 2019-05-24 | 西安理工大学 | 基于显著性的自适应阈值自然目标图像分割抽取算法 |
CN107590812B (zh) * | 2017-09-01 | 2021-07-02 | 南京农业大学 | 小麦穗部小穗识别计数方法 |
-
2018
- 2018-08-30 CN CN201811002276.8A patent/CN109145848B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109145848A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109145848B (zh) | 一种麦穗计数方法 | |
Safren et al. | Detection of green apples in hyperspectral images of apple-tree foliage using machine vision | |
CN106951836B (zh) | 基于先验阈值优化卷积神经网络的作物覆盖度提取方法 | |
CN105718945B (zh) | 基于分水岭和神经网络的苹果采摘机器人夜间图像识别方法 | |
CN109447945B (zh) | 基于机器视觉和图形处理的小麦基本苗快速计数方法 | |
CN110309781B (zh) | 基于多尺度光谱纹理自适应融合的房屋损毁遥感识别方法 | |
CN110569747A (zh) | 一种利用图像金字塔与Faster-RCNN快速统计大田水稻稻穗数的方法 | |
CN106934418B (zh) | 一种基于卷积递归网络的绝缘子红外诊断方法 | |
CN106845497B (zh) | 基于多特征融合的玉米前期图像旱情识别方法 | |
CN105184808B (zh) | 一种光场图像前后景自动分割方法 | |
CN109063754A (zh) | 一种基于OpenStreetMap的遥感影像多特征联合分类方法 | |
CN111259925B (zh) | 基于k均值聚类和宽度突变算法的田间麦穗计数方法 | |
CN109187552B (zh) | 一种基于云模型的小麦赤霉病害等级判定方法 | |
CN103164695A (zh) | 一种基于多源图像信息融合的果实识别方法 | |
Ji et al. | In-field automatic detection of maize tassels using computer vision | |
CN109871900A (zh) | 一种基于图像处理的复杂背景下苹果的识别定位方法 | |
CN106845366B (zh) | 基于图像的甘蔗覆盖度自动检测方法 | |
CN110503140A (zh) | 基于深度迁移学习与邻域降噪的分类方法 | |
CN110070545A (zh) | 一种城镇纹理特征密度自动提取城镇建成区的方法 | |
CN112053371A (zh) | 一种遥感图像中的水体提取方法和装置 | |
CN116524224A (zh) | 一种基于机器视觉的烤后烟叶类型检测方法及*** | |
CN111046838A (zh) | 一种湿地遥感信息的识别方法及装置 | |
CN106023223A (zh) | 柑橘果实大小描述及分级方法 | |
CN110047064B (zh) | 一种马铃薯疮痂检测方法 | |
CN114119634B (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 |