基于小波变换及otsu分割的农田作物行提取-韩永华.pdf

《基于小波变换及otsu分割的农田作物行提取-韩永华.pdf》由会员分享，可在线阅读，更多相关《基于小波变换及otsu分割的农田作物行提取-韩永华.pdf（8页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、第38卷第1期2016年1月电子与信息学报Journal of Electronics&Information TechnologyVdl38No1Jan2016基于小波变换及Otsu分割的农田作物行提取韩永华 汪亚明 孙麒u 赵匀蛔(浙江理工大学信息学院杭州 310018)(浙江理工大学机械与自动控制学院杭州 310018)摘要：基于机器视觉的田间车辆自动导航是农用车辆导航的热门研究方向，但含较密集杂草的农田作物行提取，目前依然是个难题。该文提出一种适用丁密集杂草农田的，主要基于频率和颜色信息的农田图像分割算法。通过小波多分辨率分解后构建的频率总量指标，利用作物产生主频信息的总量优势，结合作

2、物行的交替及最大类间方差法、颜色模型分量变换，实现农田杂草的去除，并通过最小二乘法拟合直线，实现农田作物行提取。实验表明算法能有效克服密集杂草干扰，针对480X 640像素大小图像，单幅处理时间平均为132 ms。关键词：田间导航；作物行；小波变换；最大类问方差中图分类号：TP39141 文献标识码： A 文章编号：10095896(2016)01006308DOI：1011999JEITl50421Crop Row Detection Based on肠Felet Transformation andOtsu Segmentation AlgorithmHAN Yonghua WANG Ya

3、ming SUN Qi ZHAO Yun(Department劈Electronics and Informatics，Zhejian9 ScC Tech University，Hangzhou 310018，China)”(Department of Mechanic，Zhejiang sc乒Tech University，Hangzhou 310018，China)Abstract：Visionbased agricultural vehicle navigation has become a popular research area of automated guidance，howe

4、vercrop row detection in hi【gh weeds field is still a challenging topicAn image segmentation method mainlybased on frequency and color information is proposed to remove weedsThe algorithm is based on total frequencyparameters，more total crop frequency，alternation regular of crop rows，Otsu method and

5、 color modeltransformation，The total frequency parameters are obtained from wavelet multi-resolution decompositionTheleast square method is used in fitting straight line to detect the crop rowsExperiments show that the algorithmearl effectively overcome the high weedsThe average processing time of a

6、 single 480 X 640 pixels image is 132 msKey words：Agriculture navigation；Crop rows；Wavelet transformation；Otsu1 引言田问车辆自动导航、辅助导航能有效降低农民工作强度，提高农业生产效率，减少农药中毒、皮收稿日期：201504-10；改回日期：20150913；网络出版：201511184通信作者：赵匀zhaoyunzstu，edu，cn基金项目：国家自然科学基金(61272311)，浙江省自然科学基金重点项目fLZl5F020004)，机械工程浙江省高校重中之重学科和浙江理工大学重点实

