CN101285686A - 一种农业机械导航分级定位的方法和系统 - Google Patents
一种农业机械导航分级定位的方法和系统 Download PDFInfo
- Publication number
- CN101285686A CN101285686A CNA2008101135922A CN200810113592A CN101285686A CN 101285686 A CN101285686 A CN 101285686A CN A2008101135922 A CNA2008101135922 A CN A2008101135922A CN 200810113592 A CN200810113592 A CN 200810113592A CN 101285686 A CN101285686 A CN 101285686A
- Authority
- CN
- China
- Prior art keywords
- positional parameter
- agricultural machinery
- parameter
- value
- constantly
- 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
Images
Landscapes
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种农业机械导航分级定位的方法,包括获取第一定位参数;利用卡尔曼滤波器对第一定位参数进行融合,获取第二定位参数;获取第三定位参数;利用多传感器自适应加权融合算法将第二定位参数与第三定位参数进行融合,获取目标定位参数;根据目标定位参数对农业机械进行定位。本发明还公开了一种农业机械导航分级定位的系统。本发明通过对各测量设备获取的定位参数行多次融合,可平滑DGPS定位数据,有效避免DGPS动态定位的异常结果、有效滤除测试噪声、减小系统误差,能够形成连续、稳定、精确的农业机械的位置和航向角度信息,提高了农业机械导航定位的准确性。
Description
技术领域
本发明涉及导航技术领域,特别是涉及一种农业机械导航分级定位的方法和系统。
背景技术
导航定位是农业机械导航中的首要问题,导航定位的精度直接影响到农业机械跟踪预定义路径的质量。农业机械导航定位的主要信息包括:位置信息和航向角度信息、定位方法可分为绝对定位方法和相对定位方法。DGPS(Differential Global Positioning System,差分全球定位系统)作为典型的绝对定位系统,能够全天候为接收机提供绝对位置、速度和方向等信息;但该系统需要接收四颗以上卫星才能决定车辆的位置,当受到外界干扰时,定位误差将大大增加。机器视觉属于相对定位方法,具有探测范围宽、信号丰富完整的特点,在提供导航信息的同时还可获得田间作物分布信息,其适应能力强,比较灵活,不需预先设定导航路线,非常适合行间作业的农业机械导航;但当地物标志信息不完整、外界光照条件变化时,都可能导致机器视觉导航定位系统的功能失效,并且它只能提供微观尺度的导航定位数据,不适合进行大田内宏观的路径规划及行走方向的引导。电子罗盘、加速度计、陀螺仪和里程计等都是广泛使用的传感器,短时间内具有很高的可信度;但由于随机误差和系统误差的存在,测量信号将产生飘移,一段时间后将严重影响导航性能。
可见,单一导航定位方法都存在自身难以克服的缺点,因此在国内外的研究中,多采用组合导航定位方法实现农业机械的高精度定位。目前,常见的组合导航定位方法主要有基于光纤陀螺仪或电子罗盘的DR(Dead Reckoning,航位推算)技术、基于GIS(Geographicinformation system,地理信息系统)的地图匹配技术、采用RTK-DGPS(Real Time Kinematic DGPS,实时动态DGPS)和FOG(Fiber-OpticsGyroscopes,光纤陀螺)实现精确定位的技术以及采用低价格低精度GPS(Global Positioning System,全球定位系统)融合高精度惯性传感器,通过卡尔曼滤波方法实现农业机械精确定位的技术等等。
在实现本发明过程中,发明人发现现有技术中至少存在如下问题:采用现有的组合导航定位技术时未充分考虑卡尔曼滤波自适应、难以避免滤波发散,从而造成农业机械导航定位的准确性较低,而且成本较高。
发明内容
本发明实施例要解决的问题是提供一种农业机械导航分级定位的方法和系统,以克服现有技术中农业机械导航定位的准确性较低的缺陷。
为达到上述目的,本发明实施例的技术方案提供一种农业机械导航分级定位的方法,所述方法包括以下步骤:获取第一定位参数,所述第一定位参数包括通过差分全球定位系统DGPS接收机获取的农业机械的第一位置值,通过电子罗盘获取的所述农业机械的第一航向角度值,通过对从加速度计得到的加速度值进行积分而获取的所述农业机械的当前行进速度值;利用卡尔曼滤波器对所述第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数,所述第二定位参数包括所述农业机械的第二位置值和第二航向角度值;获取第三定位参数,所述第三定位参数包括通过摄像装置获取的所述农业机械的第三位置值和第三航向角度值;利用多传感器自适应加权融合算法将所述第二定位参数与所述第三定位参数进行融合,获取目标定位参数,所述目标定位参数包括所述农业机械的目标位置值和目标航向角度值;根据所述目标定位参数对所述农业机械进行定位。
其中,在所述利用卡尔曼滤波器对第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数之前,还包括构建卡尔曼滤波器,所述构建卡尔曼滤波器的步骤具体包括:建立基于常速度的卡尔曼滤波器的状态方程和观测方程;根据所述状态方程和观测方程,建立预测方程组和校正方程组;根据所述预测方程组和校正方程组构建卡尔曼滤波器。
其中,所述状态方程为:
所述观测方程为:
Z(t)=H(t)X(t)+W(t)
其中,Z(t)为外部观测向量,H(t)为外部观测向量和状态向量之间的测量矩阵,X(t)为t时刻农业机械的状态向量,W(t)为高斯白噪声序列。
其中,所述预测方程组为:
其中,为t时刻农业机械的运动状态,Ф(t-1)为t-1时刻状态转移矩阵,为t-1时刻农业机械的运动状态,P(t|t)为t时刻滤波误差方差阵,φ(t|t)为t时刻状态转移矩阵,P(t-1|t-1)为t-1时刻滤波误差方差阵,φT(t|t)为t时刻状态转移矩阵的转置矩阵,G(t-1)为t-1时刻过程噪声分布矩阵,Q(t-1)为t-1时刻正定协方差矩阵,GT(t-1)为t-1时刻过程噪声分布矩阵的转置矩阵;
所述校正方程组为:
其中,为t时刻农业机械的运动状态,为t-1时刻农业机械的运动状态,P(t|t)为t时刻滤波误差方差阵,H(t)为t时刻外部观测向量和状态向量之间的测量矩阵,HT(t)为t时刻外部观测向量和状态向量之间的测量矩阵的转置矩阵,R-1(t)为t时刻正定协方差矩阵的逆矩阵,Z(t)为t时刻外部观测向量,P-1(t|t)为t时刻滤波误差方差阵的逆矩阵,P-1(t|t-1)为t-1时刻滤波误差方差阵。
其中,在所述卡尔曼滤波器中,利用DGPS接收机得到的第一位置值对所述状态转移矩阵进行实时更新。
其中,所述利用卡尔曼滤波器对第一位置值、第一航向角度值和当前行进速度值进行融合,具体包括:利用所述第一位置值、第一航向角度值和当前行进速度值更新所述预测方程组的矩阵;通过所述预测方程组获取下一时刻的状态估计;利用所述校正方程组对所述下一时刻的状态估计进行校正,获取当前最优估计;根据所述当前最优估计,获取当前最优估计的误差协方差矩阵。
其中,在所述利用多传感器自适应加权融合算法将第二定位参数与第三定位参数进行融合,获取目标定位参数之前,还包括以下步骤:判断所述DGPS接收机和摄像装置是否均能正常工作,如果所述DGPS接收机和摄像装置都不能正常工作,则停止导航;如果所述DGPS接收机和摄像装置中只有一个能正常工作,则判断能正常工作的装置得到的定位参数是否在误差范围内,如果是则采用单一传感器进行导航,否则停止导航;如果所述DGPS接收机和摄像装置都能正常工作,则判断能正常工作的装置得到的定位参数是否在误差范围内,如果所述DGPS接收机和摄像装置得到的定位参数都不在误差范围内,则停止导航;如果所述DGPS接收机和摄像装置得到的定位参数中只有一个在误差范围内,则采用单一传感器进行导航;如果所述DGPS接收机和摄像装置得到的定位参数都在误差范围内,则利用多传感器自适应加权融合算法将所述DGPS接收机得到的第二定位参数和摄像装置得到的第三定位参数进行融合。
其中,所述利用多传感器自适应加权融合算法将第二定位参数与第三定位参数进行融合,具体包括:
根据各测量设备获取的定位参数的方差,获取各测量设备的最优加权因子,所述最优加权因子由公式
获取;其中,Wp为最优加权因子,σ2为各测量设备获取的定位参数的方差;
根据公式
获取融合后的定位参数;其中,为融合后的定位参数,δq 2为摄像装置获取的定位参数的方差,Xp为DGPS接收机得到的定位参数值,δp 2为DGPS接收机获取的定位参数的方差,Xq为摄像装置得到的定位参数值。
本发明实施例的技术方案还提供了一种农业机械导航分级定位的系统,所述系统包括:DGPS接收机,用于获取农业机械的第一位置值;电子罗盘,用于获取所述农业机械的第一航向角度值;加速度计,用于获取所述农业机械的加速度值;摄像装置,用于获取所述农业机械的第三位置值和第三航向角度值;目标定位参数获取装置,用于对所述第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数,并将所述第二定位参数与所述第三定位参数进行融合,获取目标定位参数;定位装置,用于根据所述目标定位参数对所述农业机械进行定位。
其中,所述目标定位参数获取装置包括:卡尔曼滤波器,用于对所述第一位置值、第一航向角度值和当前行进速度值进行融合;多传感器自适应加权融合单元,用于将所述第二定位参数与所述第三定位参数进行融合。
上述技术方案仅是本发明的一个优选技术方案,具有如下优点:本发明实施例通过对各测量设备获取的定位参数行多次融合,可平滑DGPS定位数据,有效避免DGPS动态定位的异常结果、有效滤除测试噪声、减小系统误差,能够形成连续、稳定、相对精确的农业机械的位置和航向角度信息,提高了农业机械导航定位的准确性。
附图说明
图1是本发明实施例的一种农业机械导航分级定位的方法的流程图;
图2是本发明实施例的卡尔曼滤波器的工作原理图;
图3是本发明实施例的多传感器自适应加权融合算法的原理示意图;
图4是本发明实施例的多传感器自适应加权融合算法的流程图;
图5是本发明实施例的一种农业机械导航分级定位的系统的结构图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
发明实施例的一种农业机械导航分级定位的方法的流程如图1所示,所述方法包括以下步骤:
步骤s101,获取第一定位参数,所述第一定位参数包括通过差分全球定位系统DGPS接收机获取的农业机械的第一位置值,通过电子罗盘获取的所述农业机械的第一航向角度值,通过对从加速度计得到的加速度值进行积分而获取的所述农业机械的当前行进速度值。
步骤s102,构建卡尔曼滤波器。首先建立基于常速度的卡尔曼滤波器的状态方程和观测方程,然后根据所述状态方程和观测方程,建立预测方程组和校正方程组,最后根据所述预测方程组和校正方程组构建卡尔曼滤波器。
步骤s103,利用卡尔曼滤波器对所述第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数,所述第二定位参数包括所述农业机械的第二位置值和第二航向角度值。其中,利用卡尔曼滤波器对第一位置值、第一航向角度值和当前行进速度值进行融合的步骤包括:利用所述第一位置值、第一航向角度值和当前行进速度值更新所述预测方程组的矩阵;通过所述预测方程组获取下一时刻的状态估计;利用所述校正方程组对所述下一时刻的状态估计进行校正,获取当前最优估计;根据所述当前最优估计,获取当前最优估计的误差协方差矩阵。
步骤s104,获取第三定位参数,所述第三定位参数包括通过摄像装置获取的所述农业机械的第三位置值和第三航向角度值。从每一时刻采集的图像中提取导航线,得到导航线的特征参数,即导航线两端点坐标值。导航线的特征参数建立在图像坐标系下,通过坐标转换关系,可以得到导航线特征参数在世界坐标系下的对应位置值。已知导航线方程的两点坐标,根据公式θ=atan[(x2-x1)/(y2-y1)]可求得导航线与Y轴夹角,即农机与导航线之间的航向角度值。
步骤s105,利用多传感器自适应加权融合算法将所述第二定位参数与所述第三定位参数进行融合,获取目标定位参数,所述目标定位参数包括所述农业机械的目标位置值和目标航向角度值。在获取目标定位参数之前,还包括:判断所述DGPS接收机和摄像装置是否均能正常工作,如果所述DGPS接收机和摄像装置都不能正常工作,则停止导航;如果所述DGPS接收机和摄像装置中只有一个能正常工作,则判断能正常工作的装置得到的定位参数是否在误差范围内,如果是则采用单一传感器进行导航,否则停止导航;如果所述DGPS接收机和摄像装置都能正常工作,则判断能正常工作的装置得到的定位参数是否在误差范围内,如果所述DGPS接收机和摄像装置得到的定位参数都不在误差范围内,则停止导航;如果所述DGPS接收机和摄像装置得到的定位参数中只有一个在误差范围内,则采用单一传感器进行导航;如果所述DGPS接收机和摄像装置得到的定位参数都在误差范围内,则利用多传感器自适应加权融合算法将所述DGPS接收机得到的第二定位参数和摄像装置得到的第三定位参数进行融合。利用多传感器自适应加权融合算法将第二定位参数与第三定位参数进行融合的步骤包括:根据各测量设备获取的定位参数的方差,获取各测量设备的最优加权因子,所述最优加权因子由公式
获取;其中,Wp为最优加权因子,σ2为各测量设备获取的定位参数的方差;根据公式
获取融合后的定位参数;其中,为融合后的定位参数,δq 2为摄像装置获取的定位参数的方差,Xp为DGPS接收机得到的定位参数值,δp 2为DGPS接收机获取的定位参数的方差,Xq为摄像装置得到的定位参数值。
步骤s106,根据所述目标定位参数对所述农业机械进行定位。
本实施例所采用的导航传感器主要包括DGPS接收机、摄像装置采用CCD(Charge Coupled Device,电荷耦合器件)摄像头、电子罗盘和加速度计。利用DGPS接收机采集农业机械在WGS-84(WorldGeodetic System-84,1984年世界大地坐标系)下各位置点的坐标,作为初步定位信息;采用电子罗盘测量农业机械的航向角度信息;采用加速度计测其当前时刻的加速度,再通过积分得到当前时刻的速度值。将上述各传感器组合安装在农业机械上,利用数据采集软件系统动态采集数据,设定采集频率为1Hz。数据的采集、预处理和信息融合等均由软件系统实现。
本发明的融合多传感器信息的农业机械导航分级定位方法通过两级融合,逐步对初步DGPS定位数据、航向角度数据、速度数据以及机器视觉获得的定位数据进行融合,其具体实施步骤如下:
1.采用DGPS接收机测量得到农业机械的初步定位信息:
DGPS系统中,采用的是WGS-84坐标系,属于地心坐标系,为使DGPS定位数据能用于导航控制系统,需要进行高斯投影变换。
地球椭球上任一点的位置,既可以表示为,又可以表示为。两个坐标系之间的转换公式为:
式中,N为椭球的卯酉圈曲率半径,a为地球椭球长半径,b为短半径,e为椭球的第一偏心率。
N=a/W (2)
W=(1-e2sin2B)1/2 (3)
e2=(a2-b2)/a2 (4)
2.采用卡尔曼滤波器对DGPS得到的农业机械的位置值和电子罗盘得到的航向角度值及通过对加速度计得到的初始值积分得到的当前速度值进行融合:
(1)一般来说,农业机械的作业环境较为平整,工作时要求农业机械的行进速度较为缓慢并且恒定,除地头转弯不会有较大的转向动作。因此,假设农业机械的行进速度恒定且转向角加速度为0,则利用Kalman滤波器便可得到如下的状态空间描述:
X(t+1)=Ф(t)X(t)+G(t)V(t) (5)
其中,X(t)为t时刻农业机械的状态向量,且
X(t)=(x(t),y(t),v(t))T (6)
其中,x(t),y(t)为t时刻WGS-84坐标系下农业机械在高斯投影平面上的位置点坐标值,v(t)为t时刻农业机械的行进速度,Ф(t)为状态转移矩阵,G(t)为过程噪声分布矩阵,V(t)为零均值高斯过程白噪声向量,即
E[V(t)]=0,E[V(u)VT(j)]=Q(t)δkj (7)
设φ(t)为t时刻农业机械的航向角在高斯投影平面坐标下的角度分量,横轴正向为0,逆时针为正,则状态转移矩阵定义为
可利用上一时刻农业机械的航向角度和前轮转角的实际值对φ(t)矩阵进行实时更新。
(2)Kalman滤波器的基本观测方程:
Z(t)=H(t)X(t)+W(t) (8)
利用基本观测方程,将DGPS接收机得到的农业机械当前位置点坐标(x,y),速度值v作为外部观测量,则外部观测向量Z(t)=[x(t)y(t)v(t)],外部观测向量和状态向量之间的测量矩阵为H(t),是常数矩阵:
W(t)是均值为0、方差为ri 2的高斯白噪声序列。测量的噪声方差矩阵R(t)为:
r1、r2、r3分别表示DGPS接收机、电子罗盘、加速度计的测量噪声的标准差。
(3)Kalman滤波器的工作流程:
在Kalman滤波器的状态方程和观测方程的基础上,其递归差分方程组可用如下方程表示:
P-1(t|t)=P-1(t|t-1)+HT(t)R-1(t)H(t) (10)
式(9),(10)为Kalman滤波器的校正方程组,由它们可以得到
为便于算法的实现,对式(11),(12)做进一步推导,可以得到
其中,
P(t|t)=φ(t|t)P(t-1|t-1)φT(t|t)+G(t-1)Q(t-1)GT(t-1)(16)
式(15)、(16)为Kalman滤波器的预测方程组。
在上述各方程中,方程(9)为滤波方程,利用测量值Z(t)对预测的系统状态估计进行更新。方程(10)为滤波误差方差方程,利用计算出的Kalman滤波器增益对预测的滤波误差方差阵进行更新。方程(15)为Kalman滤波器的状态方程,代入t-1时刻的状态预测值,可对t时刻农业机械的运动状态进行预测。方程(16)为预报误差方差方程,可以对滤波误差方差阵进行预测。方程(14)说明农业机械的状态估计信息融合等价于状态估计与校正的融合结果之和。预测校正包括两部分,其一是各个传感器跟踪与预测差值(跟踪残留误差)的总和,每个差值由逆(预测)协方差加权;其二是整个系统预测与各个传感器预测的差值(预测残留误差)的总和,由局部逆(预测)协方差加权。
通常Kalman滤波器的应用中只用一种残留误差——跟踪残留误差来校正状态估计的融合,但本发明采用两种。跟踪残留误差反映了整个系统中的不可预测信息所带来的误差,用来校正整个系统的状态估计。但由于各传感器的状态估计与整个系统的状态估计相互关联,所以跟踪残留误差不能完全反映不可预测信息所带来的误差。因而式(14)采用预测残留误差来校正,显然比式(12)具有更明确的物理意义,且更容易实现。
(4)Kalman滤波器的工作原理:
首先将初始状态估计和初始滤波误差方差阵P(t-1)估计引入Kalman滤波器的预测方程组进行Kalman滤波。采用DGPS、电子罗盘及加速度计采集到的数据更新预测方程组的矩阵,再由预测方程组对下一时刻状态进行估计。预测过程结束后再利用Kalman滤波器进行校正,将R矩阵引入校正方程组,根据当前测量结果(针对预测方程组对下一时刻状态的估计结果),对预测方程组的状态估计进行更新,获得当前最优估计,最后解算当前最优估计的误差协方差矩阵,完成一次滤波。Kalman滤波器的工作原理图如图2所示。
3.采用多传感器自适应加权融合算法将机器视觉提供的位置值和航向角度值与初次融合得到相对精确的农业机械的位置和航向角度信息进行再次融合:
不同传感器有各自不同的加权因子,多传感器自适应加权融合算法的中心思想是在总均方误差最小这一最优条件下,根据各个传感器所得到的测量值以自适应的方式寻找各个传感器所对应的最优加权因子,使融合后的值达到最优。
多传感器自适应加权融合算法的原理示意图如图3所示,具体实现步骤如下:
(1)最优加权因子的求取
设n个传感器的方差分别为σ1 2,σ2 2,…,σn 2,所需估计的真值为X,各传感器的测量值分别为X1,X2,…,Xn,各测量值相互独立,并且是X的无偏估计;各传感器的加权因子分别为W1,W2,…,Wn,融合后的和加权因子满足:
总均方误差为:
因为X1,X2,…Xn相互独立,并且为X的无偏估计,所以:
故σ2可以写成:
总均方误差σ2是关于各加权因子的多元二次函数,因此,σ2必然存在最小值。该最小值的求取是加权因子满足公式(18)约束条件的多元函数极值的求取。根据多元函数求极值理论,可求出总均方误差最小时对应的加权因子为:
此时对应的最小均方误差为:
(2)测量方差的加窗估计算法
设DGPS接收机和机器视觉传感器的测量值分别为Xp和Xq,对应的观测误差分别为Vp和Vq,其中Vp和Vq互不相关,且均值为0。Xp和Xq的自协方差函数分别为Rpp和Rqq,互协方差函数为Rpq,第k次采样时Rpp的时间域估计值为Rpp(k),Rqq的时间域估计值为Rqq(k),Rpq的时间域估计值为Rpq(k),则由前N次采样的Rpp,Rqq和Rpq的时间域估计值可以得到如下递推公式:
随着系统长时间运行,系统参数会发生缓慢变化,噪声特征也会发生某种程度的缓变,为了在保证估计精度的前提下及时跟踪噪声的变化提出加窗方差估计方法:
当k<N时,
当k≥N时,
式中μ为采样数据的平均值,采用上述递推算式时,为实时获取μ,将得到的加权融合后的导航参数做为测量数据的μ。N为移动窗口长度。则DGPS接收机和机器视觉传感器的方差δp 2和δq 2分别为:
δp 2=Rpp-Rpq (33)
δq 2=Rqq-Rpq (34)
此时对应的最小均方误差为:
(3)多传感器自适应加权融合算法的实现:
首先判断DGPS接收机和机器视觉传感器是否均能正常工作:
1)如都不能正常工作则停止导航;
2)如只有一个正常工作,则判断其得到的定位参数是否在误差范围内,如果是则选择单一传感器导航方法,如果否则停止导航;
3)如果二者都能正常工作,仍需判断其定位数据是否在误差范围内:如两个定位数据均不在误差范围内则停止导航;如只有一个定位数据值在误差范围内则采用单一传感器导航方式;如两定位数据均在误差范围内,利用自适应加权融合估计的方法进行信息融合以提高定位的精度。
判断机器视觉传感器正常工作的标准是软件能清楚地采集图像,并对图像进行处理;判断DGPS接收机正常工作的标准是其能正常接收卫星信号,并保证能进行差分。利用设定阈值的方法判断定位数据是否满足误差范围,阈值的选择和农业机械的前轮间距有关。本实施例中农业机械的前轮间距为0.78m,因此设定误差范围为[-0.5m,0.5m]。
多传感器信息融合自适应加权融合算法不要求知道传感器测量数据的任何先验知识,只需传感器所提供的测量数据即可以自适应的方式,寻找各自所对应的最优加权因子,计算出融合后的数据,并可保证融合后总均方误差最小,使融合后的估计值达到最优,其算法流程图如图4所示。
4.三种定位方法的定位效果分析:
为验证本发明的融合多传感器信息的农业机械导航分级定位方法的稳定性,通过实验对定位效果进行测试,实验过程如下:
(1)选取一条东西方向的直线作为预定义路径,首先采用GPS4700接收机在预定义直线轨迹的两端点A,B处各进行30分钟静态定位,并对测量数据进行平均,取得两端点的精确定位结果,A点为(444729.280301,4429977.863371),B点为(444759.627519,4429977.651654),再利用解析几何方法获得该预定义路径在高斯投影平面坐标系中的几何描述;
(2)由于GPS4700的定位频率为1Hz,机器视觉传感器的定位频率约为10Hz,因此设定两传感器信息的融合频率为1Hz。实时纪录融合后的定位信息,将三种定位方法得到的定位信息进行比较,并通过统计学方法对数据进行分析;
(3)农业机械和预定义路径之间的位置偏差XTE是定位和导航精度的主要评价参数,获取准确的XTE数据是评价系统精度的基础。在获取XTE参数的基础上,对XTE进行统计分析。选取XTE的平均值|XXTE|,最大值|XXTE|max和标准差σXTE三项数据作为定位精度的评价指标。
实验结果分析:
首先对实验数据进行定性分析,然后通过统计学方法对实验数据进行定量分析。分析结果如表1所示。
表1
机器视觉(m) | 4700(m) | 融合(m) | |
|XXTE| | 0.076999 | 0.044632 | 0.004417 |
|XXTE|max | 3.51949 | 0.172502 | 0.112359 |
σXTE | 0.080704 | 0.034517 | 0.029885 |
从表1可看出融合后的统计量平均值,最大值和标准差均为最小。而且,融合算法去除了机器视觉定位中的跳变点,消除了可能出现误差的数据,提高了导航定位的稳定性和精度。同时,融合后的定位数据的波动最小,稳定性最好。
本发明实施例的一种农业机械导航分级定位的系统如图5所示,该系统包括DGPS接收机51、电子罗盘52、加速度计53、摄像装置54、目标定位参数获取装置55和定位装置56,其中目标定位参数获取装置55分别与DGPS接收机51、电子罗盘52、加速度计53、摄像装置54、定位装置56连接。
DGPS接收机51用于获取农业机械的第一位置值;电子罗盘52用于获取所述农业机械的第一航向角度值;加速度计53用于获取所述农业机械的加速度值;摄像装置54用于获取所述农业机械的第三位置值和第三航向角度值;目标定位参数获取装置55用于对所述第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数,并将所述第二定位参数与所述第三定位参数进行融合,获取目标定位参数;定位装置56用于根据所述目标定位参数对所述农业机械进行定位。
目标定位参数获取装置55包括卡尔曼滤波器551和多传感器自适应加权融合单元552,其中卡尔曼滤波器551分别与DGPS接收机51、电子罗盘52、加速度计53连接,多传感器自适应加权融合单元552分别与摄像装置54和卡尔曼滤波器551连接。
卡尔曼滤波器551用于对所述第一位置值、第一航向角度值和当前行进速度值进行融合;多传感器自适应加权融合单元552用于将所述第二定位参数与所述第三定位参数进行融合。
本发明采用卡尔曼滤波器融合DGPS得到的农业机械的当前位置值和电子罗盘得到的航向角度值及通过对加速度计得到的初始值积分得到的速度值,可平滑DGPS定位数据,有效避免DGPS动态定位的异常结果。另外,本发明提出采用多传感器自适应加权融合算法将机器视觉提供的位置值和航向角度值与初次融合得到相对精确的农业机械的位置和航向角度信息进行再次融合,获得最优的农业机械位置值和航向角度值。进一步,本发明的农业机械导航分级定位方法,将各个传感器的测量数据进行多次融合,可以有效滤除测试噪声、减小系统误差,能够形成连续、稳定、相对精确的农业机械的位置和航向角度信息
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1、 一种农业机械导航分级定位的方法,其特征在于,所述方法包括以下步骤:
获取第一定位参数,所述第一定位参数包括通过差分全球定位系统DGPS接收机获取的农业机械的第一位置值,通过电子罗盘获取的所述农业机械的第一航向角度值,通过对从加速度计得到的加速度值进行积分而获取的所述农业机械的当前行进速度值;
利用卡尔曼滤波器对所述第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数,所述第二定位参数包括所述农业机械的第二位置值和第二航向角度值;
获取第三定位参数,所述第三定位参数包括通过摄像装置获取的所述农业机械的第三位置值和第三航向角度值;
利用多传感器自适应加权融合算法将所述第二定位参数与所述第三定位参数进行融合,获取目标定位参数,所述目标定位参数包括所述农业机械的目标位置值和目标航向角度值;
根据所述目标定位参数对所述农业机械进行定位。
2、 如权利要求1所述农业机械导航分级定位的方法,其特征在于,在所述利用卡尔曼滤波器对第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数之前,还包括构建卡尔曼滤波器,所述构建卡尔曼滤波器的步骤具体包括:
建立基于常速度的卡尔曼滤波器的状态方程和观测方程;
根据所述状态方程和观测方程,建立预测方程组和校正方程组;
根据所述预测方程组和校正方程组构建卡尔曼滤波器。
4、 如权利要求2所述农业机械导航分级定位的方法,其特征在于,所述预测方程组为:
其中,为t时刻农业机械的运动状态,Φ(t-1)为t-1时刻状态转移矩阵,为t-1时刻农业机械的运动状态,P(t|t)为t时刻滤波误差方差阵,φ(t|t)为t时刻状态转移矩阵,P(t-1|t-1)为t-1时刻滤波误差方差阵,φT(t|t)为t时刻状态转移矩阵的转置矩阵,G(t-1)为t-1时刻过程噪声分布矩阵,Q(t-1)为t-1时刻正定协方差矩阵,GT(t-1)为t-1时刻过程噪声分布矩阵的转置矩阵;
所述校正方程组为:
5、 如权利要求2至4任一项所述农业机械导航分级定位的方法,其特征在于,在所述卡尔曼滤波器中,利用DGPS接收机得到的第一位置值对所述状态转移矩阵进行实时更新。
6、 如权利要求2至4任一项所述农业机械导航分级定位的方法,其特征在于,所述利用卡尔曼滤波器对第一位置值、第一航向角度值和当前行进速度值进行融合,具体包括:
利用所述第一位置值、第一航向角度值和当前行进速度值更新所述预测方程组的矩阵;
通过所述预测方程组获取下一时刻的状态估计;
利用所述校正方程组对所述下一时刻的状态估计进行校正,获取当前最优估计;
根据所述当前最优估计,获取当前最优估计的误差协方差矩阵。
7、 如权利要求1所述农业机械导航分级定位的方法,其特征在于,在所述利用多传感器自适应加权融合算法将第二定位参数与第三定位参数进行融合,获取目标定位参数之前,还包括以下步骤:
判断所述DGPS接收机和摄像装置是否均能正常工作,
如果所述DGPS接收机和摄像装置都不能正常工作,则停止导航;
如果所述DGPS接收机和摄像装置中只有一个能正常工作,则判断能正常工作的装置得到的定位参数是否在误差范围内,如果是则采用单一传感器进行导航,否则停止导航;
如果所述DGPS接收机和摄像装置都能正常工作,则判断能正常工作的装置得到的定位参数是否在误差范围内,
如果所述DGPS接收机和摄像装置得到的定位参数都不在误差范围内,则停止导航;
如果所述DGPS接收机和摄像装置得到的定位参数中只有一个在误差范围内,则采用单一传感器进行导航;
如果所述DGPS接收机和摄像装置得到的定位参数都在误差范围内,则利用多传感器自适应加权融合算法将所述DGPS接收机得到的第二定位参数和摄像装置得到的第三定位参数进行融合。
8、 如权利要求1所述消息的读写方法,其特征在于,所述利用多传感器自适应加权融合算法将第二定位参数与第三定位参数进行融合,具体包括:
根据各测量设备获取的定位参数的方差,获取各测量设备的最优加权因子,所述最优加权因子由公式
获取;其中,Wp为最优加权因子,σ2为各测量设备获取的定位参数的方差;
根据公式
获取融合后的定位参数;其中,为融合后的定位参数,δq 2为摄像装置获取的定位参数的方差,Xp为DGPS接收机得到的定位参数值,δp 2为DGPS接收机获取的定位参数的方差,Xq为摄像装置得到的定位参数值。
9、 一种农业机械导航分级定位的系统,其特征在于,所述系统包括:
DGPS接收机,用于获取农业机械的第一位置值;
电子罗盘,用于获取所述农业机械的第一航向角度值;
加速度计,用于获取所述农业机械的加速度值;
摄像装置,用于获取所述农业机械的第三位置值和第三航向角度值;
目标定位参数获取装置,用于对所述第一位置值、第一航向角度值和当前行进速度值进行融合,获取第二定位参数,并将所述第二定位参数与所述第三定位参数进行融合,获取目标定位参数;
定位装置,用于根据所述目标定位参数对所述农业机械进行定位。
10、 如权利要求9所述农业机械导航分级定位的系统,其特征在于,所述目标定位参数获取装置包括:
卡尔曼滤波器,用于对所述第一位置值、第一航向角度值和当前行进速度值进行融合;
多传感器自适应加权融合单元,用于将所述第二定位参数与所述第三定位参数进行融合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101135922A CN101285686B (zh) | 2008-05-29 | 2008-05-29 | 一种农业机械导航分级定位的方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008101135922A CN101285686B (zh) | 2008-05-29 | 2008-05-29 | 一种农业机械导航分级定位的方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101285686A true CN101285686A (zh) | 2008-10-15 |
CN101285686B CN101285686B (zh) | 2010-11-10 |
Family
ID=40058004
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008101135922A Expired - Fee Related CN101285686B (zh) | 2008-05-29 | 2008-05-29 | 一种农业机械导航分级定位的方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101285686B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102297692A (zh) * | 2011-07-12 | 2011-12-28 | 重庆邮电大学 | 用于智能轮椅在转角区域的自定位方法 |
CN102591332A (zh) * | 2011-01-13 | 2012-07-18 | 同济大学 | 用于无人驾驶汽车局部路径规划的装置及方法 |
CN101846734B (zh) * | 2009-03-26 | 2012-08-22 | 中国农业大学 | 农用机械导航定位方法、系统及农用机械工控机 |
CN102667855A (zh) * | 2009-10-19 | 2012-09-12 | Metaio有限公司 | 用于确定相机姿态以及用于识别真实环境的对象的方法 |
WO2013155919A1 (zh) * | 2012-04-16 | 2013-10-24 | 中兴通讯股份有限公司 | 一种定位方法和系统 |
CN103838210A (zh) * | 2014-02-25 | 2014-06-04 | 北京理工大学 | 突发事件现场有毒气体远程无线监测系统和方法 |
CN103929716A (zh) * | 2014-04-24 | 2014-07-16 | 黄卿 | 一种定位方法、定位信息发送方法及装置 |
CN104778722A (zh) * | 2015-03-20 | 2015-07-15 | 北京环境特性研究所 | 一种传感器的数据融合方法 |
CN104990563A (zh) * | 2015-07-22 | 2015-10-21 | 广西大学 | 车辆行驶里程计算方法 |
CN105374224A (zh) * | 2015-10-29 | 2016-03-02 | 深圳市元征科技股份有限公司 | 一种定位数据处理方法及车载终端 |
CN107576977A (zh) * | 2017-09-08 | 2018-01-12 | 北京韦加无人机科技股份有限公司 | 基于多源信息自适应融合的无人机导航系统及方法 |
CN107782304A (zh) * | 2017-10-26 | 2018-03-09 | 广州视源电子科技股份有限公司 | 移动机器人的定位方法及装置、移动机器人及存储介质 |
CN108332750A (zh) * | 2018-01-05 | 2018-07-27 | 深圳市功夫机器人有限公司 | 机器人定位方法及终端设备 |
CN108931806A (zh) * | 2017-05-25 | 2018-12-04 | 中国辐射防护研究院 | 一种用于无线定位个人辐射剂量仪的定位坐标算法 |
CN109883423A (zh) * | 2019-02-25 | 2019-06-14 | 广州市香港科大霍英东研究院 | 基于卡尔曼滤波的定位方法、系统、设备及存储介质 |
CN110430534A (zh) * | 2019-09-12 | 2019-11-08 | 北京云迹科技有限公司 | 一种定位选择方法、装置、电子设备及存储介质 |
CN111086001A (zh) * | 2019-12-25 | 2020-05-01 | 广东省智能制造研究所 | 一种足式机器人多模态感知的状态估计方法及系统 |
US10739142B2 (en) | 2016-09-02 | 2020-08-11 | Apple Inc. | System for determining position both indoor and outdoor |
CN112629529A (zh) * | 2020-12-15 | 2021-04-09 | 西安工业大学 | 无人机室内自主导航方法 |
CN112974068A (zh) * | 2021-02-11 | 2021-06-18 | 新疆农业大学 | 具有地面仿形功能的喷雾机喷杆系统 |
CN117760412A (zh) * | 2023-08-16 | 2024-03-26 | 北京路凯智行科技有限公司 | 分级更新的矿区高精地图制作装置 |
-
2008
- 2008-05-29 CN CN2008101135922A patent/CN101285686B/zh not_active Expired - Fee Related
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846734B (zh) * | 2009-03-26 | 2012-08-22 | 中国农业大学 | 农用机械导航定位方法、系统及农用机械工控机 |
US10650546B2 (en) | 2009-10-19 | 2020-05-12 | Apple Inc. | Method of providing a descriptor for at least one feature of an image and method of matching features |
US8942418B2 (en) | 2009-10-19 | 2015-01-27 | Metaio Gmbh | Method of providing a descriptor for at least one feature of an image and method of matching features |
CN102667855A (zh) * | 2009-10-19 | 2012-09-12 | Metaio有限公司 | 用于确定相机姿态以及用于识别真实环境的对象的方法 |
US10229511B2 (en) | 2009-10-19 | 2019-03-12 | Apple Inc. | Method for determining the pose of a camera and for recognizing an object of a real environment |
US10580162B2 (en) | 2009-10-19 | 2020-03-03 | Apple Inc. | Method for determining the pose of a camera and for recognizing an object of a real environment |
CN106650723A (zh) * | 2009-10-19 | 2017-05-10 | Metaio有限公司 | 用于确定相机姿态以及用于识别真实环境的对象的方法 |
US10062169B2 (en) | 2009-10-19 | 2018-08-28 | Apple Inc. | Method of providing a descriptor for at least one feature of an image and method of matching features |
US9218665B2 (en) | 2009-10-19 | 2015-12-22 | Metaio Gmbh | Method for determining the pose of a camera and for recognizing an object of a real environment |
US8837779B2 (en) | 2009-10-19 | 2014-09-16 | Metaio Gmbh | Method for determining the pose of a camera and for recognizing an object of a real environment |
CN102591332A (zh) * | 2011-01-13 | 2012-07-18 | 同济大学 | 用于无人驾驶汽车局部路径规划的装置及方法 |
CN102591332B (zh) * | 2011-01-13 | 2014-08-13 | 同济大学 | 用于无人驾驶汽车局部路径规划的装置及方法 |
CN102297692A (zh) * | 2011-07-12 | 2011-12-28 | 重庆邮电大学 | 用于智能轮椅在转角区域的自定位方法 |
CN103379619B (zh) * | 2012-04-16 | 2017-11-28 | 中兴通讯股份有限公司 | 一种定位方法和系统 |
CN103379619A (zh) * | 2012-04-16 | 2013-10-30 | 中兴通讯股份有限公司 | 一种定位方法和系统 |
WO2013155919A1 (zh) * | 2012-04-16 | 2013-10-24 | 中兴通讯股份有限公司 | 一种定位方法和系统 |
CN103838210B (zh) * | 2014-02-25 | 2017-01-04 | 北京理工大学 | 突发事件现场有毒气体远程无线监测系统和方法 |
CN103838210A (zh) * | 2014-02-25 | 2014-06-04 | 北京理工大学 | 突发事件现场有毒气体远程无线监测系统和方法 |
CN103929716A (zh) * | 2014-04-24 | 2014-07-16 | 黄卿 | 一种定位方法、定位信息发送方法及装置 |
CN104778722A (zh) * | 2015-03-20 | 2015-07-15 | 北京环境特性研究所 | 一种传感器的数据融合方法 |
CN104990563B (zh) * | 2015-07-22 | 2017-12-26 | 广西大学 | 车辆行驶里程计算方法 |
CN104990563A (zh) * | 2015-07-22 | 2015-10-21 | 广西大学 | 车辆行驶里程计算方法 |
CN105374224B (zh) * | 2015-10-29 | 2018-04-10 | 深圳市元征科技股份有限公司 | 一种定位数据处理方法及车载终端 |
CN105374224A (zh) * | 2015-10-29 | 2016-03-02 | 深圳市元征科技股份有限公司 | 一种定位数据处理方法及车载终端 |
US11859982B2 (en) | 2016-09-02 | 2024-01-02 | Apple Inc. | System for determining position both indoor and outdoor |
US10739142B2 (en) | 2016-09-02 | 2020-08-11 | Apple Inc. | System for determining position both indoor and outdoor |
CN108931806A (zh) * | 2017-05-25 | 2018-12-04 | 中国辐射防护研究院 | 一种用于无线定位个人辐射剂量仪的定位坐标算法 |
CN108931806B (zh) * | 2017-05-25 | 2024-04-12 | 中国辐射防护研究院 | 一种用于无线定位个人辐射剂量仪的定位坐标算法 |
CN107576977A (zh) * | 2017-09-08 | 2018-01-12 | 北京韦加无人机科技股份有限公司 | 基于多源信息自适应融合的无人机导航系统及方法 |
CN107576977B (zh) * | 2017-09-08 | 2020-11-03 | 北京韦加无人机科技股份有限公司 | 基于多源信息自适应融合的无人机导航系统及方法 |
CN107782304A (zh) * | 2017-10-26 | 2018-03-09 | 广州视源电子科技股份有限公司 | 移动机器人的定位方法及装置、移动机器人及存储介质 |
CN108332750A (zh) * | 2018-01-05 | 2018-07-27 | 深圳市功夫机器人有限公司 | 机器人定位方法及终端设备 |
CN109883423A (zh) * | 2019-02-25 | 2019-06-14 | 广州市香港科大霍英东研究院 | 基于卡尔曼滤波的定位方法、系统、设备及存储介质 |
CN110430534A (zh) * | 2019-09-12 | 2019-11-08 | 北京云迹科技有限公司 | 一种定位选择方法、装置、电子设备及存储介质 |
CN111086001A (zh) * | 2019-12-25 | 2020-05-01 | 广东省智能制造研究所 | 一种足式机器人多模态感知的状态估计方法及系统 |
CN112629529A (zh) * | 2020-12-15 | 2021-04-09 | 西安工业大学 | 无人机室内自主导航方法 |
CN112974068B (zh) * | 2021-02-11 | 2021-11-26 | 新疆农业大学 | 具有地面仿形功能的喷雾机喷杆系统 |
CN112974068A (zh) * | 2021-02-11 | 2021-06-18 | 新疆农业大学 | 具有地面仿形功能的喷雾机喷杆系统 |
CN117760412A (zh) * | 2023-08-16 | 2024-03-26 | 北京路凯智行科技有限公司 | 分级更新的矿区高精地图制作装置 |
CN117760412B (zh) * | 2023-08-16 | 2024-06-07 | 北京路凯智行科技有限公司 | 分级更新的矿区高精地图制作装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101285686B (zh) | 2010-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101285686B (zh) | 一种农业机械导航分级定位的方法和系统 | |
CN101082493B (zh) | 一种农业机械导航的组合定位方法 | |
CN104406605B (zh) | 机载多导航源综合导航仿真系统 | |
Gu et al. | GNSS/onboard inertial sensor integration with the aid of 3-D building map for lane-level vehicle self-localization in urban canyon | |
CN102508278B (zh) | 一种基于观测噪声方差阵估计的自适应滤波方法 | |
Jing et al. | Integrity monitoring of GNSS/INS based positioning systems for autonomous vehicles: State-of-the-art and open challenges | |
CN105509738B (zh) | 基于惯导/多普勒雷达组合的车载定位定向方法 | |
US8775063B2 (en) | System and method of lane path estimation using sensor fusion | |
CN104729506A (zh) | 一种视觉信息辅助的无人机自主导航定位方法 | |
CN104390646A (zh) | 水下潜器地形辅助惯性导航系统的位置匹配方法 | |
CN102353378A (zh) | 一种矢量形式信息分配系数的自适应联邦滤波方法 | |
Deng et al. | Large-scale navigation method for autonomous mobile robot based on fusion of GPS and lidar SLAM | |
CN102928858A (zh) | 基于改进扩展卡尔曼滤波的gnss单点动态定位方法 | |
Mu et al. | A GNSS/INS-integrated system for an arbitrarily mounted land vehicle navigation device | |
CN104677361A (zh) | 一种综合定位的方法 | |
Dawson et al. | Radar-based multisensor fusion for uninterrupted reliable positioning in GNSS-denied environments | |
Meichen et al. | Dynamic obstacle detection based on multi-sensor information fusion | |
CN106370181A (zh) | 一种高精度的组合导航定位方法 | |
US11127304B2 (en) | Method, device, and computer-readable storage medium with instructions for estimating the pose of a transportation vehicle | |
Ghanem et al. | Testing vehicle-to-vehicle relative position and attitude estimation using multiple UWB ranging | |
CN113063441B (zh) | 里程计累计推算误差的数据源纠正方法及装置 | |
Song et al. | RFID/in-vehicle sensors-integrated vehicle positioning strategy utilising LSSVM and federated UKF in a tunnel | |
Cahyadi et al. | Unscented Kalman filter for a low-cost GNSS/IMU-based mobile mapping application under demanding conditions | |
Al Hage et al. | Bounding localization errors with Student distribution for road vehicles | |
Gu et al. | Integration of 3D map based GPS positioning and on-board sensors for vehicle self-localization in urban canyon |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101110 Termination date: 20130529 |