Above-Threshold Queries of Environmental Conditions Based on Bilinear Interpolation in Wireless Sensor Networks

Wireless sensor networks can be regarded as sensor database systems, which permit users to query sensor data of interest. Among various spatial database queries, we focus the area-wise aggregate queries in the region where the sensor values are above a predefined threshold, which are summarized as above-threshold queries. In this paper, we propose a novel Bilinear Interpolation-Based (BIB) algorithm, which utilizes the bilinear interpolation to estimate the environmental variables inside a grid with the known sensor values at the vertexes, to support the above-threshold queries for regularly-deployed sensor networks and provide the closed-form solution of the above-threshold ratio. We designate experiments with both the artificially-constructed environment data and the real temperature data. Experiment results manifest that the proposed BIB algorithm shows a good performance in estimating the above-threshold ratios to support the above-threshold queries in an accurate and efficient manner.


Introduction
Wireless sensor networks are sharply emerging as one of the most significant technologies to bridge the gap between the physical world and the digital world [1][2][3]. In order to sense spatial natural phenomena in the physical world remotely, wireless sensor networks usually consist of many smart energy-efficient devices, namely sensor nodes that are equipped with inexpensive embedded sensors, processors, memories, radio communication modules, and power supplies to gather environmental conditions, ranging from temperature, precipitation, smoke, soil humidity to atmospheric pressure, to name a few. Wireless sensor networks provide us prolific spatiotemporal environment data and will play important roles in many applications, including environmental monitoring, target tracking, precision agriculture, and so forth [4][5][6][7][8][9].
Since data are constantly being sent from sensors, wireless sensor networks can be regarded as sensor database systems, such as TinyDB [10]. The ability to apply various kinds of data-querying techniques to sensors and process data is highly desirable. Various methods and techniques of data storage and data query have been proposed, which permit users to manage sensor data of interest [11,12].
Spatial database queries, which gather environmental variables, from the perspective of geography, inside a particular geographical region of interest (ROI), have become one of the primary technologies due to the spatial characteristics of the collected environmental conditions [13,14]. As one query algorithms, we first use the Kriging interpolation, which can be regarded as the gold standard to compare the interpolation weights of various algorithms. The results of theory analysis show that the proposed BIB above-threshold algorithm is an approximation of the Kriging algorithm, while the NN algorithm has a large calculation error in geographic data analysis. Then, we perform the artificially-constructed environment data experiment and real temperature data experiment to draw the boundaries and check the metrics of False Positives (FP), False Negatives (FN), and Mean Squared Error (MSE). Experiment results show that the proposed BIB algorithm has a good performance in estimating the above-threshold ratios to support the above-threshold queries in an accurate and efficient manner.
The rest of the paper is organized as follows. Section 2 discusses related works on the applications of the regular deployment sensor networks and enumerates the database management techniques for window queries and aggregate queries. Section 3 presents the system model and formulates the problem. Section 4 introduces the bilinear interpolation-based above-threshold query algorithm. Section 5 describes the experiments and discusses the results. Finally, Section 6 concludes this paper.

