CN114812564B - 基于城市动态时空风险分析的多目标无人机路径规划方法 - Google Patents
基于城市动态时空风险分析的多目标无人机路径规划方法 Download PDFInfo
- Publication number
- CN114812564B CN114812564B CN202210707989.4A CN202210707989A CN114812564B CN 114812564 B CN114812564 B CN 114812564B CN 202210707989 A CN202210707989 A CN 202210707989A CN 114812564 B CN114812564 B CN 114812564B
- Authority
- CN
- China
- Prior art keywords
- risk
- flight
- unmanned aerial
- aerial vehicle
- coefficient
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000012502 risk assessment Methods 0.000 title claims abstract description 27
- 230000003068 static effect Effects 0.000 claims abstract description 62
- 230000003116 impacting effect Effects 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 4
- 230000000873 masking effect Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明的实施例提供了一种基于城市动态时空风险分析的多目标无人机路径规划方法,涉及无人机技术领域。通过构建目标区域的密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,利用模糊动态贝叶斯网络将三个模型融合后得到无人机飞行风险地图,再设计基于风险动力学约束的混合A*算法规划飞行航线飞行风险地图。这样同时评估了无人机在面对静态建筑物、以及随不同时段的城市事件影响变化、人流信息的动态风险因素得到飞行风险地图,基于风险约束的动力学路径规划可以得到保证飞行安全、路径平滑、代价更小的飞行航线。
Description
技术领域
本发明涉及无人机技术领域,具体而言,涉及一种基于城市动态时空风险分析的多目标无人机路径规划方法。
背景技术
随着城市人口的持续增长,以及城市交通的拥堵状况的加剧,当前地面交通系统已经达到了其使用极限。由于无人机具备多功能作业、且作业效率较高、作业成本低廉等特点,使得无人机在智慧城市体系的构建过程中可以发挥极大的作用。当前人们熟知的无人机应用领域包括了:商品运输、智能监控、交通调查、空中巴士和辅助通信等等。
当前无人机的城市飞行环境具备规划范围广、空中环境复杂以及非结构化等特点。面对这样复杂的城市飞行环境,无人机路径规划方法工作主要目标就是根据城市多种风险作出响应,规划满足安全飞行、节约人力物力的多目标路径。而现有技术的路径规划,存在以下三个方面的缺陷:
首先,现有技术的城市飞行地图主要基于地面建筑物、禁飞区等静态障碍进行风险分析得到的,然而对于城市人流、车流密集区等动态的因素没有考虑其对无人机飞行的影响风险。相应地,规划得到的飞行路径由于没有考虑动态的车流、人流因素,使得无人机按照飞行路径执行任务时具有一定的风险性。
其次,现有技术中基于传输飞行地图进行无人机路径规划以避开高风险区域的方法基本不随时间变化而变化。但是,无人机城市飞行环境是动态多变的,由于人群流动、聚集、城市功能在时间维度的改变,不同时间段下不同区域的风险系数也是变化。例如:露天休闲场所、学校等区域在节假日、工作日的人群密集程度是不同。基于该城市飞行地图规划的飞行路径无法保证不同时间段内无人机的飞行安全和飞行效率。
再则,无人机自身的动力局限性、任务的时效性以及复杂的城市飞行要求限制三者之间存在矛盾。现有技术路径规划时,是在二维平面中单纯的使用加权的方法来对每个区域进行风险评估,容易导致确定航线和高度方面的代价函数不严格,使得最终的飞行路径容易出现不必要的“之”字型路线和不平滑路线。
因此,基于城市动态时空风险评估的、双向管控飞行安全和飞行效率的多目标无人机路径规划方法提出是必要且有重要意义的。
发明内容
本发明的目的在于提供一种基于城市动态时空风险分析的多目标无人机路径规划方法,以改善现有技术存在的问题。
本发明的实施例可以这样实现:
本发明提供一种基于城市动态时空风险分析的多目标无人机路径规划方法,包括:
建立目标区域所在空间的三维网格模型;所述三维网格模型包含多个网格单元;
构建所述目标区域的密集人流区域风险模型,所述密集人流区域风险模型表征无人机在不同时段飞行时的人群密度对应的坠落风险系数;
构建所述目标区域的社会属性风险模型,所述社会属性风险模型表征所述无人机在所述不同时段飞行时的地区社会属性对应的碰撞风险系数;
构建所述目标区域的静态物理障碍风险模型,所述静态物理障碍风险模型表征所述无人机在所述不同时段飞行时的静态建筑物对应的飞行风险系数;
利用模糊动态贝叶斯网络融合所述密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,获得飞行风险地图;所述飞行风险地图为所述三维网格模型中包含所述不同时段下每个所述网格单元对应的总风险系数;
根据所述飞行风险地图,设计基于风险动力学约束的混合A*算法规划所述无人机在所述目标区域的飞行航线。
在可选的实施方式中,所述构建所述目标区域的密集人流区域风险模型的步骤,包括:
获取所述目标区域的人流量信息,并基于所述人流量信息得到人群密度系数;
基于所述无人机撞击地面的动能,确定所述无人机坠落导致的撞击风险率;
根据预设系数、所述人群密度系数、所述撞击风险率,确定所述坠落风险系数。
在可选的实施方式中,所述人群密度系数的表达式为:
所述撞击风险率的表达式为:
其中,为所述无人机撞击地面的动能;为遮蔽参数,取值区间为(0,1];为所述遮蔽参数为0.5时所述撞击风险率达到50%时所需的撞击能量;为所述遮蔽参数降到0时导致撞击事故所需的撞击能量值;为所述撞击风险率;
所述坠落风险系数的表达式为:
在可选的实施方式中,所述构建所述目标区域的社会属性风险模型的步骤,包括:
计算所述目标区域中每个兴趣点的访问概率;所述目标区域中包括多个兴趣点;
获取所述目标区域的社会属性信息;所述社会属性信息包括M个地区社会属性;
基于所述访问概率和所述社会属性信息,确认所述目标区域中第m种地区社会属性的目标访问概率;
获取所述目标区域的禁飞区域和非禁飞区域各自对应的禁飞风险系数;
基于所述目标访问概率和所述禁飞风险系数,确定所述碰撞风险系数。
在可选的实施方式中,所述访问概率的表达式为:
所述目标访问概率的表达式为:
所述碰撞风险系数的表达式为:
在可选的实施方式中,所述构建所述目标区域的静态物理障碍风险模型的步骤,包括:
获取与所述三维网格模型对应的静态建筑物信息;所述静态建筑物信息表征每个所述网格单元所在位置是否存在所述静态建筑物;
基于所述静态建筑物信息,确定所述飞行风险系数。
在可选的实施方式中,所述利用模糊动态贝叶斯网络融合所述密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,获得飞行风险地图的步骤,包括:
在每个所述网格单元,根据该网格单元对应的坠落风险系数、碰撞风险系数、飞行风险系数的大小,确定所述不同时段中该网格单元对应的总风险系数。
在可选的实施方式中,所述飞行风险地图的表达式为:
其中,表示位置点所在的网格单元在时段对应的总风险系数;表示所述位置点所在的网格单元在时段对应的坠落风险系数;表示所述位置点所在的网格单元在时段对应的碰撞风险系数;表示所述位置点所在的网格单元在时段对应的飞行风险系数;
在可选的实施方式中,所述根据所述飞行风险地图,设计基于风险动力学约束的混合A*算法规划所述无人机在所述目标区域的飞行航线的步骤,包括:
获取起点和终点的位置坐标;
对所述飞行风险地图进行预处理,得到稀疏风险地图;所述稀疏风险地图中的所有网格单元均为所述无人机可飞行的区域;
基于代价函数,设计所述基于风险动力学约束的混合A*算法搜索所述起点至所述终点之间的飞行航线。
与现有技术相比,本发明实施例提供了一种基于城市动态时空风险分析的多目标无人机路径规划方法,通过构建的目标区域的密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,这三个模型从四维(三维空间和时间维度)角度出发分别考虑了在不同时段下,城市的人群密度对无人机飞行的影响、城市中不同区域各自持有的地区社会属性对无人机飞行的影响、城市中的静态建筑物对无人机飞行的影响。再利用模糊动态贝叶斯网络将三个模型融合后得到飞行风险地图,再设计基于风险动力学约束的混合A*算法规划飞行航线。同时评估了无人机在面对静态建筑物、以及随不同时段、城市事件影响变化、人流信息的动态风险因素,得到的飞行风险地图更加全面,基于风险约束的动力学路径规划可以得到保证飞行安全、路径平滑、代价更小的飞行航线。使得无人机在保证在满足复杂城市环境中飞行需求的前提下,高效地完成各种飞行任务。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1为本发明实施例提供的一种基于城市动态时空风险分析的多目标无人机路径规划方法的流程示意图之一。
图2为本发明实施例提供的一种基于城市动态时空风险分析的多目标无人机路径规划方法的流程示意图之二。
图3为本发明实施例提供的一种基于城市动态时空风险分析的多目标无人机路径规划方法的流程示意图之三。
图4为本发明实施例提供的一种基于城市动态时空风险分析的多目标无人机路径规划方法的流程示意图之四。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
需要说明的是,在不冲突的情况下,本发明的实施例中的特征可以相互结合。
随着城市人口的持续增长,以及城市交通的拥堵状况的加剧,当前地面交通系统已经达到了其使用极限。由于无人机具备多功能作业、且作业效率较高、作业成本低廉等特点,使得无人机在智慧城市体系的构建过程中可以发挥极大的作用。当前人们熟知的无人机应用领域包括了:商品运输、智能监控、交通调查、空中巴士和辅助通信等等。
在城市飞行环境中,由于其规划范围广、空中环境复杂以及非结构化的特点,规划得到无人机在城市飞行的路径的前提是需要预先得到城市飞行地图,再基于城市飞行地图进行路径规划得到无人机的飞行路径。
由于城市人口密度大,地面状况复杂,因此无人机复杂城市环境中飞行至少需要保证达到以下三种要求:(1)安全性,需要考虑地面动态流动的人群、车流以及静态建筑等影响,保证在安全的路径执行飞行任务;(2)高效性,高效使用有限的城市空域;(3)隐私性,避开需要保护隐私的禁飞区域。因此,在满足城市飞行的多种需求下,确保在复杂城市环境中无人机可顺利完成任何类型的常规任务的路径规划算法十分重要。
在现有技术中,大部分无人机路径规划时利用的城市飞行地图,主要是针对地面建筑物、禁飞区等静态障碍进行风险评估得到的,但是这样考虑的风险因素不够全面,无法得到较为全面的城市飞行地图。
而现有技术的路径规划,存在以下三方面的缺陷:
首先,现有技术的城市飞行地图主要基于地面建筑物、禁飞区等静态障碍进行风险分析得到的,然而对于城市人流、车流密集区等动态的因素没有考虑其对无人机飞行的影响风险。相应地,规划得到的飞行路径由于没有考虑动态的车流、人流因素,使得无人机按照飞行路径执行任务时具有一定的风险性。
其次,现有技术中基于传输飞行地图进行无人机路径规划以避开高风险区域的方法基本不随时间变化而变化。但是,无人机城市飞行环境是动态多变的,由于人群流动、聚集、城市功能在时间维度的改变,不同时间段下不同区域的风险系数也是变化。例如:露天休闲场所、学校等区域在节假日、工作日的人群密集程度是不同。基于该城市飞行地图规划的飞行路径无法保证不同时间段内无人机的飞行安全和飞行效率。
再则,无人机自身的动力局限性、任务的时效性以及复杂的城市飞行要求限制三者之间存在矛盾。现有技术路径规划时,是在二维平面中单纯的使用加权的方法来对每个区域进行风险评估,容易导致确定航线和高度方面的代价函数不严格,使得最终的飞行路径容易出现不必要的“之”字型路线和不平滑路线。
基于上述技术问题的发现,发明人经过创造性劳动提出下述技术方案以解决或者改善上述问题。需要注意的是,以上现有技术中的方案所存在的缺陷,均是发明人在经过实践并仔细研究后得出的结果,因此,上述问题的发现过程以及下文中本申请实施例针对上述问题所提出的解决方案,都应该是发明人在发明创造过程中对本申请做出的贡献,而不应当理解为本领域技术人员所公知的技术内容。
因此,本发明实施例提供一种基于城市动态时空风险分析的多目标无人机路径规划方法,通过构建目标区域的密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,从四维(三维空间和时间维度)角度出发分别考虑了在不同时段下,城市的人群密度对无人机飞行的影响、城市中不同区域各自持有的地区社会属性对无人机飞行的影响、城市中的静态建筑物对无人机飞行的影响,再利用模糊动态贝叶斯网络融合三个风险模型得到飞行风险地图,最后设计基于风险动力学约束的利用混合A*算法进行路径规划得到无人机的飞行航线。以下通过实施例,并配合所附附图,进行详细说明。
请参考图1,图1为本发明实施例提供的一种基于城市动态时空风险分析的多目标无人机路径规划方法的流程示意图,该方法的执行主体可以是一种电子设备,该方法包括以下步骤:
S110、建立目标区域所在空间的三维网格模型。
S120、构建目标区域的密集人流区域风险模型。
在本实施例中,密集人流区域风险模型表征无人机在不同时段飞行时的人群密度对应的坠落风险系数。结合交通流量信息,密集人流区域风险模型可以反映密集人流区域对无人机在不同时段飞行时的风险影响程度。
S130、构建目标区域的社会属性风险模型。
在本实施例中,社会属性风险模型表征无人机在不同时段飞行时的地区社会属性对应的碰撞风险系数。结合城市兴趣点(Point of Interest,POI),社会属性风险模型可以反映地区社会属性对无人机在不同时段飞行时的风险影响程度。
S140、构建目标区域的静态物理障碍风险模型。
本实施例中,静态物理障碍风险模型表征无人机在不同时段飞行时的静态建筑物对应的飞行风险系数。利用静态建筑物信息,静态物理障碍风险模型可以反映无人机飞行时受物理障碍的影响程度。
S150、利用模糊动态贝叶斯网络融合密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,获得飞行风险地图。
在本实施例中,可以设计基于模糊动态贝叶斯网络,针对每个网格单元,将三种风险系数汇总得到对应的总风险系数,绘制四维的飞行风险地图,不同时段下每个网格单元的风险成本将被计算出来并呈现在地图上,用于后续的无人机路径规划。
可以理解,飞行风险地图为三维网格模型中包含不同时段下每个网格单元对应的总风险系数。密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型三者均为从时间和空间维度出发,分别代表了三种影响因素(动态的人流量、不同的地区社会属性、静态建筑物)对无人机飞行的影响程度。而飞行风险地图中综合考量了三种影响因素的影响程度,融合得到每个网格单元对应的总风险系数。
S160、根据飞行风险地图,设计基于风险动力学约束的混合A*算法规划无人机在目标区域的飞行航线。
在飞行风险地图的基础上,可以设计基于风险动力学约束的混合A*算法(Risk-based Hybrid A*)来规划无人机的飞行航线。传统A*算法中将无人机看成一个质点,那么无人机三维笛卡尔坐标网格(x,y,z)就可以看作是它的状态向量。然而,无人机由于飞行性能问题,其在不同网格间转弯或抬高降低高度收到自身惯性影响,只选择网格中心点规划的路径不够平滑,且不能考虑惯性对无人机飞行成本代价的消耗和安全影响。
而本实施例采用风险动力学约束的混合A*算法,可以将无人机视为刚体,将其大小、朝向等纳入飞行状态考虑。
需要说明的是,在实际情况中,上述步骤S120、S130、S140之间的执行顺序不以图1所示为限,S120、S130、S140可以并列执行,也可以任意顺序先后执行。
本发明实施例提供了一种基于城市动态时空风险分析的多目标无人机路径规划方法,首先建立目标区域所在空间的包含多个网格单元三维网格模型。然后基于城市的动态时空风险分析分别构建了目标区域中的密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,再利用模糊动态贝叶斯网络将三个模型融合后得到飞行风险地图,再利用混合A*算法规划飞行航线。这样,同时基于建筑物、禁飞区这两个静态影响因素以及不同时段中变化的人流信息这一动态风险因素这三种不同的角度对无人机飞行的风险性进行分析评估,得到的飞行风险地图更加全面,路径规划得到飞行航线更能保障无人机的飞行安全,使得无人机在保证在满足复杂城市环境中飞行需求的前提下,高效地完成各种飞行任务。
三维空间模型是空间维度,在时间维度,可以将每天划分为T个时间段。以下通过介绍在T个时间段中的任意一个时间段t,无人机的三种风险模型各自对应的风险系数。
在可选的实施方式中,空中飞行无人机可能会失去控制或动力而坠落,无人机坠落在人流量密集区域,容易发生事故,例如,无人机坠落导致行人受伤或者致死。密集人流区域风险模型可以基于目标区域的地面人流、车流密度来评估无人机飞行的坠落风险系数,相应地,请参见图2,上述步骤S120的子步骤可以包括:
S121、获取目标区域的人流量信息,并基于人流量信息得到人群密度系数。
在本实施例中,人流量信息可以通过手机信令或者交通客流数据估计得到。对人流量信息利用sigmoid函数进行归一化处理可以得到人群密度系数。
可以理解,归一化处理之后,人群密度系数的取值区间在[0,1]。从三维网格模型的二维层面出发,该人群密度系数的表达式可以如下所示:
S122、基于无人机撞击地面的动能,确定无人机坠落导致的撞击风险率。
无人机坠落的动能取决于飞行高度、无人机质量和体积大小等,无人机撞击地面的动能可以表示为:
可以理解,撞击风险率可以表示无人机坠落地面导致事故的概率。撞击风险率可以表示为:
S123、根据预设系数、人群密度系数、撞击风险率,确定坠落风险系数。
在可选的实施方式中,由于城市地区中不同区域具有不同的社会属性,其风险代价和飞行规则都不相同。例如,在空间维度,城市的许多区域需要隐私保护,属于禁飞区域;在时间层面,学校、公园等地区的人流量在工作日、和休息日以及一天中的不同时段人流量存在差异,相应地,不同时段,学校、公园等区域无人机飞行的风险系数也有所不同。
因此,社会属性风险模型可以是考虑了不同的地区社会属性在不同时段对无人机飞行的影响。相应地,请参见图3,上述步骤S130的子步骤可以包括:
S131、计算目标区域中每个兴趣点的访问概率。
在本实施例中,目标区域中包括多个兴趣点。可以先获取目标区域POI数据,POI数据可以包括目标区域中的各个兴趣点的位置信息等。兴趣点的类型可以是公园、学校、社区等等。然后再计算每个兴趣点在时段的访问概率,该访问概率的表达式可以是:
因此,一个兴趣点的访问概率为周围50米范围内所有访问点的访问概率之和。
S132、获取目标区域的社会属性信息。
在本实施例中,社会属性信息包括M个地区社会属性。目标区域中,不同区域可以对应不同的地区社会属性。例如,有的地区社会属性所在区域可能属于禁飞区、禁止鸣笛区域或者限速区等等。
S133、基于访问概率和社会属性信息,确认目标区域中第m种地区社会属性的目标访问概率。
在本实施例中,每种地区社会属性所在区域的目标访问概率为该区域内的各个兴趣点的访问概率之和。第m种地区社会属性为M个地区社会属性中的任意一个。第m种地区社会属性所在区域对应的目标访问概率的表达式可以为:
S134、获取目标区域的禁飞区域和非禁飞区域各自对应的禁飞风险系数。
在本实施例中,禁飞风险系数可以是从城市空域规定信息中得到的。
S135、基于目标访问概率和禁飞风险系数,确定碰撞风险系数。
在可选的实施方式中,静态物理障碍风险模型取决于目标区域的静态建筑物,考虑了静态建筑物对无人机的飞行影响。相应地,上述步骤S140的子步骤可以包括:
S141、获取与三维网格模型对应的静态建筑物信息。
静态建筑物信息可以表征每个网格单元所在位置是否存在静态建筑物。
S142、基于静态建筑物信息,确定飞行风险系数。
飞行风险系数的表达式可以是:
在可选的实施方式中,基于上述介绍的三种风险模型各自对应的风险系数,再进一步融合可以得到飞行风险地图。在三种风险模型融合的过程中,可以利用模糊动态贝叶斯网络(MFDBN)来解决由于数据不充分和传播不准确原因造成的风险评估不确定问题,以及人群密度与时间和事件关系密切导致的风险评估需要随城市人口流动趋势进行实时调整的问题。
示例性地,可以将模糊动态贝叶斯网络与动态风险评估量化方法(DRA)结合。首先在数据不确定的情况下,同时进行不确定性随时间的传播量化和动态风险评估量化(DRA),利用动态贝叶斯网络(DBN)使用三角模糊数,以完全保留t时刻不确定性信息。
示例性地,上述步骤S150的子步骤可以包括:
S151、在每个网格单元,根据该网格单元对应的坠落风险系数、碰撞风险系数、飞行风险系数的大小,确定不同时段中该网格单元对应的总风险系数。
在本实施例中,飞行风险地图的表达式可以为:
其中,表示位置点所在的网格单元在时段对应的总风险系数;表示位置点所在的网格单元在时段对应的坠落风险系数;表示位置点所在的网格单元在时段对应的碰撞风险系数;表示位置点所在的网格单元在时段对应的飞行风险系数。
表示当坠落风险系数、碰撞风险系数、飞行风险系数均小于等于零时,对应的总风险系数为0;表示当坠落风险系数、碰撞风险系数、飞行风险系数均处于(0,1)区间内时,对应的总风险系数为;表示当存在坠落风险系数、碰撞风险系数、飞行风险系数中的任意一个大于等于1时,对应的总风险系数为1。
示例性地,请参见图4,上述步骤S160飞行风险地图可以包括以下子步骤:
S161、获取起点和终点的位置坐标。
S162、对飞行风险地图进行预处理,得到稀疏风险地图。
在本实施例中,起点和终点均位于目标区域之内。稀疏风险地图中的所有网格单元均为无人机可飞行的区域。可以利用无人机自身性能约束,剔除无效网格,得到稀疏风险地图。示例性地,无效网格可以是静态建筑物部分以及静态建筑物往上预设高度部分对应的网格单元,预设高度部分可以为1米或者2米。这样可以提高后续基于风险动力学约束的混合A*算法基于稀疏风险地图搜索飞行航线的速度。
S163、基于代价函数,设计基于风险动力学约束的混合A*算法搜索起点至终点之间的飞行航线。
在本实施例中,飞行航线可以由起点到终点之间的S个路径点连接而成。S个路径点组成路径点集合。
设计风险动力学约束的混合A*算法进行路径搜索的过程可以是:从起点开始,先计算起点所在网格单元周围的每个相邻网格单元的代价值,并在代价值最小的相邻网格单元中确定出起点后的第一个路径点;然后计算第一个路径点所在网格单元周围的每个相邻网格单元的代价值,并在代价值最小的相邻网格单元确定出第一个路径点之后的第二个路径点……以此类推,直到确定出终点前的第S个路径点。
下面给出一种根据飞行风险地图,采用混合A*算法规划无人机在目标区域的飞行航线的伪代码示例:
Begin:
2.Delete 网格单元 //预处理,得到稀疏风险地图
If <扩展过程遇到障碍物>
End
return 子网格列表
For n In子网格列表
If <n Not In ( OPEN Or CLOSED)>
Then 加入OPEN表
End
Else if <子网格 In OPEN>
Then 更新代价值
End
当前网格单元=最小代价值网格单元 In OPEN
If <当前网格单元=目标网格单元>
then返回路径
End
Else
当前网格单元 to CLOSED
End
End
Then return 路径
End
End
假设第s个路径点(可以是飞行航线上的任意一个路径点)所在的网格单元称为当前网格单元。在确定了第s个路径点后,则会从当前网格单元周围的相邻网格单元中选出代价值最小的网格单元,并在该代价值最小的网格单元中确定出第s+1个路径点。最终获得从起点到终点的飞行航线规划。
在利用代价函数计算当前网格单元的每个相邻网格单元的代价值时,对于某个相邻网格单元A,利用代价函数计算的代价值与两个方面相关:一方面是该相邻网格单元A与终点之间的直线距离,相较于当前网格单元的其余几个相邻网格单元,该直线距离越大,代价值会相对越大,反之则相对越小;另一方面是该相邻网格单元与静态建筑物的距离,相较于当前网格单元的其余几个相邻网格单元,该相邻网格单元A离静态建筑物距离越近,相应的动力学衰减系数越大,对应的代价值会相对越大,反之则相对越小。
在本实施例中,在三维空间网格中,新增了一个维度θ。假设第s个路径点到终点为一条直线,θ表示该直线与x轴方向的夹角。
在当前路径点处选择下一个路径点时(该当前路径点为路径点集合中的任意一个路径点),对当前状态进行描述,表示为:
其中,为预设的扩展步长参数,下标为代表离散的状态向量,改变和得到可作为下一个路径点的扩展子列表。扩充列表过程中,若与障碍物相撞说明网格列表与障碍物相距较近,根据子列表网格与障碍物相撞次数,加权增加子网格代价值。
需要说明的是,上述方法实施例中各个步骤的执行顺序不以附图所示为限制,各步骤的执行顺序以实际应用情况为准。
与现有技术相比,本发明实施例具有以下有益效果:
通过构建的目标区域的密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,从四维(三维空间和时间维度)角度出发分别考虑了在不同时段下,城市的人群密度对无人机飞行的影响、城市中不同区域各自持有的地区社会属性对无人机飞行的影响、城市中的静态建筑物对无人机飞行的影响。最终利用模糊动态贝叶斯网络融合这三个方面的影响,得到了目标区域的飞行风险地图,该飞行风险地图的表现形式为三维网格模型中的每个网格单元在不同时段对应的不同的总风险系数。这样,不仅考虑到了建筑物和空中禁飞区等静态影响因素的,还考虑了不同时段中变化的人流信息这个动态风险因素,得到的飞行风险地图更加全面,能够在保障无人机在城市复杂多变的地面环境下的安全飞行,动态的对无人机飞行代价进行综合智能的评估计算,减小对地面人群、车辆、建筑物等的碰撞风险。
设计基于风险动力学约束的混合A*算法,基于该飞行风险地图从四维层面进行路径规划得到无人机的飞行航线,不同于传统网格中心的飞行轨迹点选择方式,本方案针对无人机飞行性能、速度和惯性,利用基于风险动力学约束的混合A*算法,优化其路径平滑性,使路径长度、消耗的能量和路径风险同时达到最小,合并多目标获取代价小的飞行航线。这样基于风险约束的动力学路径规划可以得到保证飞行安全、路径平滑、代价更小的飞行航线,使得无人机在保证在满足复杂城市环境中飞行需求的前提下,以高效能完成各种飞行任务。
综上所述,本发明实施例提供了一种基于城市动态时空风险分析的多目标无人机路径规划方法,通过构建的目标区域的密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,从四维(三维空间和时间维度)角度出发分别考虑了在不同时段下,城市的人群密度对无人机飞行的影响、城市中不同区域各自持有的地区社会属性对无人机飞行的影响、城市中的静态建筑物对无人机飞行的影响,利用模糊动态贝叶斯网络将三个模型融合后得到飞行风险地图,再设计基于风险动力学约束的混合A*算法来规划无人机的飞行航线。这样,同时评估了无人机在面对静态建筑物、以及随不同时段、城市事件影响变化、人流信息的动态风险因素,得到的飞行风险地图更加全面,基于风险约束的动力学路径规划可以得到保证飞行安全、路径平滑、代价更小的飞行航线。
以上,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (7)
1.一种基于城市动态时空风险分析的多目标无人机路径规划方法,其特征在于,包括:
建立目标区域所在空间的三维网格模型;所述三维网格模型包含多个网格单元;
构建所述目标区域的密集人流区域风险模型,所述密集人流区域风险模型表征无人机在不同时段飞行时的人群密度对应的坠落风险系数;
构建所述目标区域的社会属性风险模型,所述社会属性风险模型表征所述无人机在所述不同时段飞行时的地区社会属性对应的碰撞风险系数;
构建所述目标区域的静态物理障碍风险模型,所述静态物理障碍风险模型表征所述无人机在所述不同时段飞行时的静态建筑物对应的飞行风险系数;
利用模糊动态贝叶斯网络融合所述密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,获得飞行风险地图;所述飞行风险地图为所述三维网格模型中包含所述不同时段下每个所述网格单元对应的总风险系数;
根据所述飞行风险地图,设计基于风险动力学约束的混合A*算法规划所述无人机在所述目标区域的飞行航线;
所述构建所述目标区域的社会属性风险模型的步骤,包括:
计算所述目标区域中每个兴趣点的访问概率;所述目标区域中包括多个兴趣点;
获取所述目标区域的社会属性信息;所述社会属性信息包括M个地区社会属性;
基于所述访问概率和所述社会属性信息,确认所述目标区域中第m种地区社会属性的目标访问概率;
获取所述目标区域的禁飞区域和非禁飞区域各自对应的禁飞风险系数;
基于所述目标访问概率和所述禁飞风险系数,确定所述碰撞风险系数;
所述访问概率的表达式为:
所述目标访问概率的表达式为:
所述碰撞风险系数的表达式为:
2.根据权利要求1所述的方法,其特征在于,所述构建所述目标区域的密集人流区域风险模型的步骤,包括:
获取所述目标区域的人流量信息,并基于所述人流量信息得到人群密度系数;
基于所述无人机撞击地面的动能,确定所述无人机坠落导致的撞击风险率;
根据预设系数、所述人群密度系数、所述撞击风险率,确定所述坠落风险系数。
4.根据权利要求1所述的方法,其特征在于,所述构建所述目标区域的静态物理障碍风险模型的步骤,包括:
获取与所述三维网格模型对应的静态建筑物信息;所述静态建筑物信息表征每个所述网格单元所在位置是否存在所述静态建筑物;
基于所述静态建筑物信息,确定所述飞行风险系数。
5.根据权利要求1所述的方法,其特征在于,所述利用模糊动态贝叶斯网络融合所述密集人流区域风险模型、社会属性风险模型、静态物理障碍风险模型,获得飞行风险地图的步骤,包括:
在每个所述网格单元,根据该网格单元对应的坠落风险系数、碰撞风险系数、飞行风险系数的大小,确定所述不同时段中该网格单元对应的总风险系数。
7.根据权利要求1所述的方法,其特征在于,所述根据所述飞行风险地图,设计基于风险动力学约束的混合A*算法规划所述无人机在所述目标区域的飞行航线的步骤,包括:
获取起点和终点的位置坐标;
对所述飞行风险地图进行预处理,得到稀疏风险地图;所述稀疏风险地图中的所有网格单元均为所述无人机可飞行的区域;
基于代价函数,设计所述基于风险动力学约束的混合A*算法搜索所述起点至所述终点之间的飞行航线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210707989.4A CN114812564B (zh) | 2022-06-22 | 2022-06-22 | 基于城市动态时空风险分析的多目标无人机路径规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210707989.4A CN114812564B (zh) | 2022-06-22 | 2022-06-22 | 基于城市动态时空风险分析的多目标无人机路径规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114812564A CN114812564A (zh) | 2022-07-29 |
CN114812564B true CN114812564B (zh) | 2022-09-20 |
Family
ID=82520311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210707989.4A Active CN114812564B (zh) | 2022-06-22 | 2022-06-22 | 基于城市动态时空风险分析的多目标无人机路径规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114812564B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116227857B (zh) * | 2023-02-15 | 2023-11-21 | 上海正弘建设工程顾问有限公司 | 一种基于大数据的工程信息智能化分析系统及方法 |
CN116523384B (zh) * | 2023-04-11 | 2024-05-14 | 中国人民解放军陆军工程大学 | 基于独立否决和联合否决的无人机效能确定方法及系统 |
CN117993227B (zh) * | 2024-04-03 | 2024-06-18 | 北京工业大学 | 人流中的超越行为的仿真方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106970648A (zh) * | 2017-04-19 | 2017-07-21 | 北京航空航天大学 | 城市低空环境下无人机多目标路径规划联合搜索方法 |
CN108196575A (zh) * | 2018-01-05 | 2018-06-22 | 湖北工业大学 | 一种无人机任务分配及路线规划方法 |
CN109839110A (zh) * | 2019-01-09 | 2019-06-04 | 浙江大学 | 一种基于快速随机搜索树的多目标点路径规划方法 |
CN112215407A (zh) * | 2020-09-23 | 2021-01-12 | 中国民航大学 | 一种三维无人机安全航路规划方法 |
CN112880684A (zh) * | 2021-01-18 | 2021-06-01 | 中国人民解放军空军工程大学 | 一种城市空间无人机安全航路规划方法 |
CN113670309A (zh) * | 2021-07-21 | 2021-11-19 | 南京航空航天大学 | 考虑安全风险与噪声影响的城市低空无人机路径规划方法 |
CN114034301A (zh) * | 2021-10-21 | 2022-02-11 | 北京航空航天大学杭州创新研究院 | 一种基于决策树的实时航路规划方法 |
CN114579892A (zh) * | 2022-05-06 | 2022-06-03 | 南京航空航天大学 | 一种基于跨城市兴趣点匹配的用户异地访问位置预测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8538673B2 (en) * | 2008-10-31 | 2013-09-17 | Czech Technical University In Prague | System and method for planning/replanning collision free flight plans in real or accelerated time |
US10789853B2 (en) * | 2016-06-10 | 2020-09-29 | ETAK Systems, LLC | Drone collision avoidance via air traffic control over wireless networks |
US11573577B2 (en) * | 2019-01-30 | 2023-02-07 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Method and system for optimal trajectory path tasking for an unmanned aerial vehicle (UAV) |
CN113932813A (zh) * | 2021-09-30 | 2022-01-14 | 北航(四川)西部国际创新港科技有限公司 | 智能空地载具、其风险地图构建及航路规划方法与装置 |
-
2022
- 2022-06-22 CN CN202210707989.4A patent/CN114812564B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106970648A (zh) * | 2017-04-19 | 2017-07-21 | 北京航空航天大学 | 城市低空环境下无人机多目标路径规划联合搜索方法 |
CN108196575A (zh) * | 2018-01-05 | 2018-06-22 | 湖北工业大学 | 一种无人机任务分配及路线规划方法 |
CN109839110A (zh) * | 2019-01-09 | 2019-06-04 | 浙江大学 | 一种基于快速随机搜索树的多目标点路径规划方法 |
CN112215407A (zh) * | 2020-09-23 | 2021-01-12 | 中国民航大学 | 一种三维无人机安全航路规划方法 |
CN112880684A (zh) * | 2021-01-18 | 2021-06-01 | 中国人民解放军空军工程大学 | 一种城市空间无人机安全航路规划方法 |
CN113670309A (zh) * | 2021-07-21 | 2021-11-19 | 南京航空航天大学 | 考虑安全风险与噪声影响的城市低空无人机路径规划方法 |
CN114034301A (zh) * | 2021-10-21 | 2022-02-11 | 北京航空航天大学杭州创新研究院 | 一种基于决策树的实时航路规划方法 |
CN114579892A (zh) * | 2022-05-06 | 2022-06-03 | 南京航空航天大学 | 一种基于跨城市兴趣点匹配的用户异地访问位置预测方法 |
Non-Patent Citations (3)
Title |
---|
Yao Ji ; Chao Dong ; Xiaojun Zhu ; Qihui Wu.Fair-Energy Trajectory Planning for Multi-Target Positioning Based on Cooperative Unmanned Aerial Vehicles.《IEEE Access》.2019,全文. * |
基于配送模式的无人机城市配送路径规划;任新惠,武彤;《科学技术与工程》;20211120;第21卷(第32期);全文 * |
无人机路径规划中的环境和威胁模型研究;高晓静等;《航空计算技术》;20130531;第43卷(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114812564A (zh) | 2022-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114812564B (zh) | 基于城市动态时空风险分析的多目标无人机路径规划方法 | |
US10706729B2 (en) | Joint search method for UAV multiobjective path planning in urban low altitude environment | |
Wen et al. | A hybrid temporal association rules mining method for traffic congestion prediction | |
Pang et al. | UAV path optimization with an integrated cost assessment model considering third-party risks in metropolitan environments | |
EP3035314B1 (en) | A traffic data fusion system and the related method for providing a traffic state for a network of roads | |
Bai et al. | UAV path planning based on improved A∗ and DWA algorithms | |
US8082102B2 (en) | Computing flight plans for UAVs while routing around obstacles having spatial and temporal dimensions | |
Pan et al. | Monkey king evolution: an enhanced ebb-tide-fish algorithm for global optimization and its application in vehicle navigation under wireless sensor network environment | |
Ahmed et al. | Vehicular traffic noise prediction and propagation modelling using neural networks and geospatial information system | |
Wang et al. | Stop-and-wait: Discover aggregation effect based on private car trajectory data | |
CN109685286B (zh) | 未知静态障碍环境下usv基于改进蚁群优化的避碰规划方法 | |
Ntakolia et al. | A swarm intelligence graph-based pathfinding algorithm (SIGPA) for multi-objective route planning | |
Sarmady et al. | Modeling groups of pedestrians in least effort crowd movements using cellular automata | |
CN111767597A (zh) | 一种城市模型验证方法、装置、设备及存储介质 | |
Pang et al. | A risk-based UAS traffic network model for adaptive urban airspace management | |
Ma et al. | Adaptive path planning method for UAVs in complex environments | |
Salehinejad et al. | Dynamic Fuzzy Logic‐Ant Colony System‐Based Route Selection System | |
Yang et al. | A novel urban emergency path planning method based on vector grid map | |
CN115355922A (zh) | 一种基于改进蚁群算法的出行路径规划方法及系统 | |
Li et al. | A reinforcement learning-based routing algorithm for large street networks | |
Hildemann et al. | 3D-flight route optimization for air-taxis in urban areas with Evolutionary Algorithms and GIS | |
Cai et al. | A novel vector-based dynamic path planning method in urban road network | |
Wu et al. | A flood-discharge-based spatio-temporal diffusion method for multi-target traffic hotness construction from trajectory data | |
Chen et al. | Urban road network partitioning based on bi-modal traffic flows with multiobjective optimization | |
Zhang et al. | Design of the fruit fly optimization algorithm based path planner for UAV in 3D environments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |