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
|
# Copyright (C) 2010, Nokia <ivan.frade@nokia.com>
# Copyright (C) 2019, Sam Thursfield <sam@afuera.me.uk>
#
# This program is free software; you can redistribute it and/or
# modify it under the terms of the GNU General Public License
# as published by the Free Software Foundation; either version 2
# of the License, or (at your option) any later version.
#
# 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 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 Street, Fifth Floor, Boston, MA
# 02110-1301, USA.
#
"""
Test the distance-calculation functions in SPARQL.
"""
import unittest as ut
import fixtures
POINT_COORDS = [
(0, 0), (1, 1), (2, 2), (3, 3), (4, 4)
]
class TestDistanceFunctions (fixtures.TrackerSparqlDirectTest):
"""
Insert some points and get the distance between them.
"""
def setUp(self):
self.counter = 0
for lat, log in POINT_COORDS:
insert = """
INSERT {
<%s> a slo:GeoLocation ;
slo:longitude %d ;
slo:latitude %d .
}
""" % ("point://test/point/" + str(self.counter), log, lat)
self.tracker.update(insert)
self.counter += 1
def tearDown(self):
for i in range(0, self.counter):
delete = """
DELETE {
<%s> a rdfs:Resource.
}
""" % ("point://test/point/" + str (i))
self.tracker.update(delete)
def get_distance_between_points(self, sum_func, id1, id2):
assert 0 <= id1 <= len(POINT_COORDS)
assert 0 <= id2 <= len(POINT_COORDS)
assert sum_func == "cartesian" or sum_func == "haversine"
query_1_to_2 = """
SELECT xsd:integer(tracker:%s-distance(?lat1,?lat2,?lon1,?lon2))
WHERE {
<point://test/point/%d> a slo:GeoLocation ;
slo:latitude ?lat1 ;
slo:longitude ?lon1 .
<point://test/point/%d> a slo:GeoLocation ;
slo:latitude ?lat2 ;
slo:longitude ?lon2 .
}
""" % (sum_func, id1, id2)
result = self.tracker.query(query_1_to_2)
return int(result[0][0])
def test_distance_cartesian_symmetry(self):
"""
setUp: Insert 5 points in the pre-defined coordinates
1. TEST: Check cartesian distance from point A to B, and from B to A
(should be the same :P)
tearDown: Remove the test points inserted before
"""
a_to_b = self.get_distance_between_points("cartesian", 1, 2)
assert a_to_b == 204601
b_to_a = self.get_distance_between_points("cartesian", 2, 1)
assert b_to_a == 204601
assert a_to_b == b_to_a
def test_distance_haversine_symmetry(self):
"""
setUp: Insert 5 points in the pre-defined coordinates
1. TEST: Check cartesian distance from point A to B, and from B to A
(should be the same :P)
tearDown: Remove the test points inserted before
"""
a_to_b = self.get_distance_between_points("haversine", 1, 2)
assert a_to_b == 157225
b_to_a = self.get_distance_between_points("haversine", 2, 1)
assert b_to_a == 157225
assert a_to_b == b_to_a
def test_distance_cartesian_proportion(self):
d_1_to_2 = self.get_distance_between_points("cartesian", 1, 2)
d_2_to_3 = self.get_distance_between_points("cartesian", 2, 3)
d_3_to_4 = self.get_distance_between_points("cartesian", 3, 4)
assert d_1_to_2 > d_2_to_3 > d_3_to_4
def test_distance_haversine_proportion(self):
d_1_to_2 = self.get_distance_between_points("haversine", 1, 2)
d_2_to_3 = self.get_distance_between_points("haversine", 2, 3)
d_3_to_4 = self.get_distance_between_points("haversine", 3, 4)
assert d_1_to_2 > d_2_to_3 > d_3_to_4
def test_distance_different(self):
d_2_to_3h = self.get_distance_between_points("haversine", 2, 3)
d_2_to_3c = self.get_distance_between_points("cartesian", 2, 3)
assert d_2_to_3h < d_2_to_3c
if __name__ == '__main__':
fixtures.tracker_test_main()
|