Related Work
Quite extensive distributed database management techniques have been proposed to support the spatial queries, mainly the window queries. Berchtold et al. [30] proposed a pyramid technique to index data in the high-dimensional database system, which is efficient for window queries. In Ref. [31], Ooi et al. presented another tunable index scheme called iMinMaxfor high-dimensional window queries, which was less complex and performed better than the pyramid technique. However, there are two major issues impeding the window queries: the lack of explicit window semantics and the low implementation efficiency. To solve this problem, Li et al. [22] established a window query technique framework by defining window semantics to avoid the intra-operator buffers and the reprocessing tuples. However, the previous works have not yet explored the spatial property in sensor network databases. Xu et al. [32] investigated a framework of spatial window query in dynamic sensor networks. To further simplify the query process, Xu et al. [33] proposed an infrastructure-free window query technique in sensor networks. Ku et al. [34] presented an approach to reduce the access latency of the spatial queries, e.g., window queries, by utilizing the neighbor information in wireless broadcast environments. The geographical location information can guide the window queries due to the spatial correlation, and the inspiration of the spatial-aware approach can also be used in wireless sensor networks. The existing database management techniques mainly focus on the geographical coordinate-type window queries, i.e., the traditional window queries, rather than the spatial region determined by environmental values. The environmental value-determined spatial region was also studied [19,23]. Note that the window queries are not only restricted to the spatial coordinators, but also applied to the environmental conditions. The latter window queries to retrieve the sensor values inside a predefined interval were discussed [16][17][18]. It is obvious that such queries degenerate into above-threshold ones if the lower bound of the query window is the threshold and the upper bound is assigned to the positive infinite.
The spatial queries combined with aggregate operation, which are closely related to the above-threshold query techniques in this paper, have also been studied and used in some environmental monitoring applications. In order to figure out the influence factors of climate changes, Wang et al. [35] established a window query algorithm to support various aggregate queries on three-dimensional time-varying datasets, such as the AVERAGE, MIN, and MAX operations.
The deployment is significant for sensor networks. A regular deployment pattern is important in some scenarios such as disaster alert [24][25][26]. In Ref. [36], Bulusu proposed a GPS (Global Positioning System)-less low-cost outdoor localization model for small devices, e.g., sensors, which are placed in a regular mesh structure. Wang et al. [37] proposed a deployment algorithm in an arbitrarily-shaped region and used the regular sensor patterns in a region with boundaries and obstacles. In Ref. [38], the typical regularly-deployed oilfield with an area of 8 × 8 km required approximately a 50-100 m distance between the sensors in seismic analysis. In Ref. [23], Aslan proposed a forest fire detection and monitoring framework and implemented it in regularly-deployed wireless sensor networks. Two popular layout models-square layout, i.e., temperature sensors are placed at the corners of squares, and hexagonal layout, i.e., sensors are placed at the corners of the hexagons-were adopted.
To estimate the environmental conditions at a given location to support spatial queries, the NN and MS algorithms are usually used. The former estimates the environmental conditions only by the nearest sensors, while the latter generates the isobars by the linear interpolation on the edges. Both are applicable to the above-threshold queries. The NN algorithm is simple, and its computation workload and local communication traffic load are small. The MS algorithm is more accurate by taking the environmental conditions on the edges into account. However, neither adequately exploits the sensor values to better the interpolation of the environmental conditions inside the region.
Thus, the idea to fully exploit the sensor values to better the interpolation inside the region inspires us to propose a novel BIB algorithm for the above-threshold queries. With small incrementation of the computation workload and the local communication traffic, the proposed BIB algorithm performs interpolation inside the region and deduces a closed-form above-threshold query solution. It is not only more accurate than the NN and MS algorithms, but also accordant with the gold standard in the geographical sense.

System Model and Problem Formulation
In this section, we present the system model and formulate the problem.

System Model
Consider a wireless sensor network deployed in a predefined square region R, consisting of N sensor nodes, e 1 , e 2 , . . . , e N . Each sensor node senses and then transfers the environment condition data to the database server. Let f (x, y) be the environmental variable at the position with the coordinate (x, y). For a given threshold T, the above-threshold (sub-)region query is used to evaluate the above-threshold region R T , which is defined as: (1) The above-threshold ratio is then defined as: where |X| is the cardinality, i.e., the area of the region X.
Suppose that the region R is split into many regular non-overlapping rectangle grids, as depicted in Figure 1. All the rectangle subregions are listed as S = {S 1 , S 2 , . . . , S m }, where m is the total number of rectangle subregions.
There are four sensor nodes distributed at four vertexes of S i . Sensors can gather environmental variables, which can be used to estimate the above-threshold ratio. Notice that adjacent rectangles may share some sensor nodes. Then, the area-wise above-threshold ratio is given by: where S Ti is defined as

Problem Formulation
Our goal is to estimate the above-threshold ratio C T . From Equation (3), we can obtain C T by summing the above-threshold area limited inside each S i . Thus, the next objective is to calculate |S Ti | according to the coordinates and the sensor values of the four vertexes in S i .
Without loss of generality and for the sake of simplicity, we only consider one of the rectangle subregions and rewrite such a region S i as S. The vertex coordinates of rectangle S are A(x A , y A ), B(x B , y B ), C(x C , y C ), and D(x D , y D ), and the sensor values at four vertexes are f A , f B , f C , and f D , respectively. Let W and H be the length and the height of S, as shown in Figure 2. Figure 2. Illustration of rectangle subregion characteristics. Now, the goal is to calculate the above-threshold area inside S with length W and height H under the condition of known coordinates and the sensor values of four vertexes. The above-threshold region S T in S is defined as:

