1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
|
#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<uint32_t>(0.5 + (x - MercatorBounds::minX)
/ (MercatorBounds::maxX - MercatorBounds::minX) * CoordSize(coordBits));
uint32_t const iy = static_cast<uint32_t>(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<int64_t>(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<double>(pt.x) * (MercatorBounds::maxX - MercatorBounds::minX)
/ CoordSize(coordBits) + MercatorBounds::minX,
static_cast<double>(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<uint64_t>(v)), coordBits);
}
pair<int64_t, int64_t> 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<int64_t, int64_t> 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);
}
|