Next Article in Journal
Phase Conductance of BiFeO3 Film
Previous Article in Journal
Enhancing Leaf Area Index Estimation for Maize with Tower-Based Multi-Angular Spectral Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

UAV Cluster Mission Planning Strategy for Area Coverage Tasks

1
College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
2
College of Artificial Intelligence, Xinjiang Vocational and Technical College of Communication, Urumqi 831401, China
*
Author to whom correspondence should be addressed.
Submission received: 17 July 2023 / Revised: 10 September 2023 / Accepted: 19 September 2023 / Published: 11 November 2023
(This article belongs to the Section Vehicular Sensing)

Abstract

:
In the context of area coverage tasks in three-dimensional space, unmanned aerial vehicle (UAV) clusters face challenges such as uneven task assignment, low task efficiency, and high energy consumption. This paper proposes an efficient mission planning strategy for UAV clusters in area coverage tasks. First, the area coverage search task is analyzed, and the coverage scheme of the task area is determined. Based on this, the cluster task area is divided into subareas. Then, for the UAV cluster task allocation problem, a step-by-step solution is proposed. Afterward, an improved fuzzy C-clustering algorithm is used to determine the UAV task area. Furthermore, an optimized particle swarm hybrid ant colony (PSOHAC) algorithm is proposed to plan the UAV cluster task path. Finally, the feasibility and superiority of the proposed scheme and improved algorithm are verified by simulation experiments. The simulation results show that the proposed method achieves full coverage of the task area and efficiently completes the task allocation of the UAV cluster. Compared with related comparison algorithms, the method proposed in this paper can achieve a maximum improvement of 21.9% in balanced energy consumption efficiency for UAV cluster task search planning, and the energy efficiency of the UAV cluster can be improved by up to 7.9%.

1. Introduction

In recent years, with the continuous development of unmanned aerial vehicle (UAV) technology, the application of UAVs in area coverage tasks such as plant protection, target search, and forest fire prevention has gradually become a reality [1,2,3,4,5]. UAVs have the advantages of high flexibility and task efficiency compared to traditional manual or mechanical operations in area coverage tasks [6,7]. Particularly in harsh environments, complex terrains, and areas with limited accessibility, using UAVs for search or agricultural operations can overcome environmental constraints, significantly reduce labor costs, and demonstrate the significant value of UAV applications [8,9]. In the process of using drones to perform tasks, a single UAV cannot meet the requirements of the task. Therefore, it is necessary to introduce multiple drones (referred to as a UAV cluster when small-scale drone operations occur or a UAV swarm when the number of drones increases) to meet the task requirements through task allocation, collaborative search, and situational awareness technologies.
Therefore, many experts and scholars have conducted research on UAV task assignments in area coverage. For example, Ref. [10] proposed an opposition-based learning parameter-adjusting harmony search algorithm to enhance the efficiency of UAV cluster task assignment. Ref. [11] presented a UAV task assignment method based on a simulated annealing strategy, which improved the task balance among UAV cluster individuals to a certain extent. Ref. [12] designed a distributed grouping cooperative dynamic task assignment method based on the extended contract net protocol, considering multiple constraints such as task writing, execution order, and environmental factors among UAV clusters, thereby optimizing the task execution efficiency. Ref. [13] addressed the problem of UAV trajectory optimization in three-dimensional complex environments and proposed a trajectory planning strategy based on an improved Harris Hawks algorithm. Ref. [14] focused on UAV participation in search and rescue scenarios and proposed an optimization strategy based on joint or decoupled optimization, which improved the search efficiency of UAVs.
We have summarized and categorized some existing research, and the results are shown in Table 1.
From Table 1, it can be seen that the aforementioned research mainly focuses on the efficiency of UAV cluster tasks, task load balance among UAV cluster individuals, and UAV trajectory planning. The methods used in these studies mainly involve mathematical modeling and optimization using swarm intelligence algorithms [18].
In addition, the experimental settings in the mentioned research are mostly in two-dimensional planes with small experimental scopes, and some studies involve only a single or a few UAVs. As research progresses and task scenarios expand, the study of UAVs in the field of area coverage tasks is becoming more diverse and complex, involving larger numbers and clusters of UAVs. In this context, the key focus of research is to design appropriate task allocation algorithms and achieve efficient task allocation and high-quality trajectory planning for UAV clusters based on relevant application scenarios [19,20].
Furthermore, the full coverage tasks of UAV clusters are mainly conducted in terrains with varying elevations, such as hills, mountains, and basins. UAV path planning needs to further consider different operating scenarios (altitude, wind speed, temperature, humidity, etc.) and make corresponding decisions for specific flight environments in order to adapt to different application situations [21]. Therefore, it is necessary to consider the impact of terrain variations on task allocation and trajectory planning for UAV clusters. Thus, it is essential to extend the research on efficient task allocation for UAV clusters in three-dimensional space beyond the two-dimensional plane.
As the task area increases and the computational complexity of three-dimensional scenarios grows, there is an urgent need to improve the efficiency of UAV cluster task planning algorithms [22]. In light of this, we first conduct a detailed analysis of area coverage tasks and provide further explanations of the research background. Then, we establish models for task area search, task area division, and UAV cluster cooperative task trajectory planning and derive the energy consumption model for UAV clusters. Next, we use an improved fuzzy C-means clustering method to divide the task area of the UAV cluster. Based on this, we propose a strategy that combines an optimized particle swarm algorithm with an ant colony algorithm for UAV cluster trajectory planning. Finally, the feasibility and effectiveness of the proposed approach are verified through relevant simulation experiments.
The schematic diagram of the proposed approach is illustrated in Figure 1.
The organization of this paper is as follows.
In Section 1, we provide a detailed analysis and summary of the research background and current status. In Section 2, we conduct a comprehensive analysis of area coverage tasks. In Section 3, we perform theoretical analysis and mathematical model derivation for UAV cluster task allocation and UAV cluster trajectory planning. In Section 4, we design relevant algorithms to solve the model optimization problems. In Section 5, we conduct simulation experiments and comparative experiments to evaluate the proposed strategies and discuss the results of these experiments. Finally, we summarize the research conducted in this paper and provide future prospects.

2. Overview of UAV Cluster Area Coverage Tasks

