基于改进算法的小麦收获机作物密度检测研究
- 格式:pdf
- 大小:238.59 KB
- 文档页数:4
基于改进算法的小麦收获机作物密度检测研究刘仲鹏1,赵建军1,刘丽娟2(1.保定学院信息技术系,河北保定071000;2.河北农业大学信息学院,河北保定071000)
摘要:针对谷物联合收获机喂入量实时调整的需求,研究成熟期作物密度与收获机喂入量间的函数关系。首先对图像进行灰度化与中值滤波处理,然后以改进的图像分割算法获得了更为准确的图像分割阈值,从而将作物的谷、叶从背景中分割出来。在此基础上,以数据拟合的方式得到小麦收获机喂入量随图像像素的变化模型,为谷物收获机实时、精准地控制喂入量提供可行的方案。关键词:小麦;联合收获机;图像处理;作物密度;改进算法中图分类号:S225.3;TP391.41文献标识码:A文章编号:1003-188X(2015)08-0049-04
0引言
谷物联合收获机的喂入量是一项关键的性能参数,也是影响收获机效率的首要因素。过低的喂入量会导致工作效率低下,降低生产率;过高的喂入量则会提升谷物的损失率,甚至导致脱粒滚筒堵塞,发生故障。在作业过程中,收获机喂入量的实时控制与调整十分重要。传统的调整模式一般是根据驾驶员的经验和感觉进行的,由于主观性太强,对驾驶员经验要求较高,难以使联合收获机维持在最佳喂入量水平。如何实现收获机对喂入量的实时、自动调节,是一个研究的热点[1-3]。本文从计算机视觉的角度进行研究,以成熟期小麦为例,通过实时获取被收获作物的冠层图像,消除图像样本外在因素的影响,提取作物冠层的图像特征,以改进的阈值分割方法获得谷、叶图像区域,并分析作物实际密度与图像区域像素数之间的函数关系,从而为收获机喂入量的智能、实时调节提供依据。1图像处理与分割1.1图像采集以华北地区比较常见的郑麦7698作为研究对象,拍摄地点为河北省保定市河北农业大学小麦教研基地。以普通尼康数码相机作为采集器件,在晴天14:00-15:00时,以相同的高度与随机的方式采集40收稿日期:2014-08-02基金项目:河北省科技计划项目(12227403);保定市科技局科学技术研究与发展指导计划项目(13ZG033)作者简介:刘仲鹏(1978-),男,河北张家口人,讲师,硕士,(E-mail)lzp_bdxy@163.com。幅图片,每幅图片面积相等。所采集的原始图像为含有R、G、B等3个分量的24位真彩图。1.2数据采集对每幅图片范围内的小麦分别进行收割、脱粒,称取其质量并记录,共计40组数据。部分数据(前5组)如表1所示。其中,将谷质量+草质量/3(kg)作为收获机的喂入量[4]。
1.3图像预处理
1)灰度处理。所采集的小麦冠层原始图像属于
24位真彩图,含有R、G、B等3个分量。为了使后续
处理计算量减少,其进行灰度化处理,通过3个步骤实现:①分别提取目标图像RGB的具体值;②求出具体的灰度值gray;③令R=B=G=gray,最终实现图像灰度化。2)中值滤波处理。在对小麦冠层图像进行处理
时,还应考虑到原始图像中可能由于一些外部干扰而掺杂的噪声成分。噪声的存在能够对图像质量造成比较直接的影响,进而降低特征的提取准确性。为此,本文在灰度预处理的基础上,继续以中值滤波对图像进行进一步的处理。中值滤波的主要原理:将图像分为不同的窗口S,窗口里的所有像素进行大小排列,将中位灰度值作为该窗口的灰度值。中值滤波已经被证实能够较为有效地避免图像细节的模糊,并滤除颗粒噪声及脉冲噪声,尽量保持图像边缘的有用信息。本研究将滤波模板定位于3×3,实现样本图像的增强。1.4改进图像分割算法
此步骤将小麦冠层从图像背景中分离出来,即进行小麦作物的谷、叶与背景的分割,以便于下一步的特征提取及分析。本文分割图像采用的是阈值法,得
·94·
2015年8月农机化研究第8期
DOI:10.13427/j.cnki.njyi.2015.08.011到的是二值图像。其基本思路是:对于图像的每一个点,高于该阈值则取1,低于该阈值则取0。考虑到小麦作物的谷与叶颜色主要为黄色及绿色,而土壤背景则为褐色,所以本研究提取图像中的颜色特征来实现分割。结合小麦的生长环境以及采集到的图像样本的特点,由于小麦冠层与图像背景之间的差异比较显著,其灰度直方图体现为双峰状。适合的分割方法包括熵最大值法、最大类间方差法及双峰法[5]等。
表1样本数据(部分)Table1Sampledata
编号总质量/kg谷物质量/kg草质量/kg(谷质量+草质量)/3/kg体积
/m3密度/kg·m-3
11.71.00.71.30.008212.521.91.20.71.40.010190.032.11.30.81.30.010210.041.81.10.71.40.009200.051.91.10.81.30.009211.1
本研究选取最大类间方差法(又名大津法),其优势在于具有较高的稳定性和较快的速度,并且常用于农林类图像处理和分割中。但经典的最大类间方差法也有一些不足之处,容易出现错分现象[6-7]。本研
究用改进的方法对其进行优化,使分割阈值更加合
理,实现更好的分割效果。传统的最大类间方差法是以灰度均值来表示目标和背景的。设目标图像幅度为M×N,且将其(x,y)处的灰度表示为f(x,y),灰度共计L个级别,则有
f(x,y)∈[0,L-1]。将某像素灰度取k的概率以
P(k)表示,可知P(k)=1MN∑
f(x,y)=k1。将目标图像的
所有像素分为2个类别:目标与背景。图像分割阈值设为t,则目标区域可表示为{f(x,y)≤t},背景区域可表示为{t<f(x,y)≤L-1}。以w0(t)表示目标区域所占比例,则w0(t)=
∑0≤k≤tp(k);以w
1(t)表示背景区域所占比例,则w1(t)
=∑t<k≤L-1p(k);以μ
0(t)表示目标区域平均值,则
μ0(t)=∑0≤k≤tkp(k)/w0(t);以μ1(t)表示背景区域平
均值,则μ1(t)=∑t<k≤L-1kp(k)/w1(t)。因此,分割阈
值为g=ArgMax0<t≤L-1{w0(t)[μ0(t)-μ]2+w1(t)
[μ
1(t)-μ]2}。
传统最大类间方差法仅关注了图像灰度分布均值,而描述图像灰度分布的参数除了均值还有方差。为此,文献[8]对传统最大类间方差法进行了优化,以
平均方差代替传统方法中的均值。其思路是结合被
处理对象的灰度直方图,计算图像背景和图像本身的方差最大值,从而获取分割阈值;而图像对象的灰度
平均方差主要体现的是被处理目标的图像均匀度。通常情况下,图像对象和图像背景往往在靠近内部的地方分布较为均匀,但在图像的边界则经常会发生灰度突变的现象。所以,通过平均方差来体现目标边缘的灰度突变,使选取阈值时能够兼顾最大的类间距与最小的类内距。算法为g=ArgMax0≤t≤L-1[σ20(t)-σ2]2[σ21(t)-σ2]
2
其中,σ20(t)=1wo(t)∑0≤k≤t[k-μ0(t)]
2p(k)表示
被处理图像的目标区域方差;σ
2
1(t)=
1
w1(t)
∑t+1≤k≤L-1[k-μ1(t)]
2
p(k)表示被处理图像的背景区域
方差;σ2=∑0≤k≤L-1[k-μ(t)]2p(k)表示被处理图像
的总方差。本研究在此算法的基础上,进一步为“类内距”分别配置一个权值,从而实现更加合理的阈值。算法为
g=Arg{Max0≤t≤L-1[σ20(t)-σ2]2[σ21(t)-σ2]2}
mσ
20(t)+nσ2
1(t)
其中,m=w0(t)
的含义是被处理图像的目标区
域所占比例;n=w1(t)
的含义是被处理图像的背景区
域所占比例。可见,在文献[8]的基础上,配置了m、n作为权值,结合目标区域和背景区域的实际像素数目,灵活调整目标与背景平均方差所占比例,从而更为合理地确定阈值[9]。
1.5图像分割
以改进的最大类间方差法对图像样本进行处理,随机从所拍摄的小麦作物图像中选取5幅作为样本,以人工阈值选取作为参照,得到的结果如图2所示。
·05·
2015年8月农机化研究第8期表2阈值选取结果Table2Thethresholdvalue
图像样本人工选取阈值经典最大类间方差法阈值及时间/s本研究改进最大类间方差法阈值及时间/s
197106/0.1002100/0.0923288105/0.098889/0.081237156/1.100272/1.024548784/0.070385/0.0625510196/1.0034102/0.0981算法的处理结果显示:本研究的改进最大类间方差法和传统算法相比,与人工选取阈值更为接近,图像样本2、3、5的阈值仅相差1个单位;而作为对比,经典最大类间方差法与人工选取阈值则相差较大。以改进算法对图像进行分割,可获得满意效果。图1所示为某样本图像的分割结果灰度图。图1图像分割效果Fig.1Imagesegmentationeffect2图像像素与小麦密度关系在获取图像样本的过程中,还应考虑到小部分被作物叶片遮蔽的麦穗,因此还应进一步计算小麦密度与像素值之间的合理关系。本研究以数据拟合的方式,选取DPS软件进行数据处理。表3所示为部分数据。以d表示小麦作物的密度值,以p表示采集图像
的像素值,拟合处理之后所得到的函数关系式为d=89.6771+0.000495p。图2所示为函数曲线。
图2小麦密度和提取图像像素值关系曲线Fig.2Relationshipsofwheatdensityandpixelvalues可见,小麦密度和提取图像像素值关系基本为直线,相关性令人满意。计算相关系数得R2=0.9025,
显著性良好,因此能够用提取图像像素值来表示小麦
密度。谷物联合收获机的喂入量与三大因素有关[10],分别是小麦密度、小麦品种以及喂入面积。由此可知,在小麦品种以及喂入面积一定的情况下,可以通过小麦密度实现喂入量的实时调节。表3地人工测算的小麦密度和提取图像像素值Table3Artificialmeasuringdensityofwheatandpixelvalues
序号像素数计算密度/kg·m-3实际密度/kg·m-3误差/kg·m-3
1233187205.1047211.60476.52149238163.5499167.84994.33156924167.3545163.7545-3.64188712183.0895193.389510.35267387222.0337216.5337-5.5
·15·
2015年8月农机化研究第8期