Completion Time Minimization for UAV-UGV-Enabled Data Collection

In unmanned aerial vehicle (UAV)-enabled data collection systems, situations where sensor nodes (SNs) cannot upload their data successfully to the UAV may exist, due to factors such as SNs’ insufficient energy and the UAV’s minimum flight altitude. In this paper, an unmanned ground vehicle (UGV)-UAV-enabled data collection system is studied, where data collection missions are conducted by a UAV and a UGV cooperatively. Two cooperative strategies are proposed, i.e., collaboration without information interaction, and collaboration with information interaction. In the first strategy, the UGV collects data from remote SNs (i.e., the SNs that cannot upload data to the UAV) as well as some normal SNs (i.e., the SNs that can upload data to the UAV), while the UAV only collects data from some normal SNs. Then, they carry the data back to the data center (DC) without interacting with each other. In the second strategy, the UGV only collects data from remote SNs, while transmitting the collected data to the UAV at a data interaction point, then the data are carried back to the DC by the UAV. There are mobile data collection nodes on the ground and in the air, and the task is to find trajectories to minimize the data collection time in the data center. A collaborative strategy selection algorithm, combining a multi-stage-based SN association and UAV-UGV path optimization algorithm, is proposed to solve the problem effectively, where techniques including convex optimization and genetic algorithm are adopted. The simulation result shows that the proposed scheme reduces the mission completion time by 36% compared with the benchmark scheme.


Introduction
Unmanned aerial vehicles (UAVs) and unmanned ground vehicles (UGVs), as mobile actuators, are being deployed in many application areas, including disaster relief [1][2][3][4][5][6], surveillance in war zones [7], and energy-efficient data collection over wireless sensor networks (WSNs) [8]. Energy consumption is one important design aspect for WSNs since battery replacement or recharging for sensor nodes (SNs) is typically inconvenient, costly, and sometimes impossible. By employing a UAV or UGV as a mobile data collector, each SN can directly send its sensed data to it.
However, new challenges arise when a UAV or UGV performs data collection missions alone. To ensure the safety of UAVs and compliance with UAV traffic rules, UAVs are usually specified with a minimum safe flight altitude and no-fly zones [9]. SNs that cannot successfully upload their data to the UAV may exist, due to factors such as SNs' insufficient energy and the UAV's minimum flight altitude. We define these SNs as remote SNs. In addition, SNs may generate data on a regular basis [10], which makes it important to minimize the mission completion time for ongoing data collection missions so that subsequent missions can begin early. Although the UGV can collect remote SNs' data nearby, it will take more time to complete the mission due to the slow movement speed. The UAV and UGV can complement each other, for example, the UAV can move quickly from one place to another and have line-of-sight (LoS) communication links with SNs.

•
We propose a UAV-UGV-enabled data collection system with two cooperative strategies. We aim to minimize the mission completion time by jointly optimizing the UAV-UGV paths, as well as the SN association while ensuring that SNs can successfully upload their data with the remaining energy.

•
To tackle the formulated problem, we propose a simple scheme in that UAV-UGV only collects data while staying near or above SNs. The problem is reduced to finding the optimal stop-over above or near the SNs and the time duration at each staying location, as well as the serving order among these locations.

•
We propose a collaborative strategy selection algorithm, combining a multi-stagebased SN association and UAV-UGV path optimization algorithm, to solve the problem effectively, where techniques including convex optimization and genetic algorithm are adopted.
The rest of the paper is organized as follows. The next section presents a literature review. In Section 3, we introduce the system model and present the problem formulation. In Section 4, we present the proposed algorithm. In Section 5, we present simulation results. Finally, conclusions and future research directions are outlined.

Literature Review
UAV-UGV collaboration systems are increasingly used in both military and civilian applications because they allow coordinated actions to leverage the strengths of each to achieve a common goal [13]. The strong complementarity between UAVs and UGVs offers many new application prospects for UAV-UGV collaboration systems. Our focus is on cooperative path planning in sensor data collection. Path planning is the basis for the research of UAV-UGV-enabled data collection [14].
Considering the heterogeneous and complementary nature of UAVs and UGVs, two types of UAV-UGV-enabled data collection approaches have been extensively studied in the existing literature. The first one is to play separate functional roles. For example, the UGV is set up as a charging station for UAVs to improve the flying range. In this case, the UAV needs to periodically visit the stationary power station acted by the UGV and replenish energy via a battery replacement system, or the UGV acts as a mobile charging station and the UAV needs to periodically rendezvous with the UGV for resupply [15].
In [16], the authors investigate the use of UGVs as mobile charging stations to solve the energy-constrained problem of UAVs. To avoid the introduction of UGVs as charging stations leading to long mission cycles, the combination of classical techniques and modern evolutionary approaches is used to minimize the distance traveled by UGVs while satisfying the UAV charging requirements and full node coverage. In [17], the path planning problem of multiple UGVs as mobile charging stations for a UAV is investigated for a large-area distribution of sensors. An FCURP-MRS model is proposed, which ensures that UGVs can charge UAVs by generating multiple UGV and UAV paths to shorten the mission cycle and minimize the travel distance of UAVs. The authors of [18] investigated a scenario where a UGV acts as a mobile charging station for multiple small UAVs by jointly optimizing the trajectories of the UGV and the UAV to minimize the maximum flight time of multiple UAVs, while ensuring that charging can be completed for multiple UAVs. The authors of [19] considered the charging of a single energy-limited UAV by a UGV that can only travel on a defined grid road, with the UAV responsible for collecting data from sensors. A two-stage algorithm was used to solve this problem. Firstly, a set of charging station locations was found using the UAV range restrictions. Secondly, a mixed-integer linear programming problem was formulated to obtain the paths.
The second approach is to play the same functional role. There are significant differences in terms of detection accuracy, detection angle, and mechanical structure. By effectively coordinating UGVs and UAVs with the same functional role, data collection, inspection, and target tracking can be achieved. When performing data collection, reconnaissance, and tracking missions, a large area can be covered quickly with the UAV, and the UGV can get closer to the target and intervene in it.
In the data collection mission, the authors of [20] proposed a data collection scheme that uses UGVs as mobile cluster heads in WSNs to assist the UAV to save data collection time and energy for SNs. The ground mobile cluster head is responsible for aggregating sensor data together and the UAV obtains the required information from it. Finally, the route of the ground mobile cluster head is optimized according to the network density. In [21], a mixed-integer linear programming (MILP) algorithm was developed by considering the communication constraints between the UAV and UGV to ensure that the data collected by the UAV can be transmitted to the UGV instantaneously, which uses a branch-and-cut strategy to plan the paths of the UAV and the UGV. In [22], a UAV-UGV cooperative environmental monitoring system is proposed, in which both the UAV and UGV are used as mobile sensors. The UAV is responsible for monitoring the airborne atmospheric quality and the UGV is responsible for monitoring the ground-based atmospheric quality, and the system has been extensively tested in field experiments. A heterogeneous system of UAVs and UGVs with the mission of monitoring the terrain on a given path is presented in [23]. Both the UAV and UGV are equipped with cameras that monitor the terrain within their field of view, making terrain mapping more efficient and mapping data more accurate. Path planning for UAV-UGV collaboration systems can be well-generalized to the path length-based traveling salesman problem (TSP) [24] and can be solved by some classical heuristic algorithms.
In summary, most researchers assume that UAVs can always successfully receive data from SNs, which is not true. To ensure the safety of UAVs and compliance with UAV traffic rules, UAVs are usually specified with a minimum safe flight altitude and no-fly zones [9]. SNs that cannot successfully upload their data to the UAV may exist, due to factors such as SNs' insufficient energy and the UAV's minimum flight altitude. We define these SNs as remote SNs. In most studies, there is no joint mission assignment for the UAV and the UGV. In this paper, a UAV-UGV-enabled data collection system is studied, where data collection missions are cooperatively conducted by a UAV and a UGV to minimize the mission completion time by jointly optimizing the selection of the cooperative strategy, the SN association, and the UAV-UGV path.

