UV-diagram: a voronoi diagram for uncertain spatial databases

The Voronoi diagram is an important technique for answering nearest-neighbor queries for spatial databases. We study how the Voronoi diagram can be used for uncertain spatial data, which are inherent in scientific and business applications. Specifically, we propose the Uncertain-Voronoi diagram (or UV-diagram), which divides the data space into disjoint “UV-partitions”. Each UV-partition \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P$$\end{document} is associated with a set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document} of objects, such that any point \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q$$\end{document} located in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P$$\end{document} has the set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document} as its nearest neighbor with nonzero probabilities. The UV-diagram enables queries that return objects with nonzero chances of being the nearest neighbor (NN) of a given point \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q$$\end{document}. It supports “continuous nearest-neighbor search”, which refreshes the set of NN objects of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q$$\end{document}, as the position of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q$$\end{document} changes. It also allows the analysis of nearest-neighbor information, for example, to find out the number of objects that are the nearest neighbors of any point in a given area. A UV-diagram requires exponential construction and storage costs. To tackle these problems, we devise an alternative representation of a UV-diagram, by using a set of UV-cells. A UV-cell of an object \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$o$$\end{document} is the extent \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e$$\end{document} for which \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$o$$\end{document} can be the nearest neighbor of any point \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q \in e$$\end{document}. We study how to speed up the derivation of UV-cells by considering its nearby objects. We also use the UV-cells to design the UV-index, which supports different queries, and can be constructed in polynomial time. We have performed extensive experiments on both real and synthetic data to validate the efficiency of our approaches.


Introduction
The Voronoi diagram, primarily designed for evaluating nearest-neighbor queries over two-dimensional spatial point [33], has raised plenty of research interest.This technique has been extended to handle different related problems, including database services in wireless broadcast environments [44,45]; high-dimensional query evaluation [7]; continuous location-based services [4,32,43]; and virus spread analysis among mobile devices [41].Conceptually, the Voronoi diagram partitions the data space into disjoint "Voronoi cells", so that all points in the same Voronoi cell have the same nearest neighbor.The task of finding the nearest neighbor of a query point is then reduced to a point query.Figure 1a illustrates a Voronoi diagram of seven points.Since the query point q is located in the Voronoi cell of O 2 , O 2 is the nearest neighbor of q.
Is it possible to use the Voronoi diagram to perform nearest-neighbor search on objects whose values are impre- cise?Data values can be uncertain for a variety of reasons.Consider a satellite image, which depicts geographical objects like airports, vehicles, and people.Using machine learning and human effort (e.g., community-based systems like Wikimapia), the location of each object on the image can be obtained.Due to the noisy transmission of satellite data, the quality of these images can be affected, and we may not be able to obtain very accurate locations.Moreover, if this location information is released to the public (e.g, for research purposes), it may need to be preprocessed for privacy purposes.In fact, recent proposals like [1,2] have suggested to represent a user's position as a larger region, in order to lower the likelihood that a user is identified at a particular site.Uncertainty is also inherent in biological data management.For example, microscopy images have been actively used to analyze the thickness of neuron layers in the retina, as well as the extent of the area of a cell.Due to factors like image resolution and measurement accuracy, it is hard to obtain exact values of the objects of interest [28,29].For this kind of data, techniques for evaluating range queries, nearest-neighbor queries, and joins have been developed.These queries return answers with probabilistic guarantees, which reflect the confidence of answers due to data uncertainty.For these applications, tools that resemble the Voronoi diagram can be potentially useful.Specifically, we would like to examine space-partitioning techniques for performing a Probabilistic Nearest-Neighbor Query (PNN).Given a query point q, a PNN returns the IDs of objects with nonzero probabilities for being the closest to q, as well as their probabilities.In the sequel, we denote the objects returned by the PNN as answer objects, and their probability values as qualification probabilities.
An uncertainty model that has been commonly used is to assume that an object O i has an "uncertainty region" and a probability distribution function (pdf).This means that the precise position of O i can only be located inside the (closed) region, with a pdf that describes the distribution of the object's position within the region.The uncertainty region can have any shape, and the pdf is arbitrary (e.g., it can be a uniform distribution, Gaussian, or a histogram).Here, we assume that O i has a two-dimensional circular uncer-tainty region.We will also explain how our solution can be extended to handle non-circular-shaped regions.We examine how the Voronoi diagram should be defined to support PNN execution.Specifically, we propose the Uncertain-Voronoi diagram (or UV-diagram), where the nearest-neighbor information of every point in the data space is recorded, based on the uncertain objects involved.The UV-diagram provides a basis for studying solutions that used the Voronoi diagram for point data.It could be interesting, for instance, to extend the solution of [44] to support uncertain data in broadcasting services.Figure 1b illustrates an example of the UV-diagram of seven uncertain objects, where the space is divided into disjoint regions called UV-partitions.Each UV-partition P is associated with a set S of one or more objects.For any point q located inside P, S is the set of answer objects of q (i.e., each object in S has a nonzero probability for being the nearest neighbor of q).The highlighted regions contain points that have two or more nearest-neighbor objects.As an example, since q 1 is inside the dashed region, O 4 has a nonzero probability for being the nearest neighbor of q 1 ; on the other hand, q 2 is located inside the dotted region, and O 6 and O 7 are the answer objects for the PNN with q 2 as the query point.Observe that the Voronoi diagram, which indexes on spatial points, is a special case of the UV-diagram, since a point can be viewed as an uncertainty region with a zero radius.Figure 1 compares the two diagrams.
The Voronoi diagram can also be used in other applications.For example, a continuous nearest-neighbor query, which constantly returns the nearest neighbor (e.g., gas station) of a moving point q (e.g., a vehicle), is a typical operation in location-based services [32,43].The Voronoi diagram supports this query; particularly, the Voronoi cell that contains the current location of q can be easily retrieved.We will illustrate how to use the UV-diagram to track the possible nearest neighbors of a moving point.Another use of the Voronoi diagram is to perform data analysis or observe interesting patterns of nearest-neighbor information.In [41], the Voronoi diagram is used to investigate the spreading pattern of bluetooth viruses among mobile users.We can also use UV-diagram to provide valuable information about these "nearest-neighbor patterns".In Fig. 1b, if the dashed region is large, it indicates that O 4 has high chance to be placed in different clusters (assuming that a nearest-neighbor clustering algorithm is used).Another interesting query is as follows: given a region R, display all UV-partitions that intersect with R, as well as the density of objects that can be the nearest neighbor in each UV-partition.Hence, a UV-diagram allows a user to visualize patterns about the nearest-neighbor information.
Challenges of constructing UV-diagram.Although the UV-diagram is useful, developing a UV-diagram is not simple.Notice that the UV-partitions are produced based on uncertainty regions, which may not be points.Unfortunately, Fig. 2 A UV-diagram for 3 uncertain objects efficient computational geometry methods for generating the Voronoi diagram (e.g., line-sweeping [19]) cannot be readily used for creating a UV-diagram, since these methods are primarily designed for spatial points, rather than uncertainty regions.Figure 2 depicts the space partition based on three uncertainty regions represented as circles.Each UV-partition (named R i , where i = 1, . . ., 7) is irregular in shape and contains different answer objects, listed on the side of the figure.In general, given a set of uncertain regions, an exponential number of UV-partitions can be created.For example, Fig. 2 shows that for three objects, there are seven UV-partitions, each of which contains one of 2 3 − 1 = 7 combinations of the three objects.Moreover, the number of edges of each UV-partition can also be exponentially large!This makes it computationally infeasible to generate and store these partitions.It is also difficult to find out which of these irregular UV-partitions contain a given query point.Indeed, our experimental results show that a brute-force approach of computing and indexing UV-partitions over 40,000 objects requires about 60 h.Therefore, a scalable method for constructing a UV-diagram is highly desirable.
Our solutions.Instead of computing UV-partitions, we have developed an alternative interpretation of the UVdiagram.For every object O i , we consider the extent a i such that O i can be the nearest neighbor of any point selected from a i .We call this extent the UV-cell of O i .We examine some basic properties of a UV-cell (e.g., its size and number of edges).We show how to represent a UV-cell as a set of objects, and develop novel methods to find this object set efficiently.For example, our batch-construction algorithm allows the UV-cells of objects that are physically close to each other to be swiftly obtained.We propose a polynomial time method for constructing an index for the UV-partitions, called the UV-index.We adopt an adaptive-grid indexing scheme, which has the advantage of adapting to different distributions of uncertain objects' positions.Our experimental results show that on both synthetic and real dataset, this index can be constructed in a much shorter time.We also explain how to use the UV-index to support different applications (e.g., PNN and nearest-neighbor pattern queries).
To summarize, our contributions are the following: -Study the UV-diagram and its basic properties; -Propose efficient algorithms for obtaining a UV-cell; -Design the UV-index; -Use the UV-index to support different queries; and -Conduct experiments on real and synthetic datasets.
The rest of the paper is as follows.Section 2 summarizes the related work.In Sect.3, we present basic concepts of the UV-diagram.In Sect.4, we study the UV-cell and its essential properties.We then explain how to represent UVcell efficiently in Sect. 5.An adaptive index based on the UV-diagram is presented in Sect.6.We present experimental results in Sect.7. Section 8 concludes the paper.