In this paper, area coverage tasks refer to tasks that require full coverage operations in a specific area, such as target search, forest fire inspection, and medication spraying in mountainous vegetation. Figure 2 illustrates the application scenarios.
This paper focuses on the research of rotary-wing UAVs. Compared to fixed-wing UAVs, rotary-wing UAVs have better maneuverability, making it easier to plan and allocate tasks for each UAV in the cluster. This enables faster and more optimal task allocation while ensuring algorithm efficiency and accuracy.
Before performing task allocation, we assume that the map of the task area has been constructed using satellite images or through UAV SLAM (simultaneous localization and mapping). Based on this, we establish a 3D coordinate system to accurately calibrate the relevant area, as shown in Figure 3.

3. Model Establishment

3.1. Coverage of UAV Search Area

In order to perform the area coverage task, UAVs are typically equipped with signal detection devices or actuators required for the task, such as infrared imaging devices, optical cameras, flame detection sensors, and spraying mechanisms. Taking infrared cameras and conventional onboard optical cameras as examples, these detection devices generally employ a cone-based detection method. The UAV serves as the vertex of the cone, and the target detection process is performed in a conical-shaped spatial area. For UAV cluster search tasks in the near-ground area, the UAV cluster carries search sensors, and the search model is illustrated in Figure 4.
Taking infrared area search as an example, the sensor carried by the UAV exhibits a conical beam search on the near-ground surface. The search radius and search coverage resolution vary at different heights.
In the process of executing the search task with an infrared sensor, the search range and search resolution are the main technical performance indicators. Due to the attenuation of infrared radiation in the medium, the search range varies for different search targets and environments. In this study, based on actual conditions, the maximum searchable height for a single UAV carrying a sensor is determined.
As for the search resolution, it should meet the minimum resolution required for imaging the target under the condition of a limited imaging field of view. Considering both imaging resolution and imaging height, the maximum searchable distance satisfying these conditions is determined. The spatial search model for a single UAV is as follows:
h s < h s m a x
S s = ( h s · t a n α s ) 2 · p i
R s = h s · t a n α s 2 · p i 4 · h s 2 · S s
where h s represents the actual height during the UAV search process, h s m a x represents the maximum value of UAV search height, indicating that the UAV cannot perform the search task when it exceeds this height; S s represents the area range that the UAV can project onto the ground from near-ground surface coverage at the current height; α s represents the cone angle during the search process; R s represents the search resolution of the sensor carried by the UAV under the current height condition.
Using the conical area coverage method, the UAV cluster in Figure 5 can achieve full coverage of the task area. The coverage result is shown in Figure 5.
In Figure 5, the color depth represents the elevation of the terrain in the actual environment (darker color indicates lower altitude). During the UAV’s search for ground targets, an optimal solution for the ideal search height needs to be obtained based on the UAV’s current position.

3.2. Division of Task Area for UAV Cluster

In the process of UAV cluster operations, it is necessary to divide the task area into sub-areas. Unlike the uniform division of general task areas, UAVs mainly encounter the following three situations when searching for task points within the task area.
In an actual UAV near-ground search, the spatial relationship between the UAV and the near-ground search area is shown in Figure 6.
As shown in the figure above, α represents the search cone angle of the UAV. When conducting search tasks, the UAV’s coverage area may be on flat terrain. However, in most cases, the terrain is not flat but has certain height variations. For different near-ground search scenarios at different heights, two additional situations can be further distinguished, as illustrated in Figure 7.
In Figure 7, α represents the search cone angle of the UAV. β 1 and β 2 are the angles between the terrain and the horizontal plane in two different scenarios. In the top left image, β 1  < α , indicating that the UAV’s movement is similar to that on flat ground. In the top right image, β 1  > α ; under the consideration of a safe distance for near-ground search, there exists a blind zone in the low-lying terrain where the UAV cannot perform the search. This is illustrated in Figure 8.
In the above figure, the brown-colored area represents the search blind zone for the UAV. In this case, the UAV needs to move directionally along the low-lying terrain to satisfy the search objectives in that area. This is illustrated in Figure 9.
For UAV cluster-based area coverage tasks, a reasonable partitioning of the task area can improve the operational efficiency of the UAV cluster and reduce search energy consumption.
Furthermore, during the process of UAV cluster-based area partitioning, since the comprehensive paths and energy consumption of the UAVs have not been fully obtained, it is necessary to initialize the task partitioning process based on ground reference points. These points should be selected as task points with a certain height that satisfies the UAV’s flight safety distance conditions.
h s m h s m a x
h s m is the height of the task points, and the maximum height is the searchable farthest altitude of the UAV. In this study, the impact of the partitioning results on the workload balance among individual UAVs within the UAV cluster is considered as the evaluation criterion for the UAV cluster partitioning method.
The process of establishing the balance model for UAV search area partitioning in this study is as follows:
h s d i s = h s m a x h s m i n
h s a v = ( h s ) / n s
x s a v = ( x s ) / n s
y s a v = ( y s i ) / n s
where h s m a x represents the maximum height within the current UAV partitioned area; h s m i n represents the minimum height within the current UAV partitioned area; h s d i s represents the maximum height difference within the current UAV partitioned area. h s represents the height of the search task points that the UAV needs to pass through; n s represents the number of task points that the UAV needs to pass through within its assigned area; h s a v represents the average task operation height of the UAV within its assigned area. x s represents the sum of the x-coordinate values of all the points that the UAV needs to pass through; x s a v represents the average x-coordinate value of the points within the spatial range that the UAV needs to pass through. y s represents the sum of the y-coordinate values of all the points that the UAV needs to pass through; y s a v represents the average y-coordinate value of the points within the spatial range that the UAV needs to pass through.
The definition of the workload difference between individual UAVs within the UAV cluster can be expressed as follows:
V s , s + 1 = i = 1 n s ( h s + 1 h s + 1 a v ) i = 1 n s ( h s h s a v ) 2 + i = 1 n s ( x s + 1 x s + 1 a v ) i = 1 n s ( x s x s a v ) 2 + i = 1 n s ( y s + 1 y s + 1 a v ) i = 1 n s ( y s y s a v ) 2 ( D s + 1 D s )
In the above equation, Vs,s+1 represents the difference in task volume between different UAV individuals, Ds is the search traversal distance value for UAV individuals. After defining the workload difference model, the task allocation in the cluster is based on the workload difference among UAVs within the task area.