Problem Description
As shown in Figure 1, one UAV and one UGV are assigned to collect data from K SNs. Some SNs cannot upload data due to their limited transmission power and the limited safe flight altitude of the UAV. Two cooperative strategies are proposed, i.e., collaboration without information interaction, and collaboration with information interaction. In the first strategy, the UGV collects data from remote SNs (i.e., the SNs that cannot upload data to the UAV) as well as some normal SNs (i.e., the SNs that can upload data to the UAV), while the UAV only collects data from some normal SNs. Then, they carry the data back to the data center (DC) without interacting with each other. In the second strategy, the UGV laboration without information interaction, and collaboration with information interaction. In the first strategy, the UGV collects data from remote SNs (i.e., the SNs that cannot upload data to the UAV) as well as some normal SNs (i.e., the SNs that can upload data to the UAV), while the UAV only collects data from some normal SNs. Then, they carry the data back to the data center (DC) without interacting with each other. In the second strategy, the UGV only collects data from remote SNs, while transmitting the collected data to the UAV at a data interaction point, then the data are carried back to the DC by the UAV. The goal is to minimize the mission completion time.

Problem Assumptions
(1) The position of SNs, the remaining energy, and the amount of data to be uploaded are known before the mission is performed. (2) To ensure the safety of the UAV, the UAV flies at a fixed safe altitude.
(3) The UGV and UAV depart from the DC with sufficient energy to complete the mission. (4) The UAV and UGV collect the data in an orthogonal frequency channel, and SNs only wake up when they are scheduled for data transmission to save energy.

Problem Models
The notations that appear in the model are provided in Table 1 . Assume that the channel between SN k s and UAV is a LoS channel. The channel gain from SN k s to UAV is: where 0 ρ denotes the channel power gain at a reference distance of 1 m, then the SN to UAV upload rate (bits/s/Hz) can be obtained as:

Problem Assumptions
(1) The position of SNs, the remaining energy, and the amount of data to be uploaded are known before the mission is performed. (2) To ensure the safety of the UAV, the UAV flies at a fixed safe altitude.
(3) The UGV and UAV depart from the DC with sufficient energy to complete the mission. (4) The UAV and UGV collect the data in an orthogonal frequency channel, and SNs only wake up when they are scheduled for data transmission to save energy.

Problem Models
The notations that appear in the model are provided in Table 1. Table 1. Model notations.

Notation Type Notation Description
The set of UAV and UGV G m = {s k |a m,k = 1, ∀k} The set of sensors associated with u m K UGV = s k E k < E min k , ∀k The set of remote SNs The movement speed of u m in time slot n H The UAV flight height The amount of data to be uploaded by the SN s k E k The remaining energy of the SN s k E min k The energy required for the SN s k to upload data to the UAV hovering above it P k The transmitting power of the SN s k P UGV The transmitting power of the UGV T m The time for the UAV to return to the DC T K The time for all data to reach the DC T The order of access to u m (ε 1 , ε 2 , . . . , ε |Gm | ) The label arrangement for u m to access SNs Decision variables a m,k Binary. If SN s k is associated with u m , a m,k = 1. }. In addition, the deployment position of SNs is known in advance, and the coordinates of SN s k are expressed as w k = (x k , y k , 0) T . Assume that the channel between SN s k and UAV is a LoS channel. The channel gain from SN s k to UAV is: where ρ 0 denotes the channel power gain at a reference distance of 1 m, then the SN to UAV upload rate (bits/s/Hz) can be obtained as: where σ 2 is additive Gaussian white noise (AWGN) power received by the receiver at the UAV and α is the path loss index. Similarly, the upload rate from SN s k to the UGV and the transmission rate from the UGV to the UAV are: The minimum energy required by the SN s k to transmit its data to the UAV, i.e., the energy consumed by the SN to upload data when the UAV is hovering above the SN is: E min k = D k P k /R hover , where R hover = log 2 1 + P k ρ 0 /σ 2 H 2 . If E k < E min , SN s k data can only be collected by the UGV.
Idle listening is one important cause of energy waste for many sensor network applications. To address this issue, the sleep/wake-up scheme has been introduced [12], which is able to adjust the ratio between the sleeping and awake time of each SN to minimize idle listening. Define a m,k ∈ {0, 1}, ∀u m ∈ M, s k ∈ K. If SN s k is awakened by u m and transmits data to u m , a m,k = 1, and inversely, a m,k = 0. Therefore, we can obtain G m = s k a m,k = 1, ∀k , ∀m. Using Γ m = ε 1 , . . . , ε |G m | to denote the label of the access order of u m , the access order can be determined as (s ε 1 , . . . , s ε |Gm | ). For example, if G m = {s 2 , s 4 , s 6 }, Γ m = {4, 2, 6}, then the u m 's access order is (s 4 , s 2 , s 6 ). To ensure that the data can be uploaded with the remaining energy, the region of feasible staying locations for u m can be obtained as:

