-
Notifications
You must be signed in to change notification settings - Fork 219
/
Copy pathcentroid_multi.cpp
193 lines (155 loc) · 7.09 KB
/
centroid_multi.cpp
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
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
// Boost.Geometry (aka GGL, Generic Geometry Library)
// Unit Test
// Copyright (c) 2007-2015 Barend Gehrels, Amsterdam, the Netherlands.
// Copyright (c) 2008-2015 Bruno Lalande, Paris, France.
// Copyright (c) 2009-2015 Mateusz Loskot, London, UK.
// Copyright (c) 2017 Adam Wulkiewicz, Lodz, Poland.
// Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
// (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
// 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)
#include <algorithms/test_centroid.hpp>
#include <boost/core/ignore_unused.hpp>
#include <boost/geometry/core/point_order.hpp>
#include <boost/geometry/strategies/cartesian/centroid_average.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
// #define REPORT_RESULTS
template <typename P>
void test_2d(bool is_integer = false)
{
typedef typename bg::coordinate_type<P>::type ct;
boost::ignore_unused<ct>();
#ifdef REPORT_RESULTS
std::cout << std::endl << "type: " << typeid(ct).name() << " size: " << sizeof(ct) << std::endl;
#endif
test_centroid<bg::model::multi_point<P> >(
"MULTIPOINT((1 1),(3 3))",
2.0, 2.0);
test_centroid<bg::model::multi_point<P> >(
"MULTIPOINT((-1 -1),(-3 -3))",
-2.0, -2.0);
if (! is_integer)
{
// Only working for floating point:
test_centroid<bg::model::multi_point<P> >(
"MULTIPOINT((1 1),(2 3),(5 0))",
2.666666666666667, 1.33333);
test_centroid<bg::model::multi_linestring<bg::model::linestring<P> > >(
"MULTILINESTRING((0 0,0 2),(1 0,1 2))",
0.5, 1.0);
// degenerated
test_centroid<bg::model::multi_linestring<bg::model::linestring<P> > >(
"MULTILINESTRING((1 1,1 1))",
1.0, 1.0);
test_centroid<bg::model::multi_polygon<bg::model::polygon<P> > >(
"MULTIPOLYGON(((1 1,1 3,3 3,3 1,1 1)),((4 1,4 3,8 3,8 1,4 1)))",
4.666666666666667, 2.0);
test_centroid<bg::model::multi_polygon<bg::model::polygon<P> > >(
"MULTIPOLYGON(((2 1.3,2.4 1.7,2.8 1.8,3.4 1.2"
",3.7 1.6,3.4 2,4.1 3,5.3 2.6,5.4 1.2,4.9 0.8,2.9 0.7,2 1.3)),"
"((10 10,10 12,12 12,12 10,10 10)))",
7.338463104108615, 6.0606722055552407);
// degenerated
test_centroid<bg::model::multi_polygon<bg::model::polygon<P> > >(
"MULTIPOLYGON(((1 1,1 1,1 1,1 1,1 1)))",
1.0, 1.0);
test_centroid<bg::model::multi_polygon<bg::model::polygon<P> > >(
"MULTIPOLYGON(((1 1)))",
1.0, 1.0);
}
// Test using real-world polygon with large (Y) coordinates
// (coordinates can be used for integer and floating point point-types)
// Note that this will fail (overflow) if centroid calculation uses float
test_centroid<bg::model::multi_polygon<bg::model::polygon<P> > >(
"MULTIPOLYGON(((426062 4527794,426123 4527731"
",426113 4527700,426113 4527693,426115 4527671"
",426133 4527584,426135 4527569,426124 4527558"
",426103 4527547,426072 4527538,426003 4527535"
",425972 4527532,425950 4527531,425918 4527528"
",425894 4527517,425876 4527504,425870 4527484"
",425858 4527442,425842 4527414,425816 4527397"
",425752 4527384,425692 4527369,425658 4527349"
",425624 4527307,425605 4527260,425598 4527213"
",425595 4527167,425582 4527125,425548 4527064"
",425535 4527027,425537 4526990,425534 4526943"
",425525 4526904,425500 4526856,425461 4526811"
",425450 4526798,425381 4526823,425362 4526830"
",425329 4526848,425298 4526883,425291 4526897"
",425268 4526923,425243 4526945,425209 4526971"
",425172 4526990,425118 4527028,425104 4527044"
",425042 4527090,424980 4527126,424925 4527147"
",424881 4527148,424821 4527147,424698 4527125"
",424610 4527121,424566 4527126,424468 4527139"
",424426 4527141,424410 4527142,424333 4527130"
",424261 4527110,424179 4527073,424024 4527012"
",423947 4526987,423902 4526973,423858 4526961"
",423842 4526951,423816 4526935,423799 4526910"
",423776 4526905,423765 4526911,423739 4526927"
",423692 4526946,423636 4526976,423608 4527008"
",423570 4527016,423537 4527011,423505 4526996"
",423480 4526994,423457 4527012,423434 4527021"
",423367 4527008,423263 4526998,423210 4526993"
",423157 4526996,423110 4526994,423071 4526984"
",423048 4526984,423032 4526994,423254 4527613"
",423889 4528156,424585 4528050,425479 4527974"
",425795 4527867,426062 4527794)))",
424530.6059719588, 4527519.619367547);
}
template <typename P>
void test_exceptions()
{
using namespace bg::model;
typedef multi_polygon<polygon<P> > multi_polygon;
typedef multi_linestring<linestring<P> > multi_linestring;
// Empty multi-polygon
test_centroid_exception<multi_polygon>("MULTIPOLYGON()");
test_centroid_exception<multi_polygon>("MULTIPOLYGON(())");
test_centroid_exception<multi_polygon>("MULTIPOLYGON((), ())");
test_centroid_exception<multi_polygon>("MULTIPOLYGON((()), ())");
test_centroid_exception<multi_polygon>("MULTIPOLYGON(((), ()))");
// Empty multi-linestring
test_centroid_exception<multi_linestring>("MULTILINESTRING()");
test_centroid_exception<multi_linestring>("MULTILINESTRING(())");
test_centroid_exception<multi_linestring>("MULTILINESTRING((), ())");
}
template <typename P>
void test_empty()
{
using namespace bg::model;
typedef multi_polygon<polygon<P> > multi_polygon;
typedef multi_linestring<linestring<P> > multi_linestring;
// Multi-linestring with empty linestring
test_centroid<multi_linestring>(
"MULTILINESTRING((), (0 0))",
0.0, 0.0);
test_centroid<multi_linestring>(
"MULTILINESTRING((0 0, 1 0), ())",
0.5, 0.0);
// Multi-polygon with empty polygon
test_centroid<multi_polygon>(
"MULTIPOLYGON((()), ((0 0)))",
0.0, 0.0);
test_centroid<multi_polygon>(
"MULTIPOLYGON(((0 0, 1 0, 1 1, 0 1, 0 0)), (()))",
0.5, 0.5);
// Multi-polygon with empty interior ring
test_centroid<multi_polygon>(
"MULTIPOLYGON(((0 0, 1 0, 1 1, 0 1, 0 0), ()))",
0.5, 0.5);
test_centroid<multi_polygon>(
"MULTIPOLYGON((()), ((0 0, 1 0, 1 1, 0 1, 0 0), ()))",
0.5, 0.5);
}
int test_main(int, char* [])
{
test_2d<bg::model::d2::point_xy<float> >();
test_2d<bg::model::d2::point_xy<double> >();
test_2d<bg::model::d2::point_xy<long int> >(true);
//test_2d<bg::model::d2::point_xy<long long> >(true);
//test_2d<bg::model::d2::point_xy<long double> >();
test_exceptions<bg::model::d2::point_xy<double> >();
test_empty<bg::model::d2::point_xy<double> >();
return 0;
}