3.3. UAV Path Planning

The UAV cluster path planning model involves deriving a mathematical model for UAV task path planning based on the determined task areas for individual UAVs.
First, the path motion of a single UAV is illustrated in Figure 10.
In Figure 10, the solid line arrow represents the direction of UAV search movement, while the dashed line arrow represents the actual movement process. During the search process, a single UAV should cover the area in a circular pattern, while the arrangement of the required task points follows a hexagonal pattern. As a result, there will be overlapping search areas during the UAV search process. In order to minimize the overlap, the UAV’s movement direction is illustrated in Figure 11.
In Figure 11, the solid line represents the direction of the UAV’s search movement, while the dashed line represents the actual trajectory. In the search process, a single UAV should ideally cover a circular area, but the arrangement of task points follows a hexagonal pattern, resulting in overlapping search areas. To optimize the search process in the task area, the UAV’s movement direction should be such that it minimizes the repeated coverage area when reaching each task point. However, in practical planning, the focus is on achieving the overall shortest path, and it may not be possible to eliminate the UAV’s actual trajectory. This study introduces an optimization condition that minimizes the repeated coverage area when each task point is reached while considering the global path optimization based on the shortest path algorithm. The specific research process for UAV’s trajectory planning model is described as follows:
Based on the criteria of achieving the fastest full coverage of task points within the area, the definition of search efficiency is as follows:
t s = i = 1 k m | p s p s + 1 | v s
where p s represents the current position of the UAV; p s + 1 represents the position of the next task point that the UAV needs to operate on; k m represents the number of task points within the UAV’s allocated area; v s represents the search speed of the UAV. t s represents the total time required for the UAV to complete the full search of the task points within its allocated area.
When a UAV completes all search tasks within a task area of size A in time t, the search efficiency can be defined as
e f s = S s t s
where S s represents the total area of UAV search, and e f s represents the efficiency value based on the fastest trajectory planning objective.
The search model based on minimizing the repeated area when arriving at new task points can be defined as
S r s = i = 0 k ( S f s S c s )
where S f s is the actual search area of the UAV; S c s is the theoretical coverage area of UAV task points; S r s is the area of repeated search generated by the UAV; and k is the number of task points already completed by the UAV.
Based on minimizing the repeated area when arriving at task points, the efficiency definition can be expressed as
e a s = S r s t s × β a
where e a s represents search efficiency for the lowest repeat area objective; β a   represents a constant coefficient that provides an equivalent calculation for two different search efficiencies. By establishing models for two search efficiencies in the UAV cluster’s trajectory planning process, the search tasks of the UAV cluster can be carried out with a balanced and efficient search for the search objectives and task areas as much as possible.

3.4. Energy Efficiency Model for the UAV Cluster

During the process of conducting area coverage search tasks, energy consumption is an important optimization criterion for the UAV cluster in addition to achieving the objective of complete coverage. Adopting a low-power search approach enables the UAV cluster to have a larger operational range during task execution. Additionally, the reduced power consumption during movement can compensate for the energy usage of the UAV’s onboard sensors, further improving search efficiency.
The energy consumption model is defined in the following process.
First, for a single UAV in the spatial task area, the distance traveled is represented by l , the length of the route, and the time spent on traveling is represented by t . The energy consumption is calculated as follows without changing the direction of movement.
E S l = l s l · E i l + l s v · E i v + t · E b
where E i l represents the energy consumption per unit distance in the horizontal direction for the UAV, E i v represents the energy consumption per unit distance in the vertical direction, l s l represents the position change in the horizontal direction, l s v represents the position change in the horizontal direction, and E b represents the energy consumption during the idle time while the UAV is in motion.
At the UAV’s turning points, we define the energy consumption for all turns as E S d , meaning that each turning event incurs an energy cost of   E S d .
For a single UAV operating within its assigned task area with a known number of turns ksd and the number of task points k , according to the trajectory planning algorithm, the distance the UAV needs to fly is l s . The energy consumption incurred by a single UAV in completing the full coverage search task within the assigned area is given by E S u .
E S u = l s l · E i l + l s v · E i v + t · E b + k s d · E S d
Furthermore, assuming that the total number of task points in the search area required by the UAV cluster is K , and the number of UAVs in the cluster is n , we have
K = k 1 + k 2 + k 3 + + k n
Then, by using the trajectory planning algorithm, we can obtain the flight paths required for UAVs to search within their respective areas, denoted as l 1 , l 2 l n . Therefore, under the current area partitioning method and trajectory planning algorithm, the total power consumption for the UAV cluster to complete the area coverage search in the task area is given by
E S u = i = 1 n l i l · E i l + i = 1 n l i v · E i v + t · n · E b + K · E S d
Finally, based on the energy consumption of the UAV cluster during the search in the target area, we can further establish a model for the energy efficiency ratio. That is, based on the known energy consumption, we can establish the UAV cluster’s energy efficiency model in terms of the time required for the UAV search.
Let the energy consumption of a single UAV to complete the search in its respective area be denoted as E S u i , the duration of the search consumption be denoted as t i . The energy efficiency of a single UAV is given by
e i = E S u i t i
For the UAV cluster, energy efficiency can be divided into balanced energy efficiency and overall energy efficiency. Balanced energy efficiency is used to characterize the differences in energy efficiency among UAVs within the cluster. Overall energy efficiency, on the other hand, compares and verifies the energy efficiency of the UAV cluster as a whole.
In the UAV cluster, the individual energy consumption of each UAV is given by E S u 1 ,   E S u 2 ,   E S u 3 ,   ,   E S u n , and the corresponding time consumed by each UAV is given by t 1 ,   t 2 ,   t 3 ,   ,   t n . The model for balanced energy efficiency is as follows
E S u a = i = 1 n E S u i n
e b = 1 i = 1 n E S u i E S u a t i / a b
where a b is constant for balanced energy efficiency.
The total energy efficiency model is as follows
e a = i = 1 n E S u i i = 1 n t i
By introducing the overall energy efficiency of the UAV cluster and the energy efficiency of individual UAVs, the algorithm performance can be comprehensively evaluated.

4. UAV Cluster Search Task Algorithm Design

