OLD | NEW |
1 // Copyright 2014 The Chromium Authors. All rights reserved. | 1 // Copyright 2014 The Chromium Authors. All rights reserved. |
2 // Use of this source code is governed by a BSD-style license that can be | 2 // Use of this source code is governed by a BSD-style license that can be |
3 // found in the LICENSE file. | 3 // found in the LICENSE file. |
4 | 4 |
5 #ifndef UI_GFX_GEOMETRY_R_TREE_H_ | 5 // Defines a hierarchical bounding rectangle data structure for Rect objects, |
6 #define UI_GFX_GEOMETRY_R_TREE_H_ | |
7 | |
8 #include <vector> | |
9 | |
10 #include "base/containers/hash_tables.h" | |
11 #include "base/gtest_prod_util.h" | |
12 #include "base/macros.h" | |
13 #include "base/memory/scoped_ptr.h" | |
14 #include "base/memory/scoped_vector.h" | |
15 #include "ui/gfx/geometry/rect.h" | |
16 #include "ui/gfx/gfx_export.h" | |
17 | |
18 namespace gfx { | |
19 | |
20 // Defines a heirarchical bounding rectangle data structure for Rect objects, | |
21 // associated with a generic unique key K for efficient spatial queries. The | 6 // associated with a generic unique key K for efficient spatial queries. The |
22 // R*-tree algorithm is used to build the trees. Based on the papers: | 7 // R*-tree algorithm is used to build the trees. Based on the papers: |
23 // | 8 // |
24 // A Guttman 'R-trees: a dynamic index structure for spatial searching', Proc | 9 // A Guttman 'R-trees: a dynamic index structure for spatial searching', Proc |
25 // ACM SIGMOD Int Conf on Management of Data, 47-57, 1984 | 10 // ACM SIGMOD Int Conf on Management of Data, 47-57, 1984 |
26 // | 11 // |
27 // N Beckmann, H-P Kriegel, R Schneider, B Seeger 'The R*-tree: an efficient and | 12 // N Beckmann, H-P Kriegel, R Schneider, B Seeger 'The R*-tree: an efficient and |
28 // robust access method for points and rectangles', Proc ACM SIGMOD Int Conf on | 13 // robust access method for points and rectangles', Proc ACM SIGMOD Int Conf on |
29 // Management of Data, 322-331, 1990 | 14 // Management of Data, 322-331, 1990 |
30 class GFX_EXPORT RTree { | 15 |
| 16 #ifndef UI_GFX_GEOMETRY_R_TREE_H_ |
| 17 #define UI_GFX_GEOMETRY_R_TREE_H_ |
| 18 |
| 19 #include "r_tree_base.h" |
| 20 |
| 21 namespace gfx { |
| 22 |
| 23 template <typename Key> |
| 24 class RTree : public RTreeBase { |
31 public: | 25 public: |
| 26 typedef base::hash_set<Key> Matches; |
| 27 |
| 28 // RTrees organize pairs of keys and rectangles in a hierarchical bounding |
| 29 // box structure. This allows for queries of the tree within logarithmic time. |
| 30 // |min_children| and |max_children| allows for adjustment of the average size |
| 31 // of the nodes within RTree, which adjusts the base of the logarithm in the |
| 32 // algorithm runtime. Some parts of insertion and deletion are polynomial |
| 33 // in the size of the individual node, so the trade-off with larger nodes is |
| 34 // potentially faster queries but slower insertions and deletions. Generally |
| 35 // it is worth considering how much overlap between rectangles of different |
| 36 // keys will occur in the tree, and trying to set |max_children| as a |
| 37 // reasonable upper bound to the number of overlapping rectangles expected. |
| 38 // Then |min_children| can bet set to a quantity slightly less than half of |
| 39 // that. |
32 RTree(size_t min_children, size_t max_children); | 40 RTree(size_t min_children, size_t max_children); |
33 ~RTree(); | 41 ~RTree(); |
34 | 42 |
35 // Insert a new rect into the tree, associated with provided key. Note that if | 43 // Insert a new rect into the tree, associated with provided key. Note that if |
36 // rect is empty, this rect will not actually be inserted. Duplicate keys | 44 // |rect| is empty, the |key| will not actually be inserted. Duplicate keys |
37 // overwrite old entries. | 45 // overwrite old entries. |
38 void Insert(const Rect& rect, intptr_t key); | 46 void Insert(const Rect& rect, Key key); |
39 | 47 |
40 // If present, remove the supplied key from the tree. | 48 // If present, remove the supplied |key| from the tree. |
41 void Remove(intptr_t key); | 49 void Remove(Key key); |
42 | 50 |
43 // Fills a supplied list matches_out with all keys having bounding rects | 51 // Fills |matches_out| with all keys having bounding rects intersecting |
44 // intersecting query_rect. | 52 // |query_rect|. |
45 void Query(const Rect& query_rect, | 53 void AppendIntersectingRecords(const Rect& query_rect, |
46 base::hash_set<intptr_t>* matches_out) const; | 54 Matches* matches_out) const; |
47 | 55 |
48 // Removes all objects from the tree. | |
49 void Clear(); | 56 void Clear(); |
50 | 57 |
51 private: | 58 private: |
52 // Private data structure class for storing internal nodes or leaves with keys | 59 friend class RTreeTest; |
53 // of R-Trees. Note that "leaf" nodes can still have children, the children | 60 friend class RTreeNodeTest; |
54 // will all be Nodes with non-NULL record pointers. | 61 |
55 class GFX_EXPORT Node { | 62 class Record : public RecordBase { |
56 public: | 63 public: |
57 // Level counts from -1 for "record" Nodes, that is Nodes that wrap key | 64 Record(const Rect& rect, const Key& key); |
58 // values, to 0 for leaf Nodes, that is Nodes that only have record | 65 virtual ~Record(); |
59 // children, up to the root Node, which has level equal to the height of the | 66 const Key& key() const { return key_; } |
60 // tree. | |
61 explicit Node(int level); | |
62 | |
63 // Builds a new record Node. | |
64 Node(const Rect& rect, intptr_t key); | |
65 | |
66 virtual ~Node(); | |
67 | |
68 // Deletes all children and any attached record. | |
69 void Clear(); | |
70 | |
71 // Recursive call to build a list of rects that intersect the query_rect. | |
72 void Query(const Rect& query_rect, | |
73 base::hash_set<intptr_t>* matches_out) const; | |
74 | |
75 // Recompute our bounds by taking the union of all children rects. Will then | |
76 // call RecomputeBounds() on our parent for recursive bounds recalculation | |
77 // up to the root. | |
78 void RecomputeBounds(); | |
79 | |
80 // Removes number_to_remove nodes from this Node, and appends them to the | |
81 // supplied list. Does not repair bounds upon completion. | |
82 void RemoveNodesForReinsert(size_t number_to_remove, | |
83 ScopedVector<Node>* nodes); | |
84 | |
85 // Given a pointer to a child node within this Node, remove it from our | |
86 // list. If that child had any children, append them to the supplied orphan | |
87 // list. Returns the new count of this node after removal. Does not | |
88 // recompute bounds, as this node itself may be removed if it now has too | |
89 // few children. | |
90 size_t RemoveChild(Node* child_node, ScopedVector<Node>* orphans); | |
91 | |
92 // Does what it says on the tin. Returns NULL if no children. Does not | |
93 // recompute bounds. | |
94 scoped_ptr<Node> RemoveAndReturnLastChild(); | |
95 | |
96 // Given a node, returns the best fit node for insertion of that node at | |
97 // the nodes level(). | |
98 Node* ChooseSubtree(Node* node); | |
99 | |
100 // Adds the provided node to this Node. Returns the new count of records | |
101 // stored in the Node. Will recompute the bounds of this node after | |
102 // addition. | |
103 size_t AddChild(Node* node); | |
104 | |
105 // Returns a sibling to this Node with at least min_children and no greater | |
106 // than max_children of this Node's children assigned to it, and having the | |
107 // same parent. Bounds will be valid on both Nodes after this call. | |
108 Node* Split(size_t min_children, size_t max_children); | |
109 | |
110 // For record nodes only, to support re-insert, allows setting the rect. | |
111 void SetRect(const Rect& rect); | |
112 | |
113 // Returns a pointer to the parent of this Node, or NULL if no parent. | |
114 Node* parent() const { return parent_; } | |
115 | |
116 // 0 level() would mean that this Node is a leaf. 1 would mean that this | |
117 // Node has children that are leaves. Calling level() on root_ returns the | |
118 // height of the tree - 1. A level of -1 means that this is a Record node. | |
119 int level() const { return level_; } | |
120 | |
121 const Rect& rect() const { return rect_; } | |
122 | |
123 size_t count() const { return children_.size(); } | |
124 | |
125 intptr_t key() const { return key_; } | |
126 | 67 |
127 private: | 68 private: |
128 // Returns all records stored in this node and its children. | 69 Key key_; |
129 void GetAllValues(base::hash_set<intptr_t>* matches_out) const; | |
130 | 70 |
131 // Used for sorting Nodes along vertical and horizontal axes | 71 DISALLOW_COPY_AND_ASSIGN(Record); |
132 static bool CompareVertical(Node* a, Node* b); | |
133 | |
134 static bool CompareHorizontal(Node* a, Node* b); | |
135 | |
136 static bool CompareCenterDistanceFromParent(Node* a, Node* b); | |
137 | |
138 // Returns the increase in overlap value, as defined in Beckmann et al as | |
139 // the sum of the areas of the intersection of each |children_| rectangle | |
140 // (excepting the candidate child) with the argument rectangle. The | |
141 // expanded_rect argument is computed as the union of the candidate child | |
142 // rect and the argument rect, and is included here to avoid recomputation. | |
143 // Here the candidate child is indicated by index in |children_|, and | |
144 // expanded_rect is the alread-computed union of candidate's rect and | |
145 // rect. | |
146 int OverlapIncreaseToAdd(const Rect& rect, | |
147 size_t candidate, | |
148 const Rect& expanded_rect); | |
149 | |
150 // Bounds recomputation without calling parents to do the same. | |
151 void RecomputeBoundsNoParents(); | |
152 | |
153 // Split() helper methods. | |
154 // | |
155 // Given two vectors of Nodes sorted by vertical or horizontal bounds, this | |
156 // function populates two vectors of Rectangles in which the ith element is | |
157 // the Union of all bounding rectangles [0,i] in the associated sorted array | |
158 // of Nodes. | |
159 static void BuildLowBounds(const std::vector<Node*>& vertical_sort, | |
160 const std::vector<Node*>& horizontal_sort, | |
161 std::vector<Rect>* vertical_bounds, | |
162 std::vector<Rect>* horizontal_bounds); | |
163 | |
164 // Given two vectors of Nodes sorted by vertical or horizontal bounds, this | |
165 // function populates two vectors of Rectangles in which the ith element is | |
166 // the Union of all bounding rectangles [i, count()) in the associated | |
167 // sorted array of Nodes. | |
168 static void BuildHighBounds(const std::vector<Node*>& vertical_sort, | |
169 const std::vector<Node*>& horizontal_sort, | |
170 std::vector<Rect>* vertical_bounds, | |
171 std::vector<Rect>* horizontal_bounds); | |
172 | |
173 // Returns true if this is a vertical split, false if a horizontal one. | |
174 // Based on ChooseSplitAxis algorithm in Beckmann et al. Chooses the axis | |
175 // with the lowest sum of margin values of bounding boxes. | |
176 static bool ChooseSplitAxis(const std::vector<Rect>& low_vertical_bounds, | |
177 const std::vector<Rect>& high_vertical_bounds, | |
178 const std::vector<Rect>& low_horizontal_bounds, | |
179 const std::vector<Rect>& high_horizontal_bounds, | |
180 size_t min_children, | |
181 size_t max_children); | |
182 | |
183 // Used by SplitNode to calculate optimal index of split, after determining | |
184 // along which axis to sort and split the children rectangles. Returns the | |
185 // index to the first element in the split children as sorted by the bounds | |
186 // vectors. | |
187 static size_t ChooseSplitIndex(size_t min_children, | |
188 size_t max_children, | |
189 const std::vector<Rect>& low_bounds, | |
190 const std::vector<Rect>& high_bounds); | |
191 | |
192 // Takes our children_ and divides them into a new node, starting at index | |
193 // split_index in sorted_children. | |
194 Node* DivideChildren(const std::vector<Rect>& low_bounds, | |
195 const std::vector<Rect>& high_bounds, | |
196 const std::vector<Node*>& sorted_children, | |
197 size_t split_index); | |
198 | |
199 // Returns a pointer to the child node that will result in the least overlap | |
200 // increase with the addition of node_rect, as defined in the Beckmann et al | |
201 // paper, or NULL if there's a tie found. Requires a precomputed vector of | |
202 // expanded rectangles where the ith rectangle in the vector is the union of | |
203 // |children_|[i] and node_rect. | |
204 Node* LeastOverlapIncrease(const Rect& node_rect, | |
205 const std::vector<Rect>& expanded_rects); | |
206 | |
207 // Returns a pointer to the child node that will result in the least area | |
208 // enlargement if the argument node rectangle were to be added to that | |
209 // nodes' bounding box. Requires a precomputed vector of expanded rectangles | |
210 // where the ith rectangle in the vector is the union of |children_|[i] and | |
211 // node_rect. | |
212 Node* LeastAreaEnlargement(const Rect& node_rect, | |
213 const std::vector<Rect>& expanded_rects); | |
214 | |
215 // This Node's bounding rectangle. | |
216 Rect rect_; | |
217 | |
218 // The height of the node in the tree, counting from -1 at the record node | |
219 // to 0 at the leaf up to the root node which has level equal to the height | |
220 // of the tree. | |
221 int level_; | |
222 | |
223 // Pointers to all of our children Nodes. | |
224 ScopedVector<Node> children_; | |
225 | |
226 // A weak pointer to our parent Node in the RTree. The root node will have a | |
227 // NULL value for |parent_|. | |
228 Node* parent_; | |
229 | |
230 // If this is a record Node, then |key_| will be non-NULL and will contain | |
231 // the key data. Otherwise, NULL. | |
232 intptr_t key_; | |
233 | |
234 friend class RTreeTest; | |
235 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeBuildHighBounds); | |
236 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeBuildLowBounds); | |
237 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeChooseSplitAxisAndIndex); | |
238 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeChooseSubtree); | |
239 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeCompareCenterDistanceFromParent); | |
240 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeCompareHorizontal); | |
241 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeCompareVertical); | |
242 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeDivideChildren); | |
243 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeLeastAreaEnlargement); | |
244 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeLeastOverlapIncrease); | |
245 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeOverlapIncreaseToAdd); | |
246 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveAndReturnLastChild); | |
247 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveChildNoOrphans); | |
248 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveChildOrphans); | |
249 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveNodesForReinsert); | |
250 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeSplit); | |
251 | |
252 DISALLOW_COPY_AND_ASSIGN(Node); | |
253 }; | 72 }; |
254 | 73 |
255 // Supports re-insertion of Nodes based on the strategies outlined in | |
256 // Beckmann et al. | |
257 void InsertNode(Node* node, int* highset_reinsert_level); | |
258 | |
259 // Supports removal of nodes for tree without deletion. | |
260 void RemoveNode(Node* node); | |
261 | |
262 // A pointer to the root node in the RTree. | |
263 scoped_ptr<Node> root_; | |
264 | |
265 // The parameters used to define the shape of the RTree. | |
266 size_t min_children_; | |
267 size_t max_children_; | |
268 | |
269 // A map of supplied keys to their Node representation within the RTree, for | 74 // A map of supplied keys to their Node representation within the RTree, for |
270 // efficient retrieval of keys without requiring a bounding rect. | 75 // efficient retrieval of keys without requiring a bounding rect. |
271 base::hash_map<intptr_t, Node*> record_map_; | 76 typedef base::hash_map<Key, Record*> RecordMap; |
272 | 77 RecordMap record_map_; |
273 friend class RTreeTest; | |
274 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeBuildHighBounds); | |
275 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeBuildLowBounds); | |
276 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeChooseSplitAxisAndIndex); | |
277 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeChooseSubtree); | |
278 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeCompareCenterDistanceFromParent); | |
279 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeCompareHorizontal); | |
280 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeCompareVertical); | |
281 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeDivideChildren); | |
282 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeLeastAreaEnlargement); | |
283 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeLeastOverlapIncrease); | |
284 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeOverlapIncreaseToAdd); | |
285 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveAndReturnLastChild); | |
286 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveChildNoOrphans); | |
287 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveChildOrphans); | |
288 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeRemoveNodesForReinsert); | |
289 FRIEND_TEST_ALL_PREFIXES(RTreeTest, NodeSplit); | |
290 | 78 |
291 DISALLOW_COPY_AND_ASSIGN(RTree); | 79 DISALLOW_COPY_AND_ASSIGN(RTree); |
292 }; | 80 }; |
293 | 81 |
| 82 template <typename Key> |
| 83 RTree<Key>::RTree(size_t min_children, size_t max_children) |
| 84 : RTreeBase(min_children, max_children) { |
| 85 } |
| 86 |
| 87 template <typename Key> |
| 88 RTree<Key>::~RTree() { |
| 89 } |
| 90 |
| 91 template <typename Key> |
| 92 void RTree<Key>::Insert(const Rect& rect, Key key) { |
| 93 scoped_ptr<NodeBase> record; |
| 94 // Check if this key is already present in the tree. |
| 95 typename RecordMap::iterator it(record_map_.find(key)); |
| 96 |
| 97 if (it != record_map_.end()) { |
| 98 // We will re-use this node structure, regardless of re-insert or return. |
| 99 Record* existing_record = it->second; |
| 100 // If the new rect and the current rect are identical we can skip the rest |
| 101 // of Insert() as nothing has changed. |
| 102 if (existing_record->rect() == rect) |
| 103 return; |
| 104 |
| 105 // Remove the node from the tree in its current position. |
| 106 record = RemoveNode(existing_record); |
| 107 |
| 108 // If we are replacing this key with an empty rectangle we just remove the |
| 109 // old node from the list and return, thus preventing insertion of empty |
| 110 // rectangles into our spatial database. |
| 111 if (rect.IsEmpty()) { |
| 112 record_map_.erase(it); |
| 113 return; |
| 114 } |
| 115 |
| 116 // Reset the rectangle to the new value. |
| 117 record->set_rect(rect); |
| 118 } else { |
| 119 if (rect.IsEmpty()) |
| 120 return; |
| 121 |
| 122 record.reset(new Record(rect, key)); |
| 123 record_map_.insert(std::make_pair(key, static_cast<Record*>(record.get()))); |
| 124 } |
| 125 |
| 126 int highest_reinsert_level = -1; |
| 127 InsertNode(record.Pass(), &highest_reinsert_level); |
| 128 } |
| 129 |
| 130 template <typename Key> |
| 131 void RTree<Key>::Clear() { |
| 132 record_map_.clear(); |
| 133 ResetRoot(); |
| 134 } |
| 135 |
| 136 template <typename Key> |
| 137 void RTree<Key>::Remove(Key key) { |
| 138 // Search the map for the leaf parent that has the provided record. |
| 139 typename RecordMap::iterator it = record_map_.find(key); |
| 140 if (it == record_map_.end()) |
| 141 return; |
| 142 |
| 143 Record* record = it->second; |
| 144 record_map_.erase(it); |
| 145 RemoveNode(record); |
| 146 |
| 147 // Lastly check the root. If it has only one non-leaf child, delete it and |
| 148 // replace it with its child. |
| 149 PruneRootIfNecessary(); |
| 150 } |
| 151 |
| 152 template <typename Key> |
| 153 void RTree<Key>::AppendIntersectingRecords( |
| 154 const Rect& query_rect, Matches* matches_out) const { |
| 155 RTreeBase::Records matching_records; |
| 156 root()->AppendIntersectingRecords(query_rect, &matching_records); |
| 157 for (RTreeBase::Records::const_iterator it = matching_records.begin(); |
| 158 it != matching_records.end(); |
| 159 ++it) { |
| 160 const Record* record = static_cast<const Record*>(*it); |
| 161 matches_out->insert(record->key()); |
| 162 } |
| 163 } |
| 164 |
| 165 |
| 166 // RTree::Record -------------------------------------------------------------- |
| 167 |
| 168 template <typename Key> |
| 169 RTree<Key>::Record::Record(const Rect& rect, const Key& key) |
| 170 : RecordBase(rect), |
| 171 key_(key) { |
| 172 } |
| 173 |
| 174 template <typename Key> |
| 175 RTree<Key>::Record::~Record() { |
| 176 } |
| 177 |
294 } // namespace gfx | 178 } // namespace gfx |
295 | 179 |
296 #endif // UI_GFX_GEOMETRY_R_TREE_H_ | 180 #endif // UI_GFX_GEOMETRY_R_TREE_H_ |
OLD | NEW |