From 700180a0fecc030ecb8b72e9f9ce1b4bb798b45c Mon Sep 17 00:00:00 2001
From: Ansis Brammanis <brammanis@gmail.com>
Date: Thu, 28 Apr 2016 18:00:39 -0700
Subject: [core] use a GridIndex for queryRenderedFeatures

---
 src/mbgl/util/grid_index.cpp | 82 ++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 82 insertions(+)
 create mode 100644 src/mbgl/util/grid_index.cpp

(limited to 'src/mbgl/util/grid_index.cpp')

diff --git a/src/mbgl/util/grid_index.cpp b/src/mbgl/util/grid_index.cpp
new file mode 100644
index 0000000000..8a2ee90213
--- /dev/null
+++ b/src/mbgl/util/grid_index.cpp
@@ -0,0 +1,82 @@
+#include <mbgl/util/grid_index.hpp>
+#include <mbgl/geometry/feature_index.hpp>
+
+#include <unordered_set>
+
+namespace mbgl {
+
+
+template <class T>
+GridIndex<T>::GridIndex(int32_t extent_, int32_t n_, int32_t padding_) :
+    extent(extent_),
+    n(n_),
+    padding(padding_),
+    d(n + 2 * padding),
+    scale(double(n) / double(extent)),
+    min(-double(padding) / n * extent),
+    max(extent + double(padding) / n * extent)
+    {
+        cells.resize(d * d);
+    };
+
+template <class T>
+void GridIndex<T>::insert(T&& t, BBox&& bbox) {
+    size_t uid = elements.size();
+
+    auto cx1 = convertToCellCoord(bbox.x1);
+    auto cy1 = convertToCellCoord(bbox.y1);
+    auto cx2 = convertToCellCoord(bbox.x2);
+    auto cy2 = convertToCellCoord(bbox.y2);
+
+    for (int32_t x = cx1; x <= cx2; x++) {
+        for (int32_t y = cy1; y <= cy2; y++) {
+            auto cellIndex = d * y + x;
+            cells[cellIndex].push_back(uid);
+        }
+    }
+
+    elements.emplace_back(t, bbox);
+}
+
+template <class T>
+std::vector<T> GridIndex<T>::query(const BBox& queryBBox) const {
+    std::vector<T> result;
+    std::unordered_set<size_t> seenUids;
+
+    auto cx1 = convertToCellCoord(queryBBox.x1);
+    auto cy1 = convertToCellCoord(queryBBox.y1);
+    auto cx2 = convertToCellCoord(queryBBox.x2);
+    auto cy2 = convertToCellCoord(queryBBox.y2);
+
+    for (int32_t x = cx1; x <= cx2; x++) {
+        for (int32_t y = cy1; y <= cy2; y++) {
+            auto cellIndex = d * y + x;
+            for (auto uid : cells[cellIndex]) {
+                if (seenUids.count(uid) == 0) {
+                    seenUids.insert(uid);
+
+                    auto& pair = elements.at(uid);
+                    auto& bbox = pair.second;
+                    if (queryBBox.x1 <= bbox.x2 &&
+                        queryBBox.y1 <= bbox.y2 &&
+                        queryBBox.x2 >= bbox.x1 &&
+                        queryBBox.y2 >= bbox.y1) {
+
+                        result.push_back(pair.first);
+                    }
+                }
+            }
+        }
+    }
+
+    return result;
+}
+
+
+template <class T>
+int32_t GridIndex<T>::convertToCellCoord(int32_t x) const {
+    return std::max(0.0, std::min(d - 1.0, std::floor(x * scale) + padding));
+}
+
+template class GridIndex<IndexedSubfeature>;
+}
-- 
cgit v1.2.1