7、验室优秀青年人才培养基金(ZSTUME01817)，训算机应用创新重点学科研究生创新研究项目(XDYl4003)，浙江理工大学科研启动基金(13032156Y)，浙江理工大学521资助计划项目Foundation Items：The National Natural Science Foundation ofChina(612723111，Zhejiang Provincial Natural Science Foundation(LZl5F020004)，The Young Researchers Foundation of ZhejiangProvincial Top Key Academ

8、ic Discipline of Mechanical Engineeringand Zhejiang SciTech University Key Laboratory(zSTuME01817)Graduate Student Innovation Research Project ofComputer Application Innovation Key Subject(XDYl4003)，Science Foundation of Zhejiang SciTech University(ZSTU)(13032156一Y)，521 Project of Zhejiang Sci-Tech

9、University肤晒伤等农业事故的发生，并能促进精细农业发展，从而得到广泛研究。相比卫星等导航方式1，基于机器视觉【2_4】的田间车辆自动导航、辅助导航，获得信息量丰富、使用灵活、体积小、能耗低，被广泛采用。作物行提取是基于机器视觉的田间导航实现的关键步骤，而杂草一直是作物行提取中难于去除的干扰，目前主要有两大类解决方法：一类，在作物行提取时去除杂草：一类，在作物行提取前去除杂草。第1类通过作物行提取算法，利用杂草分布的随机性，去除杂草干扰。如采用正、余弦模版匹配的方法【5，认为农田图像中，作物行间像素值分布和正余弦负半波匹配，但当杂草较密集时这种匹配将被打破。还有些研究者利用直线提取算法

10、的统计特性抗杂草干扰，其中Hough变换算法应用最为广泛。文献6】采用固定阈值法将采集的灰度图像二值化，万方数据电子与信息学报 第38卷然后将作物行看成长方形的盒子，盒子用若干条相邻并相交于虚交点的直线表示，通过Hough变换求取所有直线对应峰值的最大值来确定导航线位置。文献f71提出基于一点的改进Hough变换检测导航线，一定程度上克服了原始Hough变换计算量大，峰值难确定的缺点。文献f8采用子图像叠加方法减少待处理像素数，然后通过Hough变换，获得导航参数。文献9】提出了基于Hough变换和Fisher准则的多垄识别统一模型。文献101采用基于梯度的随机Hough变换方法，降低算法复杂

11、度的同时提高了算法实时性。在认为作物行间距已知的前提下，文献11禾U用TheilSen回归算法，提高作物行提取过程抗杂草干扰能力。最小二乘法112,1a1和Hough变换相比实时性要好，但提取作物行时对杂草更敏感。这类利用算法的统计性，杂草分布随机性及农作物按作物行排列的规律性去除杂草干扰的方法，当杂草较密集，尤其分布较均匀时杂草去除能力大大减弱。第2类，利用作物和杂草在颜色、外观形态等方面的差异去除杂草。如文献f141试图采用两次最大类间方差f简称Otsu)分割法，来体现杂草，作物颜色上的小差异，实现杂草的去除，但当杂草和作物颜色不可区分，或光照不匀明显，或农田中不仅含有一种杂草时，杂草去除

12、效果大打折扣；文献151尝试依据杂草和作物产生主频的差异，实现杂草的去除，但当杂草和作物主频不可区分时效果变差；文献16在认为行与行间距都相等的前提下，利用感兴趣区域像素值求和的变化体现作物和杂草的区别，当杂草分布均匀、密集时处理效果不理想，且行间距相等的限定使该方法仅能应用于强结构化，田间各处作物生长态势都相当的情况；文献17$1J用作物和杂草的高度差异去除杂草，但杂草生长是随机的，不能保证杂草在任何情况下都比作物矮。为解决上述问题，本文给出一种基于小波变换及Otsu分割的农田作物行提取方法，充分利用杂草、作物所产生颜色、频率信息的区别，并充分考虑作物按作物行交替的结构特点，来有效提高作物行

13、提取算法的抗杂草干扰能力，该方法对作物行间距是否变化、光照是否均匀、作物生长的高矮、农田所含杂草种类数都不做要求。2非绿色噪声去除拍摄俯仰角45。，相机距地面高度17 in左右，采集分辨率480640像素，天气状况为阴天时，获得的萝卜地彩色农田图像如图lfa)所示。分别采用如式(1)所示的三基色平均法及如式(2)所示的超绿变换，转换为灰度图像，转换成的灰度图，分别如图l(b)和图l(c)所示。I=(R+G+B)3 (1)I=2GRB (2)式中的R指彩色图像中红色分量灰度值，G指绿色分量灰度值，而B指蓝色分量灰度值，表示转换成的灰度图像的亮度值。式(1)在进行灰度变换过程中，各基色分量R，G，

