Next Article in Journal
The Effect of Changing the Angle of the Passenger Car Seat Backrest on the Head Trajectories of the 50th Percentile Male Dummy
Next Article in Special Issue
Detecting and Mitigating Attacks on GPS Devices
Previous Article in Journal
Wavefront Changes during a Sustained Reading Task in Presbyopic Eyes
Previous Article in Special Issue
Security Control of Cyber–Physical Systems under Cyber Attacks: A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tessellation-Based Construction of Air Route for Wireless Sensor Networks Employing UAV

Department of Computer Engineering, Dankook University, Yongin Si 16890, Republic of Korea
Sensors 2024, 24(12), 3867; https://doi.org/10.3390/s24123867
Submission received: 10 April 2024 / Revised: 30 May 2024 / Accepted: 10 June 2024 / Published: 14 June 2024

Abstract

:
In this paper, we consider a wireless sensor network consisting of an unmanned aerial vehicle (UAV) acting as a sink node and a number of sensor nodes scattered uncertainly on the ground. In the network, the UAV flies to a spatial point called point of interest and hovers to collect environmental data from neighboring sensor nodes. Then, the UAV proceeds to the next point of interest. The UAV must gather data from all the sensor nodes. On the other hand, a shorter round-trip air route of the UAV is more preferred since a battery-operated UAV needs regular recharging. To satisfy the requirement and to adhere to the recommendation as well, especially in the situation where only vague locational information about sensor nodes is available, we propose a scheme that follows three steps. First, it covers the sensor field of the wireless sensor network with three categories of hexagonal tessellations. Secondly, it establishes a point of interest at the centroid of each tile. Thirdly, it constructs an air route of the UAV, which visits every point of interest along a Hamiltonian cycle on the induced graph. Next, we develop a closed-form expression for the exact flight distance attained by the proposed scheme. For comparative evaluation, we discover some optimal schemes that minimize the flight distance by completely inspecting all patterns and corroborating the property of Hamiltonicity. The flight distance along the air route constructed by the proposed scheme is found to be only slightly longer than the flight distance yielded by an optimal scheme. Furthermore, the proposed scheme is proven to be practically valid when a common multicopter is employed as the sink node.

1. Introduction

