Next Article in Journal
Axle Temperature Monitoring and Neural Network Prediction Analysis for High-Speed Train under Operation
Previous Article in Journal
Evaluation of the Quality Supervision System for Construction Projects in China Considering the Quality Behavior Risk Transmission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cost-Effective Placement of Recharging Stations in Drone Path Planning for Surveillance Missions on Large Farms

by
Jean Louis Ebongue Kedieng Fendji
1,*,
Israel Kolaigue Bayaola
1,
Christopher Thron
2,
Marie Danielle Fendji
3 and
Anna Förster
4,*
1
Department of Computer Engineering, University Institute of Technology, University of Ngaoundéré, P.O. Box 454 Ngaoundéré, Cameroon
2
Department of Science and Mathematics, Texas A&M University Central Texas, Killeen, TX 76549, USA
3
Department of Electrical and Electronics Engineering, Faculty of Engineering and Technology, University of Buea, P.O. Box 63 Buea, Cameroon
4
Sustainable Communication Networks, University of Bremen, 8359 Bremen, Germany
*
Authors to whom correspondence should be addressed.
Submission received: 21 August 2020 / Revised: 28 September 2020 / Accepted: 9 October 2020 / Published: 12 October 2020
(This article belongs to the Section Computer)

Abstract

:
The energy limitation remains one of the biggest constraints in drone path planning, since it prevents drones from performing long surveillance missions. To assist drones in such missions, recharging stations have recently been introduced. They are platforms where the drone can autonomously land to recharge its battery before continuing the mission. However, the cost of those platforms remains a significant obstacle to their adoption. Consequently, it is important to reduce their number while planning the path of the drone. This work introduces the Single Drone Multiple Recharging Stations on Large Farm problem (SD-MRS-LF). A large farm is considered as an area of interest to cover with a set of candidate locations where recharging stations can be installed. The aim is to determine the path of the drone that minimizes the number of locations for recharging stations as well as the completion time of the surveillance mission. This path planning problem falls within the realm of computational geometry and is related to similar problems that are encountered in the field of robotics. The problem is complicated due to environmental constraints on farms such as wind speed and direction, which produce asymmetries in the optimal solution. A back-and-forth-k-opt simulated annealing (BFKSA) approach is proposed to solve the defined problem. The new approach is compared to the basic back-and-forth (BF) and a K-opt variant of the well-known simulated annealing (KSA) approach over a set of 20 random topologies in different environmental conditions. The results from computational experiments show that the BFKSA approach outperforms the others, in terms of providing feasible solutions and minimizing the number of recharges.

1. Introduction

One of the latest developments in agriculture is the use of drones. According to a recent PwC report [1], the agricultural drone market is estimated at USD 32.4 billion, taking the second largest share of the market. In, particular, drones are increasingly used in the surveillance of farms [2]. In this application, drones fly over an area of interest to take pictures that can be analyzed later for different purposes including crop yield estimation, early warning of pests and disease, disaster risk reduction, and so on. Drone can be operated in two modes: manual mode requiring a remote human intervention and automatic mode in which the path of drone is determined beforehand [3]. This last mode requires prior planning of the drone path, commonly known as coverage path planning (CPP) [4]. More precisely, CPP for drones may be defined as the determination of an optimal path allowing a drone to complete a coverage mission, while minimizing the mission’s completion time. The area of interest in agriculture is typically a farm that may be represented as a sequence of vertices that form a planar polygon. This polygon generally is decomposed into smaller cells to facilitate the coverage task. The decomposition is essentially based on the ground sampling distance (GSD), which represents the distance between pixel centers measured on the ground [5].
The major constraint in drone CPP is the small capacity of the drone’s battery, which limits the drone’s available energy and thus its flight range. For this reason, much research in the area of CPP has been focused on the development of energy-aware algorithms. One of the main ideas behind common algorithms is the assumption that a drone spends a lot of time and energy making turns [6,7]. Thus, those algorithms modify conventional trajectories such as back-and-forth [8] and HILBERT [9]. An interesting survey on coverage path planning with drones can be found in [10].
Regardless of the efficiency of algorithms, large-scale surveillance missions will require more than a single trip. For instance, the authors in [11] conducted around 200 flights to survey an area of about 29.5 ha with a battery-powered DJI Inspire 1 quadcopter. In an area of four square kilometers as it is the case in [12], a DJI Phantom 3 Professional drone requires at least 10 flights for a mission. Consequently, recharging stations have recently been introduced to assist the drone in continuous operation [13,14]. When the drone runs out of energy, it may land on a recharging station to recharge its battery. Even if it is possible for a drone to land on a moving platform [15] and perform a dynamic wireless charging [16], the movement of such a platform in many farms is limited due the lack of proper roads throughout the farm. Most of the current literature considers predefined locations for recharging stations. For example, authors in [17] use only one drone with some recharging stations located at predefined sites. They tried to minimize the mission completion time by deciding when the drone should stop for charging depending on the state of charge (SoC). Reference [18] takes an alternative approach, by considering a fleet of drones with only one recharging station located at the center of the area of interest—the aim is to cover the area of interest with a constraint on energy replenishment scheduling. Since only one recharging station is used, the surveyed area is limited to a radius of half the traveling distance of a drone, so that the drone can make a round trip to the recharging station.
To automatically survey a larger area will require at least one recharging station located at an acceptable position. The first attempt to solve this issue of determining the optimal location of recharging stations is found in [19]. Authors proposed a Drone Delivery Recharging Location Model (DDRLM) that optimizes the location of a set of recharging stations with the aim of maximizing the customer demand that can be served. In this formulation, the number of recharging stations is given beforehand.
Besides finding optimal locations for recharging stations, their number should also be minimized to reduce the cost of the overall architecture. Actually, the cost of an additional recharging station can significantly influence the cost of the overall architecture. For example, the Heisha recharging station costs about 1000 USD, while a Mavic Air drone that can be recharged using this station costs barely half the price (600 USD) [20]. This shows the necessity to reduce the number of recharging stations, especially for users with limited budgets such as low-income users.
This work introduces the Single Drone Multiple Recharging Stations in Large Farm problem (SD-MRS-LF). Unlike the work in [19] which focuses on the determination of the locations for a given number of recharging stations, this work rather tackles drone path planning with a priority on the economical point of view. In this formulation, an area of interest is given as well as a set of candidate locations within the area of interest where recharging stations may potentially be installed. The objective is two-fold: first to minimize the number of actually used recharging stations and secondly to minimize the completion time of the surveillance mission. The solution with the smallest number of recharging stations will be the best regardless the completion time, since the priority is to reduce the cost. In case more than one solution provide the smallest number of recharging stations, the one with the smallest completion time will be considered as the best. To solve this problem, a new approach has been proposed: back-and-forth-k-opt simulated annealing (BFKSA). BFKSA is a combination of the two well-known algorithms, back and forth (BF) and the K-opt variant of simulated annealing (SA). The idea behind BFKSA is based on the assumption that the quality of the initial solution provided to a SA algorithm can influence the final solution.
The rest of the paper is organized as follows: Section 2 presents the system modelling and the formulation of the Single Drone Multiple Recharging Stations in Large Farm problem. Section 3 presents and explains the Back-and-Forth-k-opt Simulated Annealing approach defined to solve the problem. Section 4 presents the simulation setup and compares the results with basic Back-and-Forth and Simulated Annealing approaches. This paper ends with a conclusion and a look at future work.