14、B所占比重一样，保持了原农田图像的明暗、颜色深浅变化；采用式(2)进行彩色图像灰度变换，能有效提高绿色植物和其它非绿色噪声的对比度，有利于非绿色噪声的去除。对图1fc)所示灰度图像进行Otsu分割，获得二值图像，提取二值图像中像素值为零的坐标，将图1(b)相同坐标位置处的像素值置零，得到灰度图像，如图2所示。图2不仅有效去除了农田中的非绿色噪声，还很好地保留了作物、杂草等绿色植物之间的明暗变化，为杂草的去除提供可能性。3基于Otsu及小波变换的绿色噪声去除采用Otsu方法对图2中保留下来的绿色植物进行分割，获得图3(a)所示二值图像。当作物和杂草颜色可区分性较强时，通过这步操作可将田间杂草(a

15、)彩色农出劁像 (1)j犟色平均法铁得的灰度图像 (c)超绿变换获得的荻度图像图1农田彩色图像到灰度图像的转换万方数据第1期 韩永华等：基于小波变换及Otsu分割的农田作物行提取 65刳2去除限绿色噪声的农灰度图像图3频率总量最多的小波分解层为非作物层的示例图去除14】，否则，也能提高密集植物的频率辨识度。这是因为密集植物容易发生叶片相互遮挡，遮挡处的颜色和未发生遮挡的颜色略有不同，当作物和杂草颜色不可分时，通过Otsu分割能将遮挡、未遮挡区域的颜色差异体现出来，从而将因遮挡发生的重叠叶片分离开，使它们成为独立体，从而能更好地通过后续的小波变换体现出叶片主频。对图3(a)进行6级Haar小波分

16、解，将分解后频率最高的细节小波分解层，称为第1级小波分解层，随着频率的降低依次将分解后的细节层，称为26级小波分解层。按式(3)依次对26级小波分解层进行计算，计算结果称为频率总量。q=瓯+吼 (3)式中f为小波多分辨率分解的第f层，h为f层水平分解系数图像，d为f层对角分解系数图像，v为f层垂直分解系数图像，瓯称为f层水平分解系数图像的频率，瓯称为f层对角分解系数图像的频率，Gh称为l层垂直分解系数图像的频率，式(3)中的c，(o=h，d，v)，由式(4)获得。=丽b善M若N尼扩式中MN为小波多分辨率分解第f层各方向细节子图像的大小，L(i，歹)为坐标点(i，J)处的小波多分辨率分解系数。将

17、计算完的频率总量按从大到小顺序排列，取频率总量最多的前3级，验证有无作物行交替，选择有作物行交替的频率总量最多的小波分解层作为作物层，进行后续分割、作物行提取操作。选频率总量最多的前3级主要是考虑到，虽多数情况下农田作物占有数量上的优势，相比其它物体会在对应小波分解层上产生更多的频率信息，但杂草特别密集时，其产生频率总量可能超过作物，或较密集杂草、作物产生的非主频信息叠加在同一小波分解层上，也可能使产生的频率总量超过作物，同时考虑到一些不可预知因素，并综合考虑算法执行实时性，选频率总量最多的前3级作为研究对象。此外小波分解的第1级主要是各种边缘信息，不能有效表达作物特征，所以计算频率总量时不考

18、虑。针对图3(a)计算的各小波分解细节层频率总量值如表l所示。表1各级小波分解层频率总量表由表1可以看出，频率总量最多的前3级依次为6，4，5级。先对频率总量最多的细节子图像进行小波重构，重构后图像如图3(b)所示。接下来就是对重构小波图像进行有无作物行交替的判定，判定过程包括两步：进行作物行交替判定的图像行的选定，针对选定图像行进行的作物行交替判定。实现流程分别如图4(a)和图4(b)所示。图4(a)中，(1)第1步对图像行像素求和生成向量日，其中参与求和的是行像素的像素值，日1为行向量，大小和图像高度(以像素为单位)相同；(2)流程图上侧分支第2步为获取向量日1极值点，以极值点所在图像行作

