1. Introduction
As a key component of advanced air mobility (AAM), unmanned aircraft systems (UAS) provide the opportunity to improve understanding of the very low-level (VLL) airspace resource [
1,
2]. Progress in the development of UAS technology is raising convenience to new heights [
3]. For instance, package delivery by UAS, or drones, has emerged to address the issue of last-mail delivery. Many conventional logistics and other start-up companies are developing drone business lines. In 2019, Antwork (a start-up company) received the first approved urban drone operation license from the Civil Aviation Administration of China (CAAC) [
4]. In 2021, SF Express (a Chinese delivery company) also received approval for its suburb delivery service. At the same time, MeiTuan (a technology-driven retail company) began to deliver take-out food using multi-rotor drones in Shenzhen, China [
5,
6]. With widespread use of UAS in various areas of industry, the number of drones (especially for recreation) is rapidly increasing. According to a statistical report published by the CAAC [
7,
8], the registered number of drones is now 832,000, an increase of 58.8 percent from last year. In the USA [
9], drones are also rapidly increasing in number and complexity. There are currently 865,505 registered drones, including 314,689 commercial drones, 538,172 recreational drones and 3644 paper registrations.
Although drones can bring greater convenience, the risks associated with their operation are gradually increasing [
10]. With continuous growth in the number of UAVs, the operational risks associated with UAVs threaten the overall security of the airspace. In recent years, the impact of unauthorized UAVs on the normal operation of airports has been growing, as illustrated by the Heathrow Airport UAV incident in 2019. As most UAVs are required to operate in very-low-altitude segregated airspace, if they crash to the ground, there is potential for severe injury to people and damage to property. To ensure the safety of people on the ground, it is crucial to precisely estimate ground risk in terms of relevant parameters (falling range, people density, etc.) before planning a UAV route.
To better understand ground risk, a review of the literature was conducted which was divided into two parts: theoretical and practical.
From a theoretical point of view, ground risk assessment is closely related to calculation of a drone’s ground impact. Many studies have estimated the degree of risk by calculating the magnitude of ground impact. Moha and Hu [
11] applied ADAMS (automatic dynamic analysis of mechanical systems) and MATLAB co-simulation methods to evaluate the impact of different UAV failure types on the crash trajectory and crash area. Andrew [
12] built a kinetic energy model to evaluate the ground risk buffer. Both these studies demonstrated that the size of the ground risk buffer is dynamic and depends mainly on the height, velocity, and mass of the operation. These two studies each provided ways to calculate ground risk. Noting that the ground impact is strongly related to descent trajectories, Anders [
13] presented a probability density function based on a second-order drag model to predict the trajectory of fall drones. Baptiste [
14,
15,
16] proposed two approaches (a DoF dynamic model and a stochastic approach) to simulate the descent states of fixed-wing drones and focused on ground impact probabilities within 95%, 99% and 99.9% limits. It was found to be possible to accurately predict the ground impact, but the approach used required large storage and computational resources, which may limit its practical application. To address these issues, Liu [
17] simplified the model and analyzed associated uncertainties, such as wind and airspeed, that can affect UAV descents.
In much of the published research, ground risk is considered alongside path planning. Many studies have focused on determining the nature of ground risk for different paths. Hu and Pang [
18,
19,
20] introduced an integrated risk assessment model to plan the path of a UAV while minimizing operational risk. These authors combined parameters such as population density, impact area, and fatality rate to calculate a risk index under different paths. Results from a case study demonstrated that the risk assessment model was able to identify high-risk areas, and that the risk map produced enabled safe UAV path planning. Stefano [
21,
22] proposed a risk map to quantify the risk associated with an operational area based on the population density, the shelter factor and no-fly zones. The studies described have helped to clarify the relationship between risk and path planning but ignore the temporal dynamics of risk. For example, assessment of risk needs to consider population density, but most population data used in current methods are static data, so that the models used cannot capture dynamic trends in population change. Therefore, the risk value determined by the model is very different from the actual risk, so that the model used is unable to meet actual operational needs. Hence, there is now an increasing focus on the spatiotemporal characteristics of risk associated with drone use. Aliaksei [
23] used spatiotemporal population density data to assess the ground risk of drone operations. This method can enable operators to assess risk more precisely at different times. With the development of artificial intelligence (AI), deep learning methods have been used extensively to predict the density of people based on spatiotemporal characteristics. Nabeela [
24] and Zhang [
25,
26] forecast citywide crowd flows using a deep learning method. The results showed that crowd flows can be predicted using information on the time period and trends by incorporating data on crowd closeness. Using the findings of previously published studies on crowd flows, a new perspective was provided enabling assessment of dynamic changes in the size of the ground risk buffer and evaluation of the corresponding ground risk by predicting the density of people.
From the perspective of practical applications, aviation bureaus (such as EASA, CAAC and FAA) have begun to publish regulations and guidelines to ensure safer, more efficient and sustainable operations [
27,
28,
29]. With respect to drone safety management, the CAAC and EASA have applied a specific operations risk assessment (SORA) methodology [
30] which was proposed by the Joint Authorities for Rulemaking on Unmanned Systems (JARUS) to evaluate and manage the operational risks of the use of drones. Using the SORA approach, risk is divided into ground risk and air risk. Ground risk is determined by parameters such as mass, the velocity of drones, the ground environment and mitigation methods (people, ground risk buffer, and so on). The SORA approach represents one of the most comprehensive methods for providing effective guidance for aviation bureaus to better manage the operation of drones. However, the SORA method is largely qualitative. For instance, the ground risk buffer (the area in which people are protected from falling drones) is one of the most important ground risk criteria because it not only determines the level of UAV safety, but also has a direct impact on operational efficiency. The standard SORA method for assignment of the ground risk buffer follows a “one-to-one” rule in which the minimal size of the ground risk buffer is equal to the height of the operation. This rule may, however, be too rigid and inflexible for actual operations. Some aviation bureaus, such as the FAA and the civil aviation department of Hong Kong, China (HKCAD) have evaluated the probability and severity of drone accidents to determine the operational risk. This method is similar to the traditional civil aviation safety management approach and is, therefore, convenient for use by safety departments as they do not need become familiar with a new method. Nevertheless, identifying all accident hazards and accurately calculating their probability of occurrence remains a challenge due to the lack of adequate operational data.
Whichever of the two methods described above is used, they each have clear drawbacks. The current method of assessing ground risk is highly reliant on a large number of experienced experts and the use of many assurance documents, which can be time-consuming and lack a logical or scientific basis. Moreover, operational risk may be mainly affected by dynamic parameters, such as people density, wind, airspeed, and weight of the drones, and the current method has difficulty calculating the size of the buffer quickly and accurately. Furthermore, the rate of utilization of airspace will be reduced if a qualitative risk evaluation method is used, with the result that the operational resources are insufficient for the increasing number of drones.
The existing methods cannot conduct real-time assessments of UAV operation risk, especially because of the large gap between the risk during assessment and the risk during actual operation due to the dynamic characteristics of some parameters. Therefore, how to accurately predict future risk remains a challenge.
Hence, we propose a deep-learning-based method to calculate the dynamic ground risk. The model is divided into three parts: we first analyze the dynamic parameters that may affect ground risk and build a kinetic model to predict the size of the ground impact. Second, we build a joint convolutional neural network and deep neural network model (C-Snet model) to predict people density and to calculate the casualty rate of a ground impact. Third, a dynamic model combining a deep learning and a kinetic model is established to predict the ground risk. The contributions of this research include the following:
- (1)
In contrast to previous methods, our approach offers a new approach to the prediction of ground risk under uncertain parameters based on a joint deep learning method and kinetic model.
- (2)
We analyze the uncertainty of the dynamic parameters and show the corresponding relationships to ground risk.
- (3)
We consider the relationship between the ground risk and the casualty rate using the model.
The remainder of this paper is structured as follows:
Section 2 analyzes the parameters that affect ground risk and builds a kinetic model to evaluate the size of the ground risk buffer.
Section 3 builds the deep learning-kinetic model to predict the corresponding people density and the size of ground impact and then calculates the operation ground risk level.
Section 4 describes the experiments undertaken to evaluate the robustness of the model and analysis of their results.
Section 5 summarizes the conclusions and discusses potential further work.
2. Parameters That Affect Ground Risk
Based on the literature review, the methods used to evaluate ground risk involve consideration of a range of parameters, including people density, shelter, ground impact and others [
31]. In this section, we first define the method for calculating the ground risk.
where
is the probaility of ground risk,
A is the area of impact point in m
2,
pop is the corresponding people density of the ground impact,
is the probability of the fallen drone, and
is the fatality rate of the people impacted by the fallen drones.
The impact of various parameters on ground risk in the context of actual operations is discussed below.
2.1. Ground Impact
The ground impact is defined as the landing point of the UAV after control is lost. This factor plays an important role in the design of a ground risk buffer. The ground impact is positively correlated with the falling range of the UAV [
32]—the further the ground impact, the larger the ground risk buffer. The size of the ground risk buffer can be modified by calculating the ground impact. The ground impact can be assessed in terms of a ballistic descent model. The ballistic descent model represents the descent of the drone following loss of control and assumes zero lift-drag. Hence, the ballistic descent of a UAV can be considered as a kind of horizontal projectile motion—free fall motion in the vertical direction and variable speed linear motion in the horizontal direction. The model is based on the standard second-order drag model.
where
m and
v are the mass and velocity of the drones,
g is the gravity, and c is a constant that encapsulates the drag coefficient.
However, the ballistic descent model may also be considered to represent a probabilistic problem due to uncertainties of wind and track error.
Track error: The GPS positioning accuracy of a UAV is affected by weather, obstacles, and its own flight control system. This can cause deviation between the actual flight path and the planned flight path. Under this circumstance, when calculating the ground risk impact of a UAV, the track error will interfere with the accuracy of the position, which makes it impossible to accurately calculate the falling range of the UAV, and results in large errors in the definition of the ground risk buffer.
The UAV track error analysis experiments undertaken included the following stages: first, a three-dimensional motion capture system for an indoor UAV, based on VICON optical position sensing, was set up in a laboratory environment. Second, a UAV was used in a flight experiment in the laboratory. Finally, the flight path position data output by the UAV flight control system was compared with the vision position data calibrated in the flight path capture system and the distributions of the initial position errors and velocity errors were analyzed.
According to Han [
31], aircraft track error follows a three-dimensional Gaussian distribution with a mean value of 0, and the track error is independent in the longitudinal, lateral, and vertical directions, as represented in the following formula:
where
P is the distribution probability column vector of the track error,
E is the aircraft track point error column vector,
M is a column vector of the three-dimensional aircraft track point error mean value, and
is the covariance matrix of the track error.
Wind: The relationship between the height and wind speed is analyzed. Because the wind speed is different at different heights, the incoming flow velocity is also different. However, at present, most of the data obtained by UAV operators are the wind speed values from ground observation stations. If the air resistance is calculated based on the wind speed, it will often cause large errors. Therefore, to ensure accurate calculation of the range of the UAV landing points, it is necessary to accurately calculate the wind speed at each altitude. Considering that the flight height of a UAV is mostly within 1000 m, this paper uses a near surface wind profile model to calculate the wind speed at different heights.
where
is the average wind speed that needs to be calculated,
h stands for the height at which the wind speed is calculated,
is the height at which the wind speed is observed,
is the observed wind speed of height
, and
n is the variation coefficient of wind speed with height.
2.2. People Density
People density is one of the most important factors in a risk assessment because it describes how many people may be involved in a drone crash and how they are distributed on the map. At present, the risk assessment for a UAV is still based on static statistical data as the basis for analysis of ground risk. Real-time assessment of UAV is not carried out according to the spatial-temporal characteristics of people density, which can lead to a significant difference between the assessed risk and the risk during actual flight. In particular, the population density is highly dynamic data, while the probability of crowd density for different times of operation is quite different. Therefore, accurate population density predictions are required for proper risk assessment. As spatial-temporal data, people density has temporal characteristics including closeness, trend, and period components and spatial characteristics including nearby and distance components [
33].
2.2.1. Temporal Characteristics
Closeness is expressed as the data for two adjacent time nodes. Crowd density data are often similar—for example, the crowd density at 9 a.m. is similar to that at 10 a.m.
Period refers to similar characteristic values for population density at certain time intervals. For example, the crowd density at 8 a.m. on Monday is similar to the crowd density at 8 a.m. on Tuesday. This similarity feature requires statistical data and analysis over time periods.
The figure above shows the people density for a whole week. From
Figure 1, the people density has the characteristic of period; the population density at the same time but on different days has strong similarity (the cosine similarity value is larger than 0.90).
Trend is expressed as the change in season or year. The density index will also change and the number of people will also exhibit a certain change trend.
Figure 2 below shows that the population density was lowest at 4:00 a.m. With elapsed time, the value continued to increase and reached a peak at 9:00 a.m.
2.2.2. Spatial Characteristics
Nearby. Defined as the inflow of region A, is affected by the outflow of nearby regions. According to geographic theory [
34,
35], a city can be divided using a grid map based on longitude and latitude where a grid denotes a region. The closer the regions are, the greater their connection.
Distance. The flow can be affected by the distance of the regions. For example, people who live far away from a central business district always take the subway or highway to work, which means that the outflow from distance areas directly affects the inflow of the central business district.
2.3. Shelter Factor
The shelter factor is defined as the degree of protection against injury to people afforded by buildings and trees on the ground after a UAV falls. Similar to the ground protection area, the shelter factor frequently takes into account the coverage provided by the buildings and trees under the path. The ground protection area is a two-dimensional region, whereas the shelter factor is a three-dimensional area. By estimating the size of nearby buildings and trees, the extent of harm that a falling UAV may cause people can be assessed. The term “shelter factor” has been defined differently in several studies. According to Guglieri [
36], there are 11 levels of ground cover from 0 to 10. Zero represents no protection and 10 is the strongest protection. However, Dalamagkidis [
37] suggested that the range of the shelter factor was between zero and positive infinity. Positive infinity denotes the strongest possible protection, whereas zero denotes no protection at all. This paper takes into account the actual computation and application requirements. As a result, the definition of the shelter factor ranges from 0 to 1.
Table 1 provides details of the shelter factor levels.
2.4. Fatality Rate
The fatality rate is defined as the probability of fallen drones impacting a person and producing fatal injuries. It is highly dependent on the impact energy and shelter factors. According to Dalamagkidis [
37], the fatality rate can be evaluated using the shelter factor, impact energy, etc. The expression for the fatality rate is as follows.
where
is the fatality rate,
is the impact energy of fallen drones,
is the shelter factor,
stands for the impact energy required for 50% mortality when
= 0.5, and
stands for the impact energy limit required to cause death when
falls to 0. According to research undertaken by the American Range Command Commission (RCC) (RCCDocument323-99 issued in 1999 and RCCDocument321-07 issued in 2007) [
19], if a 0.00454 kg (1lb) object strikes a human body, the corresponding impact kinetic energy for 10% and 90% fatalities is 50 J and 200 J. Therefore with respect to the operation of drones, we reference Hu [
19] and set
= 10
6 and
= 100 J. The larger the shelter factor, the lower the fatality rate for the same impact kinetic energy.
3. Method of Dynamic Ground Risk Assessment
In this section, a dynamic ground risk model under uncertainty parameters is established. The dynamic ground risk model has two main components: a kinetic model and a joint deep learning model. First, the size of the ground impact and the fatality rate are regarded as dynamic parameters and are calculated by the kinetic model. Second, the joint deep learning model is built to predict the people density and to evaluate the shelter factor. Finally, the dynamic ground risk model is obtained by multiplying the output of the kinetic model and the joint deep learning model.
3.1. Kinetic Method to Assess Ground Impact
The mechanical model for the quadrotor UAV coordinate system is shown in
Figure 3 below. Four evenly distributed rotors produce corresponding speed, thrust and torque [
32]. In addition, the UAV fuselage is subject to uniformly distributed gravity and air drag that need to be calculated (
Figure 4).
According to Newton’s second law, the UAV crash motion can be decomposed into vertical and horizontal directions with the following equations (Equation (5)) of motion:
where
m is the mass of the drone,
G is the gravity, and
Dx and
Dz are the air drag of the different sides, respectively.
ax and
az are the accelerated velocities of the horizontal and vertical directions, respectively. The expression can be further decomposed as follows (Equation (6)):
where
m is the mass of the drone,
Cd is the air drag coefficient, and
Ax and
Az are the frontal areas of the horizontal and vertical directions, respectively.
x and
z are the displacements of the longitudinal and vertical directions, respectively.
With initial conditions
vx0 =
v0,
vz0 = 0 and
x0 = 0,
z0 = 0, the expression for the ground impact point can be obtained:
Considering that the wind may have an impact on ground impact, the expression relative to the projection of the event point onto the ground becomes the vector [
31].
3.2. Joint Deep Learning Model to Predict People Density
In
Section 2.2, we discussed the characteristic of people density and analyzed the relationship between people density and ground risk. The reason for using a deep learning model to predict people density is as follows:
From the perspective of drone risk management, unmanned aircraft pose a high risk to the safety of people on the ground due to very low altitude operation. The operator of drones needs to assess the risks before operation. One of the key parameters for risk assessment is people density. However, how to accurately predict people density and how to calculate ground risk using an online method represent current challenges.
From the perspective of path planning, the operator needs to calculate the ground risk when planning the flight path of drones according to the regulations of the Civil Aviation Administration of China. If the people density of the ground is large, the path must be redesigned to decrease the risk. However, from geography theory, the people density is not static but has spatial and temporal properties, which means that the dynamic path must be designed according to the dynamic character of people density.
From the perspective of operation authorization, the bureau needs to know the operational environment (especially people density) of the drones; hence, a deep learning model can be utilized to better explain the ground environment to the bureau and to let the inspector know the people density and risk in the next couple of hours.
Therefore, we decided to use a deep learning model to predict the people density to solve the problem of uncertainty in the people density parameter. Using a deep learning model, the people density of an operational area can be predicted for the next few hours and the operator can predict the precise operational risk accordingly.
According to deep learning theory, an LSTM model is suitable for time series prediction, while a CNN model is suitable for discrete data prediction. Therefore, considering the spatial-temporal characteristics of population density, we decided to establish a CNN-LSTM joint model to improve the accuracy of prediction. We input spatial features, such as nearby and distance features, into the CNN model and temporal features, such as closeness, period and trend into the LSTM. Finally, the two models jointly output the predicted population density. Our model can not only take into account the space-time characteristics of population density, but can also improve the efficiency of computing through joint computing.
3.2.1. Convolutional Neural Network
The convolutional neural network sub-model consists of a convolutional layer, a pooling layer and a fully connected layer (see
Table 2 and
Figure 5) [
38,
39,
40,
41]. The pooling layer uses MaxpoLling, which can shield the unimportant parameters while maintaining the data characteristics to solve the problem of high data redundancy in the model. For convolution layer setting, because the data in this study are discrete data and are not sensitive to periodic change in time, the horizontal sliding value and vertical sliding value of the two convolution layers of the model are set to 1, and the convolution operation is performed using padding for the same 0 filling. This paper uses the Relu function as the sub-model activation function, the convolutional layer of the convolutional neural network and the depth of the linear neural network layer results of non-linear mapping to avoid the model gradient explosion and the disappearance of the gradient problem. The formula is as follows (Equation (9)):
The input parameters for the CNN are shown in Equation (10) below,
where
is the crowd density of the same area in the last week,
is the date of week to which the crowd density of the area belongs, and
is the time corresponding to the crowd density. From
Figure 6, these three parameters are input to a convolutional neural network and the predicted crowd density index is obtained using Adam’s optimization algorithm.
3.2.2. LSTM Neural Network
Long short-term memory (LSTM), a type of recurrent neural network (RNN), can learn long-term dependencies, especially in sequence-prediction problems [
42]. LSTM overcomes the drawbacks of long-term prediction instability in the RNN algorithm. Like RNN, LSTM also has a chain structure to repeat and memorize information; each structure is called a cell. However, instead of having a single neural network layer in RNN, LSTM consists of four parts: a cell state, a forget gate layer, an input gate layer, and an output gate layer (
Figure 7).
where
is the sigmoid activation, which outputs a number between 0 and 1 for each piece of information in the cell state. If the sigmoid value is 1, then the information is completely retained. Similarly, the information is completely deleted when the sigmoid value is 0.
is the input information of the cell.
is the short-term state of the former cell.
,
,
are the weight matrices of the hidden layer input to the three gate layers.
,
,
,
are the biases of the layer function. The layer concatenates the input information
and the short-term state of the former cell
using the sigmoid function
to determine whether the information should be kept.
The input parameters of the LSTM model are as follows:
where (
t − 3) is the hour’s people density, (
t − 2) the hour’s people density, and (
t − 1) the hour’s people density at each of these time points, respectively.
From
Figure 8, the LSTM model has three layers. The number of neurons in each of the three layers is 256, 128 and 64, respectively. The model uses Relu as the activation function and MSE as the loss function. The look-back value of this model is three. Finally, the model outputs the predicted people density.
3.2.3. Architecture of Joint Deep Learning Model
Some characteristics of the people density (such as nearby and distance) cannot be predicted well by the LSTM model, while the CNN model may not be suitable for representing closeness and trend in people density. Therefore, we decided to combine the CNN and LSTM models to predict people density to improve the prediction accuracy of the model.
The framework for the joint deep learning model (CS-net model) is shown in the figures (
Figure 9 and
Figure 10) below. The model is composed of a CNN component and an LSTM component. A regression function is set at the top of the model to combine the output from the CNN and LSTM components. We first put the spatial characteristics into the CNN model. In the CNN model, the parameter features use convolutional units to expand the model’s ability to extract feature information and ensure the precision of the model. At the same time, the temporal characteristics are put into the LSTM model for fitting training. The model consists of four parts. The activation function of each hidden layer is the sigmoid. Each neuron in the adjacent layer is fully connected. The CS-net model uses the mean square error (MSE) as the loss function of the top-level regression function of the model. When performing fitting training, the model needs to consider the loss functions of both the CNN and LSTM components. In addition, the model must consider the fitting rates to optimize all parameters. Joint training of a CS-net model is achieved by back-propagating the gradients from the output to both the CNN and LSTM components of the model simultaneously using mini-batch stochastic optimization. Compared with embedding learning, the CS-net model training is not independent and only requires a small number of feature types to obtain higher accuracy. All the parameters in the model are jointly trained under the mean square error (MSE) loss. Due to the combination of two modules, it is quite difficult to determine a proper learning rate to train the joint model. Therefore, we chose Adam, a stochastic gradient descent method with an adaptive step size and momentum to optimize the model. The formula (Equation (17)) for the model prediction result is:
From the time trend characteristics of crowd density, the crowd flow in the last two hours has a great influence on the crowd density in the next hour. Therefore, the input parameter of the LSTM sub-model is the surface flow in the last two hours.
where
is the activation function of the nested deep learning model,
is the LSTM transformation,
x is a feature vector of the LSTM model,
is the transpose feature’s weight of the CNN model,
is the activation function of the CNN model, and
b is the bias.
3.2.4. Hyper-Parameter Setting
The hyper-parameters in our Cs-net model include learning rate, epoch, loss function, calculation speed and neural numbers. The learning rates we set are 0.1, 0.01 and 0.001, respectively. The epochs are 100, 500, and 1000. We use mean square error as the loss function. As the Cs-net is a joint deep learning model, following LeCun [
42], we set the neural numbers to be 64, 128 and 256, respectively.
3.2.5. Evaluation of Models
To determine how the model will perform on future data, we need to evaluate the performance of the prediction model. We used a regression model performance evaluation metric to assess our model as it is based on regression algorithms. The evaluation index is as follows:
- (1)
Mean square error (MSE)
- (2)
Mean absolute error (MAE)
- (3)
R-square (
R2)
where
y is the actual output value and
is the prediction value of the model.
3.3. Fusion Model
In this section, we discuss how to fuse the kinetic method and Cs-net model to produce the dynamic ground risk model. In
Section 2, we analyze the ground risk obtained from the people density, casualty rate and probability of a fallen drone.
Hence, the structure of the model is as follows (
Figure 11): First, the input data of the model comprises three parts: people density, aircraft operation data, and map data. According to the different data characteristics, we input these three sets of data into the corresponding models. The people density and map data with spatial-temporal characteristics were input into the joint deep learning model. The aircraft operation data were input into the kinetic model. The two models are parallel due to the different input data and the two models may not affect each other. The Cs-net model is built to predict the people density using a deep learning method. Considering that the shelter factor is a special expression of the map data (the shelter factor can be calculated by the building coverage in the map data), when predicting the population density, the shelter factor is calculated using the map data accordingly. The kinetic model is applied to evaluate the size of the ground impact and the impact energy. As the fatality rate is combined with the shelter factor and the impact energy, we need to fuse the output value from the Cs-net model and the kinetic model.
where
is the probability of a fatality rate,
is the shelter factor,
is the impact energy of ground impact, and
is the kinetic coefficient.
Last, according to the ground risk definition in
Section 2, the dynamic ground risk model can be calculated using the people density output from the Cs-net model, the ground impact output from the kinetic model, the fatality rate and the probability of a fallen drone.