Related work
Data uncertainty management.Recently, researchers have proposed to consider uncertainty as a "first-class citizen" in a DBMS [13,14,18,39].Two models can be used to represent uncertain data: tuple-and attribute-uncertainty.For tuple uncertainty, each database tuple has a probability of being correct [18].Here, we assume attribute uncertainty, which represents an attribute as a range of possible values and a probability distribution function (pdf) bounded in the range [39].Common queries for attribute uncertainty include range queries [16], k-nearest-neighbors [28], skylines [25,36], and top-k queries [20].
A few works have been proposed to evaluate PNN queries over attribute uncertainty.In [15], numerical integration techniques have been presented.Probabilistic verifiers, described in [13], can generate answer objects' proba bility bounds without performing expensive integration operations.Another way to compute answer probabilities is based on sampling [24].In this paper, we focus on the efficient retrieval of answer objects.
To our understanding, the only indexing method available for nearest-neighbor search over uncertain data is to use an index like the R-tree and the grid.The R-tree is a disk-based structure that uses the minimum bounding rectangles (MBRs in short) to cluster the uncertainty regions of the objects, and organizes MBRs in a hierarchical manner [6].To evaluate PNN using the R-tree, a branch-and-prune strategy has been proposed in [15], where MBRs that may contain answer objects are traversed.However, this involves a lot of I/O cost in reading index nodes and leaf pages [13,15].Similar issues also occur with grids [31].On the other hand, retrieving answer objects from the UV-diagram is essentially a point query search: given a point q, find the objects associated with the UV-partition that contains q.Hence, a UV-diagram can support more efficient PNN search.
While it is not clear how an R-tree or grid over uncertain objects can provide pattern analysis of nearest-neighbor information (e.g., displaying the extent of a UV-partition), we will show how to use the UV-diagram to provide this information.
Other types of nearest-neighbor queries, like the "group nearest-neighbors" [26], "reverse-nearest-neighbors" [10,27], "uncertain queries" [8], and "continuous nearestneighbor queries" [12] have also been proposed.In these works, the R-tree was used to support object retrieval.It is interesting to study how the UV-diagram can be used to support the execution of these queries.In this paper, we study how to use the UV-diagram to support the execution of continuous nearest-neighbor queries.
The Voronoi diagram is an important technique for answering nearest-neighbor queries over spatial points [33].It has been extended to support other applications (e.g., [7,32,[43][44][45]).It also facilitates the analysis of spreading patterns of mobile viruses [41].In [9], the k-th order Voronoi diagram is used to evaluate a k-NN query.In [38], an index called VoR-Tree is designed to merge Voronoi diagrams into R-tree in order to answer various nearest-neighbor queries.The Voronoi diagram has also been defined for boundaries of circular objects in [23].However, these objects are not uncertain, and the method of [23] cannot be used to answer PNN queries.
Few works have studied the application of the Voronoi diagram on uncertain data.[8] Consider the "uncertain" nearest neighbor query (UNN) over spatial points.Different from PNN, the query is an uncertain region, not a query point.To evaluate a UNN, the authors propose to use a Voronoi diagram over 2D points.The portions of the Voronoi cells that overlap with the query's uncertainty region are then used to compute answer probabilities.[22] Consider the clustering of uncertain attribute data, where a Voronoi diagram is constructed for centroid points.Notice that [8,22] do not construct a Voronoi diagram for uncertain data.On the other hand, the UV-diagram is a Voronoi diagram tailored for attribute uncertainty.
In [21,37], the Voronoi diagram was modified to identify an imprecise object that is surely the nearest object of a query point q.However, the UV-diagram returns object(s) that may have chance to be the nearest neighbor of q, and can be used to answer probabilistic nearest-neighbor queries.We also study a database index for the UV-diagram, which has not been examined in these two works.
This paper is an extension of [17].Here, we improve the performance of UV-index construction by proposing batch pruning, which reduces the workload of generating UV-cells for a set of nearby objects.We provide a more detailed study of the basic properties of a UV-cell (e.g., its size and number of edges).We also examine how the UV-index can be used to answer PNN queries for a moving query point.We conduct new experiments to validate the effectiveness of these approaches.

The UV-diagram
We now present the basic notions of the UV-diagram.
We introduce the UV-cell, an alternative presentation of the UV-diagram, in Sect.3.1.We then study some applications of the UV-diagram, in Sect.3.2.

The UV-cell
As discussed before, the UV-diagram can be expensive to construct.We hence propose an alternative representation of the UV-diagram, by using UV-cells.We will later explain how the UV-cells facilitate efficient construction of the UVdiagram.Now, let O 1 , O 2 , . . ., O n be the IDs of a set O of uncertain objects, and D be a two-dimensional space that contains these objects.For simplicity, we assume that D is a square.The UV-cell is then defined as follows: has a nonzero probability to be the nearest neighbor (NN) of a point q, where q ∈ U i .
Figure 2 illustrates the UV-cells for O 1 , O 2 , and O 3 .The boundary of each UV-cell is labeled with the ID of the object.For example, the UV-cell of O 2 is a region enclosed by solidline segments.
The UV-cell can be used to recover the UV-partitions (i.e., disjoint regions of a UV-diagram).In fact, a UV-partition that contains q is the intersection of all UV-cells that contain q.This is because the objects associated with these UV-cells have nonzero qualification probabilities for q.For instance, in Fig. 2, the UV-cells of both O 1 and O 3 intersect at partitions R 5 and R 7 .This means that when q is located at any of these partitions, both O 1 and O 3 are the answer objects.Since R 7 is intersected by O 2 's UV-cell, O 2 is also associated with R 7 .Therefore, a UV-diagram is the union of all objects' UVcells.Besides, the UV-cells of all objects can be used to output which object(s) is/are the nearest neighbor of q with nonzero probabilities.
Table 1 shows the symbols used in this paper.Notice that if there is at least one uncertain object in domain D, any point in D must be covered by at least one UV-cell.In particular, if O i is the only object in domain D, then its UV-cell is exactly D.

Applications of the UV-diagram
The UV-diagram supports a number of applications.Let us now explain how to use the UV-diagram to handle the following queries: A circle centered at c with radius r dist (q, c i ) Euclidean distance between q and c i 1.The Probabilistic Nearest-neighbor (PNN) Query.This query has been mentioned in Sect. 1.To evaluate a PNN for a given point q, we can find out the UV-partition that contains q.The set A of objects associated with this partition are those that can be the nearest neighbor of q.Notice that the UV-partitions can be obtained by finding the union of all the UV-cells.For each object O i ∈ A, the probability that O i is the closest to q can be efficiently evaluated by using solutions in [13,15,24].
2. The Continuous PNN Search (CPNN), an extension of the PNN, is a query that resides in the processing server for an extensive period of time.Different from PNN, the position of a query point q changes with time [12].The objective of the CPNN is to refresh the PNN answer, when the value of q changes.This query can be used in transportation services, where q can be a moving vehicle or person, and the data can be the geographical objects retrieved from satellite images.Assuming that q reports its position to the server periodically, the UV-diagram can conveniently support CPNN.Suppose that the server receives a new position of q, say, q 1 .A simple solution is to issue a new PNN for q 1 .However, if q 1 is located in the same UV-partition that contains the old position of q, then it suffices to use the objects associated with that UVpartition to compute the query answer for q 1 .The cost of retrieving the UV-partition that contains q 1 is thus saved.
3. The UV-partition Query.The UV-diagram can also be used to retrieve the distribution and pattern information about nearest neighbors, which can be useful for analysis purposes (e.g., [41]).One such "pattern-analysis" operation is the UVpartition query.Given a region R, this query retrieves all UV-partitions inside R and the "density" of each partition R j (which is equal to the number of objects associated with R j , divided by the area of R j ).This allows a user to examine the density distribution of the nearest neighbors in his/her interested area R. 4. The UV-cell Query.This is another pattern analysis operation.Given an object O i , it returns the extent and the area of O i 's UV-cell.The query user can then obtain the area of the region where O i may be the nearest neighbor.This area can reflect the "influence" of O i (in terms of the nearest neighbor information).The shape of the UV-cell can also be displayed on the user's computer screen for further analysis.
Since the UV-diagram is expensive to construct, in Sect.6, we revisit how the above queries can be implemented by the UV-index, which is an approximate version of the UVdiagram.We next address the UV-cell in detail.

More about UV-cells
We now investigate the UV-cell, which is important for constructing the UV-index, in more details.We first present a simple method for constructing a UV-cell in Sect.4.1.We then examine the shape of a UV-cell in Sect.4.2.The number of edges of a UV-cell, and its size, is studied in Sects.4.3 and 4.4, respectively.

Constructing a UV-cell
Let us first address the relationship between a query point and UV-cells.Let p be a point in D, and let dist min (O i , p) and dist max (O i , p) be the minimum and the maximum distances of object O i from p, respectively.Figure 3 illustrates two uncertain objects, O i and O j .For any point p on the solid line shown, we require the following property to hold: We call this solid line the "UV-edge of O i with respect to O j ", denoted by E i ( j).A special property of this edge is that any point p at the region on the side of E i ( j) closer to O j has its maximum distance from O j , that is, dist max (O j , p), shorter than its minimum distance from O i , that is, dist min (O i , p).
On the other hand, if p is on the opposite side of E i ( j), then The UV-edge allows us to decide whether an object is an answer object (i.e., an object with nonzero qualification probabilities).In Fig. 3, q 0 is on the right of E i ( j), which is also closer to O j than O i .Thus, dist max (O j , q 0 ) < Fig. 3 The UV-edge dist min (O i , q 0 ).In other words, O j is always closer to q 0 than O i , and O i has no chance to be the nearest neighbor of q 0 .As another example, q 1 is on the left of E i ( j).Since dist min (O i , q 1 ) ≤ dist max (O j , q 1 ), O i has a nonzero qualification probability.Hence, given E i ( j), if the query point is on the right of E i ( j), O i can be pruned.
We can now present a simple method for constructing a UV-cell.Let us define the outside region: Definition 2 The outside region of UV-edge E i ( j), denoted by X i ( j), is the region on one side of E i ( j) such that for any point q ∈ X i ( j), O j is always closer to q than O i .We call the complement of X i ( j) the inside region, denoted by X i ( j).
Given a set of n objects O, the UV-cell U i of object O i is essentially the intersection of all other n − 1 inside regions: Definition 3 A possible region of object O i , denoted by P i , is the intersection of a set of inside regions: According to the definition, the possible region should be an area that completely covers the UV-cell of O i .An example of an object's possible region is the domain D, since D must cover any UV-cell.Here, R is the empty set.Notice that a UV-cell is also a possible region.
In Fig. 3, the outside region of the UV-edge E i ( j) is the area on the right of the solid line.Notice that since q 0 is in the outside region of E i ( j), O j is closer to q 0 than O i , and thus O i cannot be q 0 's nearest neighbor.
Given an object O i , if we know all the outside regions X i ( j) (where j = 1, . . ., n ∧ j = i), then O i 's UV-cell can be constructed by excluding all these regions from D. Algorithm 1 illustrates the basic method for constructing UV-cell for n objects.The possible region of each object O i is first initialized as the whole space (Step 2).Then, for each O j , we compute the UV-edge of O i and its corresponding outside region (Steps 4 and 5).The possible region, which contains all the points that may have O i as one of their nearest neighbors, is then "reduced" by the outside region that overlaps with it (Step 6).The UV-cell of O i is then assigned to be the final possible region (Step 8).
We now discuss Step 6 in more detail.This step uses E i ( j) to "refine" the edges of P i , that is, find the intersections of E i ( j) with P i .Specifically, for each UV-edge e of P i , we compute the intersection of E i ( j) and e.Since E i ( j) is a hyperbola, we can use the techniques in [3] to do this.The resulting intersections partition E i ( j) into some segments.For each segment s, there are two scenarios: -Case 1: s is inside P i : We refine P i by using s as one of the new edges of P i .Some existing edges of P i are removed if necessary.-Case 2: s is outside P i (except the end points of s): P i cannot be changed by s, and we do not have to do anything to handle this case.
After we have visited all the segments of E i ( j), we have found all the intersections of E i ( j) and P i .Moreover, P i is refined, and Step 6 is completed.
As an example of Step 6, consider Fig. 4a, which illustrates P i , and Fig. 4b, which shows the result of intersecting E i ( j) and P i .The segment between v 4 and v 5 is inside P i .Since Case 1 is satisfied, the existing edges between v 4 and v 5 (i.e., (v 4 , v 1 ) and (v 1 , v 5 )) are replaced by segment (v 4 ,v 5 ).On the other hand, (v 5 , v 6 ) is outside P i , and so Case 2 is satisfied.There is no need to change any edges of P i between v 5 and v 6 .The process is repeated until all the segments are visited.As shown in Fig. 4b, vertices v 1 , v 2 , and v 3 are removed from P i , while v 4 , . . ., v 9 are added to it.
Note that the order of selecting the object for refining O i 's possible region (Steps 4-6) does not affect the correctness of the algorithm; the UV-cell is produced by "shrinking" the possible regions by using the outside regions of other objects.Also, not all objects are useful in shaping a UV-cell.In Sect.5, we will explain how to prune away these objects.

The shape of a UV-cell
We now study a mathematical representation of the UV-cell.We also derive the number of UV-edges of a given UV-cell.Here, we assume that the uncertainty region of O i is a circle, with center c i and radius r i , with r i > 0. Later, we discuss the UV-cell of a "point uncertainty" (i.e., r i = 0), and also uncertainty regions that are not circle in shape.
For any point q ∈ D, we can observe from Fig. 3 that: dist max (O j , q) = dist (q, c j ) + r j (5) where (c i , r i ) denotes a circle with center c i with radius r i .Since r j > 0, dist max (O j , q) must also be positive.By substituting Eqs. 4 and 5 into Eq. 1, we have: Let the coordinates of c i and c j be (x i , y i ) and (x j , y j ).Let . Then, Eq. 6 becomes: where Essentially, Eq. 7 is a hyperbolic equation, with c i and c j as the foci, rotated by θ in an anti-clockwise sense [3]. Figure 3 illustrates that the UV-edge of O i w.r.t.O j (the solid line) is a hyperbola.Equation 7shows that a UV-cell is composed of the intersections of one or more UV-edges, which are hyperbolas.Since a hyperbola is a conic curve, an UV-edge must be concave in shape.In Fig. 2, apart from the edges of the domain space, the UV-cells of the three objects have concave edges.Note that Eq. 7 has two curves, which represent the UV-edges for each pair of objects involved.For example, in Fig. 3, the solid line is the UV-edge of O i w.r.t.O j , and the dotted line is the UV-edge of O j w.r.t.O i .
If two objects overlap, then dist (c i , c j ) < r i + r j , and in Eq. 7, b is not a real number.Physically, this means E i ( j) cannot be found, and we can treat X i ( j) as an empty region.
We now revisit Algorithm 1.
Step 4 is done using Eq. 7.
Step 5 is performed by observing that the outside region of a UV-edge must be convex in shape.To perform Step 6 (i.e., cutting the possible region by an outside region), we compute the intersections of hyperbola equations using linear algebra techniques [3], which are detailed in Appendix 9.
Let us state an interesting observation about a possible region, which we will use later.

Lemma 1 The possible region of an uncertain object is a connected region without any hole inside it.
The proof of this lemma, detailed in Appendix 10, shows that a contradiction will result if a possible region contains a hole.We next discuss the shape of the UV-cell for other kinds of uncertainty regions.
(1) Point uncertainty.Given two objects O i and O j , suppose that at least one of them has no uncertainty, i.e., r i or r j is equal to zero.There are two scenarios: -If c i = c j , without loss of generality, assume that r i = 0.
Then, E i ( j) can be obtained by Eq. 7, because all variables used in that equation are real numbers, and a, b are nonzero values.Notice that E i ( j) becomes a perpendicular line segment when r i = r j = 0. -If c i = c j , then E i ( j) does not exist.If r i = r j , Eq. 1 does not hold, and the UV-cell of O i , or U i , does not exist.If r i = r j , Eq. 1 always holds, and U i = D.
(2) Non-circular uncertainty regions.To find the UVcells for non-circular uncertainty regions, our first attempt is to derive the UV-edges for objects with rectangular uncertainty regions.As shown in Fig. 5, the UV-edge between objects O 1 and O 2 is a piecewise-quadratic line segments.This is too expensive to compute and store.Instead, for each object O i , we convert its non-circular uncertainty region to a circle, M BC(O i ), which minimally contains it.Then, we use Algorithm 1 to construct the UV-cells for these circles.We claim that M BC(O i )'s UV-cell always covers that of O i .
To understand why, notice that if some object O 1 may be q's nearest neighbor, then M BC(O 1 ) can also be q'snearest Fig. 5 A UV-Edge for rectangular regions neighbor.First, for all i = 1, . . ., n, dist min (q, O 1 ) < dist max (q, O i ).Also, it must also be located in the UV-cell of M BC(O 1 ).In other words, M BC(O i )'s UV-cell always cover that of O i .Therefore, in answering a PNN, if we found that M BC(O i ) contains q, we have to verify whether O i can be the nearest neighbor of q.In the sequel, we assume that all uncertainty regions are circular.

The number of UV-edges of a UV-cell
Let us now examine the number of UV-edges of a UV-cell.As Algorithm 1 shows, for every object O i , its UV-edge with respect to other objects is used to refine its possible region P i (Step 6).This requires computing the intersections of all edges of P i with a new UV-edge E i ( j), for some object O j .
As shown in Fig. 4b, E From this example, we can see that E i ( j), a hyperbolic curve, can have at most 2 intersections with a UV-edge of P i ; and 3 new edges can be created for P i as a result.In the worst case, the number of edges of P i increases by 3 times whenever a new edge is considered.
In general, to obtain U i , we have to take into account n −1 objects.Hence, the number of edges of the UV-cell has an (exponential) upper bound of O(3 n ).
Moreover, computing intersections between hyperbolas is complex.In our implementation, 60 h are needed to create a UV-diagram of 40,000 objects by using Algorithm 1 .We will explain how to find an efficient representation of the UV-cell, in Sect. 5. We now estimate the size of a UV-cell, under the assumption that all objects are evenly placed.We consider the hexagonal lattice model, where each object has six neighbors whose centers are equidistant from each other, with distance d 0 . 1  We assume that the uncertainty region sizes of all objects are the same, with a radius of r .Figure 6 illustrates seven objects configured in this manner.Given an object O 1 , we assume that the UV-cell U 1 of O 1 is not trimmed by the boundary of the domain space.That is, its UV-cell is solely determined by the uncertainty regions of other objects.Our goal is to find the dimension s of a square that contains U 1 .This square should be a good approximation of U 1 . 2et H (d) be a set of six objects whose uncertainty region's centers have the same distance d from that of O 1 .For example, Fig. 6, the centers of the uncertainty regions of H (d) = {O 2 , . . ., O 7 } are d units away from that of O 1 .We claim the following: Lemma 2 Let P 1,d be the possible region of O 1 generated by the objects in H (d). The length of the square, which is centered at c 1 and minimally covers P 1,d , is as follows: In the sequel, we use s(d) to denote the size of P 1,d .The proof of Lemma 2 can be found in Appendix 11.Notice that P 1,d contains U 1 .Now, observe that the centers of the six objects in H (d 0 ) form the vertices of a hexagon called H E X 1 .This hexagon is illustrated in Fig. 7.As shown in [35], a larger hexagon 1 The centers of uncertainty regions form the vertices of n hexagons, each of which has an area of Fig. 7 Illustrating O 1 and its neighbors H E X i+1 , formed by the centers of six other objects, can be obtained by rotating H E X i by π 6 radians, and scaling it by a factor √ 3. Figure 7 shows how H E X 2 and H E X 3 are generated in this way.We then obtain the following result.
3r , then the square that minimally contains U 1 has a size of s(d 0 ) obtained from Eq. 8.
The main idea of the proof is that when d 0 > 2 √ 3r , the six objects that form H E X 1 alone contribute to the edges of O 1 's UV-cell.Its details can be found in Appendix 12.
An iterative approach of finding d * .We now explain how to derive the size of a square that contains U 1 , for any value of d 0 .Our goal is to find d * among different values of d, such that the square covering P 1,d * is the smallest.We observe from the first-order derivative of s from Eq. 8 that d + = 2 √ 3r is the only inflexion point, such that s monotonously decreases when 4r

Efficient UV-cell generation
Since generating a UV-cell is inefficient, our strategy is to avoid computing it directly.Instead, we represent a UV-cell as a set of candidate reference objects (cr-objects), which can be efficiently derived.As will be discussed in Sect.6, cr-objects can be used to develop an approximate representation of the UV-diagram.Section 5.1 outlines the algorithm of yielding cr-objects.We explain the preparation phase of this algorithm as well as two techniques for finding these objects quickly, in Sects.5.2 and 5.3, respectively.Section 5.4 discusses how to derive cr-objects efficiently for a group of nearby objects.

Reference objects and candidate reference objects
Recall from Algorithm 1 that the UV-cell of an object O i , that is, U i , is the result of repeatedly subtracting the outside region of other objects (i.e., X i ( j)) from its possible region, P i .In fact, not all outside regions are useful for refining P i .In particular, if the UV-edge of O i corresponding to O j , that is, E i ( j), does not intersect with P i , then P i cannot be shrinked by X i (j).We call an object O j a reference object (or r-object) of O i , if O j defines an edge of O i 's UV-cell.We also denote F i ⊆ O to be the set of r-objects of O i .The set F i contains objects whose outside regions are responsible for defining the UV-cell of O i .In Fig. 2, for example, the set of r-objects of O 3 , that is, Given that the r-objects for each object are known, our solution (to be shown in Sect.6) can use r-objects to develop an alternative representation of the UV-diagram.This solution is much cheaper than Algorithm 1 , which requires exact UV-cells to be computed.However, finding F i itself is difficult because we do not know the UV-cell of O i .Our strategy is to find a small set C i of objects, where F i ⊆ C i .We call C i the candidate reference objects (or cr-objects in short).We next show how C i can be derived without acquiring the exact UV-cell of O i .In Sect.6, we study an indexing solution based on cr-objects.
Algorithm 2 (getcrObject(O i , S)) presents a procedure that derives the cr-object set C i for object O i , based on a set S of objects.To retrieve C i , we can simply invoke getcrObject(O i , O).In this algorithm, Step 1 (initPossibleRegion) creates a possible region P i based on a small number of objects retrieved from S. In Step 2, the "index level" pruning (or iPrune) yields a set I of objects that may contribute edges to the UV-cell.Step 3 applies "computational level" pruning (or cPrune) on I , and produces C i .Here, we assume that an R-tree index has been built on the uncertainty regions of the objects in O.Each object's information (e.g., uncertainty region and pdf) is stored in the disk.Next, we explain how to generate an initial possible region (Sect.5.2), based on which two techniques for pruning non-cr-objects are developed (Sect.5.3).

Generating a possible region
In Step 1 of Algorithm 2 , we retrieve a small number of objects, called seeds, from a set of objects S.These seeds are used to generate an "initial" possible region, using a routine similar to Steps 3-7 of Algorithm 1.This region is used by other pruning methods to produce cr-objects.
Seeds have to be selected with care.If seeds are randomly selected, a big initial region can be produced.This region may be intersected by many outside regions, resulting in a poor pruning efficiency.Ideally, we would like the initial possible region generated by these seeds to closely resemble the UV-cell.We would also prefer the number of these seeds to be small, so that the possible region can be constructed efficiently.We next present two simple steps to find "good" seeds.
Step (i).We issue a k-Nearest-Neighbor Query (k-NN) on S, by using the center c i of O i 's uncertainty region as the query point.The k objects, which are not O i and whose regions' minimum distances from c i are the shortest, are obtained.Since these objects are close to O i , we consider them to have a good chance for defining the UV-edges of U i .They are thus good candidates for being seeds.Note that if S = O, then the R-tree on O can be used to support the k-NN search.
Step (ii).Out of the k objects obtained from Step (i), we select k s seeds.These objects are chosen in way such that they are evenly spread, in order to generate a "good" possible region.In particular, we divide the domain D into k s equally sized sectors, centered at c i .For each sector, the object closest to c i is a seed.
The above method does not guarantee that all k s seeds can be found (e.g., no seeds can be found if a sector is empty).Even if this happens, however, we can still obtain an initial possible region without affecting the latter steps.This region may be larger, though.In our experiments, k s = 30, and in most cases, all seeds can be found.For each object, evaluating a k-NN query requires O(|S|) time, selecting seeds costs O(k) time, and constructing an initial region needs O(k s ) time.Hence, the cost of this step is Model-based seed selection.We can use the results in Sect.4.4 to estimate the value of k derived from Step (i).We assume that all the objects in domain D follow the hexagonal lattice model.First, we find the size s of the square that bounds the UV-cell of O i .Particularly, we check whether the condition for Theorem 1 is satisfied.If this is true, we let d min = d 0 .Otherwise, we use the iterative approach, described in Sect.4.4, to find d * , and let d min = d * .Then, we find s(d min ) by using Lemma 2. Figure 6 shows that the maximum distance of any point on the possible region P i from the center c i of O i 's uncertainty region is s 2 .If we draw a circle of radius (s − r ), centered at c i , then Theorem 3 (to be discussed in Sect.5.3) tells us that only objects located in this circle can be the reference objects F i .We can then estimate k as the expected number of objects in (c i , s − r ): We can also use the above approach in a database whose locations are uniformly distributed in D. We first compute the average uncertainty radius r a of these objects.We then suppose that all these objects have the same radius r a .We also evaluate the distance of each object from its nearest neighbor, and find the average d a of these distances.The values of the radius r and the distance d 0 of the hexagonal model are set to be r a and d a , respectively.We also compute the density ρ, which is equal to No. of objects inD Area of D . Our experiments show that this model can enhance the seed selection process for uniformly distributed data.

I-Pruning and C-Pruning
Once the possible region has been initialized, we perform I-pruning and C-pruning (Steps 2 and 3 of Algorithm 2), in order to remove objects that cannot constitute a UV-edge to the UV-cell.Let us now examine these two steps in more details.
Step 2: Index Level Pruning.To understand this step, let us consider an object O i , its possible region P i , and another object O j , which has not yet been considered for refining P i .Our goal is to establish the necessary and sufficient condition(s) for O j to have effect on the shape of P i .We first claim the following.

if and only if for every point p inside P i , dist max ( p, O j ) > dist min ( p, O i ).
Proof (If) For every p ∈ P i , p cannot be on X i ( j).If this is false, then O j is always closer to p than O i , i.e., dist max ( p, O j ) ≤ dist min ( p, O i ) (Definition 2).This violates the condition that dist max ( p, O j ) > dist min ( p, O i ).Hence, p / ∈ X i ( j), and P i = P i − X i ( j). (Only if) Suppose there exists a point p inside P i , such that dist max ( p , O j ) ≤ dist min ( p , O i ).Then O j is always closer to p than O i , and O i cannot be the nearest neighbor of p .This implies that p must be excluded from P i after O j is considered.Hence, P i cannot be equal to P i − X i ( j), resulting in a contradiction.Let b(P i ) be the boundary of P i .We have: Proof (If) Let us first show that: Suppose by contrary that the above is not correct.That is, ∃ p If we let P i be P i − X i ( j), then p ∈ X i ( j) and p / ∈ P i .From the given condition, we can see that for every p ∈ b(P i ), p / ∈ X i ( j), and p ∈ P i .Thus, P i must have a hole ( p ) inside it.However, this must not be true, according to Lemma 1.
Hence, Eq. 10 is true.Using Lemma 3, we see that P i = P i − X i ( j), and the so the "if" part is correct.
(Only if) From Lemma 3, we know that for every point p ∈ P i , dist max ( p, O j ) > dist min ( p, O i ).Since b(P i ) ⊆ P i , the "only if" part is correct.
Essentially, if we want to examine whether O j has any effect on P i , it suffices to consider the points on P i 's boundary, instead of all points in P i .We next present the following theorem, which forms the basis of I-pruning.

Theorem 3 Given an object O i with center c i and radius r i , let w be the maximum distance of P i from c i . Let C out be a circle, with center c i and radius
Proof Denote C in by a circle with center c i and radius w. Figure 8(a) illustrates O i , its possible region P i (in solid lines), C in and C out .Let us suppose on the contrary that P i is not equal to P i − X i ( j), that is, P i can be reshaped by the UV-edge of O j .Then, using Theorem 2, there must exist a point p on the boundary of P i such that: Using Eqs. 4 and 5, we have: since dist (c i , c j ) ≤ dist ( p, c j ) + dist ( p, c i ) due to the triangular inequality.Now, dist ( p, c i ) ≤ w, so Eq. 12 becomes: This implies that c j is in the circle C out , contradicting the assumption of Theorem 3. Hence, this lemma is correct.
The I-pruning method uses Theorem 3 by issuing a circular range query, centered at c i with radius 2w −r i , on the dataset O.This operation can be easily implemented by using the R-tree created for O.The range query first uses the R-tree to filter all objects that do not overlap with the range.For the remaining objects, they are removed if their centers are beyond the circular range.Hence, in this phase (Step 2 of Algorithm 2), a cost of O(n) is needed.
Step 3: Computational Level Pruning.We next discuss a method, based on distance comparison, to check whether object O j can affect the possible region of object O i .We call this method C-pruning (Step 3 of Algorithm 2).Theorem 4, discussed below, serves as the foundation of C-pruning.

Theorem 4 Given an uncertain object O i (c i , r i ) and P i 's convex hull C H (P
Proof First, the convex hull C H(P i ), which completely contains P i , must also be O i 's possible region.For every point p on C H(P i )'s boundary, suppose c j is located outside the circle ( p, dist ( p, c i )).Then, we have: ( p, c i )) for every p on C H(P i )'s boundary, O j can be safely pruned.
For convenience, let ( p, dist ( p, c i )) be a w-bound (where w = dist ( p, c i )).We also define a set S of w-bounds for every point p in U i .We now show that instead of checking all the w-bounds in S, it is only necessary to check those w-bounds constructed for the vertices of C H(P i ).Specifically, the w-bounds of the vertices must contain all other w-bounds of all points on the boundary of C H(P i ).To see this, let w k be the distance of vertex v k from O i 's center.We extend each vertex v k by the distance w k to obtain a new vertex v j (black dot in Fig. 8b).These new vertices are connected to form a polygon.We use e 1 and e 2 to represent the w-bounds (v 1 , w 1 ) and (v 2 , w 2 ), respectively.
We next show that, for any point v on C H(P i )'s edge where we let e = (v , dist (v , c i )).We draw a line c 1 c 1 , which is perpendicular with v 1 v 2 and v 1 v 2 , and intersects them at points c 1 and c 1 , respectively.As v 1 v 2 is the perpendicular bisector of c i c 1 , we see that c i c 1 is the common chord of e 1 , e 2 and e .Since e 1 or e 2 is bigger than e , e is contained by e 1 or e 2 .
Hence, to check whether O j can refine P i , we just need to check the set of w-bounds S = { (v m , dist (v m , c i ))} (where S ⊆ S).If c j is located outside all w-bounds in S , then C H(P i ) = C H(P i ) − X i ( j).Finally, since P i is completely covered by C H(P i ), P i = P i − X i ( j) must also be true.This completes the proof.
Step 3 of Algorithm 2 uses Theorem 4 to prune unqualified objects returned by I-pruning.This can be done efficiently, because only the vertices of C H(P i ) are used.Moreover, |C H(P i )| is small, since the possible region is only derived by a small number k s of seeds.The complexity of this phase is O(n).
We consider the objects that remain after this step as cr-objects (i.e., C i ).The complexity of Algorithm 2, for generating C i 's of n objects, is O(n(n + k)).