Bilinear Interpolation-Based Above-Threshold Query Algorithm
In this section, we propose the bilinear interpolation-based algorithm and describe the nearest neighbor algorithm and the marching square algorithm. Furthermore, we introduce the evaluative metrics of mean squared error, false positive, and false negative.

Coordinate Transformation and Boundary of the Above-Threshold Region
In order to simplify the calculation, we perform a translation and an affine transformation to convert the coordinate ( To obtain the original above-threshold area, the result in this subsection should be multiplied by W · H due to the area unit dx · dy = W · H · du · dv.
We utilize the bilinear interpolation to estimate the environmental variables inside S. Bilinear interpolation is a widespread interpolation algorithm for two variables on a two-dimensional grid. It is an extension of linear interpolation for one variable. The key idea of bilinear interpolation is to perform linear interpolation first in one direction and then again in the other direction. Although each step is linear, the bilinear interpolation as a whole is not linear, but rather, quadratic.
The sensor value at arbitrary point (u, v) inside S can be interpolated by the sensor values of four vertexes based on bilinear interpolation. Denotef (u, v) as the estimated environmental variable. Thus, 1) =a 11 uv + a 10 u + a 01 v + a 00 (6) where: For the four-vertex interpolation in , the estimated sensor value of arbitrary coordinate (u, v) in the square is given by: where W (u, v) is a 4 × 1 matrix, which implies the interpolation weights of four vertexes. We can conclude from Equations (6) and (7) that: Denote the estimated above-threshold regionŜ The boundary ofŜ T , which lies between the above-threshold part and the below-threshold one, is given by: The boundaryf (u, v) = T is hyperbolic if a 11 = 0 and is degraded to a straight line if a 11 = 0. Various situations for the boundaries and the above-threshold regions with the BIB algorithm are shown in Figure 3. The solid black circles on the vertexes represent the above-threshold values, and the hollow ones represent the below-threshold values. The threshold T is assigned as zero. The black lines are the boundaries wheref (u, v) = T, and the shaded areas are the above-threshold regions wherê f (u, v) ≥ T. In such a concrete case, the sensor values at points A, B, C, and D are assigned to ±3, ±2, ±4, and ±1, respectively, depending on whether the sensor values are above or below the threshold.

Above-Threshold Area and Above-Threshold Ratio
The above-threshold area |Ŝ T | can be discussed for different cases. If the boundary is hyperbolic (a 11 = 0), the expression off (u, v) can be transformed to: where:û The . Let G E denote the regionûv ≥ E, defined as: If a 11 > 0,Ŝ T is the intersection of the region G E and the regionŜ, or precisely, Without loss of generality, we only discuss the case of a 11 > 0. The illustrations of the boundaries and the above-threshold regions in theûov plane if a 11 > 0 are shown in Figure 4 We denote an auxiliary function k(E,û 0 ,v 0 ) as is the characteristic function of the set X, which satisfies: Furthermore, let us define: Denote the symbol "+" as the positive sign of E,û 0 , orv 0 and the symbol "−" as the negative one. For example, the notation "+ + +" of the symbol k +++ (E,û 0 ,v 0 ) indicates that all the independent variables E,û 0 , andv 0 are positive.
In the coordinate system (û,v), the original coordinates of four vertexes are converted into: For similarity, let us define H E (P) k(E, u P , v P ). When considering various practical situations, no matter where the square S lies, the above-threshold area |Ŝ T | is always provided by: When a 11 < 0, the above-threshold area |Ŝ T | is the total area |S|, subtracting the area |S T | for a 11 > 0. We have: If a 11 = 0, the boundary is degraded into a straight line and can be transformed to: Then, the above-threshold area |Ŝ T | can be calculated according to the following four steps: Step 3. Similarly, judge if the point B, the intersection between BC andf (u, v) = T, the point C, the intersection between CD andf (u, v) = T, the point D, and the intersection between DA and f (u, v) = T should be added to the ordered set U according to the principles in Step 1 and Step 2.
Step 4. If ordered set U = ∅, |Ŝ T | = 0. Otherwise, |Ŝ T | is equal to the area of the polygon of which the vertexes are listed in the ordered set U.
So far, the above-threshold area |Ŝ T | inside the subregion S is obtained by the bilinear interpolation for all the cases.
We propose a bilinear interpolation-based algorithm to estimate the above-threshold ratio in a predefined region, as shown in Algorithm 1.   if j in [A, B, C, D] then 20: iff (u j , v j ) ≥ T then 21: Add j to U.

