[发明专利]面向地震应急响应的多无人机调度与任务规划方法

权利要求书

1.一种面向地震应急响应的多无人机调度与任务规划方法，其特征在于，包括以下步骤：

步骤1：根据震源坐标、深度、震级、烈度参数，生成地震灾区范围，并结合居民地、道路、水系和数字高程模型数据，生成灾区重点监测区域，即任务区；

步骤2：利用路径分析方法，计算可供调遣的各无人机单位到达灾区时间，根据计算得到的时间选取前往灾区的无人机单位，即得到无人机单位调度方案；

步骤3：无人机单位根据无人机单位调度方案向灾区派遣无人机编队，然后根据所要派遣的无人机编队位置与状态，生成无人机编队前往灾区的交通方案；

步骤4：无人机编队按照交通方案前往灾区，然后根据所有派遣的无人机编队的移动情况，生成多无人机编队任务规划方案，并根据无人机编队实时的作业和交通状态，动态规划作业方案；

所述的地震灾区范围是指利用地震烈度椭圆经验公式计算出的地震烈度达到指定等级的椭圆区域；

所述的任务区域是灾区内部需监测的重点区域，包括居民居住聚集区域、陡峭路段和水系区域；

所述的路径分析方法采用最短路径分析方法，通过最短路径分析方法计算给定无人机编队当前位置到灾区、任务区的最佳路径及其最小通行时间；

所述的无人机编队是指拥有无人机的单位下属的最小作业单元，包括操作人员、无人机和数据采集部件；

所述的数字高程模型数据即DEM，是利用栅格数据结构表达地表高程的数据模型。

2.根据权利要求1所述的一种面向地震应急响应的多无人机调度与任务规划方法，其特征在于，步骤1所述的生成任务区包括以下步骤：

1)根据设定的行、列数，将灾区范围划分成若干个正交矩形单元格；

2)计算每个单元格包含或相交的居民地多边形要素，进而计算各单元格居民地密度，依次将该单元格上下左右4方向相互邻近且密度大于指定阈值的单元格所包含或相交的居民地多边形要素聚为一个任务区域；进而采用约束Delaunay三角网法计算同一任务区域居民地要素聚合边界，即得到居民地任务区；

3)利用DEM数据生成灾区坡度图，然后根据坡度大小对坡度图进行重分类，从中选取超出指定坡度值的区域；最后计算这些区域中的道路、水系范围，生成道路、水系的坡度任务区；

所述的居民地要素是指在地图上使用多边形方式表达的房屋等人工建筑实体。

所述的约束Delaunay三角网是指同一任务区域中所有居民地多边形顶点构造的且不与居民地多边形相交的Delaunay三角网；

所述的聚合边界是利用约束Delaunay三角网中包含同一任务区域所有居民地多边形的最外侧三角形边线所围成的多边形；

所述的重分类是在地理信息统中利用给定的属性和分类区间值域，从而将空间或空间实体划分为多个部分的一种方法。

所述的任务区包括步骤2)和3)中生成的居民地任务区和坡度任务区。

3.根据权利要求1中所述的一种面向地震应急响应的多无人机调度与任务规划方法，其特征在于，步骤2所述确定无人机单位调度方案，包括以下处理步骤：

1)结合全国路网数据，利用路径分析方法计算无人机单位所在地到灾区的交通时间，并生成交通时间报表；

2)选择交通时间满足应急响应要求的无人机单位并确定无人机单位调度方案，并生成无人机单位订单报表。

所述的路网数据是指将全国道路数据经数据加工、拓扑处理后的道路网络数据集。

4.根据权利要求1中所述的一种面向地震应急响应的多无人机调度与任务规划方法，其特征在于，步骤3所述的跟踪无人机编队并生成交通方案报表，包括以下处理过程：

1)无人机编队内的无人机上传编队实时信息，实时信息包括实时位置、设备状态、工作状态以及突发事件；

2)根据实时信息，重新计算各无人机编队到达灾区的交通时间，并生成交通方案报表；

步骤1)中所述的突发事件包括交通中断和无人机损毁。

步骤2)中所述的重新计算各无人机编队到达灾区的交通时间是根据无人机编队实时信息中的交通中断信息进行计算。

5.根据权利要求1中所述的一种面向地震应急响应的多无人机调度与任务规划方法，其特征在于，步骤4所述的多无人机编队任务规划方案，包括以下处理过程：

1)利用路径分析工具，以当前各无人机编队所在位置为起始结点，以无人机编队到达最近任务区以及任务区之间通行时间为权值，生成无人机-任务区覆盖路径，并在每个结点上标注无人机编队和到达任务区可能花费的时间；

2)将所有任务区覆盖路径结点，按时间点由早到晚排序，然后从排序列表中剔除重复且不是最早结点上的任务区，最后得到多无人机任务规划方案；

3)实时跟踪多无人机编队交通、任务执行状态，并在方案执行过程中，依据无人机编队上报信息，动态规划后续任务执行方案；

所述的覆盖路径是从无人机编队位置开始遍历所有任务区所需交通时间成本的一个连通子图；本发明使用的覆盖路径是从无人机编队所在位置出发，查找交通时间成本最小的任务区，再从该任务区出发，寻找距离该任务交通时间最小的且未归入覆盖路径的任务区，依此递归搜索，直到所有任务区均加入覆盖路径时结束计算。