Batch processing of cr-objects
To create the UV-index, we first find out the cr-objects for each of the n database objects.A simple way to do this is to run Algorithm 2 (i.e., getcrObject(O i , O)) for all objects O i ∈ O, as proposed in [17].However, this involves running getcrObject for n times and can be quite costly.We now present a Batch Processing algorithm (or BP in short), where the cr-objects of a group of objects are considered together.As we will show, this new algorithm allows the effort of devising an object's cr-objects to be shared by others, and consequently reduces a lot of cr-object generation overhead.
Observe that if an object O i is near to O j , then their UVcells should be similar.The cr-object set of O i , that is, C i , can then be similar to C j .The BP makes use of this principle; C i ← cPrune(P i , C G ) 6: end for of all objects in G. Its uncertainty pdf is not important here, and we assume it to be uniform.Notice that O G is only used by the BP; it will be deleted after the algorithm halts.
Step 2 invokes a slightly modified version of getcrObj ect to obtain a cr-object set C G of O G .Particularly, in Step (i) of initPossibleRegion, the k-NN search skips all objects in G. Notice that initPossibleRegion computes the possible region of an object.In Step (i) of that procedure, we obtain the seeds -objects that are useful for generating a UV-cell.In Algorithm 3, the input of getcrObject is O G , whose uncertainty region includes the uncertainty regions of all objects in G. Therefore, the uncertainty region of any object O i ∈ G overlaps with that of O G .More importantly, O i ∈ G is not useful for finding possible regions of O G , because O i does not create any UV-edge for O G 's UV-cell.We next claim the following: That is to say, any object not contained in C G cannot be an r-object of O i ∈ G.In other words, C G is a superset of r-objects for all the objects in G.The proof of this theorem, which is quite complex, is detailed in Sect.5.5.Notice that all objects in G are included in C G after the execution of Step 2. This is because in the last step of getcrObject (Algorithm 2), objects whose centers are located in the cpruning bound of O G are treated as cr-objects.Since the center of an object in G is inside O G 's c-pruning bound, it must also be a cr-object of O G .Thus, G ⊆ C G .
Steps 3-6 use C G to generate cr-objects for each object O i ∈ G. From Theorem 5, we know that an object in C G may be an r-object of O i .Thus, objects in C G can be considered as good candidates for generating an initial possible region, P i for O i .We thus pass C G to initPossibleRegion and get P i (Step 4).We then execute cPrune on C G to retrieve C i (Step 5).These two steps are repeated for all objects in G, until we obtain their cr-objects. 3he LP algorithm.We now discuss a way to use Algorithm 3 to construct cr-objects for O.The Leaf-Node Processing, or LP, performs a traversal of the Rtree that indexes O.When a leaf node, say N , is reached, BP is invoked on all objects stored in N , in order to compute their cr-objects.The algorithm terminates when all leaf nodes have been exhausted.
The LP can generate cr-objects for O quickly.This is because when the BP is called, it always uses the objects located in a leaf node.In an R-tree, the leaf node consists of a set G of objects, which are physically close to each other.Recall that the object created in Step 1 of BP (i.e., O G ) is the MBC of the uncertainty regions of objects in G. Thus, the size of O G would not be very different from those of the objects in G. Consequently, the set C G derived from Step 2 (getcrObject) should also be similar to the r-objects of G's objects.In our experiments, |C G | is much smaller than |O|.Hence, Step 4 can be carried out more efficiently than if initPossibleRegion is carried out on O for every object.
We have introduced an efficient construction method to derive the cr-object set C i for O i .We have also explained how to obtain cr-objects for O quickly.One may consider to use C i to generate the exact UV-cell of O i .However, our experiments show that |C i | may be large, and so generating the UV-cell of O i can still be costly.In Sect.6, we show how to use C i directly to construct an index for the UV-diagram.In the rest of this section, we present the proof of Theorem 5.

Proof of Theorem 5
Recall that O G is formed by a set G of objects, using Step 1 of Algorithm 3. Let P i (S) be a possible region of an object O i , constructed by using a set S ∈ O of objects.Essentially, P i (S) is the intersection of the inside regions X i (k), where O k ∈ S. Let u i be the uncertainty region of O i .We first claim the following.Figure 9 illustrates Lemma 4, which shows that P i (S) is inside P k (S).Its proof can be found in Appendix 13.1.

Lemma 5 Given objects O i and O j , if c j /
∈ P i , then ∀ p ∈ P i (S): if and only if ∃O k ∈ S, where In Fig. 10, the objects in S are shaded.The center of O j , that is, c j , is outside P i (S).Given a point p ∈ P i , Lemma 5 states that if there is an object O k ∈ S such that dist max ( p, O j ) > dist max ( p, O k ), then dist max ( p, O j ) > dist min ( p, O i ), or vice versa.Its proof can be found in Appendix 13.2.These results are used by the next lemma.

Lemma 6 Given two objects O i and O k , where u i ⊆ u k , and an object O j where c j /
∈ P k (S), if P k (S) = P k (S) − X k ( j), then P i (S) = P i (S) − X i ( j).
As shown in Fig. 9, Lemma 6 claims that given an object O j whose center is outside P k (S), if the edge E k ( j) does not affect the possible region P k (S), then E i ( j) cannot contribute to P i (S).
Proof Since P k (S) = P k (S) − X k ( j), by using Lemma 3, we have: Using the "only-if" part of Lemma 5, we have: Since u i ∈ u k , using Lemma 4, we have P i (S) ∈ P k (S).Thus, Eq. 18 becomes: 123 Using the "if" part of Lemma 5, Eq.
Using Lemma 3, Eq.20 means that P i = P i − X i ( j).Hence, the lemma is correct.
Proof of Theorem 5 Let V be the set of seeds used to construct possible region P G (V ) in Step 2 of Algorithm 3. If O j / ∈ C G , the UV-edge E G ( j) does not cut P G (V ).In other words, P G (V ) = P G (V ) − X G ( j).The I-pruning and Cpruning methods used in Step 2 also guarantee that c j is not inside P G (V ), that is, c j / ∈ P G (V ).Moreover, u i ⊆ u G .By substituting S = V and k = G, we can deduce from Lemma 6 that P i (V ) = P i (V ) − X i ( j).Now there are two cases to consider: Case 1: O j contributes an edge to P i (V ).In other words, O j ∈ V .Since an object in V is not pruned by Step 2 of Algorithm 3, V ⊆ C G , and so O j ∈ C G .However, this contradicts with the assumption that O j / ∈ C G , and so this case cannot occur.Case 2: O j does not contribute an edge to P i (V ).Since the UV-cell U i of O i must be inside P i (V ), O j cannot contribute an edge to U i .Hence, O j is not an r-object of O i , and the theorem holds.

