CN103714563B - 一种曲线形农田作业区域边界建模方法 - Google Patents
一种曲线形农田作业区域边界建模方法 Download PDFInfo
- Publication number
- CN103714563B CN103714563B CN201410002907.1A CN201410002907A CN103714563B CN 103714563 B CN103714563 B CN 103714563B CN 201410002907 A CN201410002907 A CN 201410002907A CN 103714563 B CN103714563 B CN 103714563B
- Authority
- CN
- China
- Prior art keywords
- point
- key
- key point
- abscissa
- last
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000004888 barrier function Effects 0.000 claims abstract description 7
- 238000005507 spraying Methods 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 238000013507 mapping Methods 0.000 description 4
- 239000000575 pesticide Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000007921 spray Substances 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种曲线形农田作业区域边界建模方法,采集曲线形农田作业区域边界的关键点,找出所有关键点中横坐标值最大的点和横坐标值最小的点,利用三次样条函数建立上边界模型和下边界模型。本发明可以将曲线形农田作业区域的边界以及农田作业区域内的障碍物顶点精准地记录并绘制出来,从而为农业植保机的作业人员提供准确的农田作业区域,为农业植保机的路径和航迹规划以及计算农田作业面积提供可靠的依据。
Description
技术领域
本发明涉及一种曲线形农田作业区域边界建模方法。
背景技术
随着科技的不断发展,采用现代化机械来替换人工劳动已成为各行各业的流行趋势。我国作为一个传统的农业大国,耕地面积十分广阔,然而,目前国内农田作业方面还是采用传统的人工作业方式。其中农药喷洒采取的是人工喷洒的方式,这种方式不仅效率低,而且对作业人员的身心有极大的伤害,因此,急需先进的技术来改变这一现象。农业植保机应运而生,然而,现今所有的植保机都是靠人工进行操作,对农田的喷洒区域及边界完全靠操作员的眼睛进行判断,难免产生误判、漏喷及多喷的现象。因此,如果操作员能够获得准确的农田作业区域的边界,在无人直升机进行作业时,实时的将路径和航迹在显示器上显示出来,操作员便可根据航迹判断出飞机的位置,从而很大程度的解决误判、漏喷及多喷的问题,且大大提高了农业喷洒的安全性。而这一问题的关键便在于要准确的获取农田的边界信息。因此,急需一种好的农田作业区域边界的建模方法来解决这一实际问题。
北京农业信息技术研究中心的发明专利“采集农田关键顶点测绘成图的方法”提供了一种采集农田关键顶点测绘成图的方法,通过采集四种交界顶点、两种轮廓顶点一级一种辅助顶点后,根据关键顶点间的拓扑关系实时自动校验并查看测绘结果。包括以下步骤:S1:获取GPS位置信息;S2:勾勒待测区域轮廓;S3:测绘农田关键分界顶点,并标注地块名称和注释;S4:实时校验、提示与分割农田;S5:上传测绘数据;S6:获取矢量地图。
其中,在S2中,勾勒待测区域轮廓包括:顺序标定和无序标定,所述顺序标定,以轮廓顶点的编号为顺序,依次标定待测区域轮廓;所述无序标定,是指自动标定包含全部轮廓顶点的面积最大的多边形为区域轮廓。
现有方法是以绘出农田的矢量地图为目的,而并非以画出农田作业区域边界为目的,因此,所绘出的农田轮廓包含两个部分:作业区域和非作业区域。在将此农田轮廓图用于农业植保机喷洒农药的路径和航迹规划等场合时,会出现非作业区域的多喷现象,浪费农药;在用于计算农田作业面积时,会出现作业面积的计算不准确现象;现有方法未对农田区域边界建立数学模型,从而无法为农田作业区域的路径和航迹规划提供可靠依据;没有将农田作业区域以内可能存在的障碍物顶点(如电线杆、树、信号发射塔等)标记出来,在将此农田轮廓图用于农业植保机喷洒农药的路径和航迹规划等场合时,可能导致农业植保机出现撞机或坠毁危险现象发生。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种曲线形农田作业区域边界建模方法,将曲线形农田作业区域的边界以及农田作业区域内的障碍物顶点精准地记录并绘制出来,从而为农业植保机的作业人员提供准确的农田作业区域,为农业植保机的路径和航迹规划以及计算农田作业面积提供可靠的依据。
为解决上述技术问题,本发明所采用的技术方案是:一种曲线形农田作业区域边界建模方法,该方法为:
1)从曲线形农田作业区域边界的某一个关键点D1开始,以顺时针方向采集曲线形农田作业区域边界的所有关键点,设共有m个关键点,分别为D1(x1,y1),D2(x2,y2),D3(x3,y3),…,Dm(xm,ym);建立直角坐标系,将所述m个关键点(每隔一个采样周期(GPS本身的采样频率的倒数)就自动采集一个点,当回到初始位置,即关键点D1所在位置时结束采点)设置在直角坐标系的第一象限内;
2)比较上述m个关键点的横坐标,找出横坐标值最大的关键点和最小的关键点;若有多个关键点的横坐标值均为最小值,且该多个关键点是第一个关键点D1与最后一个关键点Dm之间的关键点,则将次序最先的关键点作为横坐标值最小的关键点,若第一个关键点D1和最后一个关键点Dm的横坐标值都是最小值,则将最后一个关键点作为横坐标值最小的关键点;若有多个关键点的横坐标值均为最大值,且该多个关键点是第一个关键点D1与最后一个关键点Dm之间的关键点,则将次序最后的关键点作为横坐标值最大的关键点,若第一个关键点D1和最后一个关键点Dm的横坐标值都是最大值,则将第一个关键点D1作为横坐标值最大的关键点;将横坐标值最小的关键点Di和最大的关键点Dj分别作为所有的关键点在直角坐标系中的左端和右端;
3)以Di为起始点,Dj为终点;以顺时针方向,令Di=B0,Di+1=B1,…,Dj=Bn,将从Di到Dj的各点依次即为B0(x0,y0),K,Bn(xn,yn);建立以下上边界函数模型St(x):
其中,x∈[xt,xt+1,],t=0,1,...,n-1;xt、xt+1分别为B0、Bn之间的点Bt、Bt+1的横坐标,且Bt+1为Bt顺时针方向上的下一点;系数Mt、Mt+1通过追赶法求解;ht=xt+1-xt;yt、yt+1分别为点Bt、Bt+1的纵坐标;
4)以Di为起始点,Dj为终点;以逆时针方向,令Di=C0,Di-1=C2,…,Dj=Cr,建立以下上边界函数模型Sr(x):
其中,x∈[xr,xr+1,],r=0,1,...,m-n;xr、xr+1分别为C0、Cm-n之间的点Cr、Cr+1的横坐标,且Cr+1为Cr逆时时针方向上的下一点;系数Mr、Mr+1通过追赶法求解;hr=xr+1-xr;yr、yr+1分别为点Cr、Cr+1的纵坐标。
当所述曲线形农田作业区域内存在障碍物时,采集曲线形农田作业区域内障碍物顶点的坐标Z(x0,y0),建立以下障碍物警戒线函数方程:
一级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=a2;
二级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=b2;
三级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=c2;
其中,0<a<b<c<50。
通过追赶法求解系数Mt、Mt+1的步骤为:
1)令β1=b1,y1=d1;其中, λp=1-μp;p=1,2,...,n-1;
2)计算βq=bq-lqcq-1,yq=dq-lqyq-1;其中,q=2,3,…,n-1;
3)利用下式求解Mn-1:
4)利用下式计算Mn-2,Mn-3,L,M1:s=n-2,n-3,…,1;
5)通过下式计算M0和Mn:
与现有技术相比,本发明所具有的有益效果为:本发明可以将曲线形农田作业区域的边界以及农田作业区域内的障碍物顶点精准地记录并绘制出来,从而为农业植保机的作业人员提供准确的农田作业区域,为农业植保机的路径和航迹规划以及计算农田作业面积提供可靠的依据。
附图说明
图1为本发明一实施例曲线型边界模型示意图;
图2为本发明一实施例有障碍物的农田作业区域障碍物警戒线模型。
具体实施方式
以下结合附图详细说明本发明的具体实施方式。
1)假设采点工作者以第一个关键点开始,以顺时针方向依次采集m个关键点(以GPS的周期Ts作为采集农田作业区域边界点的采样周期),分别为D1(x1,y1),D2(x2,y2),D3(x3,y3),…,Dm(xm,ym),如图1所示;
2)比较这m个点的横坐标,找出其中横坐标最小和最大的点。假设横坐标最小的点为Di(xi,yi),横坐标最大的点为Dj(xj,yj)。(若遇到多个点的横坐标都为最小值,则把取点次序最先的点作为最小点,在此有个特例:如果D1和Dm的横坐标都相等,但是在此取Dm作为最小点;若遇到多个点的横坐标都为最大值,则把取点次序最后的点作为最大点,在此有个特例:如果Dm和D1的横坐标都相等,但是在此取D1作为最大点),那么Di(xi,yi)即为左端,Dj(xj,yj)即为右端;
3)假设由左端Di与右端Dj所确定的上边界有|j-i|+1个点,令n=|j-i|(包括左端Di与右端Dj),为便于接下来的计算,将上边界上各点(包括左端Di与右端Dj)从左至右按顺时方向依次记为B0(x0,y0),K,Bn(xn,yn)(即令Di=B0,Di+1=B1,…,Dj=Bn若当j<i,则在顺时针方向移点的过程中,如果Dn出现在Dj的左边(即xn<xj),令Dn的下一个点为D1,如果Dn=Bk(k<n)那么D1=Bk+1,D2=Bk+2,K,Dj=Bn。
则上边界函数S(x)满足:
(1)S(xt)=yt(j=0,1,…,n);
(2)S(x)在每个小区间[xt,xt+1](t=0,1,…,n-1)上是3次多项式;
(3)S(x)在[x0,xn]上有连续2阶导数;
(4)端点条件:强迫第一个和第二的三次多项式的三阶导数相等,对最后一个和倒数第二个三次多项式也做同样处理。
按求解3次样条函数的三弯矩方程方法,设S″(xt)=Mt,t=0,1,...,n,
记ht=xt+1-xt,t=0,1,...,n-1。
(1)整理三次样条函数的推导过程得μtMt-1+2Mt+λtMt+1=dt,t=1,2,...,n-1,其中, λt=1-μt,
(2)由端点条件知,
由μtMt-1+2Mt+λtMt+1=dt,t=1,2,...,n-1知μ1M0+2M1+λ1M2=d1
∴
由μtMt-1+2Mt+λtMt+1=dt,t=1,2,...,n-1知μn-1Mn-2+2Mn-1+λn-1Mn=dn-1,
∴
(3)结合(1)和(2)得矩阵
求出三弯矩阵的系数阵和右端项。
(4)另
则(3)中的矩阵变为
(5)用追赶法求出Mt,t=1,2,…,n-1。
具体过程如下:
①β1=b1,y1=d1;
②对于i=2,3,…,n-1,计算βi=bi-lici-1,yi=di-liyi-1;
③
④对于i=n-2,n-3,…,1,计算
⑤计算
(6)按下列公式求出S(x)在各个区间[xt,xt+1]上表达式St(x),其为
4)下边界函数模型
左端Di与右端Dj所确定的下边界上有n+1=m-|j-i|+1个点(包括左端Di与右端Dj),为便于后来计算,将下边界上各点(包括左端Di与右端Dj)从左至右按逆时针方向依次记为C0(x1,y1),K,Cn(xn,yn)(即令Di=C0,Di-1=C2,…,Dj=Cn,若当j>i,则在逆时针方向移点的过程中,如果D1出现在Dj的左边(即X1<Xj),那么令D1的下一个点为Dn;如果D1=Ck(k<n),那么Dn=Ck+1,Dn-1=Ck+2,...,Dj=Cn)。
下边界函数计算与上边界函数计算过程一样。
如果农田中有障碍物,则由采点工作者采集到障碍物点的坐标为Z(x0,y0),假设一级警戒线为离障碍物点a米的圆,二级警戒线为离障碍物点b米的圆,三级警戒线为离障碍物点c米的圆(50>c>b>a>0),如图2所示,图2中三个圆从内至外分别为一级警戒线、二级警戒线、三级警戒线。
则一级警戒线的函数方程为(x-x0)2+(y-y0)2=a2;
二级警戒线的函数方程为(x-x0)2+(y-y0)2=b2;
三级警戒线的函数方程为(x-x0)2+(y-y0)2=c2。
Claims (3)
1.一种曲线形农田作业区域边界建模方法,其特征在于,该方法为:
1)从曲线形农田作业区域边界的某一个关键点D1开始,以顺时针方向采集曲线形农田作业区域边界的所有关键点,设共有m个关键点,分别为D1(x1,y1),D2(x2,y2),D3(x3,y3),…,Dm(xm,ym);建立直角坐标系,将所述m个关键点设置在直角坐标系的第一象限内;
2)比较上述m个关键点的横坐标,找出横坐标值最大的关键点和最小的关键点;若有多个关键点的横坐标值均为最小值,且该多个关键点是第一个关键点D1与最后一个关键点Dm之间的关键点,则将次序最先的关键点作为横坐标值最小的关键点,若第一个关键点D1和最后一个关键点Dm的横坐标值都是最小值,则将最后一个关键点作为横坐标值最小的关键点;若有多个关键点的横坐标值均为最大值,且该多个关键点是第一个关键点D1与最后一个关键点Dm之间的关键点,则将次序最后的关键点作为横坐标值最大的关键点,若第一个关键点D1和最后一个关键点Dm的横坐标值都是最大值,则将第一个关键点D1作为横坐标值最大的关键点;将横坐标值最小的关键点Di和最大的关键点Dj分别作为所有的关键点在直角坐标系中的左端和右端;若有多个关键点的横坐标值均为最小值,其中一个关键点为第一个关键点或最后一个关键点,其余关键点为第一个关键点与最后一个关键点之间的关键点,则取从第一个关键点到最后一个关键点的闭区间内次序最先的点作为最小点;若有多个关键点的横坐标值均为最大值,其中一个关键点为第一个关键点或最后一个关键点,其余关键点为第一个关键点与最后一个关键点之间的关键点,则取从第一个关键点到最后一个关键点的闭区间内次序最后的点作为最大点;
3)以Di为起始点,Dj为终点;以顺时针方向,令Di=B0,Di+1=B1,…,Dj=Bn,将从Di到Dj的各点依次即为B0(x0,y0),...,Bn(xn,yn);建立以下上边界函数模型St(x):
其中,x∈[xt,xt+1,],t=0,1,...,n-1;xt、xt+1分别为B0、Bn之间的点Bt、Bt+1的横坐标,且Bt+1为Bt顺时针方向上的下一点;系数Mt、Mt+1通过追赶法求解;
ht=xt+1-xt,t=0,1,...,n-1;yt、yt+1分别为点Bt、Bt+1的纵坐标;
4)以Di为起始点,Dj为终点;以逆时针方向,令Di=C0,Di-1=C1,…,Dj=Cr,建立以下上边界函数模型Sr(x):
其中,x∈[xr,xr+1,],r=0,1,...,m-n;xr、xr+1分别为C0、Cm-n之间的点Cr、Cr+1的横坐标,且Cr+1为Cr逆时时针方向上的下一点;系数Mr、Mr+1通过追赶法求解;hr=xr+1-xr;yr、yr+1分别为点Cr、Cr+1的纵坐标。
2.根据权利要求1所述的曲线形农田作业区域边界建模方法,其特征在于,当所述曲线形农田作业区域内存在障碍物时,采集曲线形农田作业区域内障碍物顶点的坐标Z(x0,y0),建立以下障碍物警戒线函数方程:
一级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=a2;
二级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=b2;
三级障碍物警戒线的函数方程为(x-x0)2+(y-y0)2=c2;
其中,0<a<b<c<50。
3.根据权利要求1或2所述的曲线形农田作业区域边界建模方法,其特征在于,通过追赶法求解系数Mt、Mt+1的步骤为:
1)令β1=b1,y1=d1;其中,b2=2,...,bn-2=2, λp=1-μp;p=1,2,...,n-1;
2)计算βq=bq-lqcq-1,yq=dq-lqyq-1;其中,q=2,3,…,n-1;a2=μ2,c2=λ2,c3=λ3,...,cn-2=λn-2;
3)利用下式求解Mn-1:
4)利用下式计算Mn-2,Mn-3,…,M1:s=n-2,n-3,...,1;
5)通过下式计算M0和Mn:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410002907.1A CN103714563B (zh) | 2014-01-03 | 2014-01-03 | 一种曲线形农田作业区域边界建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410002907.1A CN103714563B (zh) | 2014-01-03 | 2014-01-03 | 一种曲线形农田作业区域边界建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103714563A CN103714563A (zh) | 2014-04-09 |
CN103714563B true CN103714563B (zh) | 2016-11-02 |
Family
ID=50407504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410002907.1A Expired - Fee Related CN103714563B (zh) | 2014-01-03 | 2014-01-03 | 一种曲线形农田作业区域边界建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103714563B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155993B (zh) * | 2014-08-11 | 2017-07-25 | 江苏恒创软件有限公司 | 一种基于无人直升机的坡地区域农作物作业方法 |
CN104503464B (zh) * | 2014-12-30 | 2017-01-18 | 中南大学 | 基于计算机的凸多边形农田无人机喷洒作业航迹规划方法 |
JP6570598B2 (ja) * | 2017-11-08 | 2019-09-04 | 本田技研工業株式会社 | 地図生成サーバおよび地図生成方法 |
CN110647141B (zh) * | 2018-06-27 | 2022-11-08 | 西安合众思壮导航技术有限公司 | 避障路径的生成方法、装置和*** |
CN110378951A (zh) * | 2018-09-10 | 2019-10-25 | 天津京东深拓机器人科技有限公司 | 确定堆垛物品支撑面积的方法和装置 |
CN111626148B (zh) * | 2020-05-09 | 2023-07-11 | 浙江数治空间规划设计有限公司 | 一种无人机农田核查方法、***、智能终端及存储介质 |
CN112233131B (zh) * | 2020-10-22 | 2022-11-08 | 广州极飞科技股份有限公司 | 一种地块分割方法、装置、设备及存储介质 |
CN115855067B (zh) * | 2023-02-23 | 2023-06-27 | 华南农业大学 | 一种曲形农田边界的路径规划方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1534274A (zh) * | 2003-03-28 | 2004-10-06 | 中国科学院力学研究所 | 一种空间三维自由曲线的测量方法 |
CN101226630A (zh) * | 2007-09-27 | 2008-07-23 | 四川虹微技术有限公司 | 基于三次样条函数的插值方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9013480B2 (en) * | 2011-10-21 | 2015-04-21 | IntegrityWare, Inc. | Methods and systems for generating and editing surfaces |
-
2014
- 2014-01-03 CN CN201410002907.1A patent/CN103714563B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1534274A (zh) * | 2003-03-28 | 2004-10-06 | 中国科学院力学研究所 | 一种空间三维自由曲线的测量方法 |
CN101226630A (zh) * | 2007-09-27 | 2008-07-23 | 四川虹微技术有限公司 | 基于三次样条函数的插值方法 |
Non-Patent Citations (1)
Title |
---|
三次样条插值函数的快速生成;李珉 等;《科技视界》;20120727(第13期);第77-79页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103714563A (zh) | 2014-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103714563B (zh) | 一种曲线形农田作业区域边界建模方法 | |
CN103679774B (zh) | 一种多边形农田作业区域边界建模方法 | |
EP3384243B1 (en) | Path planning for area coverage | |
JP5821970B2 (ja) | 算出方法、算出プログラムおよび算出装置 | |
CN103673937A (zh) | 一种双重算法模式农机田作业面积的计算***及其方法 | |
WO2019179270A1 (zh) | 植株种植数据测量方法、作业路线规划方法及装置、*** | |
CN111868566B (zh) | 一种基于定位漂移测算模型的农机作业面积测算方法 | |
CN109240284A (zh) | 一种无人驾驶农机的自主路径规划方法及装置 | |
EP3646258A1 (en) | Computer platform for controlling agricultural assets | |
CN108549409A (zh) | 一种植保无人机飞行控制方法 | |
CN107462208A (zh) | 一种农机及农机作业面积测量装置和测量方法 | |
CN105355042B (zh) | 一种基于出租车gps的道路网络提取方法 | |
CN109597420A (zh) | 一种智能割草机边界自动闭合处理***及方法 | |
CN110619649A (zh) | 作业面积确定方法及装置、终端 | |
CN105573318A (zh) | 基于概率分析的环境构建方法 | |
CN103136393A (zh) | 一种基于网格划分的区域覆盖率计算方法 | |
CN103679775B (zh) | 一种直线和曲线相结合的农田作业区域边界建模方法 | |
CN108445912A (zh) | 一种植保无人机飞行控制*** | |
CN110715657A (zh) | 一种航空施药区域全覆盖路径规划方法 | |
CN107402375A (zh) | 一种带观测时滞的室内行人定位efir数据融合***及方法 | |
CN115855067B (zh) | 一种曲形农田边界的路径规划方法 | |
CN103837196A (zh) | 一种尾矿库干滩参数监测***及方法 | |
CN106568412A (zh) | 一种不规则形状地块的面积测量方法 | |
CN103548461A (zh) | 一种基于虚拟gps的精准农业自动/手动变量施肥机 | |
CN205563277U (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161102 |
|
CF01 | Termination of patent right due to non-payment of annual fee |