Problem Formulation
The SN association and UAV-UGV path are optimized to minimize the mission completion time. The mission completion time is denoted by T K . When strategy A is used, T K = max T m , conversely, T K = T U AV . The problem can be formulated as: The objective function (6) indicates that our goal is to minimize the time for data to reach the DC. Constraint (9) indicates that any SNs can only be associated with either the UAV or UGV. Constraints (10) and (11) indicate that the energy of the SNs is guaranteed to be sufficient to upload the data. This problem is NP-hard due to the presence of integer variables and the discrete nature of Γ m .

Proposed Solution
We propose a simple scheme that UAV-UGV only collects data while staying near or above SNs, termed as staying mode (S-mode). For this mode, to find the optimized stop-over above or near the SNs and the access order among all locations, we propose an efficient algorithm by leveraging the min-max multiple traveling salesman problem (min-max m-TSP) and convex optimization techniques.
Define q ε 0 = q ε |Gm |+1 = q I . The mission completion time consists of the moving time and the collecting time. Since S-mode is used, the moving time and the collecting time are only related to the staying locations, and to minimize the moving time, both the UGV and the UAV move between staying locations at maximum speed. The moving time and the staying time can be obtained as T τ h l = D l /R m,l represents the time used by u m to collect data from SN s l , 1 ≤ l ≤ |G m |.

Two-Stage Optimization Algorithm for Strategy A
The problem is decoupled into two subproblems. In the first subproblem, the SN association {G m } and the access order {Γ m } are obtained by solving Equation (12). In the second subproblem, given the association {G m } and the access order {Γ m } obtained from the first subproblem, the staying locations, q ε l , are optimized to minimize the mission completion time.

Solving for the Association and the Access Order
In the first subproblem, given that q 0 ε l = w ε l , then the problem can be formulated as follows: min s.t. (9), (10) K UGV ∈ G UGV ∈ K The association and access order can be obtained by solving the min-max m-TSP, and although the problem is NP-hard, various algorithms exist to efficiently find high-quality approximate solutions (e.g., genetic algorithms) [25].

Optimizing the UAV-UGV Path
Since τ h l is a non-convex function, and the sum of a finite number of non-convex functions is still a non-convex function, therefore, Equation (6) is a non-convex optimization problem. To solve the problem, relaxation variables {z l } are introduced.
Since the left-hand side of Constraint (16) is a convex function, Constraint (16) is a non-convex constraint. The left-hand side of the equal sign of Constraint (16) is a convex function, which can be subjected to a first-order Taylor expansion and then obtained as a concave lower bound [26], which can be obtained by giving the local point q r ε l .
where R lb ε l is a linear function that is both convex and concave. The equation holds at q ε l = q r ε l , and R lb ε l has the same gradient as R ε l [25]. For any given local point q r ε l and the lower bound expression in Equation (17), the problem is approximated by: Since the objective function and constraint are convex constraints, Equation (20) is a convex problem, which can be efficiently solved using standard convex optimization tools such as CVX [27]. The original non-convex problem can be solved by applying the SCA technique [28] to iteratively optimize Equation (20), updating the local points q r ε l at each iteration until the objective function converges. The overall algorithm for strategy A is summarized in Algorithm 1. Given q r ε l , solve (20) and denote the optimal solution as q r+1

Three-Stage Optimization Algorithm for Strategy B
The association has been provided in strategy B, G UGV = K UGV G UAV = K − G UAV . Then, the problem is decoupled into three subproblems. The first subproblem is to obtain the access order {Γ m }. The second subproblem is to obtain the intersection location {q m c }. The access order and the intersection location are then found by iterative means. Then, the third subproblem is to optimize q ε l and {q m c }. In strategy B, the mission completion time is:

Solving for the Association and the Access Order
The first subproblem is to obtain the access order by solving the TSP given ε |G m | . It can be solved directly by a variety of algorithms, including genetic algorithms.

Solving for Intersection Location
In the second subproblem, the intersection location of the UAV and the UGV is obtained by substituting the results found in the first subproblem. Since the left-hand side of Constraint (27) Since the objective function and constraint are convex constraints, Equation (28) is a convex problem, which can be solved efficiently using standard convex optimization tools.