2. System Modelling and Problem Formulation

2.1. Basic Assumptions

The following realistic assumptions and approximations have been made to specify the problem:
  • The drone begins the mission fully charged, but a single charge does not supply enough energy for the entire mission.
  • The flying status of the drone is one of two possibilities: vertical moving (take-off/landing), or horizontal motion [21].
  • Take-off energy and landing energy is negligible compared to the energy required for flight.
  • The drone flies at a constant altitude throughout the mission (except when recharging).
  • The area to cover presents no obstacles to impede the drone’s flight. In practice, drones may fly at altitudes at 120 m which is higher than the tallest trees.
  • Weather conditions (including temperature, air pressure, relative humidity, and wind speed) are constant over the entire region, during the entire mission. Note that weather conditions such as wind and temperature can affect the energy consumption. The ambient temperature may influence the battery drain and capacity [22] while the wind may affect either positively or negatively the forward movement of the drone [23], depending on the wind speed and direction, as well as the flying status of the drone.
  • Recharging stations are located within the area to enable the drone to complete its mission. There are a number of pre-identified candidate locations within the region where charging stations may potentially be located.
  • The drone can only visit a recharging station once per mission. This assumption reflects the fact that once a recharge is performed, the station will require a period of time to recharge itself (possibly using solar energy).

2.2. Scenario Modeling

Farm areas typically can be described as planar polygons. Unlike urban areas that may contain nonflying zones or tall buildings, the aerial space over a farm is generally free of obstacles. In this scenario, an area to cover is considered as a rectangular form of size S m2 and decomposed into cells as illustrated in Figure 1a. The set of the cells is represented by A and contains n   × m elements among which some cells are candidate locations where recharging stations may conveniently be deployed. The set of candidate locations is represented by C , with C A .
The decomposition of our target area into a grid depends on the maximum image collection distance interval of the drone. This distance interval depends on the specification of the drone (the field of view angle and the resolution of the camera) and the overlap requirement during image collection. The maximal value of the distance I m a x separating image collection points (as shown in Figure 1b) is computed using (1) as proposed in [24].
I m a x   =   ( 1 r ) 2 H t a n ( θ 2 )
where   θ is the opening angle of the camera; H is the drone’s altitude, and r the overlap ratio.
Figure 2 shows the geometry used to derive (1), with I the size of a cell of the grid.
In real scenarios in farm surveillance, the minimal resolution of each collected image is an important aspect for image processing in further stages. This minimal resolution depends on the purpose of the surveillance mission. Given an image resolution ( φ x , φ y ) , the spatial resolution ρ obtained by taking a picture at an altitude H with a field of view angle θ can be computed by Equation (2).
ρ = φ x 2 H t a n ( θ 2 )
If we consider the minimal resolution ρ m required by the mission, we obtain the inequality
H φ x 2 ρ m t a n ( θ 2 )
From (3) we deduce the maximum altitude H M the drone can fly.
H M = φ x 2 ρ m t a n ( θ 2 )
The maximum altitude with respect to the image resolution requirement is generally used to minimize the time of the mission. During the remote sensing mission for monitoring sorghum growth, the authors in [25] used a ReadyMadeRC Anaconda drone equipped with a camera of 1.2 megapixels and flying at an altitude of 120 m which provided a spatial resolution of approximately 6.5 cm. In addition, the maximum altitude of the drone can be legally restricted.
Due to budget limitation only one drone is considered with a maximal capacity energy E M . Each time the drone stops at a recharging station, it recharges its battery to E M before leaving. The mission starts at a given take-off point S t and ends at a given landing point that maybe different from the take-off point.

2.3. Energy Consumption Model

Under horizontal motion, drones power expenditure may be broken into two components: one component for lift, and the second to overcome the parasitic drag that hinders its forward movement in free space [26,27]. The total power in horizontal moving is given by (5). Details on how the formula is derived are presented in [28].
P = P P + P L
with
P P = 1 2 C D F a D a v 3
P L = W 2 D a b 2 v
where
  • P   is the required total power (in watts );
  • P P is the required power to overcome the parasitic drag (in watts );
  • P L is the required power to lift the drone (in watts );
  • C D is the aerodynamic drag coefficient;
  • F a is the front facing area (in m 2 );
  • W is the total weight of the drone (in kg );
  • D a is the density of the air (in kg / m 3 );
  • b is the width of the drone (in m );
  • v is the speed of the drone relative to the wind (in m/s).
Since P P and P L are respectively proportional and inversely proportional to the speed of the drone, a tradeoff (optimum speed) that minimizes the power consumption of the drone during horizontal moving should therefore be found. To compute this optimum speed, we use (8) proposed by [28].
v o p t i m u m =   ( 2 W 2 C D F a b 2 D a 2 ) 1 4
All the parameters are the same used in (6) and (7). Several energy models have been proposed in the literature. The energy consumed during horizontal flight has been considered approximately constant [23], but it is not true in real case scenarios. The authors in [17] tried to derive a more realistic energy model based on empirical experiments. However, the model is tailored for experimental conditions in which the study has been conducted. In this paper, we will consider Equation (8) to derive an energy model. The energy E consumed during a time T and a power P is given by
E = P T
The time T can also be expressed as
T =   d v + w c o s θ
where d ,   v , w , and θ representing respectively the distance traveled, drone speed relative to the wind, the wind speed, and the relative angle between the wind and the travel direction. Considering Equations (5)–(10), we may express energy consumption in terms of drone velocity relative to wind speed as
E =   P P + P L v + w c o s θ = α v 3 + β v v + w c o s θ
where α = 1 2 C D F a D a and β = W 2 D a b 2
Taking the derivative with respect to v , we find the optimization condition for the optimal speed v o p t
( 3 α v o p t 2 β v o p t 2 ) ( w c o s θ )   + 2 α v o p t 3 2 β v o p t = 0
In the case w = 0 , this gives an optimal velocity v 0 that is expressed in Equation (13),
v 0 = ( β α ) 1 4 ,
Which agrees with (5). In the general case, we may write v o p t = z v 0 . This gives the optimum condition
( 3 α v 0 3 z 2 β v 0 1 z 2 ) ( w c o s θ ) / v 0 + 2 α v 0 3 z 3 2 β v 0 1 z   = 0
Multiplying by v 0 / β and noting that ( α β ) v 0 4 = 1 , we obtain Equation (15)
( 3 z 2 1 z 2 ) ω + 2 z 3 2 z   = 0 ,
where ω ( w c o s θ ) / v 0 . This gives a 5th order Equation (16)
2 z 5 + 3 ω z 4 2 z ω   = 0 .
Denoting the largest root of Equations (16) by z m a x , we obtain the expression of the optimal speed in Equation (17), where the parameter v 0 corresponds to the cruise speed.
v o p t = z m a x v 0
The energy consumption is then given by Equation (18)
E =   ( P P , o p t + P L , o p t ) d v o p t + w c o s θ =   ( α v o p t 3 + β v o p t ) d v o p t + w c o s θ
We use the specifications of the manufacturer to derive α and β . The specifications for the 3DR Solo drone are the following:
  • Cruise speed = 2.5 m/s,
  • Max flight time: 25 min.
  • Energy: 5200 mAh, 14.8 VDC = 77 Wh = 277,200 W-s