19、为作物行交替判定的对象，不仅能捕捉杂草区域，且能尽量避免断垄行的入选；f3)考虑到第5、第6级小波分解级重构后最小图像块大小为1616像素或3232像素，会引起行向量日1出现一系列紧邻、且数值相同或相近极值点，而这一系列极值点代表同一组图像块，因此需进行流程图上侧分支第3步，获取典型极值点的操作，即仅保留一系列紧邻极值点中位置在中间或近似中问的极值点；f4)流程图上侧分支第4步，对极值点数量进行判断，如未达到规定数量繇则进行流程图第5步，即补充极值点。考虑到典型极值点数量少，兼顾均距选取图像行，同时为避免所有极值点都选在农田结构自身具有的杂草区域，如长杂草的田头，万方数据66 电子与信息学报

20、第38卷重构的是小波分解的第5级或第6级细节子图像 锗H麓紧的极值点l 7l 极值点I补充非极值ll皇望坚笪皇I灭极值点数量是否达到要求值品第3或第4级细节子图像 一(a)进行作物行交替判定的图像行的选定计算选定图lI以蜀为阈值I 统计图像行上像行像素值卜I对图像行进卜|长度大丁等丁q的平均值l l行二值分割l l的黑色线段数量否否每份极值点数量是否达到要求值日纠；臀嬲擎刊鎏耋l雠 是 +农田结构自身引起某图像行不满足黑色线段数量递减规律?fb)针对选定图像行进行的作物行交替判定图4作物行交替判定流程图选定补充过程为：寻找相邻距离最大的两个极值点，在这两个极值点中问位置附近选取非零的非极值点补

21、充为极值点，直到极值点数量达到要求。小波分解的第2到第4级细节子图像重构后，生成的向量日1极值点较多，为减少计算量，先进行极值点筛选，图4(a)下侧分支第2步求行向量日1的均值，到第4步记录甄中保留下来的值的极值点就是用于完成这步操作；图4(a)下侧分支第5步将行向量日等分为付份，每份取S个极值点，是为了避免极值点过于集中；极值点不足S时，需通过图4(a)下侧分支第7步补充极值点：考虑到平行直线透视成像时，会在远方交于一点，导致图像远端作物行交替特征相比图像底端要弱，尤其俯仰角小时这种削弱更明显，因此优先在靠近向量甄起始位置且未补充过极值点的份中补充，且兼顾均匀取点原则，限制每份获取的总极值点

22、个数小于等于2S个。图4fb)所示作物行交替判定第1步，即计算选定图像行像素值的平均值蜀，到第2步，即以写为阈值对选定图像行进行二值分割的阂值选定方式能避免分割时不同图像行间的相互干扰，减少杂草行影响范围；第3步中，保留像素数大于等于选定阈值互的黑色线段，是为了剔除作物叶片间大的间隙古l选黑色线段数量最大值与最小_1H值差异较小且满足递减规律较I 多的小波重构层作为作物层及垄间杂草到作物之间的间隙生成的短黑色线段；第4步中的递减指黑色线段数量依次递减或和前面相邻的相同，依据图像顶到底的黑色线段数量是否满足递减来判定有无作物行交替，是考虑到：透视投影下，图像远端收入的作物行多于图像底端。图4(b

23、1所示流程图下侧第2个判别框“农田结构自身引起某图像行不满足黑色线段数量递减规律”，是指通过极值点选择仍未能避免的田间断垄、或图像中出现的田头，引起的黑色线段递减规律破坏，判别方法为：将不符合递减要求的行的位置记录下来，到其它两幅重构图上验证对应行的黑色线段数量是否也不满足当幅图像的递减规律，如不满足，则认为农田结构自身造成了这种现象。上述参数佗和参数S的选定原则为：图像采集俯仰角小时，佗选得大，S选得小，保证均匀取极值点的同时，避免出现极值点集中于图像远端；俯仰角大时，佗选得小，S选得大，此时图像远端作物行交替仍很明显，选极值点时不需考虑避开这一端，而主要考虑避开断垄的影响：份数铭取得小，能

24、减少单份中某作物行完全是断垄的可能性，从而减少极值点取在断垄处的概率。当采集农田图像俯仰角为00一400时，站选10，S选1；40。一90。时，嗨万方数据第1期 韩永华等：基于小波变换及Otsu分割的农田作物行提取 67n选5，S，选2。考虑目前田问导航常用的小俯仰角一般在20。一30。左右f5,11】，而俯仰角小时，尤其窄行作物，会有更多作物行包含进图像，从而在图像行上产生较多的黑色线段，此时以像素为单位的黑色线段长度就会变短，选定TI为图像宽度的115。通过验证，图3(b)对应行上黑色线段的数量不满足递减规律，而重构后小波分解频率总量次多的第4级重构图如图5(a)所示，满足黑色线段数量递减

25、规律，选为作物层。(a)4、波分解第4级重构图 (b)0tsu分割后的二值图图5基于Otsu及小波变换的噪声去除过程由原图1fa)可以看出，除了田头有杂草，作物田间并无杂草，理论上频率总量最多的层应有作物行交替，而事实相反，是因为萝卜地作物叶片大且密，虽在第6级大尺度辨识的小波分解级中产生较多频率，但行间距窄，有些地方甚至封垄了，导致行间产生的频率和作物叶片相似，于是在第6级大尺度辨识的小波分解级中作物行间被掩盖，其黑色线段数量的变化规律为：510(2)，480(3)，450(3)，420(2)，215(4)，185(4)，120(4)，95(3)，25(4)，5(4)。小括号外的数字表示选取