The UV-index
We now present the UV-index, an approximate version of the UV-diagram.The UV-index can be efficiently computed and stored.It also facilitates efficient query evaluation.Section 6.1 gives an overview of its structure.In Sect.6.2, we discuss how to use this index to support execution of different queries.We explain its construction process in Sect.6.3.

Structure of the UV-index
The UV-index adopts a framework similar to a quad-tree [5], in order to index the irregular and non-overlapping UVpartitions.Figure 11a illustrates this index. 4Each non-leaf node, 16 bytes each, records a pointer to each of its four child nodes, where the square region spanned by each child node is one-fourth of that of its parent.The region covered by the root node is the whole domain D. Each leaf node stores all the objects whose UV-cells overlap with the region defined for the node.To save space, a node's region is not stored, since we can easily derive the dimension of the region based on the level of the node in the tree.Also, due to approximation, a UV-cell that does not overlap with the leaf node's region may be included.However, a UV-cell that truly overlaps with the region will not be excluded.For each leaf node l, we store a linked list of disk pages, which contain tuples <ID, MBC, pointer>, where: -ID is the identity of object O i whose UV-cell may overlap with the region covered by l; -MBC is the circle that minimally bounds the uncertainty region of O i ; and -pointer stores the disk page address of the object.
We assume that all non-leaf nodes are stored in the main memory, and allocate a maximum number of M non-leaf nodes.The leaf nodes, which contain the lists of pages, are stored in the disk.Hence, M controls the amount of main memory to be used to implement the index.Next, we study how to use it to support query evaluation.