A wireless sensor network consists of sink nodes and sensor nodes. In the network, a sensor node amasses environmental information in the vicinity and wirelessly delivers it to a sink node [1]. Accordingly, a sink node collects data from neighboring sensor nodes. Wireless sensor networks have been and are expected to be deployed to various areas, often in harsh environments. For example, a wireless sensor network is used to construct a large database of environmental observations for weather forecasting [2]. A wireless sensor network is also deployed for herd management, assessing the condition of animals and reporting the data to the farm manager [3]. A wireless sensor network is further deployed to a region that is difficult to access due to extreme conditions such as high temperatures, high humidity, and high pressure [4].
Meanwhile, unmanned aerial vehicle (UAV) technologies have spread widely and become popular in various fields, including radar localization, wildfire management, agricultural monitoring, border surveillance, meteorological monitoring, and rescue missions [5]. The benefits brought by the use of UAVs, for example, flexibility in constructing a network, effectiveness in saving energy, and safety for humans, exhort wireless sensor networks to increasingly employ UAVs for collecting data from sensor nodes on the ground [5,6].
In this paper, we consider a wireless sensor network as follows. In the wireless sensor network, sensor nodes are scattered on the ground. As usual, these sensor nodes amass environmental information in the vicinity of them. However, the geographical distribution of sensor nodes is not precisely known, and rough information about the locations of sensor nodes, e.g., indistinct boundary of the region in which sensor nodes sojourn, is only available to us (Such a situation can happen due to a harsh environment in which the wireless sensor network is laid or irregular mobilities of sensor nodes comprising the wireless sensor network). To collect data from sensor nodes, the wireless sensor network employs a single UAV as a sink node. In order to support the UAV in flying over the wireless sensor network, a number of spatial points called points of interest are then established. Also, an air route is constructed for the UAV to visit these points of interest. While flying along the air route, the UAV hovers whenever it arrives at each point of interest. Then, according to a medium access control (MAC) scheme, sensor nodes near the point of interest send the environmental information that they have amassed in the vicinity of them to the lingering UAV. In such a way, the UAV, which acts as a sink node, is able to gather data from sensor nodes.
There have been several studies on the establishment of points of interest and the construction of air routes in a wireless sensor network in which a UAV plays the role of a sink node [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. In [7], heavily loaded sensor nodes were designated as special sensor nodes that a mobile sink would visit. Then, a path for the mobile sink was formed to include only the special sensor nodes. In [9], a grid of geographical points was assumed over a wireless sensor network. Then, some points that a drone would visit were selected, and a path to pass selected points was created so that the sum of flying time and hovering time was minimized. In [12], geographical points at which a mobile sink collects data were built, and a path that passes these points was designed so that the cost incurred by the mobile sink’s flying between any two points as well as the cost brought by the sensor node’s forwarding data to a point. In [16], sensor nodes were clustered in a wireless sensor network. Then, the head of each cluster that a UAV visits was selected so as to maximize the amount of efficiently transmitted data. In [19], a digital twin was created to determine the path for UAVs. A best path was empirically determined based on the experience of the digital twin in a cyberspace that resembles the real environment. Then, information about the best path was delivered to UAVs. The studies as listed above suggested distinct ways of establishing points of interest and constructing an air route for the UAV. However, most of them were commonly carried out under the assumption that considerable information about the locations of sensor nodes was available.
In the wireless sensor network under discussion, points of interest should be configured so that the UAV, which acts as a sink node, is capable of collecting data from all the sensor nodes in principle. On the other hand, it is more preferred to construct a shorter round-trip air route for the UAV in the wireless sensor network since a UAV is typically powered by a battery, which should be recharged or replaced regularly. To satisfy such a requirement and to follow such a recommendation as well, especially in the situation that scanty locational information about sensor nodes is only available to us, we propose a scheme as follows. First, we employ three categories of hexagonal tessellations [20], which bear the property of Hamiltonicity [21] underneath, to cover the sensor field of the wireless sensor network. Next, we establish a point of interest at the center of each tile belonging to the adopted tessellation so that the UAV is definitely able to collect data from all the sensor nodes in the wireless sensor network. Finally, we construct a round-trip air route along a Hamiltonian cycle, which is sought on the employed tessellation, so that the UAV can visit every point of interest while the flight distance, i.e., the length of the air route, is reduced. Then, we develop a closed-form expression of the exact flight distance attained by the proposed scheme. For comparative evaluation, we also discover an optimal scheme that minimizes the flight distance, for some radii of the sensor field by completely inspecting all patterns and corroborating the property of Hamiltonicity. Then, the proposed scheme is compared with the optimal scheme for flight distance. In addition, the proposed scheme is evaluated in the situation where a common multicopter is employed as a sink node in the wireless sensor network.
The main contributions of this paper are as follows:
  • We build three categories of hexagonal tessellations that bear the property of Hamiltonicity underneath.
  • To build an air route for the UAV, we propose a scheme which covers the sensor field with a hexagonal tessellation belonging to one of the three categories, establishes a point of interest at the centroid of each tile, and constructs an air route along a Hamiltonian cycle embedded on the hexagonal tessellation.
  • We present a closed-form expression of the flight distance yielded by the proposed scheme.
  • We find a suboptimal scheme that minimizes the number of points of interest by completely inspecting all patterns. Also, we discover an optimal scheme that minimizes the flight distance by corroborating the Hamiltonicity in suboptimal tessellation.
  • We provide a universal lower bound on the flight distance in a closed form.
The paper is organized as follows. In Section 2, we describe the wireless sensor network in which a UAV plays the role of a sink node. In Section 3, we present a scheme that covers the sensor field of the wireless sensor network with a hexagonal tessellation, establishes a point of interest at the centroid of each hexagonal tile, and constructs a round-trip air route of the UAV along a Hamiltonian cycle embedded in the hexagonal tessellation. In Section 4, we develop a close-form expression of the exact flight distance attained by the proposed scheme. Also, we obtain the flight distance yielded by a scheme, which is optimal in the sense of minimizing the flight distance, by completely inspecting all patterns and corroborating the property of Hamiltonicity as well. Further, we derive a simple lower bound on the flight distance. Section 5 is devoted to numerical examples that comparatively evaluate the proposed scheme.

2. Wireless Sensor Network Employing UAV

In this paper, we consider a wireless sensor network consisting of a UAV and a number of sensor nodes. The UAV, which plays the role of a sink node, flies over the sensor nodes and receives data from them. On the other hand, sensor nodes, which are randomly scattered on the ground, gather environmental data in the vicinity and wirelessly send the data to the UAV.
In the wireless sensor network, the geographical distribution of the sensor nodes is not precisely known, and only marginal information about the locations of sensor nodes is available to us. Assume that only a vague region that encloses all the sensor nodes is known to us in the wireless sensor network under consideration. In addition, assume that the vague region takes the shape of a circular disk with radius α . Define the sensor field of a wireless sensor network to be a region in a geometric shape that encloses all the sensor nodes in the network. Then, we set the vague region to be the sensor field of the wireless sensor network (See Figure 1). Ideally, the sensor field of the wireless sensor network can be regarded as a circular disk whose boundary forms the circumcircle of all the sensor nodes comprising the wireless sensor network. Note that sensor fields can be defined or modeled to have a variety of shapes. For example, a sensor field can be defined as being in the shape of a polygon in which all the sensor nodes are inscribed. Since vague information about the geographic distribution of sensor nodes is assumed to be available, it is appropriate to model the sensor field as a circular disk, which is one of the most general shapes.
The UAV acting as a sink node flies over the sensor field to receive data from the sensor nodes lying on the sensor field. To help the UAV efficiently navigate the sensor field for collecting data, several spatial points, which are called points of interest, are established on the sensor field (See Figure 2). Then, a round-trip air route for the UAV is constructed so that the UAV is able to visit every point of interest (See Figure 3).
The UAV flies along the pre-determined air route and visits the points of interest. As the UAV arrives at a point of interest, it hovers and receives data from the sensor nodes, which reside in a region called the coverage associated with the point of interest (See Figure 2). Note that the coverage associated with a point of interest is determined by many factors, for example, the transmission power of a sensor node, altitude of the UAV, and beamwidth of the antenna at the UAV. We set that the coverage associated with each point of interest identically takes the shape of a circular disk with radius β (See Figure 3). Such a setting can be realized by assuming that the sensor field is flat and the UAV flies at a fixed altitude. In practice, the terrain on which a wireless sensor network is deployed can be mountainous. Then, coverages can have different sizes, even if the UAV navigates at a fixed altitude. In such a case, we assume that each coverage has the minimum among the distinct sizes.
As the UAV finishes collecting data from the sensor nodes sojourning in the coverage associated with a point of interest, the UAV leaves the point of interest and flies to the next point of interest along the air route.

3. Tessellation-Based Establishment of Points of Interest and Construction of Air Route

Consider a wireless sensor network employing a UAV, as depicted in Section 2. For the UAV to collect data from the sensor nodes in the network, several points of interest are established on the sensor field, and an air route is also constructed along the points of interest a priori. As addressed in the introduction, points of interest should be established so that the UAV is able to collect data from all the sensor nodes in the wireless sensor network. Furthermore, it is more preferred to construct a shorter round-trip air route along which the UAV flies. To satisfy the requirement and to follow the recommendation as well, we propose a scheme, denoted by σ , which employs a hexagonal tessellation [20] to cover the sensor field of the wireless sensor network, establishes a point of interest at the centroid of each hexagonal tile, and constructs a round-trip air route of the UAV by finding a Hamiltonian cycle [21].

3.1. Tessellation-Based Establishment of Points of Interest

A tessellation is a collection of tiles in polygonal shapes that are mutually attached without gaps or overlaps [20]. Figure 4 shows a tessellation, which consists of tiles in two polygonal shapes.
In order to establish points of interest in the situation where vague geographical information about the sensor nodes is only available, we consider using a tessellation to cover the sensor field of the wireless sensor network and then setting up a point of interest per tile that comprises the tessellation.
There are a variety of candidate tessellations for covering the sensor field. Remind the requirement that points of interest should be established so that the UAV is able to gather data from all the sensor nodes lying in the wireless sensor network. Also, recall the recommendation that a shorter round-trip air route along which the UAV visits all the points of interest is more preferred. Thus, to cover the sensor field, a candidate tessellation is definitely favored if it consists of a small number of tiles, its tiles take appropriate shapes, and the area of each tile is as large as possible. Hereafter, while bearing the requirement and recommendation in mind, we will gradually narrow candidate tessellations for covering the sensor field down to three classes of hexagonal tessellations, i.e., tessellations consisting of regular hexagons [20]. Then, we propose using these classes of hexagonal tessellations to cover the sensor field and establishing a point of interest at the centroid of each tile.
Tessellations can be divided into monohedral tessellations, i.e., tessellations with tiles in the same polygonal shape, and non-monohedral tessellations [20]. Figure 5 illustrates a monohedral tessellation, which consists of triangular tiles. Also, note that the tessellation in Figure 4 is an example of non-monohedral tessellation.
As the sensor field is covered by a tessellation, a point of interest is established in the convex hull of each tile belonging to the tessellation. Then, the coverage associated with the point of interest, which takes the shape of a circular disk with radius β , is also built such that the point of interest is centered on the coverage. For the UAV to receive data from all the sensor nodes, each tile of the tessellation should be enclosed in the corresponding coverage. Further, to reduce the number of points of interest, it is advantageous that a coverage forms the smallest enclosing disk of the corresponding tile. A non-monohedral tessellation consists of tiles in two or more distinct shapes. In case that the sensor node is covered by such a tessellation, the coverage corresponding to a tile cannot be the smallest enclosing disk as far as the tile is not the largest tile in the tessellation. As a result, the sensor field can be covered by a monohedral tessellation, which consists of a smaller number of tiles than a non-monohedral tessellation. From now on, we thus confine our attention to monohedral tessellations for candidate tessellations to cover the sensor field.
Monohedral tessellations can be categorized into edge-to-edge monohedral tessellations, i.e., monohedral tessellations in which each tile only shares one full side with any adjacent tile, and non-edge-to-edge monohedral tessellations [20]. Figure 6 exhibits a non-edge-to-edge monohedral tessellation consisting of rectangular tiles. The tessellation shown in Figure 5 is a typical edge-to-edge monohedral tessellation.
To reduce the number of points of interest, it is more preferred to employ a tessellation in which tiles are positioned so that coverages are less overlapped. Henceforth, we thus concentrate on edge-to-edge monohedral tessellations for candidate tessellations to cover the sensor field.
A tile can take the shape of either a convex polygon or a concave polygon. Figure 7 demonstrates an edge-to-edge monohedral tessellation comprised of tiles that take the shape of a concave polygon. Note that the tessellation exhibited in Figure 5 is an edge-to-edge monohedral tessellation consisting of tiles that take the shape of a convex polygon.
For any concave tile that is enclosed in a coverage, note that there is a larger convex tile that can be enclosed in the same coverage. Since a larger tile is more favored to reduce the number of points of interest, we focus on edge-to-edge monohedral tessellations with convex tiles hereafter.
Note that tiles in the shape of any triangle can comprise an edge-to-edge monohedral tessellation. Also, tiles in the shape of any convex quadrilateral can form an edge-to-edge monohedral tessellation. On the other hand, 8 types of convex pentagons and 3 types of convex hexagons are known to be able to make up edge-to-edge monohedral tessellations [20,22]. However, there is no edge-to-edge monohedral tessellation that consists of tiles in the shape of a convex polygon with more than 6 sides [20]. Recall that a larger tile is more suitable for reducing the number of points of interest. Among the candidate polygons that can be used for edge-to-edge monohedral tessellations, note that an equilateral equiangular hexagon is the largest polygon that can be enclosed in a coverage that takes the shape of a circular disk. Henceforward, we only consider hexagonal tessellations, i.e., edge-to-edge monohedral tessellations in which tiles are in the shape of a regular (equilateral and equiangular) hexagon [20], for candidate tessellations to cover the sensor field. Note that a candidate hexagonal tessellation consists of hexagonal tiles whose sides are β in length, i.e., hexagonal tiles inscribed in coverages (See Figure 8).
Since a tessellation is more favored if it consists of a smaller number of tiles, we may use a hexagonal tessellation that covers the sensor field with the minimum number of tiles. However, searching for such a hexagonal tessellation is not tractable, especially when the sensor field is relatively large. Furthermore, a minimum number of points of interest does not necessarily invoke the shortest round-trip air route of a UAV. Thus, we would rather confine our attention to three classes of hexagonal tessellations that possess a useful property called Hamiltonicity [21], as follows.
Consider the tiles at the boundary of a hexagonal tessellation. Connect the centroids of each pair of adjacent tiles using a straight-line segment. Then, a polygon can be induced by the tessellation. Recall that the coverage associated with a point of interest takes the shape of a circular disk with radius β . For n N , let T n ( 1 ) denote a hexagonal tessellation that induces a regular hexagon whose sides are 3 n β in length. (Hereafter, we call T n ( 1 ) the hexagonal tessellation of category 1 and scale n .) For n N , let T n ( 2 ) and T n ( 3 ) denote tessellations that induce distinct isogonal hexagons. In the hexagon induced by tessellation T n ( 2 ) , the lengths of 3 sides are commonly 3 n β while the lengths of the other 3 sides are equal to 3 ( n + 1 ) β . On the other hand, tessellation T n ( 3 ) induces the hexagon in which 3 sides are commonly 3 n and the other 3 sides are identically 3 ( n + 2 ) in length (Henceforth, we call T n ( 2 ) ( T n ( 3 ) ) the hexagonal tessellation of categories 2 (category 3 ) and scale n ). Figure 9 illustrates tessellations T 2 ( 1 ) , T 2 ( 2 ) and T 2 ( 3 ) .
In addition, let T 0 ( 1 ) , T 0 ( 2 ) , and T 0 ( 3 ) denote hexagonal tessellations with 1, 3, and 6 hexagonal tiles, respectively, as illustrated in Figure 10.
From now on, we confine our attention to three categories of tessellations { T n 1 : n { 0 } N } , { T n 2 : n { 0 } N } and { T n 3 : n { 0 } N } to cover the sensor field of the wireless sensor network.
The following theorems state some properties exhibited by three categories of hexagonal tessellations when they are used to cover the sensor field. (Later, these properties will be used to derive the flight distance, i.e., the length of a round-trip air route.) First, the theorem below shows the number of tiles which comprise each of T n ( 1 ) , T n ( 2 ) and T n ( 3 ) , or equivalently, the number of points of interest established when the sensor field is covered by each of them.
Theorem 1.
Suppose that the sensor field is covered by tessellation  T n ( c )  for  { 1,2 , 3 }  and  n { 0 } N . Let  K n ( c )  denote the number of points of interest established when tessellation  T n ( c )  is used to cover the sensor field for  c { 1,2 , 3 }  and  n { 0 } N Then,
K n ( 1 ) = 3 n 2 + 3 n + 1 K n ( 2 ) = 3 n 2 + 6 n + 3 K n ( 3 ) = 3 n 2 + 9 n + 6
for  n { 0 } N .
Proof. 
A proof of Theorem 1 is given in Appendix A. □
Secondly, the following theorem reveals the maximum size of the sensor field that can be enclosed in each of T n ( 1 ) , T n ( 2 ) and T n ( 3 ) .
Theorem 2.
For  c { 1,2 , 3 }  and  n { 0 } N let  R n ( c )  denote the maximum radius of the sensor field that can be enclosed in tessellation  T n ( c ) Then,
R 0 ( 1 ) = 3 2 · β R 0 2 = β R 0 ( 3 ) = 3 · β R 2 m 1 ( 1 ) = ( 3 m 1 ) · β R 2 m 1 ( 2 ) = 9 m 2 3 m + 1 · β R 2 m 1 ( 3 ) = 3 m · β R 2 m ( 1 ) = 9 m 2 + 3 m + 1 · β R 2 m ( 2 ) = ( 3 m + 1 ) · β R 2 m ( 3 ) = 9 m 2 + 9 m + 3 · β
for  m N .
Proof. 
A proof of Theorem 2 is given in Appendix B. □
From Theorems 1 and 2, we can obtain an inequality among the numbers of points of interest when the sensor field is covered by hexagonal tessellations belonging to the three categories, respectively. Also, we can derive the relation between the maximum radii of the sensor field that can be enclosed in hexagonal tessellations belonging to the three categories, respectively, as presented in the corollary below.
Corollary 1.
For  n { 0 } N ,
K n ( 1 ) K n 2 K n 3
R n ( 1 ) R n 2 R n 3 .
Proof. 
A proof of Corollary 1 is given in Appendix C. □
Suppose that we cover the sensor field of radius α by a tessellation which belongs to one of the three categories of hexagonal tessellations. Consider the hexagonal tessellation of category Ψ α and scale Λ α , where Ψ α 1,2 , 3 and Λ α { 0 } N uniquely satisfy
K Λ α ( Ψ α ) = m i n K n ( c ) : R n ( c ) α f o r   c 1,2 , 3   a n d   n { 0 } N
for each α ( 0 , ) . Then, we propose scheme σ that uses tessellation T Λ α ( Ψ α ) to cover the sensor field and establishes a point of interest at the centroid of each tile comprising tessellation T Λ α ( Ψ α ) .

3.2. Tessellation-Based Construction of Air Route

In this section, we propose a scheme which constructs a round-trip air route of the UAV for collecting data over the sensor field of the wireless sensor network.
For category c { 1,2 , 3 } and scale n { 0 } N , we can construct a graph induced by tessellation T n ( c ) as follows. Let a node represent a tile of T n ( c ) . Let a link connect a pair of nodes if two tiles corresponding to the nodes are adjacent. Then, such sets of nodes and links form the graph induced by tessellation T n ( c ) . For c { 1,2 , 3 } and n { 0 } N , let G n ( c ) denote the graph induced by tessellation T n ( c ) . Figure 11 demonstrates graphs G 2 ( 1 ) , G 2 ( 2 ) and G 2 ( 3 ) , which are induced by tessellations T 2 ( 1 ) , T 2 ( 2 ) and T 2 ( 3 ) , respectively. Note that graph G n ( c ) is a triangular grid graph for all c { 1,2 , 3 } and n { 0 } N [23].
In a graph, a path is a sequence of distinct nodes in which there is a link between two consecutive nodes. In a graph, a cycle is a special path in which the starting node is identical to the ending node. A cycle in a graph is called Hamiltonian if the cycle visits every node in the graph. A graph is said to be Hamiltonian if a Hamiltonian cycle exists in the graph [21]. The theorem below shows that graphs G n ( 1 ) , G n ( 2 ) and G n ( 3 ) are Hamiltonian.
Theorem 3.
Graphs  G n ( 1 ) , G n ( 2 ) and G n ( 3 ) are Hamiltonian for all n { 0 } N .
Proof. 
A proof of Theorem 3 is given in Appendix D. □
A Hamiltonian cycle in G n ( 1 ) is not necessarily unique. Figure 12 illustrates exemplary Hamiltonian cycles in graphs G 2 ( 1 ) , G 2 ( 2 ) and G 2 ( 3 ) .
Recall category Ψ α and scale Λ α defined in (5). Suppose that we cover the sensor field with tessellation T Λ α ( Ψ α ) . Note that tessellation T Λ α ( Ψ α ) consists of K Λ α ( Ψ α ) tiles. Thus, graph G Λ α ( Ψ α ) which is induced by tessellation T Λ α ( Ψ α ) is comprised by K Λ α ( Ψ α ) nodes. Let { v 1 , , v K Λ α ( Ψ α ) } denote the set of nodes in graph G Λ α ( Ψ α ) . From Theorem 4, there is a Hamiltonian cycle in graph G Λ α ( Ψ α ) . Let { v π 1 , , v π K Λ α ( Ψ α ) , v π 1 } denote a Hamiltonian cycle in graph G Λ α ( Ψ α ) , where { π 1 , , π K Λ α ( Ψ α ) } is a permutation of { 1 , , K Λ α ( Ψ α ) } . For j { 1 , , K Λ α ( Ψ α ) } , let P j denote the spatial location, e.g., triple of longitude, latitude and altitude, of the point of interest which lies at the centroid of the tile corresponding to node v j . Then, we propose scheme σ that constructs a round-trip air route that starts from P π 1 and also ends at P π 1 as follows. The air route sequentially visits P π j for j { 1 , , K Λ α ( Ψ α ) } . From P π j to P π j + 1 , the air route is built as a straight-line segment for j { 1 , , K Λ α ( Ψ α ) 1 } . From P π K Λ α ( Ψ α ) to P π 1 , the air route is also constructed as a straight-line segment.
For example, suppose that radius of the sensor field α is equal to 3 β . Then, category Ψ α and scale Λ α are calculated to be 3 and 0 , respectively. Thus, tessellation T 0 ( 3 ) is used to cover the sensor field according to scheme σ . Then, a point of interest is established at the centroid of each tile in tessellation T 0 ( 3 ) (See Figure 13a). Next, graph G 0 ( 3 ) which consists of 6 nodes, denoted by v 1 , , v 6 , is induced from tessellation T 0 ( 3 ) as shown in Figure 13b. Consider sequence 1,2 , 4,5 , 6,3 , which is a permutation of 1 , 2 , 3 , 4 , 5 , 6 . Then, v 1 , v 2 , v 4 , v 5 , v 6 , v 3 , v 1 is a Hamiltonian cycle on induced graph G 0 ( 3 ) . (See Figure 13c). Finally, a round-trip air route, in which every point of interest is constructed along the Hamiltonian cycle over the sensor field, is shown in Figure 13d.

4. Analysis of Flight Distance

In Section 3, we proposed scheme σ for establishing points of interest over the sensor field and building a round-trip air route that visits every point of interest. In this section, we derive the flight distance attained by scheme σ , i.e., the length of the round-trip air route constructed according to scheme σ . For the comparative evaluation of proposed scheme σ , which will be carried out in the following section, we also consider a family of schemes that use hexagonal tessellations and analyze the flight distances achieved by such schemes as well.
In the wireless sensor network under discussion, recall that the sensor field takes the shape of circular disk with radius α and each coverage is in the shape of a circular disk with radius β . Consider a family of schemes as follows. To cover the sensor field, each scheme belonging to the family uses a hexagonal tessellation in which each hexagonal tile is inscribed in the corresponding coverage, i.e., every side of each tile is β in length. Then, the scheme establishes a point of interest at the centroid of each tile. At last, the scheme constructs the shortest round-trip air route to visit every point of interest.
For i N , consider a set of distinct hexagonal tessellations, each of which consists of i hexagonal tiles with sides of length β . For i N , let ν i denote the number of such distinct hexagonal tessellations. Note that ν 1 = ν 2 = 1 and ν 3 = 3 for example. For size i N , let S i [ p ] : p { 1 , , ν i } denote the set of such distinct hexagonal tessellations. (Hereafter, we call S i [ p ] hexagonal tessellation of pattern p and size i ). For size i N and pattern p { 1 , , ν i } , let σ i [ p ] denote the scheme which uses tessellation S i [ p ] to cover the sensor field. Then, σ i [ p ] : p 1 , , ν i , i N represents the family of schemes that we focus on.
For size i N and pattern p { 1 , , ν i } , let Q i [ p ] denote the maximum radius of the sensor field that can be enclosed in tessellation S i [ p ] . For size i N and pattern p { 1 , , ν i } , let N ( α , σ i [ p ] ) denote the number of points of interest established over the sensor field of radius α by scheme σ i [ p ] if the sensor field can be covered by tessellation S i [ p ] . Then,
N α , σ i [ p ] = i
if Q i [ p ] α . Otherwise, we set N α , σ i [ p ] = . Also, for size i N and pattern p 1 , , ν i , let D ( α , σ i [ p ] ) represent the flight distance attained by scheme σ i [ p ] , i.e., the length of the round-trip air route that visits every point of interest established over the sensor field of radius α according to scheme σ i [ p ] if the sensor field can be covered by tessellation S i [ p ] . We set D α , σ i [ p ] = if Q i [ p ] < α .
First, recall category Ψ α and scale Λ α defined in (5). The theorem below shows the flight distance attained by proposed scheme σ , i.e., the length of a round-trip air route when the sensor field is covered by tessellation T Λ α ( Ψ α ) .
Theorem 4.
Suppose that we employ scheme  σ  to establish points of interest over the sensor field of radius α  and construct a round-trip air route that visits every point of interest. Then,
D α , σ = 3 β K Λ α ( Ψ α )
where  K Λ α ( Ψ α )  is the number of points of interest established over the sensor field of radius  α  according to scheme  σ .
Proof. 
A proof of Theorem 4 is given in Appendix E. □
Secondly, consider a set of hexagonal tessellations S i [ p ] : p 1 , , ν i , i N . Focus on scheme σ i [ p ] which employs tessellation S i [ p ] for pattern p 1 , , ν i and size i N . Given the radius of the sensor field α ( 0 , ) , let Ξ α and Δ α denote pattern and size that satisfy
D ( α , σ Δ α [ Ξ α ] ) = m i n D ( α , σ i p ) : Q i [ p ] α f o r   p 1 , , ν i   a n d   i N
Then, scheme σ Δ α [ Ξ α ] invokes the minimum flight distance when the radius of the sensor field is α . Henceforth, we call scheme σ Δ α [ Ξ α ] an optimal scheme in the sense of minimizing the flight distance. Also, let σ denote a scheme that uses scheme σ Δ α [ Ξ α ] for radius α ( 0 , ) collectively. Note that optimal scheme σ Δ α [ Ξ α ] is hard to find in general.
Consider set of hexagonal tessellations S i [ p ] : p { 1 , , ν i } for size i N . Recall that Q i [ p ] denotes the maximum radius of the sensor field that can be enclosed in tessellation S i [ p ] for p { 1 , , ν i } . Let Φ i denote a pattern such that
Q i [ Φ i ] = m a x Q i [ p ] : p { 1 , , ν i }
for size i N . Then, scheme σ i [ Φ i ] , which employs tessellation Q i [ Φ i ] , satisfies
N α , σ i [ Φ i ] = m i n N α , σ j [ p ] : p 1 , , ν i , j { i , i + 1 , }
for α ( Q i 1 Φ i 1 , Q i Φ i ] , i.e., scheme σ i [ Φ i ] invokes the minimum number of points of interest among the schemes which use tessellations that can enclose the sensor field of radius α ( Q i 1 Φ i 1 , Q i Φ i ] . However, it is not guaranteed that
D α , σ i [ Φ i ] = m i n D α , σ j [ p ] : p 1 , , ν i , j { i , i + 1 , }
for α ( Q i 1 Φ i 1 , Q i Φ i ] . Thus, scheme σ i [ Φ i ] is not necessarily an optimal scheme. Hereafter, we call σ i [ Φ i ] a suboptimal scheme in the sense of minimizing the number of points of interest. Also, let σ denote a scheme that uses suboptimal scheme σ i [ Φ i ] for radius α ( Q i 1 Φ i 1 , Q i Φ i ] and size i N collectively.
Theorem 5.
For radius  α ( Q i 1 Φ i 1 , Q i Φ i ]  and size  i N suboptimal scheme  σ i [ Φ i ]  is optimal scheme is identical to optimal scheme  σ Δ α [ Ξ α ]  if the graph induced by tessellation  S i [ Φ i ]  is Hamiltonian.
Proof. 
Tessellation S i [ Φ i ] consists of i tiles. Also, the straight-line distance between two adjacent points of interest is equal to 3 β . Thus, the flight distance of a round-trip air route must be greater than or equal to 3 β i . If the graph induced by tessellation S i [ Φ i ] is Hamiltonian, there is a Hamiltonian cycle, which consists of i links. If the length of each link is commonly 3 β , the length of the Hamiltonian cycle is 3 β i . Since the flight distance is equal to the length of the Hamiltonian cycle, suboptimal scheme σ i [ Φ i ] yields the minimum flight distance among the schemes that use tessellations that can enclose the sensor field of radius α ( Q i 1 Φ i 1 , Q i Φ i ] , i.e., suboptimal scheme σ i [ Φ i ] is identical to optimal scheme σ Δ α [ Ξ α ] . □
For size i N , we can find a suboptimal scheme σ i [ Φ i ] by completely searching all ν i patterns and then calculating maximum radius Q i [ p ] for each p { 1 , , ν i } . In general, however, it is hard to obtain suboptimal scheme σ i [ Φ i ] since number patterns ν i increases drastically as size i increases and hence it is hardly tractable to identify all the patterns. Only for relatively small i , we can discover all ν i patterns, derive maximum radius derive Q i [ p ] for each pattern p { 1 , , ν i } and then find suboptimal pattern Φ i and suboptimal scheme σ i [ Φ i ] as well. The table below lists suboptimal scheme σ i [ Φ i ] for some small i . Also, the table reveals maximum radius Q i [ Φ i ] and flight distance D α , σ i [ Φ i ] achieved by the adoption of scheme σ i [ Φ i ] at α ( Q i 1 Φ i 1 , Q i Φ i ] .
Tessellation S i [ Φ i ] , which is employed by suboptimal scheme σ i [ Φ i ] in Table 1, is illustrated in Appendix F. Note that the graph induced by tessellation S i [ Φ i ] , which is adopted by suboptimal scheme σ i [ Φ i ] listed in Table 1, is Hamiltonian. Thus, from Theorem 5, we confirm that suboptimal scheme σ i [ Φ i ] is also optimal scheme s σ Δ α [ Ξ α ] for radius ( Q i 1 Φ i 1 , Q i Φ i ] .
Thirdly, a lower bound on minimum flight distance m i n D α , σ i [ p ] : p 1 , , ν i , i N can be yielded as follows. Recall that each coverage takes the shape of a circular disk with radius β . Then, the apothem of a hexagonal tile inscribed in a coverage is 3 2 β and consequently its area is equal to 3 3 2 β 2 (See Figure 8). The area of the sensor field, which takes the shape of a circular disk with radius α , is equal to π α 2 . Note that the area of tessellation S i [ p ] , which is used by scheme σ i [ p ] , must be an integral multiple of the area of a hexagonal tile. Furthermore, the area of S i [ p ] should be greater than or equal to the area of the sensor field. Thus, the number of points of interest established by scheme σ i [ p ] is bounded below as follows.
N α , σ i [ p ] 2 3 π 9 α β 2
for all p { 1 , , ν i } and i N . Consider a hexagonal tessellation, which consists of 2 3 π 9 α β 2 tiles. Assume that the graph induced by the hexagonal tessellation is Hamiltonian. Then, the minimum flight distance is equal to the length of the Hamiltonian cycle multiplied by the distance between two adjacent points of interest, which is 3 β 2 3 π 9 α β 2 . Thus, the flight distance achieved by scheme σ i [ p ] is bounded below by
D α , σ i [ p ] 3 β 2 3 π 9 α β 2
for all p { 1 , , ν i } and i N .

5. Performance Evaluation of Proposed Scheme

In this section, we evaluate proposed scheme σ by comparing proposed scheme σ and optimal scheme σ in flight distance. Also, we compare the flight distance achieved by proposed scheme σ with the lower bound expressed in (13).
Figure 14 shows the maximum radius (normalized by β ) of the sensor field that can be attained by proposed scheme σ with respect to the number of points of interest established by σ . In addition, Figure 14 exhibits the maximum radius of sensor field that can be achieved by optimal scheme σ . When the number of points of interest is small, the difference between the maximum radii yielded by proposed scheme σ and optimal scheme σ seems to be insignificant.
Figure 15 shows the numbers of points of interest established by proposed scheme σ and optimal scheme σ with respect to the radius of the sensor field (normalized by β . Also, Figure 16 exhibits the flight distances (normalized by β ) yielded by proposed scheme σ and optimal scheme σ with respect to the radius (normalized by β ) of the sensor field. When the radius of the sensor field is small, the difference between the flight distances attained by proposed scheme σ and optimal scheme σ seems not to be significant. To quantitatively evaluate proposed scheme σ further, define the average difference in flight distance, denoted by Δ ( σ , σ , ω ) , as follows.
Δ σ , σ , ω = 1 ω 0 ω D α , σ D α , σ d α
for ω ( 0 , ) . Note that optimal scheme σ is available, roughly, at α [ 0,30 β ] . When ω = 30 β , average difference Δ σ , σ , ω is calculated to be 0.588 β , which is less than 60 % of the radius of a coverage. Moreover, the average difference in flight distance is expected to dwindle as observation time ω increases since optimal scheme σ in the sense of minimizing the number of points of interest does not minimize the flight distance more frequently as ω increases.
The table below summarizes specifications of a common multicopter [24,25,26].
Consider a wireless sensor network that employs such a multicopter as a sink node. Suppose that we cover the sensor field of the wireless sensor network with a hexagonal tessellation, establish a point of interest at the centroid of each tile, and construct an air route of the multicopter to visit every point of interest. From Table 2, we obtain the radius of a coverage β = 100   m × tan 97 ° / 2 113   m . Then, the distance between two adjacent points of interest is equal to 100   m × tan 97 ° / 2 × 3 196   m . Let κ denote the number of points of interest established over the sensor field so that the multicopter can finish the round-trip flight without recharging the battery. Then, number κ should necessarily satisfy the following inequality.
κ 720   s e c × 5   m / s e c 100   m × tan 97 ° / 2 × 3 κ 18 .
By employing proposed scheme σ , we can cover the sensor field with T 1 ( 3 ) as far as radius of the sensor field α 3 × β 339   m from (2). On the other hand, by use of optimal scheme σ , we can cover the sensor field with the hexagonal tessellation with 18 tiles, which is mentioned in Table 1 and depicted in Appendix F, if radius α is up to 37 / 2 × β 344   m . Note that a limit on the flight distance is calculated by multiplying the airborne time by the speed of the UAV in (15). In practice, the speed of a UAV is influenced by many factors, e.g., altitude at which the UAV flies and the direction of wind blowing to the UAV. The airborne time highly depends on the battery that operates the UAV. The battery duration is affected by the temperature. The vertical movement of a UAV consumes more energy than horizontal movement. The expression for the maximum number of points of interest in (15) does not include all the influencing factors mentioned above. Nonetheless, it is valid and useful as a rule of thumb.
Consider two cases; Case 1, in which 4 sensor nodes are sparsely scattered in an area as shown in Figure 17a, and Case 2, in which 36 sensor nodes are densely populated in the same as shown in Figure 18a. According to proposed scheme σ , tessellation T 2 ( 1 ) is employed to enclose the sensor field, which takes the shape of a circular disk in which all the sensor nodes are inscribed in both cases. Since tessellation T 2 ( 1 ) consists of 19 tiles. Nineteen points of interest are established over the sensor field. Note that graph G 2 ( 1 ) , which is induced by tessellation T 2 ( 1 ) , is Hamiltonian. Thus, from Theorem 4, the flight distance of the air route which is constructed along a Hamiltonian cycle in G 2 ( 1 ) is equal to 3 × 19 × β 32.9 × β in both cases. (See Figure 17b and Figure 18b). For a comparative evaluation of proposed scheme σ , we consider three other methods for establishing points of interest and constructing a round-trip air route that visit every point of interest; Method 1, which establishes a point of interest on each sensor node and constructs a round-trip air route by solving a traveling salesman problem; Method 2, which was introduced in [13]; and Method 3, which was mentioned in [14]. Figure 17c and Figure 18c illustrate the air routes constructed by Method 1 in sparsely and densely populated cases, respectively. From these figures, we obtain the flight distances of the air routes to be 14 2 × β 19.8 × β and 36 2 × β 50.9 × β , respectively, in Cases 1 and 2. Figure 17d and Figure 18d exhibit the air routes constructed by Method 2 in Cases 1 and 2, respectively. From these figures, we observe that two air routes are identical in both cases, and the flight distance of these air routes is equal to 32 ×   β . Figure 17e demonstrates the air route yielded by Method 3, which is a variant of Method 1. As shown in Figure 17e, the flight distance of the air route is equal to 10 2 × β 14.1 × β . From the results above, we can conclude that either Method 1 or Method 2 does not significantly outperform proposed scheme σ . In addition to the comparison of flight distances, we should consider the following factors in the comparative evaluation. First, to use Method 1 (or Method 4), points of interest should be bijectively established on (or near) sensor nodes. Thus, precise locational information about sensor nodes should be available. On the other hand, proposed scheme σ can be applied only with vague information about the geographic distribution of sensor nodes. Secondly, the UAV always receive data from a single sensor node when Method 1 or Method 3 is used. On the other hand, the UAV occasionally has to receive data from two or more sensor nodes if the proposed scheme is employed. Thus, a scheduling-based or contention-based medium access control scheme is needed, which brings about the consumption of additional time resources for signaling and throughput degradation. Thirdly, the area in which sensor nodes lie is divided into a number of regions in order to use Method 2. Then, a region takes the shape of a square in which a coverage is inscribed. Consequently, there can be a sensor node that is never able to deliver data to the UAV at whichever point of interest the UAV is. Fourthly, the distribution of sensor nodes can be changed as time goes on. As implied in Figure 17c and Figure 18c, Method 1, which places a point of interest at each sensor node, is highly sensitive to a change occurring in the number of sensor nodes and their locations. On the other hand, the proposed scheme is relatively insensitive to such a change, even though the throughput from sensor nodes to the UAV is affected to some extent.

6. Conclusions

In this paper, we considered a wireless sensor network that employs a UAV as a sink node. The wireless sensor network under discussion also lies in a harsh environment, so the geographical distribution of sensor nodes is barely known to us. In the network, the UAV flies along a round-trip air route that visits every point of interest. Whenever the UAV arrives at a point of interest, the UAV hovers and gathers data from neighboring sensor nodes on the ground. To satisfy the requirement that the UAV should be able to receive data from all the sensor nodes and to follow the recommendation that a shorter round-trip air route of the UAV is more preferred as well, we proposed a scheme characterized by covering the sensor field of the wireless sensor network with three classes of hexagonal tessellations, establishing a point of interest at the centroid of each tile belonging to the adopted tessellation, and constructing an air route of the UAV that visits every point of interest along a Hamiltonian cycle on the graph induced by the employed tessellation. Next, we developed a closed-form expression for the exact flight distance attained by the proposed scheme. For comparative evaluation, we discovered an optimal scheme that minimizes the flight distance for a given radius of the sensor field by completely inspecting patterns and corroborating the property of Hamiltonicity. Also, we presented a universal lower bound on the flight distance. Numerical examples showed that the flight distance attained by the proposed scheme is only slightly longer compared with the flight distance yielded by an optimal scheme. Furthermore, the proposed scheme was proven to be practically valid when a common multicopter is employed as a sink node.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflicts of interest.

Appendix A

For c { 1,2 , 3 } and n { 0 } N , tessellation T n + 1 ( c ) is yielded by letting minimum number of hexagonal tiles surround tessellation T n ( c ) . Let J n ( c ) denote such minimum number of tiles. Then, K n + 1 ( c ) = K n ( c ) + J n ( c ) for c { 1,2 , 3 } and n { 0 } N .
For c = 1 , we have
K 0 1 = 1 K 1 1 = K 0 1 + 6 · 1 = 7 K 2 1 = K 1 1 + 6 · 2     = 19 .
Suppose that
J n ( 1 ) = 6 n K n 1 = 3 n 2 + 3 n + 1
for n N . Then,
K n + 1 1 = K n 1 + J n ( 1 ) = 3 n 2 + 9 n + 1 = 3 ( n + 1 ) 2 + 3 n + 1 + 1
for n N , which proves Theorem 1 for c = 1 .
For c = 2 , we have
K 0 2 = 3 K 1 2 = K 0 2 + 3 · 3 = 12 K 2 2 = K 1 2 + 3 · 5 = 27 .
Suppose that
J n ( 2 ) = 3 ( 2 n + 3 ) K n 2 = 3 n 2 + 6 n + 3
for n N . Then,
K n + 1 3 = K n 3 + J n ( 3 ) = 3 n 2 + 15 n + 18 = 3 ( n + 1 ) 2 + 9 n + 1 + 6
for n N , which proves Theorem 1 for c = 2 .
For c = 3 , we have
K 0 3 = 6 K 1 3 = K 0 3 + 6 · 2 = 18 K 2 3 = K 1 3 + 6 · 3 = 36 .
Suppose that
J n ( 3 ) = 6 ( n + 2 ) K n 3 = 3 n 2 + 9 n + 6
for n N . Then,
K n + 1 3 = K n 3 + J n ( 3 ) = 3 n 2 + 15 n + 18 = 3 ( n + 1 ) 2 + 9 n + 1 + 6
for n N , which proves Theorem 1 for c = 3 .

Appendix B

Consider the tiles that form the boundary of tessellation T n ( c ) for c { 1,2 , 3 } and n N . Among the vertices of these tiles, focus on the vertices that lie on the boundary of T n ( c ) . Let d ( v ) denote the distance between the center point of tessellation T n ( c ) , denoted by o , and vertex v on the boundary of T n ( c ) . Let v denote a vertex such that d ( v ) d ( v ) for any vertex v on the boundary of T n ( c ) . Then, the radius of the inscribed circular disk in T n ( c ) , denoted by R n ( c ) , is equal to d ( v ) .
In T n ( 1 ) , vertex v is located with respect to center point o of T n ( 1 ) as shown in Figure A1.
Figure A1. Node v in tessellation T n ( 1 ) . (a) n = 2 m 1 for m N ; (b) n = 2 m for m N .
Figure A1. Node v in tessellation T n ( 1 ) . (a) n = 2 m 1 for m N ; (b) n = 2 m for m N .
Sensors 24 03867 g0a1
From Figure A1, we obtain
R 2 m 1 ( 1 ) = 2 α + 3 m 1 α = ( 3 m 1 ) · α
when n = 2 m 1 for m N , i.e., n is an odd number and
R 2 m ( 1 ) = [ α + 3 m 1 α + 5 2 α ] 2 + ( 3 2 α ) 2 = 9 m 2 + 3 m + 1 · α
when n = 2 m for m N , i.e., n is an even number.
In T n ( 2 ) , vertex v is located with respect to center point o of T n ( 2 ) as shown in Figure A2.
Figure A2. Node v in tessellation T n ( 2 ) . (a) n = 2 m 1 for m N ; (b) n = 2 m for m N .
Figure A2. Node v in tessellation T n ( 2 ) . (a) n = 2 m 1 for m N ; (b) n = 2 m for m N .
Sensors 24 03867 g0a2
From Figure A2, we obtain
R 2 m 1 ( 2 ) = [ α + 3 m 1 α + 3 2 α ] 2 + ( 3 2 α ) 2 = 9 m 2 3 m + 1 · α
when n = 2 m 1 for m N , i.e., n is an odd number and
R 2 m ( 2 ) = α + 3 m α = ( 3 m + 1 ) · α
when n = 2 m for m N , i.e., n is an even number.
In T n ( 3 ) , vertex v is located with respect to center point o of T n ( 3 ) as shown in Figure A3.
Figure A3. Node v in tessellation T n ( 3 ) . (a) n = 2 m 1 for m N ; (b) n = 2 m for m N .
Figure A3. Node v in tessellation T n ( 3 ) . (a) n = 2 m 1 for m N ; (b) n = 2 m for m N .
Sensors 24 03867 g0a3aSensors 24 03867 g0a3b
From Figure A3, we obtain
R 2 m 1 ( 3 ) = 2 α + 3 m 1 α + α = 3 m · α
when n = 2 m 1 for m N , i.e., n is an odd number and
R 2 m ( 3 ) = [ 3 α + 3 ( m 1 ) α + 3 2 α ] 2 + ( 3 2 α ) 2 = 9 m 2 + 9 m + 3 · α
when n = 2 m for m N , i.e., n is an even number.

Appendix C

We have
K n ( 1 ) = 3 n 2 + 3 n + 1 3 n 2 + 6 n + 3 = K n ( 2 ) K n ( 2 ) = 3 n 2 + 6 n + 3 3 n 2 + 9 n + 6 = K n 3
for n { 0 } N Thus,
K n ( 1 ) K n ( 2 ) K n ( 3 )
for all n { 0 } N .
Suppose that n = 2 m 1 for m N . Then, we have
R 2 m 1 ( 1 ) = 3 m 1 = 9 m 2 6 m + 1 9 m 2 3 m + 1 = R 2 m 1 ( 2 ) R 2 m 1 ( 2 ) = 9 m 2 3 m + 1 9 m 2 = 3 m = R 2 m 1 ( 3 )
for m N . Also, we have
R 2 m ( 1 ) = 9 m 2 + 3 m + 1 9 m 2 + 6 m + 1 = 3 m + 1 = R 2 m ( 2 ) R 2 m ( 2 ) = 3 m + 1 = 9 m 2 + 6 m + 1 9 m 2 + 9 m + 3 = R 2 m ( 3 )
for m N . Thus,
R n ( 1 ) R n ( 2 ) R n ( 3 )
for all n N .

Appendix D

We will prove that the graph induced by tessellation T n ( c ) , denoted by G n ( c ) , is Hamiltonian by showing that induced graph G n ( c ) is a triangular grid graph that is connected and locally connected [23].
First, graph G n ( c ) is a triangular grid graph since G n ( c ) is induced by a hexagonal tessellation. Secondly, graph G n ( c ) is apparently connected, i.e., there is a path between any pair of nodes in G n ( c ) . Thirdly, let V n ( c ) denote the set of all the nodes in G n ( c ) . For U V n ( c ) , let H n c ( U ) denote the maximal subgraph of G n ( c ) . For v V n ( c ) , let A n ( c ) ( v ) denote the set of nodes adjacent to v , i.e., nodes to which there are links from v . For v G n ( c ) , the degree of v , i.e., the number of links incident on v , is equal to one of 3 , 4 and 6 . For v 1 , v 2 G n ( c ) , H n c ( { v 1 } A n ( c ) ( v 1 ) ) is isomorphic to H n c ( { v 2 } A n ( c ) ( v 2 ) ) as far as the degree of v 1 is equal to the degree of v 2 . Also, for v 1 , v 2 G n ( c ) , H n c ( A n ( c ) ( v 1 ) ) is isomorphic to H n c ( A n ( c ) ( v 2 ) ) as far as the degree of v 1 is equal to the degree of v 2 . Figure A4, Figure A5 and Figure A6 illustrate H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is equal to 3 , 4 , and 6 , respectively. For all v V n ( c ) , graph H n c ( A n ( c ) ( v ) ) is connected whatever the degree of node v is. Thus, graph G n ( c ) is locally connected. Therefore, graph G n ( c ) is Hamiltonian for c { 1,2 , 3 } and n N .
Figure A4. Maximal subgraphs H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is 3 . (a) H n c v A n c v ; (b) H n c ( A n ( c ) ( v ) ) .
Figure A4. Maximal subgraphs H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is 3 . (a) H n c v A n c v ; (b) H n c ( A n ( c ) ( v ) ) .
Sensors 24 03867 g0a4
Figure A5. Maximal subgraphs H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is 4 . (a) H n c v A n c v ; (b) H n c ( A n ( c ) ( v ) ) .
Figure A5. Maximal subgraphs H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is 4 . (a) H n c v A n c v ; (b) H n c ( A n ( c ) ( v ) ) .
Sensors 24 03867 g0a5
Figure A6. Maximal subgraphs H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is 6 . (a) H n c v A n c v ; (b) H n c ( A n ( c ) ( v ) ) .
Figure A6. Maximal subgraphs H n c ( { v } A n ( c ) ( v ) ) and H n c ( A n ( c ) ( v ) ) when the degree of node v is 6 . (a) H n c v A n c v ; (b) H n c ( A n ( c ) ( v ) ) .
Sensors 24 03867 g0a6

Appendix E

Consider graph G n ( c ) , which is induced by tessellation T n ( c ) . Construct a supergraph of G n ( c ) , denoted by G ~ n ( c ) , to be the complete graph built on set of nodes 1 , , K n ( c ) . Assign the value of the straight-line distance between two distinct points of interest to the link between the corresponding nodes. Then, a shortest round-trip airway that visits all the points of interest coincides with a Hamiltonian cycle on G ~ n ( c ) , if any. A Hamiltonian cycle on G ~ n ( c ) should contain only K n ( c ) links. Note that the length of a link is at least 3 β . Thus, the length of a Hamiltonian cycle on G ~ n ( c ) , if exists, is greater than or equal to 3 β K n ( c ) . Note that induced graph G n ( c ) is Hamiltonian, and the length of a Hamiltonian cycle on G n ( c ) , in which 3 β is assigned to each link, is equal to 3 β K n ( c ) . Since G n ( c ) is a spanning subgraph of G ~ n ( c ) , there exists a Hamiltonian cycle on G ~ n ( c ) whose length is equal to 3 β K n ( c ) . Therefore, the length of the shortest round-trip airway that visits every point of interest established over the sensor field when the sensor field is covered by tessellation T n ( c ) is equal to 3 β K n ( c ) .

Appendix F

Figure A7. Tessellation S i [ Φ i ] which is employed by scheme σ i [ Φ i ] in Table 1. (a) i = 1 ; (b) i = 3 ; (c) i = 4 ; (d) i = 5 ; (e) i = 6 ; (f) i = 7 ; (g) i = 9 ; (h) i = 10 ; (i) i = 12 ; (j) i = 14 ; (k) i = 16 ; (l) i = 18 .
Figure A7. Tessellation S i [ Φ i ] which is employed by scheme σ i [ Φ i ] in Table 1. (a) i = 1 ; (b) i = 3 ; (c) i = 4 ; (d) i = 5 ; (e) i = 6 ; (f) i = 7 ; (g) i = 9 ; (h) i = 10 ; (i) i = 12 ; (j) i = 14 ; (k) i = 16 ; (l) i = 18 .
Sensors 24 03867 g0a7

References

  1. Seo, H.; Park, J.K.; Choi, C.W. Timeliness-aware Anti-discrimination MAC Scheme for Wireless Passive Sensor Networks. IEIE Trans. Smart Process. Comput. 2020, 9, 151–168. [Google Scholar] [CrossRef]
  2. Seino, W.; Yoshihisa, T.; Hara, T.; Nishio, S. A Communication Protocol to Improve Fairness and Data Amount on Sensor Data Collection with a Mobile Sink. In Proceedings of the 2010 IEEE International Conference on Broadband, Wireless Computing, Communication, and Applications (BWCCA), Fukuoka, Japan, 4–6 November 2010; pp. 33–40. [Google Scholar] [CrossRef]
  3. Kwong, K.; Goh, H.; Stephen, B.; Wu, T.; Sasloglou, K.; Shen, C.; Michie, C.; Glover, I.; Andonovic, I.; Du, W. Implementation of herd management systems with wireless sensor networks. IET Wirel. Sens. Syst. 2011, 1, 55–65. [Google Scholar] [CrossRef]
  4. Chang, J.-Y.; Jeng, J.-T.; Sheu, Y.-H.; Jian, Z.-J.; Chang, W.-Y. An efficient data collection path planning scheme for wireless sensor networks with mobile sinks. EURASIP J. Wirel. Commun. Netw. 2020, 2020, 257. [Google Scholar] [CrossRef]
  5. Ho, D.-T.; Grøtli, E.I.; Sujit, P.B.; Johansen, T.A.; Sousa, J.B. Optimization of Wireless Sensor Network and UAV Data Acquisition. J. Intell. Robot. Syst. 2015, 78, 159–179. [Google Scholar] [CrossRef]
  6. Olivieri, B.; Endler, M. DADCA: An Efficient Distributed Algorithm for Aerial Data Collection from Wireless Sensors Networks by UAVs. In Proceedings of the 20th ACM International Conference on Modelling, Analysis and Simulation of Wireless and Mobile Systems (MSWIM), Miami, FL, USA, 21–25 November 2017; pp. 129–136. [Google Scholar] [CrossRef]
  7. Salarian, H.; Chin, K.-W.; Naghdy, F. An Energy-Efficient Mobile-Sink Path Selection Strategy for Wireless Sensor Networks. IEEE Trans. Veh. Technol. 2013, 63, 2407–2419. [Google Scholar] [CrossRef]
  8. Khan, A.W.; Abdullah, A.H.; Razzaque, M.A.; Bangash, J.I. VGDRA: A Virtual Grid-Based Dynamic Routes Adjustment Scheme for Mobile Sink-Based Wireless Sensor Networks. IEEE Sensors J. 2014, 15, 526–534. [Google Scholar] [CrossRef]
  9. Da Silva, R.; Nascimento, M. On Best Drone Tour Plans for Data Collection in Wireless Sensor Network. In Proceedings of the 31st Annual ACM Symposium on Applied Computing (SAC), Pisa, Italy, 4–8 April 2016; pp. 703–708. [Google Scholar] [CrossRef]
  10. Cao, H.; Yang, Z.; Li, Y.-Q. A Mobile WSN Sink Node Using Unmanned Aerial Vehicles: Design And Experiment. Int. J. Interact. Mob. Technol. 2016, 10, 64–67. [Google Scholar] [CrossRef]
  11. Wu, C.; Liu, Y.; Wu, F.; Fan, W.; Tang, B. Graph-Based Data Gathering Scheme in WSNs With a Mobility-Constrained Mobile Sink. IEEE Access 2017, 5, 19463–19477. [Google Scholar] [CrossRef]
  12. Wen, W.; Zhao, S.; Shang, C.; Chang, C.-Y. EAPC: Energy-Aware Path Construction for Data Collection Using Mobile Sink in Wireless Sensor Networks. IEEE Sensors J. 2018, 18, 890–901. [Google Scholar] [CrossRef]
  13. Kim, E.-J.; Choi, H.H.; Kwon, J.-H. Regional Density-aware Data Collection Using Unmanned Aerial Vehicle in Large-scale Wireless Sensor Networks. Sensors Mater. 2018, 30, 1735–1742. [Google Scholar] [CrossRef]
  14. Chao, F.; He, Z.; Pang, A.; Zhou, H.; Ge, J. Path Optimization of Mobile Sink Node in Wireless Sensor Network Water Monitoring System. Complexity 2019, 2019, 5781620. [Google Scholar] [CrossRef]
  15. Vera-Amaro, R.; Rivero-Ángeles, M.E.; Luviano-Juárez, A. Data Collection Schemes for Animal Monitoring Using WSNs-Assisted by UAVs: WSNs-Oriented or UAV-Oriented. Sensors 2020, 20, 262. [Google Scholar] [CrossRef] [PubMed]
  16. Zhang, F.; Wang, G.; Hu, Y.; Chen, L.; Zhu, A.-X. Design of an Integrated Remote and Ground Sensing Monitor System for Assessing Farmland Quality. Sensors 2020, 20, 336. [Google Scholar] [CrossRef] [PubMed]
  17. Park, J.; Kim, S.; Youn, J.; Ahn, S.; Cho, S. Low-Complexity Data Collection Scheme for UAV Sink Nodes in Cellular IoT Networks. IEEE Trans. Veh. Technol. 2021, 70, 4865–4879. [Google Scholar] [CrossRef]
  18. Zhu, Y.; Wang, S. Efficient Aerial Data Collection with Cooperative Trajectory Planning for Large-Scale Wireless Sensor Networks. IEEE Trans. Commun. 2021, 70, 433–444. [Google Scholar] [CrossRef]
  19. Zhao, L.; Li, S.; Guan, Y.; Wan, S.; Hawbani, A.; Bi, Y.; Guizani, M. Adaptive Multi-UAV Trajectory Planning Leveraging Digital Twin Technology for Urban IIoT Applications. IEEE Trans. Netw. Sci. Eng. 2023, 1–16. [Google Scholar] [CrossRef]
  20. Grünbaum, B.; Shephard, G. Tiles and Patterns; W. H. Freeman and Company: New York, NY, USA, 1987; ISBN -10 0486469816. [Google Scholar]
  21. Harary, F. Graph Theory; Addison-Wesley: Bosston, MA, USA, 1969; ISBN -10 0201410338. [Google Scholar]
  22. Sugimoto, T. Convex Pentagons for Edge-to-Edge Tiling, I. Forma 2012, 27, 93–103. [Google Scholar]
  23. Gordon, V.S.; Orlovich, Y.L.; Werner, F. Hamiltonian properties of triangular grid graphs. Discret. Math. 2008, 308, 6166–6188. [Google Scholar] [CrossRef]
  24. Radiansyah, S.; Kusrini, M.D.; Prasetyo, L.B. Quadcopter applications for wildlife monitoring. In IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2017; Volume 54, p. 012066. [Google Scholar] [CrossRef]
  25. Anweiler, S.; Piwowarski, D. Multicopter Platform Prototype for Environmental Monitoring. J. Clean. Prod. 2017, 155, 204–211. [Google Scholar] [CrossRef]
  26. Sidibe, A.; Loubet, G.; Takacs, A.; Ferré, G.; Ghiotto, A. Miniature drone antenna design for the detection of airliners. Int. J. Microw. Wirel. Technol. 2020, 13, 21–27. [Google Scholar] [CrossRef]
Figure 1. Sensor field of wireless sensor network which is in the shape of circular disk.
Figure 1. Sensor field of wireless sensor network which is in the shape of circular disk.
Sensors 24 03867 g001
Figure 2. Points of interest established over sensor field of wireless sensor network.
Figure 2. Points of interest established over sensor field of wireless sensor network.
Sensors 24 03867 g002
Figure 3. Coverage associated with point of interest which is in the shape of circular disk.
Figure 3. Coverage associated with point of interest which is in the shape of circular disk.
Sensors 24 03867 g003
Figure 4. Tessellation consisting of tiles in two polygonal shapes.
Figure 4. Tessellation consisting of tiles in two polygonal shapes.
Sensors 24 03867 g004
Figure 5. Monohedral tessellation consisting of triangular tiles.
Figure 5. Monohedral tessellation consisting of triangular tiles.
Sensors 24 03867 g005
Figure 6. Non-edge-to-edge monohedral tessellation consisting of rectangular tiles.
Figure 6. Non-edge-to-edge monohedral tessellation consisting of rectangular tiles.
Sensors 24 03867 g006
Figure 7. Edge-to-edge monohedral tessellation comprised by tiles that take the shape of a concave polygon.
Figure 7. Edge-to-edge monohedral tessellation comprised by tiles that take the shape of a concave polygon.
Sensors 24 03867 g007
Figure 8. Hexagonal tile inscribed in coverage taking shape of circular disk with radius β .
Figure 8. Hexagonal tile inscribed in coverage taking shape of circular disk with radius β .
Sensors 24 03867 g008
Figure 9. Three categories of hexagonal tessellations. (a) Tessellation T 2 ( 1 ) ; (b) Tessellation T 2 ( 2 ) ; (c) Tessellation T 2 ( 3 ) .
Figure 9. Three categories of hexagonal tessellations. (a) Tessellation T 2 ( 1 ) ; (b) Tessellation T 2 ( 2 ) ; (c) Tessellation T 2 ( 3 ) .
Sensors 24 03867 g009aSensors 24 03867 g009b
Figure 10. Tessellations T 0 ( 1 ) , T 0 ( 2 ) and T 0 ( 3 ) . (a) Tessellation T 0 ( 1 ) ; (b) Tessellation T 0 ( 2 ) ; (c) Tessellation T 0 ( 3 ) .
Figure 10. Tessellations T 0 ( 1 ) , T 0 ( 2 ) and T 0 ( 3 ) . (a) Tessellation T 0 ( 1 ) ; (b) Tessellation T 0 ( 2 ) ; (c) Tessellation T 0 ( 3 ) .
Sensors 24 03867 g010
Figure 11. Graphs G 2 ( 1 ) , G 2 ( 2 ) and G 2 ( 3 ) induced by tessellations T 2 ( 1 ) , T 2 ( 2 ) and T 2 ( 3 ) . (a) Graph G 2 ( 1 ) induced by tessellation T 2 ( 1 ) ; (b) Graph G 2 ( 2 ) induced by tessellation T 2 ( 2 ) ; (c) Graph G 2 ( 3 ) induced by tessellation T 2 ( 3 ) .
Figure 11. Graphs G 2 ( 1 ) , G 2 ( 2 ) and G 2 ( 3 ) induced by tessellations T 2 ( 1 ) , T 2 ( 2 ) and T 2 ( 3 ) . (a) Graph G 2 ( 1 ) induced by tessellation T 2 ( 1 ) ; (b) Graph G 2 ( 2 ) induced by tessellation T 2 ( 2 ) ; (c) Graph G 2 ( 3 ) induced by tessellation T 2 ( 3 ) .
Sensors 24 03867 g011
Figure 12. Hamiltonian cycles in graphs G 2 ( 1 ) , G 2 ( 2 ) and G 2 ( 3 ) . (a) Hamiltonian cycle in graph G 2 ( 1 ) ; (b) Hamiltonian cycle in graph G 2 ( 2 ) ; (c) Hamiltonian cycle in graph G 2 ( 3 ) .
Figure 12. Hamiltonian cycles in graphs G 2 ( 1 ) , G 2 ( 2 ) and G 2 ( 3 ) . (a) Hamiltonian cycle in graph G 2 ( 1 ) ; (b) Hamiltonian cycle in graph G 2 ( 2 ) ; (c) Hamiltonian cycle in graph G 2 ( 3 ) .
Sensors 24 03867 g012
Figure 13. Exemplary establishment of points of interest and construction of round-trip air route according to proposed scheme σ . (a) Enclosement of sensor filed with tessellation T 0 ( 3 ) and establishment of points of interest denoted by P 1 , , P 6 ; (b) inducement of graph G 0 ( 3 ) from tessellation T 0 ( 3 ) ; (c) discovery of Hamiltonian cycle v 1 , v 2 , v 4 , v 5 , v 6 , v 3 , v 1 in graph G 0 ( 3 ) ; (d) construction of air route P 1 , P 2 , P 4 , P 5 , P 6 , P 3 , P 1 over sensor field.
Figure 13. Exemplary establishment of points of interest and construction of round-trip air route according to proposed scheme σ . (a) Enclosement of sensor filed with tessellation T 0 ( 3 ) and establishment of points of interest denoted by P 1 , , P 6 ; (b) inducement of graph G 0 ( 3 ) from tessellation T 0 ( 3 ) ; (c) discovery of Hamiltonian cycle v 1 , v 2 , v 4 , v 5 , v 6 , v 3 , v 1 in graph G 0 ( 3 ) ; (d) construction of air route P 1 , P 2 , P 4 , P 5 , P 6 , P 3 , P 1 over sensor field.
Sensors 24 03867 g013
Figure 14. Maximum radii (normalized by β ) of sensor field yielded by proposed scheme σ and optimal scheme σ with respect to number of points of interest established over sensor field.
Figure 14. Maximum radii (normalized by β ) of sensor field yielded by proposed scheme σ and optimal scheme σ with respect to number of points of interest established over sensor field.
Sensors 24 03867 g014
Figure 15. Numbers of points of interest established by proposed scheme σ and optimal scheme σ with respect to the radius of the sensor field.
Figure 15. Numbers of points of interest established by proposed scheme σ and optimal scheme σ with respect to the radius of the sensor field.
Sensors 24 03867 g015
Figure 16. Flight distances (normalized by β ) yielded by proposed scheme σ and optimal scheme σ with respect to the radius (normalized by β ) of the sensor field.
Figure 16. Flight distances (normalized by β ) yielded by proposed scheme σ and optimal scheme σ with respect to the radius (normalized by β ) of the sensor field.
Sensors 24 03867 g016
Figure 17. Comparison of flight distances of air routes constructed by proposed scheme σ and Methods 1, 2, and 3. (a) Sparsely scattered sensor nodes on the ground; (b) air route constructed by proposed scheme σ ; (c) air route constructed by Method 1; (d) air route constructed by Method 2; (e) air route constructed by Method 3.
Figure 17. Comparison of flight distances of air routes constructed by proposed scheme σ and Methods 1, 2, and 3. (a) Sparsely scattered sensor nodes on the ground; (b) air route constructed by proposed scheme σ ; (c) air route constructed by Method 1; (d) air route constructed by Method 2; (e) air route constructed by Method 3.
Sensors 24 03867 g017aSensors 24 03867 g017b
Figure 18. Comparison of flight distances of air routes constructed by proposed scheme σ and Methods 1, 2, and 3. (a) Densely populated sensor nodes on the ground; (b) air route constructed by proposed scheme σ ; (c) air route constructed by Method 1; (d) air route constructed by Method 2.
Figure 18. Comparison of flight distances of air routes constructed by proposed scheme σ and Methods 1, 2, and 3. (a) Densely populated sensor nodes on the ground; (b) air route constructed by proposed scheme σ ; (c) air route constructed by Method 1; (d) air route constructed by Method 2.
Sensors 24 03867 g018aSensors 24 03867 g018b
Table 1. Scheme σ i [ Φ i ] , maximum radius Q i [ Φ i ] , and flight distance D α , σ i [ Φ i ] .
Table 1. Scheme σ i [ Φ i ] , maximum radius Q i [ Φ i ] , and flight distance D α , σ i [ Φ i ] .
Scheme
σ i [ Φ i ]
Maximum Radius
Q i [ Φ i ]
Is   Graph   Induced   by   S i [ Φ i ] Hamiltonian?Flight Distance
D α , σ i [ Φ i ]
σ 1 [ Φ 1 ] 3 / 2 · β Yes 3 · β
σ 3 [ Φ 3 ] β Yes 3 3 · β
σ 4 [ Φ 4 ] 7 / 2 · β Yes 4 3 · β
σ 5 [ Φ 5 ] ( 6 3 9 ) · β Yes 5 3 · β
σ 6 [ Φ 6 ] 3 · β Yes 6 3 · β
σ 7 [ Φ 7 ] 2 · β Yes 7 3 · β
σ 9 [ Φ 9 ] 43 3 / 36 · β Yes 9 3 · β
σ 10 [ Φ 10 ] 19 / 2 · β Yes 10 3 · β
σ 12 [ Φ 12 ] 7 · β Yes 12 3 · β
σ 14 [ Φ 14 ] 31 / 2 · β Yes 14 3 · β
σ 16 [ Φ 16 ] 31 / 11 · β Yes 16 3 · β
σ 18 [ Φ 18 ] 37 / 2 · β Yes 18 3 · β
Table 2. Specifications of multicopter.
Table 2. Specifications of multicopter.
ParameterUnitValue
AltitudeMeters100
SpeedMeters/second5
Beamwidth of receiving antennaDegrees97
Airborne timeSeconds720
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Choi, C. Tessellation-Based Construction of Air Route for Wireless Sensor Networks Employing UAV. Sensors 2024, 24, 3867. https://doi.org/10.3390/s24123867

AMA Style

Choi C. Tessellation-Based Construction of Air Route for Wireless Sensor Networks Employing UAV. Sensors. 2024; 24(12):3867. https://doi.org/10.3390/s24123867

Chicago/Turabian Style

Choi, CheonWon. 2024. "Tessellation-Based Construction of Air Route for Wireless Sensor Networks Employing UAV" Sensors 24, no. 12: 3867. https://doi.org/10.3390/s24123867

APA Style

Choi, C. (2024). Tessellation-Based Construction of Air Route for Wireless Sensor Networks Employing UAV. Sensors, 24(12), 3867. https://doi.org/10.3390/s24123867

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