3.1. Example: Sensing along the Line
Consider a one-dimensional problem where sensors are to be placed along a line in a domain . Specifically, take to be the interval . This line segment may represent a barrier in some surveillance application, or it may represent a region along which practical restrictions require all sensors to be on the same line. For whatever the practical application, for sensor placement purposes the only concern is that the sensors are to be placed in to best achieve a desired level of probabilistic coverage .
For a first example,
sensors are to be placed in
to meet a desired coverage of
Thus, there is a desire for larger coverage in a region in the middle (where
) with a lesser need for coverage outside of that region, constituting a situation with a higher priority region near the middle of the domain. For this example, let us further consider the environmental characteristics in the region to be homogeneous, with
and
. Thus, any sensor placed in the domain will observe the object of interest with probability 0.5 if the object is located within one unit of the sensor position. Applying these conditions to Equation (
10) yields a sensor density function
as shown in
Figure 1. The desired sensor density
for this case follows the shape of the desired coverage, which is expected for a uniform environment. While the shapes of
and
are the same in this case, the specific levels of
differ from those of
as they are based on not only the desired coverage
, but also the number of available sensors
and the sensor performance characteristics in the environment (
and
).
Running a sampling procedure for
sensors with the density
as shown in
Figure 1 leads to the sensor positions shown in
Figure 2, where the specific resulting sensor locations are given by the circles along the x-axis. Also shown in
Figure 2 is the resulting probabilistic coverage from the sensors, as given by
in Equation (
2), as well as the desired coverage
(for comparison). Note that the coverage obtained through sampling has a similar shape to the desired coverage; although it is generally larger because there are more sensors than required to achieve the desired coverage. The sampling procedure helps to maintain the shape in such situations. The quality of the match is clearly limited by the physical constraints on the environment and number of sensors available.
To show the quality of the solutions obtained relative to the level of match that is achievable by this specific number of sensors with these specific environmental characteristics, a direct optimization of the coverage function for
given in Equation (
2) was also performed. As this function is made up of many segments of constant levels for
, it is not differentiable and thus not amenable to gradient-based optimization approaches. The numerical optimization procedure utilized in the this paper employs a genetic algorithm metaheuristic, using a standard genetic algorithm [
26] with single-point crossover, roulette selection, and an elitist selection strategy (maintaining the top two individuals from each generation). The parameters used were a population size of 50, a mutation probability of 1/64 (corresponding to one bit of mutation per individual in the population for each generation, on average), and each sensor’s x-location
was represented with an 8 bit binary string. The algorithm was run for 1000 generations or to convergence if it converged earlier. The results of the optimized positions for this example are shown in
Figure 3, where the sensor locations are shown along with the resulting coverage. Note that the locations and resulting coverage of the sampled solution in
Figure 2 are similar to the optimized result in
Figure 3. The benefit of the sampling approach over the numerical procedure is that the sampling approach is an analytic process that can provide solutions much more rapidly than the optimization approach, while still achieving many of the features from the optimal positioning of the sensors. In
Table 1, the results of the sampling procedure are shown for the one-dimensional homogeneous environment with varying numbers of sensors
. The results shown are the quality of the match with the desired coverage function, specifically measured by the
norm as
. For comparison purposes, also included in the table are the quality of match that is optimally achievable for each number of sensors
. Note that both the sampled and optimized results show similar trends in that the match deteriorates for both very small and very large numbers of sensors, as expected.
Example 1 was a problem with uniform environmental characteristics, which is not practical. For a more realistic situation, example 2 considers the placement of
sensors for the same desired coverage
that was given in Equation (
12). However, in example 2 the environment is non-homogeneous, leading to a spatial dependency for both the sensor ranges
as well as their probabilistic performance
. These dependencies are shown in
Figure 4, and the resulting desired sensor density function
from Equation (
10) is shown in
Figure 5. Comparing
Figure 5 to
Figure 1 shows how the environmental effects have a great impact on the desired sensor density. In particular, more sensors are desired at the right side (near
) than the left side (near
) since the detection performance
for individual sensors is lower there, requiring more overlap of coverage to achieve the goal. Also, there is a “bump” in the middle of the desired coverage region (from
) that did not exist in
Figure 1, owing to the fact that the environmental characteristics can have as great an effect on where to place sensors as the desired coverage trends.
Figure 6 and
Figure 7 show the resulting placements of
sensors for example 2 that were obtained using the analytic sampling approach and the optimization approach, respectively. As in the case with environmental homogeneity, for this case the sensor locations and the resulting coverage performance are qualitatively similar. The major difference between them is a sensor to the far left (near
) in the optimized result that is not in the sampled result. This is because the relatively small number of sensors makes the sampling approach somewhat inefficient in portions of the domain with low sensor density (small
). However, that limited coverage is always going to be in a portion of the domain with lower coverage. The next example shows that this effect is not as prominent as the number of sensors increases.
For a third example, consider the same
and
of example 2 (as shown in
Figure 4) as well as the same desired coverage
from Equation (
12). However, now
sensors are placed in the region. In this dense sampling regime, the shape of the desired sensor density function
is identical to that shown in
Figure 5 for example 2, it only differs by a scaling factor of
due to the
term in Equation (
10). The resulting sampled placements and the optimized placements for example 3 are shown in
Figure 8 and
Figure 9, respectively. From this example, it is shown that the large number of sensors leads to coverage well above the desired coverage levels. The sampling approach provides a scaled version of the sampling from example 2, packing more sensors into the area of higher desired coverage and spreading out the remainder accordingly. However, the optimization approach now tries to directly match the levels of the desired coverage, leading to lower coverage in some portions of the high-coverage region (from
). In this sense, the sampling approach, in addition to being computationally much quicker than the optimization approach, may also provide solutions that are more desirable to the user (while not necessarily optimal in the
sense).
Table 2 shows the quality of the match of the results of the sampling procedure for the one-dimensional non-homogeneous environment with varying numbers of sensors
. The results are qualitatively similar to those seen in the homogeneous case in that the resulting match becomes worse for both very small and very large numbers of sensors
.
3.2. Example: Sensing in the Plane
For two-dimensional sensing, consider the situation where sensors are to be placed within a closed region
. Specifically, take
to be the square region
for these examples. Such a domain may represent an area that is to be monitored or measured for some unusual activity or concentration. The sensors under consideration are described by a range
and a probability
such that a sensor located at
will cover the disc of radius
that is centered at
with a probability
. For any practical application in such a domain, for sensor placement purposes the only concern is that the sensors are to be placed in
to best achieve some pre-defined desired level of probabilistic coverage
. For the examples that follow, the two-dimensional coverage goal is defined as follows:
as shown in
Figure 10. Note that this case has a desired nominal coverage level of 0.5 throughout most of the domain
, with a larger coverage level of 0.9 in a disc around the center, corresponding to a region of larger importance.
As a first example in this two-dimensional situation, consider a scenario where the environment is homogeneous with
and
. For this homogeneous example, the goal is to determine the best locations
to place
sensors to provide a probabilistic coverage
to best match the goal coverage
. Note that a coverage range of
implies an individual sensor coverage of area of 0.0314, which is 1/32 of the total area of
. Thus, with
total such sensors, there is not even an opportunity to cover the entire domain to the lower goal coverage level of 0.5. The question for sensor placement is to determine how much focus to put on overlapping sensors in the middle of
to achieve the desired higher coverage there versus spreading out sensors in the remainder of the domain to achieve the desired lower coverage there. Both the sampling-based placement and the optimization placement strategy provide this determination as part of their solutions. For the sampling approach, the sensor density
of Equation (
10) is computed to form the sampling distribution
as in Equation (
11). This associated CDF for the PDF
is computed and then sampled deterministically to find the sensor locations
. The resulting sampled sensor locations for this example are shown in
Figure 11, where the corresponding resulting coverage
(as computed from Equation (
2)) is also shown. From the resulting coverage, it is clear that the sampling approach provides a balance between spreading out some sensors to achieve the lower desired uniform coverage of 0.5 outside of the center of
while allowing some amount of overlap to achieve some of the larger desired coverage in the center of
.
The optimization approach for this homogeneous two-dimensional example uses the same genetic algorithm approach that was used for the one-dimensional example, where the objective function is given as in Equation (
3). In particular, the parameters
are each represented by a 16-bit binary string (8 bits for each of the two dimensions). The standard genetic algorithm [
26] is again utilized with single-point crossover, roulette selection and an elitist selection strategy. The other parameters used were a population size of 50 and a mutation rate of 1/320 (corresponding to one bit of mutation per individual in the population for each generation, on average). As in the one-dimensional examples, the genetic algorithm was run for 1000 generations, or to convergence if it converged earlier. The results of the optimization approach for the homogeneous two-dimensional example are shown in
Figure 12. As in the sampling-based result of
Figure 11, the optimization solution creates a balance between spreading out some sensors away from the center of
while allowing other sensors to provide the desired overlap in the center of
. Qualitatively, the relative split between these aspects is similar between the optimization and sampling-based approaches, and thus they provide similar levels of approximation to the desired coverage goal
that was shown in
Figure 10.
Table 3 shows the results of the
quality of match from the coverage of the sampled sensors to the coverage goal
for varying numbers of sensors
. As in the one-dimensional case, the qualitative behavior of the match for the sampling approach is similar to that of the numerically intensive optimization approach.
As a second example of the two-dimensional situation, consider a scenario with a non-homogeneous environment in which the goal is to achieve the same coverage as shown in
Figure 10. In this scenario, an individual sensor located at
has a range
that depends explicitly on its position, yet the probability
remains constant at
. This is common when the sensing modality has physical properties that are heavily dependent on local environmental conditions. Consider the variation in range to be given by the function
shown in
Figure 13. This particular function was generated by taking the four values of
at the corners of
to be
and performing a two-dimensional linear interpolation to obtain the values throughout
. For this non-homogeneous scenario, the goal of the placement procedure is to find the positions for placing
sensors in the domain
.
For the sampling procedure, the sensor density
from Equation (
10) is no longer a simple scaling of the desired coverage
, but instead is given as shown in
Figure 14. Here there is still a desire to place more overlapping sensors in the center to achieve the higher desired coverage there, but there is now a trend to place more sensors to the bottom left (of both the disc in the center as well as the overall domain) in order to compensate for the range variations as per
Figure 13. The sample-based resulting positions
and coverage
are shown in
Figure 15. As in the homogeneous case, it is clear that the sampling approach again provides a balance between spreading out sensors across
to achieve the lower coverage while allowing some overlap near the center of
to achieve the larger desired coverage there. However, as opposed to the homogeneous case, there are more sensors near the bottom left to account for the lower sensor ranges
found there.
For the optimization procedure for the non-homogeneous two-dimensional example, the same genetic algorithm procedure was used as for the homogeneous example (including the same parameter settings). The resulting optimal sensor placements
and the corresponding coverage
are shown in
Figure 16. This optimized example provides a very close match to the desired coverage
of
Figure 10. Note that this solution has the same qualitative features of the sampling-based approach, in that the balance between overlapping sensors in the middle and those spread out around the remainder of the domain
is similar to that seen in the sampling-based approach. Also, both approaches are affected similarly by the range variation that was shown in
Figure 13. While the optimization is clearly a better result, what is important here is that the sampling-based approach provides much of the qualitative features of the optimal solution, without the need to run an optimization algorithm. Thus, a rapidly computed analytical solution to finding where to position sensors has been obtained that provides performance that is close to the performance of a large-scale computational approach.
Table 4 shows the results of the
quality of match from the coverage of the sampled sensors to the coverage goal
for varying numbers of sensors
. The results are similar to those seen in the homogeneous case.