4.1. UAV Cluster Task Area Planning

The UAV cluster search task algorithm mainly consists of the search task area division of the UAV cluster and the trajectory planning algorithm design for UAVs within the corresponding areas.
The task area division of the UAV cluster is aimed at achieving the fastest coverage of the specified area under the assumption of the same search efficiency and motion state for multiple UAVs. Traditional clustering methods have difficulty in clearly clustering edge points and rely solely on the number of iterations to determine the completion of clustering tasks [23]. To address this issue, this paper proposes an improved fuzzy C-means clustering method (O-FCM) for the task area division of the UAV cluster.
The convergence process of the cluster centers is as follows:
(1)
In this algorithm, the membership function is modified by incorporating the Euclidean distance between the current traversed node and its adjacent nodes as the relevance weight. This modification enhances the similarity of node classification. The optimized membership function is as follows:
u i j = 1 ( d k = 1 c d i j d i k 2 m 1 i m i n ) σ
where d i m i n represents the Euclidean distance to the nearest node from the current node, and σ represents the inclusion of the influence factor of the nearest distance in determining the weight of the membership function. Regarding the method for selecting the nearest node, instead of calculating the Euclidean distance between the current node and all other nodes in the area, the algorithm sorts all nodes in the area in ascending order based on their positions along the X, Y, and Z axes. It then calculates the Euclidean distance between the current node and the nearest node along each axis (X, Y, and Z) separately. Finally, it selects the minimum value among the three distances as the Euclidean distance between the nearest node and the current node.
(2)
At this point, the calculation of the new centroids within the current category is performed based on the following steps:
c i = i = 1 n u i j m x i i = 1 n u i j m
Then, based on the membership degrees and coordinate values of each node within the area, the optimal coordinates of the current clustering center nodes within the area can be obtained.
(3)
The objective function value under the current condition is calculated as follows:
J z = i = 1 c J i = i = 1 c j n u i j m d i j 2
Further, the convergence condition of the current clustering algorithm is determined based on Equation (25).
J z + 1 J z < η J
when the convergence of the objective function reaches an accuracy less than η J , it can be determined that the current clustering algorithm satisfies the termination condition.
In this study, for the near-ground UAV area in space, the search coverage is a planar coverage method, which can be further equivalently represented as circular coverage of the planar area. The overlapping area around the circular area cannot be directly calculated to obtain the distance between the centers of adjacent circles in practical situations. Therefore, this study uses polygons as a replacement for circles to achieve the minimal overlapping area as much as possible. The relevant calculations are as follows:
η = 180 360 k s k s s = 360
where k s   represents the number of sides of the polygon, and s represents the number of internal angles of the polygon. Since the number of sides of a polygon should be a positive integer, solving the equation simultaneously yields a hexagonal polygon with the maximum number of sides, which corresponds to the minimum overlapping area.
By constructing a honeycomb-like continuous hexagonal coverage area, the search range of the UAV cluster is obtained. The honeycomb pattern formed by hexagons represents the actual movement coverage range for the UAVs. The common center of the circles and hexagons represents the coordinates that the UAVs need to pass through.
Here is the pseudocode for the O-FCM algorithm (See Algorithm 1):
Algorithm 1. O-FCM
1:Start the task.
2:Initialization: Initialize the clustering iteration parameter
3:Randomly generate cluster centers.
4:Repeat:
5: Recalculate the cluster center based on the current category
6: Update membership function based on current clustering
7: Regenerate the cluster based on the current cluster
8:Until:
9: The objective function satisfies the convergence condition
10:Get the cluster center and the categories that each node belongs to
11:End

4.2. UAV Trajectory Planning

