基于机器学习的小麦收获机掉头轨迹识别.pdf

来戏第54卷第9 期2023年9 月农报学业机械doi:10.6041/j.issn.1000-1298.2023.09.003基于机器学习的小麦收获机掉头轨迹识别杨丽丽1,2王新鑫1李元博1常孟帅！翟卫欣1，2吴才聪1.2（1.中国农业大学信息与电气工程学院，北京10 0 0 8 3；2.农业农村部农机作业监测与大数据应用重点实验室，北京10 0 0 8 3）摘要：识别小麦收获机运行轨迹是分析农业机械活动、提高作业效率的重要手段。本文针对小麦收获机田内作业场景，提出一种基于机器学习的收获机掉头轨迹识别算法。首先通过两步K-means聚类与三步修正识别出X形掉头轨迹点、作业异常轨迹点与作业轨迹点；为进一步从作业轨迹中分类出U形掉头轨迹点，构建了基于支持向量机模型（Supportvectormachine,SVM）的U形掉头轨迹识别算法，并对初步识别结果进行三步修正；最终识别出小麦收获机的田内X形掉头、作业异常、U形掉头与作业轨迹点，识别结果的F1值为94%，时间间隔为1 5s的数据的F1值在90%以上，实现田内轨迹的细致划分。基于去除掉头轨迹与异常轨迹后获得的有效作业轨迹，可通过距离算法计算获得农田面积，结果相比使用原始轨迹的计算误差可降低12.7 6%。该研究可为基于海量农机轨迹的作业精细化管理提供参考。关键词：小麦收获机；轨迹识别；机器学习；掉头；SVM中图分类号：S126文献标识码：A文章编号：10 0 0-12 9 8（2 0 2 3)0 9-0 0 2 7-0 8OSID:Identifying Turning Trajectories of Wheat Harvester Based onMachine LearningYANG Lili,2WANG XinxinLI YuanbolCHANG MengshuaiZHAI Weixin.2WU Caicong1,2(1.College of Information and Electrical Engineering,China Agricultural University,Beijing 100083,China2.Key Laboratory of Agricultural Machinery Monitoring and Big Data Application,Ministry of Agriculture and Rural Affairs,Beijing 100083,China)Abstract:Identifying the trajectories of wheat harvester in the field is an important means to analyze theactivities of agricultural machinery and improve the working efficiency.A machine learning basedalgorithm for recognizing the turning trajectories of wheat harvester was proposed.Identifing X-turn,abnormal working,and working trajectory through two-step K-means iterative clustering and three-stepcorrection method:the first step(M1)was performed based on the three distance features between thetrajectory segments and the cluster center of the trajectory segments.The second step(M2)was based onthe direction change of the“turning and“abnormal working trajectories.The third correction step(M3)was based on the operating characteristics to specify the start and stop positions of the turning.Inorder to further classify U-turn trajectories from working trajectories,identifying X-turn,abnormalworking,U-turn and working trajectories through SVM model and three-step correction method,firstly,the correction of U-turn boundary based on trajectory curvature(S1)was carried out.Secondly,basedon the time difference between X-turn and U-turn,the misidentification as a U-turn was corrected(S2).Thirdly,the correction was based on the change of the angle before and after the U-turn(S3).The F1-score of the four trajectories recognition results was 94%.The accuracy,recall,and F1 scores of datarecognition results at different time intervals of 1 5 s were all above 85%,indicating that the algorithmperformed well on trajectory data at 1 5 s intervals.When the time interval was extended to 10 s and15 s,the U-turn trajectory would not be recognized,indicating that the algorithm cannot be applied tooverly sparse trajectory data.The effective working trajectories were obtained after removing the X-turn收稿日期：2 0 2 3-0 3-2 00修回日期：2 0 2 3-0 5-2 2基金项目：国家精准农业应用项目（JZNYYYO01）作者简介：杨丽丽（197 4一），女，副教授，博士，主要从事计算机网络与智能信息处理研究，E-mail：l l y a n g c a u.e d u.c n通信作者：吴才聪（197 5一），男，教授，博士，主要从事农机作业大数据挖掘和农机导航与位置服务信息研究，E-mail：w u c d u.c n2023年农28机报学械业trajectories,U-turn trajectories and abnormal working trajectories of the positioning track data in a field.The error of calculating the farmland area by the distance algorithm can be reduced by 12.76%comparedwith the calculation error of using the original data.The research result can provide a reference for finemanagement of farmland operations.Key words:wheat harvester;trajectory identify;machine learning;U-turn;SVM0引言GNSS全球导航卫星系统是获取车辆轨迹数据、实现农业机械智能化管理的重要定位设备1-3。在农业机械上搭载定位设备可获得农机实时定位轨迹点，精确获取并识别农机运动轨迹，对农业机械自动化管理具有重要意义。农机轨迹数据反映农机的行为特征和田路分布等深层次信息，是农机作业调度管理4、农机作业行为分析5、农机手驾驶评价6 、农田与机耕道轨迹识别7 等研究的重要数据来源。掉头识别是提升农机工作效率、合理规划路径的重要手段之一。设计合适的农机行驶路线和掉头方式与农机的工作效率息息相关，农机作业路径规划的主要目标是高效地实现农田全覆盖作业8 ,掉头识别是路径优化的关键点之一3。相关研究表明，掉头时间可占农机在田内作业总时间的40%左右9,且U形掉头效率较高【10】,通过农机定位轨迹信息识别出农机掉头方式，进而计算出不同掉头方式下的农机作业效率，便于机手合理规划路径，提升农机作业效率。另外，掉头方式与田块形状关系密切8.1,某地区的农田中不同掉头占比可作为判断该区域农田形状是否合理的依据之一，为农田区域的合理规划提供参考。农田面积的计算也是基于农机运动轨迹的相关研究之一，距离算法是计算农田面积的常用算法12-13，对农机掉头轨迹识别后可以去除农机掉头行为导致的交错轨迹，防止面积计算过程中出现同一块区域面积的重复累加，从而提升距离算法计算农田面积的精度14-16小麦是全球粮食主产物之一17-18 ,在我国是仅次于水稻、玉米的主要粮食作物，研究小麦收获轨迹对精准农业的发展至关重要。小麦收获轨迹包含典型的X形掉头与U形掉头轨迹。本文对小麦收获机田内的X形掉头、作业异常、U形掉头与作业轨迹进行识别。对收获机田内轨迹进行细致的划分，以期为农机作业效率计算、农田面积计算和农机路径规划等研究提供参考。1数据与方法1.1数据预处理数据来自北斗农机作业大数据系统12 ，获取的数据为2 0 2 2 年6 月小麦收获季在河北、河南、山东等小麦主产区通过装载全球导航卫星系统（Globalnavigation satellite system，G NSS）的收获机，采集的产品为幅宽为2.7 5m、车身长度为6.8 m的4LZ-8E2型和4LZ-7E5型自走式谷物联合收获机小麦收获轨迹,并通过人工标注获取田内农机轨迹。每条GNSS记录包含4个参数：时间（记为t，格式为YYYY-MM-DD hh:mm:ss)、经度(World geodeticsystem，W G S8 4世界大地坐标系）、纬度（WGS84）和速度（记为u,单位：m/s）。选取时间间隔5s占比均在8 5%以上的5块已人工标注的农田，田内轨迹共1969条，作为数据集A，用于训练后续试验中的SVM模型。选取时间间隔15s占比均在8 5%以上的数据，每种时间间隔选取10 块农田内的农机轨迹数据共50 块，10 6 0 58 条数据作为数据集B，用来测试算法效果。数据集见表1。表1数据集Tab.1Datasetsituation时间间隔/s数据集总计12345A19691969B44.3072215519651 110808 865106058数据预处理包含数据格式转换、去除停歇点和去除重复点。格式转换：将原始WGS84坐标系下的经纬度转换为平面坐标系下的x、y,转换后轨迹点P：属性为（t，i，y i，），时间t；格式为YYYYMMDDhhmmss。去除停歇点：将速度小于0.5m/s的连续轨迹点视为一组停歇点5，一组停歇点只保留一个轨迹点，其t为这组停歇点中第1个轨迹点的t,x、为这组停歇点、y 的均值，u设为0。去除重复点：当存在连续轨迹点的x、y 相同时，保留第1个轨迹点，将其他轨迹点删除。1.2数据定义小麦收获机作业中的典型路线如图1所示，收获作业从农田外沿开始，绕圈向内进行收获，其形状如“回”字，定义为回形轨迹。图1a为收获机开始收获路线。收获机开人农田后沿麦田外围走直线收获至拐角处，为完成转向同时收获边角小麦需进行多次前进与倒车行为19当农田四周割出5m左右的割道后进人内圈收29杨丽丽等：基于机器学习的小麦收获机掉头轨迹识别第9 期前进方向-倒车方向割道21(a)(b)图1回形收割路径示意图Fig.1Schematics of harvesting path获2 0 ,内圈收获路线如图1b所示,空白部分为已收获完的割道，阴影部分为尚未收获的作物区域，此时收获机将以U形掉头方式进行转向,以提高掉头效率10 。O为区分不同轨迹类型,轨迹点中增加label属性。如图1所示,回形小麦收获轨迹中常采用X形与U形掉头方式19。X形掉头是指农机在田内作业过程中为转变行进方向而做出的“前进-倒车-前进”的多次进退行为，该类轨迹点对应label属性记为1,如图1a的虚线框1处；作业异常轨迹是指农机行驶时脱离原轨迹方向，发生绕行或其他行为，之后又回到原轨迹方向，该类轨迹点对应label属性记为2,如图1a的虚线框2 处；U形掉头是指无需倒车直接转向的行为，该类轨迹点对应label属性记为3,如图1b的虚线框3处；除掉头轨迹与作业异常轨迹外，其他轨迹为作业轨迹，该类轨迹点对应label属性记为0,如图1a虚线框0 处。本文针对上述4种轨迹对其进行识别。为描述农机行驶方向变化，记时间序列上相邻的3个轨迹点为P;-1、P,与P;+1,记Pi-1到P；的向量为lps-P。定义轨迹点P,的角度属性;的计算式为;=arccos(lpi-tp,lp,Pitt/Ilpi-1P,IlpPi+,1)(1)通过聚类得到的较大点称之为拐点，按时间序列计算每个轨迹点到其最近的拐点的欧氏距离，记为轨迹点P；的D_T，属性。将每个轨迹点与其时间序列上的前一个轨迹点的欧氏距离记为轨迹点P，的D_P，属性。以P:-1、P,与P+1为一组,计算每组的时间差值t;、加速度a;、速度平均值、角度平均值;、角度标准差、角度最大值max:、角度极差R、平均角度差，计算式为At;=ti+1-t;-1(2)U;+1-U;-1a;=(3)ti+1-t;-1V;-1+U;+U;+1U;=(4)3;-1+;+;+1;=(5)3(;-1-,)+(,-,)+(;+1-,)3(6)m;=max(i-1,i,i+)(7)Ra,=max(i-1,i;,i+1)-min(i-1,i,i+1)(8)iI;-;-1 I+I+1-,IA;=(9)21.3技术路线本文技术路线如图2 所示，先对农机GNSS数据进行预处理，再分别用X形掉头识别算法与U形掉头识别算法识别出4种轨迹。最后用数据集B中的50 块农田轨迹数据进行算法验证，并用距离算法比较原始轨迹数据与去除掉头和异常轨迹后数据的面积计算其精度。2掉头轨迹识别2.1X形掉头轨迹识别X形掉头轨迹识别算法包含2 个模块：聚类模K-means一次聚类SVMK-means二次聚类数据格式转换D-K-means测试集S1:基于掉头轨迹边界曲率的修正去除停歇点M1:基于轨迹片段的聚类修正1S2:基于两种掉头面积计算时间差的修正去除重复点M2:基于方向变化的修正S3:基于掉头前后角度变化的修正M3:基于掉头路线特性的修正修正修正数据预处理X形掉头识别U形掉头识别验证评价1图2技术路线Fig.2Technology roadmap农2023年机报学械业30块与修正模块。聚类模块执行2 次K-means聚类,第1次聚类的输人特征为每个轨迹点的角度属性，设定类别数k为2,初步将轨迹分为拐点与其他轨迹点。基于第1步聚类结果，输人每个轨迹点到最近拐点的距离特征D_T进行第2 步聚类，聚类后初步得到作业轨迹与X形掉头轨迹。基于轨迹片段的聚类修正M1：定义属于同一轨迹类别且在时间序列上相邻的轨迹点为一个轨迹片段,计算每个片段中轨迹点属性的标准差、平均值与最大值作为轨迹片段的特征，对相同轨迹类别的轨迹片段的特征取平均值作为该类别的簇类中心。以轨迹片段为基本单位，计算每个轨迹片段到2个簇类中心的欧氏距离、契比雪夫距离与曼哈顿距离，以投票决策的方式判定该轨迹片段属于哪一类,并将轨迹片段内所有轨迹点的label属性赋予该类别标签。通过此过程，被误识别的轨迹片段得到修正。基于方向变化的修正M2：如图1所示，X形掉头前后农机的行进方向不同，而作业异常轨迹往往是作业中途因避障而出现的绕行行为，故异常轨迹出现前后行进方向相同，根据GNSS的定位误差和小麦收获机行驶速度，定掉头前后行驶方向变化小于10 为方向未发生变化2 1-2 。根据上述特点，对误识别为X形掉头轨迹的作业异常轨迹进行修正。基于掉头路线特性的修正M3：X形掉头轨迹与作业轨迹相邻，为有效识别X形掉头轨迹，需对X形掉头轨迹起止位置进行界定：收获机在进行X形掉头时，为转变行进方向同时不遗漏作业，机身越过已收地至少半个机身位后转向2 0】，本文收获机车身长度为6.8 m,因此将每个X形掉头轨迹片段的起始拐点和终止拐点前后各3.4m的轨迹点归为X形掉头轨迹点，依此完成对X形掉头轨迹边界的界定。2.2U形掉头轨迹识别2.2.1SVM初步识别如图1b所示，在小麦收获到田地中央常采用U形掉头方式，U形掉头轨迹往往与作业轨迹相接，为区分二者,从数据集A中选取人工标注label值为0与3的轨迹点,相同label且时间相邻的3个点为一组构成493条数据。计算每组数据的t、a、i、vmx、R、特征,通过卡方检验进行单因素分析,结果如表2 所示。由表2 可知,与作业轨迹和U形掉头显著强相关的5个特征为&、max、R、Z。表2卡方检验结果Tab.2Chi-square test results特征X4.002a5.08412.3956.361 10303.277 103max1.021 104R。7.709 1034.781 103构建基于SVM的U形掉头轨迹识别模型，以493条数据的5个与U形掉头强相关的特征属性和label属性作为模型输人,将其中8 0%的数据作为训练集。使用通过人工调优和网格搜索法，确定核函数为高斯核函数，确定惩罚参数C为0.2 1,该模型在测试集上的分类准确率为97%，召回率为97%，F1值为97%，可实现作业轨迹和U形掉头轨迹的初步识别，识别结果如图3所示。S1S3前进方向PaP一QPb50S2XPU图3SVM初步识别结果示意图Fig.3Schematic of preliminary recognition results of SVM图3中红色方形为U形掉头轨迹点，黄色三角形为X形掉头轨迹点，蓝色圆形为作业轨迹点。虚框SO处为正确识别的U形掉头轨迹，虚框S1、S2与S3为3种误识别情况。2.2.2基于掉头轨迹曲率的边界修正S1在SVM初步识别的U形掉头轨迹中,存在U形掉头轨迹边界过长或过短问题，如图3的虚框S1。虚框S1中的P1、P2 为U形掉头初步识别结果的两个轨迹边界点。U形掉头轨迹是一个从直线变曲再变直的过程，为每个轨迹点构造曲率特征可以更精确地描述这种变化。(1)特征构建：定义连续3个轨迹点P;-1、P,与Pi+1的外接圆曲率为轨迹点P,的曲率。P,点坐标为（x;,y:），根据3个轨迹点坐标值计算外接圆面积S，计算式为X;-1y;+x,yi+1+X;+1yi-1-yi-Ix;-y;+1-y;x;-1S;=2(10)31杨丽丽等：基于机器学习的小麦收获机掉头轨迹识别第9 期3个轨迹点构成的三角形边长为、b、c,与三角形面积S计算外接圆半径RR.=4S.abc(11)P，曲率计算式为1Curvi(12)二R.(2)设定U形掉头边界点的曲率阈值K：取数据集A中标注的U形掉头轨迹的开始点与结束点的曲率,取最小值作为阈值K。由数据集A计算得K值为0.0 190 m-l。（3)修正：识别出U形掉头轨迹边界点的曲率值大于（或小于）K时，延长（或缩短）U形掉头轨迹至其边界点曲率值刚好小于K。修正后图3虚框S1中边界点由P1、P2 变为PI、P22.2.3基于掉头时间差修正S2因农机在X形掉头前后行驶方向的转变，与X形掉头相邻的作业轨迹段的、max、R。、特征与U形掉头类似,因此在SVM初步识别中会被误识别为U形掉头轨迹，如图3中虚框S2。通过U形掉头轨迹到最近一次的X形掉头轨迹的时间长短来判断该U形掉头轨迹是否为误识别并修正，步骤为：（1)计算最小时间差：对每段U形掉头轨迹片段寻找与其在时间序列上最近的一段X形掉头轨迹片段,计算两段轨迹片段边界点Pu与Px的最小时间差tux，计算式为Atux=min(Itu-txI)(13)式中tuU形掉头轨迹边界点P，的时间属性tx-X形掉头轨迹边界点Px的时间属性(2)设定时间阈值T:取数据集A中标注的U形掉头轨迹tux的最小值，作为时间阈值T。由数据集A计算得T为30 s。（3）修正：当tux大于等于阈值T时，保留此段U形轨迹，否则将此段轨迹点类别值label由3恢复至X形掉头识别后的类别值。修正后图3的虚框S2中U形掉头轨迹点修正为作业轨迹点。2.2.4基于掉头前后角度变化修正S3农机在麦田边角的收获中也会出现作业轨迹被误识别为U形掉头轨迹的情况，如图3中虚框S3处。因此,通过掉头轨迹前后角度的变化来判断是否为误识别并修正，步骤为：（1)计算U形掉头前后向量：U形掉头轨迹的开始点与前一个轨迹点构成向量,U形掉头轨迹的结束点与后一个轨迹点构成的向量b,计算两个向量的夹角（计算方式同式（1），如图3所示。(2)设定角度阈值D：理想状态下U形掉头轨迹前后方向发生18 0 转变2 3,即为18 0。但因实际农田形状不规则，人工操控农业机械以及GNSS定位精度等因素，U形掉头轨迹前后方向变化角度不定。取数据集A中的U形掉头轨迹前后向量夹角的最小值作为角度阈值D。由数据集A计算得D为92.2 133。(3）修正：对通过S1与S2修正后的SVM结果中,计算所有U形掉头轨迹前后的向量夹角。当小于阈值D时，此U形掉头轨迹为误识别,将其label值由3恢复至X形掉头识别算法的结果，当不小于阈值D时，保留此段U形掉头轨迹。如图3的虚框S3中U形掉头轨迹点修正为作业轨迹点。3结果与分析3.1轨迹识别结果选取一块1.7 2 hm的农田，对其进行掉头识别与修正，结果如图4所示。图4a中5个虚线框部分为误识别轨迹，对应修正过程见图4b4f。图4中红色方形为U形掉头轨迹点，蓝色圆形为作业轨迹点，黄色三角形为X形掉头轨迹点，橘黄色菱形为作业异常轨迹点。图4c4e 为X形掉头识别与修正：通过D-K-means聚类方法后得到作业轨迹与掉头轨迹，图4c中c1存在作业轨迹误识别为掉头轨迹的情况，经过基于轨迹片段的聚类修正方法M1修正后见图4c的c2;图4d中d1存在作业异常轨迹被识别为掉头轨迹的情况，经过基于方向变化的修正方法M2后见图4d的d2;图4e中el存在边界不统一的情况，经过基于收获机作业特性的修正M3后见图4e的e2。图4b、4f 为U形掉头识别与修正：通过SVM识别后初步得到U形掉头轨迹，图4b中b1存在边界长度不统一的情况，经基于掉头轨迹曲率的边界修正为图4b的b2；图4b中b2存在X形掉头轨迹相邻的轨迹点被误识别为U形掉头轨迹（黑色虚线圆框处），经基于掉头时间差的修正方法S2修正为图4b的b3；图4f中f1存在作业轨迹被误识别为U形掉头轨迹，经基于掉头前后角度变化S3修正为图 4f 的 f2。3.2轨迹识别算法评价为客观评价收获机田内轨迹识别算法的效果，选取准确率P、召回率R和F1值3种评价指标进行评价。表3为数据集B的轨迹识别算法结果。X形掉头初步识别与三步修正后F1值由55%提高到95%，表明基于聚类的X形掉头识别算法农机322023年报学业械(b)(c)(e)b162d(a)(b)d1C2d2（c）(d)f1f2(e)(f)图4掉头轨迹识别算法分步结果Fig.4Step by step result diagrams of turn trajectory recognition algorithm表3轨迹识别算法结果Tab.3Trajectory recognition algorithm results%算法PRF1值D-K-means536255D-K-means+M1536356D-K-means+M1+M2859488D-K-means+M1+M2+M3939695SVM788681SVM+S1809484SVM+S1+S2909592SVM+S1+S2+S39295943种修正方法可以有效地对收获机田内作业的X形掉头轨迹、异常轨迹与作业轨迹进行识别。在U形掉头轨迹识别中对数据集B的50 幅田内轨迹进行识别,4种轨迹SVM识别结果为8 1%，对SVM初步识别结果进行3步修正，修正前后F1值由8 1%提高到94%，表明X形掉头轨迹识别算法与U形掉头轨迹识别算法结合可以对收获机田内作业的4种轨迹进行识别。实际应用中GNSS设备采集频率不同，选取时间间隔15s的50 幅轨迹，每种时间间隔10 幅，及时间间隔为10、15s的两幅农田内轨迹数据，表4为轨迹识别结果。表4不同频率轨迹数据识别结果Tab.4Identification results of track data at different frequencies%时间作业轨迹X形掉头轨迹作业异常轨迹U形掉头轨迹间隔/sPRF1值PRF1值PRF1值PRF1值199969885979195999789949129998989596959097938898923989898949494859891939593499999995969591989497959559999999397959093919489911099969766997993939315999295489564868686000由表4可见，15s内轨迹识别结果的3种评价指标均在8 5%以上，表示算法在15s时间间隔的轨迹数据中均有良好表现。时间间隔10、15s数据的X形掉头识别结果明显较差，主要原因在于当数据变得稀疏时，轨迹点的连线将从平滑变得尖锐，当原作业轨迹点与U形掉头轨迹点的属性变大时，会被误识别为X形掉头轨迹点。被误识别为X形掉头的U形掉头轨迹点会在SVM初步识别时被去除，导致输人SVM模型的数据中无U形掉头轨迹片段,因此U形掉头轨迹识别结果为0。综上，当轨迹点时间间隔为1 5s时本算法效果较好。33杨丽丽等：基于机器学习的小麦收获机掉头轨迹识别第9 期为进一步验证算法的有效性，从数据集B中选取一块标注面积为1.7 2 hm的农田,其轨迹时间间隔为4s,分别通过距离算法计算轨迹识别前后的农田面积并进行对比。图5a为原始轨迹，包含147 9个轨迹点，距离算法计算面积为2.33hm，算法运行时间为369.67ms；图5b为X形掉头识别结果，包含X形掉头轨迹点157 个，作业异常轨迹点16 个，如图中黄色三角形和橘黄色菱形轨迹点，去除X形掉头与异常轨迹点后计算面积为2.17 hm,算法运行时间为(a)原始轨迹(b）X形掉头识别算法结果(c）U 形掉头识别算法结果(d)处理后轨迹图5收获机田内轨迹处理过程Fig.5Process of harvester in-field trajectory369.08ms，误差相比原始轨迹点面积计算误差降低9.51%；图5c为运行U形掉头识别结果，识别出U形掉头轨迹点47 个，如图中红色方形轨迹点。将两种掉头轨迹及作业异常轨迹点均去除后得到图5d，剩余轨迹点个数12 58 个，计算面积为2.11hm,算法运行时间为36 4.10 ms，面积计算误差相比标注面积仍有2 2.43%，但相比原始轨迹的计算面积误差降低12.7 6%，比只去除X形掉头轨迹和异常轨迹误差降低3.2 6%。4结论（1)通过X形掉头识别算法识别出X形掉头、作业异常和作业轨迹点，再通过U形掉头识别算法将作业轨迹点进一步识别出U形掉头轨迹点。本研究所用的50 块农田轨迹数据的综合识别结果F1值为94%，时间间隔1 5s数据的4种轨迹识别结果F1值均在90%以上。随着时间间隔增加，间隔为10 s与15s的数据识别效果变差，实际应用时可参考1 5s的时间间隔。(2)对一块1.7 2 hm的农田作业轨迹进行识别，去除掉头轨迹与异常轨迹后通过距离算法计算农田面积，相比使用原始轨迹，其面积计算误差降低12.76%。参考文献1WU Caicong,CHEN Zhibo,WANG Dongxu,et al.A cloud-based in-field fleet coordination system for multiple operationsJ.Energies,2020,13(4):775.2WANG Jie,ZHU Yuting,CHEN Zhibo,et al.Auto-steering based precise coordination method for in-field multi-operation offarm machinery J.International Journal of Agricultural and Biological Engineering,2018,11(5):174-181.3翟卫欣，王东旭，陈智博，等。无人驾驶农机自主作业路径规划方法J农业工程学报，2 0 2 1,37（16）：17.ZHAI Weixin,WANG Dongxu,CHEN Zhibo,et al.Autonomous operation path planning method for unmanned agriculturalmachineryJ.Transactions of the CSAE,2021,37(16):1-7.(in Chinese)4黄凰，陈燕燕，陈鹏宇，等.基于时间窗的农机调度技术研究进展J.中国农业科技导报，2 0 2 2,2 4（4）：9310 6.HUANG Huang,CHEN Yanyan,CHEN Pengyu,et al.Research progress of agricultural machinery scheduling technology basedon time window J.Journal of Agricultural Science and Technology,2022,24(4):93-106.(in Chinese)5JENSEN M A F,BOCHTIS D.Automatic recognition of operation modes of combines and transport units based on GNSStrajectoriesJ.IFAC Proceedings Volumes,2013,46(18):213-218.6张萍基于深度学习的农机驾驶员疲劳检测研究D淄博：山东理工大学,2 0 2 2.ZHANG Ping.Research on fatigue detection of agricultural machinery drivers based on deep learning D.Zibo:ShandongUniversity of Technology,2022.(in Chinese)7李亚硕，赵博，王长伟，等基于DBSCAN和BP_Adaboost的农机作业地块划分方法J农业机械学报,2 0 2 3,54(1）：37-44.LI Yashuo,ZHAO Bo,WANG Changwei,et al.Land division method for agricultural machinery operation based on DBSCANand BP_AdaboostJJ.Transactions of the Chinese Society for Agricultural Machinery,2023,54(1):37-44.(in Chinese)8尹彦鑫，孟志军，赵春江，等大田无人农场关键技术研究现状与展望J.智慧农业（中英文）,2 0 2 2,4（4）：1-2 5.YIN Yanxin,MENG Zhijun,ZHAO Chunjiang,et al.State-of-the-art and prospect of research on key technical for unmannedfarms of field corpJ.Smart Agriculture,2022,4(4):1-25.(in Chinese)9WU Caicong,CHEN Zhibo,WANG Dongxu,et al.Behavior modelling and sensing for machinery operations using smartphonessensor data:a case study of forage maize sowing JJ.International Journal of Agricultural and Biological Engineering,2019,12(6):66-74.10姚竟发，滕桂法，霍利民，等.联合收割机多机协同作业路径优化J农业工程学报,2 0 19,35（17）：12-18.农34机2023年报学业械YAO Jingfa,TENG Guifa,HUO Limin,et al.Optimization of cooperative operation path for multiple combine harvesterswithout conflictJJ.Transactions of the CSAE,2019,35(17):12-18.(in Chinese)11GRIFFEL L,VAZHNI K,HARTLEY S,et al.Agricultural field shape descriptors as predictors of field efficiency for perennialgrass harvesting:an empirical proof J.Computers and Electronics in Agriculture,2020,168:105088.12 吴才聪，陈瑛，杨卫中，等基于北斗的农机作业大数据系统构建J农业工程学报，2 0 2 2,38（5）：18.WU Caicong,CHEN Ying,YANG Weizhong,et al.Construction of big data system of agricultural machinery based on BeidouJ.Transactions of the CSAE,2022,38(5):1-8.(in Chinese)13 李亚硕，贾晓峰，伏伟浩，等.基于三次样条插值的农机作业面积计算J农业工程，2 0 2 2，12（3）：43-46.LI Yashuo,JIA Xiaofeng,FU Weihao,et al.Working area calculation of agricultural machinery based on cubic splineinterpolationJ.Agricultural Engineering,2022,12(3):43-46.(in Chinese)14刘卉，孟志军，王培，等基于农机空间轨迹的作业面积的缓冲区算法J.农业工程学报，2 0 15,31（7）：18 0-18 4.LIU Hui,Meng Zhijun,WANG Pei,et al.Buffer algorithms for operation area measurement based on glob