Optimizing the UAV-UGV Path and Intersection Location
Finally, the optimal staying location q ε l and the interaction point {q m c } are obtained by solving the third subproblem. The relaxation variables x and c k are introduced and a first-order Taylor expansion is performed on R UGV,U AV and R ε l at the given local points (q m c ) r and q r k . We can obtain the following problem: The objective function and constraints are convex, and the problem is convex. Equation (32) can be solved directly using convex optimization tools. The algorithm for strategy B is obtained and is shown in Algorithm 2.

Algorithm 2.
Overall algorithm for strategy B Solve TSP with a genetic algorithm to obtain the order of service {Γ m }; 8.
if T U AV < T 10.

The Collaborative Strategy Selection Algorithm
Based on the results obtained in the previous two subsections, the overall algorithm for solving Equation (6) is proposed. By solving Equation (12) to obtain T UGV and T U AV , if T UGV ≤ T UGV + (τ h UGV ) min , the UAV need not receive data from the UGV, and strategy A is adopted, where (τ h UGV ) min is the time taken by the UGV to transmit the data to the UAV hovering above the UGV; conversely, if T UGV > T U AV + (τ h UGV ) min and G ugv,k = K ugv , then strategy B is adopted. The overall algorithm for Equation (6) is shown in Algorithm 3. Execution of Algorithm 1. 5. else 6.
Execution of Algorithm 2. 7. end if

Simulation and Results
In this section, the performance of the proposed algorithm is evaluated. K = 15 SNs were uniformly and randomly distributed in a square area with a side length equal to 0.8 km. It is assumed that all SNs have the same transmit power. Assume that the remaining energy of the remote nodes is 2/3E min and the sufficient energy of the remaining SNs is 2E min . For comparison, two strategies are simulated in this paper. The main parameters are shown in Table 2. To demonstrate the effectiveness of our proposed scheme, we compared it to three benchmarks, namely the minimum staying time benchmark, the fixed interaction benchmark, and the non-synergistic benchmark.

•
In the minimum staying time benchmark, the UAV/UGV moves to the staying locations of its associated SNs with the maximum speed by solving the min-max m-TSP and stays there for a sufficient duration to complete the data collection.

•
In the fixed interaction benchmark, the interaction point area is set to be midway between the center of the remote SNs and the center of the normal SNs.