The design of the UAV cluster trajectory planning algorithm is mainly based on the collaborative mission objectives of the UAV cluster, aiming to find a comprehensive solution for UAVs’ traveling process under the given environmental background and optimization conditions. To address the drawbacks of poor convergence and susceptibility to local optima in traditional swarm intelligence algorithms, this paper proposes an optimized particle swarm optimization hybrid ant colony algorithm (PSOHAC) to efficiently achieve trajectory planning for UAV clusters.
The ant colony algorithm primarily optimizes the model based on pheromone trails along the ant paths. However, the ant colony algorithm suffers from slow convergence speed [24]. In this study, we incorporate a dynamic pheromone calculation mechanism and a pheromone trail filtering mechanism into the original ant colony algorithm. The dynamic pheromone calculation mechanism uses a dynamic weight accumulation method to accumulate pheromone trails based on the existing pheromone accumulation. The number of accumulations is related to the weight. The pheromone trail filtering mechanism eliminates paths directly when the pheromone intensity on a path is smaller than the mean of all pheromone intensities and a fixed constant during the ant colony iteration. By incorporating these two mechanisms, the convergence speed is improved while ensuring the robustness and distributed nature of the ant colony algorithm.
In this article, the individual trajectory planning of UAVs from the starting point to the destination adopts the ant colony algorithm with dynamic pheromone weights. The algorithm is outlined as follows.
In the algorithm, the ant population is set to a total of ants n a c . At each time step t a c , the path taken by an individual ant is denoted as l i j . After each ant passes through the path, if the path successfully leads to the destination, the intensity of pheromone on that path is
τ i j t k = α a c · n a c + τ i j t k 1 β a c
where τ i j t k represents the concentration of pheromone on the path from i to j at time t k . n a c represents the number of times the current path has been traversed by ants. α a c is a constant factor used in the adaptation of the ant colony algorithm’s pheromone calculation. τ i j t k 1 represents the concentration of pheromone at the previous time step. β a c is a balancing constant coefficient used to ensure the correlation between the number of times and the concentration of pheromone, while avoiding getting trapped in local optima.
The probability expression for the ant at position i to select path j in the algorithm is as follows:
p i j k = τ i j t k · η i j t k s a l l τ i s t k · η i s t k , j a l l 0 , others
where a l l represents the set of all paths that can be reached from i , the current position of the ant. η i j t k represents the visibility of the path from i   to j , which is related to the distance from the selected path to the destination.
Furthermore, the concentration of pheromone on the paths is updated after each generation of ants completes the journey from the starting point to the destination. The update equation is as follows:
τ i j t + k = ρ · τ i j t + τ i j
τ i j t + k = n i = 1 n a c τ i j n i
where ρ represents the persistence of the current pheromone on the path; τ i j t + k is the updating process of pheromone concentration on the path at time k ; n a c represents the total number of ants.
During each iteration of the ant colony, the individual ants optimize the paths and evaluate their quality based on the pheromone concentration. In the next iteration of the ant colony, further iterative calculations are performed on the high-quality paths. After completing one iteration of the ant colony optimization, the existing pheromone concentration in the optimized paths is updated, and the iteration process continues in the next generation. After completing the trajectory optimization process for all individual ants and reaching the maximum iteration limit, the algorithm concludes the optimization process of UAV individual paths. The use of dynamic pheromone concentration improves the optimization efficiency of the ant colony algorithm by reducing the computational dataset under the condition of increased computational complexity.
In the particle swarm optimization (PSO) algorithm, in each iteration, the best individual of the current generation is obtained by calculating the fitness value of all individuals. However, within each generation, the individuals have a certain degree of homogeneity in terms of fitness value and can be classified into different ranges of fitness values [25]. In this paper, the characteristics of the niche population mechanism are utilized to categorize the particle swarm population individuals in each generation [26]. Different dynamic weight proportions are assigned to the fitness values of different populations in each iteration, enhancing the comparability of the iteration of the internal groups within the particle swarm population. The specific content of the optimized particle swarm optimization algorithm is as follows.
(1)
Based on the UAV cluster collaborative task model, the updating equation of the niche particle swarm algorithm is expressed as follows:
V i n = ω V i n + c 1 r 1 ( p i n x i n ) + n m u c 2 r 2 ( p g n x i n ) x i n = x i n + V i n
where V i n represents the real-time velocity of a particle in the niche-optimized particle swarm algorithm, x i n represents the current spatial position of the particle individual. p i n represents the global best position, p g n   is the best position among all positions traversed by all particle individuals. c 1 ,   r 1 , c 2 ,   r 2 are weight constants. h is the segmentation criterion for the niche population.
(2)
After establishing the state update process for particle individuals, it is necessary to define the objective functions for particle individual optimization, niche population optimization, and the overall algorithm objective achieved by the final iteration of the particle swarm population. The objective function for particle individual optimization is defined as follows:
( F p i = C 1 + f p i ) | i = 1 n f p i K
where F p i   represents fitness function values representing individual particles; K represents the shortest spatial path achievable by a UAV cluster regardless of obstacles; C 1 represents the weighted weights of the small habitat population to which the particle fitness value belongs; f p i   represents the current fitness value of the particle in the particle population.
(3)
After achieving the optimization objectives and requirements for individual particles within the particle population, the process of dividing the particle population into small habitat populations is performed. In the small habitat particle swarm algorithm, the division rule for the small habitat population is determined based on the similarity of fitness value distribution in each generation of the particle population, as shown in the following equation:
n s u ,       F p i > F p i k i · m a x F p i n n m u ,     F p i k i · m i n F p i n F p i F p i k i · m a x F p i n n l u ,     F p i n < F p i k i · m i n F p i n
where n is the number of individuals in the current particle population, i is the specific iteration count, k is the total number of iterations, n s u is the set of excellent individuals, n m u is the set of balanced individuals, and n l u is the set of poor individuals. In the initial iterations, the screening degree of individual excellence is higher, and the number of individuals screened out is larger. As the iterations progress, the number of screened individuals decreases. Towards the end of the iteration, individuals are no longer eliminated.
(4)
After determining the optimization objectives for individual particles within the particle population and the small habitat populations, the iterative optimization process begins. It involves the dynamic data link with the ant colony algorithm, enabling the two algorithms to collaborate in the optimization process. The small habitat particle swarm algorithm ultimately determines the final optimization result based on the number of iterations. Once the specified number of iterations is reached, the algorithm terminates and outputs the computed result.
The pseudocode for the PSOHAC algorithm is as follows (see Algorithm 2):
Algorithm 2. PSOHAC
1:Start the task.
2:Initialization: Initialize individual particle data, particle swarm fitness function, niche population parameters, particle population data and ant population parameters.
3:Repeat-1:
4:    Updated the ant sport count.
5:    Calculate the state transition probability.
6:    Computational analysis of pheromone.
7:    Completion meets the iteration limit.
8:    Update the particle population properties
9:    Repeat-2:
10:        Fitness function algorithm optimization.
11:        Repeat-3:
12:            Update individual fitness of particles.
13:        Until-3:
14:            Reach the limit of quantity.
15:    Until-2:
16:        Meet the iteration limit
17:        Iterate the optimal cluster path.
18:        Update pheromone properties
19:Until-1:
20:    Iterative optimal path.
21:    Obtain optimal UAV cluster path planning.
22:End

5. Simulation Analysis

In this study, an optimized fuzzy C-means clustering algorithm and an optimized particle swarm algorithm combined with the ant colony algorithm are used to achieve task allocation for full coverage search of UAV clusters in a given area. The specific simulation parameter settings and experimental process are as follows.

5.1. Algorithm Parameter Settings

The parameters used in the algorithm are set as Table 2.

5.2. Simulation Experiment Testing

The simulation in this study is conducted in a horizontal range of 1 km × 1 km and a height range of 0 to 100 m. A 3D mountain model is generated within this range to represent the terrain. The specific illustration is presented in Section 2, Figure 3b.
First, based on the task points within the environment, simulation tests are conducted to evaluate the number of UAVs required for a complete search of the designated area.
In the task area, a higher number of UAVs will result in faster search times. However, in terms of energy consumption, the relationship between the number of UAVs and the energy consumed during the search is not straightforward. In this study, simulations are conducted to analyze the energy efficiency generated by different numbers of UAVs during the search operation. Multiple experiments are conducted under consistent conditions in terms of other parameters and derivation processes, and the average values are obtained. The test results are presented in the Figure 12.
According to Figure 12, the optimal energy consumption efficiency is achieved when the number of UAVs is eight.
In the aforementioned experimental environment, the improved fuzzy C-means clustering algorithm described in this study was applied. The spatial clustering results obtained in the mountainous planning area are depicted in Figure 13.
In the above figure, differently colored points represent the task points assigned to different UAVs. From the figure, it can be observed that the method described in this study has successfully accomplished the task allocation for the UAV cluster within the designated area. At this stage, the task area range corresponding to each UAV has been preliminarily identified. The next step is to further plan the UAV trajectories.
Using the PSOHAC algorithm described in this study, the task trajectory planning for the UAV cluster corresponding to the task area is performed. The results are shown in Figure 14.
From Figure 14, it can be observed that the PSOHAC algorithm described in this study has successfully performed trajectory planning for the respective task points of the eight UAVs. Each trajectory is smooth and without any intersections.
Furthermore, the fitness values for each of the eight UAVs are obtained, as shown in Figure 15. In the graph, the maximum difference in fitness values between UAVs is 14%, with an average difference of 9%. These results demonstrate the feasibility and effectiveness of the proposed approach described in this study.

