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
|
//
// Test Suite for C-API GEOSFrechetDistance
#include <tut/tut.hpp>
// geos
#include <geos_c.h>
#include "capi_test_utils.h"
namespace tut {
//
// Test Group
//
// Common data used in test cases.
struct test_capigeosfrechetdistance_data : public capitest::utility {};
typedef test_group<test_capigeosfrechetdistance_data> group;
typedef group::object object;
group test_capigeosfrechetdistance_group("capi::GEOSFrechetDistance");
//
// Test Cases
//
template<>
template<>
void object::test<1>()
{
geom1_ = fromWKT("LINESTRING (0 0, 100 0)");
geom2_ = fromWKT("LINESTRING (0 0, 50 50, 100 0)");
double dist;
int ret = GEOSFrechetDistance(geom1_, geom2_, &dist);
ensure_equals(ret, 1);
ensure_distance(dist, 70.7106781186548, 1e-12);
}
template<>
template<>
void object::test<2>()
{
geom1_ = fromWKT("LINESTRING (0 0, 100 0)");
geom2_ = fromWKT("LINESTRING (0 0, 50 50, 100 0)");
double dist;
int ret = GEOSFrechetDistanceDensify(geom1_, geom2_, 0.5, &dist);
ensure_equals(ret, 1);
ensure_distance(dist, 50., 1e-12);
}
// No crash with tiny densify fraction
// https://trac.osgeo.org/geos/ticket/1086
template<>
template<>
void object::test<3>()
{
double dist = 0;
geom1_ = fromWKT("LINESTRING (0 0, 3 7, 5 5)");
geom2_ = fromWKT("LINESTRING (0 0, 9 1, 2 2)");
int ret = GEOSFrechetDistanceDensify(geom1_, geom2_, 1e-40, &dist);
ensure_equals(ret, 0);
}
template<>
template<>
void object::test<5>()
{
geom1_ = fromWKT("CIRCULARSTRING (0 0, 1 1, 2 0)");
geom2_ = fromWKT("LINESTRING (1 2, 2 2)");
ensure(geom1_);
ensure(geom2_);
double dist;
ensure_equals("curved geometry not supported", GEOSFrechetDistance(geom1_, geom2_, &dist), 0);
ensure_equals("curved geometry not supported", GEOSFrechetDistance(geom2_, geom1_, &dist), 0);
}
} // namespace tut
|