From these numbers we may derive the power consumed at 2.5 m/s = 77 Wh/(25/60) = 185 Watts. From (5)–(7) we have
2.5 3 α + β 2.5 = 185
From the cruise speed, we may derive: 2.5 m / s =   ( β α ) 1 4 which implies
2.5 4 α = β
Solving Equations (19) and (20) for α ,   β gives
β = 185 2.5 2 = 231.25 ;     α = β 2.5 4 = 5.92
Finally, the energy consumption is estimated by plugging these parameters into (18)
E = ( 5.92 v o p t 3 + 231.25 v o p t ) d v o p t + w c o s θ
where
v o p t = 2.5 z m a x
and z m a x is the largest root of Equation (24)
2 z 5 + 1.2 ( w c o s θ ) z 4 2 z 0.4 w c o s θ   = 0

2.4. Problem Formulation

The Single Drone Multiple Recharging Stations in Large Farm problem (SD-MRS-LF) is defined as follows: given an area of interest decomposed into a set of cells A and a set of candidate locations for recharging station(s) C , determine a route R that minimizes the number of actually used recharging stations and the completion time of the mission. A route R may be characterized as a permutation of the elements of A . More specifically,   R = ( R 1 , R 2 ,   ,   R n × m ) with R i A ,   1 i n × m .
The objective functions are given by (25) and (26). They respectively minimize the number of recharging stations and the total traveling time.
M i n   f 1 ( R ) = i = 1 n h ( R i )
M i n   f 2 ( R ) = i = 1 n d ( R i ,   R i + 1 ) v o p t R i ,   R i + 1 + w c o s θ
with
h ( R i ) = { 1 ,     |     R i C E 0 ,       | o t h e r w i s e
  • n is the number of elements in R ;
  • R i R ;
  • d ( R i ,   R i + 1 ) is the distance between R i and R i + 1 ;
  • v o p t R i ,   R i + 1 is the optimal speed of the drone to move from point R i to point R i + 1 ;
  • w is the wind speed;
  • θ is the relative angle between the wind and the travel direction;
  • C E is the set of actually used recharging stations.

2.5. Complexity Analysis of the Problem

Let us consider the take-off point and the final landing node of the drone to be the same. The given problem can easily be reduced to the traveling salesman problem (TSP) by considering each cell to cover as a city to visit. Since each cell and recharging station can be visited once only, the problem amounts to find the best arrangement starting from the take-off point considered as the depot to the final landing point which is the same point. The given problem is therefore NP-hard with factorial time complexity O ( n ! ) .
For instance, an area of size 8 × 8 contains 64 cells corresponding to 64 !   =   1.26 × 10 89 arrangements. Evaluating such a number of arrangements using an exhaustive strategy would require   10 72 years on a CPU with a clock speed of 4 gigahertz. That is why metaheuristics are good alternatives.

3. Back-and-Forth-k-opt Simulated Annealing Approach (BFKSA)

3.1. Background

BFKSA is a combination of three algorithms: back-and-forth, k-opt, and simulated annealing. In the following subsections, we describe these three algorithms. In Section 3.2. we will describe how the three algorithms are integrated together to construct the BFKSA algorithm.

3.1.1. Back-and-Forth (BF)

BF is the simple technique mainly used in regular areas. It is inspired by ‘the way of the ox’. In fact, when an ox trails a plow and has to cover an area, it moves in a straight line along the field, turns around, and then traces a new straight path adjacent to the previous one. BF is considered as an energy-aware algorithm since it reduces the number of turns made by a drone. For instance, in Figure 3a,b, the drone makes only six turns to cover the area of interest. A BF route can be generated using Algorithm 1 or Algorithm 2.
A generated route R may not be feasible. A route is feasible if the maximum energy of the battery allows the drone to fly from the starting point or an intermediary recharging station to the closest recharging station on the route; or from the last recharging station to the final landing point. Algorithm 3 minimizes the number of recharging stations as given by (24) on a given route R when the route is feasible.
Depending on the location of the recharging stations, three cases can be observed (Figure 3):
  • The route requires the minimal energy and the minimal number of recharging stations (Figure 3a)
  • The route requires the minimal energy, but a higher number of recharging station (Figure 3b)
  • The route requires the minimal energy, but it is not feasible (Figure 3c)
In the last two cases, modifying the BF route may decrease the number or recharging stations or turn an infeasible route into a feasible one. This can be done by combining back-and-forth with metaheuristics such as simulated annealing, as described below.

3.1.2. Simulated Annealing (SA)

SA is inspired from the annealing process in metallurgy. The idea consists of generating an initial route and then performing several iterations to improve this route. A neighbor route is generated from the previous one at each iteration. The new route is accepted if it improves the value of the objective function. If it is not the case, the new route is accepted with a probability depending on the current temperature T and the difference between the new and previous values of the objective function (denoted by Δ E ). The probability in the SA algorithm generally follows the Boltzmann distribution given by e Δ E T . Since the temperature updates each time the equilibrium state is reached, the probability of accepting non-improving routes decreases.
Algorithm 1: Back and Forth Vertical (BFv)
Input: A : Area to cover
Output: R : A route;
Symmetry 12 01661 i001
Algorithm 2: Back and Forth Horizontal (BFh)
Input: A : Area to cover
Output: R : A route;
Symmetry 12 01661 i002
Algorithm 3: Recharging Station Minimization (RSM)
Input: R : a route; C : Set of Candidate RS
Output: C E : set of actual RS
Symmetry 12 01661 i003

3.1.3. K-Opt Algorithm

The basic simulated annealing randomly selects a new solution x in the neighborhood of the current solution x . This selection process presents some limitations. First, the same neighbor solution x can be evaluated several times since there is no mechanism to recognize already visited solutions. In some special cases, the algorithm can even be trapped in a local loop, switching between x and x . Secondly, by selecting only one solution in the neighborhood of x , another neighbor solution x better than x can be ignored. To mitigate these limitations, a local search in the neighborhood of a current solution x is usually performed.
One of the most widely used local search algorithms is the k-opt algorithm. It is a generalization of the 2-opt algorithm proposed by Georges Croes in 1958 to solve the traveling salesman problem. The algorithm is based on the k-opt operator and works as follows: at each step, k edges are removed from the current solution and the k partitions are therefore reconnected to form a new solution. An application of the 2-opt and the 3-opt operator on a given route R is illustrated in Figure 4.

3.2. BFKSA Algorithm

The Back-and-Forth Simulated Annealing is given in Algorithm 4.
Initial Route (line 2): The initial route is obtained using Algorithm 5 (GenerateInitialRoute). BF algorithm can consider either vertical orientation (Figure 5b) or horizontal orientation (Figure 5c). The idea in Algorithm 5 is to run BF in both directions and to consider the orientation that provide the best value for f 1 and f 2 depending on wind direction. In case there is no feasible solution, the orientation that forms the smallest angle with the wind direction is taken. In this paper, we consider three directions for the wind: (0,1), (1,1), (1,0), shown in Figure 5a.
Fitness functions (lines 3 and 11): the two objective functions namely f 1 and f 2 are evaluated using respectively Equations (3) and (4). The number of actual recharging stations in a route is reduced using Algorithm 2 (recharging station minimization).
Generating a new route (line 8–10): the process of finding a new route starts by randomly selecting three indexes that will correspond to the edges to delete on the current route. Then, by using GenerateRoutes(), new solutions are derived from the current one using a k-opt operator that deletes k edges and builds a new route. The value of k varies between 2 and 3. For a given route R, all the routes derived by applying a 2-opt or 3-opt are tested and the best is kept. Since the proposed model considers the wind direction, a reverse edge may consume a different energy quantity compared to the original edge (Figure 4). The determination of the best local route is given in Algorithm 6.
Stopping condition (line 6): the algorithm stops when the current temperature is lower than the minimal temperature or after reaching several equilibrium stages without improvement. We suppose therefore having reached at least a close to optimal route.
Equilibrium stage (line 7): If there is no improvement of the best route found so far after a certain number of iteration, we suppose having reached an equilibrium state at the current temperature and the temperature is update (line 19).
Accepting a solution (line 14 to 18): Since the problem formulation first focuses on the cost of the deployment, a route is accepted if it provides a smaller number of recharging station, or if the number of recharging station is equal to the one of the best route found so far, but with a smaller completion time. However, a non-improving route with regard to the completion time can be accepted in order to escape from local optima, but with a certain probability.
Algorithm 4: Back-and-Forth-k-opt Simulated Annealing
Input: f 1 , f 2   : Objective functions;   C : Set of Candidate RS; A : area to cover
Output: R B : best route found, D B : Corresponding time, C E B : minimal set of actual RS
Symmetry 12 01661 i004
Algorithm 5: GenerateInitialRoute
Input: f 1 , f 2   : Objective functions; A : Area to cover;     C : Set of Candidate RS; w : wind direction
Output: R : A route;
Symmetry 12 01661 i005
Accepting a solution (line 12 to 16): A route is accepted if it provides a smaller number of recharging stations, or if the number of recharging station is the equal to the one of the best route found so far, but with a smaller energy consumption. However, a non-improving route with regard to energy consumption can be accepted in order to escape from local optima, but with a certain probability.
Algorithm 6: BestLocalRoute
Input: f 1 , f 2   : Objective functions;     C : Set of Candidate RS; [ R 1 R 8 ]: Routes from 3-opt
Output: R B L : best route from 3-opt
Symmetry 12 01661 i006

4. Performance Evaluation

4.1. Parameter Settings and General Setup

Parameter values used in the simulations are summarized in Table 1. The maximum energy, optimum velocity, and altitude correspond to operating parameters of the 3DR Solo drone. Vertical velocity and acceleration were taken as 0; horizontal acceleration was taken as 0, except when the drone is approaching a cell containing a RS.
A set of 20 random topologies of sizes between 8 × 8 and 16 × 16 cells with a varying number of candidate recharging station locations was generated using a generalized Halton number generator (https://github.com/fmder/ghalton). We consider an altitude H = 120   m as adopted in [22] for two reasons. First, 3DR solo drone can be equipped with a camera able to shoot 4 K videos at 30 fps and 12MP. This is far above 1.2 MP provided by the embedded camera of the ReadyMadeRC Anaconda drone used in [22]. Secondly, although the maximum altitude H M can be higher than 120 m according to Equation (4), due to possibly legal restrictions and the maximum altitude of the drone we limit the altitude of the drone to 120 m. Using an angle of view θ = 94.4 ° , the size of a cell is around 65 m (0.4 ha) according to Equation (1). That means the total area varies between 27 ha ( 8 × 8 instances) and 108 ha ( 16 × 16 instances).
Wind direction can influence the energy consumption of the drone for higher wind speeds. In this work, the wind speed ranges from 1 to 4 m/s. In addition, we consider three prevailing wind directions: 0°, 45°, 90° (Figure 5a).
All algorithms were implemented in Python, and a UI was developed to visualize the results. Figure 6 presents a screenshot of the UI.

4.2. Simulation Results and Discussion

The proposed BFKSA is compared to the naïve Back and Forth approach and the K-Simulated Annealing (KSA). KSA is the K-opt variant of the well-known simulated annealing (KSA). The main difference with BFKSA is the fact that KSA starts with a randomly generated initial solution. All the algorithms are implemented from scratch in Python. The source code is available on GitHub at the following address: https://github.com/Suns-group/SD-MRS-LF (Supplementary Materials).
Four metrics are used to compare approaches: the ability to find a feasible solution, the number of recharging stations, the completion time of the mission, and the total energy consumption.
Simulations results are discussed based on the size of the instances and the wind speed and direction. Figure 7, Figure 8, Figure 9 and Figure 10 display the results for feasibility, number of charging stations, mission completion time, and total energy consumption respectively on 8 × 8 instances. Figure 7 shows that the higher the wind speed, the lower the ability of all algorithms to find feasible solutions for all the approaches. This is explained by the fact that the drone must at times travel against the wind, and opposing strong winds can cause very large increases in energy consumption, according to Equation (21). Figure 7 and Figure 8 show that all the approaches provide 100% feasibility and equal numbers of charging stations for wind speed less than or equal to 2 m/s. However, for higher wind speeds the BFKSA outperforms other approaches in both feasibility and reduced number of charging stations, except for wind speed 4 m/s and wind direction 45°. This may be due to the fact that the initial back-and-forth movement used by BFKSA either aligns with 90° or 0°, whereas the optimal direction to avoid wind drag should be orthogonal to the wind direction 45°. It is plausible that the initial suboptimal alignment leads to the algorithm’s convergence to a suboptimal solution. On the other hand, KSA starts with a complete random solution and has a better chance of finding optimal drone path alignments.
Figure 9 and Figure 10 show the completion time and energy consumption respectively for 8 × 8 instances for the different wind scenarios. It is known from the literature that back and forth movement provides the most efficient energy consumption and mission completion time when wind speed is negligible. However, the figures show that BFKSA can provide a slight improvement in some cases for lower wind speeds. For other wind speeds, BFKSA gives results that are close to BF and better than KS, except for power consumption with wind direction 45°.
Detailed results for the 8 × 8 case are listed in Table A1, Table A2, Table A3, Table A4, Table A5 and Table A6 in Appendix A. An examination of the tables shows that BFKSA does not provide the smallest mission completion time only in two scenarios out of 12. The tables also show that each time BF finds a feasible solution, its energy consumption is close to the one of BFKSA, except in the cases where BFKSA provides a smaller number of recharging stations.
Figure 11, Figure 12, Figure 13 and Figure 14 display the results for feasibility, number of charging stations, mission completion time, and total energy consumption respectively on 16 × 16 instances. As was the case in 8 × 8 instances, the wind speed and direction heavily influence the ability to find feasible solutions. For wind speed higher than 2 m/s and wind direction different from 0°, none of the approaches is able to find a feasible solution. This justifies the holes in Figure 12, Figure 13 and Figure 14 for wind speed higher than 2 m/s and wind direction different from 0°.
However, BFKSA manages to obtain close to 100% feasibility higher wind speeds and direction equals to 0°. In addition, Figure 12, Figure 13 and Figure 14 show that BFKSA provides the minimal number of recharging stations in all scenarios, as well as nearly matching BF in terms of completion time and energy consumption.
More detailed results for the 16 × 16 simulations are given in Table A7, Table A8, Table A9, Table A10, Table A11 and Table A12. The comparative performance for KSA is much worse for 16 × than for 8 × 8. In particular, KSA sometimes requires more than twice the number of recharging stations and the total energy required by BFKSA.
All the above results show the dominance of BFKSA over the other approaches due to its ability to find a feasible solution with a minimum number of recharging stations.

5. Conclusions

The automatic surveillance of large farms using drones requires intelligent path planning that takes into account the drone’s limited range and optimizes the locations of recharging stations. In this work, we introduced the Single Drone Multiple Recharging Stations in Large Farm problem (SD-MRS-LF). The objectives were to minimize first the number of recharging stations and then the mission completion time. An approach called BFKSA based on a combination of BF (back-and-forth), K-opt, and simulated annealing (SA) has been proposed to solve the problem. Simulation results comparing BFKSA to BF and KSA have shown the superior performance of BFKSA with respect to its ability to find feasible solutions with a minimal number of charging stations. BFKSA provides the highest percentage of feasible solutions with the smallest number of recharging stations. In certain cases, BFKSA was able to provide a mission completion time and energy consumption smaller than the values provided by BF known to be effective solution.
In this work, the drone was constrained to visit each station only once. Removing this constraint can significantly reduce the number of recharging stations, but it requires a new formulation of the problem. In addition, we suppose each recharging station having an unlimited energy. In case the drone can visit the same recharging station more than once, a charging and discharging function of the recharging station should be considered for a more realistic scenario. Moreover, other optimization techniques can be explored, and the area of interest can be extended to random fields.

Supplementary Materials

The source code is available online at: https://github.com/Suns-group/SD-MRS-LF.

Author Contributions

Conceptualization, J.L.E.K.F.; Methodology, J.L.E.K.F., C.T., and I.K.B.; Software, I.K.B. and J.L.E.K.F.; Validation, J.L.E.K.F., I.K.B., C.T., M.D.F., and A.F; Formal analysis, J.L.E.K.F. and C.T.; Writing—original draft preparation, J.L.E.K.F.; Writing—review and editing, C.T., M.D.F., and A.F.; Visualization, J.L.E.K.F., C.T., and A.F.; Supervision, A.F; Funding acquisition, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The APC was funded by the University of Bremen.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Results from Simulations

Table A1. Mean values of performance metrics on 8 × 8 instances with wind speed = 3 m/s, direction = 90°.
Table A1. Mean values of performance metrics on 8 × 8 instances with wind speed = 3 m/s, direction = 90°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst8_110010010022.25232.3339.3931.17169.09234.14170.47
Inst8_210010010022.25232.3338.9031.30169.09227.49173.24
Inst8_310010010022.20232.3339.2531.18169.09231.73170.87
Inst8_410010010022.50232.3340.5231.53169.09239.05170.67
Inst8_51009510022.16232.3338.8831.34169.09231.54171.92
Inst8_610010010022.25232.3339.6131.11169.09233.77170.16
Inst8_710010010022.30232.3339.1231.39169.09235.30170.99
Inst8_810010010022.051.9532.3338.0231.17169.09223.23170.71
Inst8_910010010022.40232.3339.7831.53169.09236.90171.13
Inst8_101009510022.47232.3339.5931.04169.09236.41170.65
Table A2. Mean values of performance metrics on 8 × 8 instances with wind speed = 3 m/s, direction = 45°.
Table A2. Mean values of performance metrics on 8 × 8 instances with wind speed = 3 m/s, direction = 45°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst8_110010010012.00129.9834.2229.28127.84193.74127.32
Inst8_210010010011.95129.9834.4029.45127.84191.31126.58
Inst8_310010010011.90129.9833.1029.46127.84185.37128.51
Inst8_410010010022129.9833.9131.31127.84192.25139.31
Inst8_510010010012129.9835.0329.11127.84195.83129.24
Inst8_610010010011.90129.9833.9829.23127.84187.46127.94
Inst8_710010010011.95129.9834.2029.61127.84194.30127.48
Inst8_810010010011.95129.9833.6529.55127.84193.49126.83
Inst8_910010010012129.9835.7329.68127.84209.81127.46
Inst8_1010010010012.10129.9836.0229.35127.84205.57128.17
Table A3. Mean values of performance metrics on 8 × 8 instances with wind speed = 3 m/s, direction = 0°.
Table A3. Mean values of performance metrics on 8 × 8 instances with wind speed = 3 m/s, direction = 0°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst8_110010010022232.3335.3931.41169.09201.00171.01
Inst8_210010010022232.3336.7331.35169.09204.80170.41
Inst8_310010010022232.3334.9031.70169.09202.64171.27
Inst8_40100100-22-37.9333.69-211.94181.85
Inst8_510010010021.90232.3336.4031.42169.09205.60171.49
Inst8_610010010021.90232.3334.1131.44169.09190.71170.56
Inst8_7100100100221.8532.3334.9331.72169.09198.22170.67
Inst8_810010010022232.3334.7731.34169.09197.91168.79
Inst8_910010010022232.3335.2632.04169.09201.67169.65
Inst8_101009010032232.3335.9133.24169.09203.91182.30
Table A4. Mean values of performance metrics on 8 × 8 instances with wind speed = 4 m/s, direction = 90°.
Table A4. Mean values of performance metrics on 8 × 8 instances with wind speed = 4 m/s, direction = 90°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst8_11007510033.07329.7735.3128.80243.69318.33253.66
Inst8_207085-3.072.94-34.6031.90-319.68264.94
Inst8_3065100-3.153-35.7331.84-323.33264.61
Inst8_4035100-2.712.95-34.6932.35-298.96278.56
Inst8_5065100-3.152.95-33.5830.84-301.53258.17
Inst8_6080100-3.062.95-36.1931.45-312.23270.15
Inst8_7090100-3.223-37.0532.45-341.63289.55
Inst8_809595-3.112.89-36.3130.92-314.37271.26
Inst8_9070100-3.072.85-36.6431.94-331.92271.32
Inst8_101004510042.892.9529.7734.3730.90243.69345.96270.49
Table A5. Mean values of performance metrics on 8 × 8 instances with wind speed = 4 m/s, direction = 45°.
Table A5. Mean values of performance metrics on 8 × 8 instances with wind speed = 4 m/s, direction = 45°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst8_107510-3.134.5-33.8025.48-301.36363.18
Inst8_20905-2.944-35.5825.96-306.98383.01
Inst8_30950-2.89--35.15--295.41-
Inst8_40500-2.70--32.68--283.64-
Inst8_507025-2.864-32.5827.33-283.07383.47
Inst8_605010-3.104-32.8425.49-304.37359.28
Inst8_709010-3.064-35.0425.87-314.40382.36
Inst8_80900-2.89--34.71--308.74-
Inst8_90800-2.88--34.18--312.63-
Inst8_100450-2.78--33.54--294.10-
Table A6. Mean values of performance metrics on 8 × 8 instances with wind speed = 4 m/s, direction = 0°.
Table A6. Mean values of performance metrics on 8 × 8 instances with wind speed = 4 m/s, direction = 0°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst8_11004510033.562.829.7735.7629.71243.69367.00247.71
Inst8_21005010033.70329.7736.0029.47243.69349.26244.48
Inst8_3045100-3.782.9-37.1831.94-377.64277.93
Inst8_400100--3--32.23--267.14
Inst8_51002510033.802.629.7735.0529.98243.69372.78246.58
Inst8_61005010033.702.829.7737.1730.06243.69369.74252.76
Inst8_71005510033.912.729.7737.5029.70243.69361.09244.90
Inst8_81006010043.83329.7737.5930.63243.69358.71256.23
Inst8_91003010033.332.929.7735.8329.97243.69325.65244.57
Inst8_1010001003-2.8529.77-29.81243.69-243.42
Table A7. Values of performance metrics on 16 × 16 instances with wind speed = 1 m/s, direction = 90°.
Table A7. Values of performance metrics on 16 × 16 instances with wind speed = 1 m/s, direction = 90°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst16_11009510058.954.55117.23213.13118.61365.61688.27370.35
Inst16_21007510059.25117.23210.77117.17365.61699.65365.39
Inst16_310010010058.754.95117.23214.71117.79365.61690.99367.52
Inst16_41006510068.855.1117.23200.31130.75407.59667.21436.48
Inst16_5100100100695117.23214.93123.34365.61692.8385.7
Inst16_610010010058.854.6117.23211.59117.51365.61681.78366.53
Inst16_71009010058.894.95117.23211.94117.37365.61689.96366.04
Inst16_810010010058.954.55117.23210.54117.93365.61677.76367.91
Inst16_910090100695117.23209.5123.24365.61679.81385.51
Inst16_101008510058.884.9117.23209.86117.61365.61682.68366.84
Table A8. Values of performance metrics on 16 × 16 instances with wind speed = 1 m/s, direction = 45°.
Table A8. Values of performance metrics on 16 × 16 instances with wind speed = 1 m/s, direction = 45°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst16_11009010058.724.35117.23211.2119.05365.61687.16371.64
Inst16_21007510058.675117.23202.92117.14365.61675.12365.32
Inst16_31009510058.744.95117.23207.91117.37365.61671.58366.12
Inst16_41006510068.465121.41201.46132.76407.59664.88442.41
Inst16_51008510068.414.8117.23203.29123.29365.61661.4385.89
Inst16_610010010058.94.8117.23214.75117.59365.61691.2366.84
Inst16_71008010058.565117.23206.86117.13365.61677.07365.28
Inst16_810010010058.74.5117.23208.23118.93365.61670.6370.98
Inst16_91009510068.845117.23209.34122.77365.61676.53384.13
Inst16_101007510059.24.85117.23208.74117.75365.61689.99367.35
Table A9. Values of performance metrics on 16 × 16 instances with wind speed = 1 m/s, direction = 0°.
Table A9. Values of performance metrics on 16 × 16 instances with wind speed = 1 m/s, direction = 0°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst16_11009510059.054.05117.23211119.53365.61685.21373.21
Inst16_21007010058.864.65117.23204.24118.39365.61682.22369.45
Inst16_31009010069.175117.23214.26123.74365.61694.76387.25
Inst16_4065100-9.085-206.48127.45-678.03398.99
Inst16_51009510058.894.3117.23211.07118.55365.61686.13369.96
Inst16_610010010059.24.65117.23215.61118.04365.61697.04368.24
Inst16_710090100594.55117.23207.57117.66365.61673.79366.93
Inst16_81009510048.954117.23212.27117.2365.61695.5365.51
Inst16_910010010068.95117.23213.94123.15365.61690.82385.23
Inst16_101009010059.064.05117.23213.61119.68365.61700.72373.65
Table A10. Values of performance metrics on 16 × 16 instances with wind speed = 2 m/s, direction = 90°.
Table A10. Values of performance metrics on 16 × 16 instances with wind speed = 2 m/s, direction = 90°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst16_1100301005115117.62205.66117.55387.9800.06387.5
Inst16_210015100611.335117.62200.25122.38387.9849.84408.39
Inst16_310040100610.55117.62207.15126.12387.9808.05423.86
Inst16_401020-118-198.72125.8-795.34606.15
Inst16_510035100610.715117.62205.72123.41387.9790.25411.09
Inst16_6100451006115117.62208.98122.62387.9840.03409.37
Inst16_7100101005105117.62201.56117.51387.9758.86387.4
Inst16_810030100510.175117.62203.88117.48387.9783.53387.43
Inst16_910015100511.334.95117.62200.54117.84387.9833.57388.05
Inst16_10100251005115117.62200.55117.49387.9825.01387.44
Table A11. Values of performance metrics on 16 × 16 instances with wind speed = 2 m/s, direction = 45°.
Table A11. Values of performance metrics on 16 × 16 instances with wind speed = 2 m/s, direction = 45°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst16_110045100511.115117.62209.6117.6387.9846.94387.8
Inst16_210051006125117.62201.71121.6387.9855.09405.19
Inst16_310055100610.735117.62204.04126.41387.9805.49425.79
Inst16_401015-118-200.06125.76-803.45615.34
Inst16_510025100610.25117.62202.68122.57387.9785.22411.15
Inst16_6100401006115117.62207.14121.7387.9820.39405.05
Inst16_710020100510.55117.62199.98117.6387.9789.66387.84
Inst16_810040100510.885117.62206117.53387.9800.92387.44
Inst16_910010100511.55117.62204.42117.49387.9856.11387.39
Inst16_1010015100510.675117.62201.68117.58387.9807.72387.69
Table A12. Values of performance metrics on 16 × 16 instances with wind speed = 2 m/s, direction = 0°.
Table A12. Values of performance metrics on 16 × 16 instances with wind speed = 2 m/s, direction = 0°.
Feasibility (%)Number of Recharging StationCompletion Time (min)Energy Consumption (Wh)
BFKSABFKSABFKSABFKSABFKSABFKSABFKSABFKSA
Inst16_110010100511.55117.62210.69117.6387.9860.01387.84
Inst16_210010100511.55117.62204.61117.55387.9921.26387.7
Inst16_310020100611.55117.62204.84124.93387.9843.27415.64
Inst16_40560-105.92-211.42125.72-776.78444.95
Inst16_510035100610.865117.62211.5122.64387.9867.79407.84
Inst16_610020100512.255117.62212.81117.48387.9900.99387.34
Inst16_710020100811.58124.53210.44124.47585.97840.68586
Inst16_810030100511.834.75117.62211.65118.12387.9873.25389.2
Inst16_910010100610.55117.62206.68123.14387.9827.71410.95
Inst16_1010001006-5117.62-123.23387.9-409.35

References

  1. Mazur, M.; Wisniewski, A.; McMillan, J. Clarity from above: PwC global report on the commercial applications of drone technology. Wars. Drone Powered Solut. PriceWater House Coopers 2016. Available online: https://www.pwc.pl/pl/pdf/clarity-from-above-pwc.pdf (accessed on 20 May 2020).
  2. Lottes, P.; Khanna, R.; Pfeifer, J.; Siegwart, R.; Stachniss, C. UAV-based crop and weed classification for smart farming. In Proceedings of the 2017 IEEE International Conference on Robotics and Automation (ICRA); IEEE: New York, NY, USA; pp. 3024–3031.
  3. Becerra, V.M. Autonomous Control of Unmanned Aerial Vehicles. Electronics 2019, 8, 452. [Google Scholar] [CrossRef] [Green Version]
  4. Choset, H. Coverage for robotics–A survey of recent results. Ann. Math. Artif. Intell. 2001, 31, 113–126. [Google Scholar] [CrossRef]
  5. Kakaes, K.; Greenword, F.; Lippincott, M.; Dosemagen, S.; Meier, P.; Wich, S.A. Drones and Aerial Observation: New Technologies for Property Rights. Hum. Rights Glob. Dev. Primer New Am. 2015, 2015, 514–519. [Google Scholar]
  6. Huang, W.H. Optimal line-sweep-based decompositions for coverage algorithms. In Proceedings of the 2001 ICRA. IEEE International Conference on Robotics and Automation (Cat. No. 01CH37164); IEEE: New York, NY, USA, 2001; Volume 1, pp. 27–32. [Google Scholar]
  7. Galceran, E.; Carreras, M. A survey on coverage path planning for robotics. Robot. Auton. Syst. 2013, 61, 1258–1276. [Google Scholar] [CrossRef] [Green Version]
  8. Jiao, Y.-S.; Wang, X.-M.; Chen, H.; Li, Y. Research on the coverage path planning of uavs for polygon areas. In Proceedings of the 2010 5th IEEE Conference on Industrial Electronics and Applications; IEEE: New York, NY, USA, 2010; pp. 1467–1472. [Google Scholar]
  9. Sadat, S.A.; Wawerla, J.; Vaughan, R. Fractal trajectories for online non-uniform aerial coverage. In Proceedings of the 2015 IEEE International Conference on Robotics and Automation (ICRA); IEEE: New York, NY, USA, 2015; pp. 2971–2976. [Google Scholar]
  10. Cabreira, T.M.; Brisolara, L.B.; Ferreira, P.R., Jr. Survey on Coverage Path Planning with Unmanned Aerial Vehicles. Drones 2019, 3, 4. [Google Scholar] [CrossRef] [Green Version]
  11. Bushaw, J.D.; Ringelman, K.M.; Rohwer, F.C. Applications of Unmanned Aerial Vehicles to Survey Mesocarnivores. Drones 2019, 3, 28. [Google Scholar] [CrossRef] [Green Version]
  12. Collins, M.D. Using a Drone to Search for the Ivory-Billed Woodpecker (Campephilus principalis). Drones 2018, 2, 11. [Google Scholar] [CrossRef] [Green Version]
  13. Choi, C.H.; Jang, H.J.; Lim, S.G.; Lim, H.C.; Cho, S.H.; Gaponov, I. Automatic wireless drone charging station creating essential environment for continuous drone operation. In Proceedings of the 2016 International Conference on Control, Automation and Information Sciences (ICCAIS); IEEE: New York, NY, USA, 2016; pp. 132–136. [Google Scholar]
  14. Raciti, A.; Rizzo, S.A.; Susinni, G. Drone charging stations over the buildings based on a wireless power transfer system. In Proceedings of the 2018 IEEE/IAS 54th Industrial and Commercial Power Systems Technical Conference (I&CPS); IEEE: New York, NY, USA, 2018; pp. 1–6. [Google Scholar]
  15. Feng, Y.; Zhang, C.; Baek, S.; Rawashdeh, S.; Mohammadi, A. Autonomous Landing of a UAV on a Moving Platform Using Model Predictive Control. Drones 2018, 2, 34. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, S.J.; Lim, G.J. A Hybrid Battery Charging Approach for Drone-Aided Border Surveillance Scheduling. Drones 2018, 2, 38. [Google Scholar] [CrossRef] [Green Version]
  17. Tseng, C.-M.; Chau, C.-K.; Elbassioni, K.M.; Khonji, M. Flight tour planning with recharging optimization for battery-operated autonomous drones. CoRR 2017, Abs170310049. [Google Scholar]
  18. Trotta, A.; Felice, M.D.; Montori, F.; Chowdhury, K.R.; Bononi, L. Joint Coverage, Connectivity, and Charging Strategies for Distributed UAV Networks. IEEE Trans. Robot. 2018, 34, 883–900. [Google Scholar] [CrossRef]
  19. Hong, I.; Kuby, M.; Murray, A.T. A range-restricted recharging station coverage model for drone delivery service planning. Transp. Res. Part. C Emerg. Technol. 2018, 90, 198–212. [Google Scholar] [CrossRef]
  20. Editor, D. Pilot-free, Heisha C200 Drone Charging Solutions for Automated Drone Applications. Droneblog. Available online: https://www.droneblog.com/2019/06/28/pilot-free-heisha-c200-drone-charging-solutions-for-automated-drone-applications/ (accessed on 2 May 2020).
  21. Hwang, M.; Cha, H.-R.; Jung, S.Y. Practical endurance estimation for minimizing energy consumption of multirotor unmanned aerial vehicles. Energies 2018, 11, 2221. [Google Scholar] [CrossRef] [Green Version]
  22. Dorling, K.; Heinrichs, J.; Messier, G.G.; Magierowski, S. Vehicle routing problems for drone delivery. IEEE Trans. Syst. Man Cybern. Syst. 2016, 47, 70–85. [Google Scholar] [CrossRef] [Green Version]
  23. Tseng, C.-M.; Chau, C.-K.; Elbassioni, K.; Khonji, M. Autonomous Recharging and Flight Mission Planning for Battery-operated Autonomous Drones. arXiv 2017, arXiv:170310049 Cs. [Google Scholar]
  24. Yang, C.-H.; Tsai, M.-H.; Kang, S.-C.; Hung, C.-Y. UAV path planning method for digital terrain model reconstruction–A debris fan example. Autom. Constr. 2018, 93, 214–230. [Google Scholar] [CrossRef]
  25. Shafian, S.; Rajan, N.; Schnell, R.; Bagavathiannan, M.; Valasek, J.; Shi, Y.; Olsenholler, J. Unmanned aerial systems-based remote sensing for monitoring sorghum growth and development. PLoS ONE 2018, 13, e0196605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tennekes, H. The Simple Science of Flight: From Insects to Jumbo Jets; MIT Press: Cambridge, MA, USA, 2009. [Google Scholar]
  27. Hill, P.G.; Peterson, C.R. Mechanics and Thermodynamics of Propulsion; Addison-Wesley Publishing Company: Boston, MA, USA, 1992; p. 764. [Google Scholar]
  28. Thibbotuwawa, A.; Nielsen, P.; Zbigniew, B.; Bocewicz, G. Energy consumption in unmanned aerial vehicles: A review of energy consumption models and their relation to the UAV routing. In Proceedings of the International Conference on Information Systems Architecture and Technology; Springer: Berlin, Germany, 2018; pp. 173–184. [Google Scholar]
Figure 1. Description of the area of interest: (a) area to cover with candidate location for recharging stations; (b) numbering of the decomposed cells.
Figure 1. Description of the area of interest: (a) area to cover with candidate location for recharging stations; (b) numbering of the decomposed cells.
Symmetry 12 01661 g001
Figure 2. Consecutive positions of drone’s camera.
Figure 2. Consecutive positions of drone’s camera.
Symmetry 12 01661 g002
Figure 3. Back-and-forth solutions on a 4 × 4 grid: (a) an instance with a BFv solution using 3 RS; (b) an instance with a BFv solution using 4 RS; (c) an instance with no feasible BF solution neither BFv (Black) nor BFh (Yellow).
Figure 3. Back-and-forth solutions on a 4 × 4 grid: (a) an instance with a BFv solution using 3 RS; (b) an instance with a BFv solution using 4 RS; (c) an instance with no feasible BF solution neither BFv (Black) nor BFh (Yellow).
Symmetry 12 01661 g003
Figure 4. Application of 2-opt and 3-opt operators on a partial instance of size 6   ×   6. (a) initial route; (b) a neighbor route by deleting edges (25,31) and (32,26); (c) a neighbor route by deleting edges (19,25), (31,32), and (26,20); (d) a neighbor route by deleting edges (19,25), (31,32), and (26,20) without reversing remaining edges.
Figure 4. Application of 2-opt and 3-opt operators on a partial instance of size 6   ×   6. (a) initial route; (b) a neighbor route by deleting edges (25,31) and (32,26); (c) a neighbor route by deleting edges (19,25), (31,32), and (26,20); (d) a neighbor route by deleting edges (19,25), (31,32), and (26,20) without reversing remaining edges.
Symmetry 12 01661 g004
Figure 5. Wind’s influence on the determination of the BF orientation. (a) Three wind’s directions are considered; (b) BF orientation for wind direction (0,1); (c) BF orientation for wind direction (1,0).
Figure 5. Wind’s influence on the determination of the BF orientation. (a) Three wind’s directions are considered; (b) BF orientation for wind direction (0,1); (c) BF orientation for wind direction (1,0).
Symmetry 12 01661 g005
Figure 6. UI for visualization. The blue cell represents the take-off and landing point. Red cells are candidate locations of recharging stations. White cells are locations that are actually used for recharging stations. (a) Drone path generated respectively by BF using two recharging stations. (b) Drone path generated respectively by BFKSA using only one recharging station.
Figure 6. UI for visualization. The blue cell represents the take-off and landing point. Red cells are candidate locations of recharging stations. White cells are locations that are actually used for recharging stations. (a) Drone path generated respectively by BF using two recharging stations. (b) Drone path generated respectively by BFKSA using only one recharging station.
Symmetry 12 01661 g006
Figure 7. Feasibility of the solutions for 8 × 8 instances with different wind speeds and directions.
Figure 7. Feasibility of the solutions for 8 × 8 instances with different wind speeds and directions.
Symmetry 12 01661 g007
Figure 8. Number of recharging stations for 8 × 8 Instances with different wind speeds and directions.
Figure 8. Number of recharging stations for 8 × 8 Instances with different wind speeds and directions.
Symmetry 12 01661 g008
Figure 9. Mission completion time for 8 × 8 instances with different wind speeds and directions.
Figure 9. Mission completion time for 8 × 8 instances with different wind speeds and directions.
Symmetry 12 01661 g009
Figure 10. Total energy consumption for 8 × 8 Instances with different wind speeds and directions.
Figure 10. Total energy consumption for 8 × 8 Instances with different wind speeds and directions.
Symmetry 12 01661 g010
Figure 11. Feasibility of the solutions for 16 × 16 Instances with different wind speeds and directions.
Figure 11. Feasibility of the solutions for 16 × 16 Instances with different wind speeds and directions.
Symmetry 12 01661 g011
Figure 12. Number of recharging station for 16 × 16 instances with different wind speeds and directions.
Figure 12. Number of recharging station for 16 × 16 instances with different wind speeds and directions.
Symmetry 12 01661 g012
Figure 13. Mission completion time for 16 × 16 instances with different wind speeds and directions.
Figure 13. Mission completion time for 16 × 16 instances with different wind speeds and directions.
Symmetry 12 01661 g013
Figure 14. Total energy consumption for 16 × 16 instances with different wind speeds and directions.
Figure 14. Total energy consumption for 16 × 16 instances with different wind speeds and directions.
Symmetry 12 01661 g014
Table 1. Parameter values
Table 1. Parameter values
ParameterValue
Maximum Energy76.96 Wh
Altitude 120 m
Angle of view94.4°
Cell size65 m
Cruise speed2.5 m/s
Wind speed1–4 m/s
Wind direction0°, 45°, 90°
Size8 × 8; 16 × 16

Share and Cite

MDPI and ACS Style

Fendji, J.L.E.K.; Bayaola, I.K.; Thron, C.; Fendji, M.D.; Förster, A. Cost-Effective Placement of Recharging Stations in Drone Path Planning for Surveillance Missions on Large Farms. Symmetry 2020, 12, 1661. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12101661

AMA Style

Fendji JLEK, Bayaola IK, Thron C, Fendji MD, Förster A. Cost-Effective Placement of Recharging Stations in Drone Path Planning for Surveillance Missions on Large Farms. Symmetry. 2020; 12(10):1661. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12101661

Chicago/Turabian Style

Fendji, Jean Louis Ebongue Kedieng, Israel Kolaigue Bayaola, Christopher Thron, Marie Danielle Fendji, and Anna Förster. 2020. "Cost-Effective Placement of Recharging Stations in Drone Path Planning for Surveillance Missions on Large Farms" Symmetry 12, no. 10: 1661. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12101661

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