5.3. Comparative Experimental Testing

To demonstrate the superiority of the improved algorithms proposed in this study, we conducted comparative experimental testing against several classical traditional algorithms under consistent parameter settings. In terms of task area partitioning, we compared the K-means clustering algorithm, the C-means clustering algorithm, and the fuzzy C-means clustering algorithm used in this study. For UAV cluster trajectory planning, we conducted comparative testing against the immune algorithm (GA) [27], ant colony optimization (ACO), particle swarm optimization (PSO) [28], and the simulated annealing (SA) algorithm [29].
First, to verify the performance of the PSOHAC algorithm, we conducted comparative experiments using the UAV trajectory planning model described in Section 3.3 as the test function. The relevant results are shown in Figure 16 and Figure 17.
Based on the comprehensive analysis of Figure 16 and Figure 17, it can be observed that the proposed PSOHAC algorithm has shown advantages in terms of convergence speed and optimization accuracy compared to the other three comparative algorithms. The dynamic pheromone calculation mechanism introduced in the ant colony algorithm, as well as the iterative process based on the small habitat theory in the particle swarm algorithm, have contributed to a 22% faster convergence speed and a 20% higher optimization accuracy.
Building upon these findings, consistent with the testing experiments in Section 5.2, further comparative experiments were conducted. The results of the comparative experimental testing are presented in the following Table 3.
From Table 3, it can be observed that our proposed method exhibits higher balanced energy consumption efficiency compared to other comparative algorithms. In other words, in the process of UAV cluster task allocation, our method achieves a more balanced distribution of tasks among individual UAVs. In terms of task allocation, our method results in a lower total task quantity compared to other comparative algorithms. Furthermore, in terms of task efficiency, our method effectively reduces the length of UAV trajectories while ensuring full coverage of the task area.

6. Conclusions

This paper addresses the issues of high-energy consumption, uneven task allocation, and low task efficiency in unmanned aerial vehicle (UAV) cluster missions for area coverage. It proposes approaches from the perspectives of task area division and UAV trajectory planning.
Pre-processing of the task area is conducted, and cluster task planning models, UAV trajectory planning models, and UAV energy efficiency models are established. For the task area division problem in UAV cluster missions, an improved fuzzy C-clustering method is proposed. It incorporates neighborhood attributes in addition to node clustering attributes to assist in accurate node clustering.
For UAV trajectory planning, a strategy combining the optimization of particle swarm optimization with ant colony optimization is proposed; the strategy focuses on multi-objective optimization around minimizing repeated area and achieving faster global planning. The strategy aims to overcome the weaknesses of classical swarm intelligence algorithms, such as poor convergence and getting trapped in local optima.
The simulation experiment results demonstrate that the proposed strategies successfully accomplish task allocation and trajectory planning for UAV clusters within the designated task area. Comparative experiments show that the proposed algorithms exhibit improvements in optimization accuracy. Additionally, the trajectory planning strategy achieves advantages in terms of convergence speed.
Compared to other methods, the proposed approach can achieve a maximum improvement of 21.9% in balanced energy consumption efficiency for UAV cluster task search planning. The overall energy efficiency of the UAV cluster can be improved by up to 7.9%. The application of cluster task allocation and path planning presented in this paper has valuable reference value. It is worth noting that the simulations may have some deviations from real-world scenarios, and further adjustments to the parameters will be made through real-world testing in future research.

Author Contributions

