CN105046229A - 一种农作物行的识别方法及装置 - Google Patents
一种农作物行的识别方法及装置 Download PDFInfo
- Publication number
- CN105046229A CN105046229A CN201510446384.4A CN201510446384A CN105046229A CN 105046229 A CN105046229 A CN 105046229A CN 201510446384 A CN201510446384 A CN 201510446384A CN 105046229 A CN105046229 A CN 105046229A
- Authority
- CN
- China
- Prior art keywords
- crop row
- unique point
- crop
- straight
- crops
- 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
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000012417 linear regression Methods 0.000 claims abstract description 30
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 239000000284 extract Substances 0.000 claims description 4
- 241000196324 Embryophyta Species 0.000 description 7
- 230000011218 segmentation Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 240000008042 Zea mays Species 0.000 description 3
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 3
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 3
- 235000005822 corn Nutrition 0.000 description 3
- 238000005286 illumination Methods 0.000 description 2
- 238000009333 weeding Methods 0.000 description 2
- 235000017060 Arachis glabrata Nutrition 0.000 description 1
- 244000105624 Arachis hypogaea Species 0.000 description 1
- 235000010777 Arachis hypogaea Nutrition 0.000 description 1
- 235000018262 Arachis monticola Nutrition 0.000 description 1
- 235000016068 Berberis vulgaris Nutrition 0.000 description 1
- 241000335053 Beta vulgaris Species 0.000 description 1
- 229920000742 Cotton Polymers 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 244000068988 Glycine max Species 0.000 description 1
- 235000010469 Glycine max Nutrition 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 240000000111 Saccharum officinarum Species 0.000 description 1
- 235000007201 Saccharum officinarum Nutrition 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 244000098338 Triticum aestivum Species 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 235000020232 peanut Nutrition 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 238000013138 pruning Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种农作物行的识别方法及装置,该方法首先通过模糊聚类方法将采集的农作物图像转化为二值图像,然后从获取的二值图像中提取感兴趣的区域ROI图像,通过横条法获取其中农作物特征点,并对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程,最后再通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。本发明的装置包括模糊聚类模块、特征点提取模块和拟合模块。本发明的方法及装置,识别农作物行准确度高,运行速度块,抗干扰能力强。
Description
技术领域
本发明属于农作物栽培技术领域,尤其涉及一种基于图像处理的农作物行的识别方法及装置。
背景技术
我国幅员辽阔,地形、气候类型复杂多样,从南到北分为热带、亚热带、温带和寒带,主要粮食作物有水稻、小麦、玉米、大豆等,经济作物有棉花、花生、油菜、甘蔗和甜菜等。然而我国人口众多,耕地面积相对较少,因此农业特别是种植业在我国的地位非常重要,关系整个国计民生。随着农业机械化和信息化的迅速发展,实现农业自动化的要求越来越迫切。
机器视觉配合大中型农业机械在农业领域的运用也越来越广泛,尤其是在视觉导航和作物识别方面,准确性和费用都得到了较大的改善。因此在农作物的自动除草、收获、施肥、修剪、耕耘等工作中,基于图像处理进行农作物行的识别显得尤为重要。
现有的农作物行的识别方法主要是利用超绿法、最大类间方差法分割图像,通过Hough变换识别农作物行。但是这些方法图像分割不准确,不能很好的区分农作物和杂草,而且计算量大,达不到实时性的要求。尤其在农作物中有大量杂草存在的情况下,不能得到理想的结果。
发明内容
本发明的目的是提供一种农作物行的识别方法及装置,以避免现有技术图像分割不准确,识别效率不高的技术问题。
为了实现上述目的,本发明技术方案如下:
一种农作物行的识别方法,所述识别方法包括:
通过模糊聚类方法将采集的农作物图像转化为二值图像;
从获取的二值图像中提取感兴趣的区域ROI图像,通过横条法获取其中农作物特征点;
对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程;
通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。
优选地,所述通过模糊聚类方法将采集的农作物图像转化为二值图像,是以采集的农作物图像作为整个样本,以像素点G通道的值所占的百分比作为样本元素,将农作物和背景的聚类中心分别初始化为0.35~0.40和0.30~0.35进行模糊聚类得到。
本发明以农作物图像作为整个样本,可以直接以农作物图像像素点的像素R\G\B三通道的值作为样本元素来进行模糊聚类。但是优选地,是采用像素点G通道的值所占的百分比作为样本元素,从而降低了聚类的维数,将农作物和背景的聚类中心分别初始化为0.35~0.40和0.30~0.35进行模糊聚类得到二值图像,提高了聚类速度。
进一步地,所述采集的农作物图像的宽为W像素,高为H像素,所述ROI图像的宽为w=W/2像素,宽为h=H/2,所述通过横条法获取其中农作物特征点,包括:
将ROI图像分为Q条宽度相同的横条,用Sp,q表示第q条横条第p列中白色像素出现的次数,其中p的取值从1到w,w为ROI图像的宽度像素;
对于第q条横条,对应有阈值uq,阈值uq为第q条横条中所有Sp,q的均值;
当Sp,q小于等于uq且Sp+1,q大于uq时,认为进入农作物行,记录此时的列坐标为p1;
当Sp,q大于等于uq且Sp+1,q小于uq时,认为离开农作物行,记录此时的列坐标为p2;
计算进入和离开农作物行时列数的差值Δ=p2-p1,如果Δ大于设定的常数d,则认为横条q上从第p1到p2的段为农作物,并取该段中点为农作物特征点;
遍历所有横条,获取ROI图像中所有的农作物特征点;
其中,所述常数d的取值范围为:W/20<d<W/15。
进一步地,所述对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程,包括:
根据特征点的分布,将特征点划分为不同的农作物行;
对于任一农作物行,设农作物行直线方程为:
y=kx+b
其中b为斜距,k为斜率,计算属于该农作物行的所有特征点到该直线的距离l:
进一步计算,属于该农作物行的所有特征点到该直线的距离平方和l′:
对上式的分子对k和b求偏导并令其为0,得到:
其中M为属于该农作物行的所有特征点的数量,r属于1~M,第r个特征点的坐标为(xr,yr),lr为第r个特征点与农作物行直线的距离,求解出上式中k和b的解分别为和将其带入直线方程得到农作物直线方程为:
进一步地,所述通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正,包括:
根据农作物行直线方程,计算属于该农作物行的特征点到该农作物行的距离,剔除距离大于设定的常数的特征点;
在剔除错误的特征点后,根据剩下的特征点再次拟合出农作物行直线方程,并再次根据新拟合出的农作物行直线方程,计算特征点到该直线的距离,剔除距离大于设定的常数的特征点;
如此循环,直到达到最大线性回归次数,或当线性回归剔除的特征点数量为0时迭代停止。
本发明还提出了一种农作物行的识别装置,所述装置包括:
模糊聚类模块,用于通过模糊聚类方法将采集的农作物图像转化为二值图像;
特征点提取模块,用于从获取的二值图像中提取感兴趣的区域ROI图像,通过横条法获取其中农作物特征点;
拟合模块,用于对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程,并通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。
本发明提出的一种农作物行的识别方法及装置,通过模糊聚类对采集图像进行分割,利用横条法确定特征点,并利用特征点在x轴上的位置判断所属农作物行,利用线性回归求出过渡农作物行,在通过多次线性回归剔除错误的特征点,得出最终的农作物行。本发明识别农作物行准确度高,运行速度块,抗干扰能力强。
附图说明
图1为本发明农作物行识别方法的流程图;
图2为本发明实施例二值图像示意图;
图3为本发明实施例ROI图像示意图;
图4为本发明实施例ROI图像中特征点示意图。
具体实施方式
下面结合附图和实施例对本发明技术方案做进一步详细说明,以下实施例不构成对本发明的限定。
本发明利的总体思路是用计算机视觉技术,对图像采集设备采集到的农作物行图像进行分析处理,来识别作物行。本实施例以玉米幼苗为例来进行说明。
如图1所示,一种农作物行的识别方法,包括如下步骤:
步骤S1、通过模糊聚类方法将采集的农作物图像转化为二值图像。
本实施例用X来表示所有样本,即采集的农作物图像的每个像素点的RBG通道的值作为一个样本元素xi,xi对应于一个像素点。将整个样本X分为农作物和背景两部分,所以聚类数c=2,即将图像中的像素点分为农作物和背景两类。
用fij表示xi属于第j类的隶属度,用vj表示第j类的聚类中心,对样本X中的每个xi进行迭代,所有样本元素计算完成依次为依次迭代,通过迭代计算出农作物和背景的聚类中心。
在迭代之前需要对聚类中心vj进行初始化。由于农作物部分在理想情况下应该是绿色,本实施例设置农作物聚类中心v1={0,255,0};剩下的背景中绿色通道相较于红色通道和蓝色通道,不需要占主导地位,设置v2={255,0,128}。
隶属度fij和和聚类中心vj的迭代公式如下所示:
其中,dij是xi到聚类中心vj的欧氏距离,λ被称作指数权重,λ>1,n为样本X的容量。
该过程反复迭代直到||fij(t+1)-fij(t)||<ε或者已经到达指定的迭代次数tmax。
迭代结束后将得到两个聚合中心:v1={V1R,V1G,V1B}和v2={V2R,V2G,V2B},由于v1是农作物的聚类中心,本实施例根据v1来计算分割阈值,分割预置=V1G/(V1R+V1G+V1B)。
对于任意像素点xi,当其RGB三通道的值满足G/(R+G+B)大于分割阈值时,将该点像素置为255,即白色,否则置为0,即黑色,从而将原图像分割为二值图像。例如图2为根据采集图像得到的一副二值图像,其中白色像素代表了农作物,黑色像素代表了背景。
在不同的情况下,通过聚类分割可以很好的提取农作物信息,相比与常用的超绿法和大津法的组合的方法相比,模糊聚类可以保存更多的细节,比如杂草,落叶等,且在不同光照条件下,都能准确的将农作物和背景区分开。
值得注意的是,在聚类结束后,对于每个聚类中心Vj,其G通道的值所占的百分比val为:
通过大量的实验数据,如表1所示:
背景v1(%) | 作物v2(%) |
0.345 | 0.403 |
0.335 | 0.347 |
0.345 | 0.386 |
0.336 | 0.353 |
0.339 | 0.412 |
0.345 | 0.386 |
表1
不难发现农作物聚类中心对应的百分比val大多在0.35和0.40之间,背景聚类中心对应的百分比val在0.30和0.35之间。而在实际聚类过程中,起决定性作用的也是G通道的值所占的百分比。为了提高计算效率,优选地,直接以G通道的值所占的百分比作为样本元素,因此样本的维数从3降到1,并在每次聚类开始前,将农作物和背景的聚类中心分别初始化为0.35~0.40和0.30~0.35,避免每次聚类迭代时的重复计算,提高聚类速度。
步骤S2、从获取的二值图像中提取感兴趣的区域ROI,通过横条法获取其中农作物特征点。
首先,从获取的二值图像中提取感兴趣的区域ROI(RegionOfInterest),对于宽为W像素,高为H像素的采集图像,一般提取的ROI的长为w=W/2,宽为h=H/2,本实施例提取的ROI如图3所示,提取的ROI至少包括一条农作物行,本实施例包括两个农作物行。
对ROI使用横条法,将ROI图像分为Q条宽度相同的横条,对于宽为w像素,高为h像素的ROI图像,将其分为Q条横条,用Sp,q表示第q条横条第p列中白色像素出现的次数,其中p的取值从1到w。
对于第q条横条,对应有阈值uq,阈值uq为第q条横条中所有Sp,q的均值:
遍历Sp,q,用以下过程提取农作物行的特征点:
(1)当Sp,q小于等于uq且Sp+1,q大于uq时,说明进入农作物行,记录此时的列坐标为p1;
(2)当Sp,q大于等于uq且Sp+1,q小于uq时,说明离开农作物行,记录此时的列坐标为p2;
(3)每次进入和离开农作物行时,计算进入和离开农作物行时列数的差值Δ=p2-p1,如果Δ大于设定的常数d,则认为横条q上从第p1到p2的段为农作物,并取该段中点为农作物特征点。
其中,d为常数,可以理解为农作物行的宽度,其取值范围为:W/20<d<W/15。
通过遍历所有横条,可以得到若干特征点,如图4所示,其中黑点为特征点。
步骤S3、对所提取的农作物特征点进行线性回归,得出农作物行直线方程。
根据特征点的分布,将特征点划分为不同的农作物行。本实施例创建的ROI里一般仅保留两条农作物行,判断所有特征点的x轴坐标,如果坐标小于宽度的一半,将该点归于左侧的农作物行,否则归于右侧。
对于任一农作物行,根据属于该农作物行的所有特征点通过线性回归得到农作物行直线方程,过程如下:
假设特征点的坐标为(x,y),农作物行直线的方程为:
y=kx+b
其中,b为斜距,k为斜率。
则特征点到该农作物行直线的距离l为:
所有特征点到直线距离的平方和l′为:
其中M为特征点的数量。
可见,如果特征点(x,y)在该直线上,则l为0,但不可能所有特征点都恰好落在该直线上,通过求k和b使l′取最小值,即对的分子的k和b求偏导并令其为0:
其中M为属于该农作物行的所有特征点的数量,r属于1~M,第r个特征点的坐标为(xr,yr),lr为第r个特征点与农作物行直线的距离,解出上式的解和将和带入直线方程,得到农作物行直线的方程为:
使得特征点基本分布在该直线两侧附近。
步骤S4、通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。
在理想状况下,此时通过步骤S3中得出的农作物行直线方程已经可以得到较为准确的农作物行,但实际情况中ROI内未必只包含两条农作物行,而且光照,天气,杂草等外界因素也会影响准确性。为了得到准确的作物行,本实施例需要剔除错误的特征点,保留正确的特征点。
在所有外界环境中,ROI中其他行的侵入和杂草的影响最大,而这些侵入的农作物行和杂草在光谱上与目标接近,所以本实施例利用点与线的位置关系来判断哪些特征点需要剔除。
本实施例引入多次线性回归的方法,通过S3,已经得到了农作物行直线方程并通过点到直线距离公式计算特征点到该直线的距离l,当l大于W/15时,剔除该点,W为图像的宽。
在剔除错误的特征点后,根据剩下的特征点再次拟合出农作物行的直线方程,并再次根据新拟合出的农作物行的直线方程,计算特征点到该直线的距离l,当l大于W/15时,剔除该点。
如此循环,直到达到最大线性回归次数,本实施例为10次,或当线性回归剔除的特征点数量为0时迭代停止。
本实施例为了提高效率,对两条农作物行分开进行线性回归,当一条农作物行上的特征点没有错误时,直接跳过对另一条农作物行进行求解。
本实施例还提出了对应于上述方法的一种农作物行的识别装置,所述装置包括:
模糊聚类模块,用于通过模糊聚类方法将采集的农作物图像转化为二值图像;
特征点提取模块,用于从获取的二值图像中提取感兴趣的区域ROI图像,通过横条法获取其中农作物特征点;
拟合模块,用于对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程,并通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。
各模块所执行的操作与上述农作物行识别方法对应,这里不再赘述。
现有技术采用Hough变换和本实施例提出的多次线性回归的对比如表2所示:
表2
表2中的数值均为平均数值。从表中可以看出,无论是在准确度还是耗时,后者均优于前者,而且当图像逐渐增多时,Hough变换的计算量增大趋势也快于本实施例的方法。对玉米幼苗进行实验后的实验数据显示,在不同天气环境的条件下,识别农作物行的准确率达到96%,误差在2°左右,运行速度在10ms以内,可以达到实时性的要求。
以上实施例仅用以说明本发明的技术方案而非对其进行限制,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种农作物行的识别方法,其特征在于,所述识别方法包括:
通过模糊聚类方法将采集的农作物图像转化为二值图像;
从获取的二值图像中提取感兴趣的区域ROI图像,通过横条法获取其中农作物特征点;
对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程;
通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。
2.根据权利要求1所述的农作物行的识别方法,其特征在于,所述通过模糊聚类方法将采集的农作物图像转化为二值图像,是以采集的农作物图像作为整个样本,以像素点G通道的值所占的百分比作为样本元素,将农作物和背景的聚类中心分别初始化为0.35~0.40和0.30~0.35进行模糊聚类得到。
3.根据权利要求1所述的农作物行的识别方法,其特征在于,所述采集的农作物图像的宽为W像素,高为H像素,所述ROI图像的宽为w=W/2像素,宽为h=H/2,所述通过横条法获取其中农作物特征点,包括:
将ROI图像分为Q条宽度相同的横条,用Sp,q表示第q条横条第p列中白色像素出现的次数,其中p的取值从1到w,w为ROI图像的宽度像素;
对于第q条横条,对应有阈值uq,阈值uq为第q条横条中所有Sp,q的均值;
当Sp,q小于等于uq且Sp+1,q大于uq时,认为进入农作物行,记录此时的列坐标为p1;
当Sp,q大于等于uq且Sp+1,q小于uq时,认为离开农作物行,记录此时的列坐标为p2;
计算进入和离开农作物行时列数的差值Δ=p2-p1,如果Δ大于设定的常数d,则认为横条q上从第p1到p2的段为农作物,并取该段中点为农作物特征点;
遍历所有横条,获取ROI图像中所有的农作物特征点;
其中,所述常数d的取值范围为:W/20<d<W/15。
4.根据权利要求1所述的农作物行的识别方法,其特征在于,所述对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程,包括:
根据特征点的分布,将特征点划分为不同的农作物行;
对于任一农作物行,设农作物行直线方程为:
y=kx+b
其中b为斜距,k为斜率,计算属于该农作物行的所有特征点到该直线的距离l:
进一步计算,属于该农作物行的所有特征点到该直线的距离平方和l′:
对上式的分子对k和b求偏导并令其为0,得到:
其中M为属于该农作物行的所有特征点的数量,r属于1~M,第r个特征点的坐标为(xr,yr),lr为第r个特征点与农作物行直线的距离,求解出上式中k和b的解分别为和将其带入直线方程得到农作物直线方程为:
5.根据权利要求4所述的农作物行的识别方法,其特征在于,所述通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正,包括:
根据农作物行直线方程,计算属于该农作物行的特征点到该农作物行的距离,剔除距离大于设定的常数的特征点;
在剔除错误的特征点后,根据剩下的特征点再次拟合出农作物行直线方程,并再次根据新拟合出的农作物行直线方程,计算特征点到该直线的距离,剔除距离大于设定的常数的特征点;
如此循环,直到达到最大线性回归次数,或当线性回归剔除的特征点数量为0时迭代停止。
6.一种农作物行的识别装置,其特征在于,所述装置包括:
模糊聚类模块,用于通过模糊聚类方法将采集的农作物图像转化为二值图像;
特征点提取模块,用于从获取的二值图像中提取感兴趣的区域ROI图像,通过横条法获取其中农作物特征点;
拟合模块,用于对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程,并通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正。
7.根据权利要求6所述的农作物行的识别装置,其特征在于,所述模糊聚类模块在通过模糊聚类方法将采集的农作物图像转化为二值图像时,是以采集的农作物图像作为整个样本,以像素点G通道的值所占的百分比作为样本元素,将农作物和背景的聚类中心分别初始化为0.35~0.40和0.30~0.35进行模糊聚类得到。
8.根据权利要求6所述的农作物行的识别装置,其特征在于,所述采集的农作物图像的宽为W像素,高为H像素,所述ROI图像的宽为w=W/2像素,宽为h=H/2,所述特征点提取模块通过横条法获取其中农作物特征点时,执行如下操作:
将ROI图像分为Q条宽度相同的横条,用Sp,q表示第q条横条第p列中白色像素出现的次数,其中p的取值从1到w,w为ROI图像的宽度像素;
对于第q条横条,对应有阈值uq,阈值uq为第q条横条中所有Sp,q的均值;
当Sp,q小于等于uq且Sp+1,q大于uq时,认为进入农作物行,记录此时的列坐标为p1;
当Sp,q大于等于uq且Sp+1,q小于uq时,认为离开农作物行,记录此时的列坐标为p2;
计算进入和离开农作物行时列数的差值Δ=p2-p1,如果Δ大于设定的常数d,则认为横条q上从第p1到p2的段为农作物,并取该段中点为农作物特征点;
遍历所有横条,获取ROI图像中所有的农作物特征点;
其中,所述常数d的取值范围为:W/20<d<W/15。
9.根据权利要求6所述的农作物行的识别装置,其特征在于,所述拟合模块对所提取的农作物特征点进行线性回归,拟合出农作物行直线方程时,执行如下操作:
根据特征点的分布,将特征点划分为不同的农作物行;
对于任一农作物行,设农作物行直线方程为:
y=kx+b
其中b为斜距,k为斜率,计算属于该农作物行的所有特征点到该直线的距离l:
进一步计算,属于该农作物行的所有特征点到该直线的距离平方和l′:
对上式的分子对k和b求偏导并令其为0,得到:
其中M为属于该农作物行的所有特征点的数量,r属于1~M,第r个特征点的坐标为(xr,yr),lr为第r个特征点与农作物行直线的距离,求解出上式中k和b的解分别为和将其带入直线方程得到农作物直线方程为:
10.根据权利要求9所述的农作物行的识别装置,其特征在于,所述拟合模块在通过多次线性回归的方法剔除错误的特征点,对农作物行直线方程进行修正时,执行如下操作:
根据农作物行直线方程,计算属于该农作物行的特征点到该农作物行的距离,剔除距离大于设定的常数的特征点;
在剔除错误的特征点后,根据剩下的特征点再次拟合出农作物行直线方程,并再次根据新拟合出的农作物行直线方程,计算特征点到该直线的距离,剔除距离大于设定的常数的特征点;
如此循环,直到达到最大线性回归次数,或当线性回归剔除的特征点数量为0时迭代停止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510446384.4A CN105046229B (zh) | 2015-07-27 | 2015-07-27 | 一种农作物行的识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510446384.4A CN105046229B (zh) | 2015-07-27 | 2015-07-27 | 一种农作物行的识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105046229A true CN105046229A (zh) | 2015-11-11 |
CN105046229B CN105046229B (zh) | 2018-11-02 |
Family
ID=54452758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510446384.4A Expired - Fee Related CN105046229B (zh) | 2015-07-27 | 2015-07-27 | 一种农作物行的识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105046229B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107067430A (zh) * | 2017-04-13 | 2017-08-18 | 河南理工大学 | 一种基于特征点聚类的麦田作物行检测方法 |
KR20180009151A (ko) * | 2016-07-18 | 2018-01-26 | 강원대학교산학협력단 | 고온가압 살균된 통옥수수 제조 시의 살균도 조절을 위한 퍼지 클러스터링을 이용한 옥수수등급 판별기술 |
CN114002951A (zh) * | 2021-09-16 | 2022-02-01 | 江苏农林职业技术学院 | 一种水稻硬地育秧模糊控制灌溉方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521596A (zh) * | 2011-12-09 | 2012-06-27 | 中国科学院长春光学精密机械与物理研究所 | 有限认知条件下农作物目标识别与分割的方法 |
CN103632381A (zh) * | 2013-12-08 | 2014-03-12 | 中国科学院光电技术研究所 | 一种利用骨架提取特征点的扩展目标跟踪方法 |
-
2015
- 2015-07-27 CN CN201510446384.4A patent/CN105046229B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102521596A (zh) * | 2011-12-09 | 2012-06-27 | 中国科学院长春光学精密机械与物理研究所 | 有限认知条件下农作物目标识别与分割的方法 |
CN103632381A (zh) * | 2013-12-08 | 2014-03-12 | 中国科学院光电技术研究所 | 一种利用骨架提取特征点的扩展目标跟踪方法 |
Non-Patent Citations (1)
Title |
---|
汪博 等: ""基于最小二乘法的早期玉米作物行检测研究"", 《浙江理工大学学报(自然科学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20180009151A (ko) * | 2016-07-18 | 2018-01-26 | 강원대학교산학협력단 | 고온가압 살균된 통옥수수 제조 시의 살균도 조절을 위한 퍼지 클러스터링을 이용한 옥수수등급 판별기술 |
KR102524360B1 (ko) * | 2016-07-18 | 2023-04-24 | 강원대학교산학협력단 | 고온가압 살균된 통옥수수 제조 시의 살균도 조절을 위한 퍼지 클러스터링을 이용한 옥수수등급 판별기술 |
CN107067430A (zh) * | 2017-04-13 | 2017-08-18 | 河南理工大学 | 一种基于特征点聚类的麦田作物行检测方法 |
CN107067430B (zh) * | 2017-04-13 | 2020-04-21 | 河南理工大学 | 一种基于特征点聚类的麦田作物行检测方法 |
CN114002951A (zh) * | 2021-09-16 | 2022-02-01 | 江苏农林职业技术学院 | 一种水稻硬地育秧模糊控制灌溉方法 |
CN114002951B (zh) * | 2021-09-16 | 2023-12-29 | 江苏农林职业技术学院 | 一种水稻硬地育秧模糊控制灌溉方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105046229B (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108846832B (zh) | 一种基于多时相遥感影像与gis数据的变化检测方法及*** | |
CN111368736B (zh) | 基于sar和光学遥感数据的水稻精细化估产方法 | |
CN108710840B (zh) | 一种农田喷药机器人视觉导航路径识别方法 | |
CN111727457B (zh) | 一种基于计算机视觉的棉花作物行检测方法、装置及存储介质 | |
CN106951836A (zh) | 基于先验阈值优化卷积神经网络的作物覆盖度提取方法 | |
CN102110227B (zh) | 基于上下文关系的多分辨率遥感图像复合分类方法 | |
CN110909679B (zh) | 冬小麦历史种植区的休耕轮作信息遥感识别方法及*** | |
CN102867115B (zh) | 一种基于模糊c均值聚类的农田划分方法 | |
CN104598908A (zh) | 一种农作物叶部病害识别方法 | |
CN110175931B (zh) | 一种大范围快速提取农作物种植面积和物候信息的方法 | |
CN114818909B (zh) | 一种基于作物长势特征的杂草检测方法和装置 | |
Ma et al. | Automatic detection of crop root rows in paddy fields based on straight-line clustering algorithm and supervised learning method | |
CN102663397B (zh) | 一种小麦出苗的自动检测方法 | |
CN107680098A (zh) | 一种甘蔗蔗节特征的识别方法 | |
CN104008551A (zh) | 一种基于可见光图像的柑橘黄龙病检测方法 | |
Quinn et al. | Modeling and monitoring crop disease in developing countries | |
CN104751122B (zh) | 一种农作物病情指数的计算方法及*** | |
Zhou et al. | Autonomous detection of crop rows based on adaptive multi-ROI in maize fields | |
CN104952070A (zh) | 一种类矩形引导的玉米田遥感图像分割方法 | |
CN108364300A (zh) | 蔬菜叶部病害图像分割方法、***及计算机可读存储介质 | |
CN105046229A (zh) | 一种农作物行的识别方法及装置 | |
CN106909881A (zh) | 基于无人机遥感影像提取玉米育种基地垄数的方法及*** | |
CN105957115A (zh) | 广义dem思想下的农作物种植结构遥感提取方法 | |
Aitchison | Classification of agricultural systems | |
CN112734119A (zh) | 一种基于空间养分聚类分析的玉米产量预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181102 Termination date: 20200727 |