26、图像行的位置，小括号内数字表示在对应行上检测到的黑色线段数量。对图5(a)所示作物层再次进行Otsu分割，去除这一层中的弱频率信息，使作物层表达的主频信息更集中，有利于进一步减少杂草干扰。当频率总量最多的前3级包含小波分解的第5级、第6级中的至少一个时，表明作物叶片较大或较密集，或兼而有之，则在更高分辨率小波分解级上保留的较高频信息更容易相互成为邻域，因此去除更高分辨率小波分解级中1个像素大小的孤立像素，不仅不会过多减少作物信息还会提高噪声去除率，但这种噪声去除操作不需对每级高分辨率小波分解级进行，仅需对重构后获得的作物层进行，此时1个像素大小的滤波操作变为：作物层为小波分解的第4级，则滤除作

27、物层Otsu分割图上等于小于8 X 8像素大小的孤立白色团块；作物层为小波分解的第3级，则滤除作物层Otsu分割图上等于小于4 X 4像素大小的孤立白色团块；作物层为小波分解的第2级，则滤除作物层Otsu分割图上等于小于2X2像素大小的孤立白色图块。对于黑色孤立团块，无论大小，在作物层中全部去除，即将其变为白色。这是因为黑色孤立团块至少4邻域为白色团块，此时黑色孤立团块如在作物行上，将其变为白色前景后，能避免将其误认为作物行交替时的黑色线段；如果黑色孤立团块是作物行间白色噪声围绕产生的，将其去除有利于做出白色噪声所在位置不存在作物行交替的判定。经过上述处理获得的作物层二值图像，如图5(b)所示

28、。为了更明显看出杂草去除效果，将图5(b)的像素值按行求和，投影成曲线，如图6(a)所示。图6(a)中，水平坐标表示所采集农田图像的第n行，针对图5(b)，0冬佗639，纵坐标表示每行像素值之和，单位为个，对于像素值要么为1，要么为0的二值图像，像素值之和也表示了每行前景像素个数。将投影曲线产生明显下跳沿，且下跳沿之后行像素值之和接近0的跳变点找出来，将跳变点之后所有行对应的像素点坐标记录下来。对图2所示图像进行Otsu二值分割，将分割后图像中和记录的像素点坐标对应的像素点的像素值置零，获得图像如图6(b)所示。由图6(b)看出杂草被有效去除。在杂草等噪声去除干净的前提下，对作物行检测算法的抗

29、噪要求大大降低，可选用易受杂草干扰，但实时性好的最小二乘法检测作物行直线，检测结果如图7所示。图中虚白线为由检测到的作物行直线得到的导航线。4实验结果与分析为验证算法去除杂草能力，对不同光照条件采集的含杂草的100余幅农田图像(其中密集杂草图片10 200 310 400 500 600I矧像行(a)行像索和投影l刘 (1，)去除尔单j疗的值分割陶蚓6朵草去除能力展示图万方数据电子与信息学报 第38卷图7导航线检测结果约占40，断垄图片约占15，约79图像包含田头，大小采集俯仰角各约占50)进行了算法验证。这些图像的大小有640480像素和480640像素两种，图1(a)所示480640像素大

