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
|
/*
Copyright (c) 2017, 2025, Oracle and/or its affiliates.
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License, version 2.0,
as published by the Free Software Foundation.
This program is designed to work with certain software (including
but not limited to OpenSSL) that is licensed under separate terms,
as designated in a particular file or component or in included license
documentation. The authors of MySQL hereby grant you an additional
permission to link the program and your derivative works with the
separately licensed software that they have either included with
the program or referenced in the documentation.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License, version 2.0, for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*/
#include "my_config.h"
#include <gtest/gtest.h>
#include <memory> // unique_ptr
#include "sql/dd/dd.h"
#include "sql/dd/impl/types/spatial_reference_system_impl.h"
#include "sql/dd/properties.h"
#include "sql/dd/types/spatial_reference_system.h"
#include "sql/gis/geometries.h"
#include "sql/gis/geometries_cs.h"
#include "sql/gis/relops.h"
namespace gis_relops_unittest {
struct Cartesian_types {
typedef gis::Cartesian_point Point;
typedef gis::Cartesian_linestring Linestring;
typedef gis::Cartesian_linearring Linearring;
typedef gis::Cartesian_polygon Polygon;
typedef gis::Cartesian_geometrycollection Geometrycollection;
typedef gis::Cartesian_multipoint Multipoint;
typedef gis::Cartesian_multilinestring Multilinestring;
typedef gis::Cartesian_multipolygon Multipolygon;
static gis::Coordinate_system coordinate_system() {
return gis::Coordinate_system::kCartesian;
}
};
struct Geographic_types {
typedef gis::Geographic_point Point;
typedef gis::Geographic_linestring Linestring;
typedef gis::Geographic_linearring Linearring;
typedef gis::Geographic_polygon Polygon;
typedef gis::Geographic_geometrycollection Geometrycollection;
typedef gis::Geographic_multipoint Multipoint;
typedef gis::Geographic_multilinestring Multilinestring;
typedef gis::Geographic_multipolygon Multipolygon;
static gis::Coordinate_system coordinate_system() {
return gis::Coordinate_system::kGeographic;
}
};
template <typename Types>
class RelopsTest : public ::testing::Test {
protected:
dd::Spatial_reference_system_impl *m_srs;
public:
RelopsTest() {
switch (Types::coordinate_system()) {
case gis::Coordinate_system::kCartesian:
// Use SRID 0.
m_srs = nullptr;
break;
case gis::Coordinate_system::kGeographic:
// EPSG 4326, but with long-lat axes (E-N).
m_srs = dynamic_cast<dd::Spatial_reference_system_impl *>(
dd::create_object<dd::Spatial_reference_system>());
m_srs->set_id(4326);
m_srs->set_name("WGS 84");
m_srs->set_created(0UL);
m_srs->set_last_altered(0UL);
m_srs->set_organization("EPSG");
m_srs->set_organization_coordsys_id(4326);
m_srs->set_definition(
"GEOGCS[\"WGS 84\",DATUM[\"World Geodetic System "
"1984\",SPHEROID[\"WGS "
"84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],"
"AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY["
"\"EPSG\",\"8901\"]],UNIT[\"degree\",0.017453292519943278,"
"AUTHORITY[\"EPSG\",\"9122\"]],AXIS[\"Lon\",EAST],AXIS[\"Lat\","
"NORTH],AUTHORITY[\"EPSG\",\"4326\"]]");
m_srs->set_description("");
m_srs->parse_definition();
break;
}
}
~RelopsTest() override { delete m_srs; }
};
typedef ::testing::Types<Cartesian_types, Geographic_types> Types;
TYPED_TEST_SUITE(RelopsTest, Types);
// The purpose of this test is to cover all type combinations, not to check if
// the results are correct.
TYPED_TEST(RelopsTest, CodeCoverage) {
typename TypeParam::Geometrycollection gc;
gc.push_back(typename TypeParam::Point(0.0, 0.0));
typename TypeParam::Linestring ls;
ls.push_back(typename TypeParam::Point(0.0, 0.0));
ls.push_back(typename TypeParam::Point(1.0, 1.0));
gc.push_back(ls);
typename TypeParam::Polygon py;
typename TypeParam::Linearring lr;
lr.push_back(typename TypeParam::Point(0.0, 0.0));
lr.push_back(typename TypeParam::Point(1.0, 0.0));
lr.push_back(typename TypeParam::Point(1.0, 1.0));
lr.push_back(typename TypeParam::Point(0.0, 1.0));
lr.push_back(typename TypeParam::Point(0.0, 0.0));
py.push_back(lr);
gc.push_back(py);
typename TypeParam::Geometrycollection gc_empty;
gc.push_back(gc_empty);
typename TypeParam::Geometrycollection gc_inner;
gc_inner.push_back(typename TypeParam::Point(0.0, 0.0));
gc.push_back(gc_inner);
typename TypeParam::Multipoint mpt;
mpt.push_back(typename TypeParam::Point(0.0, 0.0));
gc.push_back(mpt);
typename TypeParam::Multilinestring mls;
mls.push_back(ls);
gc.push_back(mls);
typename TypeParam::Multipolygon mpy;
mpy.push_back(py);
gc.push_back(mpy);
for (auto g1 : gc) {
for (auto g2 : gc) {
bool result = false;
bool is_null = false;
gis::crosses(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::disjoint(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::equals(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::intersects(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_covered_by(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_disjoint(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_equals(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_intersects(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_overlaps(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_touches(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::mbr_within(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::overlaps(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::touches(this->m_srs, g1, g2, "unittest", &result, &is_null);
gis::within(this->m_srs, g1, g2, "unittest", &result, &is_null);
}
}
}
} // namespace gis_relops_unittest
|