22:
end if 23: end if 24: if j in [AB, BC, CD, DA] then 25: if j is intersected with the linef (u, v) = T then 26: Add the intersection to the ordered set U.

Nearest Neighbor Algorithm and Marching Square Algorithm
Both the NN and MS algorithm utilize the sensor values of four vertexes in the rectangle subregion.
In the NN algorithm, the sensor value of arbitrary coordinate (u, v) in subregion S is determined by that of the nearest vertex, and the interpolation weight W NN (u, v) in Equation (8) can be expressed as: The illustrations of the boundaries and the above-threshold regions with the NN algorithm are shown in Figure 7. The MS algorithm is a computer graphics algorithm, and thus, its interpolation weight does not exist. The illustrations of the boundaries and the above-threshold regions with the MS algorithm are depicted in Figure 8. Note that there are two ambiguous situations for above-threshold calculation and boundary description. The MS algorithm is possibly considered as an approximation of the BIB algorithm, which transforms the hyperbolic boundaries by the BIB algorithm into the straight lines.
Note that in order to estimate the above-threshold coverage ratio in region R, we sum the above-threshold area limited inside each rectangle grid S i for the three algorithms. Thus, the algorithm complexities of the three mentioned algorithms are both O(M), where M is the number of grids.

OR OR
? ?

Evaluative Metrics
The MSE in mathematical statistics refers to the mean value of the squares of the differences between the estimated values and true values, which is given by: where observed t represents the exact above-threshold ratio in the tth experiment, predicted t is the estimated above-threshold ratio obtained by different algorithms in the tth experiment, and N e is the number of experiments. A lesser MSE is better. In statistical hypothesis testing, the FP and FN are known as type I and type II errors, where FP corresponds to the situation when the exact value is false, while the estimated value is true and FN is opposite. True Positive (TP) represents that the exact value and the estimated value are both true, and True Negative (TN) represents both false. In the above-threshold queries, the conditions of FP, FN, TP, and TN are listed in Table 2. The probability of FP and FN in a single experiment can be expressed as:

Experiments and Discussion
In this section, we first testify to the correctness of the proposed algorithm by theory analysis. Then, we describe the experiments with both the artificially-constructed environment data and the real temperature data. All programming and benchmarking were done in Python (Version 3.6) with the Numpy [39] and pyKrige [40] libraries running on a computer with a 2600-MHz Intel processor and 8-gigabyte memories. Furthermore, we illustrate the possible implementation frameworks in wireless sensor network and analyze the algorithm complexity. Finally, we provide the evaluation and the comparison of the NN, MS, and BIB algorithms.

Theory Analysis
In geography, the Kriging algorithm is one of the most wide-spread interpolation algorithms and can be regarded as the gold standard [41].
We use the local four-point Kriging (Kr4) interpolation to validate the performance of the proposed algorithm. Assume that the Kr4 interpolation is controlled by the four vertexes of the square. The semi-variable function is linear, and the nugget is zero. The interpolation weight W Kr4 (u, v) in Equation (8) can be calculated as: where: We denote the ith elements in , respectively, and these values are shown in Figure 9 for i = 1. We can conclude that the weights W 1 Kr4 (u, v) and W 1 BIB (u, v) are almost identical, while the weight W 1 NN (u, v) is quite different from the other two.    The differences of weights, i.e., |W 1 Figure 10. The largest weight error max Kr4 (u, v)| is less than 0.03 in the BIB algorithm, while max Kr4 (u, v)| is more than 0.5 in the NN algorithm. The same conclusion holds for the other weights W i Kr4 (u, v), W i BIB (u, v), and W i NN (u, v), i = 2, 3, 4. Therefore, the interpolation result of the BIB algorithm almost approximates that of Kr4 algorithm. Note that the above-threshold area is difficult to be obtained by the Kr4 algorithm due to the complex calculation of the integral inside a complex boundary. Moreover, the NN algorithm has a large weight error in geographic data analysis.