30、小图像相比图8(a)所示640480像素大小图像能包含更多作物行向信息，在作物行检测过程中有利于减少断垄影响。实验结果表明杂草去除率达到99，仅当杂草外观形态及颜色和作物极其接近时，杂草去除效果才会变差。这里杂草去除率指所选100余幅农田图像，经过处理，能正确获得作物行的图像所占样本的百分比；而颜色接近并不指为同种颜色，仅当杂草和作物同为绿色，且绿的深浅都不可区分时才说它们颜色相近。同时因为小波分解图像在重构时，体现光照不匀的近似级未参与重构，作物行提取过程抗光照不匀能力得到提高。41分割结果分析图8(a)为天气情况为晴天的正午采集的640x480像素大小的包含密集杂草的青菜地图像，图8(b)

31、为采用超绿变换和Otsu分割相结合，去除三基色平均法获得的灰度图像中的非绿色噪声后获得的灰度图像，图8(c)为在图8(b)所示图像的植物间再次进行Otsu分割获得的灰度图像，图8(d)是经小波分解后对选定的小波分解级，此处为小波分解的第5级进行重构获得的作物层图像，图8(e)则为对图8(d)进行Otsu分割后获得的作物二值图，去除8(e)中的黑色孤立团块后获得图8(f)，记录图8(f)中像素值为零的点的坐标，将图8(b)相同坐标处的像素值置零，获得的图像如图8(g)所示。图8(g)的存在用于直观展示杂草去除效果。从图8(c)可以看出，对图8(a)这种杂草的绿和作物的绿有较明显差异，且农田图像光

