1. 引言
我国地域广阔,地理条件复杂,自然灾害时有发生。自然灾害对国民经济和人民安全造成的影响日益严重,自然灾害不但破坏建筑,造成道路损毁,还带来通信的中断,这使得救援困难,医疗物资不能及时到位,从而导致不必要的伤亡 [1] 。而虽然自然灾害是无法避免的,但通过一个全面和高效率的灾害管理系统,可以大大减轻其影响:时间是灾害管理成功的最重要因素 [2] ,而无人机(Unmanned Aerial Vehicle, UAV)不受地形限制,道路限制,具有低成本、低损耗、可重复使用且风险小等诸多优势 [3] ,可以在第一时间将食品,饮用水,医疗药品等救援物资送达指定区域,并且可以搭载数码设备,实施对地观测,为救援人员提供道路受损信息。
无人机的航迹规划是指在适当的时间内找到满足约束条件的一条可行最优的路径。但是由于路径规划是一个NP-Hard问题,当问题规模过大时,算法很难在短时间内找到最优的路径,所以在求解大规模航迹规划时,通常采用搜索算法,启发式算法,神经网络等方法 [4] [5] [6] 求得“尽可能好”的解。文献 [7] 以最小化无人机巡视环路的飞行时间为路径规划的目标,讨论了具有遥感能力的无人机的路径规划。文献 [8] 将巡检的路径规划问题制定为扩展的旅行商问题(Travelling salesman problem, TSP)。文献 [9] 提出了几种在搜救系统中无人机离散路径规划的方法。文献 [10] 以最小化航线数量为目标提出了新的算法。文献 [11] 提出了基于旋转卡尺的算法解决起降点在路径中的覆盖路径规划问题。文献 [12] 分析了现有文献中关于覆盖路径规划问题,还介绍了通常用于评估覆盖任务成功与否的性能指标。
目前关于二,三维航迹规划研究大多是研究点到点的航迹规划,即研究出发点与目标点间的最优航迹。而当灾害发生时,使用无人机对目标道路网络进行遥感成像,收集道路受损信息要求的是以航迹上的每个点为圆心,拍摄半径为半径的圆的集合覆盖整个道路区域。
2. 问题描述及模型分析
自然灾害发生后,道路被损坏的同时,通信线路也常常被破坏。我们需要n架续航能力为LE(Long-Endurance)的无人机完成对某一道路网络的监控覆盖,为救援人员提供道路状况信息。记道路网
络为区域
,无人机的理想可视半径为r0 (不考虑遮挡的情况下,实际可
视范围与无人机和地面的相对高度差
成正比,这里设在理想相对高度差
的可视半径为理想可视半径r0),每t秒拍摄一张照片(视频也可看作是照片的集合,例如帧率为60 fps的视频可是视作每1/60秒拍摄一张照片),航迹为由方程
确定的空间曲线Fi。则第i架无人机沿着曲线Fi的
监控区域
可表示为
.
则总的监控区域
,
综上所述,本文所研究的航迹规划问题可描述为:
从所有满足
的F曲线中找到n条最优的空间曲线Fi。
目前一般的三维航迹规划问题解决的是点到点的路径优化,其搜索空间为
,
其中
分别表示待监控道路区域的经度,纬度和海拔的最小值和最大值, 是
满足拍摄分辨率的最大相对高度。
然而由
知对于
,
有
,
即
所以,本文提出的初始可行域G比一般的搜索空间Ggeneral要小。
当地形起伏不大时,无人机的相对高度变化不大,可视半径的变化可忽略不计,可视半径可看作常数,本文将三维航迹规划问题转化为二维航迹规划问题以简化计算量,从而快速获得一个较好的解。
3. 近似算法研究
为解决基于无人机的道路监控的二维航迹规划问题,我们采用栅格法进行空间的划分。将空间离散化为mgrid行ngrid列的矩形单元后,道路监控的二维航迹规划问题即变为寻找一系列相邻的单元,使得以这些单元的中心为圆心,可视半径为圆的半径的实心圆覆盖道路所在的单元。传统的算法解决点到点之间的航迹规划使用的方法是采用网格划分规划空间后,构造以当前节点为中心的九宫图作为算法的数据结构,共有8个相邻节点,航迹由相邻网格节点连接形成,并连接起始节点和目标节点。此算法的时间
复杂度大于
。可见此方法的计算量十分巨大。因此我们提出以下三阶段算法:
第一阶段:首先求出
。当问题规模较小时,以道路所在单元中心为圆心画圆确定G,将 G1 = G 所在的单元视为位置集合覆盖问题(Location Set Covering Problem, LSCP)中的候选拍摄点,建立整数规划模型。使用分支定界法求解得到目标函数的最小值。当问题规模较大时,使用二分法确定覆盖整
个道路网络所需实心圆的最少数目
,的初始搜索区间
,其中
,
,
是道路面积,
是道路长度,r是无人机的可视半径。
第二阶段:记“每一个含有
个单元的集合”为一个方案,若由
个单元中心为圆心的实心圆集合能覆盖整个道路网络。显然当无人机的航迹经过方案中
个单元的中心时,无人机能完成道路监控任务。定义不同单元间中心的距离为两个单元间的边的权重,记录下此方案的总权重。比较不同方案的总权重,使用启发式算法寻找总权重较小的方案或者遍历得到总权重最小的方案;
第三阶段最后将总权重较小的方案中的
个单元中心坐标看作旅行商问题的城市坐标,求解该旅行商问题。三阶段算法的时间复杂度主要在第一阶段求解整数规划模型,确定s种有个候选拍摄点的方案,
算法步骤:
STEP1:对区域Ggeneral使用栅格法进行空间的划分,以道路所在单元中心为圆心画圆确定G,将
所在的单元视为LSCP问题中的候选拍摄点,建立整数规划模型。使用分支定界法求解,并确定 circle_min。
STEP2:检查STEP1中求出的
个拍摄点,若存在两点之间的距离小于
,(
为无人机的速度,
为无人机拍摄两张照片所需的最小飞行间距,在一般情况下,如果
两点间的距离小于
,可以降低无人机的速度,使得无人机从
飞到时,相机能够拍照,但是降低速度无疑会使完成任务的时间变长,而灾害救援系统有时间限制),则从
删去处于更“稠密”区域的那点。返回到STEP1重新计算。如果所有点的距离均不小于
,则计算当前方案的总权重,执行STEP3。
STEP3:使用启发式算法求解出满足预定精度的方案或使用遍历算法得到最优方案。
STEP4:求解出访问STEP3中求出的最优方案中
个单元的先后顺序。
4. 算例研究
如图1所示,本文将空间划分为16 × 16个正方形单元,单元的边长记为一个单位长度,黄色单元为道路。无人机可视半径设置为3个单位长度,即
。第一阶段求解得到
。第二,三阶段使用遍历算法求解,结果如图2,黄色与白色区域为初始可行域G。红色空心圆所在的第87号,第105号,第151号和第155号单元是总权重最小的方案的四个单元。由151号单元到87号单元再到105号单元最后到155号单元即为最终航迹。值得注意的是,本文并未预设无人机的起飞点,但这并不影响算法的有效性,为了符合生活实际,只需要在第三阶段的算法中添加一个起飞点的坐标,即可得到实际问题的解决方案。
5. 结束语
时间是灾害应急救援系统中的重要因素,受灾地区的救援物资运送受到道路的影响,使用无人机遥感成像可以快速提供道路受损状况,并弥补卫星遥感成像精度不足的问题。本文将灾害应急救援系统中无人机道路监控航迹规划问题解耦为LSCP问题和TSP问题,给出了覆盖整个道路网络所需实心圆的数
目的大致范围
与一种能够减少模型求解时间的三阶段算法。
基金项目
长沙理工大学大学生创新训练计划立项项目(无人机救灾响应系统开发)。
NOTES
*第一作者。