Experiment on Artificially-Constructed Environment Data
In order to evaluate different above-threshold query algorithms, we test them with artificial data generated by the Kriging algorithm. The interpolated values by the Kriging algorithm are modeled as a Gaussian process governed by full prior covariances.
Assume that the Kriging interpolation is controlled by N c points P i , 1 ≤ i ≤ N c , randomly distributed in R, and the value f at P i is an independent random variable uniformly distributed in the interval [a, b]. Typical artificial environmental variables are generated by N c control points inside a 300 × 300 region. We denote the central 100 × 100 region as the ROI in which the sensor values are depicted in Figure 11 in a single experiment. The greater the number of the control points is, the richer the details of the variances of artificial environmental variables are. To compare different algorithms with the same number of sensors, we assume that there are in total N = (n + 1) × (n + 1) sensors and M = n × n grids.
We performed simulation experiments to contrast the three algorithms, NN, MS, and BIB, to estimate the above-threshold ratios. The simulation configuration is given in Table 3. Table 3. Simulation configuration.

Parameter Value
Side length of R 100 The real boundary is shown in Figure 12a, and the boundaries generated by the three algorithms in ROI when n = 3 are drawn in Figure 12b-d. Obviously, the boundaries by the BIB and MS algorithms are more accurate than those by the NN algorithm. Furthermore, the boundaries by the BIB algorithm are smoother than those by the other algorithms. In order to analyze the three aforementioned algorithms qualitatively, we used the metrics FP, FN, and MSE to evaluate their performance for N e experiments. The average FP, FN, and MSE versus the number of grids M for the three mentioned algorithms are shown in Figure 13.    The threshold T = f max + f min 2 , where f max and f min are the maximum and the minimum of the variables. First, all the FP, FN, and MSE metrics obtained by the three algorithms tend to a definite value, zero, along with the increasing of the number of sensors. This implies all algorithms are feasible and more sensors help to reduce the query error. Second, the NN algorithm performs worse than the others. This is mainly caused by the rough estimating of the variables inside the ROI. Since the proposed BIB algorithm is an approximation of the Kr4 algorithm and the MS algorithm is an approximation of the BIB algorithm, the MS and BIB algorithms both performed well.

Experiment of Real Temperature Data
WorldClim [42] is a set of gridded global climate data with a spatial resolution from 1-340 km 2 over 30 years, which was created by Steve Fick and Robert Hijmans. WorldClim provides various environmental variables, including minimum temperature ( • C), maximum temperature ( • C), average temperature ( • C), precipitation (mm), solar radiation (kJm −2 day −1 ), wind speed (ms −1 ) and water vapor pressure (kPa). These data can be used for mapping and spatial modeling.
The global average monthly temperature data, that is 12 sets of temperature data, obtained by WorldClim from 1970-2000, were used in the following experiments. We first randomly selected an ROI that did not contain a water surface, ranging from 42-78 • E and from 49-65 • N, to observe the boundaries. The hypsometric map of the ROI is depicted in Figure 14. Rectangle divisions were done to the ROI, and M = n × n denotes the grid numbers. The real boundary and the boundaries generated by the three algorithms in ROI when n = 20 are depicted in Figure 15.
Obviously, the boundaries by the BIB algorithm are more accurate than those by the NN and MS algorithms. Furthermore, the boundaries by the BIB algorithm are smoother than those by the other algorithms. Then, we randomly selected 20 ROIs that did not contain a water surface with a length of 10 • both in the latitude and longitude direction in Eurasia to demonstrate the universality of the proposed algorithm. We assumed that in the case that the distance between the sampling points was 30 s, i.e., the most precise spatial resolution, the average above-threshold ratio obtained by the NN algorithm was the exact above-threshold ratio. The average FP, FN, and MSE versus the number of grids M for the three mentioned algorithms in 20 ROIs are shown in Figure 16. The threshold T was assigned to the median of the temperature data. Although the effects of the results cannot be favorably compared with those in the artificial experiment due to the discontinuity points such as the mountain segregation, the shapes and the trends of the curves in the real temperature data experiment are the same as those in the artificially-constructed environment data experiment. Above all, the proposed BIB algorithm shows a good performance in estimating the above-threshold ratio and solving the above-threshold query equation within the given ROI.

