Skip to content

Commit 2fa9cba

Browse files
asinghvi17asinghvi17
committed
Add testsets and testset reader for TG
These aren't run yet but it is nice to have the code, both to confirm that GeometryOps is correct, and to confirm that TGGeometry actually works.
1 parent f6092bd commit 2fa9cba

File tree

11 files changed

+2326
-0
lines changed

11 files changed

+2326
-0
lines changed

‎test/external/jts/jts_testset_reader.jl‎

Lines changed: 123 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,123 @@
1+
using XML
2+
3+
import WellKnownGeometry
4+
import GeoFormatTypes as GFT
5+
import GeometryOps as GO
6+
import GeoInterface as GI
7+
8+
"""
9+
jts_wkt_to_geom(wkt::String)
10+
11+
Convert a JTS WKT string to a GeometryOps geometry, via WellKnownGeometry.jl and GO.tuples.
12+
13+
The reason this exists is because WellKnownGeometry doesn't work well with newlines in subsidiary geometries,
14+
so this sanitizes the input before parsing and converting.
15+
"""
16+
function jts_wkt_to_geom(wkt::String)
17+
sanitized_wkt = join(strip.(split(wkt, "\n")), "")
18+
geom = GFT.WellKnownText(GFT.Geom(), sanitized_wkt)
19+
return GO.tuples(geom)
20+
end
21+
22+
struct TestItem{T}
23+
operation::String
24+
arg1::GO.GI.Wrappers.WrapperGeometry
25+
arg2::GO.GI.Wrappers.WrapperGeometry
26+
expected_result::T
27+
end
28+
29+
Base.show(io::IO, ::MIME"text/plain", item::TestItem) = print(io, "TestItem(operation = $(item.operation), expects $(GI.trait(item.expected_result)))")
30+
Base.show(io::IO, item::TestItem) = show(io, MIME"text/plain"(), item)
31+
32+
struct Case
33+
description::String
34+
geom_a::GO.GI.Wrappers.WrapperGeometry
35+
geom_b::GO.GI.Wrappers.WrapperGeometry
36+
items::Vector{TestItem}
37+
end
38+
39+
function load_test_cases(filepath::String)
40+
doc = read(filepath, XML.Node) # lazy parsing
41+
run = only(children(doc))
42+
test_cases = Case[]
43+
for case in children(run)
44+
if tag(case) != "case"
45+
continue
46+
end
47+
push!(test_cases, parse_case(case))
48+
end
49+
return test_cases
50+
end
51+
52+
function parse_case(case::XML.Node)
53+
description = value(only(children(case.children[1])))
54+
a = jts_wkt_to_geom(value(only(children(case.children[2]))))
55+
b = jts_wkt_to_geom(value(only(children(case.children[3]))))
56+
57+
items = TestItem[]
58+
for item in children(case)[4:end]
59+
ops = children(item)
60+
for op in ops
61+
op_attrs = XML.attributes(op)
62+
operation = op_attrs["name"]
63+
arg1 = op_attrs["arg1"]
64+
arg2 = op_attrs["arg2"]
65+
expected_result = jts_wkt_to_geom(value(only(op.children)))
66+
push!(items, TestItem(operation, lowercase(arg1) == "a" ? a : b, lowercase(arg2) == "a" ? a : b, expected_result))
67+
end
68+
end
69+
return Case(description, a, b, items)
70+
end
71+
72+
73+
testfile = "/Users/anshul/.julia/dev/geo/jts/modules/tests/src/test/resources/testxml/general/TestOverlayAA.xml"
74+
75+
cases = load_test_cases(testfile)
76+
77+
using Test
78+
79+
for case in cases
80+
@testset "$(case.description)" begin
81+
for item in case.items
82+
@testset "$(item.operation)" begin
83+
result = if item.operation == "union"
84+
GO.union(item.arg1, item.arg2; target = GO.PolygonTrait())
85+
elseif item.operation == "difference"
86+
GO.difference(item.arg1, item.arg2; target = GO.PolygonTrait())
87+
elseif item.operation == "intersection"
88+
GO.intersection(item.arg1, item.arg2; target = GO.PolygonTrait())
89+
elseif item.operation == "symdifference"
90+
continue
91+
else
92+
continue
93+
end
94+
95+
finalresult = if length(result) == 0
96+
nothing
97+
elseif length(result) == 1
98+
only(result)
99+
else
100+
GI.MultiPolygon(result)
101+
end
102+
103+
if isnothing(finalresult)
104+
@warn("No result")
105+
continue
106+
end
107+
108+
if GI.geomtrait(item.expected_result) isa Union{GI.MultiPolygonTrait, GI.PolygonTrait}
109+
difference_in_areas = GO.area(GO.difference(finalresult, item.expected_result; target = GO.PolygonTrait()))
110+
if difference_in_areas > 1
111+
@warn("Difference in areas: $(difference_in_areas)")
112+
f, a, p = poly(finalresult; label = "Final result (GO)", axis = (; title = "$(case.description) - $(item.operation)"))
113+
poly!(a, item.expected_result; label = "Expected result (JTS)")
114+
display(f)
115+
end
116+
# @test difference_in_areas < 1
117+
else
118+
@test_broken GO.equals(finalresult, item.expected_result)
119+
end
120+
end
121+
end
122+
end
123+
end

