#include "indexer/cell_id.hpp" #include "geometry/pointu_to_uint64.hpp" #include "base/bits.hpp" #include "std/algorithm.hpp" namespace { inline double CoordSize(uint32_t coordBits) { return (1 << coordBits) - 1; } } m2::PointU PointD2PointU(double x, double y, uint32_t coordBits) { x = my::clamp(x, MercatorBounds::minX, MercatorBounds::maxX); y = my::clamp(y, MercatorBounds::minY, MercatorBounds::maxY); uint32_t const ix = static_cast(0.5 + (x - MercatorBounds::minX) / (MercatorBounds::maxX - MercatorBounds::minX) * CoordSize(coordBits)); uint32_t const iy = static_cast(0.5 + (y - MercatorBounds::minY) / (MercatorBounds::maxY - MercatorBounds::minY) * CoordSize(coordBits)); ASSERT_LESS_OR_EQUAL(ix, CoordSize(coordBits), ()); ASSERT_LESS_OR_EQUAL(iy, CoordSize(coordBits), ()); return m2::PointU(ix, iy); } int64_t PointToInt64(double x, double y, uint32_t coordBits) { int64_t const res = static_cast(m2::PointUToUint64(PointD2PointU(x, y, coordBits))); ASSERT_LESS_OR_EQUAL(res, 3ULL << 2 * POINT_COORD_BITS, ()); ASSERT_GREATER_OR_EQUAL(res, 0, ("Highest bits of (ix, iy) are not used, so res should be > 0.")); return res; } m2::PointD PointU2PointD(m2::PointU const & pt, uint32_t coordBits) { return m2::PointD( static_cast(pt.x) * (MercatorBounds::maxX - MercatorBounds::minX) / CoordSize(coordBits) + MercatorBounds::minX, static_cast(pt.y) * (MercatorBounds::maxY - MercatorBounds::minY) / CoordSize(coordBits) + MercatorBounds::minY); } m2::PointD Int64ToPoint(int64_t v, uint32_t coordBits) { ASSERT_LESS_OR_EQUAL(v, 3ULL << 2 * POINT_COORD_BITS, ()); return PointU2PointD(m2::Uint64ToPointU(static_cast(v)), coordBits); } pair RectToInt64(m2::RectD const & r, uint32_t coordBits) { int64_t const p1 = PointToInt64(r.minX(), r.minY(), coordBits); int64_t const p2 = PointToInt64(r.maxX(), r.maxY(), coordBits); return make_pair(p1, p2); } m2::RectD Int64ToRect(pair const & p, uint32_t coordBits) { m2::PointD const pt1 = Int64ToPoint(p.first, coordBits); m2::PointD const pt2 = Int64ToPoint(p.second, coordBits); return m2::RectD(pt1, pt2); }