Conceptualization, X.Y. and R.C.; methodology, X.Y.; software, X.Y.; validation, R.C. and Z.J.; formal analysis, R.C.; investigation, Z.J.; resources, C.R; data curation, Z.J.; writing—original draft preparation, X.Y.; writing—review and editing, R.C.; visualization, Z.J; supervision, R.C.; project administration, R.C.; funding acquisition, R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The raw/processed data required to reproduce these findings cannot be shared at this time as these data also form part of an ongoing study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, C.; Li, Y.; Cao, G.; Zhang, J. Research on Dynamic Scheduling Model of Plant Protection UAV Based on Levy Simulated Annealing Algorithm. Sustainability 2023, 15, 1772. [Google Scholar] [CrossRef]
  2. Cheng, Q.; Wang, H.; Zhu, B.; Shi, Y.; Xie, B. A Real-Time UAV Target Detection Algorithm Based on Edge Computing. Drones 2023, 7, 95. [Google Scholar] [CrossRef]
  3. Han, G.; Zhao, L.; Li, Q.; Li, S.; Wang, R.; Yuan, Q.; He, M.; Yang, S.; Qin, L. A Lightweight Algorithm for Insulator Target Detection and Defect Identification. Sensors 2023, 23, 1216. [Google Scholar] [CrossRef]
  4. Huang, J.; He, Z.; Guan, Y.; Zhang, H. Real-Time Forest Fire Detection by Ensemble Lightweight YOLOX-L and Defogging Method. Sensors 2023, 23, 1894. [Google Scholar] [CrossRef] [PubMed]
  5. Partheepan, S.; Sanati, F.; Hassan, J. Autonomous Unmanned Aerial Vehicles in Bushfire Management: Challenges and Opportunities. Drones 2023, 7, 47. [Google Scholar] [CrossRef]
  6. Ozkan, O.; Kilic, S. UAV routing by simulation-based optimization approaches for forest fire risk mitigation. Ann. Oper. Res. 2023, 320, 937–973. [Google Scholar] [CrossRef]
  7. Mukhamediev, R.I.; Symagulov, A.; Kuchin, Y.; Zaitseva, E.; Bekbotayeva, A.; Yakunin, K.; Assanov, I.; Levashenko, V. Review of some applications of unmanned aerial vehicles technology in the resource-rich country. Appl. Sci. 2021, 11, 10171. [Google Scholar] [CrossRef]
  8. Liu, K.; Zheng, J. UAV trajectory optimization for time-constrained data collection in UAV-enabled environmental monitoring systems. IEEE Internet Things J. 2022, 9, 24300–24314. [Google Scholar] [CrossRef]
  9. Mohsan, S.A.H.; Khan, M.A.; Noor, F.; Ullah, I.; Alsharif, M.H. Towards the Unmanned Aerial Vehicles (UAVs): A Comprehensive Review. Drones 2022, 6, 147. [Google Scholar] [CrossRef]
  10. Cui, Y.; Dong, W.; Hu, D.; Liu, H. The Application of Improved Harmony Search Algorithm to Multi-UAV Task Assignment. Electronics 2022, 11, 1171. [Google Scholar] [CrossRef]
  11. Huo, L.; Zhu, J.; Wu, G.; Li, Z. A Novel Simulated Annealing Based Strategy for Balanced UAV Task Assignment and Path Planning. Sensors 2020, 20, 4769. [Google Scholar] [CrossRef] [PubMed]
  12. Qin, B.; Zhang, D.; Tang, S.; Wang, M. Distributed Grouping Cooperative Dynamic Task Assignment Method of UAV Swarm. Appl. Sci. 2022, 12, 2865. [Google Scholar] [CrossRef]
  13. Zhang, R.; Li, S.; Ding, Y.; Qin, X.; Xia, Q. UAV Path Planning Algorithm Based on Improved Harris Hawks Optimization. Sensors 2022, 22, 5232. [Google Scholar] [CrossRef] [PubMed]
  14. Yanmaz, E. Joint or decoupled optimization: Multi-UAV path planning for search and rescue. Ad Hoc Netw. 2023, 138, 103018. [Google Scholar] [CrossRef]
  15. Yao, J.; Sha, Y.; Chen, Y.; Zhang, G.; Hu, X.; Bai, G.; Liu, J. IHSSAO: An Improved Hybrid Salp Swarm Algorithm and Aquila Optimizer for UAV Path Planning in Complex Terrain. Appl. Sci. 2022, 12, 5634. [Google Scholar] [CrossRef]
  16. Chowdhury, A.; De, D. RGSO-UAV: Reverse Glowworm Swarm Optimization inspired UAV path-planning in a 3D dynamic environment. Ad Hoc Netw. 2023, 140, 103068. [Google Scholar] [CrossRef]
  17. Hu, C.; Qu, G.; Zhang, Y. Pigeon-inspired fuzzy multi-objective task allocation of unmanned aerial vehicles for multi-target tracking. Appl. Soft Comput. 2022, 126, 109310. [Google Scholar] [CrossRef]
  18. Li, Q.; Xiong, H.; Ding, Y.; Song, J.; Liu, J.; Chen, Y. A Review of Unmanned Aerial Vehicle Swarm Task Assignment. In International Conference on Guidance, Navigation and Control; Springer: Singapore, 2023; pp. 6469–6479. [Google Scholar]
  19. Poudel, S.; Moh, S. Task assignment algorithms for unmanned aerial vehicle networks: A comprehensive survey. Veh. Commun. 2022, 35, 100469. [Google Scholar] [CrossRef]
  20. Sharma, A.; Shoval, S.; Sharma, A.; Pandey, J.K. Path Planning for Multiple Targets Interception by the Swarm of UAVs based on Swarm Intelligence Algorithms: A Review. IETE Tech. Rev. 2022, 39, 675–697. [Google Scholar] [CrossRef]
  21. Li, Y.; Liu, M.; Jiang, D. Application of Unmanned Aerial Vehicles in Logistics: A Literature Review. Sustainability 2022, 14, 14473. [Google Scholar] [CrossRef]
  22. Tang, J.; Duan, H.; Lao, S. Swarm intelligence algorithms for multiple unmanned aerial vehicles collaboration: A comprehensive review. Artif. Intell. Rev. 2023, 56, 4295–4327. [Google Scholar] [CrossRef]
  23. Askari, S. Fuzzy C-Means clustering algorithm for data with unequal cluster sizes and contaminated with noise and outliers: Review and development. Expert Syst. Appl. 2021, 165, 113856. [Google Scholar] [CrossRef]
  24. Wang, Y.; Han, Z. Ant colony optimization for traveling salesman problem based on parameters optimization. Appl. Soft Comput. 2021, 107, 107439. [Google Scholar] [CrossRef]
  25. Ab. Aziz, N.A.; Ibrahim, Z.; Mubin, M.; Nawawi, S.W.; Mohamad, M.S. Improving particle swarm optimization via adaptive switching asynchronous–synchronous update. Appl. Soft Comput. 2018, 72, 298–311. [Google Scholar] [CrossRef]
  26. Trappes, R.; Nematipour, B.; Kaiser, M.I.; Krohs, U.; van Benthem, K.J.; Ernst, U.R.; Gadau, J.; Korsten, P.; Kurtz, J.; Schielzeth, H.; et al. How Individualized Niches Arise: Defining Mechanisms of Niche Construction, Niche Choice, and Niche Conformance. BioScience 2022, 72, 538–548. [Google Scholar] [CrossRef]
  27. Jiacheng, L.; Lei, L. A Hybrid Genetic Algorithm Based on Information Entropy and Game Theory. IEEE Access 2020, 8, 36602–36611. [Google Scholar] [CrossRef]
  28. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 1944, pp. 1942–1948. [Google Scholar]
  29. Shao, W.; Guo, G. Multiple-Try Simulated Annealing Algorithm for Global Optimization. Math. Probl. Eng. 2018, 2018, 9248318. [Google Scholar] [CrossRef]