Using the UV-index
We now explain how to use the UV-index to support the queries that we described in Sect.3.2.
1.The PNN Query.To find the probabilistic nearest neighbors of q, we first locate the leaf node l, whose region contains q.This can be done easily by finding the grid that contains q in each index level, and traversing the index.We then retrieve the disk pages associated with l, which contains the ID and the MBC values of the objects stored in these pages.Since these objects may have their UV-cells overlap with the region of l, it is possible that q is contained in their UV-cells.Let L be the set of objects associated with l, and A be the answer objects of q.To answer a PNN, we need to retrieve A from L, where A ⊆ L. We use the method described in [15]: from the set of the MBC's of the objects in L, find d minmax , the minimum of the maximum distances of these objects from q.Any object with the minimum distance larger than d minmax is removed, since it cannot have a nonzero qualification probability.For objects that are not filtered, their probabilities are computed and returned as answers.
2. The CPNN Query maintains answers for a "moving" query point, location is periodically reported to the server.Let q 0 be the latest position of q received by the server.Let g 0 be a leaf node in the UV-index, whose region r 0 contains q 0 .We assume that the objects stored in the disk pages associated with g 0 are known.Now, suppose the new location of q, say, q 1 , is received by the server.A straightforward solution is to treat q 1 as a new PNN query, and use the PNN algorithm described above to compute the answers of q 1 .A better way is to check whether q 1 is inside r 0 .If this is true, we simply use the object set associated with g 0 to compute the answer for q 1 .This saves the effort of traversing the UV-index for q 1 .
3. The UV-Partition Query.We append a counter to each leaf node, and record the number of objects at that node.This process could be done after the UV-index is constructed.Then, a range query with range R is issued over the index, in order to find the leaf nodes whose regions overlap with R. For every leaf node whose region r overlaps with R, we compute its density, which is equal to the number of objects associated with r , divided by the area of r .The query then outputs r and its density value.
4. The UV-cell Query.Notice that if an object O i appears in a leaf node g, its UV-cell overlaps with the region of g.Hence, we can return the approximate area and the extent of O i 's UV-cell by scanning the leaf nodes associated with O i , and then summing up the total area of the regions covered by these nodes.This step can be improved by precomputing and storing the area information.For example, we can scan all the leaf-nodes once, and generate a table for each O i with its respective areas.A similar procedure can be used to support the operation of displaying the approximate shape of the UVcell on the user's screen.

Construction of the UV-index
As discussed in Sect.5, a UV-cell can be represented by a set of cr-objects, C i .We now examine how this facilitates the construction of the UV-index.
Framework.Let g be the grid node being examined, and h k (where k = 1, . . ., 4) be the four child nodes of g.We define a variable nonleafnum, which indicates the number of non-leaf nodes allocated to the index and has an initial value of 1. Originally, the root of the grid is a leaf node, whose region covered (root.region) is the domain D.
We use Algorithm 4 (InsertObj) to insert an object O i to the index.This algorithm, whose inputs are C i and node g, is a recursive procedure, where InsertObj (C i , root) is first invoked.In Step  Then, g becomes the parent of four nodes (h k ), which have been previously generated by CheckSplit.The region of each child node h k covers each of the four quarters of the region defined for g.Also, nonleafnum is incremented by a value of 1. Essentially, The information about the UV-cells previously associated with g are now represented by its child nodes, and g becomes a non-leaf node.
Decision on Splitting.When g's pages are full, either O i 's information is inserted to a new page (OVERFLOW) or split into four child nodes (SPLIT).Ideally, the region of the leaf node that covers q is completely covered by a true UV-partition.This guarantees that the set of objects returned by the UV-index is the true answer objects.The UV-index, which contains grids, is just an approximation of the UVdiagram.Apparently, the more the splitting is performed, the closer the index can resemble the actual UV-diagram, and yield better query performance.
In fact, is not always useful.Suppose that associated with 100 UV-cells.Moreover, g.region is completely covered by each of these UV-cells.Then, it is not necessary to redistribute g into four child nodes.If splitting is performed in this case, then the UV-cells associated with each child node are exactly the same.
Thus, more space is wasted to store duplicated information about the UV-cells.This can happen if the corresponding 100 objects of these UV-cells are close to each other.Then, these UV-cells have similar shapes and significant overlapping.To decide whether to split, we define split fraction, θ , as follows: which is the minimum fraction of UV-cells in one of the child nodes h k that are also in g (note that the UV-cells associated with h k must be the subset of the ones attached to g).
A small θ means that the number of UV-cells overlapping with h k .region is small compared with that of g.We now define a splitting condition of a node: where T θ ∈ [0, 1] is called the split threshold.A larger value of T θ implies a higher tendency of splitting.Algorithm 5 (CheckSplit) implements these ideas.Steps 1-3 return NORMAL if the pages of g are not full.Steps 4 and 5 return OVERFLOW if the number of non-leaf nodes allocated is higher than M. In Steps 7-16, we compute the value of θ , by creating four nodes h k (Step 7), and checking the overlap of each UV-cell with h k .region(Steps 11 and 12).If the splitting condition is satisfied (Step 17), then the SPLIT decision is returned, where Algorithm 4 (Steps 18 and 19) will assign the nodes h k to be the child nodes of g.Otherwise, the child nodes are deleted and an OVERFLOW decision is made (Steps 20 and 21).
Overlap Checking.Algorithm 6 tests whether the UV-cell of an object O i overlaps with a grid g's region r .For every crobject O k ∈ C i , if any of their corresponding outside region (X i (k)) totally contains r , then CheckOverlap returns false (Steps 1-3).Otherwise, true is returned (Step 6).To prove the correctness, we use the following lemma: Lemma 7 If region r is totally covered by X i (k), where O k ∈ C i , then r must not overlap with the UV-cell U i .
Proof We would like to show that if there exists an object O k , such that r ⊆ X i (k), then r ∩U i = φ.Let X i ( j) be the region D − X i ( j).Then U i , the UV-cell of O i , can be expressed as the intersection of all regions X i ( j), for all objects in O except O i , that is, Since r ⊆ X i (k), we have Hence, the lemma is correct.
To check whether a region r is contained in X i ( j) (Step 2), a simple way is to generate and test with the UV-edge E i ( j).This can be avoided, by carrying out a simple 4-point test.Observe that r is a square, and the UV-edge of O i with respect to O j is concave in shape.If all its four corner points are confirmed to be in X i ( j), we conclude that r ⊆ X i ( j). Figure 11b shows that the region of g 1 must not overlap with U i , since all the four corners of g are located on the outside region of one of the UV-edges.Checking whether a point is in X i ( j) is easy, because we can simply check whether the point's minimum distance from O i is larger than its maximum distance from O We thus use the four-point test Step 2.
Notice that Algorithm 6 may incorrectly judge that U i overlaps with r .Figure 11(b) shows that U i does not overlap with the region of grid g 2 .However, some corners of g 2 .regionare not contained in the outside regions of two of the UV-edges of U i .If this is true for all UV-edges of U i , then U i would be decided to be associated with g 2 !If this happens, then during query evaluation, O i will be retrieved from g 2 .This increases the execution time since O i is not in g 2 .However, query accuracy is not affected, since we can still detect that O i is not a nearest neighbor of q.In our experiments, this situation is rare, and does not have a significant effect on query evaluation time. Since ) time, mainly for performing splitting and overlap checking with four child nodes.For Algorithm 4, each UV-cell, in the worst case, needs to perform overlap and split tests with M non-leaf nodes.Hence, its total time complexity is O(Mn 2 ).The index has a maximum height of M/4, if, the data distribution is very skewed, and splitting always happens in one single quadrant.However, all non-leaf nodes, 16-byte long, can all be stored in the main memory.Thus, the tree height has little effect on query performance.

Results
We now report the results.Section 7.1 describes the experiment settings.In Sect.7.2, we discuss the results about query performance.Section 7.4 presents the results about UV-index construction.

Setup
We use both synthetic and real datasets in our experiments.For synthetic data, we use Theodoridis et al's data generator 5  to obtain 20, 40, 60, 80, and 100K objects, which are uniformly distributed in a 10K × 10K space.Each object has a circular uncertainty region with a diameter of 40 units, and a Gaussian uncertainty pdf.For each uncertainty pdf, its mean is the center of the circle, and its variance is the square of one-sixth of the uncertainty region's diameter.We represent an uncertainty pdf as 16 histogram bars, where a histogram bar records the probability that the object is in that area.We also use three real datasets of geographical objects in Germany, namely utility, roads, and rrlines, with respective sizes 17, 30, and 36K.We also test the Long Beach (or LB) dataset, which contains 53K objects. 6These objects 5 http://www.rtreeportal.org/software/SpatialDataGenerator.zip. 6http://www.rtreeportal.org/.
are represented as circles before indexing, and has the same uncertainty pdf information as that of the synthetic data.
To compare with R-tree, we use a packed R*-tree [30] to index uncertain objects.The R-tree uses 4K-byte disk pages, and has a fanout of 100.We keep all its non-leaf nodes in the main memory.
For the UV-index, each non-leaf node has four 4-byte pointers to its children.We set M, the number of non-leaf nodes in the main memory, to be 10,000.The leaf nodes of both indexes, as well as the uncertainty information about the objects, are stored in the disk.
For T θ , the splitting threshold used in constructing the UVindex, we have performed a sensitivity test.Under a wide range of T θ , the indexes only have a slight performance difference.For very small values of T θ (e.g., 0.2), however, the adaptive grid tends not to split, and degrades into long linked lists of pages.Here, we set T θ to be 1.We wrote the programs in C++ and tested them on a Core2 Duo 2.66 GHz PC.

Results on query evaluation
We first study the performance of the queries studied in Sect.3.2.We assume that the LP algorithm, presented in Sect.5.4, is used to generate the UV-index.However, as we will discuss later, the different UV-index construction methods described here has little effect on query performance.
1.The PNN Query.We first compare the PNN performance of the UV-index and the R-tree.We present the average results of 50 query points randomly selected in the data domain.We use the numerical integration method of [15] to implement the probability computation of answer objects. 7igure 12a shows the query running time (T q ) for different synthetic datasets, with the number of objects ranging from 20 to 100K.The running times of both queries increase, because with a larger dataset, potentially more objects qualify as query answers, and this increases the time for index retrieval and probability computation.Our method outperforms R-tree in all cases.For example, when |O| = 60K, the UV-index needs about 50 % of the time needed by the R-tree.
To understand why our method performs better, let us examine the traversal time of the UV-index, which is composed of the time costs for visiting non-leaf and leaf nodes.Since its non-leaf traversal time takes little time in all experiments (up to 3.9 μs), we only present the I/O overhead.In Fig. 12b, we compare the I/O performance of the UVindex and the R-tree.The UV-index requires significantly less number of I/Os than the R-tree (e.g., when |O| = 60K, the UV-index consumes about one-fifth of the I/Os needed by  the R-tree).When the R-tree is used to process a PNN query, plenty of leaf nodes needed to be retrieved.For the UV-index, we only need to look for the leaf node that contains the query point.Since the number of disk pages for each leaf node is also small, a high I/O performance can be attained.Also notice that the number of I/Os for the R-tree increases with |O|, whereas that of the UV-index is relatively stable.Figure 12c shows the time components of T q : (1) index traversal; (2) retrieval of objects' pdf; and (3) probability computation.While object retrieval and probability computation costs are similar for both indexes, the R-tree requires a higher index traversal time.This explains the difference in Fig. 12b.In Fig. 12d, we can see that the query time of both indexes increases with uncertainty region size (i.e., the radius of the uncertainty region), since the larger the region, the more probable that the corresponding object is a PNN answer.For real datasets, columns 3 and 4 of Table 2 show that the UV-index consistently attains a higher query performance than the R-tree.Again, this is because the I/O performance of the UV-index is better than that of the R-tree.
2. The UV-Partition and the UV-cell Queries.We now examine the efficiency of our index for answering the UV-partition query on our synthetic dataset.For each size of a query region R, 50 queries are generated, whose centers of R are uniformly distributed in the data domain.We can see from Fig. 12e that the retrieval time of UV-partitions (T q ) increases with the size of R, since more UV-partitions are loaded when R becomes larger.The increase is almost linear, and the query evaluation time is less than 160 ms.We have also examined the performance of the UV-cell queries on the default synthetic dataset.On average, the time for obtaining a UV-cell from the UV-index is 58.46 ms, or equivalently, 4.62 I/Os.Thus, running a UV-cell query costs little time in our experiments.
3. The CPNN query.To generate a CPNN query, we use the CanuMobiSim simulator, 8 which produces a movingpoint trajectory.The movement of a query point follows a random walk model, as detailed in [34].The location of a query point, which changes at a maximum speed of 100 units per second, is reported every second.The default "trajectory length" of a query is 60, that is, each query has 60 location  We examine two algorithms that use the UV-index to support CPNN queries.The first variant, called UV-index-n, is a naïve application of the UV-index: each time a query point is received, the UV-index is consulted once.The second one, called UV-index-e, is the enhanced version of UV-index-n, where the UV-index is only consulted if the current query point is not located in the same grid as the previous one (Sect.6.2). Figure 13a shows the evaluation time of a query over synthetic data of different sizes.As we can see, the query performance of the UV-index is at least 25 % times better than the R-tree.The reason can be explained by Fig. 13b, which shows the number of I/Os required by these methods.We observe that the I/O cost of issuing a CPNN on the UV-index is much lower than that of the R-tree.For example, when |O| = 60k, the query cost of the UV-index algorithms is about 30 % of the R-tree.We also see that UV-index-e performs better than UV-index-n.When the current query point q 1 is located in the grid g that also contains the previous query point q 0 , UV-index-e uses the objects associated with g to answer the PNN at q 1 .Thus, the effort of traversing the UV-index for q 1 can be saved.This saving is quite significant; at |O| = 60k, for instance, the number of I/Os required by UV-index-e is only 66 % of that of UV-index-n.In Fig. 13c, we examine the effect of the query trajectory length.Again, the UV-index-e performs the best among the three access methods.

Storage cost analysis
Next, we compare the sizes of R-tree and UV-index.As mentioned in Sect.7.1, for both indices, we store the non-leaf nodes in the main memory, and the leaf-nodes in the disk.The index size is the sum of the main memory and disk space required.Figure 14a compares the size of the UV-index and the R-tree.The UV-index is larger than the R-tree.While the UV-index consumes less main memory than the R-tree (Fig. 14b), it needs more disk space (Fig. 14c).Although the UV-index has a larger size than the R-tree, the UV-index provides a better query performance.Moreover, the UV-index provides functionalities that are not available by R-tree (e.g., retrieval of UV-partitions).These benefits are provided in the expense of a larger disk cost.Given the low cost of hard disk space nowadays, we believe that the extra disk space required by the UV-index is still justifiable.

Results on UV-Index Construction
We now examine several UV-index construction methods.We first study the following techniques: -Basic: a UV-cell is derived using Algorithm 1, which is then used to build the UV-index; -ICR (I-and C-pruning with Refinement): collect cr-objects through I-and C-pruning (Algorithm 2), com- We assume that the R-tree for uncertain objects is available for use by these methods.Unless stated otherwise, the modelbased seed selection and batch construction methods are not used (their effect will be examined later).For generating initial possible regions (used in IC and ICR), we set k to 300 for performing the k-NN search.Then, the domain D is divided into k s = 30 sectors to obtain the seeds.Figure 15a describes the development time (T c ) of the UVindex for the three methods.Basic increases sharply with the dataset size; handling a 40K dataset requires about 60 h.This is because constructing a UV-cell requires an exponential amount of time and numerous complex hyperbola intersections.For ICR and IC, the use of I-and C-pruning significantly reduces the number of objects examined.Their effects are shown in Figure 15(b), where p c , the pruning ratio, denotes the fraction of objects from O that has been filtered.At |O|=60k, I-pruning and C-pruning achieve a pruning ratio of 98.9 and 99.5 % respectively.Hence, a large portion of objects are removed before being considered for constructing the UV-cell.Next, we examine ICR and IC.
IC versus ICR.As shown in Fig. 15c, IC performs much better than ICR.For example, at |O| = 80K, the construction time of IC is about 10 % of that of ICR.To understand why, we analyze their time components in Fig. 15d, e. Recall the difference between the two methods is that ICR needs to find out the exact r-objects (by constructing an exact UV-cell based on the objects returned by pruning), while IC does not.For ICR, Fig. 15d shows the fraction of the construction time spent on: (1) seeds selection, (2) initial possible region computation, (3) I-and C-pruning, (4) generating r-objects, and (5) indexing UV-cells.For most datasets, ICR spends most of the time to generate exact r-objects, which is very costly.For IC, r-object is not produced (Fig. 15e).Instead, the cr-objects produced by the pruning methods are immediately passed to Algorithm 4 for indexing.The number of cr-objects generated, while larger than that of r-objects, does not increase the indexing time significantly.
In Fig. 15f, the construction time of ICR increases sharply with the objects' uncertainty region sizes.With larger uncertainty regions, it is more likely that these regions overlap with each other, making it harder to prune the objects, so that more time is needed to generate r-objects.On the other hand, IC is relatively insensitive to the change in uncertainty region sizes.
We have also measured the query times between the indexes created by IC and ICR. Figure 15g shows that the UV-index generated by the two methods is highly similar, resulting in a close query performance.The query cost of ICR is about 0.01 I/Os, or 0.13 ms, better than IC.In the sequel, we assume that IC is used.
Model-based index construction.In Sect.5.2, we have demonstrated how to use the UV-cell model (Sect.4.4) to facilitate seed selection for objects whose locations are uniformly distributed.We call the UV-index construction algorithm that employs this method as Model, and the one that does not use it as Non-model.We evaluate these two algorithms on our synthetic datasets.As we can see from Fig. 16a, Model performs better than Non-model in most cases.When |O|=80 k, about 20 % of the index construction time is saved.Figure 16b illustrates that Model is consistently better than Non-model under different uncertainty region sizes.For example, when the radius of an uncertainty region is 80, the time required by Model is about half of that of I C.
To understand why Model performs well, we compare the difference between the size S of a UV-cell estimated by our model, and its "true" size.Again, S is the length of the MBR that tightly bounds the estimated UV-cell.The true size of the UV-cell can be obtained by using Algorithm 1.Based on the vertices of this UV-cell, we obtain its minimum bounding rectangle (MBR).We use the larger length of the two dimensions of this MBR to represent the size of the UV-cell.Figure 16c shows the average size of a UV-cell under different uncertainty region sizes.The UV-cell size increases with the uncertainty region radius, since an object can be in more possible locations.This increases its chance to be a possible nearest neighbor of a query point.In this experiment, our method offers a reasonable estimation of the UV-cell's sizethe estimation error is between 4 and 12 %.This enables the selection of seeds, as well as the index construction algorithm, to be effective.Batch processing.We next examine the performance of LP, which derives cr-object based on groups of data objects (Sect.5.4).We compare LP with single, which generates a cr-object set for each data object separately.We do not use model-based seed selection in these experiments.Figure 17a shows that LP performs better than single on our synthetic datasets.At |O| = 80k, the time cost of LP is about 60 % of that of single.In LP, the cr-object set generation cost is shared among a group of objects.
We also test the performance of single and LP on larger datasets.We use the same synthetic data generator to produce two datasets that contain 0.5M and 1M objects.The 1M dataset occupies 640Mbytes.The new result, illustrated in Fig. 17b, shows that the construction performance of both single and LP increases with the dataset size in a linear manner.For the 1M dataset, LP needs 7.7 h, which is 23 % faster than single.
Figure 17c shows that when LP is used, the seed selection time of single is shortened by more than 80 %.While single generates seeds for every object individually, in LP, the seeds of every object in set G are retrieved from a set of objects C G (Step 2 of Algorithm 3).We can also see that the I-and C-pruning time required by LP is also less than single; when |O| = 60k, the improvement is over 60 %.In single, I-pruning is done for every object; in LP, I-pruning is only done once for every group.The performance gap is more profound when |O| is large, since the same domain is populated with more objects, resulting in more candidates retrieved after I-pruning.
We also examine the effect of the average uncertainty region size on the construction time.As discussed before, the larger this size, the more construction time will be needed.Figure 17d shows that LP is more stable than single.When the uncertainty region size is 60, LP needs more about 60 % time of single; when the size becomes 100, LP is 3.5 times faster than single.In Fig. 17e, we compare the query performance of the UV-indices generated by single and LP.We observe that the number of I/Os required by the two methods is the same.Their probability computation times, not shown here, are also very close.Hence, the query performance of two methods is almost the same.Next, we compare the construction time of the R-tree and UV-index, using single and LP. Figure 17f shows that the construction cost of the R-tree is less than 1 % of that of the UV-index.Hence, the R-tree introduces little overhead to the UV-index construction process.However, it improves the performance of generating the UV-index.For instance, the I-pruning phase can be executed more efficiently with the use of the R-tree.
For real datasets, LP also outperforms single (Fig. 17g).In rrline, for example, LP needs one-third of the time required by single.LP also achieves a high pruning ratio, as shown in Table 2.This explains why LP requires less time to construct a UV-index, compared with single.
Skewness.In Fig. 18, we study the effect of data skewness, by varying the variance (σ ) of the objects' mean positions.We can see that when the data are more skewed (i.e., with a smaller variance), the construction time is higher, because in a dense area where uncertainty regions have high degree of overlap, an object's UV-cell is likely small and associated with many r-objects.The LP algorithm is still more efficient than single.In the most skewed dataset that we tested (σ =1,500), LP is 33.3 % faster than single.Finally, we examine how a skewed distribution of the centers of uncertainty regions can affect our results.We obtain a 60k dataset that follows the zipfian distribution, by using the same generator that produces our uniformly distributed dataset.For the zipfian distribution, the average query I/O costs for IC and ICR are 2.48 and 2.41.Thus, the query performance of ICR is 0.07 I/Os (or 2.8 %) better than IC.Since their time difference is small (around 0.4 ms), we use IC in the rest of the experiments.
Table 3 compares these two distributions in terms of their construction and query performance, by using the batch processing (LP) technique.Observe that the construction time of the zipfian distribution is worse than the uniform distribution.In a skewed dataset, a UV-cell in a very dense area can be determined by many r-objects, and this renders lower pruning efficiency in the construction phase.However, there is only a slight query I/O difference between the two distributions, and the query performance for both distributions is almost the same.
In the same table, we study the difference between single and LP for zipfian distribution.Notice that LP requires about 42 % of time needed by single.This means that our batch processing method improves the construction performance for zipfian distribution significantly.The query performance of the UV-index constructed by LP is also slightly better (0.03 I/Os) than single.

Conclusions
The UV-diagram is a variant of the Voronoi diagram designed for uncertain data.To tackle the complexity of constructing and evaluating a UV-diagram, we introduce the concept of UV-cells and cr-objects.We study the theoretical size of a UV-cell.We propose an adaptive index for the UVdiagram, and develop efficient algorithms for building it.We also present a batch processing algorithm to further reduce the UV-index construction time.Our experiments show that this index efficiently supports PNNs and other UV-diagramrelated queries.
We plan to study the use of the UV-diagram to support other variants of probabilistic NNQs, for example, approximate NNQs [12,13]; monochromatic and bichromatic reverse-nearest-neighbor (RNN) queries [10,27,42]; and k-RNN queries [11].Another interesting problem is to design a UV-diagram such that whenever a query point is located in a UV-cell U i , we can know that the qualification probability of O i is larger than some threshold T .By using this variant of UV-diagram, we can get all the objects with qualification probability larger than T , without computing their actual probabilities.This could be beneficial to queries where a user is only interested in answers with qualification probabilities larger than T .It is also interesting to examine how the UV-diagram can support multi-dimensional data and incremental updates.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Appendix 1: Hyperbolic curve intersection
As discussed in Sect.3.1, a vertex of the UV-cell is the intersection point of two hyperbolic curves.We now outline the procedure of finding this intersection, using the method described in [3].We can represent two hyperbolic curves, C 1 and C 2 , as homogeneous conic equations: which are obtained by substituting x/z into x and y/z into y for the hyperbolas (Eq.7) of C 1 and C 2 .Next, we construct equation C λ : where λ is a real value, and C λ , a linear combination of C 1 and C 2 , is a system of hyperbolas.We then rewrite C λ in the form of ω T H ω = 0, where ω = (x, y, z) T , and Let det (H ) be the determinant of H .Our aim is to find the value(s) of λ that satisfy the characteristic equation det (H ) = 0.The real value of λ, when substituted into Eq.23, ensures that (1) there is at least one intersection between C 1 and C 2 , and (2) C λ becomes a degenerated hyperbola, in the form of two straight lines.Finally, for each of the λ found from the characteristic equation, we obtain at most four roots that simultaneously satisfy C λ and C 1 .Each root represents an intersection point of C 1 and C 2 .1. Connectivity of P i : Recall from Definition 3 that P i is the intersection of a set of inside regions.Since each inside region is a connected region (by Definition 2), P i must also be connected.2. P i cannot contain any hole inside it.Suppose by contradiction that there is a hole h inside P i , such that an arbitrary point q inside h does not have O i as its possible nearest neighbor.Figure 19 illustrates the situation.Since q must be covered by the UV-cell of some other object, let us assume that q is covered by the UV-cell of object O j .Then, dist min (q, O i ) > dist max (q, O j ), or dist (q, c i ) − r i > dist (q, c j ) + r j .
We now draw a straight line, which passes through c i and q, and intersects the boundary of P i at q .We have: dist min (q, O i )+dist (q, q ) > dist max (q, O j ) + dist (q, q ) ⇒ dist (q, c i )−r i +dist (q, q ) > dist (q, c j ) + r j + dist (q, q ) ⇒ (dist (q, c i )+dist (q, q ))−r i > (dist (q, c j ) + dist (q, q )) + r j Since dist (q, c i ) + dist (q, q ) = dist (q , c i ) and dist (q, c j ) + dist (q, q ) > dist (q , c j ), we have: dist (q , c i ) − r i > dist (q , c j ) + r j In other words, dist min (q , O i ) > dist max (q , O j ).Hence, q cannot have O i as its nearest neighbor.However, this is not possible, since q ∈ P i .Therefore, P i cannot have any hole.

Appendix 3: Size of a possible region (Lemma 2)
Here, we explain how to derive the size of a possible region, as shown in Eq. 8, Sect.4. Let us denote the six objects that have the same distance d from O 1 be {O 2 , . . ., O 7 }, as shown in Fig. 6.We consider two UV-edges E 1 (2) and E 1 (3).Let X 0 be the intersection of E 1 (2) and E 1 (3).Using Eq. 1, we have: Let X 1 be the point on O 1 such that dist (X 0 , X 1 ) = dist min (X 0 , O 1 ).Also, let X 2 (X 3 ) be the point on O 2 (O 3 ) whose distance from X 0 is the maximum between X 0 and O 2 (O 3 ).According to Eq. 24, dist (X 0 , X 1 ) = dist (X 0 , X 2 ) = dist (X 0 , X 3 ) (25) Since X 1 , X 2 and X 3 have the same distance to X 0 , they are on a circle centered at X 0 with radius R. Thus, as shown in Fig. 6, X 0 is the center of circle (X 0 , R), which is externally tangent to O 1 on X 1 , and internally tangent to O 2 (O 3 ) on X 2 (X 3 ).Therefore, Now, let the coordinates of c 1 be (c 1 .x,c 1 .y).Since c 2 c 1 X 0 = π 6 (according to [40]), we have 2 ), and X 0 = (c 1 .x,c 1 .y+ R + r ).By substituting them to Eq. 26, we have: Notice that d has to be larger than 4r √ 3 , in order for R to be positive.The dimension of the square s that bounds the possible region P 1,d is then equal to s = 2 × (R + r ).By substituting R with Eq. 27, we can obtain Eq. 8.
12 Appendix 4: Size of a UV-cell (Theorem 1) Here, we establish the condition that the possible region P 1,d 0 , formed by the six objects in H (d 0 ), is exactly the UVcell of O 1 .Recall that the centers of uncertain regions of objects in H (d 0 ), which are the closest to that of O 1 , form the vertices of a hexagon H E X 1 , as shown in Fig. 7. Now, if objects in H (d 0 ) are disregarded, then any of the object O k whose uncertainty region center is a vertex of hexagon H E X 2 must be the nearest neighbor of O 1 .Suppose that the UV-edge E i (k) cannot contribute to P 1,d 0 .Then, as all uncertainty regions are equally spaced and identical, the UVedges of other objects that are further away from H E X 1 and H E X 2 must also not change the shape of P 1,d 0 .Thus, P 1,d 0 becomes the UV-cell of O 1 , i.e., U 1 .
When does E i (k) fail to influence the shape of P 1,d 0 ?First, we calculate the minimum distance between the center of O 1 and E 1 (k), which is equal to √ 3 2 d 0 +r .We compare this with s(d 0 ) 2 , where s(d 0 ) is the size of the square that bounds

Fig. 4 a
Fig. 4 a Before checking E i ( j). b After checking E i ( j)

Fig. 6
Fig. 6 Estimating the size of a UV-cell

Fig. 8
Fig. 8 Our pruning methods it employs C i to derive C j , instead of generating C i and C j independently.Let G be a set of objects that are physically close to each other.The BP first computes a set of objects C G , a superset of C i , for every O i ∈ G.The cr-objects of objects in G are then extracted from C G .Usually, C G is smaller than the database size |O|, and thus retrieving cr-objects from C G is faster than from O. Algorithm 3 presents the BP.Given G ⊆ O, Step 1 creates a new object O G .The uncertainty region of O G is the minimum bounding circle (MBC) of the uncertainty regions Algorithm 3 BP Input: A set G of objects in O Output: cr-object set C i for each

Fig. 17
Fig. 17 Results on batch processing

Fig. 19
Fig. 19 Illustrating the proof of Lemma 1

Table 1
Notations and meanings

Algorithm 1 Generating a UV-cell Input: Uncertain
objects O = {O 1 , O 2 , . . ., O n } Output: U 1 , U 2 , . . ., U n 1, CheckOverlap investigates whether the UV-cell represented by C i overlaps with the region of grid g.If so, we check whether g is a non-leaf node.If this is true, InsertObj is called recursively (Steps i , ptr(O i )) is inserted to g's page, where ptr(O i ) is the pointer to O i 's uncertainty region and pdf. 2. OVERFLOW (Steps 12-15): g's pages are full, and a new disk page has to be associated with g, before the information about O i is inserted to the new page.3. SPLIT (Steps 16-22): g's pages are full.The page list g is removed.