‎test/external/tg/data/LICENSE.md‎

Lines changed: 21 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,21 @@
1+
Reproduced from https://github.com/tidwall/tg/blob/main/LICENSE
2+
3+
> Copyright (c) 2023 Joshua J Baker
4+
>
5+
> Permission is hereby granted, free of charge, to any person obtaining a copy
6+
> of this software and associated documentation files (the "Software"), to deal
7+
> in the Software without restriction, including without limitation the rights
8+
> to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
9+
> copies of the Software, and to permit persons to whom the Software is
10+
> furnished to do so, subject to the following conditions:
11+
>
12+
> The above copyright notice and this permission notice shall be included in
13+
> all copies or substantial portions of the Software.
14+
>
15+
> THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
16+
> IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
17+
> FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
18+
> AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
19+
> LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
20+
> OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
21+
> THE SOFTWARE.

‎test/external/tg/data/gcol.jsonc‎

Lines changed: 62 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,62 @@
1+
[
2+
{
3+
"geoms": [
4+
"GEOMETRYCOLLECTION(POLYGON((1 1,4 1,4 4,1 4,1 2.5,1 1)))",
5+
"GEOMETRYCOLLECTION(POLYGON((0 2, 1 2,1 2.5,1 3, 0 3, 0 2)))"
6+
],
7+
"dims": [ 2, 2 ],
8+
"relate": [ "FF2F11212", "FF2F11212" ],
9+
"predicates": {
10+
"equals": ["F", "F"],
11+
"intersects": ["T", "T"],
12+
"disjoint": ["F", "F"],
13+
"contains": ["F", "F"],
14+
"within": ["F", "F"],
15+
"covers": ["F", "F"],
16+
"coveredby": ["F", "F"],
17+
"crosses": ["F", "F"],
18+
"overlaps": ["F", "F"],
19+
"touches": ["T", "T"]
20+
}
21+
},
22+
{
23+
"geoms": [
24+
"GEOMETRYCOLLECTION EMPTY",
25+
{"type":"Feature","id":"b","geometry":{"type":"LineString","coordinates":[[1.5,1.8],[1.5,1.4]]}}
26+
],
27+
"dims": [ -1, 1 ],
28+
"relate": [ "FFFFFF102", "FF1FF0FF2" ],
29+
"predicates": {
30+
"equals": ["F", "F"],
31+
"intersects": ["F", "F"],
32+
"disjoint": ["T", "T"],
33+
"contains": ["F", "F"],
34+
"within": ["F", "F"],
35+
"covers": ["F", "F"],
36+
"coveredby": ["F", "F"],
37+
"crosses": ["F", "F"],
38+
"overlaps": ["F", "F"],
39+
"touches": ["F", "F"]
40+
}
41+
},
42+
{
43+
"geoms": [
44+
"GEOMETRYCOLLECTION EMPTY",
45+
"GEOMETRYCOLLECTION(POLYGON((0 2, 1 2,1 2.5,1 3, 0 3, 0 2)))"
46+
],
47+
"dims": [ -1, 2 ],
48+
"relate": [ "FFFFFF212", "FF2FF1FF2" ],
49+
"predicates": {
50+
"equals": ["F", "F"],
51+
"intersects": ["F", "F"],
52+
"disjoint": ["T", "T"],
53+
"contains": ["F", "F"],
54+
"within": ["F", "F"],
55+
"covers": ["F", "F"],
56+
"coveredby": ["F", "F"],
57+
"crosses": ["F", "F"],
58+
"overlaps": ["F", "F"],
59+
"touches": ["F", "F"]
60+
}
61+
}
62+
]

0 commit comments

Comments
 (0)