CN115294487B - Unmanned aerial vehicle risk level obtaining method, storage medium and electronic equipment - Google Patents
Unmanned aerial vehicle risk level obtaining method, storage medium and electronic equipment Download PDFInfo
- Publication number
- CN115294487B CN115294487B CN202211225542.XA CN202211225542A CN115294487B CN 115294487 B CN115294487 B CN 115294487B CN 202211225542 A CN202211225542 A CN 202211225542A CN 115294487 B CN115294487 B CN 115294487B
- Authority
- CN
- China
- Prior art keywords
- target
- area
- sub
- unmanned aerial
- aerial vehicle
- 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 35
- 238000011282 treatment Methods 0.000 claims abstract description 3
- 239000013598 vector Substances 0.000 claims description 35
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims 1
- 238000011156 evaluation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000010006 flight Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/17—Terrestrial scenes taken from planes or by drones
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/22—Image preprocessing by selection of a specific region containing or referencing a pattern; Locating or processing of specific regions to guide the detection or recognition
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Economics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Multimedia (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention relates to the field of data acquisition, in particular to an unmanned aerial vehicle risk level acquisition method, a storage medium and electronic equipment, wherein the method comprises the following steps of; determining the central point coordinate of each sub-area in the target area D as a characteristic point coordinate to obtain a characteristic point coordinate set W = (W) 1 ,w 2 ,w 3 ,...,w q ) (ii) a Obtaining a differential model m 0 *L"(t)=m 0 * g-c | L ″ (t); to w j The following treatments were carried out: obtaining w according to the differential model, the horizontal velocity distribution data, the vertical velocity distribution data and the flight height distribution data j Corresponding set of coordinates of the point of falling COOR j (ii) a According to COOR j Determining P ris (x j ,y j ) (ii) a According to P ris (x j ,y j ) Determining P ris (x j ,y j ) A first target risk level for the corresponding sub-region. Thus, the accuracy of risk level assessment for the target area can be improved.
Description
Technical Field
The invention relates to the field of data acquisition, in particular to an unmanned aerial vehicle risk level acquisition method, a storage medium and electronic equipment.
Background
Unmanned aerial vehicle has nimble practical characteristics, is used for aerial photography, air transportation or city monitoring task etc. for example by wide application in military field and civilian field at present, in order to reduce the emergence that unmanned aerial vehicle falls, can carry out risk level aassessment to unmanned aerial vehicle's flight area territory.
At present, the risk level of a flight area of an unmanned aerial vehicle is evaluated, most of the risk evaluation is carried out when the unmanned aerial vehicle collides with the ground, and the risk evaluation research of the unmanned aerial vehicle for the ground collision only considers the ground information such as population density of the flight area so as to estimate the number of casualties caused by the ground collision of the unmanned aerial vehicle, and then the risk level of the flight area is evaluated according to the estimated number of casualties.
However, since there are fewer factors to be considered in the risk assessment of the ground collision of the unmanned aerial vehicle, the accuracy of the risk level assessment of the flight area of the unmanned aerial vehicle is low.
Disclosure of Invention
Aiming at the technical problem, the technical scheme adopted by the invention is as follows:
according to an aspect of the disclosure, there is provided an unmanned aerial vehicle risk level obtaining method, including:
determining the central point coordinate of each sub-area in the target area D as a characteristic point coordinate to obtain a characteristic point coordinate set W = (W) 1 ,w 2 ,w 3 ,...,w q ),w j =(x j ,y j ) (ii) a Wherein j =1,2, 3.. Q, q is the number of subregions, w j Is the jth feature point coordinate, x j Is w j Abscissa of (a), y j Is w j The ordinate of (c).
Obtaining a differential model m corresponding to the target unmanned aerial vehicle 0 *L"(t)=m 0 * g-c | L ″ (t); the unmanned aerial vehicle control method comprises the following steps of obtaining a target unmanned aerial vehicle, wherein t is the falling duration of the target unmanned aerial vehicle, L (t), L '(t) and L' (t) are all three-dimensional vectors, L (t) is a displacement vector of the target unmanned aerial vehicle corresponding to t, g is a gravity acceleration vector of the target unmanned aerial vehicle corresponding to t, L '(t) is a speed vector of the target unmanned aerial vehicle corresponding to t, and L' (t) is an acceleration vector of the target unmanned aerial vehicle corresponding to t; c is a resistance parameter c =0.5 × ρ AC of the target unmanned aerial vehicle D And ρ is an air density constant.
To w j The following treatments were carried out:
according to differential model, L = v ″ (0) = v) 0 And L (0) = (x) j ,y j ,h 0 ) Obtaining L (t) = (L) a1j ,L b1j ,L c1j ) (ii) a Wherein, L ″ (0) = v ″ 0 And L (0) = (x) j ,y j ,h 0 ) All the initial conditions are corresponding to the differential model; l' (0) is a velocity vector of the target drone at t =0, and L (0) is a displacement vector of the target drone at t = 0; v. of 0 Is the initial velocity vector of the target drone, h 0 Is the flight altitude of the target drone; l (0) = (x) j ,y j ,h 0 ) Represents x j Is the component of L (0) along the horizontal axis, y j Is the component of L (0) along the longitudinal axis and h 0 Is the component of L (0) in the vertical direction; l is a1j Is the component of L (t) along the horizontal axis, L b1j Is the component of L (t) along the longitudinal axis, L c1j Is the component of L (t) in the vertical direction; l is a1j 、L b1j And L c1j All follow v 0 、h 0 And t.
According to L (t) = (L) a1j ,L b1j ,L c1j ) Obtaining w from the horizontal velocity distribution data, the vertical velocity distribution data and the flying height distribution data j Corresponding set of coordinates of the point of falling COOR j =((a 1j ,b 1j ),(a 2j ,b 2j ),...,(a n(j)j ,b n(j)j ) ); wherein (a) ij ,b ij ) Is w j Corresponding ith drop location coordinate, i =1,2,3 j Corresponding place of fallingThe number of coordinates.
According to COOR j Determining P ris (x j ,y j ) (ii) a Wherein, P ris (x j ,y j ) Is w j A corresponding risk value;
；P 1 is the probability of the failure of the target unmanned aerial vehicle; p 2 (a ij ,b ij ) For the target unmanned plane at (x) j ,y j ) After accident, fall to (a) ij ,b ij ) The probability of (d); p is 2 (a ij ,b ij )=N ij /(n(j)),N ij Is COOR j Is located in (a) ij ,b ij ) The number of the coordinates of the falling position in the sub-area; p is 3 (a ij ,b ij ) For the target unmanned plane at (x) j ,y j ) Get out of affairs and fall to (a) ij ,b ij ) Probability of rear impact to the target; p is 4 (a ij ,b ij ) For unmanned plane in (x) j ,y j ) Is out of affairs and falls to (a) ij ,b ij ) And the probability of switching the target object from the normal state to the abnormal state after the target object is collided.
Will P ris (x j ,y j ) Corresponding risk level as P ris (x j ,y j ) A first target risk level for the corresponding sub-region.
According to another aspect of the present disclosure, there is also provided a non-transitory computer-readable storage medium, in which at least one instruction or at least one program is stored, and the at least one instruction or the at least one program is loaded and executed by a processor to implement the above method for acquiring risk level of a drone.
According to another aspect of the present disclosure, there is also provided an electronic device comprising a processor and the non-transitory computer-readable storage medium described above.
The invention has at least the following beneficial effects:
evaluating risk level of target area of target unmanned aerial vehicle, which can be based on differential model and target non-existenceDetermining a plurality of falling point coordinates corresponding to each characteristic point coordinate according to the horizontal speed distribution data, the vertical speed distribution data and the flight height distribution data of the man-machine, and determining P according to the plurality of falling point coordinates corresponding to each characteristic point coordinate 2 (a ij ,b ij )、P 3 (a ij ,b ij ) And P 4 (a ij ,b ij ) Further, a risk value corresponding to each feature point coordinate can be determined, and at the moment, a first target risk level corresponding to each sub-region can be determined; in a word, the data with more dimensionalities need to be considered for determining the first target risk level of each sub-area, the determination process of the differential model and the risk value is obtained based on the falling principle of the unmanned aerial vehicle, the falling condition of the unmanned aerial vehicle can be accurately reflected, and the accuracy of risk level assessment on the target area can be improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a flowchart of an unmanned aerial vehicle risk level obtaining method according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without inventive step based on the embodiments of the present invention, are within the scope of protection of the present invention.
The embodiment of the invention provides an unmanned aerial vehicle risk level acquisition method, wherein the method can be completed by any one or any combination of the following steps: a terminal, a server, or other devices with processing capability, which is not limited in this embodiment of the present invention.
In the embodiment of the present invention, a server is taken as an example, and a method for acquiring a risk level of an unmanned aerial vehicle will be described below with reference to a flowchart of a method for acquiring a risk level of an unmanned aerial vehicle shown in fig. 1.
The method comprises the following steps:
step S100, determining the central point coordinate of each sub-region in the target region D as a feature point coordinate, to obtain a feature point coordinate set W = (W) 1 ,w 2 ,w 3 ,...,w q ),w j =(x j ,y j )。
Wherein j =1,2, 3.. Q, q is the number of subregions, w j Is the jth feature point coordinate, x j Is w j Abscissa of (a), y j Is w j The ordinate of (a); d has a corresponding planar rectangular coordinate system.
In a possible implementation manner, a flight area corresponding to the target unmanned aerial vehicle may be used as a target area D, a plurality of sub-areas of D are obtained by equally dividing D, and when the server determines a feature point coordinate of each sub-area, the server may obtain coordinates of a plurality of boundary points on a boundary of the sub-area, then average values of abscissa of the plurality of boundary points are obtained, and the average value of the abscissa of the plurality of boundary points is set as the abscissa of the feature point coordinate of the sub-area, then average values of ordinate of the plurality of boundary points are obtained, and the average value of ordinate of the plurality of boundary points is set as the ordinate of the feature point coordinate of the sub-area, so as to obtain a feature point coordinate set W.
Step S200, obtaining a differential model m corresponding to the target unmanned aerial vehicle 0 *L"(t)=m 0 *g-c|L´(t)|L´(t)。
The unmanned aerial vehicle control method comprises the following steps of obtaining a target unmanned aerial vehicle, wherein t is the falling duration of the target unmanned aerial vehicle, L (t), L '(t) and L' (t) are all three-dimensional vectors, L (t) is a displacement vector of the target unmanned aerial vehicle corresponding to t, g is a gravity acceleration vector of the target unmanned aerial vehicle corresponding to t, L '(t) is a speed vector of the target unmanned aerial vehicle corresponding to t, and L' (t) is an acceleration vector of the target unmanned aerial vehicle corresponding to t; c is a resistance parameter c =0.5 × ρ AC of the target unmanned aerial vehicle D ρ is the air density constant; l' is LThe first derivative of (t), | L | (t) |, is the modulus of the first derivative of L (t), and L "(t) is the second derivative of L (t).
In a possible implementation, according to newton's second law, the corresponding differential model of the target drone may be determined.
To w j The following steps S300 to S600 are performed:
step S300, according to the differential model, L = v ″ (0) = v ″ 0 And L (0) = (x) j ,y j ,h 0 ) Obtaining L (t) = (L) a1j ,L b1j ,L c1j )。
Wherein, L ″ (0) = v ″ 0 And L (0) = (x) j ,y j ,h 0 ) All the initial conditions are corresponding to the differential model; l' (0) is a velocity vector of the target drone at t =0, and L (0) is a displacement vector of the target drone at t = 0; v. of 0 Is the initial velocity vector of the target drone, h 0 Is the flight altitude of the target drone; l (0) = (x) j ,y j ,h 0 ) Denotes x j Is the component of L (0) along the horizontal axis direction of the rectangular plane coordinate system, y j Is the component of L (0) in the direction of the longitudinal axis of the rectangular plane coordinate system and h 0 Is the component of L (0) in the vertical direction; l is a1j Is the component of L (t) along the direction of the horizontal axis of the rectangular plane coordinate system, L b1j Is the component of L (t) along the longitudinal axis of the rectangular plane coordinate system, L c1j Is the component of L (t) in the vertical direction; l is a1j 、L b1j And L c1j All follow v 0 、h 0 And t, are varied.
In a possible implementation, the server may respond to the differential model by applying the corresponding initial conditions L ″ (0) = v ″ 0 And L (0) = (x) j ,y j ,h 0 ) Solving the differential model can obtain L (t) = (L) a1j ,L b1j ,L c1j )。
Step S400, according to L (t) = (L) a1j ,L b1j ,L c1j ) Obtaining w from horizontal velocity distribution data, vertical velocity distribution data and flying height distribution data j Corresponding set of coordinates of the point of falling COOR j =((a 1j ,b 1j ),(a 2j ,b 2j ),...,(a n(j)j ,b n(j)j ))。
Wherein (a) ij ,b ij ) Is w j Corresponding ith drop location coordinate, i =1,2,3 j The number of corresponding drop location coordinates.
In one possible implementation, the horizontal velocity distribution data, the vertical velocity distribution data and the flying height distribution data can be normal distribution data, historical data of a plurality of unmanned aerial vehicles can be obtained through statistics, and the historical data are manually input into the server; the server is based on L (t) = (L) a1j ,L b1j ,L c1j ) And L c1j =0, can determine w j Corresponding time of falling to the ground t 1j Displacement vector L (t) corresponding to drop location 1j )=(L a2j ,L b2j ,L c2j )；L a2j Is L (t) 1j ) Component in the direction of the transverse axis, L b2j Is L (t) 1j ) Component in the direction of the longitudinal axis, L c2j Is L (t) 1j ) Component in the vertical direction, L a2j 、L b2j And L c2j All follow v 0 And h 0 Is changed;
randomly determining w according to the horizontal velocity distribution data and the vertical velocity distribution data j Corresponding set of target velocity vectors V = (V) 1 ,v 2 ,v 3 ,...,v n(j) ) (ii) a Wherein v is i The speed vector is an ith target speed vector in the V, and the target speed vector is any speed vector corresponding to any horizontal speed in the horizontal speed distribution data and any vertical speed in the vertical speed distribution data; i =1,2, 3.. N (j), n (j) being w j The number of corresponding target velocity vectors;
randomly determining a target flying height set H = (H) according to flying height distribution data 1 ,h 2 ,h 3 ,...,h n(j) ) (ii) a Wherein h is i The ith target flying height in the H is set, and the target flying height is any flying height in the flying height distribution data;
according to V, H and L (t) 1j ) Can determine w j Corresponding target displacement vector set DIS j =(L 1j ,L 2j ,L 3j ,...,L nj )；L ij Is DIS j The ith target displacement vector of (2), L ij To satisfy v 0 =v i And h 0 =h i L (t) of 1j )；
According to DIS j To obtain w j Corresponding coordinate set COOR of falling place j =((a 1j ,b 1j ),(a 2j ,b 2j ),...,(a n(j)j ,b n(j)j ) ); wherein (a) ij ,b ij ) Is w j Corresponding ith location coordinate, a ij Is L ij Component in the direction of the transverse axis, b ij Is L ij Component in the direction of the longitudinal axis.
In another possible embodiment, the method further comprises: obtaining a wind velocity vector v wind ；
To w j The following processing is also performed: determining t 1j Corresponding wind stroke vector L wind =v wind* t 1j ；
Based on this, the above a ij Is L ij Component in the direction of the horizontal axis and L wind Sum of components in the direction of the transverse axis, b ij Is L ij Component in the longitudinal direction and L wind The sum of the components in the direction of the longitudinal axis.
Step S500, according to COOR j Determining P ris (x j ,y j )。
Wherein, P ris (x j ,y j ) Is w j A corresponding risk value; p is ris (x j ,y j ) The following conditions may be satisfied:
；
P 1 is the probability of the failure of the target unmanned aerial vehicle; p is 2 (a ij ,b ij ) For the target unmanned plane at (x) j ,y j ) After accident, fall to (a) ij ,b ij ) Probability of (P) 2 (a ij ,b ij ) Is COOR j Is located in (a) ij ,b ij ) In the place ofA ratio of the number of drop location coordinates in the sub-region to n (j); p is 3 (a ij ,b ij ) For the target unmanned plane at (x) j ,y j ) Get out of affairs and fall to (a) ij ,b ij ) Probability of rear impact to the target; p is 3 (a ij ,b ij ) The following conditions are satisfied: p 4 (a ij ,b ij ) For unmanned aerial vehicle at (x) j ,y j ) Is in failure and falls to (a) ij ,b ij ) And the probability of switching the target object from the normal state to the abnormal state after the target object is collided.
In a possible implementation mode, the target object can be a person or a tree, and correspondingly, the target object is switched from a normal state to an abnormal state, namely the person dies or the tree is knocked down, and the like; the server is based on P 1 、P 2 (a ij ,b ij )、P 3 (a ij ,b ij ) And P 4 (a ij ,b ij ) Can determine w j Corresponding risk value P ris (x j ,y j )。
In another possible embodiment, P ris (x j ,y j ) The following conditions may also be satisfied:
；
P 1 (a ij ,b ij ) For the target unmanned plane in (a) ij ,b ij ) The probability of the sub-region in which it is in trouble.
Step S600, adding P ris (x j ,y j ) Corresponding risk level as P ris (x j ,y j ) A first target risk level for the corresponding sub-region.
In one possible implementation, a plurality of risk levels are stored in the server, each candidate risk level has a corresponding risk value range, and any two different candidate risk levels have no overlapping parts in the corresponding risk value ranges; based on this, the server can determine P among several candidate risk levels ris (x j ,y j ) Corresponding waitingSelecting a risk grade as P ris (x j ,y j ) A first target risk level for the corresponding sub-region; p is ris (x j ,y j ) Within a range of risk values corresponding to the first target risk level.
Therefore, the risk level of the target area of the target unmanned aerial vehicle is evaluated, a plurality of drop point coordinates corresponding to each characteristic point coordinate can be determined based on the differential model and the horizontal speed distribution data, the vertical speed distribution data and the flight height distribution data of the target unmanned aerial vehicle, and the P is determined according to the plurality of drop point coordinates corresponding to each characteristic point coordinate 2 (a ij ,b ij )、P 3 (a ij ,b ij ) And P 4 (a ij ,b ij ) Further, a risk value corresponding to each feature point coordinate can be determined, and at the moment, a first target risk level corresponding to each sub-region can be determined; in a word, the data with more dimensionalities need to be considered for determining the first target risk level of each sub-area, the determination process of the differential model and the risk value is obtained based on the falling principle of the unmanned aerial vehicle, the falling condition of the unmanned aerial vehicle can be accurately reflected, and the accuracy of risk level assessment on the target area can be improved.
Optionally, P 1 The following conditions are satisfied:
P 1 =N acc /N fly (ii) a Wherein, N acc Total number of crashes in D for several drones flying in D, N fly The total number of flights of a plurality of unmanned planes in D.
In one possible embodiment, the target drone is any one of a plurality of drones flying in D, and the plurality of drones may each be a drone of the same model as the target drone; the server can take the historical accident probability of a plurality of unmanned aerial vehicles in D as P 1 。
In another possible embodiment, P is as described above 1 (a ij ,b ij ) Can be (a) for a plurality of drones flying in D ij ,b ij ) The ratio of the total number of accident sub-areas to the total number of flight D of the unmanned aerial vehicles。
Optionally, P 3 (a ij ,b ij ) The following conditions are satisfied:
P 3 (a ij ,b ij )=s ij *（2(r p +r f )(h p /tan(γ ij ))+π(r p +r f ) 2 ) (ii) a Wherein r is p Is the average radius of the object, r f Maximum radius of target drone, h p Is the average height of the object, gamma ij For the target unmanned plane at (x) j ,y j ) After taking accident, fall to (a) ij ,b ij ) Angle of falling on the ground, gamma ij According to L (t) = (L) a1j ,L b1j ,L c1j ) Can be determined; s ij Is (a) ij ,b ij ) The density of the target in the sub-region.
In one possible embodiment, r p Can be the average radius of the human body, h p The target density can be the human mouth density, which is the average height of the human body; r is p 、h p 、r f And s ij All can be manually input into the server, gamma ij Can be according to the above L (t) 1j )=(L a2j ,L b2j ,L c2j ) And (4) determining.
Optionally, P 4 (a ij ,b ij ) The following conditions are satisfied:
P 4 (a ij ,b ij )= 
(ii) a Where α is a first constant threshold parameter, α =36, β is a second constant threshold parameter, β =100,e ij For the target unmanned plane at (x) j ,y j ) After taking accident, fall to (a) ij ,b ij ) Kinetic energy of impact in time; p is a radical of formula ij Is (a) ij ,b ij ) The shading factor of the sub-area; the shading factor is used to represent the degree of shading of the corresponding sub-region.
In one possible embodiment, α, β and p ij All can be manually input into the server, trees in a certain subareaThe more buildings and buildings are, the larger the corresponding shading factor is, and the more vacant land of a certain sub-area is, the smaller the corresponding shading factor is; e ij =0.5*m 0 *v t1j ，v t1j =L´(t 1j )。
Optionally, there are several no-fly regions in D; based on this, the method further comprises:
taking each subarea which has a superposition part with the plurality of no-fly areas as a first subarea;
taking each sub-area except the first sub-areas as a second sub-area;
taking the highest risk grade as a second target risk grade corresponding to each first sub-area;
taking the lowest risk grade as a second target risk grade corresponding to each second sub-area;
and for each sub-area, taking the highest risk grade in the corresponding first target risk grade and the second target risk grade as the first comprehensive risk grade corresponding to the sub-area.
In a possible implementation manner, the server may take each sub-region having an overlapping portion with the plurality of no-fly regions as a first sub-region according to manually input boundary data of each no-fly region and boundary data of each sub-region; the boundary data may be a plurality of discrete coordinate points on the boundary line of the corresponding region, or continuous data for representing the boundary line of the corresponding region; the server can determine a plurality of second sub-areas, then the highest candidate risk level in the plurality of candidate risk levels is used as a second target risk level corresponding to each first sub-area, and the lowest candidate risk level in the plurality of candidate risk levels is used as a second target risk level corresponding to each second sub-area; and finally, the higher one of the first target risk level and the second target risk level corresponding to each sub-area can be used as the first comprehensive risk level corresponding to the area. Therefore, the no-fly area can be set according to actual conditions, and the unmanned aerial vehicle flies in D according to the first comprehensive risk level of each sub-area, so that the flying safety of the unmanned aerial vehicle can be improved. Wherein any two different no-fly regions do not have overlapping portions.
Optionally, the step D has a plurality of preset static areas; based on the method, a plurality of no-fly areas are determined by the following method:
acquiring the minimum value of the flight height of the target unmanned aerial vehicle;
acquiring the maximum building height of each sub-area;
taking each sub-area with the corresponding building height maximum value smaller than the flight height minimum value as a dynamic area;
and taking each dynamic area and each static area as a no-fly area.
In a possible implementation manner, the server may determine boundary data of each dynamic area according to a manually input minimum value of the flight height of the target unmanned aerial vehicle, a manually input maximum value of the building height of each sub-area, and boundary data of each sub-area; the minimum value of the flight height can be the minimum flight height after the target unmanned aerial vehicle flies stably in the current flight plan; the server can also obtain manually input boundary data of a plurality of static areas, and then the boundary data of each dynamic area and the boundary data of each static area are both used as the boundary data of the no-fly area to determine the no-fly area in D.
Optionally, a plurality of control areas are arranged in the step D; based on this, the method further comprises:
taking each sub-area which has a superposition part with the plurality of control areas as a third sub-area; each control area has a corresponding preset risk level;
taking each sub-area except the third sub-areas as a fourth sub-area;
taking the preset risk level of the control area corresponding to each third sub-area as a third target risk level;
taking the lowest risk grade as a third target risk grade corresponding to each fourth subregion;
and regarding each sub-area, taking the highest risk grade in the corresponding first target risk grade, second target risk grade and third target risk grade as a second comprehensive risk grade corresponding to the sub-area.
In a possible implementation manner, one preset risk level is the same as one candidate risk level in the plurality of candidate risk levels, and the server may take each sub-region having a portion overlapping with the plurality of control regions as a third sub-region according to manually input boundary data of each control region, the preset risk level corresponding to each control region, and the boundary data of each sub-region; the server can determine a plurality of fourth sub-areas, then the preset risk level of the control area corresponding to each third sub-area is used as a third target risk level, and the lowest candidate risk level in the plurality of candidate risk levels is used as a third target risk level corresponding to each fourth sub-area; and finally, the highest one of the first target risk level, the second target risk level and the third target risk level corresponding to each sub-area can be used as a second comprehensive risk level corresponding to the area. Therefore, the setting of the no-fly area and the control area can be carried out according to actual conditions, and then the unmanned aerial vehicle flies in D according to the second comprehensive risk level of each sub-area, so that the flying safety of the unmanned aerial vehicle can be further improved. Wherein any two different regulatory regions do not have overlapping portions.
Embodiments of the present invention also provide a non-transitory computer-readable storage medium, which may be disposed in an electronic device to store at least one instruction or at least one program for implementing a method of the method embodiments, where the at least one instruction or the at least one program is loaded into and executed by a processor to implement the method provided by the above embodiments.
Embodiments of the present invention also provide an electronic device comprising a processor and the aforementioned non-transitory computer-readable storage medium.
Embodiments of the present invention also provide a computer program product comprising program code means for causing an electronic device to carry out the steps of the method according to various exemplary embodiments of the invention described above when said program product is run on the electronic device.
Although some specific embodiments of the present invention have been described in detail by way of example, it should be understood by those skilled in the art that the above examples are for illustration only and are not intended to limit the scope of the invention. It will also be appreciated by those skilled in the art that various modifications may be made to the embodiments without departing from the scope and spirit of the invention. The scope of the invention is defined by the appended claims.
Claims (7)
1. An unmanned aerial vehicle risk level obtaining method is characterized by comprising the following steps:
determining the central point coordinate of each sub-area in the target area D as a characteristic point coordinate to obtain a characteristic point coordinate set W = (W) 1 ,w 2 ,w 3 ,...,w q ),w j =(x j ,y j ) (ii) a Wherein j =1,2, 3.. Q, q is the number of sub-regions, w j Is the jth feature point coordinate, x j Is w j Abscissa of (a), y j Is w j The ordinate of (a);
acquiring differential model m corresponding to target unmanned aerial vehicle 0 *L"(t)=m 0 * g-c | L ″ (t); wherein t is an elapsed sinking time length during which the target unmanned aerial vehicle loses an wreck and sinks, L (t), L '(t) and L "(t) are all three-dimensional vectors, L (t) is a displacement vector of the target unmanned aerial vehicle corresponding to t, g is a gravity acceleration vector of the target unmanned aerial vehicle corresponding to t, L' (t) is a velocity vector of the target unmanned aerial vehicle corresponding to t, and L" (t) is an acceleration vector of the target unmanned aerial vehicle corresponding to t; c is a resistance parameter c =0.5 ρ AC of the target unmanned aerial vehicle D ρ is an air density constant;
to w j The following treatments were carried out:
according to the differential model, L = v ″ (0) = v) 0 And L (0) = (x) j ,y j ,h 0 ) Obtaining L (t) = (L) a1j ,L b1j ,L c1j ) (ii) a Wherein, L ″ (0) = v ″ 0 And L (0) = (x) j ,y j ,h 0 ) All the initial conditions are corresponding to the differential model; a velocity vector of the target drone when L' (0) is t =0, and when L (0) is t =0A displacement vector of the target drone; v. of 0 Is the initial velocity vector, h, of the target drone 0 Is the flight altitude of the target drone; l (0) = (x) j ,y j ,h 0 ) Denotes x j Is the component of L (0) along the horizontal axis, y j Is the component of L (0) in the direction of the longitudinal axis and h 0 Is the component of L (0) in the vertical direction; l is a radical of an alcohol a1j Is the component of L (t) along the horizontal axis, L b1j Is the component of L (t) along the longitudinal axis, L c1j Is the component of L (t) in the vertical direction; l is a1j 、L b1j And L c1j Are all following v 0 、h 0 And t;
according to L (t) = (L) a1j ,L b1j ,L c1j ) Obtaining w from horizontal velocity distribution data, vertical velocity distribution data and flying height distribution data j Corresponding coordinate set COOR of falling place j =((a 1j ,b 1j ),(a 2j ,b 2j ),...,(a n(j)j,bn(j)j ) ); wherein (a) ij ,b ij ) Is w j Corresponding ith drop location coordinate, i =1,2,3 j The number of corresponding drop location coordinates;
according to COOR j Determining P ris (x j ,y j ) (ii) a Wherein, P ris (x j ,y j ) Is w j A corresponding risk value;
；P 1 the probability of the target unmanned aerial vehicle failing is set; p 2 (a ij ,b ij ) For the target drone in (x) j ,y j ) After taking accident, fall to (a) ij ,b ij ) The probability of (d); p 2 (a ij ,b ij )=N ij /(n(j)),N ij Is COOR j Is located in (a) ij ,b ij ) The number of the coordinates of the falling position in the sub-area; p is 3 (a ij ,b ij ) For the target drone in (x) j ,y j ) Is out of affairs and falls to (a) ij ,b ij ) Rear impact to targetThe ratio; p is 4 (a ij ,b ij ) For unmanned aerial vehicle at (x) j ,y j ) Is in failure and falls to (a) ij ,b ij ) And the probability that the target object is switched from a normal state to an abnormal state after the target object is collided is obtained;
will P ris (x j ,y j ) Corresponding risk level as P ris (x j ,y j ) A first target risk level for the corresponding sub-region;
P 3 (a ij ,b ij ) The following conditions are satisfied:
P 3 (a ij ,b ij )=s ij *（2(r p +r f )(h p /tan(γ ij ))+π(r p +r f ) 2 ) (ii) a Wherein r is p Is the average radius of the object, r f Is the maximum radius, h, of the target drone p Is the average height of the object, gamma ij For the target drone in (x) j ,y j ) After taking accident, fall to (a) ij ,b ij ) Angle of falling on the ground, gamma ij According to L (t) = (L) a1j ,L b1j ,L c1j ) Can be determined; s ij Is (a) ij ,b ij ) The density of the target in the sub-area;
P 4 (a ij ,b ij ) The following conditions are satisfied:
P 4 (a ij ,b ij )=
(ii) a Where α is a first constant threshold parameter, α =36, β is a second constant threshold parameter, β =100,e ij For the target drone in (x) j ,y j ) After taking accident, fall to (a) ij ,b ij ) Kinetic energy of impact in time; p is a radical of ij Is (a) ij ,b ij ) The shading factor of the sub-area; the shading factor is used to represent the degree of shading of the corresponding sub-region.
2. The method of claim 1, wherein the first and second light sources are selected from the group consisting of,characterised by the feature that P 1 The following conditions are satisfied:
P 1 =N acc /N fly (ii) a Wherein N is acc Total number of crashes in D for several drones flying in D, N fly The total number of times that the unmanned aerial vehicle flies in D.
3. The method of claim 1, wherein D has several no-fly regions therein;
the method further comprises the following steps:
taking each sub-area which has an overlapped part with the plurality of no-fly areas as a first sub-area;
taking each sub-area except a plurality of first sub-areas as a second sub-area;
taking the highest risk level as a second target risk level corresponding to each first sub-area;
taking the lowest risk level as a second target risk level corresponding to each second sub-area;
and for each sub-area, taking the highest risk grade in the corresponding first target risk grade and second target risk grade as the first comprehensive risk grade corresponding to the sub-area.
4. The method according to claim 3, wherein D has a preset number of static areas;
several of the no-fly zones are determined by the following method:
acquiring the minimum value of the flight height of the target unmanned aerial vehicle;
acquiring the maximum building height of each sub-area;
taking each sub-area with the corresponding building height maximum value smaller than the flight height minimum value as a dynamic area;
and taking each dynamic area and each static area as the no-fly area.
5. The method of claim 4, wherein D has a number of regulatory domains;
the method further comprises the following steps:
taking each sub-area which has a superposition part with the plurality of control areas as a third sub-area; each control area is provided with a corresponding preset risk level;
taking each sub-area except a plurality of the third sub-areas as a fourth sub-area;
taking the preset risk level of the control area corresponding to each third sub-area as a third target risk level;
taking the lowest risk level as a third target risk level corresponding to each fourth sub-area;
and for each sub-area, taking the highest risk grade in the corresponding first target risk grade, second target risk grade and third target risk grade as the second comprehensive risk grade corresponding to the sub-area.
6. A non-transitory computer readable storage medium having stored therein at least one instruction or at least one program, the at least one instruction or the at least one program being loaded and executed by a processor to implement the method of any one of claims 1-5.
7. An electronic device comprising a processor and the non-transitory computer readable storage medium of claim 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211225542.XA CN115294487B (en) | 2022-10-09 | 2022-10-09 | Unmanned aerial vehicle risk level obtaining method, storage medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211225542.XA CN115294487B (en) | 2022-10-09 | 2022-10-09 | Unmanned aerial vehicle risk level obtaining method, storage medium and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115294487A CN115294487A (en) | 2022-11-04 |
| CN115294487B true CN115294487B (en) | 2023-01-31 |
Family
ID=83834402
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211225542.XA Active CN115294487B (en) | 2022-10-09 | 2022-10-09 | Unmanned aerial vehicle risk level obtaining method, storage medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115294487B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116451588B (en) * | 2023-04-25 | 2024-02-27 | 中航信移动科技有限公司 | Method, medium and equipment for determining prompt information based on target object predicted track |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111177851A (en) * | 2019-12-27 | 2020-05-19 | 北航（四川）西部国际创新港科技有限公司 | Method for evaluating ground risks in unmanned aerial vehicle operation safety risk evaluation |
| CN111582740A (en) * | 2020-05-13 | 2020-08-25 | 电子科技大学 | Risk assessment system for multi-rotor unmanned aerial vehicle |
| CN111639393A (en) * | 2020-05-19 | 2020-09-08 | 广州大学 | Unmanned aerial vehicle safety early warning method for detecting multiple risk factors based on differential game |
| CN112053593A (en) * | 2020-06-23 | 2020-12-08 | 中国民用航空总局第二研究所 | Airport unmanned aerial vehicle control area planning method based on risk assessment |
| CN112132144A (en) * | 2020-07-31 | 2020-12-25 | 北京航空航天大学 | Unmanned aerial vehicle air line ground collision risk assessment method based on remote sensing image |
| CN113361174A (en) * | 2021-06-17 | 2021-09-07 | 河北科技大学 | STP model-based large unmanned aerial vehicle collision probability calculation method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA3059926A1 (en) * | 2019-05-24 | 2020-11-24 | A3 By Airbus, Llc | Risk assessment for unmanned aerial vehicles |
-
2022
- 2022-10-09 CN CN202211225542.XA patent/CN115294487B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111177851A (en) * | 2019-12-27 | 2020-05-19 | 北航（四川）西部国际创新港科技有限公司 | Method for evaluating ground risks in unmanned aerial vehicle operation safety risk evaluation |
| CN111582740A (en) * | 2020-05-13 | 2020-08-25 | 电子科技大学 | Risk assessment system for multi-rotor unmanned aerial vehicle |
| CN111639393A (en) * | 2020-05-19 | 2020-09-08 | 广州大学 | Unmanned aerial vehicle safety early warning method for detecting multiple risk factors based on differential game |
| CN112053593A (en) * | 2020-06-23 | 2020-12-08 | 中国民用航空总局第二研究所 | Airport unmanned aerial vehicle control area planning method based on risk assessment |
| CN112132144A (en) * | 2020-07-31 | 2020-12-25 | 北京航空航天大学 | Unmanned aerial vehicle air line ground collision risk assessment method based on remote sensing image |
| CN113361174A (en) * | 2021-06-17 | 2021-09-07 | 河北科技大学 | STP model-based large unmanned aerial vehicle collision probability calculation method |
Non-Patent Citations (1)
| Title |
|---|
| 无人机系统安全目标水平预估方法;张泽京等;《航空动力学报》;20180423;第33卷(第04期);1017-1024 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115294487A (en) | 2022-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109522784B (en) | Device and method for distinguishing between surmountable and non-surmountable objects | |
| DE102014114827A9 (en) | Path planning for evasive steering maneuvers in the presence of a target vehicle and surrounding objects | |
| WO2020007189A1 (en) | Obstacle avoidance notification method and apparatus, electronic device, and readable storage medium | |
| CN115294487B (en) | Unmanned aerial vehicle risk level obtaining method, storage medium and electronic equipment | |
| US11892845B2 (en) | System and method for mission planning and flight automation for unmanned aircraft | |
| DE102011100927A1 (en) | Object and vehicle detection and tracking using 3-D laser rangefinder | |
| CN112132144B (en) | Unmanned aerial vehicle air line ground collision risk assessment method based on remote sensing image | |
| US10254758B2 (en) | Object tracking by unsupervised learning | |
| EP3698227B1 (en) | Path planning for an unmanned vehicle | |
| US20230053459A1 (en) | Vehicle-based data processing method and apparatus, computer, and storage medium | |
| CN109633631A (en) | A kind of multi-functional reconnaissance radar combat duty automatic planning | |
| DE102015224796A1 (en) | Method and control unit for detecting a possible collision of an unmanned aerial vehicle with an object | |
| CN115187005B (en) | Scheduling method, device, equipment and storage medium | |
| US20220019713A1 (en) | Estimation of probability of collision with increasing severity level for autonomous vehicles | |
| DE102019216722A1 (en) | Method for locating a vehicle on a digital map | |
| CN110709791A (en) | Method and device for adjusting flying height and controlling flying of unmanned aerial vehicle | |
| CN113609784B (en) | Traffic limit scene generation method, system, equipment and storage medium | |
| US20190147749A1 (en) | System and Method for Mission Planning, Flight Automation, and Capturing of High-Resolution Images by Unmanned Aircraft | |
| CN111081061B (en) | Collision early warning method and device | |
| CN112001309A (en) | Target searching method, device and equipment based on unmanned aerial vehicle cluster and storage medium | |
| KR102314515B1 (en) | Method for providing autonomous driving service platform to be used for supporting autonomous driving of vehicles by using competitive computing and information fusion, and server using the same | |
| CN112166446A (en) | Method, system, device and computer readable storage medium for identifying trafficability | |
| CN115480582A (en) | Maneuver prediction method for LSTM-based target, electronic device and storage medium | |
| CN116483127A (en) | Unmanned aerial vehicle off-site take-off and landing method, unmanned aerial vehicle control terminal and storage medium | |
| US20210225179A1 (en) | Apparatus, method and system relating to aircraft systems |
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 |