•
In the uncooperative benchmark, the UGV is only responsible for collecting data from the remote SNs and bringing the data back itself. Figure 2 shows the UAV-UGV paths and associations for missions with different data sizes. It can be observed that as the D k increases, the UAV needs to spend more time on data collection as the SN requires more communication resources, when the advantage of the UGV arriving close to the data collection is more obvious, as expected. For example, compare Figure 2a,c. In Figure 2a, the SN has fewer data to upload, and the time difference between the UAV and the UGV for data collection is small. The mission completion time depends more on the moving time. In Figure 2c, the SN has more data to upload, the time difference between the UAV and the UGV for data collection is larger, and the mission completion time depends more on the data collection time. Therefore, when the amount of data is small, the UAV is responsible for more nodes. Otherwise, make the UGV responsible for more nodes. Similarly, compare Figure 2b,d. In Figure 2b, it is important to try to avoid the disadvantage of the UGV in terms of movement speed, and therefore to make the distance traveled by the UGV as small as possible, allowing the UAV to move further to receive the data from the UGV and then bring the data back. In Figure 2d, the UGV will collect the data in less time, compensating for the lack of movement speed, while the UAV takes time to collect the data and therefore does not need to UAV to receive the data from the UGV, and instead has the UGV wait on the return path of the UAV. Figure 3 depicts the mission completion time versus the amount of data uploaded by SNs. It can be observed that as the amount of uploaded data increases, the mission completion time increases. When D k ≤ 9 Mbits, the performance of the proposed solution is close to that of strategy B, while when D k > 9 Mbits, the performance of the proposed scheme is close to that of strategy A. This is to be expected. Combine this with the trajectory analysis in Figure 2. When D k ≤ 9 Mbits, the mission completion time is more dependent on the moving time. Comparing Figure 2a,b, to reduce the disadvantage of the slow UGV movement speed, it is important to make the distance traveled by the UGV as small as possible, so when executing strategy A, only the UGV is allowed to collect data from the remote SNs. When executing strategy B, the UGV is made to move as small a distance as possible. When D k > 9 Mbits, the amount of data to be uploaded by the SN is high, and the mission completion time is more dependent on the data collection time. Comparing Figure 2c,d, to reduce the UAV's disadvantage in data collection, the UGV has time to collect data from normal SNs after the remote SNs are collected. In strategy B, the UGV has more time to move to a position that is conducive to reducing the distance traveled by the UAV after collecting data from remote SNs. Since it is more advantageous for the UGV to be nearby for data collection, the use of strategy A to allow the UGV to collect data from more SNs yields better results. Compared with the uncooperative benchmark scheme, the proposed scheme reduces the task completion time by 36% on average.    Figure 2a,b, to reduce the disadvantage of the slow UGV movement speed, it is important to make the distance traveled by the UGV as small as possible, so when executing strategy A, only the UGV is allowed to collect data from the remote SNs. When executing strategy B, the UGV is made to move as small a distance as possible. When 9 k D > Mbits, the amount of data to be uploaded by the SN is high, and the mission completion time is more dependent on the data collection time. Comparing Figure 2c,d, to reduce the UAV's disadvantage in data collection, the UGV has time to collect data from normal SNs after the remote SNs are collected. In strategy B, the UGV has more time to move to a position that is conducive to reducing the distance traveled by the UAV after collecting data from remote SNs. Since it is more advantageous for the UGV to be nearby for data collection, the use of strategy A to allow the UGV to collect data from more SNs yields better results. Compared with the uncooperative benchmark scheme, the proposed scheme reduces the task completion time by 36% on average.   Figure 4 depicts the mission completion time versus the number of insufficient-energy SNs. Strategy A is preferred when the number of remote nodes is low when the UGV must access fewer SNs and the UGV can collect from more SNs. When the number of remote SNs is higher, strategy B is preferred, because the UGV must access more SNs, and the UAV can quickly complete its mission. Therefore, try to keep the UGV execution time as low as possible and let the UAV receive data from the UGV. As the number increases again, the amount of data collected by the UGV increases, and the time to transmit the data to the UAV exceeds that for the UGV to bring back, so strategy A is chosen as superior. Moreover, when the number of remote SNs is low, the mission completion time of strategy B tends to decrease as the number of remote SNs increases. This is because the UGV will finish collecting the data it is responsible for in a short time and then move to the interaction point to wait for the UAV. After all, the UAV is responsible for more SNs. The waiting time of the UGV is enough to collect data from normal SNs. As  Figure 4 depicts the mission completion time versus the number of insufficient-energy SNs. Strategy A is preferred when the number of remote nodes is low when the UGV must access fewer SNs and the UGV can collect from more SNs. When the number of remote SNs is higher, strategy B is preferred, because the UGV must access more SNs, and the UAV can quickly complete its mission. Therefore, try to keep the UGV execution time as low as possible and let the UAV receive data from the UGV. As the number increases again, the amount of data collected by the UGV increases, and the time to transmit the data to the UAV exceeds that for the UGV to bring back, so strategy A is chosen as superior. Moreover, when the number of remote SNs is low, the mission completion time of strategy B tends to decrease as the number of remote SNs increases. This is because the UGV will finish collecting the data it is responsible for in a short time and then move to the interaction point to wait for the UAV. After all, the UAV is responsible for more SNs. The waiting time of the UGV is enough to collect data from normal SNs. As the number of remote SNs increases and the number of SNs responsible for the UAV decreases, the time for the UAV to reach the intersection will advance, while the UGV can still reach the intersection before the UAV. As the number of remote SNs increases, the mission completion time will have a decreasing process. Compared with the uncooperative benchmark scheme, the proposed scheme reduces the task completion time by 23% on average. number of remote SNs is higher, strategy B is preferred, because the UGV must access more SNs, and the UAV can quickly complete its mission. Therefore, try to keep the UGV execution time as low as possible and let the UAV receive data from the UGV. As the number increases again, the amount of data collected by the UGV increases, and the time to transmit the data to the UAV exceeds that for the UGV to bring back, so strategy A is chosen as superior. Moreover, when the number of remote SNs is low, the mission completion time of strategy B tends to decrease as the number of remote SNs increases. This is because the UGV will finish collecting the data it is responsible for in a short time and then move to the interaction point to wait for the UAV. After all, the UAV is responsible for more SNs. The waiting time of the UGV is enough to collect data from normal SNs. As the number of remote SNs increases and the number of SNs responsible for the UAV decreases, the time for the UAV to reach the intersection will advance, while the UGV can still reach the intersection before the UAV. As the number of remote SNs increases, the mission completion time will have a decreasing process. Compared with the uncooperative benchmark scheme, the proposed scheme reduces the task completion time by 23% on average.    Figure 5 depicts the relationship between UAV energy consumption and the amount of data uploaded by SNs. When D k ≤ 9 Mbits, the energy consumption of the UAV in strategy B is smaller than that in strategy A. This is because the UAV has an advantage over the UGV in saving more time when the amount of data to be collected is less, so the UAV will collect data from more SNs. When D k > 9 Mbits, the energy consumption of the UAV in strategy B is greater than that in strategy A. This is because, when the amount of data to be collected is larger, the UGV takes less time, so it will let the UGV move further to collect data from more SNs. The solution proposed in this paper allows the UGV/UAV to take full advantage of its advantages by choosing the strategy to speed up the mission completion time at the cost of consuming more energy of the UAV or the UGV.  Figure 5 depicts the relationship between UAV energy consumption and the amount of data uploaded by SNs. When 9 k D ≤ Mbits, the energy consumption of the UAV in strategy B is smaller than that in strategy A. This is because the UAV has an advantage over the UGV in saving more time when the amount of data to be collected is less, so the UAV will collect data from more SNs. When 9 k D > Mbits, the energy consumption of the UAV in strategy B is greater than that in strategy A. This is because, when the amount of data to be collected is larger, the UGV takes less time, so it will let the UGV move further to collect data from more SNs. The solution proposed in this paper allows the UGV/UAV to take full advantage of its advantages by choosing the strategy to speed up the mission completion time at the cost of consuming more energy of the UAV or the UGV.  Figure 6 depicts the relationship between the degree of aggregation of the remote nodes and the mission completion time. Where the degree of aggregation of the remote nodes is represented by the average distance from the remote nodes to the geometric center of all the remote nodes, the result is obtained by randomly taking 20% of the energy-deficient SNs for simulation several times and then taking the average value. When the remote nodes are all distributed in a small area, the mission completion time for strategy A is shorter. As the degree of aggregation decreases, the UGV takes more time to  Figure 6 depicts the relationship between the degree of aggregation of the remote nodes and the mission completion time. Where the degree of aggregation of the remote nodes is represented by the average distance from the remote nodes to the geometric center of all the remote nodes, the result is obtained by randomly taking 20% of the energydeficient SNs for simulation several times and then taking the average value. When the remote nodes are all distributed in a small area, the mission completion time for strategy A is shorter. As the degree of aggregation decreases, the UGV takes more time to move between the remote nodes. When the time taken by the UGV to perform data collection is close to the time taken by the UAV, strategy B is chosen to allow the UAV to fly further to receive data from the UGV. Compared with the uncooperative benchmark scheme, the proposed scheme reduces the task completion time by 33% on average.  Figure 6 depicts the relationship between the degree of aggregation of the remote nodes and the mission completion time. Where the degree of aggregation of the remote nodes is represented by the average distance from the remote nodes to the geometric center of all the remote nodes, the result is obtained by randomly taking 20% of the energy-deficient SNs for simulation several times and then taking the average value. When the remote nodes are all distributed in a small area, the mission completion time for strategy A is shorter. As the degree of aggregation decreases, the UGV takes more time to move between the remote nodes. When the time taken by the UGV to perform data collection is close to the time taken by the UAV, strategy B is chosen to allow the UAV to fly further to receive data from the UGV. Compared with the uncooperative benchmark scheme, the proposed scheme reduces the task completion time by 33% on average.

Conclusions
In this paper, the data collection mission in UAV-UGV-supported WSNs was investigated to minimize the completion time by jointly optimizing the path, association, and cooperative strategy selection, while ensuring that each SN can successfully upload data with the remaining energy. Two strategies were proposed, the corresponding solution algorithms were presented, and the strategy selection algorithm was proposed. Numerical results showed that the proposed scheme achieved a significant performance improvement over the benchmark scheme. In future work, the problem of minimizing the mission completion time in the more general cases, i.e., when the UAV and the UGV can collect and transmit data while moving, will be explored.