Possible Implementation Framework and Computation/Communication Load Analysis
A naive implementation framework to support the aforementioned above-threshold queries with the BIB, NN, or MS algorithms is to collect all the sensor data and then process the data in the center server. All the positions of the nodes were assumed to be obtained with node location techniques. A possible data collection route tree to transfer the sensor data from nodes is shown in Figure 17a. In such a naive implementation, the nodes only send out their data, and the center server needs to execute the query, which is thus named the centralized query framework. The traffic loads in the sensor network of different algorithms were equal. However, the center server needed more computation for the BIB and MS algorithms to perform the interpolation than that for the NN algorithm, which in fact simply counted the number of the data above the threshold.
However, in most practical scenarios, there would be only a small proportion of the sensor values above the threshold. We herein propose a distributed implementation framework to support the aforementioned above-threshold queries.
In this distributed implementation framework to support the NN algorithm, the sensor will transfer its sensor value if and only if it is triggered by the above-threshold event, i.e., its sensor value is beyond the threshold, as shown in Figure 17b. The nodes in the route path from the triggered node to the sink node help to forward the data hop by hop. Therefore, the center server only counted the data from the triggered nodes, and the values of other nodes were assumed to be below the threshold. Now, the network traffic was significantly reduced if the above-threshold events occurred sparsely. It is more complicated to support the MS and BIB algorithms since both need local co-operation. We first grouped the sensor nodes into overlapped clusters. All four nodes in the vertexes of a square grid make a cluster. One node possibly belongs to a different cluster since it may be shared by neighbor grids. There was a cluster head for each cluster, as shown in Figure 17c. When a node, say, s 1 , is triggered by the above-threshold event, it will send its sensor value to notify the cluster heads of all the clusters to which it belongs. Each cluster head notified, say s 4 , will command its cluster members to submit their sensor values to it. Once the header collects all the sensor values in the cluster, it will perform the computation and transfer the fused result to the sink hop by hop. Therefore, the MS and BIB algorithms needed more local traffic and more computation to exchange their sensor values for the interpolation.

Comparison and Discussion
We evaluate and compare the three mentioned algorithms, NN, MS, and BIB, according to various evaluative metrics in Table 4. We can see from Figures 12 and 15 that the curve boundaries by the BIB algorithm are more accurate and smoother than the straight line boundaries by the other two algorithms. We can see from Figures 13 and 16 that all algorithms are feasible and that NN algorithm performed worse than the others. The largest weight error is less than 0.03 in the BIB algorithm, while it is more than 0.5 in the NN algorithm. Therefore, the interpolation result of the BIB algorithm almost approximates that of the Kr4 algorithm. Since the MS algorithm is possibly considered as an approximation of the BIB algorithm, the same conclusion holds for it. Both the BIB and the MS algorithms have small weight errors in geographic data analysis, while the NN algorithm have the larger one. The interpolation weights W BIB (u, v) and W NN (u, v) are given by Equations (9) and (35), respectively. The MS algorithm is a computer graphics algorithm, and thus, its interpolation weight does not exist. We can see from Figure 8 that there were two ambiguous situations for above-threshold area calculation and boundary description in the MS algorithm. For the distributed implementation framework, once the header collects all the sensor values in the cluster, it will perform the computation and transfer the fused result to the sink hop by hop in the MS and BIB algorithms. The center server only counts the data from the triggered nodes and the values of other nodes were assumed to be below the threshold in NN algorithm. Therefore, the MS and BIB algorithms need more local traffic and more computation to exchange their sensor values for the interpolation. The traffic and computation are modest in the MS and BIB algorithms due to the small above-threshold clusters. Taking all the factors into consideration, the proposed BIB algorithm has an excellent comprehensive performance compared with the other algorithms.

Conclusions
Wireless sensor networks can be regarded as sensor database systems, which permit users to query sensor data of interest. In this paper, we proposed a novel Bilinear Interpolation-Based (BIB) algorithm to support the above-threshold queries. It is mainly designed for regularly-deployed sensor networks. We utilize the bilinear interpolation to estimate the environmental variables inside a square with the known sensor values from sensor nodes located at its four vertexes. Thus, the area of the above-threshold subregion can be obtained by mathematical formulation, and finally, we proposed the closed-form solution of the above-threshold ratio.
In order to evaluate the performance of different above-threshold query algorithms, we first used the interpolation weights of Kr4 interpolation, which is commonly used in geometry research, as a golden standard to compare various algorithms. The weights of the proposed BIB algorithm are approximate to those of the Kr4 interpolation, which means BIB is valid from the perspective of geography.
We also described experiments both with the artificially-constructed environment data and with real temperature data. Experiment results manifest that the proposed BIB algorithm shows a good performance in estimating the above-threshold ratio.
In summary, with a closed-form formulation, the proposed BIB algorithm supports the above-threshold queries in an efficient and accurate manner.