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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
#include "base/SRC_FIRST.hpp"
#include "testing/testing.hpp"
#include "geometry/mercator.hpp"
#include "base/math.hpp"
#include "base/macros.hpp"
#include "base/logging.hpp"
UNIT_TEST(Mercator_Grid)
{
for (int lat = -85; lat <= 85; ++lat)
{
for (int lon = -180; lon <= 180; ++lon)
{
double const x = MercatorBounds::LonToX(lon);
double const y = MercatorBounds::LatToY(lat);
double const lat1 = MercatorBounds::YToLat(y);
double const lon1 = MercatorBounds::XToLon(x);
// Normal assumption for any projection.
TEST_ALMOST_EQUAL_ULPS(static_cast<double>(lat), lat1, ());
TEST_ALMOST_EQUAL_ULPS(static_cast<double>(lon), lon1, ());
// x is actually lon unmodified.
TEST_ALMOST_EQUAL_ULPS(x, static_cast<double>(lon), ());
}
}
}
UNIT_TEST(Mercator_DirectInferseF)
{
double const eps = 0.0000001;
double lon = 63.45421;
double x = MercatorBounds::LonToX(lon);
double lon1 = MercatorBounds::XToLon(x);
TEST_LESS(fabs(lon - lon1), eps, ("Too big round error"));
double lat = 34.28754;
double y = MercatorBounds::LatToY(lat);
double lat1 = MercatorBounds::YToLat(y);
TEST_LESS(fabs(lat - lat1), eps, ("Too big round error"));
TEST_LESS(fabs(MercatorBounds::maxX - MercatorBounds::maxY), eps, ("Non-square maxX and maxY"));
TEST_LESS(fabs(MercatorBounds::minX - MercatorBounds::minY), eps, ("Non-square minX and minY"));
}
UNIT_TEST(Mercator_ErrorToRadius)
{
double const points[] = { -85.0, -45.0, -10.0, -1.0, -0.003, 0.0, 0.003, 1.0, 10.0, 45.0, 85.0 };
double const error1 = 1.0; // 1 metre
double const error10 = 10.0; // 10 metres
for (size_t i = 0; i < ARRAY_SIZE(points); ++i)
{
for (size_t j = 0; j < ARRAY_SIZE(points); ++j)
{
double const lon = points[i];
double const lat = points[j];
m2::PointD const mercPoint(MercatorBounds::LonToX(lon), MercatorBounds::LatToY(lat));
m2::RectD const radius1 = MercatorBounds::MetresToXY(lon, lat, error1);
TEST(radius1.IsPointInside(mercPoint), (lat, lon));
TEST(radius1.Center().EqualDxDy(mercPoint, 1.0E-8), ());
m2::RectD const radius10 = MercatorBounds::MetresToXY(lon, lat, error10);
TEST(radius10.IsPointInside(mercPoint), (lat, lon));
TEST(radius10.Center().EqualDxDy(mercPoint, 1.0E-8), ());
TEST_EQUAL(m2::Add(radius10, radius1), radius10, (lat, lon));
TEST(radius10.IsPointInside(radius1.LeftTop()), (lat, lon));
TEST(radius10.IsPointInside(radius1.LeftBottom()), (lat, lon));
TEST(radius10.IsPointInside(radius1.RightTop()), (lat, lon));
TEST(radius10.IsPointInside(radius1.RightBottom()), (lat, lon));
}
}
}
UNIT_TEST(Mercator_Sample1)
{
LOG(LINFO, (MercatorBounds::XToLon(27.531491200000001385),
MercatorBounds::YToLat(64.392864299248202542)));
}
|
