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
|
#include <tut/tut.hpp>
// geos
#include <geos_c.h>
#include "capi_test_utils.h"
namespace tut {
//
// Test Group
//
struct test_geosisvalid_data : public capitest::utility {};
typedef test_group<test_geosisvalid_data> group;
typedef group::object object;
group test_geosisvalid("capi::GEOSisValid");
template<>
template<>
void object::test<1>()
{
GEOSGeometry* input = GEOSGeomFromWKT("LINESTRING (1 2, 4 5, 9 -2)");
ensure_equals(1, GEOSisValid(input));
GEOSGeom_destroy(input);
}
template<>
template<>
void object::test<2>()
{
GEOSGeometry* input = GEOSGeomFromWKT("POLYGON ((0 0, 1 0, 0 1, 1 1, 0 0))");
ensure_equals(0, GEOSisValid(input));
GEOSGeom_destroy(input);
}
// Unclosed polygon
template<>
template<>
void object::test<3>()
{
GEOSCoordSequence* shell_seq = GEOSCoordSeq_create(4, 2);
double shell_coords[] = {0,0, 0,10, 10,10, 10,0};
for (unsigned int i = 0; i < 4; i++) {
GEOSCoordSeq_setXY(shell_seq, i, shell_coords[2*i], shell_coords[2*i+1]);
}
// Unclosed ring fails in construction
// Linear ring takes ownership of coordinate sequence and
// frees it when construction fails
GEOSGeometry* shell = GEOSGeom_createLinearRing(shell_seq);
ensure(shell == nullptr);
// So we end up handing nullptr into create polygon
GEOSGeometry* polygon = GEOSGeom_createPolygon(shell, nullptr, 0);
ensure(polygon == nullptr);
// Which also returns nullptr and that goes into isvalid
char isvalid = GEOSisValid(polygon);
// Which causes an exception that we catch
ensure_equals(2, isvalid);
}
template<>
template<>
void object::test<5>()
{
input_ = fromWKT("CURVEPOLYGON (COMPOUNDCURVE( CIRCULARSTRING (0 0, 1 1, 2 0), (2 0, 1 1)))");
ensure(input_ != nullptr);
char ret = GEOSisValid(input_);
ensure_equals("error raised on curved geometry", ret, 2);
}
} // namespace tut
|
