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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
|
// Boost.Geometry
// Unit Test
// Copyright (c) 2017, Oracle and/or its affiliates.
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
// Use, modification and distribution is subject to the Boost Software License,
// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
// http://www.boost.org/LICENSE_1_0.txt)
#ifndef BOOST_GEOMETRY_TEST_SRS_CHECK_GEOMETRY_HPP
#define BOOST_GEOMETRY_TEST_SRS_CHECK_GEOMETRY_HPP
#include <geometry_test_common.hpp>
#include <boost/geometry/core/access.hpp>
#include <boost/geometry/core/coordinate_type.hpp>
#include <boost/geometry/core/tag.hpp>
#include <boost/geometry/core/tags.hpp>
#include <boost/geometry/io/wkt/read.hpp>
#include <boost/geometry/views/detail/indexed_point_view.hpp>
#include <boost/range/begin.hpp>
#include <boost/range/end.hpp>
#include <boost/range/size.hpp>
#include <boost/range/value_type.hpp>
namespace test
{
struct check_point
{
template <typename Point, typename T>
static void apply(Point const& point1, Point const& point2, T tol)
{
typename bg::coordinate_type<Point>::type
x1 = bg::get<0>(point1),
y1 = bg::get<1>(point1),
x2 = bg::get<0>(point2),
y2 = bg::get<1>(point2);
BOOST_CHECK_CLOSE(x1, x2, tol);
BOOST_CHECK_CLOSE(y1, y2, tol);
}
};
template <typename Policy = check_point>
struct check_range
{
template <typename Range, typename T>
static void apply(Range const& range1, Range const& range2, T tol)
{
size_t range1_count = boost::size(range1);
size_t range2_count = boost::size(range2);
BOOST_CHECK_EQUAL(range1_count, range2_count);
if (range1_count == range2_count)
{
apply(boost::begin(range1), boost::end(range1),
boost::begin(range2), tol);
}
}
template <typename It, typename T>
static void apply(It first1, It last1, It first2, T tol)
{
for ( ; first1 != last1 ; ++first1, ++first2)
Policy::apply(*first1, *first2, tol);
}
};
template <typename Geometry, typename Tag = typename bg::tag<Geometry>::type>
struct check_geometry_impl
{};
template <typename Point>
struct check_geometry_impl<Point, bg::point_tag>
: check_point
{};
template <typename Segment>
struct check_geometry_impl<Segment, bg::segment_tag>
{
template <typename T>
static void apply(Segment const& g1, Segment const& g2, T tol)
{
bg::detail::indexed_point_view<Segment const, 0> p1(g1);
bg::detail::indexed_point_view<Segment const, 1> p2(g1);
bg::detail::indexed_point_view<Segment const, 0> q1(g2);
bg::detail::indexed_point_view<Segment const, 1> q2(g2);
check_point::apply(p1, q1, tol);
check_point::apply(p2, q2, tol);
}
};
template <typename MultiPoint>
struct check_geometry_impl<MultiPoint, bg::multi_point_tag>
: check_range<>
{};
template <typename Linestring>
struct check_geometry_impl<Linestring, bg::linestring_tag>
: check_range<>
{};
template <typename MultiLinestring>
struct check_geometry_impl<MultiLinestring, bg::multi_linestring_tag>
: check_range< check_range<> >
{};
template <typename Ring>
struct check_geometry_impl<Ring, bg::ring_tag>
: check_range<>
{};
template <typename Polygon>
struct check_geometry_impl<Polygon, bg::polygon_tag>
{
template <typename T>
static void apply(Polygon const& g1, Polygon const& g2, T tol)
{
check_range<>::apply(bg::exterior_ring(g1), bg::exterior_ring(g2), tol);
check_range< check_range<> >::apply(bg::interior_rings(g1), bg::interior_rings(g2), tol);
}
};
template <typename MultiPolygon>
struct check_geometry_impl<MultiPolygon, bg::multi_polygon_tag>
: check_range
<
check_geometry_impl
<
typename boost::range_value<MultiPolygon>::type,
bg::polygon_tag
>
>
{};
template <typename Geometry, typename T>
inline void check_geometry(Geometry const& g1, Geometry const& g2, T tol)
{
check_geometry_impl<Geometry>::apply(g1, g2, tol);
}
template <typename Geometry, typename T>
inline void check_geometry(Geometry const& g1, std::string const& wkt2, T tol)
{
Geometry g2;
bg::read_wkt(wkt2, g2);
check_geometry_impl<Geometry>::apply(g1, g2, tol);
}
} // namespace test
#endif // BOOST_GEOMETRY_TEST_SRS_CHECK_GEOMETRY_HPP
|