32、照比较均匀、没有阴影的图像，进行两次Otsu分割，的确就像文献141所指出的能很好地去除杂草影响(经两次Otsu分割后获得的图8(c)虽仍然保留了部分杂草信息，但保留下来的多为独立的小团块，很容易去除1，但针对图1(a)进行两次Otsu分割获得的图3(a)却表明文献f141的方法失效了，导致这一结果的原因主要是图l(a)中杂草颜色和农田作物很接近，又加上图像采集光线暗淡导致颜色分辨率降低，最终使得依靠两次Otsu分割的杂草去除不再有效。针对上述100余幅农田图像，表2给出了杂草去除效果较好的Hough变换、两次Otsu分割算法【川、主要基于频率信息的方法【15和本文算法杂草去除率的对比。表2不

33、同算法杂草去除率比较()I+ i，!i捌i，。j酿 tjj芷。二j：孓：；0鬈t，；枣刮隧釜I(c)()1分割后获得的作物_值图(f)去除黑色独祝块的作物_值罔 (g)杂草去除效果图图8青菜地图像分割过程图万方数据第1期 韩永华等：基于小波变换及Otsu分割的农田作物行提取 69表2中的Hough变换指对图像进行Otsu二值分割后再进行Hough变换提取导航线；两次Otsu算法指采用文献141提出的两次Otsu分割法去除杂草，再采用最小二乘法拟合直线的过程，如牺牲算法实时性，将最小二乘法用Hough变换代替，则两次Otsu算法的杂草去除率可提高到87；基于频率信息的方法指依据杂草和作物产生主频

34、的差异，结合作物行的交替去除杂草，再采用图像旋转投影法15】提取导航线的过程。基于频率信息的方法和本文方法相比主要差别在于：(11在小波变换前及找到作物层后都未对图像进行基于颜色信息的分割操作，使产生主频和作物相近的杂草难于去除；f2)缺少依据作物浓密、叶片大小去除作物层噪声的过程，影响后续作物行提取的准确度；(3)其通过扫描小波重构层全部图像行的小波系数，依据作物行上小波系数递增、递减变换的规律性进行作物行交替判定，相比本文方法，其对小俯仰角采集图像引起的作物行宽变化及断垄敏感。42时间复杂度分析本文所提出的农田作物行提取流程，主要包括基于超绿变换和三基色平均法的灰度变换、3次Otsu分割、

35、小波分解及重构，最小二乘法导航线提取。相比一般作物行提取过程所包含的一次彩色图像到灰度图像变换、一次基于作物颜色信息的分割运算，及作物行提取运算，多出了一次灰度变换、两次基于颜色的分割运算及小波变换，但文中采用的灰度变换是目前农田车辆导航研究中常用灰度变换方法中实时性最强的，最费时的作物行提取环节因为前面分割部分噪声去除得好，可以选择实时性最好的常用算法，而较少考虑其鲁棒性。小波重构次数和作物行交替判定过程相关，最理想结果是在频率总量最多的小波分解层上就能判定存在作物行交替，最差结果是频率总量前三的小波分解层全部重构后才能判定，针对100多幅采集的农田图像实验，发现只需重构一级就能检测到作物行

36、交替的比例接近64，需重构两级的约22，需重构3级的仅约14。如果用于实验的农田图像为宏观结构化更好的机耕田，贝0fJ,波重构的总体次数将会更少。工控机采用Intel(R)Core(TM)2处理器，主频293GHz，4 G内存，运行系统为Windows XP，针对480640像素大小的图像，各算法执行的平均时间为超绿变换41 ins，三基色平均法灰度变换65 ms，3次Otsu分割332 ms，小波变换421 ms，最小二乘法46 ms。处理单幅图像平均用时约132 ms，达到导航实时性要求。43作物行提取过程实时性的进一步提高本文提出的作物行提取过程，可增加有无杂草检测环节。以图1(a)所示

37、的萝卜地图像为例，其去除背景噪声后获得了图2，计算图2行像素值的和，投影成曲线，如图9所示。坐标系横坐标表示行，纵坐标表示行像素的灰度值之和，针对图9量级为10a。对整条曲线进行扫描，获得明显跳变沿，将图像底端到距离图像底端最近的明显跳变沿间的图像作为处理对象，检测有无作物行交替；如未检测到明显跳变沿，则在整幅图像验证有无作物行交替，如有作物行交替，说明农田中无杂草。羁5蛩a纛s差z蔷O100 200 300 400 500 600图像行图9去除背景噪声的萝h地灰度图像的行像素值的和的投影图选取图像底端到明显跳变沿处的图像进行作物行交替判定，主要是考虑到：田问导航时，田头是从图像顶端渐入图像的

38、，一旦田头进入图像，只能依据图像底端的作物行进行田间作物行提取；明显跳变定义成跳变沿前后两部分行像素值和的均值相差两倍以上，该定义能避免作物自身分布或生长不匀造成的小跳变被当成田头产生的跳变；如果因为跳变选择两倍以上导致田头漏判，或者错判，会导致检测不到作物行交替，从而判定农田中有杂草，此时会进行后续小波变换等处理，不会影响作物行的正确提取。当待处理农田图像中无杂草或极少杂草时，增加的杂草有无判定操作能节省几十毫秒时间；当待处理农田图像有较多杂草时，增加的这步操作会增加不到3 mS处理时间。因文中想着重用图1(a)体现所提出的作物行检测过程的抗杂草干扰能力，因此跳过了这步操作。5结论在农田车辆

39、视觉导航图像处理过程中，杂草尤其高密度杂草干扰的去除是一大难题，经过上述超绿变换、Otsu分割、小波变换、频率总量判定的相互配合应用，不仅有效去除了杂草干扰，且通过在小波重构过程中去除近似分解层系数减少了光照不匀、阴影对作物行提取的影响，和其它作物行提取过程相比，进一步提高了作物行提取的鲁棒性。万方数据70 电子与信息学报 第38卷138】9参考文献张柏华，马红光，孙新利，等基于正交约束的导航接收机空时自适应方法J】，电子与信息学报，2015，37(4)：900906doi：101199JEITl40740ZHANG Baihua，MA Hongguang，SUN Xinli，et a1Spa

40、cetime adaptive processing technique based on orthogonalconstraint in navigation receiverJJournal of Electronics&Information Technology，2015，37(4)：900906doi：101199JEITl40740姬长英，周俊农业机械导航技术发展分析fJ，农业机械学报2014，45(9)：4454JI Changying and ZHOU JunTechnical analysis of thedevelopment of agricultural machine

41、ry navigationJTransactions o|the Chinese Society ofAgricultural Machinery2014，45(9)：4454李骏扬，金立左，费树岷，等基于多尺度特征表示的城市道路检测J】电子与信息学报，2014，36(11)：25782585doi：103724SPJ1146201400271LI Junyang，JIN Lizuo，FEI Shumin，et a1Urban roaddetection based on multi-scale feature representationJ，Journal of ElectronicsIno

42、rmation Technology，2014，36(11)：25782585doi：103724SPJ1146201400271李盛辉，田光兆，姬长英，等，自主导航农业车辆的全景视觉多运动目标识别跟踪【J农业机械学报，2015：46(1)：卜7，LU Shenhui，TIAN Guangzhao，JI Changying，et a1Multiplemoving objects tracking based on panoramic vision forautonomous navigation of agricultural vehicleJTransactions ofthe Chines

43、e Society ofAgricultural Machinery2015，46(1)：卜7KISE M， ZHANG Q， and ROVIRA M F， Astereovision-based crop rOW detection method fortractor-automated guidanceJBiosystems Engineering，200590(4)：357 367ASTRAND B and BAERVELDT A JA vision basedrowfollowing system for agricultural field machineryJMechatroni

44、es，2005，15(2)：251269Leemans V and Destain M F，A computervision basedprecision seed drill guidance assistanceJComputers andElectronics饥Agriculture，2007，59(12)：卜12BAKKER T，WOUTERS H，ASSELT K V，et a1A visionbased row detection system for sugar beetJComputers andElectronics in Agriculture，2008，60(1)：879

45、5张志斌，罗锡文，周学成，等基于Hough变换和Fisher准则的垄线识别算法【J中国图象图形学报，2007，12(12)：21642168ZHANG Zhibin，LUO Xiwen，ZHOU Xnecheng，et a1Croprows detection based Oil Hough transform and fisherdiscriminant criterion functionJJournal of Image andGraphics，2007，12(12)：21642168JI Ronghua and QI LijunA croprow detection algorithm

46、based on random hough transformationJMathematical andComputer Modelling，2011，54(34)：10161020GUERRERO J M，GUIJARR0 M，MONTALVO M，eoLAutomatic expert system based on images for accuracy croprow detection in maize fieldsJExpert Systems withApplications，2013，40(2)：656664姜国权，王志衡，赵翠君基于已知点的作物行检测方法J】应用基础与工程科

47、学学报，2013，21(5)：983990JIANG Guoquan，WANG Zhiheng、and ZHAO CuOunAnalgorithm of detecting crop rows based on knownpointsJJournal of Basic Science and Engineering，2013，21(5)：983990胡炼，罗锡文，张智刚，等株间除草装置横向偏移量识别与作物行跟踪控制J农业工程学报，2013，29(14)：8-14HU Lion，LUO Xiwen，ZHANG Zhigang，et a1Sideshiftoffset identification

48、 and control of crop row tracking forintra-row mechanical weedingJTransactions of the ChineseSociety of Agricultural Engineering，2013，29(14)：8-14MONTA【O M，PAJARES G，GUERRERO M，et。tAutomatic detection of crop rows in maize fields with highweeds pressureJExpert Systems with Applications，2012，39(15)：1188911897韩永华，汪亚明，康锋，等基于小波多分辨率分解的农田障碍物检测J农业机械学报，2013，44(6)：215221HAN Yonghua，WANG Yarning，KANG Feng，et a1Detection of obstacles in farmland based on wavelet multiresolution transformJTransactions of the Chinese Societyo，Agricultural Machinery，2013，44(6)：215221JIANG G QWANG Z Hand LI