Figure 1. Schematic Diagram of the Proposed Approach.
Figure 1. Schematic Diagram of the Proposed Approach.
Sensors 23 09122 g001
Figure 2. Schematic Diagram of Area Coverage Task Scenarios.
Figure 2. Schematic Diagram of Area Coverage Task Scenarios.
Sensors 23 09122 g002
Figure 3. Illustration of Task Area Transformation and Calibration. (a) Actual Topographic Map; (b) Transformed Task Map.
Figure 3. Illustration of Task Area Transformation and Calibration. (a) Actual Topographic Map; (b) Transformed Task Map.
Sensors 23 09122 g003
Figure 4. UAV Search Coverage Model.
Figure 4. UAV Search Coverage Model.
Sensors 23 09122 g004
Figure 5. Task Area Coverage with UAV Cone Beam. (a) 3D view; (b) 2D view.
Figure 5. Task Area Coverage with UAV Cone Beam. (a) 3D view; (b) 2D view.
Sensors 23 09122 g005
Figure 6. UAV Search Coverage Area on Flat Ground.
Figure 6. UAV Search Coverage Area on Flat Ground.
Sensors 23 09122 g006
Figure 7. UAV Search Coverage Area on Sloped Terrain.
Figure 7. UAV Search Coverage Area on Sloped Terrain.
Sensors 23 09122 g007
Figure 8. Search Blind Zone on Steep Slope with Incline Angle Greater than UAV Search Cone Angle.
Figure 8. Search Blind Zone on Steep Slope with Incline Angle Greater than UAV Search Cone Angle.
Sensors 23 09122 g008
Figure 9. Directional Search of UAV on Steep Slope.
Figure 9. Directional Search of UAV on Steep Slope.
Sensors 23 09122 g009
Figure 10. Path Motion Process of a Single UAV-1.
Figure 10. Path Motion Process of a Single UAV-1.
Sensors 23 09122 g010
Figure 11. Path Motion Process of a Single UAV-2.
Figure 11. Path Motion Process of a Single UAV-2.
Sensors 23 09122 g011
Figure 12. Relationship Between the Number of UAVs and Energy Consumption Efficiency.
Figure 12. Relationship Between the Number of UAVs and Energy Consumption Efficiency.
Sensors 23 09122 g012
Figure 13. Clustering Algorithm Results Graph. (a) 3D view; (b) 2D view.
Figure 13. Clustering Algorithm Results Graph. (a) 3D view; (b) 2D view.
Sensors 23 09122 g013
Figure 14. UAV Cluster Task Search Results. (a) 3D view; (b) 2D view.
Figure 14. UAV Cluster Task Search Results. (a) 3D view; (b) 2D view.
Sensors 23 09122 g014
Figure 15. Fitness Values for Each UAV.
Figure 15. Fitness Values for Each UAV.
Sensors 23 09122 g015
Figure 16. Three-Dimensional Trajectory Planning Results under Different Algorithms (Replacement). (a) PSOHAC; (b) PSO; (c) ACO; (d) GA.
Figure 16. Three-Dimensional Trajectory Planning Results under Different Algorithms (Replacement). (a) PSOHAC; (b) PSO; (c) ACO; (d) GA.
Sensors 23 09122 g016
Figure 17. Convergence Curves of Different Algorithms.
Figure 17. Convergence Curves of Different Algorithms.
Sensors 23 09122 g017
Table 1. Summary of Existing Research.
Table 1. Summary of Existing Research.
RefResearch ObjectivesMethod/AlgorithmNumber of UAVApplication Scenario
[10]Improve the efficiency of task allocationmulti-objective optimizationCluster2D
[11]Balanced UAV Task AssignmentSimulated Annealing-Based StrategyCluster3D
[12]Improve the efficiency of task allocationModelingCluster2D
[13]Path optimizationImproved Harris Hawks Optimization algorithmSingle3D
[14]Path optimization
Improve search efficiency
ModelingCluster2D
[15]Path optimizationHybrid Salp Swarm Algorithm
Aquila Optimizer
Single3D
[16]Path optimizationReverse Glowworm Swarm OptimizationSingle3D
[17]Balanced UAV Task Assignment
Improve the efficiency of task allocation
pigeon-inspired fuzzy multi-objective optimization algorithmCluster3D
Table 2. Constant Parameter Values and Descriptions in this Study.
Table 2. Constant Parameter Values and Descriptions in this Study.
ParameterDescriptionsValue
α s UAV Cone Beam Search Angle30°
d s f m a x UAV Maximum Safety Distance5 m
h s m a x UAV Maximum Search Altitude5 m
α d Energy-to-Distance Ratio Constant0.8
v i UAV Constant Travel Speed2.1 m/s
t i n Algorithm Initialization Time5.5
E i l Horizontal Unit Distance Energy Consumption0.2
E i v Vertical Unit Distance Energy Consumption0.5
Table 3. Comparative Experimental Results.
Table 3. Comparative Experimental Results.
Trajectory Planning MethodArea Division MethodTotal Path Length/mUAV Flight Time/sBalanced Energy Consumption EfficiencyOverall Energy Consumption Efficiency
Method used in this paperMethod used in this paper1348.678.20.71215.82
K-means clustering1346.578.90.64417.07
C-means clustering1370.481.50.6216.81
GAMethod used in this paper1343.579.90.60416.81
K-means clustering1351.280.40.58416.89
C-means clustering1380.682.10.59216.82
PSOMethod used in this paper1345.380.10.63616.90
K-means clustering1329.079.10.59616.75
C-means clustering1379.882.130.594816.86
ACOMethod used in this paper1345.280.10.63616.79
K-means clustering1341.079.80.60816.81
C-means clustering1355.380.70.61216.77
SAMethod used in this paper1342.279.90.616.53
K-means clustering1333.379.30.63216.85
C-means clustering1369.681.50.58816.84
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, X.; Chen, R.; Jiang, Z. UAV Cluster Mission Planning Strategy for Area Coverage Tasks. Sensors 2023, 23, 9122. https://0-doi-org.brum.beds.ac.uk/10.3390/s23229122

AMA Style

Yan X, Chen R, Jiang Z. UAV Cluster Mission Planning Strategy for Area Coverage Tasks. Sensors. 2023; 23(22):9122. https://0-doi-org.brum.beds.ac.uk/10.3390/s23229122

Chicago/Turabian Style

Yan, Xiaohong, Renwen Chen, and Zihao Jiang. 2023. "UAV Cluster Mission Planning Strategy for Area Coverage Tasks" Sensors 23, no. 22: 9122. https://0-doi-org.brum.beds.ac.uk/10.3390/s23229122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop