#include "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(CoordT x, CoordT 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(CoordT x, CoordT 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; } CoordPointT PointU2PointD(m2::PointU const & pt, uint32_t coordBits) { return CoordPointT( static_cast(pt.x) * (MercatorBounds::maxX - MercatorBounds::minX) / CoordSize(coordBits) + MercatorBounds::minX, static_cast(pt.y) * (MercatorBounds::maxY - MercatorBounds::minY) / CoordSize(coordBits) + MercatorBounds::minY); } CoordPointT 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) { CoordPointT const pt1 = Int64ToPoint(p.first, coordBits); CoordPointT const pt2 = Int64ToPoint(p.second, coordBits); return m2::RectD(m2::PointD(pt1.first, pt1.second), m2::PointD(pt2.first, pt2.second)); }