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
|
#include <tut/tut.hpp>
// geos
#include <geos_c.h>
#include "capi_test_utils.h"
namespace tut {
//
// Test Group
//
struct test_geosdifference_data : public capitest::utility {};
typedef test_group<test_geosdifference_data> group;
typedef group::object object;
group test_geosdifference("capi::GEOSDifference");
template<>
template<>
void object::test<1>()
{
geom1_ = fromWKT("LINESTRING (2 8, 10 8)");
geom2_ = fromWKT("LINESTRING (4 8, 6 8)");
GEOSSetSRID(geom1_, 4326);
result_ = GEOSDifference(geom1_, geom2_);
ensure(result_);
ensure_geometry_equals(result_, "MULTILINESTRING ((6 8, 10 8), (2 8, 4 8))");
ensure_equals(GEOSGetSRID(geom1_), GEOSGetSRID(result_));
}
/**
* Mixed GeometryCollection types permitted at a high-level
*/
template<>
template<>
void object::test<2>() {
GEOSGeometry* a = GEOSGeomFromWKT("GEOMETRYCOLLECTION (POINT (51 -1), LINESTRING (52 -1, 49 2))");
GEOSGeometry* b = GEOSGeomFromWKT("POINT (2 3)");
ensure(a);
ensure(b);
GEOSGeometry* ab = GEOSDifference(a, b);
GEOSGeometry* ba = GEOSDifference(b, a);
ensure(ab);
ensure(ba);
ensure_geometry_equals(ab, a);
ensure_geometry_equals(ba, b);
GEOSGeom_destroy(a);
GEOSGeom_destroy(b);
GEOSGeom_destroy(ab);
GEOSGeom_destroy(ba);
}
template<>
template<>
void object::test<3>()
{
geom1_ = fromWKT("CIRCULARSTRING (0 0, 1 1, 2 0)");
geom2_ = fromWKT("LINESTRING (1 0, 2 1)");
ensure(geom1_);
ensure(geom2_);
result_ = GEOSDifference(geom1_, geom2_);
ensure("curved geometry not supported", result_ == nullptr);
}
} // namespace tut
|
